On the role of the smallest scales of a passive scalar field in a near-wall turbulent flow

Role of the smallest diffusive scales of a passive scalar field in the near-wall turbulent flow was examined with pseudo-spectral numerical simulations. Temperature fields were analyzed at friction Reynolds number Re _τ=171 and at Prandtl numbers, Pr=1 and Pr=5.4. Results of direct numerical simulations (DNS) were compared with the under-resolved simulations where the velocity field was still resolved with the DNS accuracy, while a coarser grid was used to describe the temperature fields. Since the smallest temperature scales remained unresolved in these simulations, an appropriate spectral turbulent thermal diffusivity was applied to avoid pile-up at the higher wave numbers. In spite of coarser numerical grids, the temperature fields are still highly correlated with the DNS results, including instantaneous temperature fields. Results point to practically negligible role of the diffusive temperature scales on the macroscopic behavior of the turbulent heat transfer.     

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.     

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.     

Conjugate heat transfer study of incompressible turbulent offset jet flows

In the present case, the conjugate heat transfer involving a turbulent plane offset jet is considered. The bottom wall of the solid block is maintained at an isothermal temperature higher than the jet inlet temperature. The parameters considered are the offset ratio (OR), the conductivity ratio (K), the solid slab thickness (S) and the Prandtl number (Pr). The Reynolds number considered is 15,000 because the flow becomes fully turbulent and then it becomes independent of the Reynolds number. The ranges of parameters considered are: OR = 3, 7 and 11, K = 1–1,000, S = 1–10 and Pr = 0.01–100. High Reynolds number two-equation model (k–ε) has been used for turbulence modeling. Results for the solid–fluid interface temperature, local Nusselt number, local heat flux, average Nusselt number and average heat transfer have been presented and discussed.     

High-resolution imaging of turbulence structures in jet flames and non-reacting jets with laser Rayleigh scattering

A comparative study of the length scales and morphology of dissipation fields in turbulent jet flames and non-reacting jets provides a quantitative analysis of the effects of heat release on the fine-scale structure of turbulent mixing. Planar laser Rayleigh scattering is used for highly resolved measurements of the thermal and scalar dissipation in the near fields of CH₄/H₂/N₂ jet flames (Re _d  = 15,200 and 22,800) and non-reacting propane jets (Re _d  = 7,200–21,700), respectively. Heat release increases the dissipation cutoff length scales in the reaction zone of the flames such that they are significantly larger than the cutoff scales of non-reacting jets with comparable jet exit Reynolds numbers. Fine-scale anisotropy is enhanced in the reaction zone. At x/d = 10, the peaks of the dissipation angle PDFs in the Re _d  = 15,200 and 22,800 jet flames exceed those of non-reacting jets with corresponding jet exit Reynolds numbers by factors of 2.3 and 1.8, respectively. Heat release significantly reduces the dissipation layer curvature in the reaction zone and in the low-temperature periphery of the jet flames. These results suggest that the reaction zone shields the outer regions of the jet flame from the highly turbulent flow closer to the jet axis.     

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.     

Joint Scalar versus Joint Velocity-Scalar PDF Simulations of Bluff-Body Stabilized Flames with REDIM

Two transported PDF strategies, joint velocity-scalar PDF (JVSPDF) and joint scalar PDF (JSPDF), are investigated for bluff-body stabilized jet-type turbulent diffusion flames with a variable degree of turbulence–chemistry interaction. Chemistry is modeled by means of the novel reaction-diffusion manifold (REDIM) technique. A detailed chemistry mechanism is reduced, including diffusion effects, with N ₂ and CO ₂ mass fractions as reduced coordinates. The second-moment closure RANS turbulence model and the modified Curl’s micro-mixing model are not varied. Radiative heat loss effects are ignored. The results for mean velocity and velocity fluctuations in physical space are very similar for both PDF methods. They agree well with experimental data up to the neck zone. Each of the two PDF approaches implies a different closure for the velocity-scalar correlation. This leads to differences in the radial profiles in physical space of mean scalars and mixture fraction variance, due to different scalar flux modeling. Differences are visible in mean mixture fraction and mean temperature, as well as in mixture fraction variance. In principle, the JVSPDF simulations can be closer to physical reality, as a differential model is implied for the scalar fluxes, whereas the gradient diffusion hypothesis is implied in JSPDF simulations. Yet, in JSPDF simulations, turbulent diffusion can be tuned by means of the turbulent Schmidt number. In the neck zone, where the turbulent flow field results deteriorate, the joint scalar PDF results are in somewhat better agreement with experimental data, for the test cases considered. In composition space, where results are reported as scatter plots, differences between the two PDF strategies are small in the calculations at hand, with a little more local extinction in the joint scalar PDF results.     

Thermo-mechanical modeling of turbulent heat transfer in gas–solid flows including particle collisions

A thermo-mechanical turbulence model is developed and used for predicting heat transfer in a gas–solid flow through a vertical pipe with constant wall heat flux. The new four-way interaction model makes use of the thermal k_θ–τ_θ equations, in addition to the hydrodynamic k–τ transport, and accounts for the particle–particle and particle–wall collisions through a Eulerian/Lagrangian formulation. The simulation results indicate that the level of thermal turbulence intensity and the heat transfer are strongly affected by the particle collisions. Inter-particle collisions attenuate the thermal turbulence intensity near the wall but somewhat amplify the temperature fluctuations in the pipe core region. The hydrodynamic-to-thermal times-scale ratio and the turbulent Prandtl number in the region near the wall increase due to the inter-particle collisions. The results also show that the use of a constant or the single-phase gas turbulent Prandtl number produces error in the thermal eddy diffusivity and thermal turbulent intensity fields. Simulation results also indicate that the inter-particle contact heat conduction during collision has no significant effect in the range of Reynolds number and particle diameter studied.     

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.     

Influence of buoyancy in a mixed convection liquid metal flow for a horizontal channel configuration

This article presents the direct numerical simulation (DNS) of mixed convection turbulent heat transfer in a horizontal channel case for liquid lead. Cartesian mesh is used and the incompressible Navier-Stokes equations are discretized with highly accurate finite difference sixth-order compact schemes to perform the DNS. The influence of mixed convection in liquid metal with Prandtl number equal to 0.025 and Reynolds number equal to 4667 has been studied by varying the Richardson number (Ri = 0, 0.25, 0.50, 1.00). The obtained results are extensively analyzed and discussed in this article. In particular, large-scale circulation is observed under the influence of buoyancy. Compared to the forced convection case (Ri = 0), stronger velocity fluctuations are noticed that highlight the fact that turbulence is strongly enhanced with the increasing buoyancy. It also proves that the thermal plumes rising up from the hot wall of the channel activate the cross-stream eddies. Moreover, temperature fluctuations are found to be more homogeneously distributed with increasing buoyancy effects and mixing is more effective in the center of the channel. In addition, compared with forced convection, mixed convection has shown enlargement of the large-scale structures that only appear in the temperature field for low Prandtl number fluids. Extensive results of flow and temperature fields are analyzed and presented.     

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.     

Influence of the variable thermophysical properties on the turbulent buoyancy-driven airflow inside open square cavities

The effects of the air variable properties (density, viscosity and thermal conductivity) on the buoyancy-driven flows established in open square cavities are investigated, as well as the influence of the stated boundary conditions at open edges and the employed differencing scheme. Two-dimensional, laminar, transitional and turbulent simulations are obtained, considering both uniform wall temperature and uniform heat flux heating conditions. In transitional and turbulent cases, the low-Reynolds k − ω turbulence model is employed. The average Nusselt number and the dimensionless mass-flow rate have been obtained for a wide and not yet covered range of the Rayleigh number varying from 10³ to 10¹⁶. The results obtained taking into account variable properties effects are compared with those calculated assuming constant properties and the Boussinesq approximation. For uniform heat flux heating, a correlation for the critical heating parameter above which the burnout phenomenon can be obtained is presented, not reported in previous works. The effects of variable properties on the flow patterns are analyzed.     

Spatial resolution effects on the measurement of scalar variance and scalar gradient in turbulent nonpremixed jet flames

The effect of spatial averaging is important for scalar gradient measurements in turbulent nonpremixed flames, especially when the local dissipation length scale is small. Line imaging of Raman, Rayleigh and CO-LIF is used to investigate the effects of experimental resolution on the scalar variance and radial gradient in the near field of turbulent nonpremixed CH₄/H₂/N₂ jet flames at Reynolds numbers of 15,200 and 22,800 (DLR-A and B) and in piloted CH₄/air jet flames at Reynolds numbers of 13,400, 22,400 and 33,600 (Sandia flames C/D/E). The finite spatial resolution effects are studied by applying the Box filter with varying filter widths. The resulting resolution curves for both scalar variance and mean squared-gradient follow nearly the same trends as theoretical curves calculated from the model turbulence kinetic energy spectrum of Pope. The observed collapse of resolution curves of mean squared-gradient for nearly all studied cases implies the shape of the dissipation spectrum is approximately universal. Fluid transport properties are shown to have no effect on the dissipation resolution curve, which implies that the dissipation length scale inferred from the square gradient is equivalent to the length scale for the scalar dissipation rate, which includes the diffusion coefficient. With the Box filter, the required spatial resolution to resolve 98% of the mean dissipation rate is about one−two times of the dissipation cutoff length scale (analogous to the Batchelor scale in turbulent isothermal flows). The effects of resolution on the variances of mixture fraction, temperature, and the inverted Rayleigh signal are also compared. The ratio of the filtered variance to the true variance is shown to depend nearly linearly on the probe resolution. The inverted Rayleigh scattering signal can be used to study the resolution effect on temperature variance even when the Rayleigh scattering cross section is not constant. The experimental results also indicate that these laboratory scale turbulent jet flames have small effective Reynolds numbers, such that there is some direct interaction of the large (energy containing) and small (dissipative) scalar length scales, especially for the near field case at x/d = 7.5 of the piloted Sandia flames C/D/E.     

Heat and mass transfer study of impinging turbulent premixed flames

 Impinging jet combusting flows on granite plates are studied. A mathematical model for calculating heat release in turbulent impinging premixed flames is developed. The combustion including radiative heat transfer and local extinction effects, and flow characteristics are modeled using a finite volume computational approach. Two different eddy viscosity turbulence models, namely the standard k–ɛ and the RNG k–ɛ model with and without radiation (discrete transfer model) are assessed. The heat released predictions are compared with experimental data and the agreement is satisfactory only when both radiative heat transfer and local extinction modeling are taken into account. The results indicate that the main effect of radiation is the decrease of temperature values near the jet stagnation point and along the plate surface. Radiation increases temperature gradients and affects predicted turbulence levels independently of the closure model used. Also, the RNG k–ɛ predicts higher temperatures close the solid plate, with and without radiative heat transfer. Received on 13 November 2000 / Published online: 29 November 2001     

Spatial and Temporal Characteristics of OH in Turbulent Opposed-Jet Double Flames

Simultaneous high repetition-rate, two-point hydroxyl (OH) time-series measurements with associated PLIF/PIV measurements are employed to investigate spatio–temporal scales and flame-velocity interactions in turbulent opposed jets sustaining methane-air double flames. For a fuel-side equivalence ratio, ϕ _B  = 1.2, a rich premixed flame exists on the fuel side while a diffusion flame exists on the air side of the stagnation plane. The bulk Reynolds number (Re) and strain rate (SR) can be adjusted to generate flames at ϕ _B  = 1.2 with both well separated and completely merged flame fronts. Simultaneous PLIF/PIV measurements highlight distinct spatial OH structures of the premixed and diffusive fronts corresponding to variations in the flow field. The self-propagating tendency of the rich premixed front causes large-scale wrinkling, thereby enhancing the OH contour length by 15% as compared to the diffusive front. Two-point OH time-series measurements are implemented to quantify both spatial and temporal fluctuations via study of radial length and time scales. In general, these integral length and time scales follow similar trends and reach a minimum at the axial location of peak [OH]. In comparison to merged double flames having higher Re and SR, greater OH fluctuations are observed in the rich-premixed front as compared to the diffusive front for a well separated double flame. Because of the developing turbulence, the OH length scales exhibit reduced axial gradients across the reaction zone for higher Re in comparison to lower Re. A stochastic time-series simulation, using a state relationship based on a joint mixture fraction and progress variable, is utilized to extract estimated scalar time scales from those of measured OH. The simulations indicate that the hydroxyl fluctuations in double flames are only twice those of the underlying conserved scalar. “Turbulent Opposed-Jet Double Flames” is submitted for consideration as a full length article to Flow Turbulence and Combustion.     

Large eddy simulations of a turbulent thermal plume

Large eddy simulations of a three-dimensional turbulent thermal plume in an open environment have been carried out using a self-developed parallel computational fluid dynamics code SMAFS (smoke movement and flame spread) to study the thermal plume’s dynamics including its puffing, self-preserving and air entrainment. In the simulation, the sub-grid stress was modeled using both the standard Smagorinsky and the buoyancy modified Smagorinsky models, which were compared. The sub-grid scale (SGS) scalar flux in the filtered enthalpy transport equation was modeled based on a simple gradient transport hypothesis with constant SGS Prandtl number. The effect of the Smagorinsky model constant and the SGS Prandtl number were examined. The computation results were compared with experimental measurements, thermal plume theory and empirical correlations, showing good agreement. It is found that both the buoyancy modification and the SGS turbulent Prandtl number have little influence on simulation. However, the SGS model constant C _s has a significant effect on the prediction of plume spreading, although it does not affect much the prediction of puffing.     

A method to simultaneously image two-dimensional mixture fraction, scalar dissipation rate, temperature and fuel consumption rate fields in a turbulent non-premixed jet flame

A new imaging technique was developed that provides two-dimensional images of the mixture fraction (ξ), scalar dissipation rate (χ), temperature (T), and fuel consumption rate in a turbulent non-premixed jet flame. The new method is based on “seeding” nitric oxide (NO) into a particular carbon monoxide–air flame in which it remains passive. It is first demonstrated that the mass fraction of NO is a conserved scalar in the present carbon monoxide–air flame configuration, using both laminar flame calibration experiments and computations with full chemistry. Simultaneous planar laser-induced fluorescence (PLIF) and planar Rayleigh scattering temperature imaging allow a quantitative determination of the local NO mass fraction and hence mixture fraction in the turbulent jet flame. The instantaneous mixture fraction fields in conjunction with the local temperature fields are then used to determine quantitative scalar dissipation rate fields. Advantages of the present technique include an improved signal-to-noise ratio over previous Raman scattering techniques, improved accuracy near the stoichiometric contour because simplifying chemistry assumptions are not required, and the ability to measure ξ and χ in flames experiencing localized extinction. However, the method of measuring ξ based on the passive NO is restricted to dry carbon monoxide–air flames due to the well-controlled flame chemistry. Sample imaging results for ξ, χ, T, and are presented that show high levels of signal-to-noise while resolving the smallest mixing scales of the turbulent flowfield. The application, accuracy, and limitations of the present technique are discussed.     

Turbulent flow simulation of liquid jet emanating from pressure-swirl atomizer

The performances of three linear eddy viscosity models (LEVM) and one algebraic Reynolds stress model (ARSM) for the simulation of turbulent flow inside and outside pressure-swirl atomizer are evaluated by comparing the interface position with available experimental data and by comparing the turbulence intensity profiles at the atomizer exit. It is found that the turbulence models investigated exhibit zonal behaviors, i.e. none of the models investigated performs well throughout the entire flow field. The turbulence intensity has a significant influence on the global characteristics of the flow field. The turbulence models with better predictions of the turbulence intensity, such as Gatski-Speziale’s ARSM model, can yield better predictions of the global characteristics of the flow field, e.g. the reattachment lengths for the backward-facing step flow and the sudden expansion pipe flow, or the discharge coefficient, film thickness and the liquid sheet outer surface position for the atomizer flows. The standard k–ε model predicts stronger turbulence intensity as compared to the other models and therefore yields smaller film thickness and larger liquid sheet outer surface position. In average, the ARSM model gives both quantitatively and qualitatively better results as compared to the standard k–ε model and the low Reynolds number models.     

Flow and Mixing Fields for Transported Scalar PDF Simulations of a Piloted Jet Diffusion Flame (‘Delft Flame III’)

Numerical simulation results are presented for ‘Delft Flame III’, a piloted jet diffusion flame with strong turbulence–chemistry interaction. While pilot flames emerge from 12 separate holes in the experiments, the simulations are performed on a rectangular grid, under the assumption of axisymmetry. In the first part of the paper, flow and mixing field results are presented with a non-linear first order k–ε model, with the transport equation for ε based on a modeled enstrophy transport equation, for cold and reactive flows. For the latter, the turbulence model is applied in combination with pre-assumed β-PDF modeling for the turbulence–chemistry interaction. The mixture fraction serves as conserved scalar. Two chemistry models are considered: chemical equilibrium and a steady laminar flamelet model. The importance of the turbulence model is highlighted. The influence of the chemistry model is noticeable too. A procedure is described to construct appropriate inlet boundary conditions. Still, the generation of accurate inlet boundary conditions is shown to be far less important, their effect being local, close to the nozzle exit. In the second part of the paper, results are presented with the transported scalar PDF approach as turbulence–chemistry interaction model. A C₁ skeletal scheme serves as chemistry model, while the EMST method is applied as micro-mixing model. For the transported PDF simulations, the model for the pilot flames, as an energy source term in the mean enthalpy transport equation, is important with respect to the accuracy of the flow field predictions. It is explained that the strong influence on the flow and mixing field is through the turbulent shear stress force in the region, close to the nozzle exit.