0
Research Papers

# Simulation of Rotating Channel Flow With Heat Transfer: Evaluation of Closure ModelsOPEN ACCESS

[+] Author and Article Information
Alan S. Hsieh

Department of Aerospace Engineering Sciences,
Boulder, CO 80309

Sedat Biringen

Department of Aerospace Engineering Sciences,
Boulder, CO 80309

Alec Kucala

Department of Aerospace Engineering Sciences,
Boulder, CO 80309

1Corresponding author.

Manuscript received September 14, 2015; final manuscript received April 12, 2016; published online May 24, 2016. Assoc. Editor: Graham Pullan.

J. Turbomach 138(11), 111009 (May 24, 2016) (15 pages) Paper No: TURBO-15-1203; doi: 10.1115/1.4033463 History: Received September 14, 2015; Revised April 12, 2016

## Abstract

A direct numerical simulation (DNS) of spanwise-rotating turbulent channel flow was conducted for four rotation numbers: $Rob=0$, 0.2, 0.5, and 0.9 at a Reynolds number of 8000 based on laminar centerline mean velocity and Prandtl number 0.71. The results obtained from these DNS simulations were utilized to evaluate several turbulence closure models for momentum and heat transfer transport in rotating turbulent channel flow. Four nonlinear eddy viscosity turbulence models were tested and among these, explicit algebraic Reynolds stress models (EARSM) obtained the Reynolds stress distributions in best agreement with DNS data for rotational flows. The modeled pressure–strain functions of EARSM were shown to have strong influence on the Reynolds stress distributions near the wall. Turbulent heat flux distributions obtained from two explicit algebraic heat flux models (EAHFM) consistently displayed increasing disagreement with DNS data with increasing rotation rate.

<>

## Introduction

Turbulent channel flow subject to rotation in the spanwise direction is characterized by reduced turbulence levels near one wall and elevated turbulence levels near the opposite wall; these regions are known as the suction and pressure sides, respectively [1]. Subsequently, the symmetric profiles of mean velocity and Reynolds stress distributions in the nonrotating channel become asymmetric with respect to the channel centerline. The effects of spanwise rotation on momentum transport in turbulent channel flow have been well documented in previous DNS studies by Grundestam et al. [1] and Wu and Kasagi [2]. Rotation-induced body forces (Coriolis and centrifugal) generate a secondary cross flow, and consequently, a particular type of complex turbulent flow regime is developed with more than one mean flow gradient. The analyses of such complexities on the structure and parameterization of turbulence have relevance to engineering applications, such as gas turbine blade and rotating turbomachinery design, especially with respect to surface heat transfer and skin friction within the internal cooling passages [3,4].

Due to significantly reduced computational costs compared to DNS, Reynolds-averaged Navier–Stokes (RANS) and turbulent heat flux models are highly desirable if they can accurately parameterize turbulent transport. The strong influence of rotational forces on the Reynolds stress and turbulent heat flux transport equations imposes significant challenges for model development especially at high rotation numbers. In this work, we employ DNS to integrate the Navier–Stokes and energy equations and utilize this database to assess four RANS models proposed by (a) Reif et al. [5], (b) Speziale and Gatski [6], (c) Girimaji [7], and (d) Grundestam et al. [8]. We also evaluate two algebraic heat flux models proposed by (e) Younis et al. [9] and (f) Abe and Suga [10]. In addition, the pressure–strain functions proposed in Refs. [6] and [7] are investigated for their influence on the modeled Reynolds stress distributions.

## Numerical Method and Case Descriptions

The time-dependent, three-dimensional incompressible Navier–Stokes and energy equations were numerically integrated in a doubly periodic (in x and z directions) channel flow using a fractional-step method [11,12]. With all spatial coordinates nondimensionalized by the channel half-height δ and velocities by the laminar centerline velocity uc, the Navier–Stokes and energy equations read (in conservative form) Display Formula

(1)$∂ui∂xi=0$
Display Formula
(2)$∂ui∂t+∂uiuj∂xj=−∂p∂xi+1Rec∂2ui∂xj∂xj−RocεijkΩjΩuk$
Display Formula
(3)$∂θ∂t+∂(θuj)∂xj=1RecPr∂2θ∂xi∂xj$

where $Rec=ucδ/ν, Pr=ν/κ$, ν is the kinematic viscosity, κ is the thermal diffusivity, and the vector $u=〈u,v,w〉$ is composed of three velocity components in the x (streamwise), y (wall-normal), and z (spanwise) directions, respectively. The rotation number (or Rossby number) is defined as $Roc=2Ωδ/uc$, where Ω is the spanwise angular rotation vector. p is the nondimensional effective pressure ($p=p0−(1/8)Roc2r2$) which combines the static pressure (p0) and centrifugal force, and r represents the nondimensional distance away from the axis of rotation. Also, t is the nondimensional time and θ is the nondimensional temperature ($θ=[T−TL]/[TU−TL]$) with $TU$ and $TL$ representing dimensional temperatures on the upper and lower walls, respectively. Equation (2) is uncoupled from Eq. (3), and buoyancy effects are neglected such that present DNS thermal statistics can be compared to the results of similar heat transfer studies involving rotation where negligible buoyancy effects were also assumed. The Prandtl number (Pr) was kept constant at 0.71.

The flow geometry of the present DNS is shown in Fig. 1. Both the velocity and temperature flow fields are assumed to be statistically homogeneous in the streamwise (x) and spanwise (z) directions allowing periodic boundary conditions in those directions. No-slip conditions are imposed on the rigid channel walls (y = 0 and 2). An isothermal boundary condition on the dimensionless temperature θ ($θU=0.5$ and $θL=−0.5$) was imposed at the upper and lower walls. Temperature coupling was initiated once the turbulence field had reached a quasi-steady-state solution, and a linear distribution for the temperature field was assumed as an initial condition in the wall-normal direction.

The fractional-step time integration scheme used is a modified semi-implicit Adams–Bashforth/Crank–Nicolson method, which is second-order accurate in time and fourth-order accurate in space. Spatial derivatives were computed using fourth-order central finite differences in order to facilitate parallelization. The wall-normal direction incorporated an exponentially stretched and staggered mesh that clusters points near the two wall boundaries. All the solutions for the elliptic pressure equation and the spatially elliptic momentum equations are obtained using iterative methods available on the Portable, Extensible Toolkit for Scientific Computation Library at the Texas Advanced Computing Center (TACC). To validate the accuracy of the code, numerical solutions were compared with the hydrodynamic linear stability theory for plane channel (Poiseuille) flow. By measuring the growth rate and phase speed of a primary two-dimensional disturbance from the computation and comparing these values to the eigenvalues of the Orr–Sommerfeld equation, the accuracy of the code was assessed. Excellent agreement between the theoretical predictions and DNS results was observed (less than $1%$ deviation in amplitude), which gives confidence to the viability of the code for solving the unsteady incompressible Navier–Stokes equations. More information on the code scheme and validation may be found in Refs. [1114].

The full listing of simulations and their corresponding grid resolutions are found in the case descriptions (Table 1). Four simulation cases A–D were conducted for Rossby numbers $Rob=0$, 0.2, 0.5, and 0.9, such that models could be tested for a wide range of rotation rates. The Reynolds number based on the laminar centerline velocity $Rec=8000$ was kept constant which resulted in a friction Reynolds number $Reτ=200$ for the no-rotation case ($Rob=0$). The low Reynolds number was chosen such that reasonable comparisons could be made between the present DNS and other DNS studies of rotating turbulent channel flow with similar Reynolds numbers [1,15]. The asymmetric velocity distributions due to the rotational effects decreased the value of $Reτ$, which was calculated as an average between the two walls. Simulations were performed at constant mass flux which resulted in a bulk Reynolds number $Reb=ubδ/ν$, where $ub$ is the mean bulk velocity defined as

Display Formula

(4)$ub=12δ∫02δu¯dy$

where $u¯$ denotes a plane-averaged quantity. The bulk Reynolds number of the present DNS simulations was 3800. The domain lengths were selected as $Lx=4πδ, Ly=2δ$, and $Lz=2πδ$ for all the simulations such that two-point spatial autocorrelations in the streamwise and spanwise directions converged to zero at the largest separations.

Mesh independence was established by designating two high-resolution cases with grid numbers $nx×2ny×nz$ for case A ($Rob=0$) and $2nx×ny×2nz$ for case D ($Rob=0.9$). For both cases, the distributions of mean velocity, Reynolds stresses, and turbulent kinetic energy budgets compared very favorably to those of the original simulations, demonstrating that the selected meshes of the present DNS cases were mesh invariant. The grid spacing is also comparable to other DNS studies of spanwise-rotating turbulent channel flow [1,15].

Superscript “+” refers to nondimensionalization by the friction velocity, $uτ=ν∂u¯/∂y|wall$ or friction temperature, $Tτ=(κ/uτ)∂θ/∂y|wall$. The global friction velocity for the rotational cases is denoted as $uτ=uτs2/2+uτp2/2$, where $uτs$ and $uτp$ are the local friction velocities at the suction and pressure walls, respectively, an equivalent calculation of the global friction temperature was used for the thermal statistics. Unless otherwise specified, all the coordinate directions are nondimensionalized by δ. For all the cases, the governing equations were integrated until both the friction Reynolds number $Reτ$ and friction temperature converged availing a sufficiently long time window ($t+=tuτ2/ν≥1000$) to calculate statistics.

In fully resolved turbulent flow, energy spectra should demonstrate an energy density roll-off of several decades with increasing wavenumber and negligible aliasing at high wavenumbers [16]. Sample streamwise power spectra (E) for case A ($Rob=0$) are shown in Fig. 2. The energy spectra demonstrate sufficient energy drop-off and no aliasing at high wavenumbers in addition to proper characterization of both the inertial and dissipative energy scales.

In Fig. 3, turbulence quantities for simulation case A ($Rob=0$) are compared with the experimental hot-film data in Ref. [17] and particle-tracking velocimetry data (PTV) in Ref. [18]; a PTV distribution was not available for $u1′u2′¯+$. The present DNS results are shown to compare favourably, especially with respect to the PTV data for $(ui′ui′¯+)1/2$ and the hot-film data for $u1′u2′¯+$. In addition, good agreement was also shown for rotating simulation cases B–D with DNS results in Refs. [15] and [1]. Thermal statistics obtained for simulation cases A–D also compared favorably with the DNS results in Ref. [19].

Rotational effects on the mean velocity and temperature profiles are shown for simulation cases A–D in Figs. 4 and 5, respectively. In Fig. 4, the mean velocity profile is symmetric about the channel centerline (y = 1) for case A ($Rob=0$). With system rotation, the mean velocity distributions become asymmetric as the flow regime is separated into the pressure and suction regions. For rotational cases B–D, a laminarlike (parabolic) profile is observed near the suction wall (y = 0) which is characteristic of the suppressed turbulence in the suction region. As the flow progressively relaminarizes with increasing rotation number, the suction region expands. Near the pressure wall (y = 2), a constant gradient of $−2Ω$ is shown in the mean velocity profiles at all the rotation rates, which is consistent with previous DNS results in Refs. [1] and [15]. In Fig. 5, the thickness of the thermal diffusive sublayer, characterized by large near-wall temperature gradients, is broader near the suction wall than the pressure wall for rotating simulation cases B–D. As the rotation number increases, the size of the diffusive layer in the suction region increases and the mean temperature profile shifts toward the pressure wall [19].

In Fig. 6, maps of fluctuating streamwise vorticity $ωx′$ are shown for a y–z crossflow section at $x=2π$ for simulation cases A–C. In Fig. 6(a), there are spanwise arrays of alternating high- and low-speed streaks near both channel walls for case A ($Rob=0$), which also represent the local intensity of turbulence. In Figs. 6(b) and 6(c), the coherence of the high- and low-speed structures is disrupted near the suction wall for cases B ($Rob=0.2$) and C ($Rob=0.5$). With increasing rotation number in the suction region, a significant reduction in the number of turbulence structures is observed along with the expansion of a near-wall region of completely suppressed turbulent activity. In contrast to the suction region, the pressure region continues to demonstrate high levels of turbulent activity, and the coherence of turbulence structures is preserved for all the rotation rates.

In Fig. 7, time-averaged spanwise and wall-normal velocity vectors are shown for a y–z crossflow section at $x=2π$ for simulation cases A–C. In Fig. 7(a), no spanwise flow is observed in case A ($Rob=0$) as expected. In Figs. 7(b) and 7(c), large-scale vortical structures (roll cells) appear from the secondary mean flow advected by the Coriolis force for rotating cases B ($Rob=0.2$) and C ($Rob=0.5$). These counter-rotating pairs of roll cells, also known as Taylor–Gortler vortices [1], are observed to circulate flow throughout the pressure region. With increasing rotation number, the roll cells become more numerous and the size of the circulation region is reduced from increased relaminarization.

To illustrate the effects of rotation on wall shear stress and heat transfer, the dimensionless friction Reynolds ($Reτ$) and Nusselt (Nu) numbers for both channel walls are provided for the present simulation cases A–D in Table 2. The introduction of rotation is shown to initially decrease $Reτ$ on the suction wall while increasing $Reτ$ on the pressure wall. At higher rotation numbers, these $Reτ$ trends are shown to significantly weaken or even reverse in case D ($Rob=0.9$) for $Reτ$ on the pressure wall. These results correspond well with Ref. [1] which demonstrated that $Reτ$ on both walls trended toward convergence at high rotation numbers until the eventual full relaminarization of the flow regime. For the Nusselt number, a dimensionless number generally used to represent surface heat transfer, increasing system rotation is shown to continually decrease Nu on both channel walls. This trend was similarly observed in Ref. [19].

## RANS Turbulence Models

###### RANS Models: Overview.

The engineering approach for the calculation of turbulent flows generally consists of the use of Reynolds-averaged (time-averaged) Navier–Stokes (RANS) equations which require heuristic closure approximations to model the Reynolds stresses that appear in these equations as a result of the time-averaging process. Accurate parameterization of these quantities in terms of the fundamental variables of the problem, often with the application of an eddy (turbulent) viscosity coefficient, is essential for an accurate prediction of turbulent transport. RANS models are often incorporated into engineering applications but have difficulties especially for the modeling of complex flows with extra mean rates-of-strain. In this work, the current state of turbulence models with respect to spanwise-rotating turbulent channel flow is evaluated through the examination of four commonly used turbulence models in their ability to accurately characterize turbulent transport with rotational effects.

In this work for model testing, nonlinear eddy viscosity models are used due to their advantages over linear eddy viscosity models, such as the k–ε model [20] for complex turbulent flows. Linear eddy viscosity models invoke a linear constitutive relationship between the Reynolds stresses and mean flow straining field (Boussinesq approximation) and although these models have produced satisfactory predictions in two-dimensional thin shear flows, they perform poorly in complex turbulent flows with more than one mean flow velocity gradient such as the spanwise-rotating turbulent channel flow considered here [21]. In Fig. 8, the Reynolds stress distributions of the k–ε model were compared to the present DNS data and a sample nonlinear eddy viscosity model [8] for case D ($Rob=0.9$). In comparison to the nonlinear eddy viscosity model, the k–ε model demonstrated significantly poorer accuracy especially with respect to proper characterization of the pressure region ($1≤y≤2$). The k–ε model does produce a very similar distribution with the nonlinear eddy viscosity model for $u3′u3′¯+$ but, in general, a clear disadvantage for the linear eddy viscosity model is observed at high rotation number compared to the nonlinear model.

All of the turbulence closure approximations considered in this work invoke the nonlinear eddy viscosity hypothesis which models the Reynolds stresses as functions of the turbulent kinetic energy (k), dissipation rate (ε), and mean flow gradients in the form of the mean strain-rate ($Sij$) and rotation tensors ($Wij$). The turbulence model proposed by Reif et al. [5] (hereafter referred to as the PRDO model) is similar to the k$ε$ model in which partial differential equations (PDEs) for k and ε are solved. In addition, the PRDO model incorporates two PDEs for the temporal evolution of $v2¯$ and f, where $v2¯$ is the modeled (time-averaged) wall-normal Reynolds stress, and f is an elliptic relaxation function, which provides kinematic blocking of turbulent transport near the channel walls. Through analytical derivation for this set of four PDEs, the PRDO model produces the explicit expression for the Reynolds stresses shown in Eq. (5)Display Formula

(5)$ui′uj′¯=23kδij−2Cμ∗v2¯TSij$

where $Cμ∗$ is a function of the strain-rate ($Sij$) and rotation tensors ($Wij$) and Display Formula

(6)$T=max[kε;6(νε)1/2]$

represents a turbulent timescale.

The Speziale–Gatski (SG) and Girimaji (GI) turbulence models, proposed by Speziale and Gatski [6] and Girimaji [7], respectively, are EARSM which belong to a common subtype of nonlinear eddy viscosity models. In EARSM, an implicit algebraic model for the anisotropic tensor ($bij$) is constructed from the Reynolds stress transport equation through the modeling of individual energy balance terms and an explicit form for $bij$ is derived analytically. An expression for the Reynolds stresses is then obtained using the following relation: Display Formula

(7)$bij=ui′uj′¯2k−13δij$

The explicit expressions for the Reynolds stresses from the SG and GI models are shown in Eqs. (8) and (9), respectively. In both models, the terms η1 and η2 are equivalent to invariants of the normalized strain-rate and rotation tensors, or $SijSij$ and $WijWij$, respectively. In Eq. (8), α1 is a function of the production–dissipation ratio and in Eq. (9), the G coefficients are functions of η1 and η2Display Formula

(8)$ui′uj′¯=23kδij−2(SikSkj−13SmnSmnδij)−6(1+η12)α1k3+η12+6η22η12+6η22(Sij+SijWkj−WikSkj)$
Display Formula
(9)$ui′uj′¯=2kG1Sij+2kG2(SikWkj+WikSkj)+2kG3(SikSkj−13SmnSmnδij)+23kδij$

The model proposed by Grundestam et al. [8] (hereafter referred to as the Grundestam–Wallin–Johansson (GWJ) model) is also an EARSM derived from the Reynolds stress transport equation. Unlike the SG and GI models which use the linear SSG pressure–strain function to model the pressure–strain energy balance term [22], the GWJ model uses a nonlinear pressure–strain function in its formulation. The GWJ model proposes the explicit expression for the Reynolds stresses Display Formula

(10)$ui′uj′¯=23kδij+k∑k=110βkTijk$

where the $Tij$ tensor components are various functions of the mean strain-rate and rotation tensors. Similarly, the β coefficients depend on the tensor invariants defined as Display Formula

(11)$ΠS=SijSji, ΠΩ∗=WijWji, IV=SijWjkWki$

Most of these closure formulations have been shown to produce significant errors in predicting the mean velocity and temperature distributions for rotating turbulent channel flow (e.g., see, Ref. [8] for the GWJ model). Hence, in this work, we examine the details of these models to assess the likely sources of such errors in the modeled mean and fluctuating quantities.

###### RANS Models: Results.

The normal Reynolds stress $ui′ui′¯+$ and primary shear Reynolds stress $u1′u2′¯+$ distributions calculated using the four aforementioned RANS-based models are compared directly to the present DNS results for simulation cases A–D. For the no-rotation case A, the comparisons are shown in Fig. 9. The approximations of the GWJ and GI models are observed to be similar for all the Reynolds stress components. The SG and PRDO models demonstrate the best agreement with DNS results for the normal and shear Reynolds stress distributions, respectively. Also noted in Ref. [23] is that the SG model is observed to not satisfy realizability as the modeled distribution of $u2′u2′¯+$ contains nonphysical negative amplitudes near both channel walls. Realizability is a set of mathematical and physical principles that need to be satisfied to prevent the turbulence model from generating nonphysical results. The three realizability conditions proposed by Schumann [24] are

Display Formula

(12)$ui′ui′¯≥0$
Display Formula
(13)$ui′uj′¯2≤ui′ui′¯uj′uj′¯$
Display Formula
(14)$det(ui′uj′¯)≥0$

where det is the tensor determinant.

For case B ($Rob=0.2$), the modeled Reynolds stresses and DNS results are shown in Fig. 10 and all the four turbulence models demonstrate proper characterization of the suppressed and enhanced turbulence levels in the suction and pressure regions, respectively. However, the three EARSM-type models (GWJ, SG, and GI) display better agreement with DNS results than the PRDO model, which displayed identical distributions of $ui′ui′¯+$. The poor performance of the PRDO model is likely the result of the linear relationship between the Reynolds stress and mean strain-rate tensors (Eq. (5)). With respect to the DNS, the SG model showed the best agreement for the normal Reynolds stresses but all the four turbulence models failed to predict the Reynolds shear stress profile in the pressure region. For all the cases, the PRDO and GWJ models produce identical distributions for spanwise Reynolds stress $u3′u3′¯+$, which was unexpected due to the significant differences in the modeled PRDO and GWJ expressions shown in Eqs. (5) and (10), respectively. The SG model again fails to satisfy realizability as negative amplitudes of $u2′u2′¯+$ are observed near the pressure wall (y = 2) in Fig. 10(b).

The DNS and modeled Reynolds stresses for cases C ($Rob=0.5$) and D ($Rob=0.9$) are shown in Fig. 11 and 12, respectively. The invariance of the normal Reynolds stress distributions from the PRDO model continues to result in poor agreement with the DNS in relation to EARSM. The GWJ model demonstrates the best correspondence with DNS data for cases C and D. For case C, the modeled SG Reynolds stress distributions develop disparities with the DNS near the pressure wall as the peak modeled amplitude of $u1′u2′¯+$ is significantly overestimated in Fig. 11(d) and an unexpected decrease of $u3′u3′¯+$ amplitudes is observed in Fig. 11(c). For case D, these issues worsen and the modeled SG distributions demonstrate significant deviations from the DNS for all the Reynolds stress components except $u2′u2′¯+$. For cases C and D, the SG model is shown to satisfy the previously violated realizability condition as the modeled distributions of $u2′u2′¯+$ in Figs. 11(b) and 12(b) remain non-negative throughout the entire channel.

For no-rotation, all the four tested turbulence models produced Reynolds stress distributions in close agreement with the present DNS results. For rotating simulation cases B–D, the EARSM (GWJ, SG, and GI) had significantly better agreement with the DNS than the PRDO model, a different type of nonlinear eddy viscosity model. This observation supports the current trend of turbulence model design for rotational favoring EARSM-type closures [8]. One notable deficiency of EARSM is generally manifested in the near-wall regions where significant deviations of $ui′uj′¯$ from the DNS distributions are observed. It is therefore instructive to examine the parameterization used for the pressure–strain term as this term is responsible for intercomponent energy transfer.

###### EARSM: Pressure–Strain Modeling.

The pressure–strain term is of paramount importance as the correlations between pressure and rate-of-strain fluctuations play a dominant role in intercomponent energy transfer [21] within the Reynolds stress equation Display Formula

(15)$Dui′uj′¯Dt=Pij+Πij−εij+Cij+DijT$

The terms on the right-hand side of Eq. (15) represent, respectively, the production ($Pij$), pressure–strain ($Πij$), dissipation ($εij$), Coriolis ($Cij$), and diffusion terms ($DijT$).

Both the SG and GI EARSM apply the linear Speziale, Sarkar, and Gatski (SSG) pressure strain function (Eq. (16)), albeit with different sets of numerical coefficients as their algebraic formulations differ Display Formula

(16)$Πij=−C1εbij+C2kSij+C3k(bikSkj+Sikbkj)−23C3bmnSmnδij−C4k(bikWkj−Wikbkj)$

In Figs. 13 and 14, the distributions of the pressure–strain closure formula for the SG and GI models are compared to the pressure–strain distributions directly extracted from the present DNS data to evaluate possible sources of error specifically for cases A ($Rob=0$) and B ($Rob=0.2$).

Figure 13 displays results for case A ($Rob=0$), the modeled pressure–strain function values for both models are similar to the DNS results in the interior of the channel. However, near both walls where most of the intercomponent energy transfer takes place, the modeled pressure–strain amplitudes are drastically overestimated for both the SG and GI models. This discrepancy is not unexpected as the SSG function is a single-point closure which captures only the local effects of the pressure–strain correlations, neglecting nonlocal effects such as wavevector information [25]. In Fig. 14, similar comparisons are presented for case B ($Rob=0.2$). It is observed that the modeled pressure–strain distributions have better agreement with the DNS data with increasing rotation number near the suction wall (y = 0) and the only major disparity occurs near the pressure wall with suppressed amplitudes from relaminarization in the suction region.

In Figs. 13 and 14, it is observed that $Π11$ is a loss term for $u1′u1′¯+$ and $Π22$ and $Π33$ are gain terms for their corresponding Reynolds normal stresses. This is expected because the main contribution of the pressure–strain correlation term is to isotropize (return-to-isotropy) the normal Reynolds stress components. In nonrotating turbulent channel flow, the production term contributes only to $u1′u1′¯+$, so the pressure–strain tensor isotropizes the normal Reynolds stresses through a negative contribution to $u1′u1′¯$ and positive contributions to $u2′u2′¯+$ and $u3′u3′¯+$. It is also observed that these pressure–strain models extract too much energy from $u1′u1′¯+$ and transfer excess energy to $u2′u2′¯+$ and $u3′u3′¯+$, which is directly reflected in the corresponding modeled Reynolds stress distributions in Figs. 9 and 10.

## Heat Transfer Models

###### Heat Transfer Models: Overview.

Heat transfer models parameterize turbulent heat fluxes which are correlations between velocity and temperature fluctuations. The following two heat flux models that we consider in this work are used frequently in engineering calculations [26]. The YWL heat flux model proposed by Younis et al. [9] is an EAHFM. For EAHFM, the transport equation for $ui′θ′¯$ may be written as Display Formula

(17)$Dui′θ′¯Dt=Diθ+Piθ+Giθ+Πiθ+εiθ$

The terms $Diθ, Piθ, Giθ, Πiθ$, and $εiθ$ represent heat flux diffusion, production of heat flux through shear forces, production of heat flux through rotational forces, pressure–temperature gradient correlations, and viscous dissipation, respectively.

The $ui′θ′¯$ expression obtained from the YWL model is given in Eq. (18). Here, the coefficients $C2t, C3t$, and $C4t$ are absolute constants while $C1t$ is a wall function dependent on invariants of the anisotropic stress tensor and Ω is the angular rotation rate vector Display Formula

(18)$−ui′θ′¯=C1tk2ε∂θ¯∂xi+∂θ¯∂xj[C2tkεui′uj′¯+C3tk3ε2(∂Ui¯∂xj+εmjiΩm)]+C4tk2ε2∂θ¯∂xj[ui′uk′¯(∂Uj¯∂xk+εjmkΩm)+uj′uk′¯(∂Ui¯∂xk+εimkΩm)]$

The heat flux model proposed by Abe and Suga [10] (hereafter referred to as the SA model) is also an EAHFM but unlike the YWL model, the SA model incorporates nonlinear Reynolds stress terms in its model expressions given in Eq. (19). The terms $Cθ$ and τ represent a wall function and time scale ($k/ε$), respectively. The expressions for tensors $αij$ and $σij$ are shown in Eqs. (20) and (21), respectively. $Ωij$ is the mean rotation tensor and the coefficients $cα0, cα1, cσ0, cσ1$, and $cσ2$ are composed of weight functions and the normalized second invariant of the mean strain-rate tensor Display Formula

(19)$ui′θ′¯=−Cθkτ(σij+αij)∂θ¯∂xj$
Display Formula
(20)$αij=cα0τΩij+cα1τ(Ωilul′uj′¯/k+Ωilul′uj′¯/k)$
Display Formula
(21)$σij=cσ0δij+cσ1ui′uj′¯/k+cσ2ui′ul′¯ul′uj′¯/k2$

###### Heat Transfer Models: Results.

Similar to previous RANS model testing, the YWL and SA modeled turbulent heat flux distributions are directly compared to DNS results for simulation cases A–D. In Figs. 15(a) and 15(b), the DNS and model distributions of the streamwise and wall-normal turbulent heat fluxes, $u1′θ′¯+$ and $u2′θ′¯+$, respectively, are shown for the no-rotation case A ($Rob=0$). Good agreement is displayed with the DNS results for both $u1′θ′¯+$ and $u2′θ′¯+$. In Figs. 15(c) and 15(d), inconsistencies develop in the YWL and SA modeled heat flux distributions for case B ($Rob=0.2$). Both models show good agreement with the DNS for $u1′θ′¯+$ in the suction region but overestimate the peak amplitudes near the pressure wall (y = 2). In addition, neither model approximates the linear distribution of $u2′θ′¯+$ at the channel center shown in the DNS results.

In Fig. 16, the modeled heat flux distributions from the YWL and SA models are shown and compared with the DNS results for cases C ($Rob=0.5$) and D ($Rob=0.9$). In Fig. 16(a), a nonphysical drop-off in $u1′θ′¯+$ distributions from both the SA and YWL models is observed for case C near y = 0.6. In Fig. 16(c), this discrepancy worsens at a higher rotation number in case D ($Rob=0.9$). Increasing departure from DNS data was also observed for increasing rotation number in the modeled distributions of $u2′θ′¯+$. With an increase in Rossby number, the DNS amplitudes of $u2′θ′¯+$ are shown to become more suppressed in the suction region but the modeled distributions instead display significant amplitude increases of $u2′θ′¯+$ throughout the channel. Through comparison of turbulent heat flux distributions with DNS results for system rotation, both heat transfer models displayed significant deviations that increased with increasing rotation number. As pressure–strain correlations were shown to influence EARSM accuracy, it is instructive to examine the parameterization used for the pressure–temperature gradient correlations in EAHFM.

## EAHFM: Pressure–Temperature Gradient Correlation Modeling

In this work, EAHFM modeling of the pressure–temperature gradient correlation term Display Formula

(22)$Πiθ=p′∂θ′/∂xi¯$

in the transport equation for the turbulent heat fluxes (Eq. (17)) is investigated for error contributions to the corresponding modeled profiles of $ui′θ′¯$. In Sec. 3.3, correlations of pressure fluctuations were demonstrated to produce significant contributions to turbulence dynamics and Liu and Lun [19] noted that $Pi2θ$ was a major source term to $u2′θ′¯$. The SA model uses the pressure–scalar correlation model proposed in Ref. [27] Display Formula

(23)$Πiθ=−C1θεkui′θ¯−C2θuk′θ¯∂Ui¯∂xk$

where ($C1θ, C2θ$) = (3, 0.4). The YWL model uses the same pressure–scalar correlation model with additional corrective terms Display Formula

(24)$Πiθ=−C1θεkui′θ′¯−C2θuk′θ′¯∂Ui¯∂xk−C1θwPiθ−C2θwC2θPiθninkf$

where ($C1θ, C2θ, C1θw, C2θw$) = (3, 0.5, 0.5, 0.5), and $Piθ$ represents the production rate of heat flux through mean velocity and temperature gradients Display Formula

(25)$Piθ=−(uk′θ′¯∂Ui¯∂xk+uk′ui′¯∂θ¯∂xk)$

n and f denote a unit vector and wall-damping function, respectively.

In Fig. 17, the modeled pressure–temperature gradient distributions from the YWL and SA models are compared to DNS results for selected cases A ($Rob=0$) and B ($Rob=0.2$). For case A, $Π1θ$ and $Π2θ$ are gain terms for the corresponding turbulent heat fluxes, in agreement with Ref. [19]. In Figs. 17(a) and 17(b), both models demonstrate significant overestimations of $Π1θ$ in the near-wall regions, similar to the tested pressure–strain models. There is significantly better agreement between the modeled and DNS distributions for $Π2θ$, and the SA model shows exceptional accuracy.

For case B, $Π1θ$ continues to be a source term for $u1′θ′¯+$ but $Π2θ$ changes to become a loss term for $u2′θ′¯+$ in the suction region. This feature is not captured by the YWL and SA models, resulting in significant model disparity from DNS results in the suction region. This major error in the pressure–temperature gradient correlation term is shown to contribute to inaccuracies in the corresponding modeled heat flux distributions in Fig. 15(d). $Π2θ$ is incorrectly characterized by both models as a gain term for case B, resulting in significant overestimations of modeled $u2′θ′¯+$ in the suction region. The examination of error contributions from $Πij$ and $Πiθ$ functions to their corresponding models demonstrated the profound effects of pressure correlations on turbulent transport in nonrotating and rotating turbulent channel flow.

## Conclusions

For spanwise-rotating turbulent channel flow, the Reynolds stress distributions produced from a linear and nonlinear eddy viscosity model were compared to demonstrate improved accuracy for the nonlinear model over the linear model. In a comparison of four nonlinear eddy viscosity models with DNS data, EARSM were the most compatible with DNS results in modeling the Reynolds stresses for turbulent channel flow subject to spanwise rotation. The Speziale–Gatski (SG) model was shown to be the most compatible for zero and low rotation numbers but displayed significant deviations near the pressure wall at high rotation numbers. The GWJ model showed the best agreement with the DNS data at high rotation numbers. Heat flux models were in good agreement with DNS data for the no-rotation case but with system rotation, the models deviated from the DNS, increasing at higher rotation rates.

The pressure–strain models of two EARSM (Girimaji and SG) were shown to have significant disagreements with the DNS data in the near-wall regions and the pressure–temperature gradient models of two EAHFM (YWL and SA) demonstrated inaccurate characterization of the suction region with system rotation. The errors in the modeled contributions from these terms resulted in degeneration of the predictive capabilities of their respective closure models. These errors contributed to inaccurate Reynolds stress amplitudes in the near-wall regions and an inaccurate modeled distribution shape for wall-normal turbulent heat flux in EARSM and EAHFM, respectively. The present results indicate that correct characterization of pressure fluctuations is a crucial factor in both EARSM and EAHFM design for spanwise-rotating turbulent channel flow.

## Acknowledgements

This work was supported in part by NSF Grant IDR No. 1131802 and by the Air Force Office of Scientific Research under Grant No. FA9550-15-0495. The authors acknowledge the NSF Extreme Science and Engineering Discovery Environment (XSEDE) and Texas Advanced Computing Center (TACC) at the University of Texas at Austin for providing the computational resources that have contributed to the results reported within this paper under Grant No. TG-CTS120002. The authors thank Cesar Galan funded through the University of Colorado Undergraduate Research Opportunity Program for his contributions to this work.

## Nomenclature

• bij =

anisotropic tensor, $ui′uj′¯/2k−13δij$

• Cij =

Coriolis

• $Diθ$ =

diffusion (heat flux transport)

• $DijT$ =

diffusion budget term (heat flux transport)

• $Giθ$ =

production from rotational forces (heat flux transport)

• k =

turbulent kinetic energy, $12ui′ui′¯$

• Nu =

Nusselt number, $qwδ/κΔT$

• p =

nondimensional effective pressure, $p0−(1/8)Roc2r2$

• Pij =

production

• $Piθ$ =

production from shear forces (heat flux transport)

• p =

quantities on the pressure wall

• p0 =

static pressure

• Pr =

Prandtl number, $ν/κ$

• qw =

mean heat flux at channel walls, $κ∂θ¯/∂y|wall$

• r =

nondimensional distance from axis of rotation

• Rec =

centerline Reynolds number, $ucδ/ν$

• $Reτ$ =

friction Reynolds number, $uτδ/ν$

• Rob =

bulk Rossby number, $2Ωδ/ub$

• Roc =

centerline Rossby number, $2Ωδ/uc$

• s =

quantities on the suction wall

• Sij =

mean strain-rate tensor, $12(∂Ui/∂xj+∂Uj/∂xi)$

• t =

nondimensional time

• $Tτ$ =

friction temperature

• ub =

mean bulk velocity

• uc =

mean laminar centerline velocity

• $uτ$ =

friction velocity

• Ui =

mean velocity in xi direction [U, V, W]

• $ui′$ =

fluctuating velocity in xi direction [$u′, v′, w′$]

• $ui′ui′¯$ =

Reynolds stress

• $ui′θi′¯$ =

turbulent heat flux

• ν =

kinematic viscosity

• Wij =

mean rotation tensor, $12(∂Ui/∂xj−∂Uj/∂xi)+εijkΩk$

• xi =

Cartesian coordinate [x, y, z]

• $ΔT$ =

temperature difference between upper/lower walls

• δ =

channel half-height

• δij =

Kronecker delta

• ε =

turbulent dissipation rate

• $εijk$ =

cyclic permutation tensor

• εij =

dissipation

• $εiθ$ =

viscous dissipation (heat flux transport)

• θ =

nondimensional temperature

• κ =

thermal diffusivity

• τ =

turbulent time scale, $k/ε$

• Πij =

pressure-strain

• $Πiθ$ =

pressure–temperature gradient correlations (heat flux transport)

• Ωi =

angular rotation vector

## References

Grundestam, O. , Wallin, S. , and Johansson, A. , 2008, “ Direct Numerical Simulations of Rotating Turbulent Channel Flow,” J. Fluid Mech., 598, pp. 177–199.
Wu, H. , and Kasagi, N. , 2004, “ Effects of Arbitrary Directional System Rotation on Turbulent Channel Flow,” Phys. Fluids, 16(4), pp. 979–990.
Acharya, S. , Sethuraman, E. , and Nikotopoulos, S. , 2012, “ Mass and Heat Transfer in Rotating, Smooth, High Aspect Ratio Coolant Channels With Curved Walls,” ASME J. Turbomach., 131(2), p. 021002.
Lakehal, D. , 2002, “ Near-Wall Modeling of Turbulent Convective Heat Transport in Film Cooling of Turbine Blades With the Aid of Direct Numerical Simulation Data,” ASME J. Turbomach., 124(3), pp. 485–498.
Reif, B. P. , Durbin, P. , and Ooi, A. , 1999, “ Modeling Rotational Effects in Eddy-Viscosity Closure,” Int. J. Heat Fluid Flow, 20(6), pp. 563–573.
Speziale, C. , and Gatski, T. , 1993, “ On Explicit Algebraic Stress Models for Complex Turbulent Flows,” J. Fluid Mech., 254, pp. 59–78.
Girimaji, S. , 1996, “ Fully Explicit and Self-Consistent Reynolds Stress Model,” Theor. Comput. Fluid Dyn., 8(6), pp. 387–402.
Grundestam, O. , Wallin, S. , and Johansson, A. , 2005, “ An Explicit Algebraic Reynolds Stress Model Based on a Nonlinear Pressure Strain Rate Model,” Int. J. Heat Fluid Flow, 26(5), pp. 732–745.
Younis, B. , Weigand, B. , and Laqua, A. , 2012, “ Prediction of Turbulent Heat Transfer in Rotating and Nonrotating Channels With Wall Suction and Blowing,” ASME J. Heat Transfer, 134(7), p. 071702.
Abe, K. , and Suga, K. , 2000, “ Towards the Development of a Reynolds-Averaged Algebraic Turbulent Scalar-Flux Model,” Int. J. Heat Fluid Flow, 22(1), pp. 19–29.
Waggy, S. , Biringen, S. , and Sullivan, P. , 2013, “ Direct Numerical Simulation of Top-Down and Bottom-Up Diffusion in the Convective Boundary Layer,” J. Fluid Mech., 724, pp. 581–606.
Waggy, S. , Kucala, A. , and Biringen, S. , 2014, “ Parallel Implementation of a Navier–Stokes Solver: Turbulent Ekman Layer Direct Simulation,” Int. J. Comput. Methods, 11(5), pp. 443–460.
Waggy, S. , Marlatt, S. , and Biringen, S. , 2011, “ Direct Numerical Simulation of the Turbulent Ekman Layer: Instantaneous Flow Structures,” J. Thermophys. Heat Transfer, 25(2), pp. 309–318.
Waggy, S. , Biringen, S. , and Kucala, A. , 2015, “ Wake Effects on Turbulent Transport in the Convective Boundary Layer,” Geophys. Astrophys. Fluid Dyn., 109(5), pp. 465–479.
Kristofferson, R. , and Andersson, H. , 1993, “ Direct Simulations of Low Reynolds Number Turbulent Flow in a Rotating Channel,” J. Fluid Mech., 253, pp. 163–197.
Kim, J. , Moin, P. , and Moser, R. , 1987, “ Turbulence Statistics in Fully Developed Channel Flow at Low Reynolds Number,” J. Fluid Mech., 177, pp. 133–166.
Kreplin, H. , and Eckelmann, M. , 1979, “ Behavior of the Three Fluctuating Velocity Components in the Wall Region of a Turbulent Channel Flow,” Phys. Fluids, 22(7), pp. 1233–1239.
Kasagi, N. , and Nishino, K. , 1991, “ Probing Turbulence With Three-Dimensional Particle-Tracking Velocimetry,” Exp. Therm. Fluid Sci., 4(5), pp. 601–612.
Liu, N. , and Lun, X. , 2007, “ Direct Numerical Simulation of Spanwise Rotating Turbulent Channel Flow With Heat Transfer,” Int. J. Numer. Methods Fluids, 53(11), pp. 1689–1706.
Launder, B. , and Shama, B. , 1974, “ Application of the Energy-Dissipation Model of Turbulence to the Calculation of Flow Near a Spinning Disc,” Lett. Heat Mass Transfer, 1(2), pp. 131–138.
Launder, B. , Reece, G. , and Rodi, W. , 1975, “ Progress in the Development of a Reynolds Stress Turbulent Closure,” J. Fluid Mech., 68(3), pp. 537–566.
Speziale, C. , Sarkar, S. , and Gatski, T. , 1991, “ Modeling the Pressure-Strain Correlation of Turbulence: An Invariant Dynamical System Approach,” J. Fluid Mech., 227, pp. 245–272.
Rung, T. , and Thiele, F. , 1999, “ On the Realizability of Nonlinear Stress-Strain Relationships for Reynolds Stress Closures,” Flow Turbul. Combust., 60(4), pp. 333–359.
Schumann, U. , 1977, “ Realizability of Reynolds-Stress Turbulence Models,” Phys. Fluids, 20(5), pp. 721–725.
Mishra, A. , and Girimaji, S. , 2013, “ Intercomponent Energy Transfer in Incompressible Homogeneous Turbulence: Multi-Point Physics and Amenability to One-Point Closures,” J. Fluid Mech., 731, pp. 639–681.
Li, X. , Ren, J. , and Jiang, H. , 2014, “ Film Cooling Modeling of Turbine Blades Using Algebraic Anisotropic Turbulence Models,” ASME J. Turbomach., 136(11), p. 111006.
Launder, B. , 1975, “ On the Effects of a Gravitational Field on the Turbulent Transport of Heat and Momentum,” J. Fluid Mech., 67(03), pp. 569–581.
View article in PDF format.

## References

Grundestam, O. , Wallin, S. , and Johansson, A. , 2008, “ Direct Numerical Simulations of Rotating Turbulent Channel Flow,” J. Fluid Mech., 598, pp. 177–199.
Wu, H. , and Kasagi, N. , 2004, “ Effects of Arbitrary Directional System Rotation on Turbulent Channel Flow,” Phys. Fluids, 16(4), pp. 979–990.
Acharya, S. , Sethuraman, E. , and Nikotopoulos, S. , 2012, “ Mass and Heat Transfer in Rotating, Smooth, High Aspect Ratio Coolant Channels With Curved Walls,” ASME J. Turbomach., 131(2), p. 021002.
Lakehal, D. , 2002, “ Near-Wall Modeling of Turbulent Convective Heat Transport in Film Cooling of Turbine Blades With the Aid of Direct Numerical Simulation Data,” ASME J. Turbomach., 124(3), pp. 485–498.
Reif, B. P. , Durbin, P. , and Ooi, A. , 1999, “ Modeling Rotational Effects in Eddy-Viscosity Closure,” Int. J. Heat Fluid Flow, 20(6), pp. 563–573.
Speziale, C. , and Gatski, T. , 1993, “ On Explicit Algebraic Stress Models for Complex Turbulent Flows,” J. Fluid Mech., 254, pp. 59–78.
Girimaji, S. , 1996, “ Fully Explicit and Self-Consistent Reynolds Stress Model,” Theor. Comput. Fluid Dyn., 8(6), pp. 387–402.
Grundestam, O. , Wallin, S. , and Johansson, A. , 2005, “ An Explicit Algebraic Reynolds Stress Model Based on a Nonlinear Pressure Strain Rate Model,” Int. J. Heat Fluid Flow, 26(5), pp. 732–745.
Younis, B. , Weigand, B. , and Laqua, A. , 2012, “ Prediction of Turbulent Heat Transfer in Rotating and Nonrotating Channels With Wall Suction and Blowing,” ASME J. Heat Transfer, 134(7), p. 071702.
Abe, K. , and Suga, K. , 2000, “ Towards the Development of a Reynolds-Averaged Algebraic Turbulent Scalar-Flux Model,” Int. J. Heat Fluid Flow, 22(1), pp. 19–29.
Waggy, S. , Biringen, S. , and Sullivan, P. , 2013, “ Direct Numerical Simulation of Top-Down and Bottom-Up Diffusion in the Convective Boundary Layer,” J. Fluid Mech., 724, pp. 581–606.
Waggy, S. , Kucala, A. , and Biringen, S. , 2014, “ Parallel Implementation of a Navier–Stokes Solver: Turbulent Ekman Layer Direct Simulation,” Int. J. Comput. Methods, 11(5), pp. 443–460.
Waggy, S. , Marlatt, S. , and Biringen, S. , 2011, “ Direct Numerical Simulation of the Turbulent Ekman Layer: Instantaneous Flow Structures,” J. Thermophys. Heat Transfer, 25(2), pp. 309–318.
Waggy, S. , Biringen, S. , and Kucala, A. , 2015, “ Wake Effects on Turbulent Transport in the Convective Boundary Layer,” Geophys. Astrophys. Fluid Dyn., 109(5), pp. 465–479.
Kristofferson, R. , and Andersson, H. , 1993, “ Direct Simulations of Low Reynolds Number Turbulent Flow in a Rotating Channel,” J. Fluid Mech., 253, pp. 163–197.
Kim, J. , Moin, P. , and Moser, R. , 1987, “ Turbulence Statistics in Fully Developed Channel Flow at Low Reynolds Number,” J. Fluid Mech., 177, pp. 133–166.
Kreplin, H. , and Eckelmann, M. , 1979, “ Behavior of the Three Fluctuating Velocity Components in the Wall Region of a Turbulent Channel Flow,” Phys. Fluids, 22(7), pp. 1233–1239.
Kasagi, N. , and Nishino, K. , 1991, “ Probing Turbulence With Three-Dimensional Particle-Tracking Velocimetry,” Exp. Therm. Fluid Sci., 4(5), pp. 601–612.
Liu, N. , and Lun, X. , 2007, “ Direct Numerical Simulation of Spanwise Rotating Turbulent Channel Flow With Heat Transfer,” Int. J. Numer. Methods Fluids, 53(11), pp. 1689–1706.
Launder, B. , and Shama, B. , 1974, “ Application of the Energy-Dissipation Model of Turbulence to the Calculation of Flow Near a Spinning Disc,” Lett. Heat Mass Transfer, 1(2), pp. 131–138.
Launder, B. , Reece, G. , and Rodi, W. , 1975, “ Progress in the Development of a Reynolds Stress Turbulent Closure,” J. Fluid Mech., 68(3), pp. 537–566.
Speziale, C. , Sarkar, S. , and Gatski, T. , 1991, “ Modeling the Pressure-Strain Correlation of Turbulence: An Invariant Dynamical System Approach,” J. Fluid Mech., 227, pp. 245–272.
Rung, T. , and Thiele, F. , 1999, “ On the Realizability of Nonlinear Stress-Strain Relationships for Reynolds Stress Closures,” Flow Turbul. Combust., 60(4), pp. 333–359.
Schumann, U. , 1977, “ Realizability of Reynolds-Stress Turbulence Models,” Phys. Fluids, 20(5), pp. 721–725.
Mishra, A. , and Girimaji, S. , 2013, “ Intercomponent Energy Transfer in Incompressible Homogeneous Turbulence: Multi-Point Physics and Amenability to One-Point Closures,” J. Fluid Mech., 731, pp. 639–681.
Li, X. , Ren, J. , and Jiang, H. , 2014, “ Film Cooling Modeling of Turbine Blades Using Algebraic Anisotropic Turbulence Models,” ASME J. Turbomach., 136(11), p. 111006.
Launder, B. , 1975, “ On the Effects of a Gravitational Field on the Turbulent Transport of Heat and Momentum,” J. Fluid Mech., 67(03), pp. 569–581.

## Figures

Fig. 3

Verification of DNS data for case A (Rob=0): (a) (u1′u1′¯+)1/2, (b) (u2′u2′¯+)1/2, (c) (u3′u3′¯+)1/2, and (d) u1′u2′¯+. Black: DNS; red: Kasagi and Nishino [18]; and square: Kreplin and Eckelmann [17].

Fig. 2

One-dimensional energy spectra along the streamwise axis for case A (Rob=0). Solid: Euu; dash-dot: Evv; dotted: Eww; and dashed: Kolmogorov −5/3 spectrum.

Fig. 1

Geometry of DNS computational domain for turbulent channel flow with rotation in the spanwise direction

Fig. 4

Mean velocity distributions with 2Ω lines for full simulation cases A–D. Black: case A (Rob=0); red: case B (Rob=0.2); blue: case C (Rob=0.5); and magenta: case D (Rob=0.9). −−−: 2Ω lines.

Fig. 6

Instantaneous ωx′ map for a y–z section at x=2π: (a) case A (Rob=0), (b) case B (Rob=0.2), and (c) case C (Rob=0.5)

Fig. 5

Mean temperature distributions for full simulation cases A–D. Black: case A (Rob=0); red: case B (Rob=0.2); blue: case C (Rob=0.5); and magenta: case D (Rob=0.9).

Fig. 7

Time-averaged v and w velocity vectors for a y–z section at x=2π: (a) case A (Rob=0), (b) case B (Rob=0.2), and (c) case C (Rob=0.5)

Fig. 10

Modeled Reynolds stress profiles for case B (Rob=0.2): (a) u1′u1′¯+, (b) u2′u2′¯+, (c) u3′u3′¯+, and (d) u1′u2′¯+. Black: DNS; blue: PRDO; green: SG; red: GI; and magenta: GWJ.

Fig. 11

Modeled Reynolds stress profiles for case C (Rob=0.5): (a) u1′u1′¯+, (b) u2′u2′¯+, (c) u3′u3′¯+, and (d) u1′u2′¯+. Black: DNS; blue: PRDO; green: SG; red: GI; and magenta: GWJ.

Fig. 12

Modeled Reynolds stress profiles for case D (Rob=0.9): (a) u1′u1′¯+, (b) u2′u2′¯+, (c) u3′u3′¯+, and (d) u1′u2′¯+. Black: DNS; blue: PRDO; green: SG; red: GI; and magenta: GWJ.

Fig. 8

Modeled Reynolds stress profiles for case D (Rob=0.9): (a) u1′u1′¯+, (b) u2′u2′¯+, (c) u3′u3′¯+, and (d) u1′u2′¯+. Black: DNS; red: linear eddy viscosity model [20]; and blue: nonlinear eddy viscosity model [8].

Fig. 9

Modeled Reynolds stress profiles for case A (Rob=0): (a) u1′u1′¯+, (b) u2′u2′¯+, (c) u3′u3′¯+, and (d) u1′u2′¯+. Black: DNS; blue: PRDO; green: SG; red: GI; and magenta: GWJ.

Fig. 13

Modeled pressure–strain budget term profiles for case A (Rob=0): (a) Π11, (b) Π22, (c) Π33, and (d) Π12. Black: DNS; red: SG; and blue: GI.

Fig. 14

Modeled pressure–strain budget term profiles for case B (Rob=0.2): (a) Π11, (b) Π22, (c) Π33, and (d) Π12. Black: DNS; red: SG; and blue: GI.

Fig. 15

Modeled turbulent heat flux profiles for cases A (Rob=0) and B (Rob=0.2): (a) u1′θ′¯+ (case A), (b) u2′θ′¯+ (case A), (c) u1′θ′¯+ (case B), and (d) u2′θ′¯+ (case B). Black: DNS; red: YWL; and blue: SA.

Fig. 16

Modeled turbulent heat flux profiles for cases C (Rob=0.5) and D (Rob=0.9): (a) u1′θ′¯+ (case C), (b) u2′θ′¯+ (case C), (c) u1′θ′¯+ (case D), and (d) u2′θ′¯+ (case D). Black: DNS; red: YWL; and blue: SA.

Fig. 17

Modeled pressure–temperature gradient profiles for cases A (Rob=0) and B (Rob=0.2): (a) Π1θ (case A), (b) Π2θ (case A), (c) Π1θ (case B), and (d) Π2θ (case B). Black: DNS; red: YWL; and blue: SA.

## Tables

Table 1 Case descriptions and initial conditions
Table 2 Friction Reynolds ($Reτ$) and Nusselt (Nu) numbers for present DNS cases A–D. Subscripts s and p denote the suction (y = 0) and pressure walls (y = 2), respectively.

## Errata

Some tools below are only available to our subscribers or users with an online account.

### Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related Proceedings Articles
Related eBook Content
Topic Collections