A New Insight into the Convective Boundary Condition

The steady mixed convection boundary layer flows over a vertical surface adjacent to a Darcy porous medium and subject respectively to (i) a prescribed constant wall temperature, (ii) a prescribed variable heat flux, $q_\mathrm{w} =q_0 x^{-1/2}$ q w = q 0 x ? 1 / 2 , and (iii) a convective boundary condition are compared to each other in this article. It is shown that, in the characteristic plane spanned by the dimensionless flow velocity at the wall ${f}^{\prime }(0)\equiv \lambda $ f ′ ( 0 ) ≡ λ and the dimensionless wall shear stress $f^{\prime \prime }(0)\equiv S$ f ′ ′ ( 0 ) ≡ S , every solution $(\lambda , S)$ ( λ , S ) of one of these three flow problems at the same time is also a solution of the other two ones. There also turns out that with respect to the governing mixed convection and surface heat transfer parameters $\varepsilon $ ε and $\gamma $ γ , every solution $(\lambda , S)$ ( λ , S ) of the flow problem (iii) is infinitely degenerate. Specifically, to the very same flow solution $(\lambda , S)$ ( λ , S ) there corresponds a whole continuous set of values of $\varepsilon $ ε and $\gamma $ γ which satisfy the equation $S=-\gamma (1+\varepsilon -\lambda )$ S = ? γ ( 1 + ε ? λ ) . For the temperature solutions, however, the infinite degeneracy of the velocity solutions becomes lifted. These and further outstanding features of the convective problem (iii) are discussed in the article in some detail.     

Linear Instability of a Horizontal Thermal Boundary Layer Formed by Vertical Throughflow in a Porous Medium: The Effect of Local Thermal Nonequilibrium

In this paper we investigate the onset of convection in a saturated porous medium where uniform suction into a horizontal and uniformly hot bounding surface induces a stationary thermal boundary layer. Particular attention is paid to how the well-known linear stability characteristics of this boundary layer are modified by the presence of local thermal nonequilibrium effects. The basic conduction state is determined and it is found that the boundary layer forms two distinct regions when the porosity is small or when the conductivity of the fluid is small compared with that of the solid. A linearised stability analysis is performed which results in an ordinary differential eigenvalue problem for the critical Darcy–Rayleigh number as a function of the wave number and the two nondimensional parameters, $H$ H and $\gamma $ γ , which are associated with local thermal nonequilibrium. This eigenvalue problem is solved numerically by first approximating the equations by fourth order compact finite differences, and then the critical Rayleigh number is computed iteratively using the inverse power method and minimised over the wavenumber. The variation of the critical Rayleigh number and wavenumber with $H$ H and $\gamma $ γ is presented. One of the unusual effects of local thermal nonequilibrium is that there exists a parameter regime within which the neutral curve is bimodal.     

Mixed Convection Boundary-Layer Flow Over a Vertical Surface Embedded in a Porous Material Subject to a Convective Boundary Condition

The mixed convection boundary-layer flow on one face of a semi-infinite vertical surface embedded in a fluid-saturated porous medium is considered when the other face is taken to be in contact with a hot or cooled fluid maintaining that surface at a constant temperature $T_\mathrm{{f}}$ . The governing system of partial differential equations is transformed into a system of ordinary differential equations through an appropriate similarity transformation. These equations are solved numerically in terms of a dimensionless mixed convection parameter $\epsilon $ and a surface heat transfer parameter $\gamma $ . The results indicate that dual solutions exist for opposing flow, $\epsilon <0$ , with the dependence of the critical values $\epsilon _\mathrm{{c}}$ on $\gamma $ being determined, whereas for the assisting flow $\epsilon >0$ , the solution is unique. Limiting asymptotic forms for both $\gamma $ small and large and $\epsilon $ large are also discussed.     

Dynamical Parallelepipeds in Minimal Systems

For a topological dynamical system $(X,T)$ ( X , T ) and $d\in \mathbb N $ d ∈ N , the associated dynamical parallelepiped $\mathbf{Q}^{[d]}$ Q [ d ] was defined by Host–Kra–Maass. For a minimal distal system it was shown by them that the relation $\sim _{d-1}$ ～ d ? 1 defined on $\mathbf{Q}^{[d-1]}$ Q [ d ? 1 ] is an equivalence relation; the closing parallelepiped property holds, and for each $x\in X$ x ∈ X the collection of points in $\mathbf{Q}^{[d]}$ Q [ d ] with first coordinate $x$ x is a minimal subset under the face transformations. We give examples showing that the results do not extend to general minimal systems.     

Conjugate Natural Convection in a Porous Enclosure Sandwiched by Finite Walls Under the Influence of Non-uniform Heat Generation and Radiation

Conjugate natural convection in a square porous enclosure sandwiched by finite walls under the influence of non-uniform heat generation and radiation is studied numerically in the present article. The horizontal heating is considered, where the vertical walls heated isothermally at different temperatures, while the horizontal walls are kept adiabatic. The Darcy model is used in the mathematical formulation for the porous layer and finite difference method is applied to solve the dimensionless governing equations. The governing parameters considered are the ratio of wall thickness to its width $(0.02 \le D \le 0.3)$ ( 0.02 ≤ D ≤ 0.3 ) , the wall to porous thermal conductivity ratio $(0.1 \le k_\mathrm{r} \le 10.0)$ ( 0.1 ≤ k r ≤ 10.0 ) , the internal heating $(0 \le \gamma \le 5)$ ( 0 ≤ γ ≤ 5 ) , and the local heating exponent parameters $(1 \le \lambda \le 20)$ ( 1 ≤ λ ≤ 20 ) . It is found that the average Nusselt number on the hot and cold interfaces increases with increasing the radiation intensity. Very high non-uniformity heating does not affect the average Nusselt number at very thick walls.     

Local Dynamics at Focal Points Coupled with Elliptic Directions

We investigate local dynamics of a 4-dimensional system with small and slowly varying time periodic forcing. By assuming the unperturbed system is autonomous and has a fixed point with eigenvalues $(0,0,i,-i)$ ( 0 , 0 , i , ? i ) , we study homoclinic, subharmonic solutions and Hopf bifurcation in a $O(\epsilon )$ O ( ? ) neighborhood of the fixed point, where $\epsilon $ ? is the perturbation parameter.     

Generalized Resolvent Estimates of the Stokes Equations with First Order Boundary Condition in a General Domain

In this paper, we prove unique existence of solutions to the generalized resolvent problem of the Stokes operator with first order boundary condition in a general domain ${\Omega}$ of the N-dimensional Eulidean space ${\mathbb{R}^N, N \geq 2}$ . This type of problem arises in the mathematical study of the flow of a viscous incompressible one-phase fluid with free surface. Moreover, we prove uniform estimates of solutions with respect to resolvent parameter ${\lambda}$ varying in a sector ${\Sigma_{\sigma, \lambda_0} = \{\lambda \in \mathbb{C} \mid |\arg \lambda| < \pi-\sigma, \enskip |\lambda| \geq \lambda_0\}}$ , where ${0 < \sigma < \pi/2}$ and ${\lambda_0 \geq 1}$ . The essential assumption of this paper is the existence of a unique solution to a suitable weak Dirichlet problem, namely it is assumed the unique existence of solution ${p \in \hat{W}^1_{q, \Gamma}(\Omega)}$ to the variational problem: ${(\nabla p, \nabla \varphi) = (f, \nabla \varphi)}$ for any ${\varphi \in \hat W^1_{q', \Gamma}(\Omega)}$ . Here, ${1 < q < \infty, q' = q/(q-1), \hat W^1_{q, \Gamma}(\Omega)}$ is the closure of ${W^1_{q, \Gamma}(\Omega) = \{ p \in W^1_q(\Omega) \mid p|_\Gamma = 0\}}$ by the semi-norm ${\|\nabla \cdot \|_{L_q(\Omega)}}$ , and ${\Gamma}$ is the boundary of ${\Omega}$ . In fact, we show that the unique solvability of such a Dirichlet problem is necessary for the unique existence of a solution to the resolvent problem with uniform estimate with respect to resolvent parameter varying in ${(\lambda_0, \infty)}$ . Our assumption is satisfied for any ${q \in (1, \infty)}$ by the following domains: whole space, half space, layer, bounded domains, exterior domains, perturbed half space, perturbed layer, but for a general domain, we do not know any result about the unique existence of solutions to the weak Dirichlet problem except for q =  2.     

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.     

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.     

A Comment on the Flow and Heat Transfer Past a Permeable Stretching/Shrinking Surface in a Porous Medium: Brinkman Model

An analytical solution is presented for the boundary-layer flow and heat transfer over a permeable stretching/shrinking surface embedded in a porous medium using the Brinkman model. The problem is seen to be characterized by the Prandtl number $Pr$ , a mass flux parameter $s$ , with $s>0$ for suction, $s=0$ for an impermeable surface, and $s<0$ for blowing, a viscosity ratio parameter $M$ , the porous medium parameter $\Lambda $ and a wall velocity parameter $\lambda $ . The analytical solution identifies critical values which agree with those previously determined numerically (Bachok et al. Proceedings of the fifth International Conference on Applications of Porous Media, 2013) and shows that these critical values, and the consequent dual solutions, can arise only when there is suction through the wall, $s>0$ .     

Invariant Parallels, Invariant Meridians and Limit Cycles of Polynomial Vector Fields on Some 2-Dimensional Algebraic Tori in \mathbb R ^3

We consider the polynomial vector fields of arbitrary degree in $\mathbb R ^3$ R 3 having the 2-dimensional algebraic torus $$\begin{aligned} \mathbb T ^2(l,m,n)=\{(x,y,z)\in \mathbb R ^3 : (x^{2l}+y^{2m}-r^2)^2+z^{2n}-1=0\}, \end{aligned}$$ T 2 ( l , m , n ) = { ( x , y , z ) ∈ R 3 : ( x 2 l + y 2 m - r 2 ) 2 + z 2 n - 1 = 0 } , where $l,m$ l , m , and $n$ n positive integers, and $r\in (1,\infty )$ r ∈ ( 1 , ∞ ) , invariant by their flow. We study the possible configurations of invariant meridians and parallels that these vector fields can exhibit on $\mathbb T ^2(l,m,n)$ T 2 ( l , m , n ) . Furthermore, we analyze when these invariant meridians or parallels are limit cycles.     

Fracture Propagation Driven by Fluid Outflow from a Low-Permeability Aquifer

Deep saline aquifers are promising geological reservoirs for $\mathrm {CO}_2$ CO 2 sequestration if they do not leak. The absence of leakage is provided by the caprock integrity. However, $\mathrm {CO}_2$ CO 2 injection operations may change the geomechanical stresses and cause fracturing of the caprock. We present a model for the propagation of a fracture in the caprock driven by the outflow of fluid from a low-permeability aquifer. We show that to describe the fracture propagation, it is necessary to solve the pressure diffusion problem in the aquifer. We solve the problem numerically for the two-dimensional domain and show that, after a relatively short time, the solution is close to that of one-dimensional problem, which can be solved analytically. We use the relations derived in the hydraulic fracture literature to relate the width of the fracture to its length and the flux into it, which allows us to obtain an analytical expression for the fracture length as a function of time. Using these results we predict the propagation of a hypothetical fracture at the In Salah $\mathrm {CO}_2$ CO 2 injection site to be as fast as a typical hydraulic fracture. We also show that the hydrostatic and geostatic effects cause the increase of the driving force for the fracture propagation and, therefore, our solution serves as an estimate from below. Numerical estimates show that if a fracture appears, it is likely that it will become a pathway for $\mathrm {CO}_2$ CO 2 leakage.     

Non-uniqueness of solutions in asymptotically self-similar shock reflections

The present study addresses the self-similar problem of unsteady shock reflection on an inclined wedge. The start-up conditions are studied by modifying the wedge corner and allowing for a finite radius of curvature. It is found that the type of shock reflection observed far from the corner, namely regular or Mach reflection, depends intimately on the start-up condition, as the flow “remembers” how it was started. Substantial differences were found. For example, the type of shock reflection for an incident shock Mach number $M=6.6$ and an isentropic exponent $\gamma =1.2$ changes from regular to Mach reflection between $44^\circ $ and $45^\circ $ when a straight wedge tip is used, while the transition for an initially curved wedge occurs between $57^\circ $ and $58^\circ $ .     

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.     

The Onset of Prandtl–Darcy–Prats Convection in a Horizontal Porous Layer

We consider the effect of finite Prandtl–Darcy numbers of the onset of convection in a porous layer heated isothermally from below and which is subject to a horizontal pressure gradient. A dispersion relation is found which relates the critical Darcy–Rayleigh number and the induced phase speed of the cells to the wavenumber and the imposed Péclet and Prandtl–Darcy numbers. Exact numerical solutions are given and these are supplemented by asymptotic solutions for both large and small values of the governing nondimensional parameters. The classical value of the critical Darcy–Rayleigh number is $4\pi ^2$ 4 π 2 , and we show that this value increases whenever the Péclet number is nonzero and the Prandtl–Darcy number is finite simultaneously. The corresponding wavenumber is always less than $\pi $ π and the phase speed of the convection cells is always smaller than the background flux velocity.     

The Onset of Convection in a Layered Porous Medium with Vertical Throughflow

An analytical investigation of the effect of vertical throughflow on the onset of convection in a composite porous medium consisting of two horizontal layers has been made. The cases of iso-flux and iso-temperature boundaries are both investigated. The critical Rayleigh number depends on a Péclet number $Q$ , a permeability ratio $K_{r}$ , a thermal conductivity ratio $k_{r}$ , and a depth ratio $\delta $ . For the case of small $Q$ an approximate solution is obtained, which shows that in general throughflow has a stabilizing effect whose magnitude may be increased or decreased by the heterogeneity. This solution is supplemented by an asymptotic solution valid for large $Q.$      

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.     

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^*}$ .