Stability analysis of thermosolutal convection in a horizontal porous layer using a thermal non-equilibrium model

A stability analysis is carried out to investigate the onset of thermosolutal convection in a horizontal porous layer when the solid and fluid phases are not in a local thermal equilibrium, and the solubility of the dissolved component depends on temperature. To study how the reaction and thermal non-equilibrium affect the double-diffusive convection, the effects of scaled inter-phase heat transfer coefficient H and dimensionless reaction rate k on thermosolutal convection are discussed . The critical Rayleigh number and the corresponding wave number for the stability and overstability convections are obtained. Specially, asymptotic analysis for both small and large values of H and k is presented, and the corresponding asymptotic solutions are compared with numerical results. At last, a nonlinear stability analysis is presented to study how H and k affect the Nusselt number.     

Linear and non-linear double diffusive convection in a rotating porous layer using a thermal non-equilibrium model

Double diffusive convection in a fluid-saturated rotating porous layer heated from below and cooled from above is studied when the fluid and solid phases are not in local thermal equilibrium, using both linear and non-linear stability analyses. The Darcy model that includes the time derivative and Coriolis terms is employed as momentum equation. A two-field model that represents the fluid and solid phase temperature fields separately is used for energy equation. The onset criterion for stationary, oscillatory and finite amplitude convection is derived analytically. It is found that small inter-phase heat transfer coefficient has significant effect on the stability of the system. There is a competition between the processes of thermal and solute diffusions that causes the convection to set in through either oscillatory or finite amplitude mode rather than stationary. The effect of solute Rayleigh number, porosity modified conductivity ratio, Lewis number, diffusivity ratio, Vadasz number and Taylor number on the stability of the system is investigated. The non-linear theory based on the truncated representation of Fourier series method predicts the occurrence of subcritical instability in the form of finite amplitude motions. The effect of thermal non-equilibrium on heat and mass transfer is also brought out.     

The effect of local thermal non-equilibrium on conduction in porous channels with a uniform heat source

We examine the effect of local thermal non-equilibrium on the steady state heat conduction in a porous layer in the presence of internal heat generation. A uniform source of heat is present in either the fluid or the solid phase. A two-temperature model is assumed and analytical solutions are presented for the resulting steady-state temperature profiles in a uniform porous slab. Attention is then focussed on deriving simple conditions which guarantee local thermal equilibrium.     

Bounds on convective heat transport in a rotating porous layer

Using the background field variational method, bounds on convective heat transport in a rotating porous layer heated from below are derived from the primitive equations. The enhancement of heat transport beyond the minimal conduction value (the Nusselt number Nu) is bounded in terms of the dimensionless temperature difference across the layer (the Rayleigh number Ra) according to

This rigorous upper bound shows that rotation has a retarding effect on convective heat transport.     

Stability analysis of soret-driven double-diffusive convection of Maxwell fluid in a porous medium

Stability analysis of double-diffusive convection for viscoelastic fluid with Soret effect in a porous medium is investigated using a modified-Maxwell-Darcy model. We use the linear stability analysis to investigate how the Soret parameter and the relaxation time of viscoelastic fluid effect the onset of convection and the selection of an unstable wavenumber. It is found that the Soret effect is to destabilize the system for oscillatory convection. The relaxation time also enhances the instability of the system. The effects of Soret coefficient and relaxation time on the heat transfer rate in a porous medium are studied using the nonlinear stability analysis, the variation of the Nusselt number with respect to the Rayleigh number is derived for stationary and oscillatory convection modes. Some previous results can be reduced as the special cases of the present paper.     

Bifurcation analysis for thermal convection in a rotating porous layer

The linear and weakly nonlinear thermal convection in a rotating porous layer is investigated by constructing a simplified model involving a system of fifth-order nonlinear ordinary differential equations. The flow in the porous medium is described by Lap wood-Brinkman-extended Darcy model with fluid viscosity different from effective viscosity. Conditions for the occurrence of possible bifurcations are obtained. It is established that Hopf bifurcation is possible only at a lower value of the Rayleigh number than that of simple bifurcation. In contrast to the non-rotating case, it is found that the ratio of viscosities as well as the Darcy number plays a dual role on the steady onset and some important observations are made on the stability characteristics of the system. The results obtained from weakly nonlinear theory reveal that, the steady bifurcating solution may be either sub-critical or supercritical depending on the choice of physical parameters. Heat transfer is calculated in terms of Nusselt number.     

Onset of convection in a horizontal porous channel with uniform heat generation using a thermal nonequilibrium model

This paper considers the onset of free convection in a horizontal fluid-saturated porous layer with uniform heat generation. Attention is focused on cases where the fluid and solid phases are not in local thermal equilibrium, and where two energy equations describe the evolution of the temperature of each phase. Standard linearized stability theory is used to determine how the criterion for the onset of convection varies with the inter-phase heat transfer coefficient, H, and the porosity-modified thermal conductivity ratio, γ. We also present asymptotic solutions for small values of H. Excellent agreement is obtained between the asymptotic and the numerical results.     

Thermal stretching in two-phase porous media: Physical basis for Maxwell model

An alternate yet general form of the classical effective thermal conductivity model (Maxwell model) for two-phase porous materials is presented, serving an explicit thermo-physical basis. It is demonstrated that the reduced effective thermal conductivity of the porous media due to non-conducting pore inclusions is caused by the mechanism of thermal stretching, which is a combination of reduced effective heat flow area and elongated heat transfer distance (thermal tortuosity).     

MHD boundary-layer flow of an upper-convected Maxwell fluid in a porous channel

Two-dimensional magnetohydrodynamic (MHD) boundary layer flow of an upper-convected Maxwell fluid is investigated in a channel. The walls of the channel are taken as porous. Using the similarity transformations and boundary layer approximations, the nonlinear partial differential equations are reduced to an ordinary differential equation. The developed nonlinear equation is solved analytically using the homotopy analysis method. An expression for the analytic solution is derived in the form of a series. The convergence of the obtained series is shown. The effects of the Reynolds number Re, Deborah number De and Hartman number M are shown through graphs and discussed for both the suction and injection cases.     

Check on a nonlocal Maxwellian theory of sound propagation in fluid-saturated rigid-framed porous media

A macroscopic nonlocal theory of sound propagation in homogeneous rigid-framed porous media permeated with a viscothermal fluid has been recently proposed in this journal. It accounts for the first time for the full temporal and spatial dispersion effects, independently of the nature of the microgeometry. In this paper this new Maxwellian theory is validated in the case of sound propagation in cylindrical circular tubes, by showing that it matches exactly the long-known direct Kirchhoff–Langevin’s solutions.     

Combined effect of thermal modulation and rotation on the onset of stationary convection in a porous layer

The stability of a fluid-saturated horizontal rotating porous layer subjected to time-periodic temperature modulation is investigated when the condition for the principle of exchange of stabilities is valid. The linear stability analysis is used to study the effect of infinitesimal disturbances. A regular perturbation method based on small amplitude of applied temperature field is used to compute the critical values of Darcy–Rayleigh number and wavenumber. The shift in critical Darcy–Rayleigh number is calculated as a function of frequency of modulation, Taylor number, and Darcy–Prandtl number. It is established that the convection can be advanced by the low frequency in-phase and lower-wall temperature modulation, where as delayed by the out-of-phase modulation. The effect of Taylor number and Darcy–Prandtl number on the stability of the system is also discussed. We found that by proper tuning of modulation frequency, Taylor number, and Darcy–Prandtl number it is possible to advance or delay the onset of convection.     

On the convective instability of a thermal boundary layer

When the surface temperature of a liquid is a harmonic function of time with a frequency, a temperature wave propagates into the liquid. The amplitude of this wave decreases exponentially with distance from the surface. The temperature oscillation is essentially concentrated in a layer of the order of (2/)^1/2, where x is the thermal conductivity of the liquid (thermal boundary layer). Depending on the phase, at certain positions below the surface the temperature gradient is directed downwards and if its magnitude is sufficiently large (the magnitude is a function of the amplitude and frequency of the surface oscillations) the liquid can become unstable with respect to the onset of convection. In that case the convective motion may spread beyond the initial unstable layer. For low frequencies the stability condition can be derived from the usual static Rayleigh criterion, on the basis of the Rayleigh number and the average temperature gradient of the unstable layer. This quasi-static approach, used by Sal'nikov [1], is appropriate to those cases in which the period of the temperature oscillations is much larger than the characteristic time of the perturbations. But when these times are of the same order, the problem must be analyzed in dynamic terms. The stability problem must then be formulated as a problem of parametricresonance excitation of velocity oscillations due to the action of a variable parameter-the temperature gradient.In an earlier work [2] we considered the problem of the stability of a horizontal layer of liquid with a periodically varying temperature gradient. It was assumed that the thickness of the layer was much smaller than the penetration depth of the thermal wave, so that the temperature gradient could be assumed to be independent of position. In the present work we consider the opposite case, in which the liquid layer is assumed to be much larger than the penetration depth, i. e., a thermal boundary layer can be defined. The temperature gradient at equilibrium, which is a parameter in the equations determining the onset of perturbations, is here a periodic function of time and a relatively complicated function of the depth coordinate z. The periodic oscillations are solved by the Fourier method; the equations for the amplitudes are solved by the approximate method of KarmanPohlhausen.The authors are grateful to L. G. Loitsyanskii for helpful criticism.     

Thermal instability in a rotating anisotropic porous layer saturated by a viscoelastic fluid

Linear and non-linear thermal instability in a rotating anisotropic porous medium, saturated with viscoelastic fluid, has been investigated for free-free surfaces. The linear theory is being related to the normal mode method and non-linear analysis is based on minimal representation of the truncated Fourier series analysis containing only two terms. The extended Darcy model, which includes the time derivative and Coriolis terms has been employed in the momentum equation. The criteria for both stationary and oscillatory convection is derived analytically. The rotation inhibits the onset of convection in both stationary and oscillatory modes. A weak non-linear theory based on the truncated representation of Fourier series method is used to find the thermal Nusselt number. The transient behaviour of the Nusselt number is also investigated by solving the finite amplitude equations using a numerical method. The results obtained during the analysis have been presented graphically.     

Linear and nonlinear double diffusive convection in a rotating sparsely packed porous layer using a thermal non-equilibrium model

Double diffusive convection in a fluid-saturated rotating porous layer is studied when the fluid and solid phases are not in local thermal equilibrium, using both linear and nonlinear stability analyses. The Brinkman model that includes the Coriolis term is employed as the momentum equation. A two-field model that represents the fluid and solid phase temperature fields separately is used for the energy equation. The onset criterion for stationary, oscillatory, and finite amplitude convection is derived analytically. It is found that small inter-phase heat transfer coefficient has significant effect on the stability of the system. There is a competition between the processes of thermal diffusion, solute diffusion, and rotation that causes the convection to set in through either oscillatory or finite amplitude mode rather than stationary. The effect of solute Rayleigh number, porosity modified conductivity ratio, Lewis number, diffusivity ratio, Vadasz number, and Taylor number on the stability of the system is investigated. The nonlinear theory based on the truncated representation of Fourier series method predicts the occurrence of subcritical instability in the form of finite amplitude motions. The effect of thermal non-equilibrium on heat and mass transfer is also brought out. Some of the convection systems previously reported in the literature is shown to be special cases of the system presented in this study.     

Weakly nonlinear stability analysis of triple diffusive convection in a Maxwell fluid saturated porous layer

The weakly nonlinear stability of the triple diffusive convection in a Maxwell fluid saturated porous layer is investigated. In some cases, disconnected oscillatory neutral curves are found to exist, indicating that three critical thermal Darcy-Rayleigh numbers are required to specify the linear instability criteria. However, another distinguishing feature predicted from that of Newtonian fluids is the impossibility of quasi-periodic bifurcation from the rest state. Besides, the co-dimensional two bifurcation points are located in the Darcy-Prandtl number and the stress relaxation parameter plane. It is observed that the value of the stress relaxation parameter defining the crossover between stationary and oscillatory bifurcations decreases when the Darcy-Prandtl number increases. A cubic Landau equation is derived based on the weakly nonlinear stability analysis. It is found that the bifurcating oscillatory solution is either supercritical or subcritical, depending on the choice of the physical parameters. Heat and mass transfers are estimated in terms of time and area-averaged Nusselt numbers.     

Coupling model for unsteady MHD flow of generalized Maxwell fluid with radiation thermal transform

This paper introduces a new model for the Fourier law of heat conduction with the time-fractional order to the generalized Maxwell fluid. The flow is influenced by magnetic field, radiation heat, and heat source. A fractional calculus approach is used to establish the constitutive relationship coupling model of a viscoelastic fluid. We use the Laplace transform and solve ordinary differential equations with a matrix form to obtain the velocity and temperature in the Laplace domain. To obtain solutions from the Laplace space back to the original space, the numerical inversion of the Laplace transform is used. According to the results and graphs, a new theory can be constructed. Comparisons of the associated parameters and the corresponding flow and heat transfer characteristics are presented and analyzed in detail.     

Analysis of matching conditions at the boundary surface of a fluid-saturated porous solid and a bulk fluid: the use of Lagrange multipliers

The compatibility conditions matching macroscopic mechanical fields at the contact surface between a fluid-saturated porous solid and an adjacent bulk fluid are considered. The general form of balance equations at that discontinuity surface are analyzed to obtain the compatibility conditions for the tangent and normal components of the velocity and the stress vector fields. Considerations are based on the procedure similar to that used in the phenomenological thermodynamics for derivation of constitutive relations, where the entropy inequality and the concept of Lagrange multipliers are applied. This procedure made possible to derive the compatibility conditions for the viscous fluid flowing tangentially and perpendicularly to the boundary surface of the porous solid and to formulate the generalized form of the so called slip condition for the fluid velocity field, postulated earlier by Beavers and Joseph, J. Fluid. Mech. 30, 197–207 (1967). PACS 47.55.Mh Communicated by Y.D. Shikhmurzaev     

Finite element analysis of convective heat transfer in porous media

The finite element method is used to analyse convective heat transfer in a porous medium. Convection past a vertical surface embedded in the medium and convection in a confined porous medium enclosure are analysed using the above method. The results are compared with those available in the literature and the agreement is found to be good. The method is applicable for two-dimensional analysis in a porous body of any arbitrary shape. The restriction of the boundary layer assumption is relaxed.     

Stability of triple diffusive convection in a viscoelastic fluid-saturated porous layer

The triple diffusive convection in an Oldroyd-B fluid-saturated porous layer is investigated by performing linear and weakly nonlinear stability analyses. The condition for the onset of stationary and oscillatory is derived analytically. Contrary to the observed phenomenon in Newtonian fluids, the presence of viscoelasticity of the fluid is to degenerate the quasiperiodic bifurcation from the steady quiescent state. Under certain conditions, it is found that disconnected closed convex oscillatory neutral curves occur, indicating the requirement of three critical values of the thermal Darcy-Rayleigh number to identify the linear instability criteria instead of the usual single value, which is a novel result enunciated from the present study for an Oldroyd-B fluid saturating a porous medium. The similarities and differences of linear instability characteristics of Oldroyd-B, Maxwell, and Newtonian fluids are also highlighted. The stability of oscillatory finite amplitude convection is discussed by deriving a cubic Landau equation, and the convective heat and mass transfer are analyzed for different values of physical parameters.