1. Introduction
As fluids flow over solid walls, the transfer of internal energy and momentum across the fluid–solid interface is expected to be similar, assuming identical boundary conditions and driving sources (e.g. pressure gradient and heat source) and comparable molecular momentum and thermal diffusivities. This similarity is known as the Reynolds analogy, named after Reynolds’ pioneering observations on pipe-flow heat transfer (Reynolds Reference Reynolds1874). Although the ideal conditions for this analogy are rare in practical applications, understanding its breakdown in simplified scenarios can improve predictive heat-transfer models and flow-control strategies, particularly when aiming for dissimilar control of momentum and heat transfer (Hasegawa & Kasagi Reference Hasegawa and Kasagi2011). This work focuses on cases where temperature differences are small enough that buoyancy effects are negligible, allowing temperature to be treated as a passive scalar field. Thus, temperature and scalar field will be used interchangeably.
An analogy between momentum and heat transfer might suggest a direct similarity between the velocity and temperature fields. However, this comparison is formally inappropriate due to their different nature: velocity is a three-dimensional solenoidal vector field, while temperature is a scalar field. As pointed out by Abe, Antonia & Kawamura (Reference Abe, Antonia and Kawamura2009), a more sound approach consists in seeking similarity between the turbulent kinetic energy
$\overline {\kappa } = \overline {u_i^\prime u_i^\prime }/2$
, where
$u_i$
denotes the
$i$
th velocity component, and the scalar-fluctuation variance
$\overline {\gamma } = \overline {\theta ^{\prime 2}}/2$
(a prime superscript indicates temporal fluctuations, whereas an over-bar stands for time averaging). These two quantities, both scalar fields governed by similar equations (Corrsin Reference Corrsin1952), are good candidates to assess momentum and scalar similarities within the flow field.
Their qualitative similarity is well-established for flows that closely adhere to the Reynolds analogy, such as smooth-wall turbulent channel flows (Kim & Moin Reference Kim and Moin1989; Kasagi, Tomita & Kuroda Reference Kasagi, Tomita and Kuroda1992; Antonia, Abe & Kawamura Reference Antonia, Abe and Kawamura2009; Pirozzoli, Bernardini & Orlandi Reference Pirozzoli, Bernardini and Orlandi2016). In contrast, departures from analogy can manifest themselves through dissimilar distributions of velocity and scalar fluctuations. For instance, for turbulent flows over surface roughness, Leonardi et al. (Reference Leonardi, Orlandi, Djenidi and Antonia2015); Peeters & Sandham (Reference Peeters and Sandham2019); and Hantsis & Piomelli (Reference Hantsis and Piomelli2022) have shown that roughness acts differently on the momentum and scalar turbulent fluxes, thereby modifying the production of velocity and scalar fluctuations within the roughness sublayer.
Although departures from the analogy can be observed through differing distributions of
$\overline {\kappa }$
and
$\overline {\gamma }$
, linking this observation to a quantitative measure of the Reynolds analogy remains challenging. In wall-bounded flows, the analogy is measured by the Reynolds-analogy factor
$RA=2C_h/C_f$
, with
$RA=1$
indicating perfectly similar momentum and scalar transfer (here
$C_f$
and
$C_h$
are the friction coefficient and Stanton number, respectively). Given that the total momentum and scalar fluxes at the wall, and thus
$C_f$
and
$C_h$
, must be related to the total power input required to sustain the flow, inclusive of mean and turbulent contributions, it seems short-sighted to seek a direct connection between
$C_f$
and
$C_h$
and turbulent energy fields (e.g.
$\overline {\kappa }$
and
$\overline {\gamma }$
).
Based on this argument, we propose to consider the Reynolds-analogy problem in terms of the total mechanical energy and half of the squared scalar field,
$K=u_iu_i/2$
and
$\Gamma =\theta ^2/2$
, respectively. Their governing equations can be written in the form
\begin{align} \frac {\partial \mathit {K}}{\partial t} + u_i \frac {\partial B}{\partial x_i} &= \nu \frac {\partial ^2 B}{\partial x_i \partial x_i} - \nu \omega _i \omega _i + f_i u_i, \end{align}
\begin{align} \frac {\partial \mathit {\Gamma }}{\partial t} + u_i \frac {\partial \mathit {\Gamma }}{\partial x_i} &= \alpha \frac {\partial ^2 \mathit {\Gamma }}{\partial x_i \partial x_i} - \alpha \frac {\partial \theta }{\partial x_i} \frac {\partial \theta }{\partial x_i} + q\theta , \end{align}
where
$B=K+(p/\rho)$
is the total mechanical energy,
$p,\, \rho$
,
$\nu$
, and
$\alpha$
are, respectively, the pressure, fluid density, kinematic viscosity and thermal diffusivity,
$\omega _i$
is the
$i$
th vorticity component, and
$f_i$
and
$q$
are, respectively, external volume-force and scalar-source terms. Equation (1.1a
) is obtained by taking the inner product of the velocity with the momentum equation, and recalling that
$p,_{ii}/\rho = \omega _i \omega _i/2 - s_{ij}s_{ij}$
(Bradshaw & Koh Reference Bradshaw and Koh1981). The notation
$(\cdot),_i$
indicates differentiation in the
$i$
direction. Here
$s_{ij}=\ ({1}/{2}) ( u_{i,j} + u_{j,i} )$
is the strain-rate tensor. Similarly, (1.1b
) is obtained by multiplying the transport equation for
$\theta$
by
$\theta$
itself.
Equations (1.1a
) and (1.1b
) highlight a key difference in the evolution of
$K$
and
$\Gamma$
: the point-wise time rate of change of
$K$
depends on the transport and molecular diffusion of total mechanical energy
$B$
, while that of
$\Gamma$
depends solely on the transport and diffusion of
$\Gamma$
itself. For statistically stationary flows, time averaging the equations yields two steady-state advection–diffusion equations for
$\overline {B}$
and
$\overline {\Gamma }$
, in which the terms
$\nu \overline {\omega _i \omega _i}$
and
$\alpha \overline {\theta ,_i\theta ,_i}$
play the role of dissipation rates of
$\overline {B}$
and
$\overline {\Gamma }$
, respectively. Note that, according to (1.1a
), the former occurs through the action of the local dissipation rate of kinetic energy,
$\nu \omega _i \omega _i$
(Wu, Zhou & Fan Reference Wu, Zhou and Fan1999). For statistically stationary flows,
$\nu \overline {\omega _i \omega _i}$
amounts also to the local dissipation rate of the time-averaged total mechanical energy,
$\overline {B}$
; this equivalence holds because the balance of
$\overline {B}$
differs from that of
$\overline {K}$
only through the pressure transport terms of convection and molecular diffusion. The integral form of (1.1a
) and (1.1b
) establishes a link between the total power necessary to sustain the flow and field quantities associated with the mechanical and scalar energy. This connection can be further extended to include the total momentum and scalar fluxes at the wall, thus determining a bridge to the Reynolds-analogy factor.
This work explores this approach to study fully developed turbulent channel flow over homogeneous random surface roughness. Surface roughness impacts scalar and momentum transfer depending on the relative size of roughness elements to the viscous sublayer thickness imposed by the flow. In the hydraulically smooth regime, where roughness elements are significantly smaller than the viscous sublayer thickness, no roughness effects are observed. Heat and momentum transfer occur primarily via molecular diffusion across the channel walls, maintaining similar transfer mechanisms. As the roughness size approaches that of the viscous sublayer, the flow enters the transitionally rough regime. Here, the momentum transfer departs from an analogous behaviour due to the momentum loss caused by the form-drag induced by roughness elements, which lacks an equivalent mechanism for the heat transfer. Ultimately, in the fully rough regime, the form-drag becomes the predominant contribution to momentum loss, whereas there is not yet unanimous consensus on a scalar transfer law in this regime (see e.g. the recent review of Hantsis & Piomelli Reference Hantsis and Piomelli2024).
In the present work, direct numerical simulations (DNS) are used to replicate two selected cases presented by Thakkar, Busse & Sandham (Reference Thakkar, Busse and Sandham2017); Thakkar (Reference Thakkar2017); Thakkar, Busse & Sandham (Reference Thakkar, Busse and Sandham2018) and Peeters & Sandham (Reference Peeters and Sandham2019). In this context, (1.1a ) and (1.1b ) provide a novel perspective on the breaking of the Reynolds analogy, enabling both qualitative and quantitative assessments of the phenomenon.
This work is organized as follows: in § 2, (1.1a ) and (1.1b ) are specialized for the flow under investigation; the details about the simulation set-ups and the numerical procedure are presented in § 3; results are reported in § 4, while concluding remarks are summarized in § 5.
2. Problem formulation
We consider a fully developed turbulent channel flow with mean streamwise velocity directed along the
$x=x_1$
direction. The wall-normal direction is denoted by
$z=x_3$
. No-slip and non-penetration boundary conditions are applied at the stationary walls of the channel, where the passive scalar
$\theta$
is fixed to a constant value
$\theta _w$
. To have analogous wall boundary conditions for the velocity and the passive scalar field, we set
$\theta _w=0$
. Upon non-dimensionalization using viscous units, the time-averaged equations (1.1a
) and (1.1b
) become
\begin{align} \overline {u_i^+ \frac {\partial B^+}{\partial x_i^+}} &= \frac {\partial ^2 \overline {B}^+}{\partial x_i^+ \partial x_i^+} - \overline {\omega _i^+ \omega _i^+} + \frac {1}{Re_\tau } \delta _{1i} \overline {u_i}^{\;+}, \end{align}
\begin{align} \frac{1}{Pr}\overline{u_i^+ \frac{\partial \mathit{\Gamma}^+}{\partial x_i^+}} &= \frac{1}{Pr^2}\frac{\partial^2 \overline{\mathit{\Gamma}}^+}{\partial x_i^+ \partial x_i^+} - \frac{1}{Pr^2}\overline{\frac{\partial \theta^+}{\partial x_i^+} \frac{\partial \theta^+}{\partial x_i^+}} + \frac{1}{Pe_\tau}\overline{\theta}^{\;+}, \end{align}
where
$Re_\tau =u_\tau \delta /\nu$
,
$Pe_\tau =Pr Re_\tau$
. The friction velocity is
$u_\tau =(\tau _w/\rho )^{1/2}$
,
$\tau _w$
being the mean total wall-shear stress and
$\delta$
the half-channel height. Note that the Prandtl number,
$Pr=\nu /\alpha$
, can be replaced by the Schmidt number in case
$\theta$
represented a species-concentration field. The friction passive scalar is defined as
$\theta _\tau =\alpha q_w/u_\tau$
, where
$q_w$
is the mean scalar-flux at the wall,
$q_w = \ ({1}/{A})\int _{A_w}\theta ,_i n_i \textrm{d}A$
(where
$A_w$
and
$A$
are, respectively, the wetted and plan area of the channel wall, and
$n_i$
indicates the unit vector normal to the wall surface). The
$ + $
superscript is used to denote quantities non-dimensionalised using
$\delta_{\nu}=\nu /u_{\tau}$
,
$u_{\tau}$
and
$ \theta_\tau $
.
The mean bulk velocity
$U_b^+$
and scalar
$\Theta _b^+$
can be computed by integrating, respectively, (2.1a
) and (2.1b
) over the channel volume
$\mathcal {V}$
and dividing the result by the channel volume. Namely,
\begin{align} U_b^+ & = \frac {Re_\tau }{\mathcal {V}^+} \int _{\mathcal {V}^+} \overline {\omega _i^+ \omega _i^+} \textrm{d} \mathcal {V}^+, \end{align}
\begin{align} \Theta_b^+ & = \frac{Pe_\tau}{Pr^2 \mathcal{V}^+}\int_{\mathcal{V}^+} \overline{\frac{\partial \theta^+}{\partial x_i^+} \frac{\partial \theta^+}{\partial x_i^+}} d \mathcal{V}^+. \end{align}
Equation (2.2a
) states the integral balance between the total mechanical power exerted by the flow, i.e.
$U_b^+$
, and the integral dissipation of mechanical energy on the right-hand side of the equation. Analogously, (2.2b
) is interpreted as the integral balance between the mean bulk scalar,
$\Theta _b^+$
, and the integral dissipation of
$\overline {\Gamma }$
. Note that these expressions are equivalent to the integral energy and scalar balances presented, for smooth-wall flows, by Laadhari (Reference Laadhari2007); Abe & Antonia (Reference Abe and Antonia2016, Reference Abe and Antonia2017); Wei (Reference Wei2018). Yet, (2.2a
) and (2.2b
) remain valid even in the case of non-smooth-wall flows and regardless of the decomposition into mean and turbulent flow contributions.
Expressions (2.2a
) and (2.2b
) can be inserted in the definitions of the friction coefficient
$C_f = 2/U_b^{+2}$
and of the Stanton number
$C_h=1/(U_b^+ \Theta _b^+)$
, which, in turn, can be used to compute the Reynolds-analogy factor
$RA = 2C_h / C_f$
:
\begin{equation} RA = \frac{U_b^+}{\Theta_b^+} = Pr \left( \int_{\mathcal{V}^+} \overline{\omega_i^+ \omega_i^+} d \mathcal{V}^+\right)\times \left( \int_{\mathcal{V}^+} \overline{\frac{\partial \theta^+}{\partial x_i^+}} \overline{\frac{\partial \theta^+}{\partial x_i^+}} d \mathcal{V}^+ \right)^{-1}. \end{equation}
This expression relates the Reynolds-analogy factor to quantities defined within the bulk of the flow and, as such, can be instrumental in unveiling flow phenomena that lead to breaking the analogy. It must be emphasized that, locally, departures from the Reynolds analogy can manifest themselves through dissimilar distributions of all the terms in (2.1a
) and (2.1b
); however, from an intergal perspective, (2.3) shows that a dissimilar scalar and momentum transfer can be quantified globally through the volume integrals of the dissipation rates of
$\overline {B}$
and
$\overline {\Gamma }$
.
3. Case set-ups and numerical procedure
The present numerical simulations replicate two cases of the rough-wall fully developed turbulent channel-flow simulations presented by Thakkar et al. (Reference Thakkar, Busse and Sandham2017), Thakkar (Reference Thakkar2017), Thakkar et al. (Reference Thakkar, Busse and Sandham2018) and Peeters & Sandham (Reference Peeters and Sandham2019). In particular, in one case the flow is transitionally rough, but close to smooth-wall conditions, with
$k^+=15$
(here
$k^+=k u_\tau /\nu$
, and
$k$
is the mean peak-to-trough height of the roughness topography), whereas for the second case, the flow is close to fully rough conditions, with
$k^+=90$
. The roughness topography is obtained from a scan of a real grit-blasted surface and is statistically homogeneous in the wall-parallel directions
$x_1=x$
and
$x_2=y$
. A sample of the roughness topography is shown in figure 1. The physical size of the roughness is
$k=\delta /12$
, for the
$k^+=15$
case, and
$k=\delta /6$
, for the
$k^+=90$
case. The friction Reynolds number is
$Re_\tau =180$
for the former case and
$Re_\tau =540$
for the latter. An additional smooth-wall channel flow at
$Re_\tau =180$
is also simulated. For all the cases, a unit Prandtl number is considered.
Sample of the roughness topography for the case
$k^+=15$
(a). Details of the computational grid:
$k^+=15$
case (b);
$k^+=90$
case (c).
The numerical integration of the incompressible Navier–Stokes equations is carried out using the flow solver
$\texttt{Nek5000}$
(Fisher, Lottes & Kerkemeier Reference Fisher, Lottes and Kerkemeier2019), which is based on the spectral-element method proposed by Patera (Reference Patera1984). Surface roughness is introduced into the simulations using a body-fitted mesh, obtained by deforming, in the wall-normal direction
$z$
, the spectral elements of a plane-channel box having the bottom and top walls at, respectively,
$z=0$
and
$z=2\delta$
. The mesh deformation is symmetric with respect to the channel centre line and varies linearly with the wall-normal distance from the bottom wall of the non-deformed channel. As a result, the mean height of the deformed channel is
$2\delta$
. Figures 1(b) and 1(c) illustrate the details of the body-conforming computational meshes used in the two rough-wall cases.
The non-deformed computational domain has size
$L_x\times L_y\times 2\delta$
in the streamwise, spanwise and wall-normal directions, respectively. Correspondingly, the number of elements in the three coordinate directions are denoted by
$n_x$
,
$n_y$
and
$n_z$
. A uniform distribution of elements is used in the streamwise and spanwise directions, whereas, in the wall-normal direction, elements are clustered near the channel walls. The domain size, the number of elements in each direction, the polynomial degree of the solution, and the total number of degrees of freedom (DOF) are reported in table 1. The total number of DOF corresponds to the number of grid points in the physical space and is given by
$n_x\times n_y \times n_z \times (P+1)^3$
, since each element is divided into
$P+1$
grid points. The points distribution along each element direction is non-uniform and follows the Gauss–Legendre–Lobatto distribution of quadrature nodes. This can be appreciated in the panels (b) and (c) of figure 1 for the case
$P=7$
.
Simulation parameters. Here,
$P$
is the polynomial degree of the solution, whereas DOF indicates the total number of degrees of freedom.
The present numerical set-up is validated against the reference data of Thakkar (Reference Thakkar2017) and Peeters & Sandham (Reference Peeters and Sandham2019). As an example, figure 2 shows a comparison of the mean streamwise-velocity and temperature profiles with the available reference data. In the figure and throughout the text, angular brackets denote intrinsic averaging, i.e. spatial averaging along the homogeneous channel directions performed in fluid regions only; namely,
$\left \lt \cdot \right \gt = ( \int _{A_f} \cdot \, \textrm {d}A_f )/A_f$
, where
$A_f$
indicates the plan area covered by the fluid at a given wall-normal height. Time-statistics reported in this work have been collected for different time integration intervals,
$T_{avg}$
, depending on the considered case. For each simulated case,
$T_{avg}$
, expressed in terms of external time units,
$\delta /U_b$
, is reported in table 1. For the
$k^+=90$
case, a significantly shorter integration time interval was used, in comparison with the other cases, to save on computational resources. For this case, the convergence of the statistics was checked by comparing the results for
$T_{avg}=140.9\delta /U_b$
and
$T_{avg}=94.7\delta /U_b$
and no significant differences were observed in the acquired statistics.
Comparison with data from Thakkar (Reference Thakkar2017) and Peeters & Sandham (Reference Peeters and Sandham2019). (a) Mean streamwise-velocity profiles; (b) mean scalar profiles. Black dashed lines indicate smooth-wall data at
$Re_\tau =180$
.
4. Results
Compared with a smooth-wall flow at matched
$Re_\tau$
, a rough-wall flow exhibits increased momentum and scalar transfer which manifest themselves as a downward shift of the logarithmic regions of the mean streamwise-velocity and scalar profiles (see, for instance, the review of Chung et al. Reference Chung, Hutchins, Schultz and Flack2021). This can be appreciated in figure 2, where the mean velocity and scalar profiles are shown for all the simulated cases. From an integral perspective, the roughness effect results in a decrease in the mean bulk velocity and scalar,
$U_b^+$
and
$\Theta _b^+$
, respectively. The latter, computed using (2.2a
) and (2.2b
), and the Reynolds-analogy factor (2.3) are reported in table 1. As expected, as the roughness Reynolds number,
$k^+$
, increases, so does the departure from the Reynolds analogy. This is shown by the faster decrease of
$U_b^+$
compared with
$\Theta _b^+$
for increasing
$k^+$
. Correspondingly, the Reynolds-analogy factor decreases, from
$RA=1$
for the smooth-wall case, to
$RA=0.8$
for the
$k^+=90$
case. As expected, while surface roughness promotes momentum and scalar exchange, the former is much favoured compared with the latter. Further, (2.2a
) indicates that
$U_b^+$
and
$\Theta _b^+$
represent, respectively, the total mechanical and scalar power realized by the flow for a specified mean pressure gradient and heat source. Hence, at fixed
$Re_\tau$
, an increase in
$k^+$
shows as a decrease in the total power that can be dissipated within the flow.
4.1. Mechanical- and scalar-energy budgets
The intrinsic average of the transport, viscous diffusion, dissipation and power-input terms of (2.1a
) and (2.1b
) are shown in figure 3. The definition of each term is reported in table 2. The power-input terms,
$\Pi$
and
$\Pi _\theta$
, correspond to the mean streamwise-velocity and scalar profiles scaled by the Reynolds and Péclet numbers, respectively. Hence, the observations, made in reference to the mean profiles in figure 2, hold when comparing
$\Pi$
and
$\Pi _\theta$
. Figure 3 further allows to relate the power input to the other terms in the budgets. For all the cases, towards the channel centre line, the viscous diffusion and the dissipation terms become negligible compared with the other terms, and
$\Pi \approx T$
and
$\Pi _\theta \approx T_\theta$
. Closer to the wall, all terms in the budgets become relevant and each case needs to be addressed separately.
Mechanical- and scalar-energy budgets. (a) Smooth wall,
$Re_\tau =180$
; (b)
$k^+=15$
,
$Re_\tau =180$
; (c)
$k^+=90$
,
$Re_\tau =540$
. 
,
$-\varepsilon$
; 
,
$-\varepsilon _\theta$
; 
,
$-T$
; 
,
$-T_\theta$
; 
,
$D$
; 
,
$D_\theta$
; 
,
$\Pi$
; 
,
$\Pi _\theta$
. The vertical black dashed lines in (b) and (c) indicate the boundaries of the roughness canopy.
For the smooth-wall case, an almost perfect overlap is observed everywhere between the terms in the mechanical- and the scalar-energy equations. It should be noted that, at higher Reynolds numbers, such a perfect correspondence might be lost. As shown by Abe & Antonia (Reference Abe and Antonia2017), at
$Pr=1$
, and for sufficiently high Reynolds numbers,
$U_b^+ = 2.54 \log (Re_\tau ) + 2.41$
and
$\Theta _b^+ = 2.18 \log (Re_\tau ) + 4$
. These two relations predict diverging
$U_b^+$
and
$\Theta _b^+$
as the Reynolds number increases, with the scalar transfer becoming more effective than momentum transfer for sufficiently large Reynolds numbers (see below, figure 4
a and the related discussion). Hence, differences can appear also in the budget terms shown in the figure.
For the
$k^+=15$
case, figure 3(b) shows that roughness produces changes in the mechanical- and scalar-energy budget terms only well within the roughness canopy, leaving the distributions practically unaltered farther from the wall. Within the roughness canopy, the viscous diffusion is almost exclusively balanced by the dissipation term for both the mechanical and scalar budgets. Further,
$\varepsilon _\theta$
and
$D_\theta$
are amplified compared with
$\varepsilon$
and
$D$
. This scenario differs for the
$k^+=90$
case in figure 3(c). Within the roughness canopy, the largest differences occur, whereas, away from the wall, the terms of the mechanical- and scalar-energy budgets tend to overlap, with the exception of
$T_\theta$
and
$\Pi _\theta$
which remain larger in magnitude than
$T$
and
$\Pi$
, respectively. Focusing on the mechanical-energy budget, the transport term
$T$
penetrates deep in the roughness sublayer , where it exceeds in magnitude the viscous diffusion
$D$
. Concurrently, the mechanical-energy dissipation exceeds, in magnitude, the joint contribution of
$T$
and
$D$
, thus realizing a slightly negative power input for
$z\lesssim -0.045\delta$
. The same cannot occur for the scalar power input
$\Pi _\theta$
, for the scalar
$\theta$
, and thus
$\Pi _\theta$
, is bounded to vary in the interval
$ [0, \, 1 ]$
. Consistently, within the roughness canopy,
$\varepsilon_\theta$
retains a smaller magnitude compared with the joint contribution of
$T_\theta$
and
$D_\theta$
, therefore realizing a small, positive,
$\Pi _\theta$
. Interestingly, the profiles of
$T_\theta$
,
$D_\theta$
and
$\varepsilon _\theta$
are qualitatively similar in shape to the corresponding mechanical-energy budget terms, even though they are larger in magnitude.
Definition of time- and space-averaged terms of the budget equations for
$\overline {B}$
and
$\overline {\Gamma }$
.
(a) Mean bulk velocity (blue) and scalar (red) as functions of the friction Reynolds number. Smooth wall, circle markers; rough wall
$k^+=15$
, cross markers; rough wall
$k^+=90$
, square markers. (b) Mean and turbulent contributions to the volume-averaged dissipation rates. 
,
$Re_\tau \mathcal{E}^m$
; 
,
$Re_\tau \mathcal{E}^t$
; 
,
$Pe_\tau \mathcal{E}_\theta^m/(Pr^2)$
; 
,
$Pe_\tau \mathcal{E}_\theta^t/(Pr^2)$
.
4.2. Mechanical and scalar dissipation
The bulk mechanical and scalar dissipation, denoted respectively by
$\mathcal {E}$
and
$\mathcal {E}_\theta$
, are defined as the volume average of the dissipation rates
$\overline {\omega _i^+\omega _i^+}$
and
$\overline {\theta ,_i^+ \theta ,_i^+}$
. Their relationship to the mean bulk velocity and scalar is given by (2.2a
) and (2.2b
). Figure 4(a) shows the computed values of
$U_b^+$
and
$\Theta _b^+$
as functions of the friction Reynolds number. The plot shows also the logarithmic distributions expected for smooth-wall channel flow (Abe & Antonia Reference Abe and Antonia2016, Reference Abe and Antonia2017). The increase in momentum and scalar transfer introduced by the roughness causes the rough-wall data to lie below the smooth-wall asymptotes. The downward shift of
$U_b^+$
and
$\Theta _b^+$
might rest upon different mean and turbulent flow contributions. This is verified by using the Reynolds decomposition to split
$\mathcal {E}$
and
$\mathcal {E}_\theta$
into their mean,
$\mathcal {E}^m$
,
$\mathcal {E}_\theta ^m$
, and turbulent contributions,
$\mathcal {E}^t$
,
$\mathcal {E}_\theta ^t$
, respectively.
The contributions to
$U_b^+$
and
$\Theta _b^+$
of the mean and turbulent field are shown in figure 4(b) for the smooth- and rough-wall cases. Note that the smooth-wall data at
$Re_\tau =540$
are computed using the asymptotic relationships given by Abe & Antonia (Reference Abe and Antonia2016, Reference Abe and Antonia2017). At
$k^+=15$
, the mean-flow contribution exceeds that due to the stochastic turbulent field. It is interesting to observe that
$\mathcal {E}^t$
,
$\mathcal {E}_\theta ^t$
attain similar values, which are also close to the smooth-wall data at the same Reynolds number. Conversely,
$\mathcal {E}^m$
and
$\mathcal {E}_\theta ^m$
are reduced from the corresponding smooth-wall values, and
$\mathcal {E}_\theta ^m$
is approximately
$10\%$
larger than
$\mathcal {E}^m$
. Compared with the smooth-wall predictions at
$Re_\tau =540$
, the rough-wall data at
$k^+=90$
shows that
$\mathcal {E}^m$
suffers the largest decrease, whereas
$\mathcal {E}^t$
, still being significantly lower than its corresponding smooth-wall counterpart, exceeds
$\mathcal {E}^m$
. On the other hand,
$\mathcal {E}_\theta ^m$
and
$\mathcal {E}_\theta ^t$
appear similar to each other.
Interestingly, figure 4(b) shows that, for both rough-wall cases, the largest contribution to the breaking of the Reynolds analogy appears to be caused by the mean-field: at
$k^+=15$
,
$\mathcal {E}_\theta ^m\approx 1.10\mathcal {E}^m$
and
$\mathcal {E}_\theta ^t\approx \mathcal {E}^t$
; at
$k^+=90$
,
$\mathcal {E}_\theta ^m\approx 1.96\mathcal {E}^m$
and
$\mathcal {E}_\theta ^t\approx 0.93 \mathcal {E}^t$
. Even though the limited amount of data are not sufficient to confirm this trend, it is instructive to further investigate the mean-flow contribution to
$U_b^+$
and
$\Theta _b^+$
. To do so, the time-averaged flow field is further decomposed in its time- and space-averaged parts and its dispersive (or wake) contribution
$\overline {q} = \left \lt \overline {q}\right \gt + \overline {q}^{\prime \prime }$
, where
$(\cdot )^{\prime \prime }$
indicates the dispersive fluctuation of
$\overline {q}$
. The contributions to
$\mathcal {E}^m$
and
$\mathcal {E}_\theta ^m$
of, respectively,
$\overline {\omega _i}^{\prime \prime +}\overline {\omega _i}^{\prime \prime +}$
and
$\overline {\theta ,_i}^{\prime \prime +}\overline {\theta ,_i}^{\prime \prime +}$
are shown in figure 4(b) as hatched areas laid on top of the
$\mathcal {E}^m$
and
$\mathcal {E}_\theta ^m$
boxes. At
$k^+=15$
, the dispersive components of the dissipation rates contribute little to
$\mathcal {E}^m$
and
$\mathcal {E}_\theta ^m$
and their magnitude is very similar between the mechanical and scalar dissipation rates. As expected, this scenario changes significantly at
$k^+=90$
, as the dispersive contribution to
$\mathcal {E}^m$
and
$\mathcal {E}_\theta ^m$
becomes the prominent one and
$\overline {\theta ,_i}^{\prime \prime +}\overline {\theta ,_i}^{\prime \prime +}$
significantly exceeds
$\overline {\omega _i}^{\prime \prime +}\overline {\omega _i}^{\prime \prime +}$
.
Dispersive dissipation rates for the
$k^+=15$
(a,b) and
$k^+=90$
(c,d) cases. Panels (a,c) show
$\overline {\omega _i}^{\prime \prime +}\overline {\omega _i}^{\prime \prime +}$
; (b,d)
$\overline {\theta ,_i}^{\prime \prime +}\overline {\theta ,_i}^{\prime \prime +}$
. Dashed lines represent lines of constant
$\overline {\omega _y}^{\prime \prime +}=\pm 0.046$
(a,c) and
$\overline {\theta ,_z}^{\prime \prime +}=\pm 0.046$
(b,d). Blue and red colours indicate, respectively, negative and positive values. In panels (c,d), a black dotted line represents the isocontour line of zero mean streamwise velocity.
Figure 5 compares the distributions of
$\overline {\omega _i}^{\prime \prime +}\overline {\omega _i}^{\prime \prime +}$
and
$\overline {\theta ,_i}^{\prime \prime +}\overline {\theta ,_i}^{\prime \prime +}$
within typical roughness patterns for the two rough-wall cases . High shear near roughness peaks leads to intense local heat and momentum transfer, where both fluxes exhibit a similar behaviour, as noted by Zhong, Hutchins & Chung (Reference Zhong, Hutchins and Chung2023). Locally, the strong, positive, velocity and scalar gradients exceed significantly their spatial mean at the same wall-normal height, thus resulting in intense positive dispersive fluctuations
$\overline {\omega _i}^{\prime \prime +}$
and
$\overline {\theta ,_i}^{\prime \prime +}$
. Correspondingly, the dispersive dissipation fields,
$\overline {\omega _i}^{\prime \prime +}\overline {\omega _i}^{\prime \prime +}$
and
$\overline {\theta ,_i}^{\prime \prime +}\overline {\theta ,_i}^{\prime \prime +}$
, attain their maxima. In figure 5, these regions are identified by isocontour lines for positive values of
$\overline {\omega _y}^{\prime \prime +}$
(figure 5
a,c) and
$\overline {\theta },_z^{\prime \prime +}$
(figure 5
b,d). They grow attached to the windward side of roughness hills and develop, exposed to high shear, up to the hill’s crest, where they eventually detach from the roughness contour. Such regions will be named attached, exposed regions (AER) in the following. In AER, an apparent similarity is observed between
$\overline {\omega _i}^{\prime \prime +}\overline {\omega _i}^{\prime \prime +}$
and
$\overline {\theta ,_i}^{\prime \prime +}\overline {\theta ,_i}^{\prime \prime +}$
.
Large roughness elements tend to shelter their downstream regions from the incoming flow, thereby generating pockets of slow-moving, recirculating flow, which entrain outer fluid towards the bottom of the roughness valleys (Forooghi, Stripf & Frohnapfel Reference Forooghi, Stripf and Frohnapfel2018). These regions are characterized by nearly uniform scalar distributions and smaller scalar gradients compared with their mean value at the same wall-normal height. This results in negative dispersive fluctuations of the scalar gradient,
$\overline {\theta ,_i}^{\prime \prime +}$
, which, in turn, increase the local dispersive dissipation rate
$\overline {\theta ,_i}^{\prime \prime +}\overline {\theta ,_i}^{\prime \prime +}$
. Hereafter, these regions will be referred to as well-mixed regions (WMR). The WMR are also characterized by negative dispersive fluctuations of the velocity gradient components, and hence of
$\overline {\omega _i}^{\prime \prime +}$
, due to the slow-moving, possibly recirculating, flow. In figure 5, WMR are identified by isocontour lines for negative values of
$\overline {\omega _y}^{\prime \prime +}$
(figure 5
a,c) and
$\overline {\theta },_z^{\prime \prime +}$
(figure 5
b,d). It can be seen how WMR extend throughout roughness valleys between large roughness elements. At
$k^+=15$
, WMR exhibit very similar distributions of dispersive dissipation rates. On the contrary, within these regions, significant differences are evident for the
$k^+=90$
case, where
$\overline {\theta ,_i}^{\prime \prime +}\overline {\theta ,_i}^{\prime \prime +}$
appears significantly larger than
$\overline {\omega _i}^{\prime \prime +}\overline {\omega _i}^{\prime \prime +}$
.
The spatial distributions of
$\overline {\omega _i}^{\prime \prime +}\overline {\omega _i}^{\prime \prime +}$
and
$\overline {\theta ,_i}^{\prime \prime +}\overline {\theta ,_i}^{\prime \prime +}$
represent an indication of how similarly the dissipation rates of
$\overline {B}$
and
$\overline {\Gamma }$
respond to the spatial changes imposed by the roughness topography. For instance, figures 5(a) and 5(b) clearly highlight differences in the distributions of
$\overline {\omega _i}^{\prime \prime +}\overline {\omega _i}^{\prime \prime +}$
and
$\overline {\theta ,_i}^{\prime \prime +}\overline {\theta ,_i}^{\prime \prime +}$
within a WMR. Understanding the physical reasons for such dissimilarities requires further detailed analysis of local flow features. In the specific example of figures 5(a) and 5(b), the observed differences are likely related to the presence of separated, recirculating flow in the large roughness valley, in agreement with observations in the literature (for instance, see Spalart & Coleman Reference Spalart and Coleman1997; MacDonald, Hutchins & Chung Reference MacDonald, Hutchins and Chung2019) . The zero-mean-streamwise-velocity contour line in the figure emphasizes the back-flow at the bottom of the valley.
Although the visualization of
$\overline {\omega _i}^{\prime \prime +}\overline {\omega _i}^{\prime \prime +}$
and
$\overline {\theta ,_i}^{\prime \prime +}\overline {\theta ,_i}^{\prime \prime +}$
can provide a qualitative identification of flow regions where the dissipation rates behave dissimilarly, it does not provide an indication of the net contribution that such flow regions have on
$\mathcal {E}^m$
and
$\mathcal {E}_\theta ^m$
. However, this information can be obtained by assessing the mean dissipation rates,
$\overline {\omega _i}^{+}\overline {\omega _i}^{+}$
and
$\overline {\theta ,_i}^{+}\overline {\theta ,_i}^{+}$
, over the identified flow regions. In particular, the relative contributions of AER and WMR to
$\mathcal {E}^m$
and
$\mathcal {E}_\theta ^m$
is estimated by conditional integration of the mean dissipation rates.
The AER and WMR contributions to
$Re_\tau \mathcal{E}^m$
and
$Pe_\tau \mathcal{E}_\theta^m/(Pr^2)$
. (a)
$k^+=15$
; (b)
$k^+=90$
. 
, AER; 
, WMR; 
, sum of AER and WMR contributions; a black solid outline indicates
$Re_\tau \mathcal{E}^m$
and
$Pe_\tau \mathcal{E}_\theta^m/(Pr^2)$
.
This is achieved by integrating
$\overline {\omega _i}^{+}\overline {\omega _i}^{+}$
and
$\overline {\theta ,_i}^{+}\overline {\theta ,_i}^{+}$
over AER and WMR. The former were identified as the regions over which
$\overline {\theta },_z^{\prime \prime +} \gt 0.046$
, whereas for the latter, the condition
$\overline {\theta },_z^{\prime \prime +} \lt -0.046$
was used. This particular choice was made based on the visual inspection of the flow and upon noting that these two conditions correlate well with the presence of AER and WMR, as visualized in figure 5, throughout the roughness canopy.
Figure 6 shows the result of the conditional integration of the dissipation rates, normalized by the empty channel volume. For both cases, the largest contribution to
$\mathcal {E}^m$
and
$\mathcal {E}_\theta ^m$
is introduced by AER, whereas the contribution of WMR, despite being significant at
$k^+=15$
, becomes marginal at
$k^+=90$
. Interestingly, the contribution of WMR is very similar for the mechanical and scalar dissipation rates at both roughness Reynolds numbers. Namely, the dissimilarity between
$\mathcal {E}^m$
and
$\mathcal {E}_\theta ^m$
can be attributed to the differences displayed by the dissipation rates in AER alone. Note that the joint contribution of AER and WMR accounts for a large part of
$\mathcal {E}^m$
and
$\mathcal {E}_\theta ^m$
, especially at
$k^+=90$
, where it represents more than
$90\%$
of
$\mathcal {E}^m$
and
$\mathcal {E}_\theta ^m$
. The residual represents the contribution of the mean-flow dissipation rates that takes place in the core of the channel.
5. Conclusions
The present work explored a novel perspective for the assessment of the Reynolds analogy in turbulent, forced convection. In particular, it is shown that
$B$
(the total mechanical energy, sum of the kinetic energy and pressure work) and
$\Gamma$
(one-half of the squared scalar field) share a similar structure in their time-averaged governing equations. For fully developed channel-flow configurations, the integral balance of total mechanical and scalar energy allows expressing the mean bulk velocity
$U_b^+$
and scalar
$\Theta _b^+$
in terms of the volume averages of the dissipation fields
$\overline {\omega _i^+\omega _i^+}$
and
$\overline {\theta ,_i^+ \theta ,_i^+}$
, respectively. In turn, the Reynolds-analogy factor (the ratio between the total mechanical and scalar energy realized by the flow, i.e.
$U_b^+/\Theta _b^+$
) can be related to the dissipation fields.
This energetic perspective on the Reynolds-analogy problem enables the identification of flow regions where departures from the analogy occur and provides a means for quantifying their net effect on the momentum and scalar transfer. This approach was showcased on fully developed turbulent channel flow over homogeneous surface roughness, a well-known instance of the failure of the Reynolds analogy. The study presented DNS results of two rough-wall turbulent channel flows at
$k^+=15$
and
$k^+=90$
. In addition, a smooth-wall channel flow at
$Re_\tau =180$
was also considered. For all the cases, a unit Prandtl number was considered.
The present data show that similarity exists between the time- and space-averaged wall-normal distributions of the mechanical and scalar budget-equation terms. For the smooth-wall case, each term overlaps in both equations. For the rough-wall cases, roughness amplifies significantly
$\overline {\Gamma }$
budget-equation terms while maintaining similar shapes of the corresponding terms in the
$\overline {B}$
budget equation.
An increase in
$k^+$
determines a decrease in the total power exerted by the flow; equivalently, the total energy dissipated within the flow decreases as the roughness Reynolds number increases. For the present data, roughness affects more significantly the volume-averaged mean dissipation rates,
$\mathcal {E}^m$
and
$\mathcal {E}_\theta ^m$
, compared with the turbulent contributions
$\mathcal {E}^t$
and
$\mathcal {E}_\theta ^t$
. While at
$k^+=15$
the latter are practically unchanged, compared with their smooth-wall equivalent at matched
$Re_\tau$
, at
$k^+=90$
, the relative change, compared with the smooth-wall reference, of
$\mathcal {E}^m$
,
$\mathcal {E}_\theta ^m$
exceeds that of
$\mathcal {E}^t$
,
$\mathcal {E}_\theta ^t$
. Further, at both roughness Reynolds numbers, the largest departure from Reynolds analogy is observed for the mean-flow dissipation contributions, whereas
$\mathcal {E}^t$
and
$\mathcal {E}_\theta ^t$
remain similar. It is noted that further data are necessary to explore higher Reynolds number scenarios, whereby the turbulent contributions to the dissipation rates possibly start to play a dominant role in shaping the asymptotic behaviours of
$U_b^+$
and
$\Theta _b^+$
.
The analysis of the mean dissipation rates reveals the flow regions that contribute the most to the mean scalar and momentum transfer within the roughness sublayer. These are AER, which develop on the windward side of roughness hills, and WMR, which occur within the valleys of the roughness topography. In both cases, AER contribute most significantly to
$\mathcal {E}^m$
and
$\mathcal {E}_\theta ^m$
, while the net contribution of WMR becomes negligible at
$k^+=90$
. Interestingly, the WMR contribution is nearly identical for both mechanical and scalar dissipation rates for the two investigated roughness Reynolds numbers. This indicates that the observed dissimilarity between
$\mathcal {E}^m$
and
$\mathcal {E}_\theta ^m$
arises primarily from differences in the dissipation rates within AER alone.
This analysis highlights the role that AER and WMR have, in different flow regimes, on momentum and scalar transfer over rough walls, in agreement with observations in the literature (Zhong et al. Reference Zhong, Hutchins and Chung2023). Future studies should expand upon the present findings by examining a broader range of roughness Reynolds numbers and roughness topographies; in this context, the presented energetic perspective can be valuable for exploring asymptotic scalar transfer laws in the fully rough regime.
Funding
This work was supported by the German Research Foundation (DFG) under the Collaborative Research Centre TRR150 (project number 237267381). The simulations were performed on the HPE Apollo (Hawk) supercomputer at the High Performance Computing Center Stuttgart (HLRS) under grant number 44198.
Declaration of interests
The authors report no conflict of interest.
Data availability statement
The data that support the findings of this study are openly available in the KITopen repository at doi:10.35097/s3jx9n5yhzz2f24g.
Author contributions
Conceptualization, F.S.; methodology, F.S.; software, F.S.; validation, F.S.; formal analysis, F.S., D.G., U.P. and B.F.; investigation, F.S.; resources, B.F.; data curation, F.S.; writing – original draft preparation, F.S.; writing – review and editing, F.S., D.G., U.P. and B.F.; visualization, F.S.; supervision, D.G., U.P. and B.F.; project administration, B.F.; funding acquisition, B.F.
Open access
