Time evolution of double-diffusive convection in a vertical cylinder with radial temperature and axial solutal gradients

The characteristics of transient double-diffusive convection in a vertical cylinder are numerically simulated using a finite element method. Initially the fluid in the cavity is at uniform temperature and solute concentration, then constant temperature and solute concentration, which are lower than their initial values, are imposed along the sidewall and bottom wall, respectively. The time evolution of the double-diffusive convection is investigated for specific parameters, which are the Prandtl number, Pr = 7, the Lewis number, Le = 5, the thermal Grashof number, Gr_T = 10⁷, and the aspect ratio, A = 2, of the enclosure. The objective of the work is to identify the effect of the buoyancy ratio (the ratio of solutal Grashof to thermal Grashof numbers: N = Gr_S/Gr_T) on the evolution of the flow field, temperature and solute field in the cavity. It is found that initially the fluid near the bottom wall is squeezed by the cold flow from the sidewall, a crest of the solute field forms and then pushed to the symmetry line. In the case of N > 0, a domain with higher temperature and weak flow (dead region) forms on the bottom wall near the symmetry line, and the area of dead region increases when N varies from 0.5 to 1.5. More crests of the solute field are formed and the flow near the bottom wall fluctuates continuously for N < 0. The frequency of the fluctuation increases when N varies from −0.5 to −1.5. Corresponding to the variety of the thermal and solutal boundary layers, the average rates of heat transfer (Nu) at the sidewall remain almost unchanged while the average rates of mass transfer (Sh) at the bottom wall change much in the cases of N = 1, 0, −1.     

Oscillatory double-diffusive mixed convection in a two-dimensional ventilated enclosure

By starting from a steady flow configuration based on the work of Deng et al. [Qi-Hong Deng, Jiemin Zhou, Chi Mei, Yong-Ming Shen, Fluid, heat and contaminant transport structures of laminar double-diffusive mixed convection in a two-dimensional ventilated enclosure, Int. J. Heat Mass Transfer 47 (2004) 5257–5269], a numerical investigation was conducted to analyse the unsteady double-diffusive mixed convection in two-dimensional ventilated room due to heat and contaminant sources. Owing to the large number of parameters, the results are reported only for a constant buoyancy ratio N equal to 1. The flow is found to be oscillatory for a fixed Reynolds number (700 ⩽ Re ⩽ 1000) when the Grashof number is varied in a wide range (10³ ⩽ Gr ⩽ 10⁶). Results of the simulations show that the onset of the oscillatory indoor airflow occurs for couples (Re, Gr) values that can be correlated as Re = aGr^b.     

Soret effect on the boundary layer flow regime in a vertical porous enclosure subject to horizontal heat and mass fluxes

The combined thermo- and double-diffusive convection in a vertical tall porous cavity subject to horizontal heat and mass fluxes was investigated analytically and numerically using the Darcy model with the Boussinesq approximation. The investigation focused on the effect of Soret diffusion on the boundary layer flow regime. The governing parameters were the thermal Rayleigh number, R_T, the Lewis number, Le, the buoyancy ratio, N, the Soret parameter, M, which characterized the Soret effect, and the aspect ratio of the enclosure, A_r. The results demonstrated the existence of a boundary layer flow solution for which the Soret parameter had a strong effect on the heat and mass transfer characteristics. For M ≠ 1 and M ≠ −1/Le, the profiles of the vertical velocity component, v, temperature, T, and solute concentration, S, exhibited boundary layer behaviors at high Rayleigh numbers. Furthermore, as R_T increased, the temperature and solute concentration became vertically and linearly stratified in the core region of the enclosure. The thermo-diffusion effect on the boundary layer thickness, δ, was discussed for a wide range of the governing parameters. It was demonstrated analytically that the thickness of the boundary layer could either increase or decrease when the Soret parameter was varied, depending on the sign of the buoyancy ratio. The effect of R_T on the fluid flow properties and heat and mass transfer characteristics was also investigated.     

Unsteady mixed convection flow from a moving vertical plate in a parallel free stream: Influence of heat generation or absorption

An unsteady mixed convection flow over a moving vertical plate in a parallel free stream is considered to investigate the combined effects of buoyancy force and thermal diffusion in presence of heat generation or absorption. The unsteadiness is introduced by the time dependent free stream velocity as well as by the moving plate velocity. The governing boundary layer equations are transformed into a non-dimensional form by a special group of non-similar transformations. The resulting system of coupled non-linear partial differential equations is solved by an implicit finite difference scheme in combination with the quasi-linearization technique. Computations are performed and numerical results are displayed graphically to illustrate the influence of the buoyancy (mixed convection) parameter, Prandtl number, the ratio of free stream velocity to the composite reference velocity and heat generation or absorption parameter on the velocity and temperature profiles. The numerical results for the local skin-friction coefficient and local Nusselt number are also presented. Present results are compared with previously published work and are found to be in excellent agreement. It is found that in presence of buoyancy force (λ > 0), the velocity profile exhibits velocity overshoot 80% more for lower Prandtl number (Pr = 0.7) as compared to the magnitude of the velocity overshoot for higher Prandtl number (Pr = 7.0).     

Numerical study of turbulent double-diffusive natural convection in a square cavity by LES-based lattice Boltzmann model

Turbulent double-diffusive natural convection in a square cavity represents numerous important problems in practice as well as in fundamental. However up to date the study on it is quite sparse and most previous studies just focus on laminar regime. To the best knowledge of the present authors, only several k–? models were developed to investigate turbulent double-diffusive convection and there is no attempt to use Large Eddy Simulation (LES). In order to deepen our knowledge on turbulent double-diffusive convection in a square cavity, we propose a novel LES-based lattice Boltzmann (LB) model to simulate such turbulent convectional flow. Previous LES-based LB models can be recovered from the present model. We find that the symmetry of the fluid circulation becomes broken since the Rayleigh number Ra = 10⁸, although the asymmetry is more clear when Ra ? 10¹⁰. More important, in the present study we find the power-law relationship among the Nusselt (Nu), the ratio of buoyancy forces (N) and the Rayleigh number (Ra) still exists in turbulent regime. The formula among them can be concluded as Nu = a × (Ra × ∣1 ? N∣)^b + c. The values of parameters a, b and c are given in this work.     

Influence of periodicity of sinusoidal bottom boundary condition on natural convection in porous enclosure

In this paper, natural convective flow within a rectangular enclosure has been investigated numerically. All the walls of the enclosure are adiabatic except the bottom wall, which is partially heated and cooled by sinusoidal temperature profile. Both situations: medium is hydro-dynamically isotropic and anisotropic are considered. The governing equations are written under assumption of Brinkman-extended non-Darcy model, including material derivative, and then solved by numerically using spectral element method (SEM). Main emphasize is given on effect of periodicity parameter (N) on local heat transfer rate (Nu_x) as well as flow mechanism in the enclosure. The result shows that, as the periodicity is decreased on increasing N, the absolute value of Nu_x at the bottom left corner point increases in both situations. For odd values of N, the local heat transfer profile is symmetric about the line X = 0.5. The entire flow is governed by two type convective cells: (i) rotating clockwise (ii) rotating anticlockwise. Furthermore, increasing of N increases the multiple cellular structures in the form of N + 1. The relative impact of media permeability on Nu_x at (0, 0) shows that enhancement of media permeability along Y-direction is much more effective than the same along X-direction. In contrast to isotropic porous media, the effect of the permeability and permeability orientation angle on the flow is complex and often non intuitive. In particular the present analysis shows that, different periodicity of temperature boundary condition has significant effect on the flow mechanism and consequently on the heat transfer rate for both situations.     

Jet impingement interaction with cross flow in horizontal porous layer under thermal non-equilibrium conditions

In the present article the jet impingement cooling of heated portion of a horizontal surface immersed in a thermally non-equilibrium porous layer is considered for investigation numerically with the presence of a cross flow. The mathematical model is derived for steady, two-dimensional laminar flow based on Darcy model and two-energy equation for fluid and solid phases. A parametric study is carried out by varying the following parameters: cross flow to jet flow velocity ratio parameter (0 ⩽ M ⩽ 1); porosity scaled thermal conductivity ratio parameter (0.1 ⩽ K_r ⩽ 1000); heat transfer coefficient parameter (0.1 ⩽ H ⩽ 1000); Péclet number (1 ⩽ Pe ⩽ 1000) and Rayleigh number (10 ⩽ Ra ⩽ 100). The total average Nusselt number is defined based on the overall thermal conductivity, which is assumed to be the arithmetic mean of the porosity scaled thermal conductivity of the fluid and solid phases. The total average Nusselt number as well as the average Nusselt number for both fluid and solid phases is presented for different governing parameters. It is found that the presence of a weak cross flow in a jet impinging jet may degrade the heat transfer. The results show that the average Nusselt number calculated from the thermal equilibrium model are the maximum possible values and these values can be reproduced by large values of H × K_r.     

Heatlines based natural convection analysis in tilted isosceles triangular enclosures with linearly heated inclined walls: effect of various orientations

Natural convection in isosceles triangular enclosures with various configurations (case 1 — inverted, case 2 — straight and case 3 — tilted) is studied via heatline analysis for linear heating of inclined walls. Detailed analysis and comparison for various base angles (φ = 45°, 60°) of triangular enclosures have been carried out for a range of fluids (Pr = 0.015 − 1000) within Ra = 10³ − 10⁵ using Galerkin finite element method. The heat flow distributions indicate conduction dominant heat transfer at low Ra (Ra = 10³) for case 1 and case 2 whereas in case 3, convective heat flow is observed due to high buoyancy force. As Ra increases, enhanced thermal mixing is observed at the core of the cavity. Wall to wall heat transfer occurs at walls AB and AC due to linear heating boundary condition in all the cases. Although the distributions of fluid flow and heat flow are qualitatively similar for φ = 45° and 60°, the intensity of fluid flow and heat flow decreases as φ increases. Strength of fluid flow and heat flow circulation cells is found to be higher in case 3 for identical parameters. Results show that upper side wall (AC) for case 3 exhibits higher heat transfer rates whereas heat transfer rates for walls AB and AC are the same for case 1 and case 2. Also Nu_AB is higher for case 2 followed by case 1 and case 3 at the middle portion of wall AB. Thus to achieve high heat transfer from fluid to wall at the central region, case 2 and case 3 configurations may be recommended at high Ra (Ra = 10⁵) and Pr, irrespective of φ.     

Magnetoconvection in an enclosure of water near its density maximum with Soret and Dufour effects

In this paper, unsteady double-diffusive magnetoconvection of water in an enclosure with Soret and Dufour effects around the density maximum has been numerically investigated. The right vertical wall has constant temperature, θ_c, while left vertical wall is θ_h, with θ_h > θ_c. The concentration in right wall is maintained lower than left wall (c_h > c_c). The remaining horizontal walls are adiabatic. The governing equations are solved by control volume method using SIMPLE algorithm with QUICK scheme. Representative results illustrating the effects of the thermal Rayleigh number, Hartmann number, the direction of magnetic field, density inversion parameter, buoyancy ratio, Schmidt number, and Soret and Dufour parameters on the contour maps of the fluid flow, temperature and concentration as well as the profile of velocity at mid-section of the enclosure are reported. In addition, numerical results for the average Nusselt and Sherwood numbers are presented for various parametric conditions and discussed.     

Heat transfer in rotating scale-roughened trapezoidal duct at high rotation numbers

An experimental study of heat transfer in a radially rotating trapezoidal duct with two bevel walls roughened by deepened scales is performed with cooling applications to gas turbine rotor blades. Laboratory scale heat transfer data along the centerlines of two scale-roughened walls is generated within the parametric ranges of 7500 ? Re ? 15,000, 0 ? Ro ? 1.8 and 0.13 ? Δρ/ρ ? 0.42. No previous study has examined the heat transfer in a rotating scale-roughened channel and the present Ro range extends considerably from other researches to date. A selection of experimental data illustrates the individual and interactive impacts of Re, Ro and buoyancy number (Bu) on local heat transfer through which the manners of isolated and synergetic influences of Coriolis force and rotating buoyancy on heat transfer are examined. Local Nusselt number ratios between rotating and static channels on the stable (leading) and unstable (trailing) scale-roughened surfaces with Ro varying from 0.1 to 1.8 fall in the ranges of 0.8–2 and 1.1–2.5, respectively. Rotating buoyancy effects are weakened as Ro increases that impair local heat transfer for the present test configurations. Empirical heat transfer correlations for developed flow regions on two scale-roughened surfaces are derived that permit the evaluation of interactive and individual effects of Re, Ro and Bu on local heat transfer.     

History recovery and source identification of multiple gaseous contaminants releasing with thermal effects in an indoor environment

Thermal transport and transient dispersion of pollutants emitted from two discrete strips within the displacement ventilation enclosure have been modeled numerically. Following the full numerical simulation of turbulent air flows, the inverse determinations of multiple pollutant sources were conducted by the use of quasi reversibility methodology. Direct simulation together with the turbulent streamlines and turbulent heatlines demonstrate that the enclosure flow pattern, enclosure air thermal level and heat transfer potential will depend on the interactions of external forced flow and thermal buoyancy driven flows, i.e., Reynolds number (2 × 10³ ? Re ? 10⁴) and Grashof number (10⁶ ? Gr ? 10¹⁰). In subsequent forward time and backward time modeling of airborne pollutant transports, temporal evolutions of enclosure average concentration and pollutant exhaust are shown to depend on the supplying velocity (Re), thermal plume (Gr), pollutant diffusivity (0.1 ? Sc ? 2), and the pitch between both sources (0.2H ? d_PSL = d_PSR ? 0.7H). Reverse time modeling of airborne spread has demonstrated that increasing the spread rate and the concentration sensitivity of airborne pollutants will facilitate the identification of pollutant sources.     

Direct numerical simulation of wall-normal rotating turbulent channel flow with heat transfer

Direct numerical simulation of wall-normal rotating channel flow with heat transfer has been performed for the rotation number N_τ from 0 to 0.1, the Reynolds number 194 based on the friction velocity of non-rotating case and the half-height of the channel, and the Prandtl number 1. The objective of this study is to reveal the effects of rotation on the characteristics of turbulence and heat transfer. Some statistical turbulence and heat transfer quantities, including the mean velocity, temperature and their fluctuations, turbulent heat fluxes, and turbulence structures, are investigated. Based on the present calculated results, two typical rotation regimes are identified. When 0 < N_τ < 0.06, the turbulence statistics correlated with the spanwise velocity fluctuation are enhanced since the shear rate of spanwise mean flow induced by Coriolis force increases; however, the other statistics are suppressed. When N_τ > 0.06, all the turbulence statistics are suppressed significantly. To elucidate the effects of rotation on the turbulent heat transfer, the budget terms in the transport equation of turbulent heat fluxes are analyzed. Remarkable change of the direction of near-wall streak structures of the velocity and temperature fluctuations, nearly in alignment with the absolute mean flow direction, is revealed. An attempt to evaluate the mean spacing and the direction of streaky structures near the wall has been examined based on the two-point correlations of the velocity and temperature fluctuations.     

Experimental study of thermal behaviors in a rectangular channel with baffle of pores

This experimental study attempts to explore the local heat transfer in rectangular channel with baffles, and analyzes the experimental results of baffles with different heights and pores in the event of five Reynolds numbers and three heating quantities. Apart from increasing the perturbation of flow field, the channel's flow field with baffles, which is similar to a backward-facing step flow field, is very helpful to heat transfer. To obtain an optimized baffle and increase the perturbation of flow field, this experiment employed baffles with five heights (H = 10–50 mm) and different numbers of pores (N = 1–3), as well as heat flux: Q = 40–100 l/min, Reynolds number: 702–1752, and heating quantity: q_in = 90–750 W/m². In addition to measurement of overall temperature distribution, emphasis is also placed on analysis of local heat transfer coefficient. Furthermore, heat transfer distribution of channel can be applied to explain how the baffles of pores have an influence upon backward-facing step flow field, shear layer, recirculation region, reattachment region and redeveloped boundary layer. Finally, some empirical formulas derived form experimental results may provide a reference for future design.     

Velocity field and turbulence effects on heat transfer characteristics from surfaces with V-shaped ribs

The flow field of smooth surfaces and surfaces with V-shaped ribs (V-SR) was studied experimentally with a Laser-Doppler Anemometry (LDA) system. In addition, heat transfer characteristics were experimentally investigated. Heat transfer results from these surfaces under impingement of a circular jet array (5 × 3) using an infrared thermal imaging technique are presented.The velocity profiles were measured at Reynolds number of 10,000 and at H/d equal to 3 and 12. For each H/d position, profiles were collected from x/d = 0 to 6 axial locations. The heat transfer data were obtained at Reynolds numbers equal to 2000, 6000, and 10,000. Along the target plate, different boundary layer profiles were obtained for smooth and V-SR plates at H/d = 3 and 12. Positions of maximum radial and axial velocities and turbulence intensities have been determined for smooth and V-SR plates. For low jet-to-plate spacings, the production of turbulence kinetic energy is higher for the V-SR surfaces as compared to smooth surfaces. For H/d = 3, the radial velocities are higher for the V-SR surfaces as compared to smooth surfaces but for H/d = 12, the radial velocities are not nearly changed all x/d locations. The heat transfer results have also been compared with those of a smooth surface under the same flow conditions to determine the enhancement in the heat transfer coefficient from x/d = 0 to 3 locations. In these locations, the Nusselt numbers are higher for the V-SR surfaces as compared to smooth surfaces. The locations of the peaks and the minima are influenced by cross flow velocities which in turn depend on jet-to-plate spacing and V-SR arrangements. For all results, the Nusselt numbers at the stagnation points decrease with increase in H/d.     

Effect of optical properties on oscillatory hydromagnetic double-diffusive convection within semitransparent fluid

The effect of radiative heat transfer on the hydromagnetic double-diffusive convection in two-dimensional rectangular enclosure is studied numerically for fixed Prandtl, Rayleigh, and Lewis numbers, Pr = 13.6, Ra = 10⁵, Le = 2. Uniform temperatures and concentrations are imposed along the vertical walls while the horizontal walls are assumed to be adiabatic and impermeable to mass transfer. The influences of the optical thickness and scattering albedo of the semitransparent fluid on heat and mass transfer with and without magnetic damping are depicted. When progressively varying the optical thickness, multiple solutions are obtained which are steady or oscillatory accordingly to the initial conditions. the mechanisms of the transitions between steady compositionally dominated flow and unsteady thermally dominated flow are analyzed.     

Effect of fin thickness on the air-side performance of wavy fin-and-tube heat exchangers under dehumidifying conditions

This study investigated the effect of fin thickness on the air-side performance of wavy fin-and-tube heat exchangers under dehumidifying conditions. A total of 10 samples were tested with associated fin thickness (δ_f) of 0.115 mm and 0.25 mm, respectively. For a heat exchanger with two rows (N = 2) and fin pitch F_p of 1.41 mm, the effect of fin thickness on the heat transfer coefficient is more pronounced. The heat transfer coefficients for δ_f = 0.25 mm is about 5–50% higher than those for δ_f = 0.115 mm whereas the pressure drop for δ_f = 0.25 mm is about 5–20% higher. The unexpected difference in heat transfer coefficient subject to fin thickness is attributable to better interactions between the directed main flow and the swirled flow caused by the condensate droplet for δ_f = 0.25 mm. The maximum difference in heat transfer coefficients for N = 2 and F_p = 2.54 mm subject to the influence of fin thickness is reduced to about 20%, and there is no difference in heat transfer coefficient when the frontal velocity is above 3 m/s. For N ⩾ 4 and F_p = 2.54 mm, the influence of fin thickness on the heat transfer coefficients diminishes considerably. This is because of the presence of tube row, and the unsteady/vortex shedding feature at the down stream of wavy channel. Based on the present test results, a correlation is proposed to describe the air-side performance for wavy fin configurations, the mean deviations of the proposed heat transfer and friction correlations are 7.9% and 7.7%, respectively.     

Numerical research of nature convective heat transfer enhancement filled with nanofluids in rectangular enclosures

Heat transfer enhancement utilizing nanofluids in a two-dimensional enclosure is investigated for various pertinent parameters. The Khanafer's model is used to analyze heat transfer performance of nanofluids inside an enclosure taking into account the solid particle dispersion. Transport equations are model by a stream function-vorticity formulation and are solved numerically by finite-difference approach. Based upon the numerical predictions, the effects of Rayleigh number (Ra) and aspect ratio (AR) on the flow pattern and energy transport within the thermal boundary layer are presented. The diameter of the nanoparticle d_p is taken as 10 nm in nanofluids. The buoyancy parameter is 10³ ≤ Ra ≤ 10⁶ and aspect ratios (AR) of two-dimensional enclosure are 1/2, 1, 2. Results show that increasing the buoyancy parameter and volume fraction of nanofluids cause an increase in the average heat transfer coefficient. Finally, the empirical equation was built between average Nusselt number and volume fraction.     

Enhancement of heat transfer in turbulent channel flow over dimpled surface

A systematic numerical investigation of heat transfer in turbulent channel flow over dimpled surface is conducted. Both symmetric (or spherical) and asymmetric dimple with different depth ratios (h/D) and skewness (Dx and Dz) are considered for a series of Reynolds numbers Re_2H (based on bulk velocity and full channel height) between 4000 and 6000 while Prandtl number Pr is fixed at 0.7. It is found that the optimum dimple configuration for enhancing heat transfer measured in terms of the volume goodness factor is obtained for the case of asymmetric dimple with a depth ratio of h/D = 15% and stream-wise skewness of Dx = 15%. The heat transfer capacity in terms of Nusselt number is significantly increased, while the associated pressure loss is kept almost to the same level as the symmetric dimple with the same depth ratio. The present study also suggests that the heat transfer enhancement is closely related to ejection with counter-rotating flow, intensified secondary flow and vortex structures at the downstream rim of asymmetric dimple. All these findings suggest that a carefully designed asymmetric dimpled surface presents a viable means of enhancing heat transfer compared to the symmetric dimple.     

Convective diffusion from strip-like micro-probes into colloidal suspensions

Full equations of convective diffusion are solved numerically for a strip-like (2D) electrodiffusion friction probe in a stream of microdisperse liquid, assuming a non-linear near-to-wall velocity profile ranging from simple shear flow (p = 1) to ideal slip (p = 0). The range of generalized Peclet number H from H = 0.01 (almost pure spatial diffusion) to H = 100 (diffusion layer with negligible longitudinal diffusion) covers all cases of possible experimental relevance. The main result is expressed as a relative deviation of actual total diffusion flux N from its diffusion-layer approximation N_DLA, Ψ = N/N_DLA ? 1.     

Mixed convection flow and heat transfer across a square cylinder under the influence of aiding buoyancy at low Reynolds numbers

In this paper, mixed convection flow and heat transfer around a long cylinder of square cross-section under the influence of aiding buoyancy are investigated in the vertical unconfined configuration (Reynolds number, Re = 1–40 and Richardson number, Ri = 0–1). The semi-explicit finite volume method implemented on the collocated grid arrangement is used to solve the governing equations along with the appropriate boundary conditions. The onset of flow separation occurs between Re = 1–2, between Re = 2–3 and between Re = 3–4 for Ri = 0, 0.5 and 1, respectively. The flow is found to be steady for the range of conditions studied here. The friction, pressure and total drag coefficients are found to increase with Richardson number, i.e., as the influence of aiding buoyancy increases drag coefficients increase at the constant value of the Reynolds number. The temperature field around the obstacle is presented by isotherm contours at the Prandtl number of 0.7 (air). The local and average Nusselt numbers are calculated to give a detailed study of heat transfer over each surface of the square cylinder and an overall heat transfer rate and it is found that heat transfer increases with increase in Reynolds number and/or Richardson number. The simple expressions for the wake length and average cylinder Nusselt number are obtained for the range of conditions covered in this work.