Effect of rotation on the onset of double diffusive convection in a Darcy porous medium saturated with a couple stress fluid

The effect of rotation on the onset of double diffusive convection in a horizontal couple stress fluid-saturated porous layer, which is heated and salted from below, is studied analytically using both linear and weak nonlinear stability analyses. The extended Darcy model, which includes the time derivative and Coriolis terms, has been employed in the momentum equation. The onset criterion for stationary, oscillatory and finite amplitude convection is derived analytically. The effect of Taylor number, couple stress parameter, solute Rayleigh number, Lewis number, Darcy–Prandtl number, and normalized porosity on the stationary, oscillatory, and finite amplitude convection is shown graphically. It is found that the rotation, couple stress parameter and solute Rayleigh number have stabilizing effect on the stationary, oscillatory, and finite amplitude convection. The Lewis number has a stabilizing effect in the case of stationary and finite amplitude modes, with a destabilizing effect in the case of oscillatory convection. The Darcy–Prandtl number and normalized porosity advances the onset of oscillatory convection. A weak nonlinear theory based on the truncated representation of Fourier series method is used to find the finite amplitude Rayleigh number and heat and mass transfer. The transient behavior of the Nusselt number and Sherwood number is investigated by solving the finite amplitude equations using Runge–Kutta method.     

On the linear stability analysis of a Maxwell fluid with double-diffusive convection

The problem of double-diffusive convection and cross-diffusion in a Maxwell fluid in a horizontal layer in porous media is re-examined using the modified Darcy–Brinkman model. The effect of Dufour and Soret parameters on the critical Darcy–Rayleigh numbers is investigated. Analytical expressions of the critical Darcy–Rayleigh numbers for the onset of stationary and oscillatory convection are derived. Numerical simulations show that the presence of Dufour and Soret parameters has a significant effect on the critical Darcy–Rayleigh number for over-stability. In the limiting case some previously published results are recovered.     

An analytical study of linear and nonlinear double diffusive convection in a fluid saturated anisotropic porous layer with Soret effect

The double diffusive convection in a horizontal anisotropic porous layer saturated with a Boussinesq fluid, which is heated and salted from below in the presence of Soret coefficient is studied analytically using both linear and nonlinear stability analyses. The normal mode technique is used in the linear stability analysis while a weak nonlinear analysis based on a minimal representation of double Fourier series method is used in the nonlinear analysis. The generalized Darcy model including the time derivative term is employed for the momentum equation. The critical Rayleigh number, wavenumber for stationary and oscillatory modes and frequency of oscillations are obtained analytically using linear theory. The effect of anisotropy parameters, solute Rayleigh number, Soret parameter and Lewis number on the stationary, oscillatory, finite amplitude convection and heat and mass transfer are shown graphically.     

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.     

Porous convection and thermal oscillations

A theory is developed for thermal convection in a fluid saturated porous material when the temperature may propagate as a thermal wave. In particular, we are interested in the mechanism of thermal oscillation and so allow for Guyer?CKrumhansl effects but employ a heat flux equation developed by Christov and Morro. The instability mechanism is investigated in complete detail and it is shown that stationary convection is likely to prevail under normal terrestrial conditions, but if the thermal relaxation time is sufficiently large there is a possible parameter range which allows for oscillatory convection. However, the presence of the Guyer?CKrumhansl terms has the effect of damping the oscillatory convection and returning the instability mechanism to one of stationary convection.     

Oscillatory convection and the Cattaneo law of heat conduction

The problem of thermal convection is investigated for a layer of fluid when the heat flux law of Cattaneo is adopted. The boundary conditions are those appropriate to two fixed surfaces. It is shown that for small Cattaneo number the critical Rayleigh number initially increases from its classical value of 1707.765 until a critical value of the Cattaneo number is reached. For Cattaneo numbers greater than this critical value a notable Hopf bifurcation is observed with convection occurring at lower Rayleigh numbers and by oscillatory rather than stationary convection. The aspect ratio of the convection cells likewise changes.     

Effect of magnetic field dependent viscosity on ferroconvection in the presence of dust particles

This paper deals with the theoretical investigation of the effect of magnetic field dependent (MFD) viscosity on the thermal convection in a ferromagnetic fluid in the presence of dust particles. For a flat ferromagnetic fluid layer contained between two free boundaries, the exact solution is obtained using a linear stability analysis and a normal mode analysis method. For the case of stationary convection, dust particles always have a destabilizing effect, whereas the MFD viscosity has a stabilizing effect on the onset of convection. In the absence of MFD viscosity, the destabilizing effect of magnetization is depicted but in the presence of MFD viscosity, non-buoyancy magnetization may have a destabilizing or a stabilizing effect on the onset of convection. The critical wave number and critical magnetic thermal Rayleigh number for the onset of stationary convection are also determined numerically for sufficiently large values of buoyancy magnetization parameter M ₁. Graphs have been plotted by giving numerical values to the parameters to depict the stability characteristics. It is 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 ferromagnetic fluid heated from below in the absence of dust particles. The oscillatory modes are introduced due to the presence of the dust particles, which were non-existent in their absence. A sufficient condition for the non-existence of overstability is also obtained.     

Magnetohydrodynamic stationary and oscillatory convective stability in a mushy layer during binary alloy solidification

For the case of solidification of a bottom cooled binary alloy, the magnetohydrodynamic stationary and oscillatory convective stability in the mushy layer is investigated analytically using normal mode linear stability analysis. In the limit of large Stefan number (S_t), a near–eutectic approximation with large far field temperature is considered in the present research. To ascertain the instability in the mushy layer, the strength of the superimposed magnetic field is so chosen that it corresponds to a given mush Hartmann number (Ha_m) of the problem. The results are presented for various values of mush Hartmann numbers in the range, 0 ≤ Ha_m ≤ 50. The critical Rayleigh number for stationary convection shows a linear relationship with increasing Ha_m. The magnetohydrodynamic effect imparts a stabilizing influence during stationary convection. In comparison to that of the stationary convective mode, the oscillatory mode appears to be critically susceptible at higher values of β (β = S_t/℘² ϒ², ℘ is the compositional ratio, ϒ = 1 + S_t/℘), and vice versa for lower β values. Analogous to the behavior for stationary convection, the magnetic field also offers a stabilizing effect in oscillatory convection and thus influences global stability of the mushy layer. Increasing magnetic strength shows reduction in the wavenumber and in the number of rolls formed in the mushy layer.     

Direct Resonance in Double-Diffusive Systems

In thermohaline convection and a number of other systems, both direct and oscillatory modes of instability are possible. Should these modes coalesce at some point in parameter space, nonlinear instabilities are likely to be more dramatic in the neighborhood of such a point than elsewhere. A finite-amplitude evolution equation describing such events is derived here by the method of multiple scales. The results are compared with those obtained previously for model systems and by other methods. The direct resonance point of view is found to provide some new insights. A generalization to allow for slow spatial modulation of the amplitude is given and, among other possibilities, solitary-wave solutions are obtainable.     

Stationary statistical properties of Rayleigh‐Bénard convection at large Prandtl number

This is the third in a series of our study of Rayleigh‐Bénard convection at large Prandtl number. Here we investigate whether stationary statistical properties of the Boussinesq system for Rayleigh‐Bénard convection at large Prandtl number are related to those of the infinite Prandtl number model for convection that is formally derived from the Boussinesq system via setting the Prandtl number to infinity. We study asymptotic behavior of stationary statistical solutions, or invariant measures, to the Boussinesq system for Rayleigh‐Bénard convection at large Prandtl number. In particular, we show that the invariant measures of the Boussinesq system for Rayleigh‐Bénard convection converge to those of the infinite Prandtl number model for convection as the Prandtl number approaches infinity. We also show that the Nusselt number for the Boussinesq system (a specific statistical property of the system) is asymptotically bounded by the Nusselt number of the infinite Prandtl number model for convection at large Prandtl number. We discover that the Nusselt numbers are saturated by ergodic invariant measures. Moreover, we derive a new upper bound on the Nusselt number for the Boussinesq system at large Prandtl number of the form which asymptotically agrees with the (optimal) upper bound on Nusselt number for the infinite Prandtl number model for convection. © 2007 Wiley Periodicals, Inc.     

Thermal convection problem of micropolar fluid subjected to hall current

Thermal instability of a micropolar fluid layer heated from below in the presence of hall currents is investigated. Using the appropriate boundary conditions on the boundary surfaces of the fluid layer, the frequency equation is derived and then critical Rayleigh number is determined. It is found that hall current parameter has destabilizing effect on the system. For specific values of parameters, oscillatory convection in observed in the system. The behavior of Rayleigh number with wavenumber is also computed for different values of various parameters. The results of some earlier workers have been reduced as a special case from the present problem.     

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.     

On optimization of the RBF shape parameter in a grid-free local scheme for convection dominated problems over non-uniform centers

Global optimization techniques exist in the literature for finding the optimal shape parameter of the infinitely smooth radial basis functions (RBF) if they are used to solve partial differential equations. However these global collocation methods, applied directly, suffer from severe ill-conditioning when the number of centers is large. To circumvent this, we have used a local optimization algorithm, in the optimization of the RBF shape parameter which is then used to develop a grid-free local (LRBF) scheme for solving convection–diffusion equations. The developed algorithm is based on the re-construction of the forcing term of the governing partial differential equation over the centers in a local support domain. The variable (optimal) shape parameter in this process is obtained by minimizing the local Cost function at each center (node) of the computational domain. It has been observed that for convection dominated problems, the local optimization scheme over uniform centers has produced oscillatory solutions, therefore, in this work the local optimization algorithm has been experimented over Chebyshev and non-uniform distribution of the centers. The numerical experiments presented in this work have shown that the LRBF scheme with the local optimization produced accurate and stable solutions over the non-uniform points even for convection dominant convection–diffusion equations.     

Linear and nonlinear stability analysis of binary viscoelastic fluid convection

The linear and weakly nonlinear stability analysis of the quiescent state in a viscoelastic fluid subject to vertical solute concentration and temperature gradients is investigated. The non-Newtonian behavior of the viscoelastic fluid is characterized using the Oldroyd model. Analytical expressions for the critical Rayleigh numbers and corresponding wave numbers for the onset of stationary or oscillatory convection subject to cross diffusion effects is determined. A stability diagram clearly demarcates non-overlapping regions of finger and diffusive instabilities. A Lorenz system is obtained in the case of the weakly nonlinear stability analysis. The effect of Dufour and Soret parameters on the heat and mass transports are determined and discussed. Due to consideration of dilute concentrations of the second diffusing component the route to chaos in binary viscoelastic fluid systems is similar to that of single-component (thermal) viscoelastic fluid systems.     

Long-wave oscillatory Marangoni instability in a horizontal liquid layer

The long-wave instability in the problem of thermocapillary convection in a horizontal layer with a free deformable boundary and a solid bottom is investigated. The transcendental equation for the main asymptotic term of the spectral parameter is written in explicit form. The main attention is paid to investigating oscillatory instability. For the frequency of neutral oscillations, simple transcendental equations are obtained that contain the Prandtl and Biot numbers. In a number of cases, exact solutions are indicated. Explicit formulae are given for the main asymptotic term of the Marangoni number. In the case of a non-heat-conducting solid wall, the relation between the critical values of the parameters for inverse Prandtl numbers is found. It is shown that, for different Prandtl numbers, the asymptotic values are in good agreement with the numerical values.     

On the stability of mixed convective motions in a vertical layer with wave-like boundaries

The stability of the stationary and oscillatory convective motions which develop in a vertical layer with periodically curved boundaries is studied for the case of longitudinal fluid injection. The amplitude of the boundary undulations and the flow of fluid along the layer are both assumed to be small, and methods of perturbation theory are used. The characteristic properties of the incremental spectrum of the spatially periodic motions are studied and the most dangerous types of perturbations as well as the forms of the stability regions are determined.

Theoretical investigations of the effect of spatial inhomogeneity of the boundary conditions on the stability of convection were sparse, and they deal mainly with horizontal layers of fluid /1–3/. Stationary, spatially periodic motions in a vertical layer with curved boundaries were investigated in /4/ for the case of free convection (when the flow was closed), and their stability was investigated in /5/. It was established that the presence of a small but finite flow of fluid along the layer leads to an increase in the number of different modes of flow, and to the appearance of non-stationary convective motions in the region near the threshold.     

Turing pattern outside of the Turing domain

There are two simple solutions to reaction-diffusion systems with limit-cycle reaction kinetics, producing oscillatory behaviour. The reaction parameter mu gives rise to a 'space-invariant' solution, and mu versus the ratio of the diffusion coefficients gives rise to a 'time-invariant' solution. We consider the case where both solution types may be possible. This leads to a refinement of the Turing model of pattern formation. We add convection to the system and investigate its effect. More complex solutions arise that appear to combine the two simple solutions. The convective system sheds light on the underlying behaviour of the diffusive system.     

Nontrivial stationary points of two-species self-structured communities

The two-species model of self-structured stationary biological communities proposed by U. Dieckmann and R. Law is considered. A way of investigating the system of integro-differential equations describing the model equilibrium is developed, nontrivial stationary points are found, and constraints on the model parameter space resulting in similar stationary points are studied. The results are applied to a number of widely known biological scenarios.     

Continuous Dependence on the Heat Source in Resonant Porous Penetrative Convection

I study the structural stability for a problem in a porous medium when the density of saturating liquid is a nonlinear function of temperature and an internal heat source is present. It has been shown that for this problem when one considers thermal convection in a plane infinite layer then resonance may occur between internal layers that arise. A key parameter is the internal heat source and its presence may lead to oscillatory instability inducing resonance. Therefore, in this paper, I analyze the general structural stability problem of continuous dependence on the heat source itself for a model of nonisothermal flow in a porous medium of Forchheimer type, in a general three‐dimensional domain.     

Stochastic-periodic homogenization of Maxwell’s equations with linear and periodic conductivity

Stochastic-periodic homogenization is studied for the Maxwell equations with the linear and periodic electric conductivity. It is shown by the stochastic-two-scale convergence method that the sequence of solutions to a class of highly oscillatory problems converges to the solution of a homogenized Maxwell equation.