Direct numerical simulation of forced thermal convection in square ducts up to $Re_\tau \approx 2000$

TL;DR

This study uses direct numerical simulation to analyze turbulent heat transfer in square ducts at high Reynolds numbers, confirming the validity of the hydraulic diameter and revealing detailed flow and temperature behaviors.

Contribution

It provides the first high-Reynolds-number DNS data for turbulent square duct flow with heat transfer, validating the hydraulic diameter and analyzing turbulence and heat flux modeling.

Findings

Mean temperature exhibits logarithmic layers with a different slope than velocity.

Secondary motions contribute about 5% to heat transfer.

Traditional linear models for turbulent heat flux show shortcomings.

Abstract

We carry out direct numerical simulation (DNS) of flow in a turbulent square duct by focusing on heat transfer effects, considering the case of unit Prandtl number. Reynolds numbers up to $Re_\tau \approx 2000$ are considered which are much higher than in previous studies, and which yield clear scale separation between inner- and outer-layer dynamics. Close similarity between the behavior of the temperature and the streamwise velocity fields is confirmed as in previous studies related to plane channels and pipes. Just like the mean velocity, the mean temperature is found to exhibit logarithmic layers as a function of the nearest wall, however with a different slope. The most important practical implication is the validity of the traditional hydraulic diameter as the correct reference length for reporting heat transfer data, as we rigorously show here. Temperature and velocity fluctuations also have similar behavior, but apparently logarithmic growth of their inner-scaled peak variances is not observed here unlike in canonical wall-bounded flows. Analysis of the split contributions to the heat transfer coefficient shows that mean cross-stream convection associated with secondary motions is responsible for about $5\%$ of the total. Finally, we use the DNS database to highlight shortcomings of traditional linear closures for the turbulent heat flux, and show that substantial modeling improvement may be in principle obtained by retaining at least the three terms in the vector polynomial integrity basis expansion.