Density Maximum Effect on Double-diffusive Natural Convection in a Porous Cavity with Variable Wall Temperature

The effect of density maximum of water on double-diffusive natural convection in a two-dimensioned cavity filled with a water saturated isotropic porous medium is studied numerically. The horizontal walls of the cavity are insulated. The opposing vertical walls are kept at different temperatures θ _h (linearly varies with height) and θ _c (θ _c ≤ θ _h ). The concentration levels at cold wall and hot wall are, respectively, c ₁ and c ₂ with c ₁ > c ₂. Brinkman-Forchheimer extended Darcy model is used to investigate the average heat and mass transfer rates. The non-dimensional equations for momentum, energy, and concentration are solved by finite volume method with power law scheme for convection and diffusion terms. The results are presented in the form of streamlines, isotherms, and isoconcentration lines for various values of Grashof numbers, Schmidt number, porosity, and Darcy numbers. It is observed that the density maximum of water has profound effect on the thermosolutal convection. The effects of different parameters on the velocity, temperature, and species concentrations are also shown graphically.     

Onset of Buoyancy-driven Instability in Porous Media Melted from Below

When a nonhomogeneous solid is melting from below, convection may be induced in a thermally–unstable melt layer. In this study, the onset of buoyancy-driven convection during time-dependent melting is investigated by using similarly transformed disturbance equations. The critical Darcy–Rayleigh numbers based on the melt-layer thickness, Ra _H,c, are found numerically for various conditions. For small superheats, the present predictions show that Ra _H,c is located between 27.1 and 4π ² and it approaches the well-known results of the original Horton–Rogers–Lapwood problem. However, for high superheats, it is dependent on the phase change rate λ and the relation of Ra _H,c λ = 25.89 is shown.     

Stability of an Initially, Stably Stratified Fluid Subjected to a Step Change in Temperature

The onset of convective instability in an initially quiescent, stably stratified fluid layer between two horizontal plates is analyzed with linear theory. The bottom boundary is heated suddenly from below, subjected to a step change in surface temperature. The critical time t _c to mark the onset of Rayleigh-Bénard convection is predicted by propagation theory. This theory uses the length scaled by , where α denotes thermal diffusivity. Under the normal mode analysis the dimensionless disturbance equations are obtained as a function of τ(=αt/d ²) and ζ(=Z/), where d is the fluid layer depth and Z is the vertical distance. The resulting equations are transformed to self-similar ones by using scaling and finally fixing τ as τ_c under the frame of coordinates τ and ζ. For a given γ, Pr and τ_c, the minimum value of Ra is obtained from the marginal stability curve. Here γ denotes the temperature ratio to represent the degree of stabilizing effect, Pr is the Prandtl number and Ra is the Rayleigh number. With γ=0, the minimum Ra value approaches the well-known value of 1708 as τ_c increases. However, it is inversely proportional to τ_c ^3/2 as τ_c decreases. With increasing γ, the system becomes more stable. It is interesting that in the present system, propagation theory produces the stability criteria to bound the available experimental data over the whole domain of time. Received 5 November 2001 and accepted 29 March 2002 Published online: 2 October 2002 RID="*" ID="*" This work has been supported by both SK Chemicals Co. Ltd. and LG Chemical Ltd., Seoul under the Brain Korea 21 Project of the Ministry of Education. Communicated by H.J.S. Fernando     

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     

Nonsimilarity solutions for non-Darcy mixed convection from horizontal surfaces in a porous medium

Non-Darcy mixed convection in a porous medium from horizontal surfaces with variable surface heat flux of the power-law distribution is analyzed. The entire mixed convection regime is divided into two regions. The first region covers the forced convection dominated regime where the dimensionless parameter ζ _f =Ra^* _x /Pe² _x is found to characterize the effect of buoyancy forces on the forced convection with K ^′ U _∞/ν characterizing the effect of inertia resistance. The second region covers the natural convection dominated regime where the dimensionless parameter ζ _n =Pe _x /Ra^*1/2 _x is found to characterize the effect of the forced flow on the natural convection, with (K ^′ U _∞/ν)Ra^*1/2 _x /Pe _x characterizing the effect of inertia resistance. To obtain the solution that covers the entire mixed convection regime the solution of the first regime is carried out for ζ _f =0, the pure forced convection limit, to ζ _f =1 and the solution of the second is carried out for ζ _n =0, the pure natural convection limit, to ζ _n =1. The two solutions meet and match at ζ _f =ζ _n =1, and R ^* _h =G ^* _h . Also a non-Darcy model was used to analyze mixed convection in a porous medium from horizontal surfaces with variable wall temperature of the power-law form. The entire mixed convection regime is divided into two regions. The first region covers the forced convection dominated regime where the dimensionless parameter ξ _f =Ra _x /Pe _x ^3/2 is found to measure the buoyancy effects on mixed convection with Da _x Pe _x /ɛ as the wall effects. The second region covers the natural convection dominated region where ξ _n =Pe _x /Ra _x ^2/3 is found to measure the force effects on mixed convection with Da _x Ra _x ^2/3/ɛ as the wall effects. Numerical results for different inertia, wall, variable surface heat flux and variable wall temperature exponents are presented. Received on 8 July 1996     

A Summary of New Predictive High Frequency Thermo-Vibrational Models in Porous Media

In this paper, we consider the effect of mechanical vibration on the onset of convection in porous media. The porous medium is saturated either by a pure fluid or by a binary mixture. The importance of a transport model on stability diagrams is presented and discussed. The stability threshold for the Darcy–Brinkman case in the Ra _Tc -R and k _c -R diagrams is presented (where Ra _Tc , k _c and R are the critical Rayleigh number, the critical wave number and the vibration parameters, respectively). It is shown that there is a significant deviation from the Darcy model. In the thermo-solutal case with the Soret effect, the influence of vibration on the reduction of multi-cellular convection is emphasized. A new analytical relation for obtaining the threshold of mono-cellular convection is derived. This relation shows how the separation factor Ψ is related to the controlling parameters of the problem, Ψ = f (R, ε*, Le), when the wave number k → 0. The importance of vibrational parameter definition is highlighted and it is shown how, by using a proper definition for vibrational parameter, we may obtain compact relationship. It is also shown how this result may be used to increase component separation.     

Soret effect and slow mass diffusion as a catalyst for overstability in Marangoni-Bénard flows

We study the onset of time dependent Marangoni-Bénard convection in binary mixtures subject to Soret effect by numerical computation of linear instability thresholds in infinite fluid layers and two-dimensional boxes. The calculations are done for positive Marangoni numbers (Ma > 0) and negative Marangoni Soret parameters S _M = –(D _S γ _c )/(Dγ _T ) where D _S and D are the Soret and mass diffusion coefficients, respectively, and γ _T , γ _c are the first derivatives of the surface tension with respect to temperature and concentration. Our purpose is to understand why for particular choices of Prandtl and Schmidt numbers, the increase of the stabilizing solutal contribution leads to a decrease of the critical temperature difference, a phenomenon already reported by Chen & Chen [5] and Skarda et al. [12] For various choices of Prandtl and Schmidt numbers we analyze the evolution of the critical Marangoni number Ma _c , critical wavenumber k _c and angular frequency ω _c with S _M and compute the corresponding eigenvectors. We next propose a physical mechanism which explains how the stabilizing solutal contribution acts as a catalyst for overstability. Finally, we extend our results to two dimensional boxes of small aspect ratio.     

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 a fluid-saturated porous medium contained in a vertical cylinder heated from below by forced convection

A linear stability analysis determining the critical Rayleigh number R _c for onset of convection in a bounded vertical cylinder containing a fluid-saturated porous medium is performed for insulated sidewalls, isothermal top surface, and bottom surface heated by forced convection. This Newtonian heating of the bottom surface involves a Biot number Bi that allows consideration of the continuum of boundary conditions ranging from constant heat flux, with global minimum R _min=27.096 found as Bi→0, to isothermal, with global minimum R _min=4π² found as Bi→ ∞. In both cases and for most cylinder aspect ratios, incipient convection sets in as an asymmetric mode, though islands of aspect ratio exist where the onset mode is symmetric. Sample three-dimensional renderings of disturbance temperature distributions showing preferred modes at onset of convection for fixed Bi are provided and an analytical fit to R _min as a function of Bi is given.     

Streamwise evolution of a high-Schmidt-number passive scalar in a turbulent plane wake

The concentration fluctuation c of diluted fluorescein dye, a high-Schmidt-number passive scalar (Sc=ν/D ≈ 2000, ν and D are the fluid momentum and dye diffusivities, respectively), is measured in the wake of a circular cylinder using a single-point laser-induced fluorescence (SPLIF) technique. The streamwise decay rate of the mean and rms values of c is slow in comparison to that of θ, the temperature fluctuation for which the molecular Prandtl number Pr=ν/κ is about 0.7 (κ is the thermal diffusivity). The comparison between mean and rms distributions of c and θ highlights the combined role the Reynolds and Schmidt numbers play in terms of dispersing the scalar. The streamwise evolution of the probability density functions (pdfs) of c and θ suggest that while p(θ) is approximately Gaussian in the intermediate wake (x/d ≈ 80), p(c) is strongly non-Gaussian, and depends on both x/d and Re. The skewness of c is larger than that of θ along the wake centreline. Arguably, the asymmetry of p(c) reflects the relatively strong organisation of the large-scale motion in the far-wake. Received: 27 July 2000/Accepted: 22 December 2000     

Influence of inclination angle on buoyancy-driven convection in triangular enclosure filled with a fluid-saturated porous medium

The present investigation deals with the numerical analysis of steady-state laminar buoyancy-driven convection in an inclined triangular enclosure filled with fluid saturated porous media using the Darcy law equation. One wall of the enclosure is isothermally heated and the other is cooled, while the remaining wall is adiabatic. The effect of inclination angle on natural convection is investigated by varying the angle of inclination (φ) between 0° and 360°. The governing transformed equations are solved numerically using a finite-difference method. Obtained results are shown in the form of streamlines, isotherms, mean Nusselt numbers and dimensionless stream function for different values of the Rayleigh number Ra in the range 100 ≤ Ra ≤ 1,000. It is found that the values of the maximum and minimum mean Nusselt number are reached for φ = 330° and φ = 210° , respectively. However, the lowest flow strength is formed at φ = 240° for all values of Ra.     

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

The effect of vertical heterogeneity of permeability, on the onset of convection in a horizontal layer of a saturated porous medium, uniformly heated from below but with a non-uniform basic temperature gradient resulting from vertical throughflow, is studied analytically using linear stability theory. It is found that, to first order, a linear variation of the reciprocal of permeability with depth has no effect on the critical value of the Rayleigh number Ra _c based on the harmonic mean of the permeability, but a quadratic variation increasing in the upwards direction leads to a reduction in Ra _c.     

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.     

Prediction of natural convection flow using network model and numerical simulations inside enclosure with distributed solid blocks

Steady state natural convection of a fluid with Pr ≈ 1 within a square enclosure containing uniformly distributed, conducting square solid blocks is investigated. The side walls are subjected to differential heating, while the top and bottom ones are kept adiabatic. The natural convection flow is predicted employing the nondimensional volumetric flow rate (Q_max^* Q_{\max }^{*} ) by using a network model and also using numerical simulations. For identical solid and fluid thermal conductivities (i.e. k _s  = k _f ), a parametric study of the effect of number of blocks (N ²), gap size (δ) and enclosure Rayleigh number (Ra) on Q_max^* Q_{\max }^{*} is performed using the two approaches. Network model predictions are observed to agree well with that from the simulations until Raδ³ ~ 12. Considering the enclosure with blocks as a porous medium, for a fixed enclosure Ra number, increasing the number of blocks for a fixed volumetric porosity leads to a decrease in enclosure permeability, which in turn reduces the flow rate. When the number of blocks is fixed, and for a given Ra number, the flow rate increases as the porosity increases by widening the gap between the blocks.     

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 γ.     

Nusselt-Rayleigh correlations for free convection in 2D air-filled parallelogrammic enclosures with isothermal active walls

In the present study Nu-Ra-α correlations are proposed to calculate the steady-state natural convection heat transfer taking place in 2D air-filled cavities of parallelogrammic section. The thermal conditions and the dimensions of the enclosures permit to cover a large range of Rayleigh numbers, 1.7 × 10³  ≤ Ra ≤ 3.0 × 10⁹, suitable for diverse engineering applications. The two active walls of the cavities are kept vertical and isothermal at hot and cold temperatures T _h and T _c respectively. Separated by a horizontal distance H, they have the same height H and are connected by a closed adiabatic channel whose upper and lower walls can be inclined at an angle α with respect to the horizontal, varying between −60° to +60°. That gives rise to a conducting or insulating cavity, in the convective sense of the term (diode cavity). A computational model based on the finite volume method is used to solve the governing equations. The large number of treated configurations led to propose Nu-Ra-α correlations for large ranges of Ra and α which can be applied to many engineering areas. The results of this numerical study have been successfully compared with calculated and measured available data.     

Laminar thermally developing flow inside right-angularly triangular ducts

An analytical approach based on the generalized integral transform technique is presented, for the solution of laminar forced convection within the thermal entry region of ducts with arbitrarily shaped cross-sections. The analysis is illustrated through consideration of a right triangular duct subjected to constant wall temperature boundary condition. Critical comparisons are made with results available in the literature, from direct numerical approaches. Numerical results for dimensionless average temperature and Nusselt numbers are presented for different apex angles.Nomenclature a,b sides of right triangular duct - A _c cross-sectional area of duct - c _p specific heat of fluid - D _h =4A _c /p hydraulic diameter, with P the wet perimeter - h(z) heat transfer coefficient at duct wall - k thermal conductivity - Pe=c _p D _h /k Peclet number - T(x, y, z) temperature distribution - T _o inlet temperature - T _w prescribed wall temperature - u(x, y); U(X, Y) dimensional and dimensionless velocity profile - average flow velocity - x; X dimensional and dimensionless normal coordinate (Fig. 1) - x ₁(y); X ₁(Y) dimensional and dimensionless position at irregular boundary (Fig. 1) - y; Y dimensional and dimensionless normal coordinate (Fig. 1) - z; Z dimensional and dimensionless axial coordinate Greek letters side of right triangular duct in X direction (dimensionless) - side of right triangular duct in Y direction (dimensionless) - density of fluid - (X, Y, Z) dimensionless temperature distribution - * apex angle of triangular duct (Fig. 1) - ** apex angle of triangular duct (Fig. 1)     

Thermal stability of composite superconducting tape under the effect of a two-dimensional dual-phase-lag heat conduction model

Thermal stability of composite superconducting tape subjected to a thermal disturbance is numerically investigated under the effect of a two-dimensional dual-phase-lag heat conduction model. It is found that the dual-phase-lag model predicts a wider stable region as compared to the predictions of the parabolic and the hyperbolic heat conduction models. The effects of different design, geometrical and operating conditions on superconducting tape thermal stability were also studied.a conductor width, (m) - A conductor cross sectional area of, (m²) - As conductor aspect ratio, (a/b) - b conductor thickness, (m) - Bi Biot number - B dimensionless disturbance Intensity - C heat capacity, (J m^–3 K^–1) - D disturbance energy density, (W m^–3) - f volume fraction of the stabilizer in the conductor - g(T) steady capacity of the Ohmic heat source, (W m^–3) - g_max Ohmic heat generation with the whole current in the stabilizer, (W m^–3) - G_max dimensionless maximum Joule heating - h convective heat transfer coefficient, (W m^–2 K^–1) - J current density, (A m^–2) - k thermal conductivity of conductor, (W m^–1 K^–1) - q conduction heat flux vector, (W m^–2) - Q dimensionless Joule heating - R relaxation times ratio (_T/2_q) - t rime, (s) - T temperature, (K) - T_c critical temperature, (K) - T_c1 current sharing temperature, (K) - T_i initial temperature, (K) - T_o ambient temperature, (K) - x, y co-ordinate defined in Fig. 1, (m) - thermal diffusivity (m² s^–1) - dimensionless time - _i dimensionless duration time - dimensionless y-variable - _o superconductor dimensionless thickness - dimensionless temperature - _c1 dimensionless current sharing temperature - 1 dimensionless maximum temperature - dimensionless disturbance energy - numerical tolerance - x width of conductor subjected to heat disturbances, (m) - y thickness of conductor subjected to heat disturbances, (m) - dimensionless x-variable - _o superconductor dimensionless width - stabilizer electrical resistivity, () - _q relaxation time of heat flux, (s) - _T relaxation time of temperature gradient, (s) - i initial - sc current sharing - max maximum - o ambient     

Average and extremal properties of heat transfer and shear stress on a wall surface in Rayleigh–Bénard convection

In this study surface-averaged and extremal properties of heat transfer and shear stress on the upper wall surface of Rayleigh–Bénard convection are numerically examined. The Prandtl number was raised up to 10³, and the Rayleigh number was changed between 10⁴ and 10⁷. As a result, average Nusselt number Nu and shear rate τ/Pr depends on Pr, Ra, and the entire numerical results are distributed between two correlation equations corresponding to small and large Pr. The small and large Pr equations are closely related to steady and unsteady flow regimes, respectively. Nevertheless, a single relation τ/Pr ~ Nu ^3.0 exists to explain the entire results. Similarly the change of local maximal properties Nu _max and τ _max/Pr depends on Pr, Ra, and these values are also distributed between two correlation equations corresponding to small and large Pr cases. Despite such complicated dependence we can obtain a correlation equation as a form of τ _max/Pr ~ Nu _max^2.6, which has not been obtained theoretically.