Heat Transfer Modulation by Inertial Particles in Particle-Laden Turbulent Channel Flow

2018 ◽  
Vol 140 (11) ◽  
Author(s):  
Caixi Liu ◽  
Shuai Tang ◽  
Yuhong Dong ◽  
Lian Shen

We study the effects of particle-turbulence interactions on heat transfer in a particle-laden turbulent channel flow using an Eulerian–Lagrangian simulation approach, with direct numerical simulation (DNS) for turbulence and Lagrangian tracking for particles. A two-way coupling model is employed in which the momentum and energy exchange between the discrete particles and the continuous fluid phase is fully taken into account. Our study focuses on the modulations of the temperature field and heat transfer process by inertial particles with different particle momentum Stokes numbers (St), which in a combination of the particle-to-fluid specific heat ratio and the Prandtl number results in different particle heat Stokes numbers. It is found that as St increases, while the turbulent heat flux decreases due to the suppression of wall-normal turbulence velocity fluctuation, the particle feedback heat flux increases significantly and results in an increase in the total heat flux. The particle thermal feedback effect is illustrated using the instantaneous structures and statistics of the flow and temperature fields. The mechanisms of heat transfer modulation by inertial particles are investigated in detail. The budget of turbulent heat flux is examined. Moreover, by taking advantage of the ability of numerical simulation to address different momentum and heat processes separately, we investigate in detail the two processes of particles affecting heat transfer for the first time, namely the direct effect of particle thermal feedback to the fluid (i.e., heat feedback) and the indirect effect of the modulation of turbulent velocity field induced by the particles (i.e., momentum feedback). It is found that the contribution of heat transfer from turbulent convection is reduced by both heat and momentum feedback due to the decrease of the turbulent heat flux. The contribution of heat transfer from particle transport effects is barely influenced by the momentum feedback, even if St is large and is mainly affected by the heat feedback. Our results indicate that both heat and momentum feedback are important when the particle inertia is large, suggesting that both feedback processes need to be taken into account in computation and modeling.

2015 ◽  
Vol 776 ◽  
pp. 512-530 ◽  
Author(s):  
S. Leonardi ◽  
P. Orlandi ◽  
L. Djenidi ◽  
R. A. Antonia

Direct numerical simulations (DNS) are carried out to study the passive heat transport in a turbulent channel flow with either square bars or circular rods on one wall. Several values of the pitch (${\it\lambda}$) to height ($k$) ratio and two Reynolds numbers are considered. The roughness increases the heat transfer by inducing ejections at the leading edge of the roughness elements. The amounts of heat transfer and mixing depend on the separation between the roughness elements, an increase in heat transfer accompanying an increase in drag. The ratio of non-dimensional heat flux to the non-dimensional wall shear stress is higher for circular rods than square bars irrespectively of the pitch to height ratio. The turbulent heat flux varies within the cavities and is larger near the roughness elements. Both momentum and thermal eddy diffusivities increase relative to the smooth wall. For square cavities (${\it\lambda}/k=2$) the turbulent Prandtl number is smaller than for a smooth channel near the wall. As ${\it\lambda}/k$ increases, the turbulent Prandtl number increases up to a maximum of 2.5 at the crests plane of the square bars (${\it\lambda}/k=7.5$). With increasing distance from the wall, the differences with respect to the smooth wall vanish and at three roughness heights above the crests plane, the turbulent Prandtl number is essentially the same for smooth and rough walls.


2016 ◽  
Vol 138 (11) ◽  
Author(s):  
Alan S. Hsieh ◽  
Sedat Biringen ◽  
Alec Kucala

A direct numerical simulation (DNS) of spanwise-rotating turbulent channel flow was conducted for four rotation numbers: Rob=0, 0.2, 0.5, and 0.9 at a Reynolds number of 8000 based on laminar centerline mean velocity and Prandtl number 0.71. The results obtained from these DNS simulations were utilized to evaluate several turbulence closure models for momentum and heat transfer transport in rotating turbulent channel flow. Four nonlinear eddy viscosity turbulence models were tested and among these, explicit algebraic Reynolds stress models (EARSM) obtained the Reynolds stress distributions in best agreement with DNS data for rotational flows. The modeled pressure–strain functions of EARSM were shown to have strong influence on the Reynolds stress distributions near the wall. Turbulent heat flux distributions obtained from two explicit algebraic heat flux models (EAHFM) consistently displayed increasing disagreement with DNS data with increasing rotation rate.


Author(s):  
Kyoungyoun Kim ◽  
Radhakrishna Sureshkumar

A direct numerical simulation (DNS) of viscoelastic turbulent channel flow with the FENE-P model was carried out to investigate turbulent heat transfer mechanism of polymer drag-reduced flows. The configuration was a fully-developed turbulent channel flow with uniform heat flux imposed on both walls. The temperature was considered as a passive scalar. The Reynolds number based on the friction velocity (uτ) and channel half height (δ) is 125 and Prandtl number is 5. Consistently with the previous experimental observations, the present DNS results show that the heat-transfer coefficient was reduced at a rate faster than the accompanying drag reduction rate. Statistical quantities such as root-mean-square temperature fluctuations and turbulent heat fluxes were obtained and compared with those of a Newtonian fluid flow. Budget terms of the turbulent heat fluxes were also presented.


2014 ◽  
Vol 753 ◽  
pp. 360-401 ◽  
Author(s):  
R. Vicquelin ◽  
Y. F. Zhang ◽  
O. Gicquel ◽  
J. Taine

AbstractThe role of radiative energy transfer in turbulent boundary layers is carefully analysed, focusing on the effect on temperature fluctuations and turbulent heat flux. The study is based on direct numerical simulations (DNS) of channel flows with hot and cold walls coupled to a Monte-Carlo method to compute the field of radiative power. In the conditions studied, the structure of the boundary layers is strongly modified by radiation. Temperature fluctuations and turbulent heat flux are reduced, and new radiative terms appear in their respective balance equations. It is shown that they counteract turbulence production terms. These effects are analysed under different conditions of Reynolds number and wall temperature. It is shown that collapsing of wall-scaled profiles is not efficient when radiation is considered. This drawback is corrected by the introduction of a radiation-based scaling. Finally, the significant impact of radiation on turbulent heat transfer is studied in terms of the turbulent Prandtl number. A model for this quantity, based on the new proposed scaling, is developed and validated.


2016 ◽  
Vol 804 ◽  
pp. 646-687 ◽  
Author(s):  
Ryoichi Kurose ◽  
Naohisa Takagaki ◽  
Atsushi Kimura ◽  
Satoru Komori

Turbulent heat transfer across a sheared wind-driven gas–liquid interface is investigated by means of a direct numerical simulation of gas–liquid two-phase turbulent flows under non-breaking wave conditions. The wind-driven wavy gas–liquid interface is captured using the arbitrary Lagrangian–Eulerian method with boundary-fitted coordinates on moving grids, and the temperature fields on both the gas and liquid sides, and the humidity field on the gas side are solved. The results show that although the distributions of the total, latent, sensible and radiative heat fluxes at the gas–liquid interface exhibit streak features such that low-heat-flux regions correspond to both low-streamwise-velocity regions on the gas side and high-streamwise-velocity regions on the liquid side, the similarity between the heat-flux streak and velocity streak on the gas side is more significant than that on the liquid side. This means that, under the condition of a fully developed wind-driven turbulent field on both the gas and liquid sides, the heat transfer across the sheared wind-driven gas–liquid interface is strongly affected by the turbulent eddies on the gas side, rather than by the turbulent eddies and Langmuir circulations on the liquid side. This trend is quite different from that of the mass transfer (i.e. $\text{CO}_{2}$ gas). This is because the resistance to heat transfer is normally lower than the resistance to mass transfer on the liquid side, and therefore the heat transfer is controlled by the turbulent eddies on the gas side. It is also verified that the predicted total heat, latent heat, sensible heat and enthalpy transfer coefficients agree well with previously measured values in both laboratory and field experiments. To estimate the heat transfer coefficients on both the gas and liquid sides, the surface divergence could be a useful parameter, even when Langmuir circulations exist.


2016 ◽  
Vol 798 ◽  
pp. 299-349 ◽  
Author(s):  
Bo Sun ◽  
Sudheer Tenneti ◽  
Shankar Subramaniam ◽  
Donald L. Koch

Fluctuations in the gas-phase velocity can contribute significantly to the total gas-phase kinetic energy even in laminar gas–solid flows as shown by Mehrabadi et al. (J. Fluid Mech., vol. 770, 2015, pp. 210–246), and these pseudo-turbulent fluctuations can also enhance heat transfer in gas–solid flow. In this work, the pseudo-turbulent heat flux arising from temperature–velocity covariance, and average fluid-phase conduction during convective heat transfer in a gas–solid flow are quantified and modelled over a wide range of mean slip Reynolds number and solid volume fraction using particle-resolved direct numerical simulations (PR-DNS) of steady flow through a random assembly of fixed isothermal monodisperse spherical particles. A thermal self-similarity condition on the local excess temperature developed by Tenneti et al. (Intl J. Heat Mass Transfer, vol. 58, 2013, pp. 471–479) is used to guarantee thermally fully developed flow. The average gas–solid heat transfer rate for this flow has been reported elsewhere by Sun et al. (Intl J. Heat Mass Transfer, vol. 86, 2015, pp. 898–913). Although the mean velocity field is homogeneous, the mean temperature field in this thermally fully developed flow is inhomogeneous in the streamwise coordinate. An exponential decay model for the average bulk fluid temperature is proposed. The pseudo-turbulent heat flux that is usually neglected in two-fluid models of the average fluid temperature equation is computed using PR-DNS data. It is found that the transport term in the average fluid temperature equation corresponding to the pseudo-turbulent heat flux is significant when compared to the average gas–solid heat transfer over a significant range of solid volume fraction and mean slip Reynolds number that was simulated. For this flow set-up a gradient-diffusion model for the pseudo-turbulent heat flux is found to perform well. The Péclet number dependence of the effective thermal diffusivity implied by this model is explained using a scaling analysis. Axial conduction in the fluid phase, which is often neglected in existing one-dimensional models, is also quantified. As expected, it is found to be important only for low Péclet number flows. Using the exponential decay model for the average bulk fluid temperature, a model for average axial conduction is developed that verifies standard assumptions in the literature. These models can be used in two-fluid simulations of heat transfer in fixed beds. A budget analysis of the mean fluid temperature equation provides insight into the variation of the relative magnitude of the various terms over the parameter space.


Author(s):  
Vera Papp ◽  
Andrea Pucciarelli ◽  
Medhat Sharabi ◽  
Walter Ambrosini

This work proposes simulations of heat transfer under supercritical pressure conditions showing improvements with respect to previous works. This is obtained by the introduction of the Algebraic Heat Flux Model (AHFM) for evaluating the turbulent heat flux in turbulence production terms, using the in-house code THEMAT and the STAR-CCM+ code. The first code makes use of the AHFM also in the energy balance equations, while for the commercial code simplifying assumptions are considered in the implementations. Custom sets of parameters for every condition of inlet temperature and internal diameter are tuned in some cases, driven by the opinion that a single set of parameters cannot be suitable in every flow conditions, considering the complexity of the variables that concur in the heat transfer deterioration phenomenon. The AHFM model gives promising results with new sets of parameters in order to model the deterioration and the recovery phases because of its term related to the variance of temperature.


Author(s):  
I. Otic´ ◽  
G. Gro¨tzbach

Results of direct numerical simulation (DNS) of turbulent Rayleigh-Be´nard convection for a Prandtl number Pr = 0.025 and a Rayleigh number Ra = 105 are used to evaluate the turbulent heat flux and the temperature variance. The DNS evaluated turbulent heat flux is compared with the DNS based results of a standard gradient diffusion turbulent heat flux model and with the DNS based results of a standard algebraic turbulent heat flux model. The influence of the turbulence time scales on the predictions by the standard algebraic heat flux model at these Rayleigh- and Prandtl numbers is investigated. A four equation algebraic turbulent heat flux model based on the transport equations for the turbulent kinetic energy k, for the dissipation of the turbulent kinetic energy ε, for the temperature variance θ2, and for the temperature variance dissipation rate εθ is proposed. This model should be applicable to a wide range of low Prandtl number flows.


Sign in / Sign up

Export Citation Format

Share Document