Three-Dimensional Hydrostatic Curved Channel Flow Simulations Using Non-Staggered Triangular Grids

Non-staggered triangular grids have many advantages in performing river or ocean modeling with the finite-volume method. However, horizontal divergence errors may occur, especially in large-scale hydrostatic calculations with centrifugal acceleration. This paper proposes an unstructured finite-volume method with a filtered scheme to mitigate the divergence noise and avoid further influencing the velocities and water elevation. In hydrostatic pressure calculations, we apply the proposed method to three-dimensional curved channel flows. Approximations reduce the numerical errors after filtering the horizontal divergence operator, and the approximation is second-order accurate. Numerical results for the channel flow accurately calculate the velocity profile and surface elevation at different Froude numbers. Moreover, secondary flow features such as the vortex pattern and its movement along the channel sections are also well captured.

Stable channel flow with spanwise heterogeneous surface temperature

Recent studies have demonstrated that large secondary motions are excited by surface roughness with dominant spanwise length scales of the order of the flow's outer length scale. Inspired by this, we explore the effect of spanwise heterogeneous surface temperature in weakly to strongly stratified closed channel flow (at $Ri_\tau =120$ , 960; $Re_\tau = 180$ , 550) with direct numerical simulations. The configuration consists of equally sized strips of high and low temperature at the lower and upper boundaries, while an overall stable stratification is induced by imposing an average temperature difference between the top and bottom. We consider the influence of the width of the strips ( ${\rm \pi} /8 \leq \lambda /h \leq 4{\rm \pi} $ ), Reynolds number, stability and upper boundary condition on the mean flow structure, skin friction and heat transfer. Results indicate that secondary flows are excited, with alternating high- and low-momentum pathways and vortices, similar to the patterns induced by spanwise heterogeneous surface roughness. We find that the impact of the surface heterogeneity on the outer layer depends strongly on the spanwise heterogeneity length scale of the surface temperature. Comparison to stable channel flow with uniform temperature reveals that the heterogeneous surface temperature increases the global friction coefficient and reduces the global Nusselt number in most cases. However, for the high-Reynolds cases with $\lambda /h \geq {\rm \pi} /2$ , we find a reduction of the friction coefficient. At stronger stability, the vertical extent of the vortices is reduced and the impact of the heterogeneous temperature on momentum and heat transfer is smaller.

Direct numerical simulation of turbulent mass transfer at the surface of an open channel flow

We present direct numerical simulation results of turbulent open channel flow at bulk Reynolds numbers up to 12 000, coupled with (passive) scalar transport at Schmidt numbers up to 200. Care is taken to capture the very large-scale motions which appear already for relatively modest Reynolds numbers. The transfer velocity at the flat, free surface is found to scale with the Schmidt number to the power ‘ $-1/2$ ’, in accordance with previous studies and theoretical predictions for uncontaminated surfaces. The scaling of the transfer velocity with Reynolds number is found to vary, depending on the Reynolds number definition used. To compare the present results with those obtained in other systems, we define a turbulent Reynolds number at the edge of the surface-influenced layer. This allows us to probe the two-regime model of Theofanous et al. (Intl J. Heat Mass Transfer, vol. 19, 1976, pp. 613–624), which is found to correctly predict that small-scale vortices significantly affect the mass transfer for turbulent Reynolds numbers larger than 500. It is further established that the root mean square of the surface divergence is, on average, proportional to the mean transfer velocity. However, the spatial correlation between instantaneous surface divergence and transfer velocity tends to decrease with increasing Schmidt number and increase with increasing Reynolds number. The latter is shown to be caused by an enhancement of the correlation in high-speed regions, which in turn is linked to the spatial distribution of surface-parallel vortices.