Nusselt number and friction factor in thermally stratified turbulent channel flow under Non-Oberbeck–Boussinesq conditions

In stably stratified turbulence, computations under Oberbeck–Boussinesq (OB) hypothesis of temperature-independent fluid properties may lead to inaccurate representation of the flow field and to wrong estimates of momentum/heat transfer coefficients. This is clearly assessed here comparing direct numerical simulations of stratified turbulence under OB conditions to simulations under NOB (Non-Oberbeck–Boussinesq) conditions of temperature-dependent fluid viscosity and thermal expansion coefficient. Compared to the OB case, NOB conditions may induce local flow relaminarization with significant variations (up to 30%) of heat and momentum transfer coefficients. Together with DNS results, we propose a phenomenological model (based on turbulent bursts) for heat transfer prediction in stratified turbulence under OB and NOB conditions. Implications of NOB assumptions on mixing efficiency (i.e. flux Richardson number Ri_f) and turbulent Prandtl number (Pr_t) are also discussed. These results are of specific importance in RANS modelling, where the condition Pr_t = 1 is usually assumed (Reynolds analogy). Although this assumption is valid in some situations (i.e. boundary layer, pipe flow) there is uncertainty about its validity for stably-stratified turbulence. We demonstrate that this assumption is inaccurate when NOB effects become significant.     

Effect of thermo‐solutal stratification on recirculation flow patterns in a backward‐facing step channel flow

This paper presents results on the combined effect of thermo‐solutal buoyancy forces on the recirculatory flow behavior in a horizontal channel with backward‐facing step and the ensuing impact on heat and mass transfer phenomena. The governing equations for double diffusive mixed convection are represented in velocity–vorticity form of momentum equations, velocity Poisson equations, energy and concentration equations. Galerkin's finite‐element method has been employed to solve the governing equations. Recirculatory flow fields with heat and mass transfer are simulated for opposing and aiding thermo‐solutal buoyancy forces by assuming suitable boundary conditions for energy and concentration equations. The effect of Richardson number (0.1?Ri?10) and buoyancy ratio (?10?N?10) on the recirculation bubble and Nusselt and Sherwood numbers are studied in detail. For Richardson number greater than unity, distinct variations in the gradients of Nusselt number and Sherwood number with buoyancy ratio are observed for flow regimes with opposing and aiding buoyancy forces. Copyright © 2009 John Wiley & Sons, Ltd.     

Computation of conjugate heat transfer in the turbulent mixed convection regime in a vertical channel with multiple heat sources

This paper presents the results of a comprehensive numerical study to analyze conjugate, turbulent mixed convection heat transfer from a vertical channel with four heat sources, uniformly flush-mounted to one of the channel walls. The results are presented to study the effect of various parameters like thermal conductivity of wall material (k _s), thermal conductivity of flush-mounted discrete heat source (k _c), Reynolds number of fluid flow (Re _s), modified Richardson number (Ri ⁺) and aspect ratio (AR) of the channel. The standard k-ε turbulence model, modified by including buoyancy effects with physical boundary conditions, i.e. without wall functions, has been used for the analysis. Semi-staggered, non-uniform grids are used to discretise the two dimensional governing equations, using finite volume method. A correlation, encompassing a wide range of parameters, is developed for the non-dimensional maximum temperature (T ^*) using the asymptotic computational fluid dynamics (ACFD) technique.     

Numerical study on heat/mass transfer by shear/thermal/solutal-convection in a cavity under counteracting buoyancies

This study provides physical insights into a lid-driven square cavity filled with a mixture of a solvent vapor and non-condensable gas, subjected to the vertically parallel thermal and solutal gradients. The top lid is maintained at constant speed while bottom lid and the other two walls are kept fixed. Zero heat and mass fluxes are imposed on the vertical side walls. The transport equations are solved numerically through a pressure-correction-based iterative algorithm (SIMPLE) with the QUICK scheme for convective terms. The diffusivities of heat and salt are assumed to be equal throughout this investigation. The essential details of flow, temperature and concentration fields are presented for the opposing buoyancy forces ratio (B < 0) with special attention being given for the values of parameters for which either the flow inside the cavity is operated by the mechanically induced convection; or the flow structure inside the cavity is akin to a single-diffusive thermal or solutal convection. The variations of average rates of heat and mass transfer are uniform with the Reynolds number, while the variations of these quantities against the solutal Richardson number (Ri _C ) and thermal Richardson number (Ri _T ) point to the existence of the local minimum/maximum. Finally, two linear relations between Ri _C and Ri _T at a constant sliding speed are proposed to identify the above points of high and low transport phenomena and justified with the exhibition of the flow structures inside the cavity.     

Experimental study of thermal mixing layer using variable temperature hot-wire anemometry

The buoyancy effects on the development of the thermal mixing layer downstream from a horizontal separating plate were studied by comparing stable and unstable counter-gradient configurations. In this study, the novel experimental technique called parameterizable constant temperature anemometer, proposed by Ndoye et al. (Meas Sci Technol 21(7):075401, 2010), was improved to make possible the simultaneous measurement of temperature and two velocity components with an x-wire probe. The buoyancy effects on the flow are discussed through the transport equations of turbulent kinetic energy and temperature variance. In view of the low Richardson numbers at stake (Ri _f < 0.03), the buoyancy forces appeared logically to be quantitatively negligible compared to the main driving forces, but such a low-energy forcing mechanism was in fact sufficient in unstable configurations to increase the shear stress and the expansion rate of the mixing layer significantly, both phenomena being associated with enhanced production of turbulence.     

Instability threshold of a negatively buoyant fountain

Experimental simulations were carried out to investigate the onset of instability in negatively buoyant fountains by injecting glycerin–water mixtures into silicon oil. The transition from a stable to an unstable fountain structure is primarily governed by the Richardson number, and to a lesser extent, Reynolds number, viscosity ratio, Weber number and vent geometry. Transition nominally occurs at a Ri = 1.0. For a fountain issuing from a cylindrical pipe, the major effect of the Reynolds number is in determining whether or not the fountain is laminar or turbulent. The Reynolds number effect can be largely accounted for by basing a corrected Richardson number on the root mean square of the mean velocity. Viscosity ratio deviating from unity has the effect of stabilizing the flow structure and thereby reducing the transition Richardson number. Similarly, interfacial tension stabilizes the flow pattern resulting in a trend of increasing transition Richardson number with increasing Weber number. The results are valid in rectangular vents if the Richardson number and Reynolds number are based on the hydraulic diameter.     

Numerical prediction of turbulent boundary layer noise from a sharp-edged flat plate

An efficient hybrid uncorrelated wall plane waves–boundary element method (UWPW-BEM) technique is proposed to predict the flow-induced noise from a structure in low Mach number turbulent flow. Reynolds-averaged Navier-Stokes equations are used to estimate the turbulent boundary layer parameters such as convective velocity, boundary layer thickness, and wall shear stress over the surface of the structure. The spectrum of the wall pressure fluctuations is evaluated from the turbulent boundary layer parameters and by using semi-empirical models from literature. The wall pressure field underneath the turbulent boundary layer is synthesized by realizations of uncorrelated wall plane waves (UWPW). An acoustic BEM solver is then employed to compute the acoustic pressure scattered by the structure from the synthesized wall pressure field. Finally, the acoustic response of the structure in turbulent flow is obtained as an ensemble average of the acoustic pressures due to all realizations of uncorrelated plane waves. To demonstrate the hybrid UWPW-BEM approach, the self-noise generated by a flat plate in turbulent flow with Reynolds number based on chord Re_c = 4.9 × 10⁵ is predicted. The results are compared with those obtained from a large eddy simulation (LES)-BEM technique as well as with experimental data from literature.     

Heat transfer in turbulent pipe flow based on a new mixing length model

A new model for the heat transfer in turbulent pipe flow is presented based on a modified form of the mixing length theory developed by Cebeci [1] for boundary layer flow problems. The model predicts the velocity and temperature distributions and the Nusselt number for fluids with low, medium and high Prandtl numbers (Pr=.02 to 15) and fits the available experimental data very accurately for values of Reynolds number exceeding 10⁴. Expressions for the eddy conductivity and for the turbulent Prandtl number are presented and shown to be dependent upon the Reynolds number, the Prandtl number, and the distance from the tube wall.     

壁面加热湍流相干结构的子波分析实验研究

通过在平板湍流边界层沿流向固壁表面平行放置若干条通电加热的金属细丝,在平板表面形成沿展向周期性分布的温度场,利用该温度场引起的空气热对流,在湍流边界层近壁区域产生一组沿湍流边界层展向周期分布的大尺度流向涡结构,改变了平板湍流边界层中不同尺度结构及其能量分布。采用对壁湍流多尺度结构的子波分析表明,在湍流边界层近壁区域产生规则的流向涡结构将壁湍流各种尺度湍涡结构不规则的脉动有序地组织起来,抑制了壁湍流各种尺度湍涡结构脉动,特别抑制了能量最大尺度湍涡结构的脉动,减小由于湍流脉动引起的在湍流边界层法向和展向的动量和能量损耗,从而减小了湍流的阻力。     

Large eddy simulation of a developing turbulent boundary layer at a low Reynolds number

A spectral code has been used to simulate a developing turbulent boundary layer at low Reynolds number Re_θ (based on free stream velocity and momentum thickness) between 353 and 576. The starting field was generated by allowing a step change of temperature to diffuse outwards from one wall in a fully developed channel flow. The thermal boundary layer so created was conditionally sampled to convert it into a momentum boundary layer with an irrotational free stream region, a process which is justified by appeal to experiments. This initial field was allowed to develop until the momentum boundary layer thickness δ₉₉₅ had grown to about 1·5 times its original thickness. The results of the simulation have been compared with a wide range of experimental data. The outcome of this comparison is generally very satisfactory; the main trends of the experiments are well reproduced and our simulation supplements and extends the existing sets of experimental data. The simulation also gives pressure statistics which cannot be obtained experimentally. In particular, it gives the contribution of pressure diffusion to the balance equations for the Reynolds stress and indicates the error produced by omitting this term.     

The influence of the periodic disturbance on the local heat transfer in separated and reattached flow

A numerical study based on the large eddy simulation methodology was made of heat transfer in locally disturbed turbulent separated and reattached flow over a backward facing step. The local disturbance was given to the flow by a sinusoidally blowing/suction of the fluid into a separated shear layer. The Reynolds number was fixed at 33,000 and Richardson number at 0.5. The disturbance frequency was varied in the range 0 ≤ St ≤ 2, where St is the Strouhal number of disturbance. The obtained results revealed the existence of an optimum perturbation frequency value, St = 0.25, in terms of the reduced reattachment length. At this frequency the heat transfer is significantly enhanced in the recirculation zone. The influence of the frequency and the amplitude of disturbance, in the maximum heat transfer positions and the maximum local Nusselt number, is analysed.     

Mixed convection heat transfer from three heated square cylinders in cross-flow at low Reynolds numbers

This paper presents the effects of cross buoyancy and Prandtl number on the flow and heat transfer characteristics around three equal isothermal square cylinders arranged in a staggered configuration within an unconfined medium. Transient two-dimensional numerical simulations are performed with a finite volume code based on the SIMPLEC algorithm in a collocated grid system. The pertinent dimensionless parameters, such as Reynolds, Prandtl and Richardson numbers are considered in the range of 1 ≤ Re ≤ 30, 0.7 ≤ Pr ≤ 100 and 0 ≤ Ri ≤ 1. The representative streamlines, vortex structures and isotherm patterns are presented and discussed. In addition, the overall drag and lift coefficients and average Nusselt numbers are determined to elucidate the effects of Reynolds, Prandtl and Richardson numbers on flow and heat transfer. The flow is observed to be steady for all the ranges of parameters considered. The drag coefficient is found to decrease with Re (for Ri = 0) and Ri at low Pr, whereas it increases with Pr at higher Ri. The lift coefficient decreases with Ri at low Pr and increases with Pr at higher Ri. The time and surface average cylinder Nusselt number is found to increase monotonically with Re as well as Pr while it remains almost insensitive to Ri at low Pr.     

Analyzing the influence of compressibility on the rapid pressure–strain rate correlation in turbulent shear flows

The influence of compressibility on the rapid pressure–strain rate tensor is investigated using the Green’s function for the wave equation governing pressure fluctuations in compressible homogeneous shear flow. The solution for the Green’s function is obtained as a combination of parabolic cylinder functions; it is oscillatory with monotonically increasing frequency and decreasing amplitude at large times, and anisotropic in wave-vector space. The Green’s function depends explicitly on the turbulent Mach number M _t , given by the root mean square turbulent velocity fluctuations divided by the speed of sound, and the gradient Mach number M _g , which is the mean shear rate times the transverse integral scale of the turbulence divided by the speed of sound. Assuming a form for the temporal decorrelation of velocity fluctuations brought about by the turbulence, the rapid pressure–strain rate tensor is expressed exactly in terms of the energy (or Reynolds stress) spectrum tensor and the time integral of the Green’s function times a decaying exponential. A model for the energy spectrum tensor linear in Reynolds stress anisotropies and in mean shear is assumed for closure. The expression for the rapid pressure–strain correlation is evaluated using parameters applicable to a mixing layer and a boundary layer. It is found that for the same range of M _t there is a large reduction of the pressure–strain correlation in the mixing layer but not in the boundary layer. Implications for compressible turbulence modeling are also explored.      

Simulation of a Large-Eddy-Break-up Device (LEBU) in a Moderate Reynolds Number Turbulent Boundary Layer

A well-resolved large eddy simulation (LES) of a large-eddy break-up (LEBU) device in a spatially evolving turbulent boundary layer is performed with, Reynolds number, based on free-stream velocity and momentum-loss thickness, of R e _θ ≈ 4300. The implementation of the LEBU is via an immersed boundary method. The LEBU is positioned at a wall-normal distance of 0.8 δ (δ denoting the local boundary layer thickness at the location of the LEBU) from the wall. The LEBU acts to delay the growth of the turbulent boundary layer and produces global skin friction reduction beyond 180δ downstream of the LEBU, with a peak local skin friction reduction of approximately 12 %. However, no net drag reduction is found when accounting for the device drag of the LEBU in accordance with the towing tank experiments by Sahlin et al. (Phys. Fluids 31, 2814, 1988). Further investigation is performed on the interactions of high and low momentum bulges with the LEBU and the corresponding output is analysed, showing a ‘break-up’ of these large momentum bulges downstream of the LEBU. In addition, results from the spanwise energy spectra show consistent reduction in energy at spanwise length scales for \(\lambda _{z}^{+} > 1000\) independent of streamwise and wall-normal location when compared to the corresponding turbulent boundary layer without LEBU.     

Turbulent diffusion of a buoyant layer at a wall

Summary When a light fluid is injected at a steady rate at the roof of a tunnel in which there is a turbulent main flow of a heavier fluid, the turbulent diffusion of the light layer may be considerably reduced due to buoyancy. For large Richardson numbers turbulent mixing ceases altogether.The equations of motion and diffusion were solved by introducing an eddy diffusivity which is dependant on the Richardson number. Experiments were made on brine (floor) layers in a water flow, and on methane (roof) layers in an air flow. Results were essentially in agreement with theory.The motion and mixing of the layers depend mainly on the inclination of the tunnel and on a dimensionless combination of main-flow velocity, gravity, relative density difference, volume input rate of layer fluid, and tunnel width. Values of the dimensionless parameter are suggested to overcome the effects of buoyancy on mixing, and to prevent layers from moving up a slope against the main flow.     

Measurement of interfacial distortions at a stratified entrainment interface

A laser beam scanning technique was used to measure the rms amplitude of the largest wave (interfacial distortion) that is present at the entrainment interface of a stratified fluid. Both linearly-stratified and two-layer fluid systems, subjected to entrainment by oscillating-grid induced shear-free turbulence, are considered. The measurements were found to be in general agreement with a theory due to Long (1978), which predicts /D _e Ri ^–3/4 , where D _e is the mean thickness of the mixed layer and Ri is the appropriately defined bulk Richardson number at the interface.     

Opposing mixed convection and its interaction with radiation inside eccentric horizontal cylindrical annulus

In the present investigation, the coupled phenomenon of opposing mixed convection and radiation within differentially heated eccentric horizontal cylindrical annulus has been numerically simulated. The radiation transfer contributed from the participating medium is obtained by solving the nonlinear integro‐differential radiative transfer equation using discrete ordinate method. The participating gray medium is considered to be emitting, absorbing and isotropically scattering. The walls of the annulus are considered to be opaque, diffuse and gray. In the study it has been observed that the Richardson number ‘Ri’ has a small effect on the total Nusselt number ‘Nu’ in mixed convection heat transfer with or without radiation. From the present investigation it is found that substantial changes occur in isotherms as well as in flow patterns, when the Richardson number is allowed to vary in the range of 0.01–1. The influence of radiative parameters on the interaction phenomenon has been delineated through isotherm and streamline pattern. Copyright © 2008 John Wiley & Sons, Ltd.     

The variation of heat transfer in a two‐sided lid‐driven differentially heated square cavity with nanofluids containing carbon nanotubes for physical properties of fluid dependent on temperature

The behavior of nanofluids containing cylindrical nanoparticles are investigated numerically inside a two‐sided lid‐driven differentially heated square cavity to gain insight into the convective recirculation and flow processes induced by a nanofluid. The physical properties of the base fluid such as viscosity, thermal conductivity and thermal expansion coefficient are, respectively, assumed to be temperature independent (taking the mean temperature of the left and right walls) and temperature dependent. A model is developed to analyze the behavior of nanofluids taking into account the nanoparticle volume fraction whereas the transport equations are solved numerically with finite volume approach using SIMPLEC algorithm. The left and right moving walls are maintained at different constant temperatures while the upper and bottom walls are thermally insulated. The directions of the moving walls were considered in a way that the force and natural convections aid each other. The governing parameter Richardson number was 0.1<Ri<50.0 but due to space constraints only the results for 0.1<Ri<10.0 from fluid flow are presented. It was found that the temperature dependency of physical properties at different Richardson numbers and nanoparticle volume fractions affects the fluid flow and heat transfer in the cavities. Finally, comparisons between the behaviors of the average Nusselt number at the left wall for two cases are presented. Copyright © 2011 John Wiley & Sons, Ltd.     

A von Neumann–Smagorinsky turbulent transport model for stratified shear flows

A simple subgrid turbulent diffusion model based on an analogy to the von Neumann–Richtmyer artificial viscosity is explored for use in modelling mixing in turbulent stratified shear flow. The model may be more generally applicable to multicomponent turbulent hydrodynamics and to subgrid turbulent transport of momentum, composition and energy. As in the case of the von Neumann artificial viscosity and many subgrid-scale models for large-eddy simulation, the turbulent diffusivity explicitly depends on the grid size and is not based on a quantitative model of the unresolved turbulence. In order to address the issue that it is often not known a priori when and where a flow will become turbulent, the turbulent diffusivity is set to zero when the flow is expected to be stable on the basis of a Richardson/Rayleigh–Taylor stability criterion, in analogy to setting the von Neumann artificial viscosity to zero in expanding flows. One-dimensional predictions of this model applied to a simple shear flow configuration are compared to those obtained using a K–ε model. The density and velocity profiles predicted by both models are shown to be very similar.     

Laboratory experiments of a jet impinging on a stratified interface

 The entrainment rates of vertical and inclined jets impinging on a stratified interface are measured in water tank experiments. At moderate Richardson number, the entrainment rate of the vertical jet is proportional to Ri ^-1/2, independent of Reynolds number. The inclined jets are tilted at 15° from the vertical. In one case, the jet nozzle is rotated about a vertical axis, so that the inclined jet precesses, while in the other, it is stationary. The inclined jets entrain at a rate proportional to Ri ^-3/2, whether precessing or not. This behavior is consistent with a new model of stratified entrainment which accounts for vortex persistence. Received: 15 October 1996/Accepted: 19 December 1996