Effect of particles on turbulent thermal field of channel flow with different Prandtl numbers

The direct numerical simulation(DNS) of heat transfer in a fully developed non-isothermal particle-laden turbulent channel flow is performed.The focus of this paper is on the modulation of the particles on turbulent thermal statistics in the particle-laden flow with three Prandtl numbers(P r = 0.71,1.5,and 3.0) and a shear Reynolds number(Reτ = 180).Some typical thermal statistics,including normalized mean temperature and their fluctuations,turbulent heat fluxes,Nusselt number and so on,are analyzed.The results show that the particles have less effects on turbulent thermal fields with the increase of Prandtl number.Two reasons can explain this.First,the correlation between fluid thermal field and velocity field decreases as the Prandtl number increases,and the modulation of turbulent velocity field induced by the particles has less influence on the turbulent thermal field.Second,the heat exchange between turbulence and particles decreases for the particle-laden flow with the larger Prandtl number,and the thermal feedback of the particles to turbulence becomes weak.     

DNS and Modelling of Passive Scalar Transport in Turbulent Channel Flow with a Focus on Scalar Dissipation Rate Modelling

A direct numerical simulation of turbulent channel flow with an imposed mean scalar gradient is analyzed with a focus on passive scalar flux modelling and in particular the treatment of the passive scalar dissipation equation. The Prandtl number is 0.71 and the Reynolds number based on the wall friction velocity and the channel half width is 265. Budgets are presented for the passive scalar variance and its dissipation rate, as well as for the individual scalar flux components. These form a basis for a discussion of modelling issues related to explicit algebraic scalar flux modelling. This revised version was published online in July 2006 with corrections to the Cover Date.     

Turbulence Budgets in Buoyancy-affected Vertical Backward-facing Step Flow at Low Prandtl Number

The paper investigates buoyancy impact on the vertical flow over a backward-facing step at low Prandtl number by Direct Numerical Simulation. In particular, the very low Prandtl number of liquid sodium, 0.0088, is considered in the regime of mixed convection, i.e. for Richardson numbers below unity. The effects of buoyancy on mean flow, heat transfer and turbulence are assessed. Buoyancy is found to attenuate recirculation and, consequently, increase heat transfer. Turbulence is decreased in the attached boundary layer for moderate buoyancy impact but surpasses the levels found in forced convection at the largest Richardson number investigated. Beyond the mean flow and second moments, the budgets of turbulent kinetic energy, Reynolds shear stress, temperature variance, and turbulent heat flux components are studied and related to the alterations in the mean field quantities. Due to scale separation, production and dissipation nearly balance for temperature variance while this is not the case for turbulent kinetic energy. Similar findings for the turbulent heat fluxes show that the correlation between temperature and pressure gradient is the most important contribution to the budget aside from production and dissipation. In addition to the physical insight into this flow, the data presented may be used for the validation and improvement of turbulence models for liquid metal flows.     

An application of the k-ɛ model with variable Prandtl number to heat transfer computations in air flows

The two-equation `low Reynolds number' k-? model of turbulence with a set of universal constants suggested by Launder and Sharma is modified in the present paper. The variability of the turbulent Prandtl number Pr_t in the energy equation is assumed along with a change of a constant in the dissipation term of the turbulent kinetic energy equation. The turbulent heat transfer is computed for an air flow in a circular pipe for the Reynolds number within the range of 10⁴?4. The modification considerably improves the agreement between the numerical results and the experiment data published in the available literature.     

Measurement and decomposition of periodic flow structures downstream of a test turbine

The temperature dissipation rate inferred from the balance of $\overline{\theta^{2}}/2$ budget is used for the purpose of studying different methods employed to directly measure dissipation. The terms involved in the budget equation of temperature variance are measured with laser Doppler velocimetry and cold-wire thermometry used simultaneously. This study focuses on the centerline of a turbulent round jet, in the far field, at high Reynolds number (x/D = 30, Re _D  = 1.5 × 10⁵ and Re _λ  = 548). Particular attention is devoted to statistical convergence of second- and third-order moments of velocity and temperature fluctuations. Temperature dissipation obtained by Taylor’s hypothesis and radial temperature derivative spectra confirm local isotropy. A high level of low wave number content is reported for the longitudinal derivative spectra, probably due to transverse mode spectral aliasing and noise contamination for small wire separation. A parallel is drawn between finite difference formulations and the behavior of the autocorrelation coefficient for small wire separations. The temperature dissipation estimates found are close to the budget reference value, but spectral analysis cast doubts on the validity of the streamwise derivative obtained with a pair of probes.     

The influence of pipe length on thermal statistics computed from DNS of turbulent heat transfer

We present results from direct numerical simulation of turbulent heat transfer in pipe flow at a bulk flow Reynolds number of 5000 and Prandtl numbers ranging from 0.025 to 2.0 in order to examine the effect of streamwise pipe length (πδ ≡ πD/2 ? L ? 12πδ) on the convergence of thermal turbulence statistics. Various lower and higher order thermal statistics such as mean temperature, rms of fluctuating temperature, turbulent heat fluxes, two-point auto and cross-correlations, skewness and flatness were computed and it is found that the value of L required for convergence of the statistics depends on the Prandtl number: larger Prandtl numbers requires comparatively shorter pipe length for convergence of most of the thermal statistics.     

The influence of pipe length on thermal statistics computed from DNS of turbulent heat transfer

We present results from direct numerical simulation of turbulent heat transfer in pipe flow at a bulk flow Reynolds number of 5000 and Prandtl numbers ranging from 0.025 to 2.0 in order to examine the effect of streamwise pipe length (πδ ≡ πD/2 ⩽ L ⩽ 12πδ) on the convergence of thermal turbulence statistics. Various lower and higher order thermal statistics such as mean temperature, rms of fluctuating temperature, turbulent heat fluxes, two-point auto and cross-correlations, skewness and flatness were computed and it is found that the value of L required for convergence of the statistics depends on the Prandtl number: larger Prandtl numbers requires comparatively shorter pipe length for convergence of most of the thermal statistics.     

DNS of turbulent channel flow with conjugate heat transfer: Effect of thermal boundary conditions on the second moments and budgets

Budgets of turbulent heat fluxes and temperature variance obtained from the Direct Numerical Simulation of an incompressible periodic channel flow with a Reynolds number of 150 (based on friction velocity) and a Prandtl number of 0.71 are presented and analysed for four cases: locally imposed temperature at the wall (constant Dirichlet), locally imposed heat flux (constant Neumann), heat exchange coefficient (Robin) and 3D conjugate heat transfer. The dissipation rate associated with the temperature variance is strongly impacted by the thermal boundary condition. For non-conjugate cases, a straightforward analytical analysis establishes the connection between the boundary condition, the temperature variance and the wall-normal part of the thermal dissipation rate at the wall. For the conjugate case, the two-point correlations of the thermal field in the solid domain confirms the existence of very large scale thermal structures.     

Flow and heat transfer in a pipe with a fin attached to inner wall

Enhancement of heat transfer to the fluid can be done by turbulence promoters such as attached fins to the pipe walls. In this study, the flow field and the heat transfer rates were numerically investigated in a pipe with an internally attached fin. Numerical simulations were conducted for four different types of fluids and for different fin heights and locations, and as the Reynolds number was varied, the effects of the fin on Nusselt number and friction factors were investigated. For all the Reynolds numbers considered in this study, the effect of fin location on the heat transfer rate and friction factor was negligible. As the fin height was increased, the mean Nusselt number and the friction factor also increased in the turbulent flow regimes. For low Prandtl number fluids (Pr = 0.011), the main heat transfer mode is conduction, and hence the mean Nusselt number slightly affected the flow rates.     

Numerical simulation of heat transfer in a pipe with non-homogeneous thermal boundary conditions

Direct numerical simulations of heat transfer in a fully-developed turbulent pipe flow with circumferentially-varying thermal boundary conditions are reported. Three cases have been considered for friction Reynolds number in the range 180–360 and Prandtl number in the range 0.7–4. The temperature statistics under these heating conditions are characterized. Eddy diffusivities and turbulent Prandtl numbers for radial and circumferential directions are evaluated and compared to the values predicted by simple models. It is found that the usual assumptions made in these models provide reasonable predictions far from the wall and that corrections to the models are needed near the wall.     

Passive heat transfer in a turbulent channel flow simulation using large eddy simulation based on the lattice Boltzmann method framework

In this paper, a large eddy simulation based on the lattice Boltzmann framework is carried out to simulate the heat transfer in a turbulent channel flow, in which the temperature can be regarded as a passive scalar. A double multiple relaxation time (DMRT) thermal lattice Boltzmann model is employed. While applying DMRT, a multiple relaxation time D3Q19 model is used to simulate the flow field, and a multiple relaxation time D3Q7 model is used to simulate the temperature field. The dynamic subgrid stress model, in which the turbulent eddy viscosity and the turbulent Prandtl number are dynamically computed, is integrated to describe the subgrid effect. Not only the strain rate but also the temperature gradient is calculated locally by the non-equilibrium moments. The Reynolds number based on the shear velocity and channel half height is 180. The molecular Prandtl numbers are set to be 0.025 and 0.71. Statistical quantities, such as the average velocity, average temperature, Reynolds stress, root mean square (RMS) velocity fluctuations, RMS temperature and turbulent heat flux are obtained and compared with the available data. The results demonstrate great reliability of DMRT–LES in studying turbulence.     

LES-based filter-matrix lattice Boltzmann model for simulating turbulent natural convection in a square cavity

In this paper, a novel thermal filter-matrix lattice Boltzmann model based on large eddy simulation (LES) is proposed for simulating turbulent natural convection. In this study, the Vreman subgrid-scale eddy-viscosity model is introduced into the present framework of LES to accurately predict the flow in near-wall region. Two dimensional numerical simulations of natural convection in a square cavity were performed at high Rayleigh number varying from 10⁷ to 10¹⁰ with a fixed Prandtl number of Pr = 0.71. The influences of the higher-order terms upon the present results at high Rayleigh numbers are examined, taking Ra = 10⁷ and 10⁸ as the example, revealing that the proper minimization of the higher-order terms can improve numerical accuracy of present model for high Rayleigh convective flow. For the turbulent convective flow, the time-averaged quantities in the median lines are presented and compared against those available results from previous studies. The general structure of turbulent boundary layers is well predicted. All numerical results exhibit good agreement with the benchmark solutions available in the previous literatures.     

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.     

A comprehensive assessment of the Reynolds Analogy in predicting heat transfer in turbulent wall-bounded shear flows

Heat transfer modeling plays a major role in design and optimization of modern and efficient thermal-fluid systems. However, currently available models suffer from a fundamental shortcoming: their development is based on the general notion that accurate prediction of the flow field will guarantee an appropriate prediction of the thermal field, known as the Reynolds Analogy. This investigation presents a comprehensive assessment of the capability of the Reynolds Analogy in predicting turbulent heat transfer when applied to turbulent shear flows of fluids with different Prandtl numbers. It turns out that the Reynolds Analogy is able to provide acceptable results for first order statistics only when fluids with Prandtl number close to unity are considered. Further, it is shown that unsteady simulations could provide acceptable results on second order statistics concerning fluids with different Prandtl numbers, if appropriate grid design/resolution is provided that allows to resolve essential dynamics of the thermal field. However, accurate prediction of higher order statistics close to solid surface requires more advanced heat transfer models that can provide accurate information on thermal time scales, in case the grid is too coarse to support accurate resolving of the essential thermal dynamics in these regions.     

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.     

Passive heat transfer in a turbulent channel flow simulation using large eddy simulation based on the lattice Boltzmann method framework

In this paper, a large eddy simulation based on the lattice Boltzmann framework is carried out to simulate the heat transfer in a turbulent channel flow, in which the temperature can be regarded as a passive scalar. A double multiple relaxation time (DMRT) thermal lattice Boltzmann model is employed. While applying DMRT, a multiple relaxation time D3Q19 model is used to simulate the flow field, and a multiple relaxation time D3Q7 model is used to simulate the temperature field. The dynamic subgrid stress model, in which the turbulent eddy viscosity and the turbulent Prandtl number are dynamically computed, is integrated to describe the subgrid effect. Not only the strain rate but also the temperature gradient is calculated locally by the non-equilibrium moments. The Reynolds number based on the shear velocity and channel half height is 180. The molecular Prandtl numbers are set to be 0.025 and 0.71. Statistical quantities, such as the average velocity, average temperature, Reynolds stress, root mean square (RMS) velocity fluctuations, RMS temperature and turbulent heat flux are obtained and compared with the available data. The results demonstrate great reliability of DMRT–LES in studying turbulence.     

Investigation of transitional and turbulent heat and momentum transport in a rotating cavity

The paper gives the results of the DNS/LES which was performed to investigate the transitional and turbulent non-isothermal flows within a rotor/stator cavity. Computations were performed for the cavity of aspect ratio L = 2–35, Rm = 1.8 and for rotational Reynolds numbers up to 290000. The main purpose of the investigations was to analyze the influence of aspect ratio and Reynolds number on the flow structure and heat transfer. The numerical solution is based on a pseudo-spectral Chebyshev–Fourier–Galerkin collocation approximation. The time scheme is semi-implicit second-order accurate, which combines an implicit treatment of the diffusive terms and an explicit Adams–Bashforth extrapolation for the non-linear convective terms. In the paper we analyze distributions of the Reynolds stress tensor components, the turbulent heat flux tensor components, Nusselt number distributions and the turbulent Prandtl number and other structural parameters, which can be useful for modeling purposes. Selected results are compared with the experimental data obtained for single heated rotating disk by Elkins and Eaton (2000).     

A low Reynolds number variant of partially-averaged Navier-Stokes model for turbulence

A low Reynolds number (LRN) formulation based on the Partially Averaged Navier-Stokes (PANS) modelling method is presented, which incorporates improved asymptotic representation in near-wall turbulence modelling. The effect of near-wall viscous damping can thus be better accounted for in simulations of wall-bounded turbulent flows. The proposed LRN PANS model uses an LRN k-ε model as the base model and introduces directly its model functions into the PANS formulation. As a result, the inappropriate wall-limiting behavior inherent in the original PANS model is corrected. An interesting feature of the PANS model is that the turbulent Prandtl numbers in the k and ε equations are modified compared to the base model. It is found that this modification has a significant effect on the modelled turbulence. The proposed LRN PANS model is scrutinized in computations of decaying grid turbulence, turbulent channel flow and periodic hill flow, of which the latter has been computed at two different Reynolds numbers of Re = 10,600 and 37,000. In comparison with available DNS, LES or experimental data, the LRN PANS model produces improved predictions over the standard PANS model, particularly in the near-wall region and for resolved turbulence statistics. Furthermore, the LRN PANS model gives similar or better results - at a reduced CPU time - as compared to the Dynamic Smagorinsky model.     

Pr数对湍流被动标量输运的影响

采用直接数值模拟的方法,研究分子Pг数对湍流被动标量输运的影响,并提供充分的证据证明,湍流Pг数明显依赖于分子Pг数.在算例中,湍流雷诺平均PгT数与分子Pг数的倒数呈线性关系;湍流亚格子Pгt数与分子Pг数的关系较为复杂,在分子Pг数为1附近时,湍流亚格子Pгt数出现极小值.     

Eine asymptotische Analyse des Wärmeüberganges im Einlaufbereich von turbulenten Kanal- und Rohrströmungen

The basic equations for turbulent entrance flow are deduced from an asymptotic expansion of the Navier-Stokes equations and the thermal energy equation forRe→∞. Together with a turbulence model they can be solved numerically. Solutions are independent of the Reynolds and Prandtl number. Based on theses solutions, the skin friction and heat transfer as well as velocity and temperature profiles can be determined for finite Reynolds numbers and Prandtl numbersO (1).