Passive scalar transport in rotating turbulent channel flow Geert Brethouwer 1 1 Linné FLOW Centre, KTH Mechanics, SE-144 Stockholm, Sweden (Dated: October 31, 217) arxiv:171.1931v1 [phsics.flu-dn] 3 Oct 217 Abstract Passive scalar transport in turbulent channel flow subject to spanwise sstem rotation is studied b direct numerical simulations. The Renolds number Re = U b h/ν is fixed at 2 and the rotation number Ro = 2Ωh/U b is varied from to 1.2, where U b is the bulk mean velocit, h the half channel gap width and Ω the rotation rate. The scalar value is constant but different at the two walls, leading to stead scalar transport across the channel. The rotation causes an unstable channel side with relativel strong turbulence and turbulent scalar transport, and a stable channel side with relativel weak turbulence or laminar-like flow, weak turbulent scalar transport but large scalar fluctuations and steep mean scalar gradients. The distinct turbulent-laminar patterns observed at certain Ro on the stable channel side induce similar patterns in the scalar gradient field. The main conclusions of the stud are that rotation reduces the similarit between the scalar and velocit field and that the Renolds analog for scalar-momentum transport does not hold for rotating turbulent channel flow. This is shown b a reduced correlation between velocit and scalar fluctuations, and a strongl reduced turbulent Prandtl number of less than.2 on the unstable channel side awa from the wall at higher Ro. On the unstable channel side, scalar scales become larger than turbulence scales according to spectra and the turbulent scalar flux vector becomes more aligned with the mean scalar gradient owing to rotation. Budgets in the governing equations of the scalar energ and scalar fluxes are presented and discussed as well as other statistics relevant for turbulence modelling. 1
1. INTRODUCTION In this paper, I present a numerical stud of passive scalar transport in turbulent channel flow subject to spanwise sstem rotation. A passive scalar, a scalar that does not influence the flow, can represent e.g. a contaminant or small temperature variation. The present investigation is therefore connected to heat and mass transfer in rotating flows. Examples of these can be found in man industrial apparatus such as turbo machiner, separators and chemical reactors. Brethouwer (25) and Kassinos et al. (27) investigated passive scalar transport in homogeneous turbulent shear flow subject to rotation about an axis normal to the mean shear plane b direct numerical simulations (DNS). The found that rotation has a large impact on the scalar fluctuation intensit and direction of turbulent scalar transport. With a transverse mean scalar gradient, scalar transport is predominantl in the streamwise direction if the flow is non-rotating, but if the sstem rotation exactl counteracts the rotation b the mean shear the scalar transport is nearl aligned with the mean scalar gradient. Brethouwer (25) observed that rotation can strongl reduce the turbulent Prandtl number. Man of these rotation effects are also shown b rapid distortion theor (Brethouwer 25). A generic case to investigate the influence of sstem rotation on mass or heat transfer in wall flows is turbulent plane channel flow including a passive scalar subject to spanwise sstem rotation. The effect of spanwise rotation on turbulent channel flow has alread been thoroughl investigated numericall (see e.g. Kristoffersen & Andersson 1993, Grundestam et al. 28, Xia et al. 216) and experimentall (see e.g. Johnston et al. 1972). Most of the numerical studies are limited to Re τ = u τ h/ν 194, where u τ is the friction velocit, h the half channel gap width, ν the viscosit, but recentl I have extended the stud of spanwise rotating channel flow to Re = U b h/ν = 31 6 and a wide range of rotation numbers Ro = 2Ωh/U b (Brethouwer 217). Here, U b is the bulk mean velocit and Ω the sstem rotation rate. The main conclusions of these studies can be summarized as follows. Turbulence and especiall the wall-normal velocit fluctuations are augmented on the channel side were the sstem rotation is anti-cclonic, i.e. in the direction opposite to the rotation induced b the mean shear, while the are damped on the channel side were the rotation is cclonic at moderate Ro (Xia et al. 216). These channel sides are from now on called the unstable and stable sides, respectivel. Large streamwise roll cells, 2
sometimes called Talor-Görtler vortices, are observed on the unstable channel sides at low to moderate Ro (Liu & Lu 27, Dai et al. 216, Brethouwer 217). The flow on the stable side relaminarizes partl or full at moderate Ro whereas in the limit of ver high Ro the whole channel relaminarizes (Grundestam et al. 28, Wallin et al. 213, Brethouwer 217). A ke feature of spanwise rotating channel flow is the development of linear part in the mean streamwise velocit profile on the unstable side where the absolute mean vorticit is nearl zero (Grundestam et al. 28, Xia et al. 216). Another remarkable feature is the occurrence of a linear instabilit in a range of Re and Ro, which leads to recurring strong bursts of turbulence on the stable channel side (Brethouwer et al. 214, Brethouwer 216). Passive scalar transport in non-rotating turbulent channel flow has been investigated extensivel as well, see e.g. the DNS studies b Kawamura et al. (1998, 1999) and Johansson & Wikström (1999). The velocit and scalar fields show a quite high degree of similarit in several aspects, especiall near the wall (Abe & Antonia 29, Antonia et al. 29, Dharmarathne et al. 216). Pirozzoli et al. (216) carried out DNS of passive scalar transport in turbulent channel flow for Re τ up to 488 and Prandtl number P r = ν/α between.2 and 1, where α is the scalar diffusivit. The observed that the mean scalar profiles follow a logarithmic law and that large-scale structures are present in the scalar field in the outer laer. The turbulent Prandtl number was found to be close to.85 in a large part of the channel, suggesting that the Renolds analog for scalar-momentum transfer is valid in a non-rotating turbulent channel flow. Scalar transport in rotating channel flow has been much less examined. Matsubara & Alfredsson (1996) have experimentall investigated the momentum and heat transfer in a laminar rotating channel flow. The found that streamwise roll cells, present at higher rotation rates, have a profound effect on heat transfer, and concluded that the Renolds analog is not valid since the Nusselt number changes but not the skin friction under the influence of rotation. Nagano & Hattori (23) and Liu & Lu (27) carried out DNS of passive scalar transport in spanwise rotating turbulent channel flow at Re τ = 15 and 194, respectivel, for Ro τ = 2Ωh/u τ 7.5 corresponding to Ro.5. The observed a reduced turbulent scalar transport on the stable channel side, resulting in a relativel steep mean scalar gradient and strong scalar fluctuations. Liu & Lu (27) also observed that the streamwise turbulent scalar transport is reduced b rotation, like in rotating homogeneous turbulent shear flow (Brethouwer 25), and that streamwise roll cells significantl contribute to scalar transport, 3
FIG. 1. Flow configuration and coordinate sstem. in particular at low rotation rates. The Nusselt number showed a moderate decrease with rotation rate in their DNS. Wu & Kasagi (24) studied scalar transport in turbulent channel flow at Re = 228 subject to sstem rotation with varing directions b DNS. Also in their DNS, rotation has a large impact on the direction and rate of turbulent scalar transport. Yang et al. (211) observed that the velocit and scalar field are quite strongl correlated on the unstable side but much less so on the stable side in DNS of spanwise rotating channel flow at Re = 2666. The objective of this stud is to investigate passive scalar transport in plane turbulent channel flow subject to spanwise sstem rotation b DNS at higher Re and for a wider range of Ro than in previous studies in order to obtain a better understanding of turbulent heat and mass transfer in rapidl rotating wall-bounded flows. I will show that rotation weakens the similarit between the scalar and velocit field and that the Renolds analog is not necessaril valid in rotating channel flow. With m stud, I aim to assist the development of models for transport and mixing in rotating turbulent flows. Several models for turbulent scalar and heat transfer in rotating wall flows have been proposed (Nagano & Hattori 23, Hattori et al. 29, Müller et al. 215, Hsieh et al. 216) and data on scalar transport in rotating channel flow ma be used to validate them. 2. NUMERICAL PROCEDURE I have carried out DNS of plane turbulent channel flow with a passive scalar subject to rotation about the spanwise direction. The flow geometr and coordinate sstem are shown in figure 1. The incompressible Navier-Stokes equations (Brethouwer 216) are solved 4
TABLE I. DNS parameters: N x, N and N z are the number of modes in the streamwise, wallnormal and spanwise direction, respectivel. Ro Re τ Re u τ Re s τ N x N N z 1 1 1 256 385 192.15 976 117 825 234 385 1728.45 8 964 594 248 361 1536.65 7 851 55 192 321 144.9 544 677 365 1536 257 1152 1.2 423 51 326 1152 217 864 with a pseudo-spectral code, as in Brethouwer (216, 217), with Fourier expansions in the streamwise x- and spanwise z-direction and Chebshev polnomials in the wall-normal -direction (Chevalier et al. 27). Together with the Navier-Stokes, the code solves the advection-diffusion equation for a passive scalar Θ t + U Θ = 1 Re P r 2 Θ (1) where U is the dimensionless velocit and Θ the scalar value. In the x- and z-direction, periodic boundar conditions are used for the the scalar and velocit, and no-slip conditions for the velocit at the walls. Further, Θ = at one wall at = 1 and Θ = 1 at the other wall at = 1, where is made non-dimensional with h. The scalar is thus kept at constant but different values at the wall, like in Johansson & Wikström (1999) and Nagano & Hattori (23). In the statisticall stationar state the mean scalar fluxes are equal at both walls. In the DNS the flow rate and thus Re was kept constant at 2 and P r =.71. The domain size was 8πh, 2h and 3π in the x-, - and z-direction, respectivel, and the spatial resolution was similar as in other DNS of turbulent channel flow (Lee & Moser 215). The rotation number Ro was varied from (no rotation) to a quite high value of 1.2. The parameters of the DNS are listed in table I. The friction velocit is defined as u τ = [u 2 τu/2 + u 2 τs/2] 1/2, where u τu and u τs are the friction velocit of unstable and stable channel side, respectivel (Grundestam et al. 28). Re u τ and Re s τ are Renolds numbers based on u τu and u τs, respectivel. Before the statistics were collected, I ran the DNS for a sufficientl long time to reach 5
1.8 1.5 1.2.9.6.3 (a) -1 -.5.5 1 2.8 2.4 2 1.6 1.2.8.4 (c) -1 -.5.5 1 3.2 2.8 2.4 2 1.6 1.2.8.4 (b) -1 -.5.5 1 1.5 -.5-1 -1.5 (d) -1 -.5.5 1 FIG. 2. Profiles of (a) U/U b, (b) u +, (c) v +, and (d) uv +. ( ) Ro =, ( ) Ro =.15, ( ) Ro =.45, ( ) Ro =.65, ( ) Ro =.9. ( ) Ro = 1.2. a statisticall stationar state with a constant mean scalar flux. I experienced that it can take a long time before the scalar field reaches such a state when the channel is rotating. 3. FLOW FIELD Before I present and discuss the DNS results on the scalar transport I will briefl discuss some one-point flow statistics. A more comprehensive stud of rotating channel flow is presented in Brethouwer (216, 217) and other publications cited in the Introduction. Note that in this paper I onl show results for Re = 2 whereas in Brethouwer (217) I often show results for other Re. Here below, U denotes the mean streamwise velocit and u, v and w the streamwise, wall-normal and spanwise velocit fluctuations, respectivel. An overbar implies averaging over time and homogeneous x- and z-directions. Figure 2.(a) shows that all profiles of U/U b for Ro > are skewed and have a linear slope part with du/d 2Ω on the unstable side. In the DNS = 1 and 1 correspond to the wall on the unstable and stable channel sides, respectivel. Figure 2.(b) and (c) show the 6
profiles of the root-mean-square of the streamwise and wall-normal velocit fluctuations u + respective v + scaled b u τ. On the stable side u + is strongl reduced b rotation and the sharp peak near the wall disappears. The maximum of v + grows, respectivel, decas monotonicall with Ro on the unstable and stable channel side while the normalized Renolds shear stress, uv + = uv/u 2 τ, decas with Ro on the stable side (figure 2.d). Both v + and uv + are small or nearl zero on the stable channel side if Ro.9. When Ro is further raised turbulence becomes also weak on the unstable side and if Ro 2.4 turbulent momentum transfer is negligible and the flow approaches a laminar state (Brethouwer 217). In the DNS Re is constant but Re τ varies and decreases with Ro, see table I, owing to the reduced Renolds stresses. 4. VISUALIZATIONS OF THE SCALAR FIELD AND TIME SERIES One-dimensional velocit spectra at Re = 2 and Ro =,.15,.45 and.9 and visualizations of the instantaneous flow field are presented in Brethouwer (217). These spectra and visualizations reveal the presence of large and long streamwise counter-rotating roll cells on the unstable channel side at Ro =.15 with a spanwise size of πh/2. The roll cells leave an imprint on the scalar field visualized in figure 3. Whereas no clear coherent structures can be seen in the scalar field in the non-rotating channel (figure 3.a), at Ro =.15 narrow streaks with low scalar values are observed in the x-z plane on the unstable channel side in the outer laer (figure 3.b). These are approximatel aligned with the flow direction and caused b updrafts between the counter rotating roll cells coming from the lower wall at = 1. With increasing Ro roll cells tend to become smaller and less obvious, accordingl, the streaks with low scalar values become less coherent and the spanwise distance between them diminishes (figure 3.c, d) and at Ro =.9 roll cells are hardl perceptible (figure 3.e). These observations are consistent with the energ spectra presented in Brethouwer (217). Figure 4 shows the instantaneous scalar gradient on the wall at = 1 on the stable channel side. The fluctuations in the scalar gradient impl that the flow is full turbulent here if Ro.15 (figure 4.a, b), but at higher Ro it starts to partl relaminarize on the stable channel side. At Ro =.45 oblique banded patterns develop with alternating turbulent and laminar-like flow (Brethouwer 217). A corresponding pattern emerges in the scalar gradient field with oblique banded regions with strong scalar gradient fluctuations 7
(a) (b) (c) (d) (e) FIG. 3. Visualizations of the instantaneous scalar field in an x-z plane at (a) Ro = and =.6, (b) Ro =.15 and =.6, (c) Ro =.45 and =.8, (d) Ro =.65 and =.8 and (e) Ro =.9 and =.85. Dark colours correspond to low scalar values. induced b the turbulence and regions where these gradient fluctuations are mostl absent since the flow is locall laminar-like (figure 4.c). Similar oblique turbulent-laminar patterns have been observed in several flow tpes, see for example, Duguet et al. (29), Brethouwer et al. (212) and Deusebio et al. (214). A distinguishing feature is that the patterns do not span the whole channel in the present case but are onl found on the stable side. At a higher 8
(a) (b) (c) (d) (e) FIG. 4. Visualizations of the instantaneous scalar gradient at = 1 at (a) Ro =, (b) Ro =.15, (c) Ro =.45, (d) Ro =.65 and (e) Ro =.9. Dark colours correspond to high scalar gradients. Ro =.65 the patterns with turbulence and strong scalar gradient fluctuations become less coherent, as seen in figure 4.(d). The also appear to become more unstead and their size varies in time, as indicated b other visualizations (not shown) and temporal variations in the mean scalar gradient and skin friction at = 1, shown later. When Ro.9 fluctuations are weak on the stable side and small-scale fluctuations in the scalar gradient field cannot be observed (figure 4.e), impling that near-wall turbulence is suppressed b the sstem 9
(a) (b) (c) FIG. 5. Visualizations of the instantaneous scalar field in a -z-plane at (a) Ro =, (b) Ro =.45 and (c) Ro =.9. Red and blue colours correspond to low and high scalar values, respectivel. rotation. Instantaneous plots of the scalar field in a -z-plane at Ro =,.45 and.9 are shown in figure 5. At Ro =.45, the visualization reveals large-scale plume-like structures in the scalar field (figure 5.b), which are likel caused b the large streamwise roll cells. At higher Ro and Ro =, the roll cells are absent or much smaller (Brethouwer 217) and accordingl no or less large plumes are seen in the scalar field (figure 5.a, c). Large jumps in the scalar field are observed when Ro > in the core region near the border of stable and unstable channel side where the mean scalar gradient is large, as discussed later. It is important to remark that at Re = 2 a linear instabilit of Tollmien-Schlichtinglike waves occurs when Ro.9 (Brethouwer 216). This instabilit causes recurring bursts of turbulence mostl confined to the stable channel side and its mechanism is examined in Brethouwer et al. (214) and Brethouwer (216). The instabilit and strong bursts naturall affect the scalar field, as shown b the time series in figure 6. In this figure, time series of the volume averaged turbulent kinetic energ K m and root-mean-square of the scalar fluctuations θ m are presented as well as the wall shear stress τ ws at = 1 and wall temperature gradients (dθ/d) u and (dθ/d) s at = 1 and = 1, respectivel. The stress and gradients are averaged over the wall and the time t is scaled b h/u b. When Ro.45 all these quantities are approximatel constant in time (figures 6.a) but at Ro =.65, τ ws and (dθ/d) s displa noticeable variations caused b the growing and shrinking of the turbulent areas on the stable channel side, as discussed before. Dai et al. (216) have observed similar variations 1
1.5 1.5 1 1.5.5 (a) 1 3 4 t (b) 1.5 1.8 1 1.2.5.6 8 16 24 (d) t 2 4 6 8 15 2 t 2.4 5 1 t Time series at (a) R =.45, (b) R =.65, (c) R =.9 and (d) R = 1.2 of the volume averaged Km ( (dθ/d)s ( 2 2 (c) FIG. 6. ) and θm ( ), and wall gradients τws ( ), (dθ/d)u ( ) and ). τws is normalized with the time-averaged wall shear on both wall whereas the other quantities are normalized with their respective time-averaged values. of τws in their DNS and attributed these to the dnamics and changes of the streamwise roll cells. At Ro =.9 and 1.2 the time series of θm and (dθ/d)s show simultaneousl with Km and τws recurring sharp peaks with a time interval of about 7t, whereas (dθ/d)u and the wall shear stress at = 1 (not shown) sta approximatel constant. This shows that the bursts of turbulence triggered b the linear instabilit causes sharp and significant but relativel short surges in the scalar fluctuations and scalar flux at the wall on the stable channel side. Such recurring bursts are also observed at higher Ro. However, when Ro 2.1 and the turbulence is weak in the whole channel (Brethouwer 217) the bursts trigger sharp spikes in the scalar flux not onl at the wall at = 1 but also at the other wall (not shown). Although interesting, I will not further explore these recurring bursts in the DNS at Ro.9 but focus in the rest of the paper in these DNS on the relativel calm periods between the bursts. From a modelling point of view these calm periods are more relevant 11
1.8.6 Θ.4.2 (a) -1 -.5.5 1 7 6 5 4 θ + 3 2 1 (b) -1 -.5.5 1 FIG. 7. Profiles of (a) Θ and (b) θ +. Lines as in figure 2. and give a better fundamental insight into the influence of rotation on turbulent scalar transport. Including the bursts would add further complexit to the problem. I therefore compute the scalar statistics from the statistics collected during the calm periods between the bursts. In Brethouwer (217) I explain in more detail how I exclude the periods with the bursts from the statistics. 5. PASSIVE SCALAR TRANSPORT In this section, I will discuss the basic statistics of the scalar field and scalar transport. Here below, Θ is used to denote the mean scalar value, i.e. averaged over time and x- and z-directions, and θ the scalar fluctuation. A superscript + implies, unless stated otherwise, scaling in terms of viscous wall units ν and u τ and θ τ = Q w /u τ where Q w = α(dθ/d) w is the mean scalar flux at the wall, which is equal at both walls in the statisticall stationar state. Figure 7.(a) shows profiles of Θ at different Ro. The mean scalar value Θ goes down monotonicall with Ro on the unstable channel side, whereas on the stable side the mean scalar gradient is steep in rotating channel flow, much steeper than on the unstable side. This can be understood b considering the mean scalar transport across the channel in the stead state α dθ d vθ = Q w. (2) On the stable channel side, vθ is strongl diminished in rotating channel flow, as shown 12
25 2 15 U + 1 5 (a) 1 1 1 1 2 1 3 + 2 16 12 Θ + 8 4 (b) 1 1 1 1 2 1 3 + FIG. 8. Profiles of (a) U + and (b) Θ + as function of +. The straight green and blue dashed lines in (b) are explained in the text. Other lines as in figure 2. later, which naturall implies a large αdθ/d to maintain the scalar flux balance. Profiles of the root-mean-square of the scalar fluctuations θ + are shown in figure 7.(b). At Ro = the profile of θ + has not onl peaks near both walls caused b intense near-wall turbulence but also a maximum at the centre as a result of a large dθ/d and, consequentl, a high production of scalar fluctuations here, as shown later. When.15 Ro.65, θ + including the near-wall peak on the unstable channel side declines with Ro and grows on the stable channel side. This growth is caused b the large dθ/d and related large production of scalar fluctuations on the stable side. Note that scalar fluctuations are strong although turbulence is weakened b rotation on the stable channel side. At higher Ro, θ + declines in the near-wall region of the stable side. Instead, θ + has a pronounced maximum near the position where dθ/d strongl increases going from the unstable towards the stable channel side. Quite similar trends regarding the profiles of Θ and θ + are also observed at lower Renolds numbers (Nagano & Hattori 23, Liu & Lu 27). Figure 8 shows the mean velocit U + and scalar Θ + profiles on the unstable channel side as function of the distance to the wall in wall units +. The mean velocit profile at Ro = follows approximatel a log-law behaviour between the near-wall and centre region, but a similar log-law behaviour is not observed when Ro > and cannot be expected from theoretical arguments since the velocit profile depends on Ro. Previous studies have shown that in the overlap region of non-rotating channel flow the Θ + -profile has a similar log-law behaviour as the velocit (Kawamura et al. 1999, Pirozzoli et al. 216). Accordingl, in the 13
present DNS at Ro =, the mean scalar profile between the near-wall and centre region approximatel matches Θ + = 1 log + + C θ (3) κ θ with κ θ =.42 and C θ = 2.8 given b the straight green dashed line (figure 8.b). The present value of κ θ =.42 is similar to the values.43 and.46 found b Kawamura et al. (1999), respectivel, Pirozzoli et al. (216). When Ro.45 the profiles of Θ + on the unstable side awa from the wall approximatel also follow the logarithmic profile (3) but with κ θ = 2. and decreasing values of C θ with Ro, given b the straight blue dashed lines in figure 8.(b), despite the absence of a similar log-law behaviour in the U + -profile. A closer inspection shows that this log-law region partl overlaps with the region where the mean velocit profile is approximatel linear and du/d 2Ω. Whether this log-law behaviour of the Θ + -profiles with the same logarithmic slope in rotating channel flows is just a coincidence is not et clear and there is obvious explanation for this behaviour. I will now consider turbulent scalar transport. Figure 9.(a) and (b) show the profiles of the streamwise and wall-normal turbulent scalar fluxes uθ + respective vθ + in wall units. In non-rotating channel flow uθ + has, except around the centreline, a larger magnitude than vθ + (Johansson & Wikström 1999) meaning that the turbulent scalar flux is mainl aligned with the streamwise direction and the same applies to the stable channel side in the rotating cases. When Ro rises, uθ + including its near-wall peak monotonicall declines on the unstable side while vθ + stas near unit, meaning that the turbulent scalar flux vector turns towards the wall-normal direction and is nearl aligned with the wall-normal direction awa from the wall in the region where du/d 2Ω. On the stable channel side, vθ + declines with Ro and is ver small for Ro.9 when the flow relaminarizes there, whereas uθ + stas large for Ro.65 but declines at higher Ro and its near-wall peak disappears (figure 9). The scalar transport on the stable channel side is thus mainl diffusive at higher Ro. These results are consistent with previous DNS of scalar transport in rotating channel flow at lower Re and with a smaller range of Ro (Nagano & Hattori 23, Liu & Lu 27). The are also consistent with DNS and rapid distortion theor of homogeneous turbulent shear flow subject to spanwise sstem rotation. In that case, the turbulent scalar flux vector is mostl aligned with the streamwise direction when the rotation is absent, or cclonic as on the stable channel side. When the rotation is anti-cclonic as on the unstable channel side the turbulent scalar flux vector becomes more aligned with the wall-normal direction and is 14
8 6 4 2 uθ + -2-4 -6-8 (a) -1 -.5.5 1 -.5 vθ + -1 (b) -1 -.5.5 1 1.5 wθ + (c) -1 -.5.5 1 FIG. 9. Profiles of (a) uθ + and (b) vθ +. Lines as in figure 2. (c) Profile of wθ + at ( ) Ro =.45 Ro =.45 and ( ) Ro =.65. almost full aligned with the mean scalar gradient when the mean shear du/d equals 2Ω, where Ω is the imposed rotation (Brethouwer 25, Kassinos et al. 27), like in the present case. In non-rotating turbulent channel flow the mean spanwise turbulent scalar flux wθ + is naturall zero. But figure 9.(c) shows that when Ro =.45 or.65 and oblique turbulent patterns are present in the near-wall region of the stable side (figure 4.c, d) wθ + is significant in this region. In these cases, uw is also non-zero. The oblique patterns induce thus a noticeable spanwise turbulent momentum and scalar transport near the wall. Figure 1 shows the correlation coefficients ρ uθ = uθ/(u θ ) and ρ vθ = vθ/(v θ ), where a prime denotes root-mean-square values. Consistent with Abe & Antonia (29) and Pirozzoli et al. (216), ver high streamwise flux correlations ρ uθ are found at Ro = near the wall owing the ver high similarit between the streamwise velocit and scalar field 15
1.75.5.25 ρ uθ -.25 -.5 -.75-1 (a) -1 -.5.5 1 1.8.6 ρ vθ.7.6.5.4.3.2.1 (b) -1 -.5.5 1 1.8.6 ρ uθ.4 ρ vθ.4.2.2 (c) 5 1 15 2 + (d) 5 1 15 2 + FIG. 1. Profiles of (a,c) ρ uθ and (b,d) ρ vθ, (c) and (d) show a close-up near the wall on the unstable side. Lines as in figure 2. (figure 1.c), while in the outer laer, except around the centreline, ρ uθ is lower et still large (figure 1.a). A large correlation means here large positive or large negative values of ρ uθ. In the rotating cases, the magnitude of ρ uθ remains large on the stable channel side, but on the unstable side it declines rapidl in the outer laer with Ro and even becomes negative at high Ro (figure 1.a). Close to the wall on the unstable side, ρ uθ declines when Ro.65 but remains quite large (figure 1.c). Thus, rotation reduces the similarit between the u- and θ-field on the unstable side, especiall in the outer laer. Yang et al. (211) came to a similar conclusion for DNS of rotating channel flow at much lower Re. The high near-wall correlation at Ro = is motivated b the similarit between the governing equations for u and θ with the main difference being the pressure gradient term (Abe & Antonia 29). However, if Ro > an additional Coriolis term appears in the former, which diminishes the similarit between the governing equations and, accordingl, the correlation between u and 16
3.5 3 2.5 u + 2 max θ 1.5 rms + 1.5.3.6.9 1.2 Ro FIG. 11. The maximum near-wall value of u + ( ) and θ + ( ) on the unstable (blue and red lines/smbols) and stable channel side (green lines/smbols) as function of Ro. Blue lines and smbols: using wall units based on u τ. Red lines and smbols: using wall units based on u τu. Green lines and smbols: using wall units based on u τs. θ. B contrast, ρ vθ is quite insensitive to Ro on the unstable channel side awa from the wall with values between.4 and.5 whereas it declines with Ro on the stable side (figure 1.b). On the other hand, near the wall on the unstable side, ρ vθ increases considerabl with Ro and becomes quite large (figure 1.d). The behaviour and magnitude of ρ vθ with Ro follows here closel the correlation coefficient for the Renolds shear stress uv/(u v ) (not shown) which also becomes large near the wall at higher Ro. These high correlations might be related to the formation of elongated streamwise near-wall vortices in rotating channel flow on the unstable side (Yang & Wu 212). To conclude, the alignment of the scalar flux with the mean scalar gradient with Ro on the unstable side is related to the reduced correlation between streamwise velocit and scalar fluctuations whereas the reduced turbulent scalar fluxes on the stable side are caused b both declining v-θ correlations and Renolds stresses. Figure 11 shows the near-wall peak values of u + and θ + as function of Ro. The velocit and scalar fluctuations on the unstable side are either scaled b u τ and θ τ = Q w /u τ, respectivel, or b u τu and θ τu = Q w /u τu, respectivel, while on the stable channel side the are scaled b u τs and θ τs = Q w /u τs, respectivel. The high similarit between the u and θ field in the near-wall region in non-rotating channel flow (Abe & Antonia 29) is reflected 17
b the closeness of the peak values of u + and θ +, although the peak value of u + is slightl higher than that of θ +, as a consequence of a smaller than unit P r (Pirozzoli et al. 216). In rotating channel flow the difference between the peak values of u + and θ + is considerable if the scaling is based on u τ and θ τ, but if the scaling of the peak values on the unstable and stable channel side is based on u τu, θ τu and u τs, θ τs, respectivel, the differences are much smaller. In the latter case, the peak values of u + and θ + on the unstable channel side show a ver similar downward trend with Ro. The trend on the stable channel side is non-monotonic, which ma be related to the appearance of turbulent-laminar patterns at Ro =.45 and.65. Rotation reduces thus the similarit between the u- and θ-field on the unstable side, but their near-wall peak values are still highl correlated and displa a similar scaling in terms of wall units. 6. BUDGET EQUATIONS In this section, the budgets of the governing equations of the scalar energ K θ = 1 2 θθ and uθ and vθ are considered to obtain insights into the generation of scalar fluctuations and fluxes. Budgets for scalar transport in non-rotating channel have been studied b e.g. Johansson & Wikström (1999) and in spanwise rotating channel flow b Nagano & Hattori (23) and Liu & Lu (27), albeit at low Renolds numbers. In the present stead-state case, the governing equation of the scalar energ K θ reads = vθ Θ } {{ } P θ 1 vθ 2 2 } {{ } Dθ t where P θ represents production, D t θ and Dm θ + α 2 K θ 2 }{{} D m θ and ε θ dissipation. Figure 12 shows the budgets P + θ, Dt+ θ and ε + θ.15,.65 and 1.2. The molecular diffusion D m θ α θ θ, (4) x k x }{{ k} ε θ turbulent and molecular diffusion, respectivel, in wall units at Ro =, is not shown since it is small. The turbulent diffusion Dθ t appears onl significant near the wall and at Ro = 1.2 in the region awa from the wall where K θ and P θ are large, impling that there is primaril a balance between P θ and ε θ. Accordingl, at Ro = and Ro.9 the ratio P θ /ε θ is near unit in the outer region awa from the border between the stable and unstable channel side (figure 13.b). However, P + θ and ε + θ are both ver small on the unstable side awa from the wall when.15 Ro.65 (figure 13.a) since the mean scalar gradient is small. In these cases the 18
.2.2.1.1 -.1 (a) -1 -.5.5 1.2 -.1 (b) -1 -.5.5 1.2.1.1 -.1 (c) -1 -.5.5 1 -.1 (d) -1 -.5.5 1 FIG. 12. Budgets of the K θ -balance equation in wall units at (a) Ro =, (b) Ro =.15, (c) Ro =.65 and (d) Ro = 1.2. ( ) P + θ, ( ) ε+ θ, ( ) Dt+ θ..2 1.5 1 ε + θ.2 P θ ε θ.5.2 (a) -1 -.5.5 1 (b) -1 -.5.5 1 FIG. 13. Profiles of (a) ε + θ and (b) P θ/ε θ. Lines as in figure 2. 19
ratio P θ /ε θ deviates significantl from unit (figure 13.b) impling that Dθ t is significant. This disbalance between P θ and ε θ and the considerable diffusion of K θ are likel caused b the large-scale streamwise roll cells being present at these rotation rates but which are much less dominant at higher Ro (Brethouwer 217). Liu & Lu (27) show that these roll cells indeed have a large impact on scalar transport in rotating channel flows. The maximum value of P + θ close to the walls is independent of Ro and close to its theoretical value P r/4 (Johansson & Wikström 1999). At Ro =.15 and even more so at Ro =.65, both ε + θ P + θ are large on the stable channel side (figure 12.b, c) where also the mean scalar gradient and K θ are large (figure 7). At higher Ro, the peaks of ε + θ and P + θ on the stable side move awa from the wall to the position, going from the unstable towards the stable side, where dθ/d steepens while vθ is still large (figure 12.d). Accordingl, the maximum of θ + at high Ro moves awa from the wall at high Ro (figure 7.b). Moving closer to the wall on the stable side, vθ and consequentl P + θ and ε + θ diminish rapidl. The respective governing equations of uθ and vθ read in the present case = uv Θ U vθ uvθ + ( αu θ ) u + νθ + p θ u θ (α + ν) + ρ x x }{{}}{{}}{{} k x k } 2Ωvθ {{}}{{} C 1 P 1 D1 t D1 m Π 1 = v 2 Θ }{{} P 2 v2 θ }{{} D t 2 pθ }{{} D p 2 and (5) + ( αv θ ) v + νθ + p θ v θ (α + ν) ρ x }{{} k x k } 2Ωuθ {{}}{{} C 2 D2 m Π 2 where P 1, P 2 represent production, D t 1, D t 2 turbulent diffusion, D p 2 pressure diffusion, and D m 1, D m 2 molecular diffusion. The sum of pressure scalar-gradient correlation and diffusive and viscous dissipation, Π 1 and Π 2, is considered because diffusive and viscous dissipation are generall considered to be small and therefore often added to the pressure scalar-gradient correlation term in turbulence modelling (Wikström et al. 2). The Coriolis force leads to the additional terms C 1 and C 2 in the governing equations of the scalar fluxes. Figure 14 shows the budgets P + 1, Π + 1, D t+ 1 and C + 1 of the governing equation (5) of uθ in wall units at Ro =,.15,.65 and 1.2. The molecular diffusion D m 1 (6) is again not shown because of its smallness. At Ro =, P + 1 is in the whole channel mainl balanced b Π + 1 and the same applies to the near wall region at Ro =.15 (figure 14.a, b). When.6.4 all budget terms are small in the latter case, but a closer inspection reveals that the Coriolis 2
.45.3.15 -.15 -.3 -.45 (a) -1 -.5.5 1.3.2.1 -.1 -.2 (c) -1 -.5.5 1.6.4.2 -.2 (b) -1 -.5.5 1.3.15 -.15 -.3 (d) -1 -.5.5 1 FIG. 14. Budgets of the uθ-balance equation in wall units at (a) Ro =, (b) Ro =.15, (c) Ro =.65 and (d) Ro = 1.2. ( ) P + 1, ( ) Π+ 1, ( ) Dt+ 1, ( ) C+ 1. term C + 1 is of the same order as P + 1 and Π + 1 and negative since vθ is negative. At a higher Ro =.65, P + 1 and Π + 1 are large and approximatel balance each other in the near-wall region on the unstable side and a large part of the stable side (figure 14.c). On the other hand, in the outer region of the unstable side, C + 1 is significant and balanced b P + 1 and contributes to a small uθ in this part of the channel (figure 9.a). When Ro is raised to 1.2 turbulence disappears on the stable side and all budget terms are small there (figure 14.d). In the major part of the unstable channel side D t+ θ and Π + 1 are small but P + 1 and C + 1 are both large and balance each other. The balance of P + 1 and C + 1 in rotating channel flow on the unstable channel side can be understood b considering the production term more closel. This term has two contributions, namel uv( Θ/ ) and vθ( U/ ), see equation (5). The first part is small on the unstable side awa from the wall because Θ/ is small. The second part vθ( U/ ) is large but is nearl balanced b C 1 = vθ2ω since 21
.15.1.5 -.5 -.1 (a) -1 -.5.5 1.3.2.1 -.1 (b) -1 -.5.5 1.3.2.15.1 -.1 -.2 -.15 (c) -1 -.5.5 1 -.3 (d) -1 -.5.5 1 FIG. 15. Budgets of the vθ-balance equation in wall units at (a) Ro =, (b) Ro =.15, (c) Ro =.65 and (d) Ro = 1.2. ( C + 2. ) P + 2, ( ) Π+ 2, ( ) Dt+ 2 + Dp+ 2, ( ) Π+ 2 + Dp+ 2, ( ) U/ 2Ω on the unstable side. The term induced b the Coriolis force has thus a major impact on the streamwise scalar flux and contributes to the alignment between the turbulent scalar flux vector and mean scalar gradient in rotating channel flow. Figure 15 shows the budgets P + 2, Π + 2, the sum D t+ 2 + D p+ 2 and C + 2 in the governing equation (6) of vθ in wall units at Ro =,.15,.65 and 1.2. The molecular diffusion D m 2 again small. At Ro = the dominant contributions come from P + 2 is and Π + 2, except ver near the walls and near the centreline where diffusion D t+ 2 + D p+ 2 is significant. At Ro =.15 diffusion and C + 2 are significant as well as P + 2 and Π + 2 in the main part of the channel. The Coriolis term C + 2 contributes to the production of negative vθ on the unstable side but counteracts it on the stable channel side since uθ changes sign. At a higher Ro =.65, C + 2 is important and large on the stable side and near the wall on the unstable side but 22
small awa from the wall since uθ is small. In a large part of the unstable channel side there is thus primaril a balance between P + 2 and Π + 2. At a higher Ro = 1.2, P + 2 is nearl constant awa from the wall on the unstable side and balanced b Π + 2 and C + 2 while it rapidl declines approaching the wall on the stable channel side where both Π + 2 and C + 2 are large (figure 15.d). Note that C 2 is positive on the stable side in rotating channel flow and thus contributes to the reduced wall-normal turbulent scalar transport on that side. On the stable channel side, both Π + 2 and diffusion dominated b D p+ 2 displa large variations at Ro =.65 and 1.2. B contrast, the profile of the sum Π + 2 + D p+ 2 = ((θ p/ )/ρ) +, which is also displaed in figure 15, is much smoother and does not show extreme values near the wall in rotating channel flow. This suggests that it might be easier to model the latter term instead of Π 2 and D p 2 separatel. Figure 16 shows in more detail the production terms of turbulent kinetic energ P K and P θ near the wall on the unstable side at Ro = and.65. Furthermore, the figure shows the sum P 2 + C 2, which can be regarded as a total production of vθ, and P 1θ = uv( Θ/ ) and P 1u = vθ(2ωvθ U/ ), which are the production owing to the mean scalar gradient, respective, the production b mean shear together with Coriolis term contribution in the governing equation (5) of uθ. The latter term, P 1u, is small awa from the wall on the unstable side of a rotating channel when U/ 2Ω, as explained before. All terms and the distance to the wall are in wall units. The production terms are premultiplied b to accentuate the behaviour awa from the wall. Here, u uτ is used instead of u τ for the scaling at Ro =.65 because this velocit scale is more appropriate on the unstable side. A DNS of scalar transport with P r = 1 in a non-rotating channel flow (not presented here) shows that the profiles of PK, P θ, P 1θ and P 1u including the maxima all collapse near the wall, see also Pirozzoli et al. (216). Figure 16.(a) shows that in the present case with P r =.71 the profiles do not collapse when Ro = since the maximum of PK approaches 1/4 whereas the maximum of Pθ is P r/4 as explained before. The maxima of P 1θ and P 1u are in between these two previous maxima and approximatel equal to each other, showing that production owing to the mean scalar gradient and mean shear are equall important. All these maxima are nearl independent of Ro (compare figure 16.a and b). The profiles of P1θ and P 1u approximatel collapse for 2 (figure 16.b), also at other Ro, while the do not collapse if the Coriolis term is not included in P1u, which is somewhat remarkable. 23
4 3 2 4 3 2 1 1 (a) 2 1 1 1 (b) 2 1 1 6 FIG. 16. Premultiplied production terms in wall units on the unstable side at (a) Ro = and (b) Ro =.65. ( ) P K, ( ) P θ, ( ) P 1θ, ( ) P 1u, ( ) (P 2 + C 2 ). 7. EFFICIENCY OF SCALAR TRANSPORT In this section, some quantities characterising the efficienc of scalar transport are discussed; some of them are directl relevant for engineering. An important quantit in engineering is the Nusselt number. It is here defined as the ratio of wall-normal scalar flux in the present turbulent channel flows and laminar channel flow, i.e. Nu = 2h α Θ Q w, (7) where Θ is the imposed scalar difference at the walls. With this definition, Nu = 1 for laminar channel flow. However, Nu does not distinguish between the scalar transport on the stable and unstable channel sides. To examine the effectivit of scalar transport on both channel sides, I have also computed a Nusselt number defined similarl as in some previous studies, e.g. Pirozzoli et al. (216), to account for difference in the mean scalar gradients on both sides. For the unstable channel side it is defined as Nu u = 5 8 Q δ w α Θ w Θ m. (8) Here, δ = w where is the position where the total shear stress, i.e. the sum of the viscous and turbulent shear stress, is zero, and w the wall position and Θ w the scalar value at the wall on the unstable channel side. I consider the place separating the stable and 24
45 35 25 Nu 15 5 (a).3.6.9 1.2 Ro 2St/C f 2.5 2 1.5 1.5 (b).3.6.9 1.2 Ro FIG. 17. (a) Nu (blue line/smbols), Nu u (red line/smbols) and Nu s (green line/smbols). (b) 2St/C f on the unstable (red line/smbols) and stable channel side (green line/smbols). unstable channel sides. Further Θ m = 1 δ w UΘ U m d, U m = 1 δ w Ud. (9) The Nusselt number for the stable channel side, Nu s, is defined in the same wa with the integrals in (9) from to w where w and Θ w are now the wall position and the scalar value at the wall on the stable channel side. The factor 5/8 in (8) and in the similar expression for Nu s ensures that Nu u = Nu s = 1 if the flow is laminar. The main difference between the expressions for Nu and Nu u, Nu s is that Θ/2h is replaced b Θ w Θ m /δ as an effective mean scalar gradient to characterize both channel sides. Figure 17(a) shows Nu together with Nu u and Nu s as function of Ro. Nu slightl grows at first but then declines rapidl with Ro, impling that rotation drasticall inhibits the cross channel scalar transport. The behaviour of N u in the present case is different from that in a rotating channel flow at a lower Re τ = 194 when Nu monotonicall decreases with Ro (Liu & Lu 27). The high values of Nu u for.15 Ro.9 and rapid decline of Nu s with Ro reflect the rapid and slow turbulent scalar transport on the unstable and stable channel sides, respectivel. Both Nu and Nu s approach unit at high Ro meaning that diffusive scalar transport becomes significant. Figure 17(b) shows the ratio of the Stanton number to skin friction coefficient 2St/C f where St = Q w U m T m, 25 C f = 2u2 τ U m. (1)
Here, U m, T m and u τ are either computed for the unstable or the stable channel side according to equation (9). If the Renolds analog for the momentum and scalar transport is valid 2St/C f should be near unit (Abe & Antonia 217), which indeed it is in the non-rotating channel (figure 17.b). On the stable channel side 2St/C f stas near unit while on the unstable channel side it grows with Ro and clearl deviates from unit, suggesting that the Renolds analog is valid for the stable channel side but not for the unstable side where scalar transport is relativel rapid. The ratio of scalar to turbulence time-scale r = K θ/ε θ K/ε (11) is in some turbulence models assumed to be a constant. Antonia et al. (29) found that in non-rotating turbulent channel flow with a passive scalar with P r =.71, r is indeed nearl constant; its value is between.5 and.6 in the outer region. In rotating homogeneous shear flow, however, r varies with the imposed sstem rotation rate (Brethouwer 25). Figure 18(a) shows the profiles of r in rotating channel flow. The variation of r across the channel in the present DNS at Ro = is larger than in the DNS b Antonia et al. (29), which is the consequence of a different forcing of the scalar field in those two DNS. Nevertheless, in the outer region but outside the centre region r is close to that observed b Antonia et al. (29). In rotating channel flow, r at first grows with Ro and varies between.7 and.85 at Ro =.15, but then declines somewhat on the unstable channel side and varies there between approximatel.6 and.7. Thus, the time-scale ratio r displas some variations with Ro but these are not ver large. A measure of turbulence to scalar fluctuation intensit is given b the dimensionless parameter B = K/dU/d Kθ /dθ/d (12) Figure 18(b) shows the profiles of B in the present DNS. The values of B of approximatel 1.2 to 1.3 at Ro = in the outer regions.95.6 and.4.95 are the same as those observed b Antonia et al. (29) in their DNS of non-rotating turbulent channel flow at similar Re. In the centre region, B is much larger in the present DNS since the mean velocit gradient is small, unlike the fluctuations and mean scalar gradient. With increasing Ro, however, B declines strongl in the outer region of the unstable channel side to less than.4, meaning that rotation augments scalar fluctuations relative to velocit fluctuations. The 26
1.2 1.6 1.8 1.2 r.6 B.8.4.2.4 (a) -1 -.5.5 1 (b) -1 -.5.5 1 1.2 1.8 P r T.6.4.2 (c) -1 -.5.5 1 FIG. 18. Profiles of (a) r, (b) B and (c) P r T. Lines as in figure 2. same trend is observed in rotating homogeneous shear flow with a passive scalar (Brethouwer 25). Especiall when the mean shear du/d equals 2Ω, as in the absolute mean vorticit region in the present DNS, scalar fluctuations are relativel strong, although neither in the governing equation of K nor K θ there is a direct influence of rotation. The small values of B in rotating channel flow, however, suggest that P θ is relativel large compared to the production of K. An important parameter in turbulence modelling, often assumed to be a constant, is the ratio of turbulent viscosit to scalar diffusivit given b the turbulent Prandtl number P r T = ν T α T = uv dθ/d vθ du/d. (13) It can also be considered as a relative measure of the production term of K to K θ. Figure 18(c) shows the profiles of P r T obtained from the present DNS. At Ro =, P r T.85 in the outer region, consistent with the values computed b Antonia et al. (29) and Pirozzoli 27
et al. (216) for non-rotating channel flow. On the other hand, in rotating channel flow P r T is considerabl smaller in the outer region of the unstable channel side; it is less than.2 for Ro.45. Also on the stable channel side at Ro =.15, P r T is smaller than in the non-rotating case awa from the wall. These results show, like the ratio 2St/C f in figure 17.(b), that the Renolds analog for scalar-momentum transfer does not hold for rotating channel flow and that turbulent scalar transport is ver efficient compared to momentum transfer on the unstable channel side. The low value of P r T also implies that the production of K θ is relativel strong to that of K, consistent with the small values of B observed in rotating channel flow. Small values of P r T are also observed in rotating homogeneous shear flow when the imposed shear du/d = 2Ω, where Ω is the imposed sstem rotation rate (Brethouwer 25), consistent with the present results. The small values of B and P r T suggest that the scalar and velocit field are dissimilar in rotating channel flow, as alread indicated b the correlation coefficients in figure 1. This dissimilarit is further explored b spectra in the next section. 8. SPECTRA Premultiplied one-dimensional spanwise spectra of turbulent kinetic energ k z E K (k z ) and scalar variance k z E θθ (k z ) as well as premultiplied one-dimensional streamwise spectra of turbulent kinetic energ k x E K (k x ) and scalar variance k x E θθ (k x ) at Ro =,.15,.45 and.9 are shown in figure 19 as function of the spanwise and streamwise non-dimensional wave length λ z and λ x, respectivel. Here, k x and k z are the streamwise and spanwise wave number, respectivel. The streamwise and spanwise spectra are normalized such that E β (k x )dk x = E β (k z )dk z = 1, (14) respectivel (β stands for K or θθ), and are computed at three different wall-normal positions on the unstable channel side, i.e. near the wall at 1, at 1 and in the outer region at.55. A superscript implies scaling b ν/u τu, which is appropriate for the spectra on the unstable channel side. At Ro =, k z E θθ (k z ) shows near the wall a peak at λ z 1, corresponding to the turbulent near-wall ccle, while in the outer region it has a peak at scales larger than λ z = h (figure 19.a), like the scalar spectra presented b Pirozzoli et al. (216) for non-rotating channel 28
.3.8.6 k z E K (k z ) k z E θθ (k z ).2.1 k x E K (k x ) k x E θθ (k x ).4.2 (a) 1 1 1 2 1 3.4 λ z (b) 1 2 1 3 1 4.8 λ x k z E K (k z ) k z E θθ (k z ).3.2.1 k x E K (k x ) k x E θθ (k x ).6.4.2 (c) 1 1 1 2 1 3 1 4.3 λ z (d) 1 2 1 3 1 4.8 λ x k z E K (k z ) k z E θθ (k z ).2.1 k x E K (k x ) k x E θθ (k x ).6.4.2 (e) 1 1 1 2 1 3.3 λ z (f) 1 2 1 3 1 4.8 λ x k z E K (k z ) k z E θθ (k z ).2.1 k x E K (k x ) k x E θθ (k x ).6.4.2 (g) 1 1 1 2 1 3 λ z (h) 1 2 1 3 1 4 λ x FIG. 19. Premultiplied one-dimensional spanwise and streamwise spectra as function of λ z respective λ x at Ro = (a, b), Ro =.15 (c, d), Ro =.45 (e, f) and Ro =.9 (g, h). Left column shows the spanwise spectra k z E K (k z ) (red lines) and k z E θθ (k z ) (black lines). Right column shows streamwise spectra k x E K (k x ) (red lines) and k x E θθ (k x ) (black lines). The spectra are at 1 (dotted lines), 1 (dashed lines) and.55 (solid lines). The vertical dashed line indicates scales with λ x = λ z = 1. 29