Steady thermocapillary-buoyant convection in a shallow annular pool. Part 1: Single layer fluid

This paper examines the steady thermocapillarybuoyant convection in a shallow annular pool subjected to a radial temperature gradient. A matched asymptotic theory is used to obtain the asymptotic solutions of the flow and thermal fields in the case of small aspect ratios,which is defined as the ratio of the layer thickness to the gap width. The flow domain is divided into the core region away from the cylinder walls and two end regions near each cylinder wall. Asymptotic solutions are obtained in the core region by solving the core and end flows separately and then joining them through matched asymptotic expansions. For the system of silicon melt,the asymptotic solutions are compared with the results of numerical simulations. It is found that the two kinds of solutions have a good agreement in the core region for a small aspect ratio. With the increase of aspect ratio,the applicability of the present asymptotic solutions decreases gradually.     

Stability of Gravity Driven Convection in a Cylindrical Porous Layer Subjected to Vibration

The linear stability theory is used to investigate analytically the effects of gravity modulation on convection in a homogeneous cylindrical porous layer heated from below. The linear stability results show that increasing the frequency of vibration stabilizes the convection. In addition the aspect ratio of the porous cylinder is shown to influence the stability of convection for all frequencies analysed. It was also observed that only synchronous solutions are possible in cylindrical porous layers, with no transition to subharmonic solutions as was the case in Govender (2005a) [Transport Porous Media 59(2), 227–238] for rectangular layers or cavities.     

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.     

Convection regime flow in a vertical slot: continuum of solutions from capped to open ends

In a recent publication Bühler (Heat Mass Transfer 39:631–638, 2003) reported new results for conduction regime flow between vertical differentially-heated walls that provide a continuum of solutions between capped and open ends. In this paper we extend Bühler’s work to realize a continuum of solutions of convection regime flow using empirical results for the vertical temperature gradient that develops in tall aspect ratio geometries. The mass flux is determined analytically for this three-parameter family of solutions. Identical viscous and thermal boundary layers exist at the opposing walls when the cavity is capped. However, as the flow evolves to one with open ends, there is an intensification (attenuation) of the boundary layers near the hot (cold) walls. In the limit corresponding to an open-ended cavity, the boundary layer at the cold wall vanishes altogether.     

Entropy generation in parallel plate microchannels

In this study the entropy generation in microchannels in microdevices induced by the transient laminar forced convection in the combined entrance region between two parallel plates has been investigated numerically. The study considers the microscales in the region of Kn < 0.001. The effects of aspect ratio, Reynolds number, Prandtl number, Brinkman number, and the motion of the lower plate on the entropy generation during the simultaneously developing flow in a parallel-plates channel are investigated. The obtained results addressing all cases are thoroughly in good agreement with the expectations that the entropy generation has its highest value at channel with the smallest aspect ratio at counter motion of the lower plate with the highest Re, Pr and Br/Ω values considered in the problem. An erratum to this article can be found at     

Variable permeability effect on convection in binary mixtures saturating a porous layer

The Darcy Model with the Boussinesq approximation is used to study natural convection in a shallow porous layer, with variable permeability, filled with a binary fluid. The permeability of the medium is assumed to vary exponentially with the depth of the layer. The two horizontal walls of the cavity are subject to constant fluxes of heat and solute while the two vertical ones are impermeable and adiabatic. The governing parameters for the problem are the thermal Rayleigh number, R _T, the Lewis number, Le, the buoyancy ratio, φ, the aspect ratio of the cavity, A, the normalized porosity, ε, the variable permeability constant, c, and parameter a defining double-diffusive convection (a = 0) or Soret induced convection (a = 1). For convection in an infinite layer, an analytical solution of the steady form of the governing equations is obtained on the basis of the parallel flow approximation. The onset of supercritical convection, or subcritical, convection are predicted by the present theory. A linear stability analysis of the parallel flow model is conducted and the critical Rayleigh number for the onset of Hopf’s bifurcation is predicted numerically. Numerical solutions of the full governing equations are found to be in excellent agreement with the analytical predictions.     

Free convection from a vertical permeable circular cone with non-uniform surface heat flux

 In this paper, the problem of laminar free convection from a vertical permeable circular cone maintained with non-uniform surface heat flux is considered. The governing boundary layer equations are reduced non-similar boundary layer equations with surface heat flux proportional to x ⁿ (where x is the distance measured from the leading edge). The solutions of the reduced equations are obtained by using three distinct solution methodologies; namely, (i) perturbation solution for small transpiration parameter, ξ, (ii) asymptotic solution for large ξ, and (iii) the finite difference solutions for all ξ. The solutions are presented in terms of local skin-friction and local Nusselt number for smaller values of Prandtl number and heat flux gradient and are displayed in tabular form as well as graphically. Effects of pertinent parameters on velocity and temperature profiles are also shown graphically. Solutions obtained by finite difference method are also compared with the perturbation solutions for small and large ξ and found to be in excellent agreement. Received on 1 October 1999     

Combined Laminar Mixed Convection and Surface Radiation using Asymptotic Computational Fluid Dynamics (ACFD)

This paper reports the use of the technique of combining asymptotics with computational fluid dynamics (CFD), known as asymptotic computational fluid dynamics (ACFD), to handle the problem of combined laminar mixed convection and surface radiation from a two dimensional, differentially heated lid driven cavity. The fluid under consideration is air, which is radiatively transparent, and all the walls are assumed to be gray and diffuse and having the same hemispherical, total emissivity (ɛ). The computations have been performed on FLUENT 6.2. The full radiation problem (i.e. all the walls are radiatively black corresponding to ɛ = 1) is first taken up and the method of “perturbing and blending” is used wherein, first, limiting solutions of natural and forced convection are perturbed, to obtain correlations for the weighted average convective Nusselt numbers for the full radiation case. These correlations are then blended suitably in order to obtain a composite correlation for the weighted average convective Nusselt number that is valid for the entire mixed convection range, i.e., 0 ≤ Ri ≤ ∞. This correlation is then expanded in terms of ɛ to obtain an expression for the average convective Nusselt number that is valid for any ɛ in the range 0 ≤ ɛ ≤ 1. In so far as radiation heat transfer is concerned, using asymptotic arguments, a new weighted average radiation Nusselt number is defined such that this quantity can be expanded just in terms of ɛ. Hence, by the use of ACFD, the number of solutions required to obtain reasonably accurate correlations for both the convective and radiative heat transfer rates and hence the total heat transfer rate (Nu _total = Nu _C + Nu _R), is substantially reduced. More importantly, the correlations for convection and radiation are asymptotically correct at their ends. The effect of secondary variables like aspect ratio and the case of unequal wall emissivities can also be included without significant additional effort.     

Double-diffusive convection for a heated cylinder submerged in a salt-stratified fluid layer

Double-diffusive convection due to a cylindrical source submerged in a salt-stratified solution is numerically investigated in this study. For proper simulation of the vortex generated around the cylinder, a computational domain with irregular shape is employed. Flow conditions depend strongly on the thermal Rayleigh number, Ra _T , and the buoyancy ratio, R _ρ. There are two types of onset of instability existing in the flow field. Both types are due to either the interaction of the upward temperature gradient and downward salinity gradient or the interaction of the lateral temperature gradient and downward salinity gradient. The onset of layer instability due to plume convection is due to the former, whereas, the onset of layer instability of layers around the cylinder is due to the latter. Both types can be found in the flow field. The transport mechanism of layers at the top of the basic plume belongs to former while that due to basic plume and layer around the cylinder are the latter. The increase in Ra _T reinforces the plume convection and reduces the layer numbers generated around the cylinder for the same buoyancy ratio. For the same Ra _T , the increase of R _ρ suppresses the plume convection but reinforces the layers generated around the cylinder. The profiles of local Nusselt number reflects the heat transfer characteristics of plume convection and layered structure. The profiles of averaged Nusselt number are between the pure conduction and natural convection modes and the variation is due to the evolution of layers. Received on 13 September 1996     

Stability of double-diffusive convection induced by selective absorption of radiation in a fluid layer

Linear and nonlinear stability analyses were performed on a fluid layer with a concentration-based internal heat source. Clear bimodal behaviour in the neutral curve (with stationary and oscillatory modes) is observed in the region of the onset of oscillatory convection, which is a previously unobserved phenomenon in radiation-induced convection. The numerical results for the linear instability analysis suggest a critical value γ _c of γ, a measure for the strength of the internal heat source, for which oscillatory convection is inhibited when γ > γ _c . Linear instability analyses on the effect of varying the ratio of the salt concentrations at the upper and lower boundaries conclude that the ratio has a significant effect on the stability boundary. A nonlinear analysis using an energy approach confirms that the linear theory describes the stability boundary most accurately when γ is such that the linear theory predicts the onset of mostly stationary convection. Nevertheless, the agreement between the linear and nonlinear stability thresholds deteriorates for larger values of the solute Rayleigh number for any value of γ.     

Exact solutions for the flow of second grade fluid in annulus between torsionally oscillating cylinders

The velocity field and the associated shear stress corresponding to the torsional oscillatory flow of a second grade fluid, between two infinite coaxial circular cylinders, are determined by means of the Laplace and Hankel transforms. At time t = 0, the fluid and both the cylinders are at rest and at t = 0 + , cylinders suddenly begin to oscillate around their common axis in a simple harmonic way having angular frequencies ω 1 and ω 2 . The obtained solutions satisfy the governing differential equation and all imposed initial and boundary conditions. The solutions for the motion between the cylinders, when one of them is at rest, can be obtained from our general solutions. Furthermore, the corresponding solutions for Newtonian fluid are also obtained as limiting cases of our general solutions.     

Steady mixed convection boundary-layer flow over a vertical flat surface in a porous medium filled with water at 4°C: variable surface heat flux

The steady mixed convection boundary-layer flow over a vertical impermeable surface in a porous medium saturated with water at 4°C (maximum density) when the surface heat flux varies as x ^m and the velocity outside the boundary layer varies as x ^(1+2m)/2, where x measures the distance from the leading edge, is discussed. Assisting and opposing flows are considered with numerical solutions of the governing equations being obtained for general values of the flow parameters. For opposing flows, there are dual solutions when the mixed convection parameter λ is greater than some critical value λ _c (dependent on the power-law index m). For assisting flows, solutions are possible for all values of λ. A lower bound on m is found, m > −1 being required for solutions. The nature of the critical point λ _c is considered as well as various limiting forms; the forced convection limit (λ = 0), the free convection limit (λ → ∞) and the limits as m → ∞ and as m → −1.     

Natural Convection Induced by a Centrifugal Force Field in a Horizontal Annular Porous Layer Saturated with a Binary Fluid

The Darcy Model with the Boussinesq approximation is used to study natural convection in a horizontal annular porous layer filled with a binary fluid, under the influence of a centrifugal force field. Neumann boundary conditions for temperature and concentration are applied on the inner and outer boundary of the enclosure. The governing parameters for the problem are the Rayleigh number, Ra, the Lewis number, Le, the buoyancy ratio, j{\varphi } , the radius ratio of the cavity, R, the normalized porosity, e{\varepsilon } , and parameter a defining double-diffusive convection (a = 0) or Soret induced convection (a = 1). For convection in a thin annular layer (R → 1), analytical solutions for the stream function, temperature and concentration fields are obtained using a concentric flow approximation and an integral form of the energy equation. The critical Rayleigh number for the onset of supercritical convection is predicted explicitly by the present model. Also, results are obtained from the analytical model for finite amplitude convection for which the flow and heat and mass transfer are presented in terms of the governing parameters of the problem. Numerical solutions of the full governing equations are obtained for a wide range of the governing parameters. A good agreement is observed between the analytical model and the numerical simulations.     

Penetrative convection in sublayer of water including density inversion

Numerical solutions of stability and convective flow in an infinite horizontal water layer, including density inversion, have been obtained using a finite element code. The evolution of the temperature field and flow pattern near the onset of convection are studied in detail. It is known that natural convection develops primarily in the lower unstably stratified layer. Of interest is the penetration of the convection rolls into the upper stably stratified layer and concurrent liquid entrainment as a function of the increasing Rayleigh number at different aspect ratios. Individual convection rolls may grow and expand before splitting up into two roll cells. It is shown that changing the aspect ratio influences critical Rayleigh number, flow symmetry, flow pattern, and transitions between flow patterns. Numerical results on heating from above or from below, agree well with available results in the literature. A correlation to predict critical Rayleigh numbers is given for the case of heating from above.     

Non-Darcy mixed convection from a vertical plate in saturated porous media-variable surface heat flux

Nonsimilarity solutions for non-Darcy mixed convection from a vertical impermeable surface embedded in a saturated porous medium are presented for variable surface heat flux (VHF) of the power-law form. The entire mixed convection region is divided into two regimes. One region covers the forced convection dominated regime and the other one covers the natural convection dominated regime. The governing equations are first transformed into a dimensionless form by the nonsimilar transformation and then solved by a finite-difference scheme. Computations are based on Keller Box method and a tolerance of iteration of 10⁻⁵ as a criterion for convergence. Three physical aspects are introduced. One measures the strength of mixed convection where the dimensionless parameter Ra^* _x /Pe^3/2 _x characterizes the effect of buoyancy forces on the forced convection; while the parameter Pe _x /Ra^*2/3 _x characterizes the effect of forced flow on the natural convection. The second aspect represents the effect of the inertial resistance where the parameter K′U _∞/ν is found to characterize the effect of inertial force in the forced convection dominated regime, while the parameter (K′U _∞/ν)(Ra^*2/3 _x /Pe _x ) characterizes the effect of inertial force in the natural convection dominated regime. The third aspect is the effect of the heating condition at the wall on the mixed convection, which is presented by m, the power index of the power-law form heating condition. Numerical results for both heating conditions are carried out. Distributions of dimensionless temperature and velocity profiles for both Darcy and non-Darcy models are presented. Received on 26 May 1997     

Effect of the Source Term on Steady Free Convection Boundary Layer Flows over an Vertical Plate in a Porous Medium. Part I

The Darcy free convection boundary layer flow over a vertical flat plate is considered in the presence of volumetric heat generation/absorption. In the present first part of the paper it is assumed that the heat generation/absorption takes place in a self-consistent way, the source term q ^′′′≡ S of the energy equation being an analytical function of the local temperature difference T − T _∞. In a forthcoming second part, the case of the externally controlled source terms S = S(x,y ) will be considered. It is shown that due to the presence of S, the physical equivalence of the up- and downflows gets in general broken, in the sense that the free convection flow over the upward projecting hot plate (“upflow”) and over its downward projecting cold counterpart (“downflow”) in general become physically distinct. The consequences of this circumstance are examined for different forms of S. Several analytical solutions are given. Some of them describe algebraically decaying boundary layers which can also be recovered as limiting cases of exponentially decayingones. This asymptotic phenomenon is discussed in some detail.     

The onset of buoyancy-driven convection in a ferromagnetic fluid saturated porous medium

The onset of buoyancy-driven convection in an initially quiescent ferrofluid saturated horizontal porous layer in the presence of a uniform vertical magnetic field is investigated. The Brinkman-Lapwood extended Darcy equation with fluid viscosity different from effective viscosity is used to describe the flow in the porous medium. The lower boundary of the porous layer is assumed to be rigid-paramagnetic, while the upper paramagnetic boundary is considered to be either rigid or stress-free. The thermal conditions include fixed heat flux at the lower boundary, and a general convective–radiative exchange at the upper boundary, which encompasses fixed temperature and fixed heat flux as particular cases. The resulting eigenvalue problem is solved numerically using the Galerkin technique. It is found that increase in the Biot number Bi, porous parameter σ, viscosity ratio Λ, magnetic susceptibility χ, and decrease in the magnetic number M ₁ and non-linearity of magnetization M ₃ is to delay the onset of ferroconvection in a porous medium. Further, increase in M ₁, M ₃, and decrease in χ, Λ, σ and Bi is to decrease the size of convection cells.     

A class of higher order compact schemes for the unsteady two‐dimensional convection–diffusion equation with variable convection coefficients

A class of higher order compact (HOC) schemes has been developed with weighted time discretization for the two‐dimensional unsteady convection–diffusion equation with variable convection coefficients. The schemes are second or lower order accurate in time depending on the choice of the weighted average parameter μ and fourth order accurate in space. For 0.5?μ?1, the schemes are unconditionally stable. Unlike usual HOC schemes, these schemes are capable of using a grid aspect ratio other than unity. They efficiently capture both transient and steady solutions of linear and nonlinear convection–diffusion equations with Dirichlet as well as Neumann boundary condition. They are applied to one linear convection–diffusion problem and three flows of varying complexities governed by the two‐dimensional incompressible Navier–Stokes equations. Results obtained are in excellent agreement with analytical and established numerical results. Overall the schemes are found to be robust, efficient and accurate. Copyright © 2002 John Wiley & Sons, Ltd.     

Steady mixed convection boundary layer flow over a vertical flat plate in a porous medium filled with water at 4°C: case of variable wall temperature

The problem of steady mixed convection boundary layer flow over a vertical impermeable flat plate in a porous medium saturated with water at 4°C (maximum density) when the temperature of the plate varies as x ^m and the velocity outside boundary layer varies as x ^{2 m} , where x measures the distance from the leading edge of the plate and m is a constant is studied. Both cases of the assisting and the opposing flows are considered. The plate is aligned parallel to a free stream velocity U _∞ oriented in the upward or downward direction, while the ambient temperature is T _∞ = T _m (temperature at maximum density). The mathematical models for this problem are formulated, analyzed and simplified, and further transformed into non-dimensional form using non-dimensional variables. Next, the system of governing partial differential equations is transformed into a system of ordinary differential equations using the similarity variables. The resulting system of ordinary differential equations is solved numerically using a finite-difference method known as the Keller-box scheme. Numerical results for the non-dimensional skin friction or shear stress, wall heat transfer, as well as the temperature profiles are obtained and discussed for different values of the mixed convection parameter λ and the power index m. All the numerical solutions are presented in the form of tables and figures. The results show that solutions are possible for large values of λ and m for the case of assisting flow. Dual solutions occurred for the case of opposing flow with limited admissible values of λ and m. In addition, separation of boundary layers occurred for opposing flow, and separation is delayed for the case of water at 4°C (maximum density) compared to water at normal temperature.     

The criterion of the existence or inexistence of transverse shock wave at wedge supported oblique detonation wave

A simplified theoretic method and numerical simulations were carried out to investigate the characterization of propagation of transverse shock wave at wedge supported oblique detonation wave.After solution validation,a criterion which is associated with the ratio Φ (u 2 /u CJ) of existence or inexistence of the transverse shock wave at the region of the primary triple was deduced systematically by 38 cases.It is observed that for abrupt oblique shock wave (OSW)/oblique detonation wave (ODW) transition,a transverse shock wave is generated at the region of the primary triple when Φ < 1,however,such a transverse shock wave does not take place for the smooth OSW/ODW transition when Φ > 1.The parameter Φ can be expressed as the Mach number behind the ODW front for stable CJ detonation.When 0.9 < Φ < 1.0,the reflected shock wave can pass across the contact discontinuity and interact with transverse waves which are originating from the ODW front.When 0.8 < Φ < 0.9,the reflected shock wave can not pass across the contact discontinuity and only reflects at the contact discontinuity.The condition (0.8 < Φ < 0.9) agrees well with the ratio (D ave /D CJ) in the critical detonation.