Direct numerical simulation of turbulent heat transfer in annuli: Effect of heat flux ratio

Fully developed turbulent flow and heat transfer in a concentric annular duct is investigated for the first time by using a direct numerical simulation (DNS) with isoflux conditions imposed at both walls. The Reynolds number based on the half-width between inner and outer walls, δ=(r₂-r₁)/2, and the laminar maximum velocity is Re_δ=3500. A Prandtl number Pr=0.71 and a radius ratio r^*=0.1 were retained. The main objective of this work is to examine the effect of the heat flux density ratio, q^*=q₁/q₂, on different thermal statistics (mean temperature profiles, root mean square (rms) of temperature fluctuations, turbulent heat fluxes, heat transfer, etc.). To validate the present DNS calculations, predictions of the flow and thermal fields with q^*=1 are compared to results recently reported in the archival literature. A good agreement with available DNS data is shown. The effect of heat flux ratio q^* on turbulent thermal statistics in annular duct with arbitrarily prescribed heat flux is discussed then. This investigation highlights that heat flux ratio has a marked influence on the thermal field. When q^* varies from 0 to 0.01, the rms of temperature fluctuations and the turbulent heat fluxes are more intense near the outer wall while changes in q^* from 1 to 100, lead to opposite trends.     

Direct numerical simulation of compressible turbulent flows

This paper reviews the authors＇ recent studies on compressible turbulence by using direct numerical simulation （DNS）,including DNS of isotropic（decaying） turbulence, turbulent mixing-layer,turbulent boundary-layer and shock/boundary-layer interaction.Turbulence statistics, compressibility effects,turbulent kinetic energy budget and coherent structures are studied based on the DNS data.The mechanism of sound source in turbulent flows is also analyzed. It shows that DNS is a powerful tool for the mechanistic study of compressible turbulence.     

Direct numerical simulation of a turbulent plane Couette-Poiseuille flow with zero-mean shear

Direct numerical simulations of a turbulent Couette-Poiseuille flow with zero-mean-shear at the moving wall (SL-flow) is performed to examine flow features compared to those for a turbulent pure Poiseuille flow (P-flow). Profiles of the streamwise mean velocity, indicator function and ratio of production to dissipation show that the logarithmic region is significantly elongated for the SL-flow compared to that for the P-flow at a similar Reynolds number. In addition, the magnitudes of the Reynolds stresses are found to be larger in both inner and outer layers for the SL-flow than those for the P-flow. The spanwise spectra of the production term in the turbulent kinetic energy equation are examined to provide a structural basis for explaining the statistical behaviors. In addition, because the growth of the energy-containing motions extends to the outer layer further for the SL-flow due to the presence of a positive mean shear throughout the entire wall layer, the self-similar behavior of the energy balance between the production and transport terms with respect to the self-similar wavenumber is found far from the wall. We also find the increase in the number of uniform momentum zones in the SL-flow, revealing the hierarchical distribution of the energy-containing eddies which are composed of multiple uniform momentum zones. These coherent motions lead to the elongation of the logarithmic region for the SL-flow. Finally, investigation of the turbulent energy transfer process in a spectral domain for the SL-flow demonstrates importance of outer layer very-long structures, and these structures attribute to the energy transport process in an entire flow field.     

Direct numerical simulation of turbulent heat transfer in a square duct with transverse ribs mounted on one wall

Turbulent heat transfer in a ribbed square duct of three different blockage ratios are investigated using direct numerical simulation (DNS). The results of ribbed duct cases are compared with those of a heated smooth duct flow. It is observed that owing to the existence of the ribs and confinement of the duct, organized secondary flows appear as large streamwise-elongated vortices, which intensely interact with the rib elements and four sidewalls and have profound influences on the transport of momentum and thermal energy. This study also shows that the drag and heat transfer coefficients are highly sensitive to the rib height. It is observed that as the rib height increases, the impinging effect of the flow on the windward face of the rib strengthens, leading to enhanced rates of turbulent mixing and heat transfer. The influence of sidewalls and rib height on the turbulence structures associated with temperature fluctuations are analyzed based on multiple tools such as vortex swirling strengths, temporal auto-correlations, spatial two-point cross-correlations, joint probability density functions (JPDF) between the temperature and velocity fluctuations, statistical moments of different orders, and temperature spectra.     

Direct numerical simulation of stable and unstable turbulent thermal boundary layers

This paper presents direct numerical simulations (DNS) of stable and unstable turbulent thermal boundary layers. Since a buoyancy-affected boundary layer is often encountered in an urban environmental space where stable and unstable stratifications exist, exploring a buoyancy-affected boundary layer is very important to know the transport phenomena of the flow in an urban space. Although actual observation may qualitatively provide the characteristics of these flows, the relevant quantitative turbulent quantities are very difficult to measure. Thus, in order to quantitatively investigate a buoyancy-affected boundary layer in detail, we have here carried out for the first time time- and space-developing DNS of slightly stable and unstable turbulent thermal boundary layers. The DNS results show the quantitative turbulent statistics and structures of stable and unstable thermal boundary layers, in which the characteristic transport phenomena of thermally stratified boundary layers are demonstrated by indicating the budgets of turbulent shear stress and turbulent heat flux. Even though the input of buoyant force is not large, the influence of buoyancy is clearly revealed in both stable and unstable turbulent boundary layers. In particular, it is found that both stable and unstable thermal stratifications caused by the weak buoyant force remarkably alter the structure of near-wall turbulence.     

Direct numerical simulation of a turbulent Couette-Poiseuille flow,part 2: Large- and very-large-scale motions

A direct numerical simulation dataset of a fully developed turbulent Couette-Poiseuille flow is analyzed to investigate the spatial organization of streamwise velocity-fluctuating u-structures on large and very large scales. Instantaneous and statistical flow fields show that negative-u structures with a small scale on a stationary bottom wall grow throughout the centerline due to the continuous positive mean shear, and they penetrate to the opposite moving wall. The development of an initial vortical structure related to negative-u structures on the bottom wall into a large-scale hairpin vortex packet with new hairpin vortices, which are created upstream and close to the wall, is consistent with the auto-generation process in a Poiseuille flow (Zhou et al., J. Fluid Mech., vol. 387, 1999, pp. 353–396). Although the initial vortical structure associated with positive-u structures on the top wall also grows toward the bottom wall, the spatial development of the structure is less coherent with weak strength due to the reduced mean shear near the top wall, resulting in less turbulent energy on the top wall. The continuous growth of the structures from a wall to the opposite wall explains the enhanced wall-normal transport of the streamwise turbulent kinetic energy near the centerline. Finally, an inspection of the time-evolving instantaneous fields and conditional averaged flow fields for the streamwise growth of a very long structure near the centerline exhibits that a streamwise concatenation of adjacent large-scale u-structures creates a very-large-scale structure near the channel centerline.     

Direct numerical simulation of a particle-laden low Reynolds number turbulent round jet

Direct numerical simulation method is used for the investigating of particle-laden turbulent flows in a spatially evolution of low Reynolds number axisymmetric jet, and the Eulerian–Lagrangian point-particle approach is employed in the simulation. The simulation uses an explicit coupling scheme between particles and the fluid, which considers two-way coupling between the particle and the fluid. The DNS results are compared well with experimental data with equal Reynolds number (R_e = 1700). Our objects are: (i) to investigate the correlation between the particle number density and the fluctuating of fluid streamwise velocity; (ii) to examine whether the three-dimensional vortex structures in the particle-laden jet are the same as that in the free-air jet and how the particles modulate the thee-dimensional vortex structures and turbulence properties with different Stokes number particles; (iii) to discover the particle circumferential dispersion with different Stokes number particles. Our findings: (i) all the particles, regardless of their particle size, tend to preferentially accumulate in the region with large-than-mean fluid streamwise velocity; (ii) the small Stokes number particles take an important part in the modulation of three-dimensional vortex structures, but for the intermediate and larger sized particles, this modulation effect seems not so apparent; (iii) the particle circumferential dispersion is more effective for the smaller and intermediate sized particles, especially for the intermediate sized particles.     

Direct numerical simulation of particle interaction with ejections in turbulent channel flows

Direct numerical simulations (DNS) of incompressible turbulent channel flows coupled with Lagrangian particle tracking are performed to study the characteristics of ejections that surround solid particles. The behavior of particles in dilute turbulent channel flows, without particle collisions and without feedback of particles on the carrier fluid, is studied using high Reynolds number DNS (Re = 12,500). The results show that particles moving away from the wall are surrounded by ejections, confirming previous studies on this issue. A threshold value separating ejections with only upward moving particles is established. When normalized by the square root of the Stokes number and the square of the friction velocity, the threshold profiles follow the same qualitative trends, for all the parameters tested in this study, in the range of the experiments. When compared to suspension thresholds proposed by other studies in the Shields diagram, our simulations predict a much larger value because of the measure used to characterize the fluid and the criterion chosen to decide whether particles are influenced by the surrounding fluid. However, for intermediate particle Reynolds numbers, the threshold proposed here is in fair agreement with the theoretical criterion proposed by Bagnold (1966) [Bagnold, R., 1966. Geological Survey Professional Paper, vol. 422-1]. Nevertheless, further studies will be conducted to understand the normalization of the threshold.     

DNS of heat transfer in turbulent and transitional channel flow obstructed by rectangular prisms

Direct numerical simulation (DNS) of heat transfer in a channel flow obstructed by rectangular prisms has been performed for Re_τ = 80–20, where Re_τ is based on the friction velocity, the channel half width and the kinematic viscosity. The molecular Prandtl number is set to be 0.71. The flow remains unsteady down to Re_τ = 40 owing to the disturbance induced by the prism. For Re_τ = 30 and 20, the flow results in a steady laminar flow. In the vicinity of the prism, the three-dimensional complex vortices are generated and heat transfer is enhanced. The Reynolds number effect on the time-averaged vortex structure and the local Nusselt number are investigated. The mechanism of the heat transfer enhancement is discussed. In addition, the mean flow parameters such as the friction factor and the Nusselt number are examined in comparison with existing DNS and experimental data.     

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.     

Effect of coarse particles on the heat transfer in a particle-laden turbulent boundary layer

The temperature distribution in particle-laden turbulent flow, in a flume, was investigated both by DNS and experimentally. Simulations were performed at Re=171 and Pr=5.4 in order to study the interaction between the particle motion and flow turbulence. Two-way coupling was used to obtain various turbulence statistics, the grid resolution was sufficiently fine to resolve all essential turbulent scales. The effect of particle diameter on momentum, heat transfer and particle deposition was considered. The details of particle-turbulence interaction depend on the particle Stokes number and the particle Reynolds number.

The spatial structures of instantaneous flow and temperature fields were visualized. Low frequency small oscillations of deposited particles were observed. It was found that these small deviations from the initial position, caused strong changes in the instantaneous temperature field near the particle.

The experiments provided details of the temperature field on the heated wall close to the particle. In the front of the particle, a sharp increase in heat transfer coefficient was observed. The experimental results agree well with the computational predictions.     

Gas-solid interaction force from direct numerical simulation (DNS) of binary systems with extreme diameter ratios

nding the large particles becomes less in case ofa bidisperse mixture, as compared to a monodisperse system with the same volume fraction. We further investigated this discrepancy by calculating the volume per particle by means of Voronoi tessellation.     

Direct numerical simulation of turbulent wake behind a surface-mounted square cylinder

In this research, direct numerical simulation has been performed to study the turbulent wake behind a wall-mounted square cylinder with aspect ratio 4 and Reynolds number 12 000 (based on the free-stream velocity and obstacle side length) in a developing boundary layer. Owing to the relatively high Reynolds number and high aspect ratio of the cylinder tested, the wake is wide spread behind the cylinder and exhibits complex and energetic vortex motions. The lateral and tip vortex shedding patterns at different frequencies, coherent structures downstream of the obstacle, the production rate and distribution of turbulent kinetic energy, and the instantaneous pressure distribution in the wake region have been thoroughly investigated. In order to validate the numerical results, the first- and second-order flow statistics obtained from the simulations have been carefully compared against available wind-tunnel measurement data.     

Coherent structures in trailing-edge cooling and the challenge for turbulent heat transfer modelling

The present paper tests the capability of a standard Reynolds-Averaged Navier–Stokes (RANS) turbulence model for predicting the turbulent heat transfer in a generic trailing-edge situation with a cutback on the pressure side of the blade. The model investigated uses a gradient-diffusion assumption with a scalar turbulent-diffusivity and constant turbulent Prandtl number. High-fidelity Large-Eddy Simulations (LES) were performed for three blowing ratios to provide reliable target data and the mean velocity and eddy viscosity as input for the heat transfer model testing. Reasonably good agreement between the LES and recent experiments was achieved for mean flow and turbulence statistics. The LES yielded coherent structures which were analysed, in particular with respect to their effect on the turbulent heat transfer. For increasing blowing ratio, the LES replicated an also experimentally observed counter-intuitive decrease of the cooling effectiveness caused by the coherent structures becoming stronger. In contrast, the RANS turbulent heat transfer model failed in predicting this behaviour and yielded significantly too high cooling effectiveness. It is shown that the model cannot predict the strong upstream and wall-directed turbulent heat fluxes caused by large coherent structures, which were found to be responsible for the counter-intuitive decrease of the cooling effectiveness.     

Multiscale simulation of flow and heat transfer of nanofluid with lattice Boltzmann method

Based on the lattice Boltzmann (LB) approach, a novel hybrid method has been proposed for getting insight into the microscale characteristics of the multicomponent flow of nanofluid. In this method, the whole computational domain is divided into two regions in which different-sized meshes are involved for simulation (fine mesh and coarse mesh). The multicomponent LB method is adopted in the fine mesh region, and the single-component LB approach is applied to the coarse mesh region where the nanofluid is treated as a mixed single-component fluid. The conservation principles of mass, momentum and energy are used to derive a hybrid scheme across the different scaled regions. Numerical simulation is carried out for the Couette flow and convective heat transfer in a parallel plate channel to validate the hybrid method. The computational results indicate that by means of the present method, not only the microscopic characteristics of the nanofluid flow can be simulated, but also the computational efficiency can be remarkably improved compared with the pure multicomponent LB method.     

Direct numerical simulation of turbulent heat transfer over fully resolved anisotropic porous structures

Lattice Boltzmann direct numerical simulations of turbulent heat transfer over and inside anisotropic porous media are performed. This study considers turbulent plane channel flows whose bottom walls are made from the porous media at the bulk Reynolds number of 2900 with isothermal and conjugate heat transfer wall conditions. Four different porous walls are considered. They are walls with only the wall-normal permeability, with the wall-normal and spanwise permeabilities, with the wall-normal and streamwise permeabilities, and with the isotropic wall-normal, spanwise and streamwise permeabilities. The porosity of the porous walls ranges from 0.6 to 0.8. Discussions on the effects of the anisotropic permeability on turbulent thermal fields are carried out by the instantaneous flow visualizations and the statistical quantities. In particular, temperature fluctuations, turbulent and dispersion heat fluxes are examined both inside and outside the porous walls. Finally, the heat transfer performance is discussed considering the effects of the anisotropic permeability.     

Direct numerical simulation of turbulent pipe flow

Fully developed turbulent pipe flow at low Re-number is studied by means of direct numerical simulation (DNS). In contrast to many previous DNS's of turbulent flows in rectangular geometries, the present DNS code, developed for a cylindrical geometry, is based on the finite volume technique rather than being based on a spectral method. The statistical results are compared with experimental data obtained with two different experimental techniques. The agreement between numerical and experimental results is found to be good which indicates that the present DNS code is suited for this kind of numerical simulations.     

Effects of the roughness height in turbulent boundary layers over rod- and cuboid-roughened walls

Direct numerical simulations (DNSs) of spatially developing turbulent boundary layers (TBLs) over two-dimensional (2D) rod-roughened walls and three-dimensional (3D) cuboid-roughened walls are conducted to investigate the effects of the roughness height on the flow characteristics in the outer layer. The rod elements are periodically aligned along the downstream direction with a pitch of p_x/θ_in = 12, and the cuboid elements are periodically staggered with a pitch of p_x/θ_in = 12 and p_z/θ_in = 3, where p_x and p_z are correspondingly the streamwise and spanwise pitches of the roughness and θ_in is the momentum thickness at the inlet. The first surface roughness is placed 80θ_in downstream from the inlet, leading to a step change from a smooth to rough surface. The rod and cuboid roughness height (k) is varied in the range of 0.1 ≤ k/θ_in ≤ 1.8 (13 ≤ δ/k ≤ 285), respectively (δ is the boundary layer thickness), and the Reynolds number based on the momentum thickness (θ) is varied in the range of Re_θ = 300 ~ 1400. For each case, the self-preservation form of the velocity-defect and the turbulent Reynolds stresses is achieved along the downstream direction. As the roughness height increases, the roughness function (ΔU⁺) extracted from the mean velocity profiles increases, although the velocity-defect profiles for the rough-wall cases show good agreement with the profile from the smooth-wall case. The magnitude of the Reynolds stresses in the outer layer increases with an increase of k/δ. The outer layer similarity between the flows over the rough and smooth-walls is found when δ/k ≥ 250 and 100 for the 2D rod and 3D cuboid, respectively. The continuous increase of the Reynolds stresses in the outer layer with an increase of k/δ is explained by a large population of very long structures over the rough-wall flows. Because the characteristic width of the structures increases continuously with an increase of k/δ for the rod and cuboid roughness, a wide width of the structures leads to frequent spanwise merging between adjacent structures. The active spanwise merging events with an increase of k/δ increase the streamwise coherence of the structures with the appearance of significant meandering.