Direct Numerical Simulation of Turbulent Heat Transfer Across a Mobile, Sheared Gas-Liquid Interface

2003 ◽  
Vol 125 (6) ◽  
pp. 1129-1139 ◽  
Author(s):  
D. Lakehal ◽  
M. Fulgosi ◽  
G. Yadigaroglu ◽  
S. Banerjee
Keyword(s):  
Heat Transfer ◽  
Numerical Simulation ◽  
Direct Numerical Simulation ◽  
Gas Phase ◽  
Liquid Interface ◽  
Turbulent Heat Transfer ◽  
Turbulent Heat ◽  
Flow Data ◽  
Prandtl Numbers ◽  
The Impact

The impact of interfacial dynamics on turbulent heat transfer at a deformable, sheared gas-liquid interface is studied using Direct Numerical Simulation (DNS). The flow system comprises a gas and a liquid phase flowing in opposite directions. The governing equations for the two fluids are alternately solved in separate domains and then coupled at the interface by imposing continuity of velocity and stress. The deformations of the interface fall in the range of capillary waves of waveslope ak=0.01 (wave amplitude a times wavenumber k), and very small phase speed-to-friction velocity ratio, c/u*. The influence of low-to-moderate molecular Prandtl numbers Pr on the transport in the immediate vicinity of the interface is examined for the gas phase, and results are compared to existing wall-bounded flow data. The shear-based Reynolds number Re* is 171 and Prandtl numbers of 1, 5, and 10 were studied. The effects induced by changes in Pr in both wall-bounded flow and over a gas-liquid interface were analyzed by comparing the relevant statistical flow properties, including the budgets for the temperature variance and the turbulent heat fluxes. Overall, Pr was found to affect the results in very much the same way as in most of the available wall flow data. The intensity of the averaged normal heat flux at high Prandtl numbers is found to be slightly greater near the interface than at the wall. Similar to what is observed in wall flows, for Pr=1 the turbulent viscosity and diffusivity are found to asymptote with z+3, where z+ is the distance to the interface, and with z+n, where n>3 for Pr=5 and 10. This implies that the gas phase perceives deformable interfaces as impermeable walls for small amplitude waves with wavelengths much larger than the diffusive sublayers. Moreover, high-frequency fluctuating fields are shown to play a minor role in transferring heat across the interface, with a marked filtering effect of Pr. A new scaling law for the normalized heat transfer coefficient, K+ has been derived with the help of the DNS data. This law, which could be used in the range of Pr=1 to 10 for similar flow conditions, suggests an approximate Pr−3/5 relationship, lying between the Pr−1/2 dependence for free surfaces and the Pr−2/3 law for immobile interfaces and much higher Prandtl numbers. A close inspection of the transfer rates reveals a strong and consistent relationship between K+, the frequency of sweeps impacting the interface, the interfacial velocity streaks, and the interfacial shear stress.

Direct Numerical Simulation of the Flow in a Bare Rod Bundle at Different Prandtl Numbers

2019 ◽  
Vol 141 (12) ◽  
Author(s):  
Jonathan K. Lai ◽  
Giacomo Busco ◽  
Elia Merzari ◽  
Yassin A. Hassan
Keyword(s):  
Heat Transfer ◽  
Numerical Simulation ◽  
Reynolds Number ◽  
Prandtl Number ◽  
Direct Numerical Simulation ◽  
Frequency Analysis ◽  
Turbulent Heat Transfer ◽  
Turbulent Heat ◽  
Rod Bundle ◽  
Prandtl Numbers

Abstract A direct numerical simulation (DNS) of bare rod bundles with a low pitch-to-diameter ratio is performed with heat transfer at different Prandtl numbers. Turbulence statistics for temperature and velocity as well as the turbulent budgets have been collected. High-fidelity simulations are performed with the spectral element method (SEM) using Nek5000, a highly scalable code. To pertain to industrial-related flows, a rod bundle model is based on Hooper and Wood's (Hooper, J. D., and Wood, D., 1984, “Fully Developed Rod Bundle Flow Over a Large Range of Reynolds Number,” Nucl. Eng. Des., 83(1), pp. 31–46) experimental setup. Both wall normalized velocity profile and turbulent kinetic energy are validated with a Reynolds number of 22,600. Kolmogorov length scales and time scales are calculated to be within the simulation's spatial–temporal resolution. Moreover, gap vortices and coherent structures are quantified by using Lambda2 vortex criterion, frequency analysis, and two-point correlation. Heat transfer statistics are discussed with a constant heat flux for six different Prandtl numbers ranging from 2 to 0.002. This range shows significantly different characteristics in temperature for both mean and variance. Mean temperature profiles in the subchannel center are very sensitive to the Prandtl number when it becomes small. It is also found that the location of the local maxima for the variance of temperature fluctuations becomes very sensitive at larger Prandtl numbers. The temperature frequency analysis reveals a shift to lower frequencies for low Prandtl numbers. The DNS results provided in this work will contribute as benchmark for the improvement and development of existing and new turbulent heat transfer models at different Prandtl number regimes.

Direct numerical simulation of turbulent heat transfer across a sheared wind-driven gas–liquid interface

2016 ◽  
Vol 804 ◽  
pp. 646-687 ◽  
Author(s):  
Ryoichi Kurose ◽  
Naohisa Takagaki ◽  
Atsushi Kimura ◽  
Satoru Komori
Keyword(s):  
Heat Transfer ◽  
Numerical Simulation ◽  
Mass Transfer ◽  
Heat Flux ◽  
Streamwise Velocity ◽  
Liquid Interface ◽  
Turbulent Heat Transfer ◽  
Turbulent Heat ◽  
Transfer Coefficients ◽  
Turbulent Eddies

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.

Sign in / Sign up

Export Citation Format

Share Document