## Abstract

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.

1.
Komori
,
S.
,
Nagaosa
,
R.
, and
Murakami
,
Y.
,
1993
, “
Turbulence Structure and Mass Transfer Across a Sheared Air-Water Interface in Wind-Driven Turbulence
,”
J. Fluid Mech.
,
249
, p.
161
161
.
2.
Rashidi
,
M.
,
Hetsroni
,
G.
, and
Banerjee
,
S.
,
1992
, “
Wave-Turbulence Interaction in Free-Surface Channel Flows
,”
Phys. Fluids A
,
4
, p.
2727
2727
.
3.
Danckwerts
,
P. V.
,
1951
, “
Significance of Liquid Film Coefficients in Gas Absorption
,”
Ind. Eng. Chem.
,
43
, p.
1460
1460
.
4.
Kasagi
,
N.
,
Tomita
,
Y.
, and
Kuroda
,
A.
,
1992
, “
Direct Numerical Simulation of Passive Scalar Field in a Turbulent Channel Flow
,”
ASME J. Heat Transfer
,
114
, p.
598
598
.
5.
Kawamura
,
H.
,
Ohsake
,
K.
,
Abe
,
H.
, and
Yamamoto
,
K.
,
1998
, “
DNS of Turbulent Heat Transfer in Channel Flow With Low to Medium-High Prandtl Number
,”
Int. J. Heat Fluid Flow
,
19
, p.
482
482
.
6.
Lyons
,
S. L.
,
Hanratty
,
T. J.
, and
McLaughlin
,
J. B.
,
1991
, “
Large-Scale Computer Simulation of Fully Developed Turbulent Channel Flow With Heat Transfer
,”
Int. J. Numer. Methods Fluids
,
13
, p.
999
999
.
7.
Na
,
Y.
,
Papavassiliou
,
D. V.
, and
Hanratty
,
T. J.
,
1999
, “
Use of Direct Numerical Simulation to Study the Effect of Prandtl Number on Temperature Fields
,”
Int. J. Heat Fluid Flow
,
20
, p.
187
187
.
8.
Campbell
,
J. A.
, and
Hanratty
,
T. J.
,
1983
, “
Mechanisms of Turbulent Mass Transfer at Solid Boundaries
,”
AIChE J.
,
29
, p.
221
221
.
9.
,
M. J.
,
,
E.
, and
Hanratty
,
T. J.
,
1986
, “
Computer Simulation of Turbulent Mass Transfer at a Mobile Interface
,”
AIChE J.
,
32
, p.
1108
1108
.
10.
Calmet
,
I.
, and
Magnaudet
,
J.
,
1997
, “
Large-Eddy Simulation of High-Schmidt Number Mass Transfer in a Turbulent Channel Flow
,”
Phys. Fluids
,
9
, p.
438
438
.
11.
Piller
,
M.
,
Nobile
,
E.
, and
Hanratty
,
T. J.
,
2002
, “
DNS Study of Turbulent Transport at Low Prandtl Numbers in a Channel Flow
,”
J. Fluid Mech.
,
458
, p.
419
419
.
12.
Lombardi
,
P.
,
De Angelis
,
V.
, and
Banerjee
,
S.
,
1996
, “
Direct Numerical Simulation of Near-Interface Turbulence in Coupled Gas-Liquid Flow
,”
Phys. Fluids
,
8
, p.
1643
1643
.
13.
De Angelis
,
V.
,
Lombardi
,
P.
, and
Banerjee
,
S.
,
1997
, “
Direct Numerical Simulation of Turbulent Flow Over a Wavy Wall
,”
Phys. Fluids
,
9
, p.
2429
2429
.
14.
De Angelis, V., Lombardi, P., Andreussi, P., and Banerjee, S., 1997, “MicroPhysiscs of Scalar Transfer at Air-Water Interfaces,” in Proceedings of Wind-Over-Wave Couplings: Perspectives and Prospects, S. G. Sajjadi, N. H. Thomas, and J. C. R. Hunt, eds., Oxford University Press.
15.
Fulgosi
,
M.
,
Lakehal
,
D.
,
Banerjee
,
S.
, and
De Angelis
,
V.
,
2003
, “
Direct Numerical Simulation of Turbulence in a Sheared Air-Water Flow With Deformable Interface
,”
J. Fluid Mech.
,
482
, p.
319
319
.
16.
Delhaye
,
J. M.
,
1974
, “
Jump Conditions and Entropy Sources in Two-Phase Systems. Local Instant Formulation
,”
Int. J. Multiphase Flow
,
1
, p.
305
305
.
17.
Lakehal
,
D.
,
Meier
,
M.
, and
Fulgosi
,
M.
,
2002
, “
Interface Tracking Towards the Direct Simulation of Heat and Mass Transfer in Multiphase Flows
,”
Int. J. Heat Fluid Flow
,
23
, p.
242
242
.
18.
Temam, R., 1979, Navier-Stokes Equations, North-Holland, Amsterdam.
19.
Tiselj
,
I.
,
Bergant
,
R.
,
Mavko
,
B.
,
Bajsic
,
I.
, and
Hetsroni
,
G.
,
2001
, “
DNS of Turbulent Heat Transfer in Channel Flow With Heat Conduction in the Solid Wall
,”
ASME J. Heat Transfer
,
123
, p.
849
849
.
20.
Phillips, O. M., 1997, The Dynamics of the Upper Ocean, Cambridge University Press, Cambridge, UK.
21.
Tiselj
,
I.
,
Pogrebnyak
,
E.
,
Li
,
C.
,
Mosyak
,
A.
, and
Hetsroni
,
G.
,
2001
, “
Effect of Wall Boundary Condition on Scalar Transfer in a Fully Developed Turbulent Flume
,”
Phys. Fluids
,
13
(
4
), p.
1028
1028
.
22.
Papavassiliou
,
D. V.
, and
Hanratty
,
T. J.
,
1997
, “
Transport of a Passive Scalar in a Turbulent Channel Flow
,”
Int. J. Heat Mass Transfer
,
40
, p.
1303
1303
.
23.
Banerjee
,
S.
,
1971
, “
A Note on Turbulent Mass Transfer at High Schmidt Numbers
,”
Chem. Eng. Sci.
,
26
, p.
989
989
.
24.
Lam
,
K.
, and
Banerjee
,
S.
,
1992
, “
On the Condition of Streaks Formation in a Bounded Turbulent Flow
,”
Phys. Fluids A
,
4
, p.
306
306
.
25.
Nagano, Y., 2002, “Modelling Heat Transfer in Near-Wall Flows,” in Closure Strategies for Turbulent and Transitional Flows, B. Launder and N. Sandham, eds., Cambridge University Press.