A nonlinear stability analysis for rotating magnetized ferrofluid heated from below saturating a porous medium

A nonlinear (energy) stability analysis is performed for a rotating magnetized ferrofluid layer heated from below saturating a porous medium, in the stress-free boundary case. By introducing a generalized energy functional, a rigorous nonlinear stability result for a thermoconvective rotating magnetized ferrofluid is derived. The mathematical emphasis is on how to control the nonlinear terms caused by magnetic body force. It is found that the nonlinear critical stability magnetic thermal Rayleigh number does not coincide with that of linear instability analysis, and thus indicates that the subcritical instabilities are possible. However, it is noted that, in case of non-ferrofluid, global nonlinear stability Rayleigh number is exactly the same as that for linear instability. For lower values of magnetic parameters, this coincidence is immediately lost. The effect of magnetic parameter, M ₃, medium permeability, D _a , and rotation, , on subcritical instability region has also been analyzed. It is shown that with the increase of magnetic parameter, M ₃, and Darcy number, D _a , the subcritical instability region between the two theories decreases quickly while with the increase of Taylor number, , the subcritical region expands. We also demonstrate coupling between the buoyancy and magnetic forces in the presence of rotation in nonlinear energy stability analysis as well as in linear instability analysis.      

A nonlinear stability analysis of a rotating double-diffusive magnetized ferrofluid

A nonlinear stability result for a double-diffusive magnetized ferrofluid layer rotating about a vertical axis for stress-free boundaries is derived via generalized energy method. The mathematical emphasis is on how to control the nonlinear terms caused by magnetic body and inertia forces. The result is compared with the result obtained by linear instability theory. The critical magnetic thermal Rayleigh number given by energy theory is slightly less than those given by linear theory and thus indicates the existence of subcritical instability for ferrofluids. For non-ferrofluids, it is observed that the nonlinear critical stability thermal Rayleigh number coincides with that of linear critical stability thermal Rayleigh number. For lower values of magnetic parameters, this coincidence is immediately lost. The effect of magnetic parameter, M₃, solute gradient, S₁, and Taylor number, T_A1, on subcritical instability region have been analyzed. We also demonstrate coupling between the buoyancy and magnetic forces in the nonlinear stability analysis.     

Effect of dust particles on a layer of micropolar ferromagnetic fluid heated from below saturating a porous medium

The problem of the effect of dust particles on the thermal convection in micropolar ferromagnetic fluid saturating a porous medium subject to a transverse uniform magnetic field has been investigated theoretically. Linear stability analysis and normal mode analysis methods are used to find an exact solution for a flat micropolar ferromagnetic fluid layer contained between two free boundaries. In case of stationary convection, the effect of various parameters like medium permeability, dust particles, non-buoyancy magnetization, coupling parameter, spin-diffusion parameter and micropolar heat conduction parameter are analyzed. For sufficiently large values of magnetic parameter M₁, the critical magnetic thermal Rayleigh number for the onset of instability is determined numerically and results are depicted graphically. It is also observed that the critical magnetic thermal Rayleigh number is reduced solely because the heat capacity of clean fluid is supplemented by that of the dust particles. The principle of exchange of stabilities is found to hold true for the micropolar ferromagnetic fluid saturating a porous medium heated from below in the absence of micropolar viscous effect, microinertia and dust particles.     

Electrohydrodynamic instability of a rotating layer of a viscoelastic fluid heated from below

The stability of a rotating layer of viscoelastic dielectric liquid (Walters liquid B) heated from below is considered. Linear stability theory is used to derive an eigenvalue system of ten orders and exact eigenvalue equation for a neutral instability is obtained. Under somewhat artificial boundary conditions, this equation can be solved exactly to yield the required eigenvalue relationship from which various critical values are determined in detail. Critical Rayleigh heat numbers and wavenumber for the onset of instability are presented graphically as function of the Taylor number for various values of electric Rayleigh number and the elastic parameters.     

Stability of plane parallel shear flow in a rotating layer heated from below

The discussion of stability of plane parallel shear flow in an infinite rotating layer heated from below requires a mathematical analysis of this problem in dependence on four parameters. These are the Reynolds- and Rayleigh-number, controlling the strength of the shear flow and the heating power, respectively, the Prandtl-number, which measures the relative influence of viscosity and thermal conductivity, and the rotation rate of the layer. After discussing some physical background, possible applications and laboratory experiments two major problems are addressed: i) To find out the cases where unconditional (global) stability up to criticality takes place. In these situations theory makes the clearest predictions and coincidence between experiments and mathematical theory can be expected. ii) To prove that the (monotonic) energy-stability limit is assumed by 2-dimensional (with respect to the spatial variables) perturbations. The solution of this variational problem shows that in certain situations the critical perturbations are 2-dimensional. In these situations, at least, the stability problem is completely solved. Received May 21, 1996 / Accepted December 17, 1996     

Bifurcation analysis of a viscoelastic fluid heated from below

Convection of a viscoelastic fluid in a square domain heated from below is investigated for the case of nondeformable free surfaces. To describe the rheological behavior of the fluid the generalized Oldroyd model is used. A weakly nonlinear analysis is performed in order to determine the character of branching for both the monotonic and oscillatory modes. We also perform a reduction of the boundary value-problem to the set of nonlinear amplitude equations. The analysis of this dynamic system demonstrates the onset and competition of five convection modes.     

Nonlinear rotating convection in a sparsely packed porous medium

We investigate linear and weakly nonlinear properties of rotating convection in a sparsely packed Porous medium. We obtain the values of Takens–Bogdanov bifurcation points and co-dimension two bifurcation points by plotting graphs of neutral curves corresponding to stationary and oscillatory convection for different values of physical parameters relevant to rotating convection in a sparsely packed porous medium near a supercritical pitchfork bifurcation. We derive a nonlinear two-dimensional Landau–Ginzburg equation with real coefficients by using Newell–Whitehead method [16]. We investigate the effect of parameter values on the stability mode and show the occurrence of secondary instabilities viz., Eckhaus and Zigzag Instabilities. We study Nusselt number contribution at the onset of stationary convection. We derive two nonlinear one-dimensional coupled Landau–Ginzburg type equations with complex coefficients near the onset of oscillatory convection at a supercritical Hopf bifurcation and discuss the stability regions of standing and travelling waves.     

A numerical investigation of the nonlinear wave stability of vertical thermal boundary layer flow in a porous medium

A two-dimensional, nonlinear, time-dependent, elliptic, numerical method coupled with an appropriate coordinate transformation is used to investigate the stability of free convection induced by an isothermally heated semi-infinite surface embedded in a fluid-saturated porous medium. It is found that the basic boundary layer flow is stable even to large amplitude disturbance for nondimensional distances of up to 1024 from the leading edge of the heated surface.     

Linear,hydrodynamic flow in a rotating saturated porous medium

Linear, steady, axisymmetric flow of a homogeneous fluid in a rigid, bounded, rotating, saturated porous medium is analyzed. The fluid motions are driven by differential rotation of horizontal boundaries. The dynamics of the interior region and vertical boundary layers are investigated as functions of the Ekman number E(=v/ΩL ²) and rotational Darcy 3 numberN(=kΩ/v) which measures the ratio between the Coriolis force and the Darcy frictional term. IfN≫E ^−1/2, the permeability is sufficiently high and the flow dynamics are the same as those of the conventional free flow problem with Stewartson'sE ^1/3 andE ^1/4 double layer structure. For values ofN≤E ^−1/2 the effect of porous medium is felt by the flow; the Taylor-Proudman constraint is no longer valid. ForN≪E ^−1/3 the porous medium strongly affects the flow; viscous side wall layer is absent to the lowest order and the fluid pumped by the Ekman layer, returns through a region of thicknessO(N ⁻¹). The intermediate rangeE ^−1/3≪N≪E ^−1/2 is characterized by double side wall layer structure: (1)E ^1/3 layer to return the mass flux and (ii) (NE)^1/2 layer to adjust the interior azimuthal velocity to that of the side wall. Spin-up problem is also discussed and it is shown that the steady state is reached quickly in a time scaleO(N).     

Numerical solution of a heated subsidence mound problem in a porous medium

Two different numerical models are constructed to solve a two dimensional subsidence mound problem heated along the moving wet/dry interface. One numerical model is based on cartesian coordinates while the other is based on polar coordinates. In both approaches coordinate transformations are used that render the interface stationary. The problem involves a system of three coupled equations; an elliptic equation for a stream function, a parabolic equation for the temperature and a non-linear equation for the boundary location. Good agreement is found between the results of both methods. Graphic results are presented for the decay of a subsidence mound for different values of the various parameters in the model problem.     

Stability of a layer of dipolar fluid heated from below

The equations of Bleustein and Green [2] are formulated in a way suitable to describe the convective instability which occurs when a layer of dipolar fluid is heated from below. The linear instability boundary is shown to coincide with the nonlinear stability curve and the critical Rayleigh numbers describing this boundary are found; in particular, the non-dimensional micro-length is found to always stabilize.     

A remark on the porous medium equation with nonlinear source

In the present paper, we obtain a new a priori estimate of the solution of the initial-boundary value problem for the porous medium equation with nonlinear source and formulate the conditions guaranteeing the global solvability of this problem.     

A nonlinear two phase fluid flow through a porous medium in presence of a well

We study a flow of fresh and salt water in a two dimensional axially symmetric coastal aquifer with a well on the central axis. The flow is governed by a nonlinear Darcy's law. We also show the behaviour of the solution when the out flow of salt water at well goes to 0. Received May 1999     

On critical Fujita exponents for the porous medium equation with a nonlinear boundary condition

We establish the critical Fujita exponents for the solution of the porous medium equation u_t=Δu^m, x∈R₊^N, t>0, subject to the nonlinear boundary condition −∂u^m/∂x₁=u^p, x₁=0, t>0, in multi-dimension.     

Separation of variables of a generalized porous medium equation with nonlinear source

This paper considers a general form of the porous medium equation with nonlinear source term: u_t=(D(u)u_xⁿ)_x+F(u), n≠1. The functional separation of variables of this equation is studied by using the generalized conditional symmetry approach. We obtain a complete list of canonical forms for such equations which admit the functional separable solutions. As a consequence, some exact solutions to the resulting equations are constructed, and their behavior are also investigated.     

Existence of weak solutions for fractional porous medium equations with nonlinear term

We study the following fractional porous medium equations with nonlinear term $\{\begin{array}{cc}{u}_t+{(-\Delta )}^{\sigma /2}\left({\left|u\right|}^{m-1}u\right)+g\left(u\right)=h\text{,}& \text{in}\Omega \times {\mathbb{R}}^{+}\text{,}\\ u(x,t)=0\text{,}& \text{in}\partial \Omega \times {\mathbb{R}}^{+}\text{,}\\ u(x,0)=u_0\text{,}& \text{in}\Omega \text{.}\end{array}$ The authors in de Pablo et al. (2011) and de Pablo et al. (2012) established the existence of weak solutions for the case $g\left(u\right)\equiv 0$. Here, we consider the nonlinear term $g$ is without an upper growth restriction. The nonlinearity of $g$ leads to the invalidity of the Crandall–Liggett theorem, which is the critical method to establish the weak solutions in de Pablo et al. (2011) and de Pablo et al. (2012). In addition, because of $g$ does not have an upper growth restriction, we have to apply the weak compactness theorem in an Orlicz space to prove the existence of weak solutions by using the Implicit Time Discretization method.     

Symmetrization for linear and nonlinear fractional parabolic equations of porous medium type

We establish symmetrization results for the solutions of the linear fractional diffusion equation _∂tu+(−Δ)^σ/2u=f${\partial }_{t}u+{(-\mathrm{\Delta })}^{\sigma /2}u=f$ and its elliptic counterpart hv+(−Δ)^σ/2v=f$hv+{(-\mathrm{\Delta })}^{\sigma /2}v=f$, h>0$h>0$, using the concept of comparison of concentrations. The results extend to the nonlinear version, _∂tu+(−Δ)^σ/2A(u)=f${\partial }_{t}u+{(-\mathrm{\Delta })}^{\sigma /2}A(u)=f$, but only when the nondecreasing function A:_R+→_R+$A:{\mathbb{R}}_{+}\to {\mathbb{R}}_{+}$ is concave. In the elliptic case, complete symmetrization results are proved for B(v)+(−Δ)^σ/2v=f$B(v)+{(-\mathrm{\Delta })}^{\sigma /2}v=f$ when B(v)$B(v)$ is a convex nonnegative function for v>0$v>0$ with B(0)=0$B(0)=0$, and partial results hold when B is concave. Remarkable counterexamples are constructed for the parabolic equation when A is convex, resp. for the elliptic equation when B   is concave. Such counterexamples do not exist in the standard diffusion case σ=2$\sigma =2$.     

Helical flows of a heated generalized Oldroyd-B fluid subject to a time-dependent shear stress in porous medium

This paper presents an analysis for helical flows of a heated generalized Oldroyd-B fluid subject to a time-dependent shear stress in porous medium, where the motion is due to the longitudinal time-dependent shear stress and the oscillating velocity in boundary. The exact solutions are established by using the sequential fractional derivatives Laplace transform coupled with finite Hankel transforms in terms of generalized G function. Moreover, the effects of various parameters (relaxation time, fractional parameter, permeability and porosity) on the flow and heat transfer are analyzed in detail by graphical illustrations.