1. Introduction
Flows over curved surfaces, involving unsteady separation and reattachment in space and time, occur in numerous engineering applications, such as engine nacelles, curved ducts, bluff bodies and so on. For such complex flows, it is quite challenging to accurately predict the separation and reattachment behaviour at affordable computational cost (Larsson et al. Reference Larsson, Kawai, Bodart and Bermejo-Moreno2016; Bose & Park Reference Bose and Park2018). In particular, the separation location fluctuates spatially and temporally, and strongly affects the downstream flow structures including the subsequent reattachment process.
Flow over a periodic arrangement of smoothly contoured two-dimensional hills (ERCOFTAC test case 81) has in recent years become one of the most widely used test cases to investigate the physics of turbulent separated flows over curved surfaces, as well for validation of computational fluid dynamics codes and turbulence models (Rodi, Bonnin & Buchal Reference Rodi, Bonnin and Buchal1995). The configuration of periodic channel flow was originally proposed by Almeida, Durao & Heitor (Reference Almeida, Durao and Heitor1993), and then modified by Mellen, Fröhlich & Rodi (Reference Mellen, Fröhlich, Rodi, Deville and Owens2000) to be more suitable for numerical simulations. The history of this case is found in the review paper by Breuer et al. (Reference Breuer, Peller, Rapp and Manhart2009). The periodic hills geometry has been investigated both numerically and experimentally over a wide range of Reynolds numbers. Mellen et al. (Reference Mellen, Fröhlich, Rodi, Deville and Owens2000) investigated the periodic hill flow at $Re_h = 7100$ (based on hill height
$h$ and bulk mean velocity through the hill crest
$U_b$) using wall-modelled and wall-resolved large-eddy simulations (WMLES and WRLES), and the impact of different subgrid-scale (SGS) models and grid qualities was assessed. A similar investigation was carried out by Temmerman et al. (Reference Temmerman, Leschziner, Mellen and Fröhlich2003) at
$Re_h = 10\,595$ using WMLES, which combined different SGS models and wall-law functions. It was demonstrated that the flow features are surprisingly more sensitive to the wall model than the SGS model, and the classical wall models developed for attached flows are not satisfactory for this flow. Fröhlich et al. (Reference Fröhlich, Mellen, Rodi, Temmerman and Leschziner2005) performed highly resolved (only the bottom wall) incompressible LES at
$Re_h = 10\,595$ using two different second-order finite-volume discretizations with two different SGS models, i.e. the dynamic Smagorinsky model (Germano et al. Reference Germano, Piomelli, Moin and Cabot1991) and the ‘wall-adapted local eddy-viscosity’ model (Ducros, Nicoud & Poinsot Reference Ducros, Nicoud, Poinsot and Baines1998). A detailed analysis of structural characteristics of this flow configuration has revealed a number of interesting features. Breuer et al. (Reference Breuer, Peller, Rapp and Manhart2009) presented a comprehensive review of direct numerical simulations (DNS) and WRLES performed to date, comparing experimental data with numerical results using DNS up to
$Re_h = 5600$ and WRLES up to
$Re_h = 10\,595$. Rapp & Manhart (Reference Rapp and Manhart2011) experimentally investigated this flow at four Reynolds numbers
$(5600 \leqslant Re_h \leqslant 37\,000)$ with two-dimensional particle image velocimetry and one-dimensional laser Doppler anemometry measurements. The streamwise periodicity of the flow and sidewall effects were evaluated, together with investigations of Reynolds number effects. Kähler, Scharnowski & Cierpka (Reference Kähler, Scharnowski and Cierpka2016) repeated the experiment
$(Re_h=8000,33\,000)$ with high-resolution particle image velocimetry and particle tracking velocimetry techniques. This high-resolution measurement makes possible a precise analysis of the near-wall flow features. Recently, Krank, Kronbichler & Wall (Reference Krank, Kronbichler and Wall2018) presented DNS at
$Re_h=5600$ and
$10\,595$ using spectral incompressible discontinuous Galerkin schemes. Although many DNS and LES have been performed for
$Re_h \leqslant 10\,595$, there has been somewhat limited research for larger Reynolds numbers due to extremely high-resolution requirements in the near-wall region. To the best of our knowledge, only two hybrid Reynolds-averaged Navier–Stokes (RANS)/LES have been performed at
$Re_h=37\,000$ (Chaouat & Schiestel Reference Chaouat and Schiestel2013; Mokhtarpoor, Heinz & Stoellinger Reference Mokhtarpoor, Heinz and Stoellinger2016). For wall-bounded turbulent flows, a tenable solution for investigating higher-Reynolds-number cases is to employ the wall modelling approach since the mesh resolution requirement of WMLES scales linearly with increasing
$Re$ (Choi & Moin Reference Choi and Moin2012).
In the past four decades, several wall models have been proposed for canonical flows in simple geometries (Schumann Reference Schumann1975; Grötzbach Reference Grötzbach1987; Piomelli et al. Reference Piomelli, Ferziger, Moin and Kim1989; Marusic, Kunkel & Porté-Agel Reference Marusic, Kunkel and Porté-Agel2001; Piomelli & Balaras Reference Piomelli and Balaras2002). The reader is referred to the review paper by Bose & Park (Reference Bose and Park2018) for recent developments of wall model techniques. However, there are a couple of primary challenges when it comes to flow over curved surfaces. First, most wall models follow the equilibrium stress assumption and imply a logarithmic-law profile in the near-wall region, which breaks down when turbulent boundary layers are subjected to strong adverse pressure gradients leading to separation, extra strain due to curvature, etc. (Diurno Reference Diurno2001; Bose & Park Reference Bose and Park2018). Meanwhile, some enhanced wall-function models, which take into account the streamwise pressure gradient, have been developed to simulate flow over periodic hills (Breuer, Kniazev & Abel Reference Breuer, Kniazev and Abel2007; Manhart, Peller & Brun Reference Manhart, Peller and Brun2008; Duprat et al. Reference Duprat, Balarac, Métais, Congedo and Brugière2011), but the free parameters in these models restrict their application in general. Second, most wall modelling strategies including detached eddy simulations (DES) fall into the hybrid RANS/LES methodology in complex geometries, which solves the RANS equations in the inner layer and provides wall shear stress boundary conditions for the outer LES region (Cabot & Moin Reference Cabot and Moin1999; Piomelli & Balaras Reference Piomelli and Balaras2002; Kawai & Asada Reference Kawai and Asada2013; Park & Moin Reference Park and Moin2016). This hybrid method is not only sensitive to the choice of the RANS model and its associated model coefficients, but also causes the so-called ‘scale disparity’ problem on the nominal interfaces between the RANS and LES regions (Germano Reference Germano2004; Piomelli Reference Piomelli2008). Alternatively, Chung & Pullin (Reference Chung and Pullin2009) proposed the virtual wall model, which dynamically couples the outer resolved region with the inner wall region, and offers a slip velocity boundary condition for the filtered velocity field on the ‘virtual’ wall. This wall model has been successfully deployed in canonical flows without separation (Inoue & Pullin Reference Inoue and Pullin2011; Saito, Pullin & Inoue Reference Saito, Pullin and Inoue2012), and then extended by Cheng, Pullin & Samtaney (Reference Cheng, Pullin and Samtaney2015) to simulate flat-plate turbulent boundary layer flows with separation and reattachment. Recently, this virtual wall model was extended to generalized curvilinear coordinates by Gao et al. (Reference Gao, Zhang, Cheng and Samtaney2019) and utilized in WMLES for flow past airfoils. The same framework is adopted and tested in the present simulations.
Turbulent boundary layer flow over a curved surface, involving separation and reattachment, generates larger pressure fluctuations than that in the equilibrium boundary layer. Wall-pressure fluctuations play a key role in a variety of engineering applications, such as flow-induced panel flutter and structural vibration, aircraft cabin noise and hydroacoustics of underwater vehicles (Blake Reference Blake1970). Many investigations of wall-pressure fluctuations beneath a turbulent boundary layer have been performed in the past several decades, including zero pressure gradient turbulent boundary layer (Bradshaw Reference Bradshaw1967; Willmarth Reference Willmarth1975; Farabee & Casarella Reference Farabee and Casarella1991); flat-plate turbulent boundary layers with adverse pressure gradient (Mabey Reference Mabey1972; Simpson, Ghodbane & McGrath Reference Simpson, Ghodbane and McGrath1987; Na & Moin Reference Na and Moin1998a,Reference Na and Moinb; Abe Reference Abe2017); and turbulent flows over a backward- or forward-facing step (Farabee & Casarella Reference Farabee and Casarella1986; Ji & Wang Reference Ji and Wang2012; Awasthi et al. Reference Awasthi, Devenport, Glegg and Forest2014; Doolan & Moreau Reference Doolan and Moreau2016). Presently, for turbulent flow in a channel with streamwise periodic constrictions, we analyse the pressure fluctuations by relating these to the mean pressure in the separation bubble and the development of the mixing layer.
In the present investigation, we emphasize three main objectives. First, the extended virtual wall model developed by Gao et al. (Reference Gao, Zhang, Cheng and Samtaney2019) is applied in a periodic channel flow (both bottom and top walls). All the WMLES results are validated with experimental data wherever available. Some WRLES results are also utilized for verifications, especially the pressure and skin friction coefficients which are not reported in the experiments but are important in separation and reattachment. Based on these verifications and validations, the effects of Reynolds number on the skin friction coefficient, separation bubble size and pressure fluctuations are analysed. Second, the details of unsteady separation in this flow have not been reported in the past. Recent work by Cheng, Pullin & Samtaney (Reference Cheng, Pullin and Samtaney2017, Reference Cheng, Pullin and Samtaney2018) in WRLES of flow past a smooth and grooved cylinder at subcritical and supercritical Reynolds numbers emphasized the role of unsteady separation and the dynamics of separation bubbles in the phenomenon of the drag crisis. They attributed the drag crisis, mainly due to a large change in form drag, to the topology change induced by the movement of the location of curvature-controlled large-scale separation. In the present study, the periodic hill channel may be considered as a combination of flat and curved surfaces, and this geometric complexity would result in rich flow physics associated with separation and reattachment. The last main objective comes from the fact that almost all of the previous investigations paid little attention to the flat top wall of the channel. The empirical friction law and universal logarithmic law, which are well captured in plane channel flows, are also evaluated from the top wall of the channel.
The paper is organized as follows. In § 2, we describe the physical set-up, followed by a discussion of the governing equations, and briefly discuss wall and SGS models employed (details are relegated to appendices). In § 3, the WMLES results are validated and verified with experimental and WRLES results. After that, in §§ 4 and 5, we provide new insights based on LES results up to $Re_h=10^5$. The flow near the bottom wall is analysed in § 4, emphasizing on the separation/reattachment behaviour. Here we investigate scaling relations for the skin friction coefficient and velocity profiles in separation zones. Instantaneous skin friction lines at three
$Re_h$ are compared, with focus on the existence of a small separation bubble near the top of the hill. Further in § 5, the flow at the top wall is characterized using the empirical friction law and universal logarithmic law which are proposed for planar channel flow. Finally, the conclusions are drawn in § 6.
2. Physical and numerical set-up, equations and models
In this section, the flow configuration is described, and following that we present the essential set of equations including the boundary conditions at the virtual wall required to perform WMLES in a wall-bounded domain. The detailed derivations for the wall model in the generalized curvilinear coordinates were given by Gao et al. (Reference Gao, Zhang, Cheng and Samtaney2019). We include the details of the wall model, the SGS model and the numerical methods in appendices for the sake of completeness.
2.1. Flow configuration
We perform LES of turbulent flow past periodically constructed hills on the bottom wall in a channel with a flat wall on the top. The physical geometry and simulation domain are illustrated in figure 1. The shape of the hill is defined in the form of a polynomial, taken from the experimental study by Almeida et al. (Reference Almeida, Durao and Heitor1993). This is also a standard test case for an ERCOFTAC/IAHR workshop, and has been widely used to validate various numerical schemes and physical models. It should be noted that geometry G1 (figure 1a) was adopted in WRLES by Fröhlich et al. (Reference Fröhlich, Mellen, Rodi, Temmerman and Leschziner2005) for $Re_h=10\,595$ and geometry G2 (figure 1b) was the region of focus in the experimental investigations of Kähler et al. (Reference Kähler, Scharnowski and Cierpka2016). These two geometries are essentially identical since the hills are constructed periodically in the
$x$ direction. For the sake of simplicity of verifications and validations, G1 is adopted in the case of
$Re_h=10\,595$, while G2 is used for the cases of
$Re_h=33\,000$ and
$Re_h=10^5$.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_fig1.png?pub-status=live)
Figure 1. Sketch of the hill geometry adopted in the present paper. (a) Geometry G1: simulation geometry adopted in the present WMLES and WRLES of Fröhlich et al. (Reference Fröhlich, Mellen, Rodi, Temmerman and Leschziner2005) for $Re_h=10\,595$; (b) geometry G2: simulation geometry adopted in the present WMLES and WRLES and experimental domain of Kähler et al. (Reference Kähler, Scharnowski and Cierpka2016) for
$Re_h=33\,000$.
2.2. Governing equations
To compute the flow in this set-up, assuming constant (unity) density, we solve the filtered incompressible Navier–Stokes equations on a body-fitted curvilinear grid. In generalized curvilinear coordinates, the equations of motion are
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn1.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn2.png?pub-status=live)
where $(\xi ^1,\xi ^2,\xi ^3)=(\xi ,\eta ,\zeta )$ denote the generalized curvilinear coordinates;
$U^m$ (the volume flux normal to the surface of constant
$\xi ^m$) and
$F_i^m$ are given by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn3.png?pub-status=live)
where $(x_1,x_2,x_3)=(x,y,z)$ are the Cartesian coordinates with corresponding velocity components
$(u_1,u_2,u_3)=(u,v,w)$,
$p$ is the pressure,
$\nu$ is the kinematic viscosity,
$J^{-1}$ is the Jacobian of the transformation and
$ {G}^{mn}$ is the mesh skewness tensor. It should be noted that the
$\zeta$ direction is congruent with the
$z$ direction since the spanwise geometry is uniform in the present research. Applying a nominal filter to the incompressible Navier–Stokes equations, the filtered LES equations are written below in terms of the resolved velocity field:
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn4.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn5.png?pub-status=live)
where tildes denote filtered quantities and $ {T}_{ij} = \widetilde {u_iu_j} - \tilde u_i \tilde u_j$ is the SGS stress tensor. The stretched spiral vortex model is adopted to compute
$ {T}_{ij}$, details of which are in appendix A.
2.3. Boundary conditions
For the simulations, periodic boundary conditions are prescribed in the spanwise and streamwise directions. The spanwise extent ($L_z$) is of importance in order to obtain reliable and physically reasonable results. To ensure this criterion, the two-point correlations in the spanwise direction have to decay to sufficiently small values in the half-width of the domain size chosen. Based on the investigations by Mellen et al. (Reference Mellen, Fröhlich, Rodi, Deville and Owens2000), a spanwise extension of the computational domain of
$L_z=4.5h$ is used in all computations presented. This spanwise domain extent was also used in the investigation by Fröhlich et al. (Reference Fröhlich, Mellen, Rodi, Temmerman and Leschziner2005) and many other DNS/LES studies (Breuer et al. Reference Breuer, Peller, Rapp and Manhart2009; Krank et al. Reference Krank, Kronbichler and Wall2018). It represents a well-balanced compromise between spanwise extent and spanwise resolution.
For boundary conditions on the solid walls (both bottom and top walls), the no-slip boundary condition is specified at the actual wall in WRLES, and a Dirichlet boundary condition for the velocity, derived from the wall model, is specified at the virtual wall in WMLES. Similar to the turbulent plane channel flow case, the non-periodic behaviour of the pressure field is accounted for by inclusion of a mean pressure gradient as a source term in the streamwise momentum equation. Two alternatives are possible: either the pressure gradient is fixed which results in a fluctuating mass flux or the mass flux is held constant which requires adjustment of the mean pressure gradient in time. Since a fixed Reynolds number can only be guaranteed by a fixed mass flux, the second option is chosen using the method proposed by Benocci & Pinelli (Reference Benocci, Pinelli, Rodi and Ganic1990).
2.4. Wall model
The virtual wall model in generalized curvilinear coordinates developed by Gao et al. (Reference Gao, Zhang, Cheng and Samtaney2019), coupled with the stretched vortex SGS model, has been strongly verified and validated in various airfoil flows. We briefly describe the essential idea of the wall model with details relegated to appendix B. It should be noted that, consistent with other approaches involving body-fitted mesh computations, the computational mesh is constrained to be locally orthogonal to the solid walls for wall-normal averaging, and the wall-normal coordinate is denoted as $y_n$. As shown in figure 2, the distance
$h'$ is typically chosen as the distance of the first grid point from the solid wall and
$h_0$ is the height of the virtual wall that is further discussed below.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_fig2.png?pub-status=live)
Figure 2. Sketch of the near-wall velocity components. The dashed line above the solid wall denotes the virtual wall, and the blue point denotes the centre of the first grid cell off the solid wall.
We define the magnitude of the resultant velocity $\tilde q$ and velocity angle
$\theta$ on the wall-parallel plane as
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn6.png?pub-status=live)
where $\tilde {u}_s$ is the streamwise velocity parallel to the solid wall (as shown in figure 2), given by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn7.png?pub-status=live)
where $\theta _w$ denotes the local angle between the solid wall and the
$x$ coordinate.
The wall-normal gradient of $\tilde q$, i.e.
$\eta _0 (\xi , \zeta , t)$, and the associated friction velocity
$u_\tau$ are defined as
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn8.png?pub-status=live)
Based on the near-wall filtering and wall-normal averaging approach, the following ordinary differential equation (ODE) for $\eta _0$ can then be derived (see appendix B for details):
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn9.png?pub-status=live)
where
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn10.png?pub-status=live)
Detailed expressions for $F_\xi$,
$F_\zeta$ and
$M$ are given by (B 10), (B 12) and (B 14), respectively; an approximate analytical solution to (2.9) is given in appendix B.
Once $\eta _0(\xi ,\zeta ,t)$ is known, the velocity angle
$\theta (\xi , z, t)$ is estimated as
$\arccos (\tilde u_s|_{h'}/\tilde q|_{h'})$ from the first grid cell of the resolved LES field, an approximation justified based on the work of Cheng et al. (Reference Cheng, Pullin and Samtaney2015) (turbulent boundary layer flow with separation and reattachment) and Gao et al. (Reference Gao, Zhang, Cheng and Samtaney2019) (separated flow past airfoils). The local wall shear stress components may then be computed as
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn11.png?pub-status=live)
Here $\mu =\rho \nu$ is the dynamic viscosity and
${\boldsymbol {\tau }}_w \equiv (\tau _{w,s}, \tau _{w,z})$ is the LES representation of the surface stress vector. Above, we make the approximation that the velocity angle
$\theta$ is constant within the first grid cell,
$0 \leqslant y_n \leqslant h$. Cheng et al. (Reference Cheng, Pullin and Samtaney2015) proposed an algebraic model for
$\theta$ in turbulent boundary layer simulations and concluded that there is little difference between the constant velocity angle model and the algebraic model. In the present paper, the constant velocity angle model is adopted for simplicity.
Finally, we present a slip velocity $\tilde q|_{h_0}$ on the virtual wall as
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn12.png?pub-status=live)
where $h_0^+=u_{\tau } h_0/\nu$ and
$h^+_\nu$ is the intercept between the linear and logarithmic components in the law of the wall. Experimental research shows that the outer edge of the viscous sublayer is located at
$h^+_\nu \approx 11$, which is approximately equivalent to the offset (
$=5.0$) in the classical logarithmic law. This empirical value is adopted by Chung & Pullin (Reference Chung and Pullin2009) and Inoue & Pullin (Reference Inoue and Pullin2011), and also by Cheng et al. (Reference Cheng, Pullin and Samtaney2015) and Gao et al. (Reference Gao, Zhang, Cheng and Samtaney2019) in modelling the boundary layer flows with separation and reattachment. In the present case,
$h^+_\nu = 11$ is used as an empirical parameter in the wall model. Both
$h_0^+$ and
$h^+_\nu$ are fixed in all the simulations presented, the same values being used in our previous work, and hence these parameters may not be considered as ‘tunable’ parameters. In the attached region (
$\tau _{w,s}>0$), the linear–logarithmic relation is essentially the same as that of Chung & Pullin (Reference Chung and Pullin2009), which is derived from the stretched-vortex SGS model (see appendix A) and the Kármán-like constant
$\mathscr {K}_1$ is dynamically computed. In the separated region (
$\tau _{w,s} \leqslant 0$), the log-like relation is no longer valid and Cheng et al. (Reference Cheng, Pullin and Samtaney2015) proposed a linear relationship which appears to work reasonably well in regions of flow separation. Here, we follow the linear law of Cheng et al. (Reference Cheng, Pullin and Samtaney2015).
Based on the above formulation, the wall model can be summarized as follows. In the near-wall region, (2.9) is solved for $\eta _0$, in which the coefficients on the right-hand side are approximated with the resolved LES field at the first grid cell, i.e.
$h'=h_0+{\rm \Delta} y_n/2$ (the choice of
$h_0$ is given in appendix B). Equation (2.12) is then used to compute the resultant velocity
$\tilde q|_{h_0}$ on the virtual wall with the streamwise and spanwise velocity components given by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn13.png?pub-status=live)
The contribution of the wall-normal velocity component $\tilde u_n|_{h_0}$ to
$\tilde u$ and
$\tilde v$ is assumed to be small comparing with
$\tilde u_s|_{h_0}$, and we use
$\tilde u_n|_{h_0}=0$. Finally, the slip velocity boundary condition on the virtual wall
$y_n=h_0$ is
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn14.png?pub-status=live)
with the spanwise velocity component $\tilde w|_{h_0}$ given by (2.13b).
2.5. Summary of numerical cases
Three cases, as summarized in table 1, are considered. The energy-conservative fourth-order finite-difference scheme is used for spatial discretizations, and the discretized governing equations are solved using a semi-implicit fractional step method (see appendixC for details). For $Re_h=10\,595$, only WMLES is performed and compared with both WRLES from Fröhlich et al. (Reference Fröhlich, Mellen, Rodi, Temmerman and Leschziner2005) and experimental data from Rapp & Manhart (Reference Rapp and Manhart2011). To check the effects of mesh resolution on the WMLES, a mesh convergence study of this case is presented in appendix D. For
$Re_h=33\,000$, neither the pressure nor skin friction coefficients were measured in the experimental investigation of Kähler et al. (Reference Kähler, Scharnowski and Cierpka2016); thus we perform both WRLES and WMLES. Accurate predictions of the separation and reattachment depend on the skin friction coefficient, which is also challenging in terms of wall modelling especially for flow over complex geometries. The WRLES results are utilized to verify the WMLES results for some quantities that were not reported in the experiment for
$Re_h=33\,000$. Furthermore, the WRLES (R2) case not only resolves the flow at the lower wall but also resolves the upper wall by a DNS-like representation. We briefly remark that the computational expense of WMLES is approximately 12 times that of WMLES per time interval with the same CPU cores. The grid spacings in the wall-normal direction (i.e.
$\eta$ direction) on both the bottom and top walls are illustrated in figure 3, and the mesh size ratios are also given in table 1. Note that we employ a higher resolution in the spanwise direction than the streamwise direction for case M1. For cases M2 and M3 the resolution in terms of (
${\rm \Delta} \xi ^+/ {\rm \Delta} \eta ^+_b$,
${\rm \Delta} z^+/ {\rm \Delta} \eta ^+_b$) is (5.1, 7.0) for M2 and (4.1, 7.0) for M3, respectively. We believe that the disparity in resolution between the spanwise and streamwise directions is not severe within the context of WMLES.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_fig3.png?pub-status=live)
Figure 3. Dimensionless wall-normal grid spacings in wall units ${\rm \Delta} \eta ^+_{(b,t)}$: (a) bottom wall and (b) top wall.
, WMLES for
$Re_h=10\,595$;
, WRLES for
$Re_h=33\,000$;
, WMLES for
$Re_h=33\,000$;
, WMLES for
$Re_h=10^5$. Note that the right-hand
$y$ axis is for WRLES.
Table 1. Summary of the performed numerical cases. The ‘$+$’ superscript indicates the expression of the mesh quantities using wall units and the subscripts ‘
$b$’ and ‘
$t$’ denote the bottom and top walls, respectively. Time
$t_r=9h/U_b$ is the typical flow-through time and
$t_{a}$ is the total simulation time.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_tab1.png?pub-status=live)
3. Numerical results for statistically averaged quantities
In this section, we present several time- and spanwise-averaged quantities that are used to verify and validate the present WMLES/WRLES code: these are the skin friction coefficient, $C_f$; the pressure coefficient,
$C_p$; the normalized velocity profiles in
$x$ and
$y$ directions,
$\bar u/U_b$ and
$\bar v/U_b$; and the normalized Reynolds stresses profiles,
$\overline {u^{\prime }u^{\prime }}/U_b^2$,
$\overline {v^{\prime }v^{\prime }}/U_b^2$ and
$\overline {u^{\prime }v^{\prime }}/U_b^2$. The Reynolds stress containing SGS corrections to the resolved flow is calculated as
$\overline {u'_iu_j'} = \overline {\tilde {u}'_i\tilde {u}_j'} + \overline{{T}}_{ij}$ (Inoue & Pullin Reference Inoue and Pullin2011). The skin friction coefficient and pressure coefficient are computed as
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn15.png?pub-status=live)
where $\bar {\tau }_{w,s}$ is the local streamwise (parallel to the actual wall) mean wall shear stress,
$\bar {p}$ is the local pressure and
$p_{r}$ is a reference wall pressure taken at the centre of the hill crest. In WMLES,
$\bar {\tau }_{w,s}$ is computed from the wall model; in WRLES, the third-order accuracy one-sided velocity-derivative method (Cheng et al. Reference Cheng, Pullin and Samtaney2017) is used for the wall-normal differentiation. It should be noted that all the flow configurations are rescaled to geometryG1 (see figure 1a) for simplicity of investigation, and this compatible coordinate system will be used in the later sections if not mentioned specifically.
3.1. Reynolds number
$Re_h=10\,595$
The skin friction coefficient $C_f$ and pressure coefficient
$C_p$ for
$Re_h=10\,595$ on both the bottom and top walls are shown in figure 4. We compare our results with the highly resolved LES results from Fröhlich et al. (Reference Fröhlich, Mellen, Rodi, Temmerman and Leschziner2005) for verification (note that Fröhlich etal. (Reference Fröhlich, Mellen, Rodi, Temmerman and Leschziner2005) did not report
$C_f$ at the top wall). The present WMLES shows good agreement with the reference LES data, although slightly smaller
$C_f$ on the leeward side of the hill within the separation zone, i.e.
$x/h \approx 1.0$–
$2.0$. The
$C_f$ plot shows the primary separation point at approximately
$x/h=0.22$, with reattachment at
$x/h=4.59$. This compares well with the experimental results of Rapp & Manhart (Reference Rapp and Manhart2011), who reported the reattachment point at
$x/h=4.21$. The predicted values also compare well with the WRLES (Fröhlich et al. Reference Fröhlich, Mellen, Rodi, Temmerman and Leschziner2005) values of
$x/h=0.20$ and
$x/h=4.56$ (see figure 4b), respectively, for the separation and reattachment points. Two very small recirculation regions, one on the hill top,
$x/h \approx 0$, and the other in the post-reattachment zone on the windward side of hill,
$x/h \approx 7.0\text {--}7.4$, for
$Re_h > 200$ are reported by Breuer et al. (Reference Breuer, Peller, Rapp and Manhart2009), and also confirmed for
$Re_h =10\,595$ by WRLES (Mellen et al. Reference Mellen, Fröhlich, Rodi, Deville and Owens2000; Fröhlich et al. Reference Fröhlich, Mellen, Rodi, Temmerman and Leschziner2005) and DNS (Diosady & Murman Reference Diosady and Murman2014; Krank et al. Reference Krank, Kronbichler and Wall2018). However, the experimental results at
$Re_h=8000$ and
$33\,000$ from Kähler et al. (Reference Kähler, Scharnowski and Cierpka2016) do not confirm the existence of the local separation bubble in these regions, and this was attributed to an insufficient duration of time for averaging. These very small separation regions are also absent in the present WMLES, and are further discussed in § 4.5.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_fig4.png?pub-status=live)
Figure 4. (a) Pressure coefficient $C_p$ and (b) skin friction coefficient
$C_f$ along the bottom and top walls for
$Re_h=10\,595$. The solid and dash-dotted curves are the corresponding results at the bottom and top walls from the present WMLES. In (a), the circles and squares denote
$C_p$ on the bottom and top walls from WRLES by Fröhlich et al. (Reference Fröhlich, Mellen, Rodi, Temmerman and Leschziner2005); in (b), the circles denote
$C_f$ on the bottom wall from WRLES by Fröhlich et al. (Reference Fröhlich, Mellen, Rodi, Temmerman and Leschziner2005). The bold grey line denotes the geometry of the channel with streamwise periodic hills. A horizontal dashed line
$C_f=0$ in (b) is shown for reference: zero crossing of this line indicates separation/reattachment.
Figure 5 shows the normalized mean velocity profiles ($\bar u/U_b$ and
$\bar v/U_b$) in
$x$ and
$y$ directions along vertical lines at ten different locations, i.e.
$x/h=0.05, 0.5,1,2,3,4,5,6,7,8$. Comparisons with experimental data from Rapp & Manhart (Reference Rapp and Manhart2011) show excellent agreement. For profiles in the separation zone, i.e.
$x/h=0.5,1,2,3,4$, negative back flow is noted near the bottom wall (see figure 5a), and the peak value of the back-flow velocity occurs at
$x/h=2$, close to the centre of the main separation zone. The steep velocity gradient near the top wall is also captured by the present wall model, which is absent in WRLES by Fröhlich et al. (Reference Fröhlich, Mellen, Rodi, Temmerman and Leschziner2005) due to low mesh resolution near the top wall in their simulation. Since the total mass flux through the channel was fixed, the poor prediction of
$\bar u/U_b$ near the top wall has a certain influence on the overall flow development (Breuer et al. Reference Breuer, Peller, Rapp and Manhart2009). The peak value of
$\bar v/U_b$ is observed at
$x/h=8$ (see figure 5b) due to strong acceleration on the windward side of the hill (see also figure 4a for the pressure distributions).
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_fig5.png?pub-status=live)
Figure 5. Normalized mean velocity profiles in the $x$ direction (a) and
$y$ direction (b) for
$Re_h=10\,595$. From left to right, the mean velocity profiles are located at
$x/h=0.05,0.5,1,2,3,4,5,6,7,8$. For the sake of clarity these profiles are shifted by
$2$ and
$0.5$ for
$\bar u/U_b$ and
$\bar v/U_b$, respectively.
$\circ$, Experiment by Rapp & Manhart (Reference Rapp and Manhart2011);
, present WMLES.
Figure 6 shows the normalized Reynolds stresses profiles ($\overline {u^{\prime }u^{\prime }}/U_b^2$,
$\overline {v^{\prime }v^{\prime }}/U_b^2$ and
$\overline {u^{\prime }v^{\prime }}/U_b^2$) at the same locations where profiles of mean velocity were presented above. Both
$\overline {v^{\prime }v^{\prime }}/U_b^2$ and
$\overline {u^{\prime }v^{\prime }}/U_b^2$ agree quite well with the experimental and WRLES results. However, deviations are visible for
$\overline {u^{\prime }u^{\prime }}/U_b^2$ in the range
$x/h=4$–
$6$. Considering only the post-reattachment region at
$x/h=5$, the measured peak values of
$\overline {u^{\prime }u^{\prime }}/U_b^2$ are at most
$8\,\%$ higher than in the present WMLES. Nevertheless, the locations of the peak values compare well with both experimental and WRLES results. In the separation zone,
$x/h=0.5$–
$3$, these comparisons are very satisfactory. In the vicinity of the top wall, the steep variations of
$\overline {u^{\prime }u^{\prime }}/U_b^2$ are captured by the present WMLES but absent in WRLES of Fröhlich et al. (Reference Fröhlich, Mellen, Rodi, Temmerman and Leschziner2005). The maximum
$\overline {v^{\prime }v^{\prime }}/U_b^2$ and
$\overline {u^{\prime }v^{\prime }}/U_b^2$ occurs in the separation zone, inside the constant
$C_p$ region (see figure 4a), in which flow deceleration and acceleration occur alternately.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_fig6.png?pub-status=live)
Figure 6. Normalized Reynolds stress profiles: (a) $\overline {u^{\prime }u^{\prime }}/U_b^2$, (b)
$\overline {v^{\prime }v^{\prime }}/U_b^2$ and (c)
$-\overline {u^{\prime }v^{\prime }}/U_b^2$ for
$Re_h=10\,595$. From left to right, the Reynolds stress profiles are located at
$x/h=0.05,0.5,1,2,3,4,5,6,7,8$, shifted by
$0.1$ along the abscissa.
$\circ$, Experiment by Rapp & Manhart (Reference Rapp and Manhart2011);
, present WMLES;
, WRLES from Fröhlich et al. (Reference Fröhlich, Mellen, Rodi, Temmerman and Leschziner2005).
3.2. Reynolds number
$Re_h=33\,000$
The skin friction coefficient $C_f$ and pressure coefficient
$C_p$ for
$Re_h=33\,000$ on both the bottom and top walls are shown in figure 7. The present WRLES (case R2) is used to verify the performance of the wall model. Similar to that for
$Re_h=10\,595$, a nearly constant pressure plateau is observed in the
$C_p$ distribution from the separation point to the centre of the separation zone (
$x/h \approx 2$). The skin friction coefficient
$C_f$ in the present WMLES is slightly smaller than that from WRLES in the separation zone (
$x/h \approx 1.5\text {--}2.5$). Overall, the present WMLES results agree well with the WRLES results. The separation and reattachment points from the present WMLES
$C_f$ plot occur at approximately
$x/h=0.27$ and
$x/h=3.94$, respectively. These predictions are comparable to the corresponding values of
$x/h=0.27$ and
$x/h=3.96$ from the present WRLES (see figure 7b). These also match well with the experimental results of Kähler et al. (Reference Kähler, Scharnowski and Cierpka2016), who reported the separation point at
$x/h=0.34 \pm 0.05$ and the reattachment point at
$x/h=3.80\pm 0.05$. Compared with the
$Re_h=10\,595$ case, the peak value of
$C_f$ in both attached and separated zones decreases, confirming the observations by Breuer et al. (Reference Breuer, Peller, Rapp and Manhart2009). Furthermore, the undulations in
$C_f$ at the beginning of the main separation bubble seem to be related to geometry and not much affected by Reynolds number effects. The very small separation zone at the hill top is detected in the region from
$x/h=8.932$ to
$9.985$ in the present WRLES, but absent in the WMLES, as was also the case for
$Re_h=10\,595$. The very small separation bubble at the foot of thewindward side of the hill does not occur in either simulation from an inspection of the
$C_f$ plot.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_fig7.png?pub-status=live)
Figure 7. (a) Pressure coefficient $C_p$ and (b) skin friction coefficient
$C_f$ along the bottom and top walls for
$Re_h=33\,000$.
$\circ$,
$C_p$ and
$C_f$ on the bottom wall from the present WRLES;
$\square$,
$C_p$ and
$C_f$ on the top wall from the present WRLES;
,
$C_p$ and
$C_f$ on the bottom wall from the present WMLES;
,
$C_p$ and
$C_f$ on the top wall from the present WMLES. The bold grey line denotes the geometry of the channel with streamwise periodic hills. A horizontal dashed line
$C_f=0$ in (b) is shown for reference: zero crossing of this line indicates separation/reattachment.
Figure 8 shows the normalized mean velocity profiles in $x$ and
$y$ directions along vertical lines at
$17$ different locations, i.e.
$x/h \in [0,9]$, with a constant gap distance of
$0.5$. Comparisons with experimental data from Kähler et al. (Reference Kähler, Scharnowski and Cierpka2016) show good agreement, although the
$\bar v/U_b$ profile at
$y/h \approx 1.0$ (in the separated shear layer) is slightly smaller than the measured values after the separation point (see figure 8b). Overall, the shape of the mean velocity profiles is similar to that for
$Re_h=10\,595$ at the same locations, which may be construed as agreement with the
$Re$-independent behaviour suggested by Kähler et al. (Reference Kähler, Scharnowski and Cierpka2016). For profiles in the separation zone, i.e.
$x/h=0.5\text {--}3.5$, a negative back-flow velocity is clearly observed near the bottom wall (see figure 8a), and the peak value of the back-flow velocity occurs at
$x/h=2.0$, close to the centre of the main separation zone. The steep velocity gradient near the top wall is also captured by the present wall model (see also figure 8a). The peak value of
$\bar v/U_b$ is visible at
$x/h=8.5$ (see figure 8b) due to strong acceleration on the windward side of the hill, similar to the case for
$Re_h=10\,595$. The magnitude of this extremum is larger than that in the
$Re_h=10\,595$ case, and the vertical location is shifted towards the bottom wall. This is in accordance with the trend described by Breuer et al. (Reference Breuer, Peller, Rapp and Manhart2009) and Kähler et al. (Reference Kähler, Scharnowski and Cierpka2016). The downward flow (
$\bar {v}/U_b <0$) enhances the momentum transfer from the upper half of the channel to the separation zone, which is fed by the outer high-momentum fluid, and thus promotes earlier reattachment compared with the
$Re_h=10\,595$ case.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_fig8.png?pub-status=live)
Figure 8. Normalized mean velocity profiles in the $x$ direction (a) and
$y$ direction (b) for
$Re_h=33\,000$. From left to right, the mean velocity profiles are located at
$x/h \in [0,9]$ with gap distance of 0.5, shifted by
$2$ and
$0.5$ for
$\bar u/U_b$ and
$\bar v/U_b$, respectively.
$\circ$, Experiment by Kähler et al. (Reference Kähler, Scharnowski and Cierpka2016);
, present WMLES;
, present WRLES.
Figure 9 shows the normalized Reynolds stress profiles at the same monitoring locations. For flow outside the developing shear layer, i.e. $y/h<0.5$ and
$y/h>1.5$, both WMLES and WRLES results compare well with the experimental data from Kähler etal. (Reference Kähler, Scharnowski and Cierpka2016), but WRLES overpredicts the
$\overline {u^{\prime }u^{\prime }}/U_b^2$ profiles in the vicinity of the bottom wall, especially in the post-separation zone (see figure 9a). The overshoot of the Reynolds stresses after the hill top (
$x/h>0$) is also captured well by the simulations. In the separated shear layer region, both WMLES and WRLES overpredict the three Reynolds stress components, but WRLES matches better with the peak values of these profiles. Bose & Park (Reference Bose and Park2018) argued that the proper resolution of the separated shear layer is critical in WMLES, and Kähler et al. (Reference Kähler, Scharnowski and Cierpka2016) also reported that the unresolved mixing process in the simulations and insufficient averaging time (time for ensemble averaging in the LES and DNS should be 12–48 times longer) could also affect the flow statistics. The contribution of the SGS stress to the total Reynolds stress is shown in appendix E, and we found that the fraction of the SGS stress is larger for
$Re_h=33\,000$ than
$Re_h=10\,595$, especially in the separated shear layer. Nevertheless, the reason behind these discrepancies should be explored further. The peak from the experiment displays a narrower distribution than the present LES. The maximum
$\overline {u^{\prime }u^{\prime }}/U_b^2$ occurs upstream of the primary separation point, similar to
$Re_h=10\,595$, because we expect larger fluctuations in the region where the transition from attached to separated flow states occurs. The maximum
$\overline {v^{\prime }v^{\prime }}/U_b^2$ and
$\overline {u^{\prime }v^{\prime }}/U_b^2$ is also located in the constant
$C_p$ region, similar to the
$Re_h=10\,595$ case. These qualitative observations are also in accordance with the experimental research of Kähler et al. (Reference Kähler, Scharnowski and Cierpka2016).
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_fig9.png?pub-status=live)
Figure 9. Normalized Reynolds stress profiles: (a) $\overline {u^{\prime }u^{\prime }}/U_b^2$, (b)
$\overline {v^{\prime }v^{\prime }}/U_b^2$ and (c)
$-\overline {u^{\prime }v^{\prime }}/U_b^2$ for
$Re_h=33\,000$. From left to right, the Reynolds stress profiles are located at
$x/h \in [0,9]$ with gap distance of 0.5, shifted by
$0.1$ along the abscissa.
$\circ$, Experiment by Kähler et al. (Reference Kähler, Scharnowski and Cierpka2016);
, present WMLES;
, present WRLES.
4. Reynolds number effects at the bottom wall: separation and reattachment
In the previous section, we noted that the our WMLES results are in reasonable agreement with both experimental and WRLES results. Based on these verifications and validations, we consider another case with $Re_h=10^5$ (numerical details in table 1), which hitherto is the highest
$Re_h$ case in periodic hill channel flows. We focus on the mean skin friction coefficients indicative of separation and reattachment, separation bubble size and pressure fluctuations that are related to the mean pressure and the developing mixing layer after the crest of the hill.
4.1. Mean separation and reattachment: distribution of
$C_f$
4.1.1. Skin friction coefficient
$C_f$ and its maximum value
The skin friction coefficient $C_f$ for
$Re_h=10^5$ on both the bottom and top walls is shown in figure 10(a). The separation and reattachment points from the present WMLES
$C_f$ plot occur at approximately
$x/h=0.26$ and
$x/h=3.57$, respectively. Two aspects of the skin friction coefficient are assessed quantitatively:
$(1)$ the peak value of the skin friction coefficient (
$C_{f,max}$) on the bottom wall and
$(2)$ the location of the mean reattachment point (
$x_r/h$).
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_fig10.png?pub-status=live)
Figure 10. (a) Skin friction coefficient $C_f$ along the bottom and top walls for
$Re_h=10^5$. The solid and dash-dotted lines denote the WMLES results at the bottom and top walls. The bold grey line denotes the geometry of the channel with streamwise periodic hills. A horizontal dashed line
$C_f=0$ is shown for reference: zero crossing of this line indicates separation/reattachment. (b) The peak value of
$C_f$ for
$Re_h$ ranging from
$700$ to
$10^5$:
$\circ$, DNS/WRLES from Breuer et al. (Reference Breuer, Peller, Rapp and Manhart2009);
$\times$, WRLES from Fröhlich et al. (Reference Fröhlich, Mellen, Rodi, Temmerman and Leschziner2005);
$\square$, the present WMLES results;
, scaling equation (4.1).
The peak values of $C_f$ for different
$Re_h$
$(700 \leqslant Re_h \leqslant 10^5)$ are summarized in figure 10(b). For
$700 \leqslant Re_h \leqslant 5600$, the data are from DNS results of Breuer et al. (Reference Breuer, Peller, Rapp and Manhart2009); for
$Re_h=10\,595$, the data are from the WRLES results of Breuer et al. (Reference Breuer, Peller, Rapp and Manhart2009) and Fröhlich et al. (Reference Fröhlich, Mellen, Rodi, Temmerman and Leschziner2005) and the present WMLES result. It is found that the power-law fit of the data for
$700 \leqslant Re_h \leqslant 10\,595$ follows
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn16.png?pub-status=live)
The predicted values of $C_{f,max}$ from WMLES at
$Re_h =33\,000$ and
$10^5$ match well with this power-law scaling equation (4.1), as shown in figure 10(b).
4.1.2. Skin friction coefficient
$C_f$ inside the separation zone
In the separation zone, Adams, Johnston & Eaton (Reference Adams, Johnston and Eaton1984) proposed a ‘laminar-like’ skin friction law of the form $C_{f,U_N} \varpropto Re_N^{-1}$. Here
$C_{f,U_N}$ is the local skin friction coefficient normalized by
$\frac {1}{2} \rho U_N^2$, i.e.
$C_{f,U_N}=2|\tau _{w,s}|/\rho U_N^2$, where
$U_N$ is the magnitude of the maximum back-flow velocity. Parameter
$Re_N$ is the wall-layer Reynolds number based on
$U_N$ and
$N$, where
$N$ is the distance of
$U_N$ away from the wall. The values of
$U_N$ and
$N$ can be determined locally by a search in the wall-normal direction. Le, Moin & Kim (Reference Le, Moin and Kim1997) proposed another correlation form based on the DNS results of backward-facing step flow at
$Re_h=5100$:
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn17.png?pub-status=live)
This formula does not quite follow the ‘$-1$ slope’ of Adams et al. (Reference Adams, Johnston and Eaton1984). Numerical data from the present simulations are shown in figure 11, based on which another correlation shows better fit, i.e.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn18.png?pub-status=live)
4.2. Separation bubble size
The locations of separation ($x_s/h$) and reattachment (
$x_r/h$) at different
$Re_h$ are summarized in figure 12. Since numerical methods and physical models (SGS model and wall model) affect the locations of separation and reattachment, results from previous numerical studies are also included in the graph. It is seen that the locations of reattachment move towards the upstream direction with increasing
$Re_h$ due to stronger mixing (Breuer et al. Reference Breuer, Peller, Rapp and Manhart2009; Kähler et al. Reference Kähler, Scharnowski and Cierpka2016). The length of the separation bubble in numerical simulations (DNS and LES) is larger than that in the experiment (see also table 2), and Kähler et al. (Reference Kähler, Scharnowski and Cierpka2016) attributed this to an unresolved turbulent mixing process, i.e. the turbulence level in the mixing layer is lower than in the experiment. Kähler et al. (Reference Kähler, Scharnowski and Cierpka2016) proposed a power scaling law for the location of reattachment based on the experimental data, i.e.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn19.png?pub-status=live)
The estimated reattachment for $Re_h = 10^5$ occurs at
$x_r/h=3.73$, which is 4.3 % larger than the WMLES results (
$x_r/h=3.57$) as shown in figure 10(b).
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_fig12.png?pub-status=live)
Figure 12. Summary of (a) separation and (b) reattachment locations for different Reynolds numbers with various approaches: , WMLES from Mellen et al. (Reference Mellen, Fröhlich, Rodi, Deville and Owens2000);
$\times$, WMLES from Temmerman et al. (Reference Temmerman, Leschziner, Mellen and Fröhlich2003);
$\square$, WRLES from Fröhlich et al. (Reference Fröhlich, Mellen, Rodi, Temmerman and Leschziner2005);
$\blacktriangle$, DNS from Peller & Manhart (Reference Peller and Manhart2006);
$\star$, DES and LES from Šarić et al. (Reference Šarić, Jakirlić, Breuer, Jaffrézic, Deng, Chikhaoui, Fröhlich, Von Terzi, Manhart, Peller, Cancès and Gerbeau2007);
$+$, WRLES from Hickel, Kempe & Adams (Reference Hickel, Kempe and Adams2008);
$\blacktriangledown$, WRLES from Ziefle, Stolz & Kleiser (Reference Ziefle, Stolz and Kleiser2008);
$\diamond$, DNS and WRLES from Breuer et al. (Reference Breuer, Peller, Rapp and Manhart2009);
, WMLES from Mokhtarpoor et al. (Reference Mokhtarpoor, Heinz and Stoellinger2016);
$\ast$, DNS from Krank et al. (Reference Krank, Kronbichler and Wall2018);
$\bullet$, experiment results from Rapp & Manhart (Reference Rapp and Manhart2011);
$\blacksquare$, experiment results from Kähler et al. (Reference Kähler, Scharnowski and Cierpka2016);
$\circ$, WM/WRLES from the present research;
, scaling equation (4.4) given by Kähler et al. (Reference Kähler, Scharnowski and Cierpka2016).
Table 2. Comparison of bubble size with experimental and LES results: centre coordinate of the recirculation zone ($x_c$,
$y_c$), length of the bubble (
$L_b$) and height of the bubble (
$H_b$).
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_tab2.png?pub-status=live)
The quantitative comparisons (bubble centre location, $x_c/h, y_c/h$; height and length of the bubble,
$H_b/h,L_b/h$) that characterize the extent of the separation bubbles are listed in table 2. Here it should be noted that the centre location is determined from the streamline plot, i.e. the centre of the closed streamlines,
$H_b/h$, is evaluated from the maximum height of the curve enveloping the separation bubble, and
$L_b/h$ is determined from the separation and reattachment points, i.e.
$L_b=x_r-x_s$. Rapp & Manhart (Reference Rapp and Manhart2011) experimentally found the height of the recirculation zone independent of the Reynolds number for
$Re_h \geqslant 10\,595$. From the present study, the height of the separation bubble decreases gradually with increasing
$Re_h$ for
$Re_h \geqslant 10\,595$. The length of the separation bubble for
$Re_h \geqslant 10\,595$ also decreases with increasing
$Re_h$ but is stronger than the variation of
$H_b$ (see table 2). Similar behaviour of the bubble length with increasing Reynolds number has also been observed in the experimental research of the backward-facing step flow by Armaly et al.(Reference Armaly, Durst, Pereira and Schönung1983) and also DNS of flow over a bump by Mollicone et al. (Reference Mollicone, Battista, Gualtieri and Casciola2017). Overall, the present numerical results compare well with the experimental results of Kähler et al. (Reference Kähler, Scharnowski and Cierpka2016) and WRLES of Fröhlich et al. (Reference Fröhlich, Mellen, Rodi, Temmerman and Leschziner2005), although the bubble length is
$6.6\,\%$ larger than the experimental data for
$Re_h=33\,000$.
4.3. Scaling of streamwise velocity profiles in the separation zone
According to the analysis of separated recirculating flow by Simpson (Reference Simpson1983) and Devenport & Sutton (Reference Devenport and Sutton1991), it is suggested that the characteristic scales for velocity and length within the recirculation region are, respectively, the magnitude of the maximum back-flow velocity $U_N$ and its distance away from the wall
$N$. The Simpson equation (Simpson Reference Simpson1983), covering the region
$0<y_n/N<1$, gives
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn20.png?pub-status=live)
where $A$ is an empirical constant,
$\bar u_s$ is the mean streamwise velocity (parallel to the bottom wall) along the wall-normal direction and
$y_n$ is the wall-normal coordinate. Simpson (Reference Simpson1983) suggested
$A=0.3$ based on the experimental research of a separating turbulent boundary layer with reattachment, while Dianat & Castro (Reference Dianat and Castro1989) found better correlations with
$A=0.235$ in another adverse pressure gradient turbulent boundary layer experiment. Furthermore, Devenport & Sutton (Reference Devenport and Sutton1991) adopted
$A$ ranging from 0.242 to 0.426 to scale the back-flow velocity profiles in a pipe with sudden expansion.
To compare with the above scaling formula, isolines of mean back-flow velocity in the recirculation region from both WMLES and WRLES are shown in figure 13. The size of the recirculation region decreases with increasing $Re_h$. The maximum back-flow velocities from the present WMLES compare well with those from WRLES. The streamwise mean velocity profiles along wall-normal lines at four locations (see figure 13) are extracted and plotted in figure 14, with
$x/h=1.0$ (
$P1$),
$x/h=1.5$ (
$P2$),
$x/h=2.0$ (
$P3$) and
$x/h=2.5$ (
$P4$). Qualitatively, the shapes of the scaled velocity profiles in the recirculation region are well described by (4.5). It is found that a unique value of
$A$ could not fit all the velocity profiles well at different locations, although
$Re_h=10^5$ shows better collapse with
$A=0.3$ than the other two cases. For
$Re_h = 10\,595$ and
$33\,000$,
$A=0.3$ suggested by Simpson (Reference Simpson1983) compares favourably with the scaled velocity profiles just at
$x/h=1.5$ (
$P2$), a profile for which there are sufficient mesh points for
$y_n<N$. On the other hand, for
$P3$ and
$P4$, the mesh points from the wall to the maximum back-flow velocity are fewer, and this may lead to an unfavourable comparison. The largest deviation with
$A=0.3$ occurs at
$x/h=1.0$, which is just downstream of the separation point and closer to the zero-velocity points. This is also reported in DNS of backward-facing step flow by Le et al. (Reference Le, Moin and Kim1997). The largest best-fit value of
$A$ (0.6 for M1; 0.45 and 0.4 for M2 and R2; 0.35 for M3) occurs at
$x/h=1.0$ (
$P1$), while the smallest best-fit value of
$A$ occurs at
$x/h=1.5$ (0.3 for M1; 0.25 for M2, R2 and M3).
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_fig13.png?pub-status=live)
Figure 13. Isolines of mean back-flow velocity in the recirculation region: (a) $Re_h=10\,595$ (WMLES), (b)
$Re_h=33\,000$ (WMLES), (c)
$Re_h=33\,000$ (WRLES) and (d)
$Re_h=10^5$ (WMLES). Velocity profiles probed at
$P1$–
$P4$ (wall-normal dashed lines) are extracted for scaling in figure 14.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_fig14.png?pub-status=live)
Figure 14. Scaling of the mean streamwise velocity profiles in the back-flow region for (a) $Re_h=10\,595$ (WMLES), (b)
$Re_h=33\,000$ (WMLES), (c)
$Re_h=33\,000$ (WRLES) and (d)
$Re_h=10^5$ (WMLES):
$\circ$,
$x/h=1.0$
$(P1)$;
$\square$,
$x/h=1.5$
$(P2)$;
$\diamond$,
$x/h=2.0$
$(P3)$;
,
$x/h=2.5$
$(P4)$, where
$P1$–
$P4$ are wall-normal lines as illustrated in figure 13. Scaled velocity profiles from (4.5) given by Simpson (Reference Simpson1983):
,
$A=0.3$;
,
$A=0.6$ in (a),
$A=0.45$ in (b),
$A=0.4$ in (c) and
$A=0.35$ in (d);
,
$A=0.45$ in (a),
$A=0.25$ in (b–d).
4.4. Pressure fluctuations
Figure 15 shows the root mean square (r.m.s.) value of pressure fluctuations in the entire channel (contour plot) and near-wall domain (wall pressure fluctuations on the bottom wall, i.e. $p_{w,rms}$). In the region around the separation point, the maximum pressure fluctuations occur very close to the wall. Inside the separation bubble, the pressure fluctuations are significantly enhanced away from the wall, and the maximum pressure fluctuations occur in the separated shear layer and turn around the bubbles (see figure 15a–c). This was also observed by Na & Moin (Reference Na and Moin1998b) in DNS of flat-plate turbulent boundary layer flows with separation and reattachment, and they attributed this to the movement of those vortical structures from upstream of separation. After the reattachment, the flow starts to redevelop and the location of maximum r.m.s. pressure fluctuations moves towards the wall, and the magnitude is much smaller than that in the separated shear layer. It should be noted that the area of high-intensity pressure fluctuations
$(\!p_{rms}/\rho U_b^2 \geqslant 0.06)$ is decreased with increasing
$Re_h$, which shows a similar trend to the separation bubble size (see also table 2). This is likely due to the dominant interaction between turbulence and mean shear in this region, which is an ‘effective’ source for pressure fluctuations (Willmarth Reference Willmarth1975; Simpson Reference Simpson1996; Hu, Reiche & Ewert Reference Hu, Reiche and Ewert2017).
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_fig15.png?pub-status=live)
Figure 15. Contours of $p_{rms}/\rho U_b^2$ superimposed on time- and spanwise-averaged streamlines: (a)
$Re_h=10\,595$, (b)
$Re_h=33\,000$ and (c)
$Re_h=10^5$. Distributions of the normalized wall pressure fluctuations: (d)
$p_{w,rms}/\rho U_b^2$, (e)
$p_{w,rms}/-\rho \overline {u'v'}_{max}$ and (f)
$p_{w,rms}/\rho \overline {v'v'}_{max}$.
,
$Re_h=10\,595$;
,
$Re_h=33\,000$;
,
$Re_h=10^5$.
To further assess the effects of Reynolds number on the wall pressure fluctuations, some scaling relations are utilized to understand the pressure-producing mechanisms. The wall pressure fluctuations normalized by the reference dynamic pressure $\rho U_b^2$ and local maximum Reynolds stress
$\rho \overline {u'v'}_{max}$ and
$\rho \overline {v'v'}_{max}$ are shown in figures 15(d)–15(f), respectively. There is less variation of wall pressure fluctuations when normalized by
$\rho {\overline {u'_iu'_j}}_{max}$ than by
$\rho U_b^2$. Many experimental and numerical investigations have shown that the local maximum Reynolds shear stress is a better scale for normalizing wall pressure fluctuations in separated turbulent boundary layers (Simpson et al. Reference Simpson, Ghodbane and McGrath1987; Ji & Wang Reference Ji and Wang2012; Abe Reference Abe2017).
Figure 15(d) shows the the r.m.s. value of wall pressure fluctuations normalized by the reference dynamic pressure $p_{w,rms}/\rho U_b^2$. Two local peaks could be noted: one at the separation point (
$p_{w,rms}/\rho U_b^2 = 0.06 \sim 0.07$) and the other close to the reattachment (
$p_{w,rms}/\rho U_b^2 = 0.05 \sim 0.06$), i.e.
$0.7 \leqslant (x-x_s)/L_b \leqslant 0.8$. The magnitudes of these peak values increase gradually with
$Re_h$. In the region of incipient detachment, the wall pressure fluctuations decrease rapidly due to the existence of the constant-pressure region in the bubble (see figures 4a and 7a). Downstream of the wall pressure fluctuation peak near the reattachment, the r.m.s. pressure decreases up to approximately half reattachment lengths. Mabey (Reference Mabey1972) summarized pressure fluctuations of step-induced separation and reattaching flows and found that the peak value of the wall pressure fluctuations occurs near the reattachment, i.e.
$p_{w,rms}/\rho U_b^2 \approx 0.06$ for the backward-facing step flow. Kiya & Sasaki (Reference Kiya and Sasaki1983) also reported a local peak value of
$p_{w,rms}/\rho U_b^2 \approx 0.06$ near the reattachment for separated flow along the side of a blunt flat plate. This value is close to our results, which is insensitive to fairly wide changes in Reynolds number and different types of separation bubbles. Na & Moin (Reference Na and Moin1998a) proposed that these peaks are partly due to the wandering of the reattachment location (further discussion of the near-wall streamlines is presented in § 4.5) and the wide variation of turbulence structures impinging on the wall at reattachment.
Figures 15(e) and 15(f) show the the r.m.s. values of wall pressure fluctuations normalized by the local maximum Reynolds stress $p_{w,rms}/\rho \overline {u'v'}_{max}$ and
$p_{w,rms}/\rho \overline {v'v'}_{max}$, respectively. Good collapse, which seems to be independent of
$Re_h$, is obtained from the separation point until some distance beyond the reattachment, i.e.
$0 \leqslant (x-x_s)/L_b \leqslant 1.5$. In experimental and numerical investigations of flat-plate turbulent boundary layer flows involving separation and reattachment, Simpson et al. (Reference Simpson, Ghodbane and McGrath1987) and Na & Moin (Reference Na and Moin1998b), respectively, found that the wall pressure fluctuations normalized by
$\rho \overline {u'v'}_{max}$ lead to near plateau, i.e.
$p_{w,rms}/-\rho \overline {u'v'}_{max} = 2.5$–
$3$, in the separated region, while reducing to be about
$1.8$ in the near-reattachment region. Ji & Wang (Reference Ji and Wang2012) pointed out that
$\rho \overline {v'v'}_{max}$ is a better scale in the reattached region (slightly less than
$L_b$) rather than
$\rho \overline {u'v'}_{max}$ in both forward- and backward-facing step flows. Moreover, they found that
$p_{w,rms}/\rho \overline {v'v'}_{max} \approx 1.35$ near the reattachment and farther downstream, when the step height is sufficiently large to generate a strong separated shear layer. Abe (Reference Abe2017) performed DNS of flat-plate turbulent boundary layer flows with large adverse and favourable pressure gradients, and found that
$p_{w,rms}/\rho \overline {v'v'}_{max} \approx 1.2$ near the reattachment independently of the Reynolds number and pressure gradient. In the present simulations, we find that both
$p_{w,rms}/-\rho \overline {u'v'}_{max}$ and
$p_{w,rms}/\rho \overline {v'v'}_{max}$ approach unambiguously some constant value in the vicinity of the reattachment (
$0.9 \leqslant (x-x_s)/L_b \leqslant 1.5$), i.e.
$p_{w,rms}/-\rho \overline {u'v'}_{max} \approx 1.7$ and
$p_{w,rms}/\rho \overline {v'v'}_{max} \approx 1.1$ independently of
$Re_h$ (see figure 15e,f). The reason for this difference is possibly associated with the convective nature of a hill-induced separation bubble near the reattachment where the interactions of the separated shear layer and turbulence, rather than the strong pressure gradient, play a dominant role in the wall pressure fluctuations.
4.5. Instantaneous skin friction fields: on instantaneous bubbles
The mean streamwise skin friction coefficients along the bottom wall provide useful information for interpreting this flow involving separation and reattachment behaviour. However, the instantaneous near-wall flow structure is more complex, and should be analysed more carefully.
Figure 16 depicts the instantaneous near-wall flow field (spanwise-averaged skin friction, $C_f$; skin friction vector field,
$\boldsymbol {C_f}=(C_{f},C_{f,z})$ (
$C_f$ and
$C_{f,z}$ denote the streamwise and spanwise skin friction coefficients, respectively), also referred to as skin friction lines; and probability density function (PDF) of the velocity angle,
$\theta$) at
$t=456.6h/U_b$ and
$t=522.9h/U_b$ for
$Re_h=10\,595$. These two time instants are chosen to correspond to a flow state where separation occurs near the hill crest or not. The instantaneous spanwise-averaged
$C_f$ is fluctuating in both attached and separated zones at most locations except
$x/h=8.2$–
$8.4$, which is different from the smooth time- and spanwise-averaged
$C_f$ profile. At
$t=456.6h/U_b$, the first separation occurs at
$x/h \approx 0.2$ (see
$\boldsymbol {C_f}$ vector field in figure 16a), also visible in the PDF(
$\theta$) with some portion of reverse flow (
$\theta \geqslant {\rm \pi}/2$). However, the first separation point detected from the spanwise-averaged
$C_f$ plot is slightly larger than that. The flow in the mean attached zone also experiences separation,
$x/h=6.0\text {--}7.6$ and
$x/h>8.9$ from the plots of
$\boldsymbol {C_f}$ and PDF(
$\theta$). The spanwise-averaged
$C_f$ fluctuates around
$C_f=0$ from
$x/h=6.6$ to
$x/h=7.5$, which implies separation and attachment in this region in spite of narrower separation zone due to the spanwise-averaging approach. The mean separation zone on the windward side of the hill,
$x/h \approx 7.0\text {--}7.4$, as reported for other DNS and LES (Mellen et al. Reference Mellen, Fröhlich, Rodi, Deville and Owens2000; Fröhlich et al. Reference Fröhlich, Mellen, Rodi, Temmerman and Leschziner2005; Breuer et al. Reference Breuer, Peller, Rapp and Manhart2009; Diosady & Murman Reference Diosady and Murman2014; Krank et al. Reference Krank, Kronbichler and Wall2018), is visible in the instantaneous field. This also confirms the high level of spanwise velocity fluctuations in the post-reattachment zone as reported by Breuer et al. (Reference Breuer, Peller, Rapp and Manhart2009), which results in the ‘splatting’ of large-scale eddies originating from the shear layer and convecting towards the windward slope. They also reported a very small counterclockwise rotating structure with positive wall shear stress for
$Re_h \leqslant 1400$ in the region
$x/h \approx 0.6\text {--}0.8$ at the leeward side of the hill, and this could also be observed in the present instantaneous field although the size is much smaller. Meanwhile, the reported very small separation zone on the hill crest is also found in figure 16(b). At
$t=522.9h/U_b$, the spanwise-averaged
$C_f$ is fluctuating around
$C_f=0$ from
$x/h=3.6$ to
$x/h=6.9$, followed by a very small separation zone at
$x/h \approx 7.0\text {--}7.5$. In the acceleration zone,
$x/h=8.0\text {--}8.9$, the skin friction lines are more uniform than in other regions even with wall curvature. For flow close to the mean locations of separation (
$x/h=0.22$) and reattachment (
$x/h=4.59$), the flow varies strongly with velocity angle covering the range from
$0$ to
${\rm \pi} $, which implies that the classic wall model assuming logarithmic law may not be suitable for such type of flow.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_fig16.png?pub-status=live)
Figure 16. The near-wall flow field on the bottom wall from WMLES for $Re_h=10\,595$: (a)
$t=456.6h/U_b$ and (b)
$t=522.9h/U_b$. The top, middle and bottom part of each panel show, respectively, the spanwise-averaged skin friction
$C_f$ (the dashed line denotes
$C_f=0$), the skin friction lines and the PDF of the velocity angle
$\theta$ along the
$x$ direction.
Figure 17 shows the instantaneous near-wall flow field at $t=186.7h/U_b$ and
$t=267.9h/U_b$ for
$Re_h=33\,000$ (WRLES). Similar to the
$Re_h=10\,595$ case, the spatially fluctuating behaviour of the spanwise-averaged
$C_f$ with higher intensity is observed. The very small separation bubble at the hill top detected from the present WRLES is visible in figure 17(a),
$x/h \approx -0.15$ to
$0.15$. At
$t=186.7h/U_b$, the very small counterclockwise rotating structure with positive wall shear stress ranges from
$x/h \approx 0.15$ to
$x/h \approx 0.22$, with larger size than that for
$Re_h=10\,595$. Excluding this very small structure, flow separation occurs from the hill top to
$x/h \approx 5.3$. At
$t=267.9h/U_b$, the very small separation bubble on the windward side of the hill is observed, ranging from
$x/h \approx 6.8$ to
$x/h \approx 7.0$, from the spanwise-averaged
$C_f$ in figure 17(b). This is also visible in the plots of
$\boldsymbol {C_f}$ and PDF(
$\theta$), but covering a slightly larger range. In the acceleration zone,
$x/h \approx 7.8\text {--}8.8$, the skin friction lines are more uniform than in other regions. For flow close to the mean locations of separation (
$x/h=0.27$) and reattachment (
$x/h=3.96$), the variations of flow directions are much stronger than in other regions, similar to that observed for
$Re_h=10\,595$.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_fig17.png?pub-status=live)
Figure 17. The near-wall flow field on the bottom wall from WRLES for $Re_h=33\,000$: (a)
$t=186.7h/U_b$ and (b)
$t=267.9h/U_b$. The top, middle and bottom part of each panel show, respectively, the spanwise-averaged skin friction
$C_f$ (the dashed line denotes
$C_f=0$), the skin friction lines and the PDF of the velocity angle
$\theta$ along the
$x$ direction.
Figure 18 shows the instantaneous near-wall flow field at $t=307.4h/U_b$ and
$t=326.9h/U_b$ for
$Re_h=10^5$. It is also noted that the spanwise-averaged
$C_f$ fluctuates along the bottom wall, even in the region close to the
$C_f$ peak. For this high-
$Re_h$ case, a major feature that is different from the previous two cases is the scale separation of the near-wall flow. At both instants, small-scale reversal flows are more concentrated, leaving sufficient space for detached flow development. In the post-reattachment region (
$x/h > 3.57$), the portion of reverse flow is smaller than for the other two cases, as shown in figure 18(a,b). Thus, the separation/reattachment point in the spanwise-averaged
$C_f$ plot is more clear. This is not observed in lower-
$Re_h$ cases, where most of the downhill region is at the vicinity of separation. Small-scale reverse-flow structures on the top and windward side of the hill are also observed from the skin friction portraits and PDF(
$\theta$), missing in the spanwise-averaged
$C_f$ plot, but the portion of reverse flow is much smaller than for
$Re_h=10\,595$ and
$33\,000$.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_fig18.png?pub-status=live)
Figure 18. The near-wall flow field on the bottom wall from WMLES for $Re_h=10^5$: (a)
$t=307.4h/U_b$ and (b)
$t=326.9h/U_b$. The top, middle and bottom part of each panel show, respectively, the spanwise-averaged skin friction
$C_f$ (the dashed line denotes
$C_f=0$), the skin friction lines and the PDF of the velocity angle
$\theta$ along the
$x$ direction.
5. Flow at the top wall: comparisons with plane channel flows
In this section, we discuss the similarities and differences between the plane channel flow and the flow field in the vicinity of the flat top wall. The mean flow field in this region is averaged in the streamwise direction since it is gradually varying along the top wall.
5.1. Skin friction
$\overline {C}_f$
The first quantity of interest is the skin friction $\overline {C}_f$ on the wall top, which may be described by some empirical friction law. Two types of formulas are adopted for comparisons: the logarithmic skin friction law (Durand Reference Durand1935),
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn21.png?pub-status=live)
where $\mathscr {K}_1=0.40$ and
$C=-0.03$ as suggested by Schultz & Flack (Reference Schultz and Flack2013), and the power law given by Dean (Reference Dean1978),
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn22.png?pub-status=live)
where $Re_m$ is based on the channel depth (
$d$) and bulk mean velocity (
$U_b$), i.e.
$Re_m = U_b d/\nu$. In plane channel flows, the channel depth is twice the channel half-height. However, in the present case, the channel is partly occupied by a submerged hill, and
$d=H-h=2.035h$ is adopted here. Figure 19 shows the LES-predicted
$\overline {C}_f$, together with other experimental and numerical results (DNS and LES) for plane channel flows. The power law (5.2) matches slightly better, with maximum relative error of
$5.4\,\%$ at
$Re_h=10\,595$, while logarithmic law (5.1) produces maximum relative error of
$9.4\,\%$. For
$Re_h=33\,000$, both formulas are in good agreement with the predicted
$\overline {C}_f$. The
$Re_h = 10^5$ case shows better agreement with the logarithmic law, in accordance with other investigations at higher
$Re_m$. The deviations from the logarithmic law may be attributed to the separated shear layer after the hill crest, which exerts a drag effect on the top wall. Nevertheless, the skin friction on the top wall seems to follow the relationship with
$Re_m$ as in plane channel flows.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_fig19.png?pub-status=live)
Figure 19. Variation of skin friction coefficient with $Re_m$:
$\diamond$, DNS from Moser, Kim & Mansour (Reference Moser, Kim and Mansour1999);
, DNS from Del Álamo et al. (Reference Del Álamo, Jiménez, Zandonade and Moser2004);
$\circ$, experimental data from Schultz & Flack (Reference Schultz and Flack2013);
$\square$, DNS from Lee & Moser (Reference Lee and Moser2015);
, LES from Giovanetti, Hwang & Choi (Reference Giovanetti, Hwang and Choi2016);
$\bullet$, results from the present simulations;
, logarithmic skin friction law (5.1);
, power skin friction law (5.2).
5.2. Mean velocities and turbulent intensities
The streamwise mean velocity in a plane channel follows a logarithmic profile with distance from the wall, $y$:
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn23.png?pub-status=live)
where $\bar {u}^+=u/u_{\tau }$,
$y^+=yu_{\tau }/\nu$ and
$B$ is a parameter that depends on the roughness of the wall. Townsend (Reference Townsend1976) proposed the attached-eddy hypothesis, which leads to a logarithmic profile for the streamwise turbulence intensity:
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn24.png?pub-status=live)
where ${\overline {u^{\prime }u^{\prime }}}^+={\overline {u^{\prime }u^{\prime }}}/{u^2_{\tau }}$ and
$\delta$ is the boundary layer thickness (or pipe radius, or channel half-height); here
$\delta =d/2$ is chosen to be compatible with the skin friction analysis. One of the difficulties in clearly demarcating the logarithmic region is that the streamwise mean velocity profile deviates slowly from (5.3), as argued by Marusic et al. (Reference Marusic, Monty, Hultmark and Smits2013). There are a variety of estimates for the bounds of the logarithmic region (Wei et al. Reference Wei, Fife, Klewicki and McMurtry2005; Eyink Reference Eyink2008; Klewicki, Fife & Wei Reference Klewicki, Fife and Wei2009; Marusic et al. Reference Marusic, Monty, Hultmark and Smits2013), and most of them relate the bound to the boundary layer thickness or
$Re_{\tau }$ (
$=u_{\tau }\delta /\nu$). We choose an estimate of the lower bound as
$y^+=2.6Re_{\tau }^{1/2}$, suggested by Klewicki et al. (Reference Klewicki, Fife and Wei2009), and the upper bound at
$y^+=0.15Re_{\tau }$, which is adopted by Marusic et al. (Reference Marusic, Monty, Hultmark and Smits2013) in analysis of wall turbulence.
Figures 20$(Re_h=10\,595, Re_{\tau }=597)$, 21
$(Re_h=33\,000, Re_{\tau }=1540)$ and 22
$(Re_h=10^5,Re_{\tau }=4000)$ show the scaling of the mean velocity and turbulence intensity profiles from the top wall, together with the DNS datasets from plane channel flows for comparisons. For an interval of
$y^+$, the match between the profiles of both the mean flow and turbulence intensities and the logarithmic scaling (5.3) and (5.4) is good for both Reynolds numbers. As mentioned above, the bound of the logarithmic region for the mean flow is difficult to discern; however, the streamwise turbulence intensity deivates more abruptly than the mean flow outside the estimated logarithmic region (
$2.6Re_{\tau }^{1/2} \leqslant y^+ \leqslant 0.15Re_{\tau }$). It is evident that the bound estimate of the logarithmic region is slightly smaller than actual bounds (see also figures 20–22), especially for
$Re_h=10\,595$. One possible reason is that the estimates based on high-Reynolds-number flows (
$Re_{\tau }>20\,000$) may not be good enough for the present cases.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_fig20.png?pub-status=live)
Figure 20. Scaling of the streamwise-averaged (a) mean velocity profile and (b) turbulence intensity profile from the top wall for $Re_h=10\,595$:
$\circ$, present WMLES;
, DNS from Moser et al. (Reference Moser, Kim and Mansour1999), at
$Re_{\tau }=590$;
, DNS from Del Álamo et al. (Reference Del Álamo, Jiménez, Zandonade and Moser2004), at
$Re_{\tau }=950$. The solid diagonal line in (a) denotes the generalized logarithmic profile defined in (5.3) with
$\mathscr {K}_1=0.384$,
$B=4.6$ from Inoue & Pullin (Reference Inoue and Pullin2011), and the solid diagonal line in (b) denotes the distribution given in (5.4) with
$A_1=1.74$,
$B_1=0.82$. The dashed vertical lines indicate the bounded region
$2.6Re_{\tau }^{1/2} \leqslant y^+ \leqslant 0.15Re_{\tau }$.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_fig21.png?pub-status=live)
Figure 21. Scaling of the streamwise-averaged (a) mean velocity profile and (b) turbulence intensity profile from the top wall for $Re_h=33\,000$:
$\square$, present WMLES;
$\circ$, present WRLES. The dot-dashed lines denote DNS from Lee & Moser (Reference Lee and Moser2015) at
$Re_{\tau }=1000$ (
) and
$Re_{\tau }=2000$ (
). The solid diagonal line in (a) denotes the generalized logarithmic profile defined in (5.3) with
$\mathscr {K}_1=0.384$,
$B=4.6$ from Inoue & Pullin (Reference Inoue and Pullin2011), and the solid diagonal line in
$(b)$ denotes the distribution given in (5.4) with
$A_1=0.88$,
$B_1=4.56$. The dashed vertical lines indicate the bounded region
$2.6Re_{\tau }^{1/2} \leqslant y^+ \leqslant 0.15Re_{\tau }$.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_fig22.png?pub-status=live)
Figure 22. Scaling of the streamwise-averaged (a) mean velocity profile and (b) turbulence intensity profile from the top wall for $Re_h=10^5$:
$\circ$, present WMLES;
, DNS from Bernardini, Pirozzoli & Orlandi (Reference Bernardini, Pirozzoli and Orlandi2014), at
$Re_{\tau }=4000$;
, DNS from Lee & Moser (Reference Lee and Moser2015), at
$Re_{\tau }=5200$. The solid diagonal line in (a) denotes the generalized logarithmic profile defined in (5.3) with
$\mathscr {K}_1=0.384$,
$B=4.6$ from Inoue & Pullin (Reference Inoue and Pullin2011), and the solid diagonal line in (b) denotes the distribution given in (5.4) with
$A_1=1.26$,
$B_1=3.02$. The dashed vertical lines indicate the bounded region
$2.6Re_{\tau }^{1/2} \leqslant y^+ \leqslant 0.15Re_{\tau }$.
The $Re_h=10\,595$ case is closer to plane channel flow at
$Re_{\tau }=590$ (Moser et al. Reference Moser, Kim and Mansour1999) as inferred from
$\overline {C}_f$ plot in figure 19. The mean velocity profiles are in good agreement with this plane channel flow and (5.3) for
$y^+ \leqslant 180$; however, the streamwise turbulence intensities are more close to
$Re_{\tau }=950$ (Del Álamo et al. Reference Del Álamo, Jiménez, Zandonade and Moser2004), and curve fitting using (5.4) gives
$A_1=1.74$, larger than the Townsend–Perry constant
$A_1 = 1.26$ (Marusic etal. Reference Marusic, Monty, Hultmark and Smits2013) for
$Re_{\tau }>20\,000$. For
$Re_h=33\,000$, it is close to
$Re_{\tau }=1000$ (Lee & Moser Reference Lee and Moser2015). The mean velocity profiles compare well with plane channel flows for
$y^+ \leqslant 400$, while the streamwise turbulence intensity is even larger than
$Re_{\tau }=2000$ (Lee & Moser Reference Lee and Moser2015), and
$A_1=0.88$ smaller than the recommended value. For
$Re_h=10^5$, it is close to
$Re_{\tau }=4000$ (Bernardini et al. Reference Bernardini, Pirozzoli and Orlandi2014), but the streamwise turbulence intensity is larger, similar to
$Re_h=33\,000$. The mean velocity profiles compare well with plane channel flows (
$y^+ \leqslant 650$) with the peak location of the streamwise turbulence intensity slightly off from the plane channel flow. Moreover, the streamwise turbulence intensity is even larger than
$Re_{\tau }=5200$ (Lee & Moser Reference Lee and Moser2015), and the Townsend–Perry constant
$A_1=1.26$ fits well the LES data. These results would imply that
$A_1$ is not a universal constant at insufficiently high Reynolds number. The hill in the channel acts as a perturbation source, which strongly enhances the turbulence intensity, and the flow far away from the top wall is also affected by the developed shear layer. These would also change the eddy structures and thus the associated parameters in the eddy intensity functions (Marusic & Monty Reference Marusic and Monty2019). These would also explain why the mean flow and streamwise turbulence intensity far away from the top wall deviate sharply from plane channel flows. Overall, the flow in the region close to the top wall behaves similarly to plane channel flows, following the classic logarithmic scaling (5.3) and (5.4). However, the hill in the channel perturbs the entire flow field, and the separated shear layer transfers these effects to the upper flow region.
6. Conclusion
We have presented results from LES of turbulent flow in a channel constricted by streamwise periodically distributed hills at $Re_h=10\,595,33\,000$ and
$10^5$. To overcome the high computational cost of resolving the small-scale turbulent scales in the near-wall region (both top and bottom walls), we couple the generalized virtual wall model (near-wall region) and the stretched spiral vortex model (outer LES resolved region) to formulate the WMLES framework. The numerical methodology is based on fourth-order spatial differencing, with a multigrid Poisson solver for the pressure utilizing Gauss–Seidel line smoothers.
The WMLES results for $Re_h=10\,595$ and
$33\,000$ are verified and validated against both WRLES and experimental data, including the pressure and skin friction coefficients, mean velocity and Reynolds stress profiles. We note that capturing separation and reattachment in complex flows still remains a challenge for the turbulence simulation community. In our WMLES, the locations of separation and reattachment are in good agreement with experiments. Furthermore, comparisons of the mean velocity profiles with experiments show good agreement, while profiles of Reynolds stress components are in reasonable agreement with experiments. For
$Re_h=10^5$, the largest
$Re_h$ case to the best of our knowledge, the locations of reattachment compare well with the scaling law given by Kähler et al. (Reference Kähler, Scharnowski and Cierpka2016), and the peak values of the skin friction coefficient are in good agreement with the best-fit scaling law based on the DNS/LES results from
$Re_h=700$ to
$10\,595$.
Based on these verifications and validations, we also investigated the effects of Reynolds number inside the recirculation zone. We find that the normalized skin friction coefficient follows ‘laminar-like’ behaviour, which is similar to the correlation in backward-facing step flows. The centre coordinate and length of the separation bubble behind the hill decrease with $Re_h$, while the dependence of the bubble height on the Reynolds number is not that strong. Meanwhile, the mean velocity profiles inside the recirculation zone are compared with the Simpson equation, and we find that a unique empirical constant
$A$ does not fit all the velocity profiles at different locations, especially near the separation point. Furthermore, the r.m.s. value of pressure fluctuations is examined in the entire channel and the maximum pressure fluctuations occur in the separated shear layer. The wall pressure fluctuations normalized by the reference dynamic pressure,
$\rho U_b^2$, indicate local peaks near the separation and reattachment point, i.e.
$p_{w,rms}/\rho U_b^2 \approx 0.06$, although with slight increase with
$Re_h$. The wall pressure fluctuations show less of a variation around the recirculation zone when scaled by the local maximum Reynolds stress
$-\rho \overline {u'v'}_{max}$ and
$\rho \overline {v'v'}_{max}$. This normalization produces some constant value in the vicinity of the reattachment, i.e.
$p_{w,rms}/-\rho \overline {u'v'}_{max} \approx 1.7$ and
$p_{w,rms}/\rho \overline {v'v'}_{max} \approx 1.1$, independently of
$Re_h$.
We examined the details of unsteady behaviour of separation and reattachment by examining the surface streamlines. An often-ignored feature in flows over periodic constrictions is the flow field near the top wall. We have compared the flow field at the top wall with planar channel flows. The skin friction at the top wall is in agreement with previously proposed empirical relations. Also, the mean velocity profiles and turbulent intensities are compared with planar channel flows, showing the existence of a logarithmic layer. The periodic hills act as perturbations to planar channel flows and the shear layer from the crest of the hill has an effect in terms of drag on the top wall.
Acknowledgements
The Cray XC40 Shaheen II at KAUST was used for all simulations reported. This research was partially supported under KAUST OCRF URF/1/1394-01 and under baseline research funds of R.S.
Declaration of interests
The authors report no conflict of interest.
Appendix A. Stretched vortex SGS model
The stretched vortex SGS model has been previously widely deployed in LES of wall-bounded turbulent flows by Pullin and co-workers. Here, for the sake of completeness, we present the essential features of this structure-based SGS model, which assumes that the turbulent fine scales are composed of tube-like structures with spiral vortices (Lundgren Reference Lundgren1982). In each computational cell, the ensemble dynamics is dominated by a stretched vortex with orientation $\boldsymbol {e}^\nu$, taken from a delta-function probability density function (Misra & Pullin Reference Misra and Pullin1997). Consequently, the SGS stress tensor
${T}_{ij}$ is modelled as (Chung & Pullin Reference Chung and Pullin2009)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn25.png?pub-status=live)
where $K$ is the SGS kinetic energy,
$k_c ={\rm \pi} /{\rm \Delta} _c$ is the cut-off wavenumber (
${\rm \Delta} _c = ({\rm \Delta} _x {\rm \Delta} _y {\rm \Delta} _z)^{1/3}$) and
$E(k)$ is the SGS energy spectrum given by Lundgren (Reference Lundgren1982). The integration of the energy spectrum gives
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn26.png?pub-status=live)
where $\Gamma$ is the incomplete gamma function,
$\tilde {a}=e^v_i e^v_j \tilde {{S}}_{ij}$ is the stretching felt along the subgrid vortex axis imposed by the resolved scales and
$\tilde {{S}}_{ij}$ is the resolved SGS strain tensor. The composite parameter
$\mathcal {K}_0'$ is obtained dynamically by structure–function matching at the grid-scale cut-off (Voelkl, Pullin & Chan Reference Voelkl, Pullin and Chan2000), i.e.
$\mathcal {K}_0' = \langle F_2 \rangle / \langle Q( \kappa _c, d) \rangle$, where
$\langle .\rangle$ denotes a local-averaging operator computed from the neighbouring
$26$ points,
$F_2$ is the second-order velocity structure function from the resolved LES field and the calculation of
$Q( \kappa _c, d)$ is similar to that in Voelkl et al. (Reference Voelkl, Pullin and Chan2000) and Chung & Pullin (Reference Chung and Pullin2009) with
$\kappa _c=r/{\rm \Delta} _c$ and
$r$ the distance from the neighbouring point to the vortex axis.
Appendix B. Wall modelling in complex geometry
Starting with the Navier–Stokes equations in the generalized curvilinear coordinates ((2.2) and (2.3)), we apply near-wall filtering along with the assumption of inner scaling to derive an ODE governing the local wall-normal velocity gradient, and a slip Dirichlet boundary condition for velocity at a virtual wall.
B.1. Derivation of the wall shear stress
Following Chung & Pullin (Reference Chung and Pullin2009), we define two near-wall filters in the near-wall region:
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn27.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn28.png?pub-status=live)
where (B 1) and (B 2) define the wall-parallel filter and the wall-normal averaging filter, respectively. Applying the wall-parallel filter (B 1) to the momentum equations (2.2), we obtain
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn29.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn30.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn31.png?pub-status=live)
where $(u,v,w)=(u_1,u_2,u_3)$ denote the velocity components in Cartesian coordinates.
We assume $\tilde q$, defined in (2.6a), follows inner scaling, i.e.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn32.png?pub-status=live)
where $F(y_n^+)$ is an unknown function of the normal distance from the solid wall in wall units. The following ODEs can then be derived:
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn33.png?pub-status=live)
Applying the wall-normal averaging filter (B 2) to the governing equation of $\tilde q$ results in
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn34.png?pub-status=live)
where $\tilde q|_{h'} = u_\tau F(h^+)$ is the resultant velocity at a distance
$h'$ from the solid wall (see figure 2). It should be noted that (B 6) is an exact consequence of (B 2) and (B\,4a,b). Moreover, an explicit form of
$F(y_n^+)$ is not required owing to cancellation.
From the definition of $\tilde q$, (2.6a,b), we have
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn35.png?pub-status=live)
From the definition of $\tilde u_s$, (2.7) and combing (B 3a) and (B 3b), we can obtain the streamwise momentum equation along the wall:
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn36.png?pub-status=live)
Then combining (B 3c), (B 7) and (B 8), we can obtain
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn37.png?pub-status=live)
where
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn38.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn39.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn40.png?pub-status=live)
For the purpose of modelling, we make the following two approximations within the first grid cell $0 \leqslant \eta \leqslant h'$:
(i) the velocity angle
$\theta$, i.e.
$\tilde u_s/ \tilde q$ and
$\tilde w/ \tilde q$, is constant;
(ii) the Jacobian of the transformation, i.e.
$\sqrt g$, is constant.
Then we can reduce $F_\eta$ to be
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn41.png?pub-status=live)
where
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn42.png?pub-status=live)
If the first wall layer is forced to be orthogonal, then the terms with $ {G}^{ij} (i\neq j)$ can be neglected. If the spanwise geometry is uniform, then the terms with
$\partial \xi ^i/\partial z\ (i=1,2)$ can be neglected.
Terms $F_\xi$,
$F_\zeta$ and
$M$ feature in the right-hand side of (B 17) (see below) and are unknown, which are estimated by resolved-scale quantities at the first grid point
$(y_n=h')$ above the solid wall. For example, the first term on the right-hand side of (B 10) is approximated by LES resolved-scale values at
$y_n=h'$ as
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn43.png?pub-status=live)
where
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn44.png?pub-status=live)
where $ {T}_{xx}$,
$ {T}_{xy}$ and
$ {T}_{xz}$ are the SGS stress tensor components. The other terms in
$F_\xi$ (B 10),
$F_\zeta$ (B 12) and
$M$ (B 14) are approximated in a similar manner.
Equations (B 6) and (B 7) can be algebraically manipulated to derive an ODE for $\eta _0$:
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn45.png?pub-status=live)
where
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn46.png?pub-status=live)
Rewriting the ODE we have
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn47.png?pub-status=live)
If $C_1/C_2$ is weakly dependent on
$t$, (B 19) can be solved analytically:
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn48.png?pub-status=live)
where
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn49.png?pub-status=live)
in which $t_0$ is the current time and
${\rm \Delta} t$ is the time interval in the simulations.
B.2. Slip velocity boundary conditions
Once $\eta _0$ is solved using (B 20), we complete the wall model with a slip velocity at a raised virtual wall plane located at
$\eta =h_0$, which scales with the boundary layer thickness but remains small, i.e.
$h_0 \leqslant 0.1\delta$. Typically,
$h_0$ is chosen as a small fraction of the near-wall cell size. In previous studies, both in channel flow by Chung & Pullin (Reference Chung and Pullin2009) and in turbulent boundary layer flows (zero and adverse pressure gradients) by Inoue & Pullin (Reference Inoue and Pullin2011) and Cheng et al. (Reference Cheng, Pullin and Samtaney2015),
$h_0=0.18{\rm \Delta} \eta$ is recommended. Their verifications and validations capture the near-wall flow physics well, and we follow this choice in the present research. Finally, the resultant slip velocity
$\tilde q|_{h_0}$ on the virtual wall can be modelled as in (2.12).
Appendix C. Numerical method
The governing equations (2.4) and (2.5) are discretized as
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn50.png?pub-status=live)
where ${\delta }/{\delta \xi ^m}$ represents the energy-conservative fourth-order finite-difference operator (Morinishi et al. Reference Morinishi, Lund, Vasilyev and Moin1998),
$C_i$ and
$S_i$ represent the convective term and SGS term and
$D_i$ and
$R_i$ are discrete operators for the viscous term and the pressure gradient term, respectively. These quantities are
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_eqn51.png?pub-status=live)
where the convective term is computed in the skew-symmetric form to minimize the aliasing error (Zang Reference Zang1991; Morinishi et al. Reference Morinishi, Lund, Vasilyev and Moin1998). The fractional step method (Zang, Street & Koseff Reference Zang, Street and Koseff1994; Zhang et al. Reference Zhang, Cheng, Gao, Qamar and Samtaney2015) is used to solve the governing equations. This method follows the predictor–corrector procedure, and the pressure Poisson equation is solved using the multigrid method with line-relaxed Gauss–Seidel as a smoother. The code is parallelized using standard MPI protocol. To achieve near-optimal load balancing, the mesh is divided into blocks of equal size and each of them is assigned to a uniqueprocessor. All the simulations are performed on the Shaheen-Cray XC40 at KAUST. The DNS version of the code (without the SGS stress terms) with the same method is described in Zhang et al. (Reference Zhang, Cheng, Gao, Qamar and Samtaney2015) for flow past an airfoil. The WRLES and WMLES versions of the same code were previously applied successfully in flow past cylinders (Cheng et al. Reference Cheng, Pullin and Samtaney2017, Reference Cheng, Pullin and Samtaney2018) and different airfoil shapes at various angles of attack and Reynolds numbers ranging from $10^4$ to
$2.1\times 10^6$ (Gao et al. Reference Gao, Zhang, Cheng and Samtaney2019), respectively.
Appendix D. The effect of mesh resolution
To evaluate the effects of mesh resolution on WMLES, we performed another WMLES with a coarse mesh ($N_{\xi } \times N_{\eta } \times N_z = 96 \times 48 \times 48$) at
$Re_h=10\,595$. The wall-normal mesh size
${\rm \Delta} \eta$ is doubled (both on bottom and top walls), and correspondingly the height of the virtual wall is doubled since we have fixed
$h_0= 0.18{\rm \Delta} \eta$. Bose & Park (Reference Bose and Park2018) summarized WMLES results in complex geometries, i.e. different airfoils and periodic hill configurations, and found that the proper resolution of the separated shear layer is important.
To show convergence, in figure 23 we plot both the time- and spanwise-averaged skin friction coefficient on the bottom wall and ${\overline {u^{\prime }u^{\prime }}}/{U^2_b}$ profile near the centre of the separation zone (
$x/h=2$). It is evident from the
$C_f$ plot that the WMLES with the coarser mesh results in reattachment at a smaller
$x/h$ location although the separation occurs at almost the same location, and the peak value of
$C_f$ at
$x/h \approx 8.6$ is slightly larger than for the WRLES and the WMLES with the finer mesh. The largest discrepancies occur inside the separation zone. The same behaviour can be observed in the
${\overline {u^{\prime }u^{\prime }}}/{U^2_b}$ profile (see figure 23b), where we note that the WMLES with the coarser mesh produces the largest values in the core region of the separation bubble and smallest values above the separation bubble. Overall, these comparisons indicate that the present wall model can capture the separation bubble well even with the coarse mesh, and mesh resolution should be judiciously chosen to resolve the separated shear layer.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_fig23.png?pub-status=live)
Figure 23. (a) Skin friction coefficient $C_f$ along the bottom and (b)
${\overline {u^{\prime }u^{\prime }}}/{U^2_b}$ near the centre of the separation bubble (
$x/h=2$) for
$Re_h=10\,595$:
, WRLES from Fröhlich et al. (Reference Fröhlich, Mellen, Rodi, Temmerman and Leschziner2005);
, present WMLES with fine mesh;
, present WMLES with coarse mesh (case M1). In (a), the bold grey line denotes the geometry of the channel with streamwise periodic hills and the horizontal dashed line is the
$C_f=0$ line plotted for reference.
Appendix E. The contribution of the SGS stress to the Reynolds stress
To evaluate the contribution of the SGS stress to the total Reynolds stress, in figure 24 we plot the fraction of the mean SGS stress components $\overline {{T}}_{xx}$ that contributes to
${\overline {u^{\prime }u^{\prime }}}$ (
$\overline {u'u'} = \overline {\tilde {u}'\tilde {u}'} + \overline {{T}}_{xx}$). Two locations are specified,
$x/h=2$ (separation zone) and
$x/h=6$ (attached zone). In the attached zone,
$\overline {{T}}_{xx}/{\overline {u^{\prime }u^{\prime }}}$ is smaller than that in the separation zone, especially in the separated shear layer after the hill crest. Compared with the
$Re_h=10\,595$ case, the SGS stress contribution is larger in the other two higher-
$Re_h$ cases. This may also potentially explain why the the Reynolds stress comparisons for
$Re_h=33\,000$ with experimental data are not as good as for
$Re_h=10\,595$. The resolution in the separated shear layer is important for resolving the mixing process, i.e. sharp variation of the Reynolds stress in the separation zone, as mentioned in § 3.2.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20200813172526632-0418:S0022112020005121:S0022112020005121_fig24.png?pub-status=live)
Figure 24. The contribution of the SGS stress to the Reynolds stress $\overline {{T}}_{xx}/{\overline {u^{\prime }u^{\prime }}}$ at (a)
$x/h=2$ and (b)
$x/h=6$ for three
$Re_h$ cases:
,
$Re_h=10\,595$;
,
$Re_h=33\,000$;
,
$Re_h=10^5$.