Leray-α Regularization of the Smagorinsky-Closed Filtered Equations for Turbulent Jets at High Reynolds Numbers

The article reports on blending of the Leray-α regularization with the conventional Smagorinsky subgrid-scale closure as an option for large-eddy-simulation of turbulent flows at very high Reynolds number on coarse meshes. The model has been tested in the self-similar far-field region of a jet at a range of Reynolds numbers spanning over two decades (4×10³, 4×10⁴ and 4×10⁵) on two very coarse meshes of 2×10⁵ and 3×10⁴ mesh cells. The results are compared with the well-resolved DNS for $Re_D=4\times 10^3$ on 15 million cells and experimental data for higher Re numbers. While the pure Leray-α can fail badly at high Re numbers on very coarse meshes, a blending of the two strategies by adding a small amount of extra-dissipation performs well even at a huge jet Reynolds number of $Re_D=4\times 10^5$ on a very coarse mesh (2×10⁵ cells), despite the ratio of the typical mesh spacing to the Kolmogorov length exceeding 300. It is found that the main prerequisite for successful LES, both for the classic Smagorinsky and the blended Leray-α/Smagorinsky model, is to resolve the shear-length $L_s=\sqrt{\varepsilon/{\cal S}^3}$ (where ${\cal S}$ is the shear-rate modulus), defined by the constraint Δ/L _s ?<?1, where Δ is the typical mesh-cell size. For the mixed Leray-α/Smagorinsky model the regularization parameter should also be related to the shear-length rather than the local mesh size or Reynolds number, for which we propose a guide criterion α?=?0.15÷0.3 L _s .     

Effect of an impinging shock wave on a partially opened door

The behavior of an aluminum door hanging at the exit of an open shock tube at different angles, from 5 $^\circ $ to 85 $^\circ $ , and thereby providing partially open space for the exiting flow, was investigated experimentally. Experiments were conducted with an incident shock wave Mach number of $M_\mathrm{is}=1.1$ impinging on the partially opened door. Both pressure measurements in the vicinity of the door, on its center and inside the shock tube, and schlieren visualization were undertaken for studying the door movement and its maximum opening angle relative to its initial position. It was found that for an initial opening angle smaller than 25 $^\circ $ the door opened completely while for larger angles its motion is marginal. In addition, for an initial door opening angle of about 10 $^\circ $ the lowest pressures were recorded inside the shock tube behind the evolving waves after exiting of the incident shock wave. The present experimental results may be useful to numerical studies of fluid–structure interactions, e.g., in designing safety valves in jet engines. Such a device is needed for preventing rupture in the case when a sudden overpressure pulse is generated inside the aircraft engine compartment.     

Global Compensated Compactness Theorem for General Differential Operators of First Order

Let A ₁(x, D) and A ₂(x, D) be differential operators of the first order acting on l-vector functions ${u= (u_1, \ldots, u_l)}$ in a bounded domain ${\Omega \subset \mathbb{R}^{n}}$ with the smooth boundary ${\partial\Omega}$ . We assume that the H ¹-norm ${\|u\|_{H^{1}(\Omega)}}$ is equivalent to ${\sum_{i=1}^2\|A_iu\|_{L^2(\Omega)} + \|B_1u\|_{H^{\frac{1}{2}}(\partial\Omega)}}$ and ${\sum_{i=1}^2\|A_iu\|_{L^2(\Omega)} + \|B_2u\|_{H^{\frac{1}{2}}(\partial\Omega)}}$ , where B _i  = B _i (x, ν) is the trace operator onto ${\partial\Omega}$ associated with A _i (x, D) for i = 1, 2 which is determined by the Stokes integral formula (ν: unit outer normal to ${\partial\Omega}$ ). Furthermore, we impose on A ₁ and A ₂ a cancellation property such as ${A_1A_2^{\prime}=0}$ and ${A_2A_1^{\prime}=0}$ , where ${A^{\prime}_i}$ is the formal adjoint differential operator of A _i (i = 1, 2). Suppose that ${\{u_m\}_{m=1}^{\infty}}$ and ${\{v_m\}_{m=1}^{\infty}}$ converge to u and v weakly in ${L^2(\Omega)}$ , respectively. Assume also that ${\{A_{1}u_m\}_{m=1}^{\infty}}$ and ${\{A_{2}v_{m}\}_{m=1}^{\infty}}$ are bounded in ${L^{2}(\Omega)}$ . If either ${\{B_{1}u_m\}_{m=1}^{\infty}}$ or ${\{B_{2}v_m\}_{m=1}^{\infty}}$ is bounded in ${H^{\frac{1}{2}}(\partial\Omega)}$ , then it holds that ${\int_{\Omega}u_m\cdot v_m \,{\rm d}x \to \int_{\Omega}u\cdot v \,{\rm d}x}$ . We also discuss a corresponding result on compact Riemannian manifolds with boundary.     

Regularity and Asymptotic Behavior of Nonlinear Stefan Problems

We study the following nonlinear Stefan problem $$\left\{\begin{aligned}\!\!&u_t\,-\,d\Delta u = g(u) & &\quad{\rm for}\,x\,\in\,\Omega(t), t > 0, \\ & u = 0 \, {\rm and} u_t = \mu|\nabla_{x} u|^{2} &&\quad {\rm for}\,x\,\in\,\Gamma(t), t > 0, \\ &u(0, x) = u_{0}(x) &&\quad {\rm for}\,x\,\in\,\Omega_0,\end{aligned} \right.$$ where ${\Omega(t) \subset \mathbb{R}^{n}}$ ( ${n \geqq 2}$ ) is bounded by the free boundary ${\Gamma(t)}$ , with ${\Omega(0) = \Omega_0}$ , μ and d are given positive constants. The initial function u ₀ is positive in ${\Omega_0}$ and vanishes on ${\partial \Omega_0}$ . The class of nonlinear functions g(u) includes the standard monostable, bistable and combustion type nonlinearities. We show that the free boundary ${\Gamma(t)}$ is smooth outside the closed convex hull of ${\Omega_0}$ , and as ${t \to \infty}$ , either ${\Omega(t)}$ expands to the entire ${\mathbb{R}^n}$ , or it stays bounded. Moreover, in the former case, ${\Gamma(t)}$ converges to the unit sphere when normalized, and in the latter case, ${u \to 0}$ uniformly. When ${g(u) = au - bu^2}$ , we further prove that in the case ${\Omega(t)}$ expands to ${{\mathbb R}^n}$ , ${u \to a/b}$ as ${t \to \infty}$ , and the spreading speed of the free boundary converges to a positive constant; moreover, there exists ${\mu^* \geqq 0}$ such that ${\Omega(t)}$ expands to ${{\mathbb{R}}^n}$ exactly when ${\mu > \mu^*}$ .     

A-Priori Direct Numerical Simulation Modelling of Co-variance Transport in Turbulent Stratified Flames

Three-dimensional Direct Numerical Simulations of statistically planar turbulent stratified flames at global equivalence ratios <???>?=?0.7 and <???>?=?1.0 have been carried out to analyse the statistical behaviour of the transport of co-variance of the fuel mass fraction Y _F and mixture fraction ξ (i.e. $\widetilde{Y_F^{\prime\prime} \xi ^{\prime\prime}}={\overline {\rho Y_F^{\prime\prime} \xi^{\prime\prime}} } \Big/ {\overline \rho })$ for Reynolds Averaged Navier Stokes simulations where $\overline q $ , $\tilde{q} ={\overline {\rho q} } \big/ {\overline \rho }$ and $q^{\prime\prime}= q-\tilde{q}$ are Reynolds averaged, Favre mean and Favre fluctuation of a general quantity q with ρ being the gas density and the overbar suggesting a Reynolds averaging operation. It has been found that existing algebraic expressions may not capture the statistical behaviour of $\widetilde{Y_F^{\prime\prime} \xi^{\prime\prime}}$ with sufficient accuracy in low Damköhler number combustion and therefore, a transport equation for $\widetilde{Y_F^{\prime\prime} \xi^{\prime\prime}}$ may need to be solved. The statistical behaviours of $\widetilde{Y_F^{\prime\prime} \xi^{\prime\prime}}$ and the unclosed terms of its transport equation (i.e. the terms originating from turbulent transport T ₁ , reaction rate T ₄ and molecular dissipation $\left( {-D_2 } \right))$ have been analysed in detail. The contribution of T ₁ remains important for all cases considered here. The term T ₄ acts as a major contributor in <???>?=?1.0 cases, but plays a relatively less important role in <???>?=?0.7 cases, whereas the term $\left( {-D_2 } \right)$ acts mostly as a leading order sink. Through an a-priori DNS analysis, the performances of the models for T ₁ , T ₄ and $\left( {-D_2 } \right)$ have been addressed in detail. A model has been identified for the turbulent transport term T ₁ which satisfactorily predicts the corresponding term obtained from DNS data. The models for T ₄ , which were originally proposed for high Damköhler number flames, have been modified for low Damköhler combustion. Predictions of the modified models are found to be in good agreement with T ₄ obtained from DNS data. It has been found that existing algebraic models for $D_2 =2\overline {\rho D\nabla Y_F^{\prime\prime} \nabla \xi^{\prime\prime}} $ (where D is the mass diffusivity) are not sufficient for low Damköhler number combustion and therefore, a transport equation may need to be solved for the cross-scalar dissipation rate $\widetilde{\varepsilon }_{Y\xi } ={\overline {\rho D\nabla Y_F^{\prime\prime} \nabla \xi^{\prime\prime}} } \big/ {\overline \rho }$ for the closure of the $\widetilde{Y_F^{\prime\prime} \xi^{\prime\prime}}$ transport equation.     

Elongational rheology of NIPAM-based hydrogels

Hydrogels of different composition based on the copolymerization of N-isopropyl acrylamide and surfmers of different chemical structure were tested in elongation using Hencky/real definitions for stress, strain, and strain rate, offering a more scientific insight into the effect of deformation on the properties. In a range between $\dot {\varepsilon }=10$ and 0.01 s $^{-1}$ , the material properties are independent of strain rate and show a very clear strain hardening with a “brittle” sudden fracture. The addition of surfmer increases the strain at break $\varepsilon _{\mathrm {H}}^{\max }$ and at the same time leads to a failure of hyperelastic models. The samples can be stretched up to Hencky strains $\varepsilon _{\mathrm {H}}^{\max }$ between 0.6 and 2.5, depending on the molecular structure, yielding linear Young’s moduli E $_{0}$ between 2,700 and 39,000 Pa. The strain-rate independence indicates an ideal rubberlike behavior and fracture in a brittle-like fashion. The resulting stress at break $\sigma _{\textrm max}$ can be correlated with $\varepsilon _{\mathrm {H}}^{\max } $ and $E_{0}$ as well as with the solid molar mass between the cross-linking points $M_{\mathrm {c}}^{\textrm {solids}} $ , derived from $E_{0}$ .     

On the Stokes and Oseen Problems with Singular Data

We consider the steady Stokes and Oseen problems in bounded and exterior domains of ${\mathbb{R}^n}$ of class C ^k-1,1 (n = 2, 3; k ≥ 2). We prove existence and uniqueness of a very weak solution for boundary data a in ${W^{2-k-1/q,q} (\partial\Omega)}$ . If ${\Omega}$ is of class ${C^\infty}$ , we can assume a to be a distribution on ${\partial\Omega}$ .     

Drainage of Capillary-Trapped Oil by an Immiscible Gas: Impact of Transient and Steady-State Water Displacement on Three-Phase Oil Permeability

In a previous paper (Dehghanpour et al., Phys Rev E 83:065302, 2011a), we showed that relative permeability of mobilized oil, $k_\mathrm{ro}$ , measured during tertiary gravity drainage, is significantly higher than that of the same oil saturation in other tests where oil is initially a continuous phase. We also showed that tertiary $k_\mathrm{ro}$ strongly correlates to both water saturation, $S_\mathrm{w}$ , water flux (water relative permeability), $k_\mathrm{rw}$ , and the change in water saturation with time, $\mathrm{d}S_\mathrm{w}/\mathrm{d}t$ . To develop a model and understanding of the enhanced oil transport, identifying which of these parameters ( $S_\mathrm{w},\,k_{\mathrm{rw}}$ , or $\mathrm{d}S_\mathrm{w}/\mathrm{d}t$ ) plays the controlling role is necessary, but in the previous experiments these could not be deconvolved. To answer the remaining question, we conduct specific three-phase displacement experiments in which $k_{\mathrm{rw}}$ is controlled by applying a fixed water influx, and $S_\mathrm{w}$ develops naturally. We obtain $k_{\mathrm{ro}}$ by using the saturation data measured in time and space. The results suggest that steady-state water influx, in contrast to transient water displacement, does not enhance $k_{\mathrm{ro}}$ . Instead, reducing water influx rate results in excess oil flow. Furthermore, according to our pore scale hydraulic conductivity calculations, viscous coupling and fluid positioning do not sufficiently explain the observed correlation between $k_{\mathrm{ro}}$ and $S_{\mathrm{w}}$ . We conclude that tertiary $k_{\mathrm{ro}}$ is controlled by the oil mobilization rate, which in turn is linked to the rate of water saturation decrease with time, $\mathrm{d}S_\mathrm{w}/\mathrm{d}t$ . Finally, we develop a simple model which relates tertiary $k_{\mathrm{ro}}$ to transient two-phase gas/water relative permeability.     

On Formation of a Locally Self-Similar Collapse in the Incompressible Euler Equations

The paper addresses the question of the existence of a locally self-similar blow-up for the incompressible Euler equations. Several exclusion results are proved based on the L ^p -condition for velocity or vorticity and for a range of scaling exponents. In particular, in N dimensions if in self-similar variables ${u \in L^p}$ and ${u \sim \frac{1}{t^{\alpha/(1+\alpha)}}}$ , then the blow-up does not occur, provided ${\alpha > N/2}$ or ${-1 < \alpha \leq N\,/p}$ . This includes the L ³ case natural for the Navier–Stokes equations. For ${\alpha = N\,/2}$ we exclude profiles with asymptotic power bounds of the form ${ |y|^{-N-1+\delta} \lesssim |u(y)| \lesssim |y|^{1-\delta}}$ . Solutions homogeneous near infinity are eliminated, as well, except when homogeneity is scaling invariant.     

Double-Diffusive Natural Convection in Anisotropic Porous Medium Bounded by Finite Thickness Walls: Validity of Local Thermal Equilibrium Assumption

Double-diffusive natural convection in fluid-saturated porous medium inside a vertical enclosure bounded by finite thickness walls with opposing temperature, concentration gradients on vertical walls as well as adiabatic and impermeable horizontal ones has been performed numerically. The Darcy model was used to predict fluid flow inside the porous material, while thermal fields are simulated based on two-energy equations for fluid and solid phases on the basis of a local thermal non-equilibrium model. Computations have been performed for different controlling parameters such as the buoyancy ratio $N$ , the Lewis number Le, the anisotropic permeability ratio $R_\mathrm{p}$ , the fluid-to-solid thermal conductivity ratio $R_\mathrm{c}$ , the interphase heat transfer coefficient $\mathcal{H}$ , the ratio of the wall thickness to its height $D$ , the wall-to-porous medium thermal diffusivity ratio $R_\mathrm{w}$ , and the solid-to-fluid heat capacity ratio $\gamma $ . Thus, the effects of the controlling parameters on heat and mass transfer characteristics are discussed in detail. Moreover, the validity domain of the local thermal equilibrium (LTE) assumption has been delimited for different set of the governing parameters. It has been shown that Le has a noticeable significant effect on fluid temperature profiles and that higher $N$ values lead to a significant enhancement in heat and mass transfer rates. Moreover, for higher $\mathcal{H}, R_\mathrm{c}$ , $R_\mathrm{p}, R_\mathrm{w}$ , or $D$ values and/or lower $\gamma $ values, the solid and fluid phases tend toward LTE.     

An experimental study on turbulence modification in the near-wall boundary layer of a dilute gas-particle channel flow

Turbulence modifications of a dilute gas-particle flow are experimentally investigated in the lower boundary layer of a horizontal channel by means of a simultaneous two-phase PIV measurement technique. The measurements are conducted in the near-wall region with y ⁺?<?250 at Re _τ (based on the wall friction velocity u _τ and half channel height h)?=?430. High spatial resolution and small interrogation window are used to minimize the PIV measurement uncertainty due to the velocity gradient near the wall. Polythene beads with the diameter of 60?μm (d _p ⁺ ?=?1.71, normalized by the fluid kinematic viscosity ν and u _τ) are used as dispersed phase, and three low mass loading ratios (Φ _m ) ranging from 10^?4 to 10^?3 are tested. It is found that the addition of the particles noticeably modifies the mean velocity and turbulent intensities of the gas-phase, as well as the turbulence coherent structures, even at Φ _m ?=?0.025?%. Particle inertia changes the viscous sublayer of the gas turbulence with a smaller thickness and a larger streamwise velocity gradient, which increases the peak value of the streamwise fluctuation velocity ( $ u_{\text{rms}}^{ + } $ ) of the gas-phase with its location shifting to the wall. Particle sedimentation increases the roughness of the bottom wall, which significantly increases the wall-normal fluctuation velocity ( $ v_{\text{rms}}^{ + } $ ) and Reynolds shear stress ( $ - \langle u^{ \prime } v^{\prime } \rangle^{ + } $ ) of the gas-phase in the inner region of the boundary layer (y ⁺?<?10). Under effect of particle–wall collision, the Q2 events (ejections) of the gas-phase are slightly increased by particles, while the Q4 events (sweeps) are obviously decreased. The spatial scale of the coherent structures near the wall shrinks remarkably with the presence of the particles, which may be attributed to the intensified crossing-trajectory effects due to particle saltation near the bottom wall. Meanwhile, the $ v_{\text{rms}}^{ + } $ and $ - \langle u^{ \prime } v^{\prime } \rangle^{ + } $ of the gas-phase are significantly reduced in the outer region of the boundary layer (y ⁺?>?20).     

Semianalytical Solution for \text{ CO}_{2} Plume Shape and Pressure Evolution During \text{ CO}_{2} Injection in Deep Saline Formations

The injection of supercritical carbon dioxide ( $\text{ CO}_{2})$ in deep saline aquifers leads to the formation of a $\text{ CO}_{2}$ rich phase plume that tends to float over the resident brine. As pressure builds up, $\text{ CO}_{2}$ density will increase because of its high compressibility. Current analytical solutions do not account for $\text{ CO}_{2}$ compressibility and consider a volumetric injection rate that is uniformly distributed along the whole thickness of the aquifer, which is unrealistic. Furthermore, the slope of the $\text{ CO}_{2}$ pressure with respect to the logarithm of distance obtained from these solutions differs from that of numerical solutions. We develop a semianalytical solution for the $\text{ CO}_{2}$ plume geometry and fluid pressure evolution, accounting for $\text{ CO}_{2}$ compressibility and buoyancy effects in the injection well, so $\text{ CO}_{2}$ is not uniformly injected along the aquifer thickness. We formulate the problem in terms of a $\text{ CO}_{2}$ potential that facilitates solution in horizontal layers, with which we discretize the aquifer. Capillary pressure is considered at the interface between the $\text{ CO}_{2}$ rich phase and the aqueous phase. When a prescribed $\text{ CO}_{2}$ mass flow rate is injected, $\text{ CO}_{2}$ advances initially through the top portion of the aquifer. As $\text{ CO}_{2}$ is being injected, the $\text{ CO}_{2}$ plume advances not only laterally, but also vertically downwards. However, the $\text{ CO}_{2}$ plume does not necessarily occupy the whole thickness of the aquifer. We found that even in the cases in which the $\text{ CO}_{2}$ plume reaches the bottom of the aquifer, most of the injected $\text{ CO}_{2}$ enters the aquifer through the layers at the top. Both $\text{ CO}_{2}$ plume position and fluid pressure compare well with numerical simulations. This solution permits quick evaluations of the $\text{ CO}_{2}$ plume position and fluid pressure distribution when injecting supercritical $\text{ CO}_{2}$ in a deep saline aquifer.     

Random Data Cauchy Theory for the Generalized Incompressible Navier–Stokes Equations

In this paper, we consider the generalized Navier?CStokes equations where the space domain is ${\mathbb{T}^N}$ or ${\mathbb{R}^N, N\geq3}$ . The generalized Navier?CStokes equations here refer to the equations obtained by replacing the Laplacian in the classical Navier?CStokes equations by the more general operator (???) ^?? with ${\alpha\in (\frac{1}{2},\frac{N+2}{4})}$ . After a suitable randomization, we obtain the existence and uniqueness of the local mild solution for a large set of the initial data in ${H^s, s\in[-\alpha,0]}$ , if ${1 < \alpha < \frac{N+2}{4}, s\in(1-2\alpha,0]}$ , if ${\frac{1}{2} < \alpha\leq 1}$ . Furthermore, we obtain the probability for the global existence and uniqueness of the solution. Specially, our result shows that, in some sense, the Cauchy problem of the classical Navier?CStokes equation is local well-posed for a large set of the initial data in H ^?1+, exhibiting a gain of ${\frac{N}{2}-}$ derivatives with respect to the critical Hilbert space ${H^{\frac{N}{2}-1}}$ .     

Well-Posedness of Nematic Liquid Crystal Flow in {L^{3}_{\rm uloc}(\mathbb{R}^3)}

In this paper, we establish the local well-posedness for the Cauchy problem of a simplified version of hydrodynamic flow of nematic liquid crystals in ${\mathbb{R}^3}$ for any initial data (u ₀, d ₀) having small ${L^{3}_{\rm uloc}}$ -norm of ${(u_{0}, \nabla d_{0})}$ . Here ${L^{3}_{\rm uloc}(\mathbb{R}^3)}$ is the space of uniformly locally L ³-integrable functions. For any initial data (u ₀, d ₀) with small ${\|(u_0, \nabla d_0)\|_{L^{3}(\mathbb{R}^3)}}$ , we show that there exists a unique, global solution to the problem under consideration which is smooth for t > 0 and has monotone deceasing L ³-energy for ${t \geqq 0}$ .     

Spatiotemporal complexity of a three-species ratio-dependent food chain model

In this paper, we investigate the complex dynamics of a ratio-dependent spatially extended food chain model. Through a detailed analytical study of the reaction–diffusion model, we obtain some conditions for global stability. On the basis of bifurcation analysis, we present the evolutionary process of pattern formation near the coexistence equilibrium point $(N^*,P^*,Z^*)$ via numerical simulation. And the sequence cold spots $\rightarrow $ stripe–spots mixtures $\rightarrow $ stripes $\rightarrow $ hot stripe–spots mixtures $\rightarrow $ hot spots $\rightarrow $ chaotic wave patterns controlled by parameters $a_1$ or $c_1$ in the model are presented. These results indicate that the reaction–diffusion model is an appropriate tool for investigating fundamental mechanism of complex spatiotemporal dynamics.     

Bounded and Unbounded Motions in Asymmetric Oscillators at Resonance

We consider the boundedness and unboundedness of solutions for the asymmetric oscillator $$\begin{aligned} x''+ax^+-bx^-+g(x)=p(t), \end{aligned}$$ where $x^+=\max \{x,0\},x^-=\max \{-x,0\}, a$ and $b$ are two positive constants, $ p(t)$ is a $2\pi $ -periodic smooth function and $g(x)$ satisfies $\lim _{|x|\rightarrow +\infty }x^{-1}g(x)=0$ . We establish some sharp sufficient conditions concerning the boundedness of all the solutions and the existence of unbounded solutions. It turns out that the boundedness of all the solutions and the existence of unbounded solutions have a close relation to the interaction of some well-defined functions $\Phi _p(\theta )$ and $\Lambda (h)$ . Some explicit conditions are given for the boundedness of all the solutions and the existence of unbounded solutions. Unlike many existing results in the literature where the function $g(x)$ is required to be a bounded function with asymptotic limits, here we allow $g(x)$ be unbounded or oscillatory without asymptotic limits.     

Hadamard Variational Formula for the Green’s Function of the Boundary Value Problem on the Stokes Equations

For every ${\varepsilon > 0}$ , we consider the Green’s matrix ${G_{\varepsilon}(x, y)}$ of the Stokes equations describing the motion of incompressible fluids in a bounded domain ${\Omega_{\varepsilon} \subset \mathbb{R}^d}$ , which is a family of perturbation of domains from ${\Omega\equiv \Omega_0}$ with the smooth boundary ${\partial\Omega}$ . Assuming the volume preserving property, that is, ${\mbox{vol.}\Omega_{\varepsilon} = \mbox{vol.}\Omega}$ for all ${\varepsilon > 0}$ , we give an explicit representation formula for ${\delta G(x, y) \equiv \lim_{\varepsilon\to +0}\varepsilon^{-1}(G_{\varepsilon}(x, y) - G_0(x, y))}$ in terms of the boundary integral on ${\partial \Omega}$ of ${G_0(x, y)}$ . Our result may be regarded as a classical Hadamard variational formula for the Green’s functions of the elliptic boundary value problems.     

An Updated Portrait of Transition to Turbulence in Laminar Pipe Flows with Periodic Time Dependence (A Correlation Study)

Transition to turbulence in axially symmetrical laminar pipe flows with periodic time dependence classified as pure oscillating and pulsatile (pulsating) ones is the concern of the paper. The current state of art on the transitional characteristics of pulsatile and oscillating pipe flows is introduced with a particular attention to the utilized terminology and methodology. Transition from laminar to turbulent regime is usually described by the presence of the disturbed flow with small amplitude perturbations followed by the growth of turbulent bursts. The visual treatment of velocity waveforms is therefore a preferred inspection method. The observation of turbulent bursts first in the decelerating phase and covering the whole cycle of oscillation are used to define the critical states of the start and end of transition, respectively. A correlation study referring to the available experimental data of the literature particularly at the start of transition are presented in terms of the governing periodic flow parameters. In this respect critical oscillating and time averaged Reynolds numbers at the start of transition; Re _os,crit and Re _ta,crit are expressed as a major function of Womersley number, $\sqrt {\omega ^\prime } $ defined as dimensionless frequency of oscillation, f. The correlation study indicates that in oscillating flows, an increase in Re _os,crit with increasing magnitudes of $\sqrt {\omega ^\prime } $ is observed in the covered range of $1<\sqrt {\omega ^\prime } <72$ . The proposed equation (Eq. 7), ${\rm{Re}}_{os,crit} ={\rm{Re}}_{os,crit} \left( {\sqrt {\omega ^\prime } } \right)$ , can be utilized to estimate the critical magnitude of $\sqrt {\omega ^\prime }$ at the start of transition with an accuracy of ±12?% in the range of $\sqrt {\omega ^\prime } <41$ . However in pulsatile flows, the influence of $\sqrt {\omega ^\prime }$ on Re _ta,crit seems to be different in the ranges of $\sqrt {\omega ^\prime } <8$ and $\sqrt {\omega ^\prime } >8$ . Furthermore there is rather insufficient experimental data in pulsatile flows considering interactive influences of $\sqrt {\omega ^\prime } $ and velocity amplitude ratio, A ₁. For the purpose, the measurements conducted at the start of transition of a laminar sinusoidal pulsatile pipe flow test case covering the range of 0.21<?A ₁?<0.95 with $\sqrt {\omega ^\prime } <8$ are evaluated. In conformity with the literature, the start of transition corresponds to the observation of first turbulent bursts in the decelerating phase of oscillation. The measured data indicate that increase in $\sqrt {\omega ^\prime } $ is associated with an increase in Re _ta,crit up to $\sqrt {\omega ^\prime } =3.85$ while a decrease in Re _ta,crit is observed with an increase in $\sqrt {\omega ^\prime } $ for $\sqrt {{\omega }'} >3.85$ . Eventually updated portrait is pointing out the need for further measurements on i) the end of transition both in oscillating and pulsatile flows with the ranges of $\sqrt {\omega ^\prime } <8$ and $\sqrt {\omega ^\prime } >8$ , and ii) the interactive influences of $\sqrt {\omega ^\prime } $ and A ₁ on Re _ta,crit in pulsatile flows with the range of $\sqrt {\omega ^\prime } >8$ .     

Global Regularity for the Two-Dimensional Anisotropic Boussinesq Equations with Vertical Dissipation

This paper establishes the global in time existence of classical solutions to the two-dimensional anisotropic Boussinesq equations with vertical dissipation. When only vertical dissipation is present, there is no direct control on the horizontal derivatives and the global regularity problem is very challenging. To solve this problem, we bound the derivatives in terms of the ${L^\infty}$ -norm of the vertical velocity v and prove that ${\|v\|_{L^{r}}}$ with ${2\leqq r < \infty}$ does not grow faster than ${\sqrt{r \log r}}$ at any time as r increases. A delicate interpolation inequality connecting ${\|v\|_{L^\infty}}$ and ${\|v\|_{L^r}}$ then yields the desired global regularity.