Direct Numerical Simulation of a Turbulent Boundary Layer with Strong Pressure Gradients
TThis draft was prepared using the LaTeX style file belonging to the Journal of Fluid Mechanics Direct Numerical Simulation of a TurbulentBoundary Layer with Strong PressureGradients
Riccardo Balin † , and K. E. Jansen Ann and H. J. Smead Aerospace Engineering Sciences, University of Colorado Boulder,Boulder, CO 80309, USA(Received xx; revised xx; accepted xx)
The turbulent boundary layer over a Gaussian shaped bump is computed by directnumerical simulation (DNS) of the incompressible Navier-Stokes equations. The two-dimensional bump causes a series of strong pressure gradients alternating in rapidsuccession. At the inflow, the momentum thickness Reynolds number is approximately1 ,
000 and the boundary layer thickness is 1 / Key words:
1. Introduction
Turbulent boundary layers undergoing pressure gradients and separation have beenthe subject of a large number of studies due to their ubiquity in science and engineering.Indeed, a deep understanding of their complex physics is necessary for the accurateprediction and design of many engineering systems. While still limited to low Reynoldsnumbers and small domains, the continuous improvement in computational resources hasmade investigation by direct numerical simulation (DNS) increasingly popular. DNS havenot only been of great complement to experimental work in the study of boundary layerphysics, but have also provided detailed and dense high fidelity data for the evaluation † Email address for correspondence: [email protected] a r X i v : . [ phy s i c s . f l u - dyn ] O c t R. Balin and K. E. Jansen and improvement of all lower fidelity turbulence models, from large eddy simulation(LES) to Reynolds Averaged Navier-Stokes (RANS) closures (Abe et al. et al. et al. et al. (2012); Abe (2017); Kitsios et al. (2016); Coleman et al. (2018), performingsimulations of the same kind at increasingly higher values of the momentum thicknessReynolds number Re θ .Flat plate boundary layer flows provide valuable insight into smooth-body shallowseparation caused by continued APG effects, which is a well-known deficiency of RANSclosures. However, they lack two important characteristics that are often present in en-gineering applications, namely strong favorable pressure gradients, in particular a strongFPG upstream of the APG region, and strong streamline curvature. The effects of thesetwo phenomena on boundary layer turbulence were initially investigated experimentallywith flows over bumps and hills. Tsuji & Morikawa (1976) reported the breakdown of thelogarithmic law in the FPG region and the formation of internal layers at the locationswhere the pressure gradient changed sign. Baskaran et al. (1987, 1991) conducted adetailed investigation of streamwise pressure gradient and curvature effects for the flowover a convex hill. A decoupling of the inner and outer layers was observed, with thelatter being isolated and behaving as a free shear layer. Curvature effects were mostpronounced in the outer layer, affecting the turbulence structure and rapidly decreasingthe Reynolds shear stress. Conversely, changes in pressure gradient initially appeared inthe inner region. Webster et al. (1996) conducted experiments of the flow over a circularbump. They reported a significant deviation of the velocity in the FPG region above thelogarithmic law and suggested that the early stages of a relaminarization process weretaking place.Numerical simulations of boundary layers over hills and bumps have also been per-formed, however most employed LES due to the large computational cost required bythese types of flows. Wu & Squires (1998) simulated the experimental bump flow ofWebster et al. (1996) with wall-resolved LES at a Reynolds number almost three timeslower, and confirmed the significant departure above the logarithmic law in the FPGregion. This geometry was also studied by Cavar & Meyer (2011) and Matai & Durbin(2019) with LES producing similar results. Matai & Durbin (2019) expanded the studyby considering a family of bumps of increasing height. For all cases, a large departureof the velocity above the logarithmic law was observed in the FPG as well as a plateauor a rapid oscillation in the skin friction at the start of the APG. They showed that thenon-dimensional favorable pressure gradient ∆ p = νρu τ ∂p∂s (1.1)exceeded the value of -0.018 identified by Patel (1965) as the start of a relaminarizationprocess. Note that in (1.1), ν is the kinematic viscosity, ρ is the fluid density, p is the staticpressure at the wall, u τ is the friction velocity, and s is the streamwise direction. Uzun& Malik (2018) performed a wall-resolved LES study of the NASA wall-mounted hump NS of a Boundary Layer with Strong Pressure Gradients et al. (2006). Bothnumerical and experimental data showed the presence of a plateau in the skin frictioncoefficient profile on the upstream side of the bump where the FPG was strong. It wasnoted that in this region the boundary layer was undergoing a relaminarization processaccording to the acceleration parameter based on the edge velocity U e K = νU e ∂U e ∂s (1.2)surpassing the limit of 3 × − (Narasimha & Sreenivasan 1979), but full relaminarizationwas not achieved.Although often not the focus of the aforementioned studies on bump and hill flows, allreported a region of strong FPG with large deviations from standard turbulent boundarylayer behavior. Some of them even mentioned signs of relaminarization (also known asreverse transition). Patel (1965) was among the first to report on these effects, noting adeviation of the streamwise velocity profile above the logarithmic law for large enoughnegative values of the pressure gradient parameter ∆ p . They supposed the breakdownof the law was due to the process of reversion to laminar flow and proposed a tentativecritical value of ∆ p = − . ∆ p did not describe the near-wall flow completely. In a later study, Patel & Head (1968) expanded the initial workand proposed the use of a different quantity to measure the effects of a favorable pressuregradient on the inner layer. They argued that the non-dimensional shear stress gradient ∆ τ = ναρu τ (1.3)is more universal and identified a critical value of − .
009 for the beginning of thedeparture above the logarithmic law and reverse transition. Note that in (1.3), α isthe gradient of the total shear stress in the wall-normal direction n across the viscoussublayer such that τ = τ w + αn . They further argued that the acceleration parameter K is not an adequate quantity since it is not able to describe the near-wall flow directly,and based on the analysis on ∆ τ , they proposed a critical value on the pressure gradientof ∆ p = − . ∆ τ to − . ∆ p = − .
025 and K = 3 × − as further indication of relaminarization, as well as alocal minimum in the shape factor. These critical values on the pressure gradient weresupported by a DNS of sink flows (Spalart 1986). Warnack & Fernholz (1998) analyzed theReynolds stresses of a boundary layer during reverse transition. They observed a decreaseof the shear and streamwise stresses when normalized by the local friction velocity thatpersisted until the local minimum in the Reynolds number Re θ , and then a sharp rise dueto retransition to turbulence. Similar behavior was shown for the shear and streamwisestress production when normalized by inner units.Despite the APG region and separation receiving most of the attention in the literature,particularly in the context of improving RANS closures, it will become apparent through-out this paper that a complete understanding and accurate modeling of the FPG region,which may include a relaminarization process, are critical for the prediction of certainaerodynamic flows such as the one studied here. However, current RANS and near-wallmodels for LES perform adequately in mild FPG, but often significantly overpredict the R. Balin and K. E. Jansen skin friction and shear stress in strong FPG (Matai & Durbin 2019; Balin et al. ≈
20 boundary layer thicknesses δ ) causesthe onset of relaminarization and a significant weakening of the near-wall turbulence.The boundary layer, however, does not relaminarize completely and stays intermittentlyturbulent. At the peak of the bump, the weakened near-wall turbulence experiences asudden enhancement in intensity due to partial retransition, which in turn leads to anatypical skin friction response and a more resilient boundary layer. The DNS was designedto focus on the part of the flow leading up to incipient separation rather than downstreamof it, with additional emphasis on the FPG region, and the discussion of results reflectsthis choice.It must be mentioned that the present DNS is part of a larger joint study withother research groups and a recent publication on this flow already exists. However,key differences are present between this paper and Uzun & Malik (2020), making themdistinct reports of the flow. In particular, this paper focuses on the strong FPG effects anddiscusses the flow physics involved in light of the body of literature on relaminarizationand strongly accelerated flows, while Uzun & Malik (2020) takes a more holistic approachwith a greater emphasis on the incipient separation and subsequent recovery.This paper is organized as follows. § § §
2. Numerical Setup
Problem Definition
The flow computed in this study is the turbulent boundary layer over the prismaticextrusion of a two-dimensional (2D) Gaussian shaped bump. The surface is defined by(2.1), which depends on the height parameter h and the length parameter x , and isshown by the black curve on the lower surface of the domain in figure 1. y ( x ) = h exp (cid:16) − (cid:0) x/x (cid:1) (cid:17) (2.1)Note that the x coordinate is aligned with the freestream flow far upstream of the bump,the y coordinate is vertical and normal to the freestream, and the z coordinate is alignedwith the spanwise direction. Moreover, (2.1) defines the entire lower surface of the domain,meaning that there is no flat-plate region on either side of the bump and the curvatureis everywhere continuous. It is important to mention that the geometry selected for this NS of a Boundary Layer with Strong Pressure Gradients Figure 1.
The black curves outline the domain of the bump flow, while the green curves showthe boundary layer thickness on both no-slip walls predicted by preliminary RANS. The dottedvertical line marks the location of the inflow to the DNS.
DNS is exactly the centerline of a three-dimensional (3D) bump developed at The BoeingCompany (Slotnik 2019) and studied experimentally at the University of Washington(Williams et al. h/L = 0 .
085 and x /L = 0 . L = 3 ft is thelength of the square cross-section of the wind tunnel used for the 3D bump experiments.These values were selected in order to obtain the desired pressure gradients and separationusing preliminary RANS simulations of the 3D bump (Slotnik 2019). The flow studiedhas a Reynolds number of Re L = 1 . × , corresponding to Re h = 85 ,
000 whenmeasured against the bump height. The freestream velocity is U ∞ = 16 .
40 m/s, which atstandard sea level conditions results in a small enough Mach number ( M ∞ = 0 . Re θ = 1 , / et al. et al. y/L = 0 .
5. The streamwise location of theorigins of the top and bottom boundary layers is at the leading edges of the no-slip wallsat x/L = − et al. x direction equal to U ∞ to be correctly imposed at the inflow ( x/L ≈ − . Solution Approach
While figure 1 describes the entire RANS flow domain, only a fraction of it wasincluded in the DNS. The inflow was moved downstream to x/L = − . R. Balin and K. E. Jansen & Shur 1997; Shur et al. et al. (2018) was applied to correct the underpredictionof the skin friction coefficient by the SA model at low Re θ . The spanwise period ofthe DNS domain was set to 4 . δ in , where δ in is the inflow 99.5% boundary layerthickness. Due to the significant growth of the boundary layer on the downstream sideof the bump as shown in figure 1, this period may introduce some confinement effects inthis part of the solution (Coleman et al. et al. (2020), wherein it was showedthat, while the near-wall region is poorly predicted by the wall model, the outer layerReynolds stresses and boundary layer thickness in the region of interest of the floware captured accurately. Since the spanwise period affects the largest scales of thedomain, which are well captured in both simulations, this exercise is representative ofthe confinement effects experienced by the DNS. No significant differences were observedin the flow and turbulent quantities of interest (skin friction, pressure gradient, velocity,and Reynolds stresses) obtained with the two WMLES. Consequently, the solution in theregion of interest of this DNS is considered to be free of any confinement effects from thespanwise period chosen.The boundary conditions enforced in the DNS were as follows. The bump surfacewas treated as a no-slip wall. The top surface was modeled as an inviscid wall offsetby the RANS predicted displacement thickness described above with zero transpiration(zero velocity component normal to the surface) and zero traction. At the outflow, weakenforcement of zero pressure was applied along with zero traction. Effects from thisboundary condition on the interior domain are contained within a streamwise distanceof one local boundary layer thickness and thus did not affect the upstream solution. Atthe inflow, the synthetic turbulence generator (STG) of Shur et al. (2014) was selectedto introduce unsteady flow into the domain, which has been shown to produce realisticturbulence a short distance downstream of the inlet for both wall-modeled LES (Shur et al. et al. x/L = − .
8. The mean Reynolds stress and velocity profiles requiredby the STG method were extracted from an additional RANS simulation at the inflowlocation of x/L = − .
6. This simulation used the same turbulence model (SARC-lowRe)and setup described above with only the following exception to the domain describedby figure 1. The top boundary was no longer flat and modeled with a no-slip condition,but instead it was slanted and modeled as an inviscid wall to match the DNS domain.This was done intentionally to extract profiles along the entire height of the inflow planeconsistent with the DNS domain and boundary conditions. Moreover, it was used to verifythat the constriction effects were appropriately captured by the slanted upper surface,and that was indeed the case.The computational grid used for the DNS of the bump was structured with a total of554 million points. It was designed with spacing ∆s + = 15, ∆z +max < ∆n +1 = 0 .
1, and
NS of a Boundary Layer with Strong Pressure Gradients Figure 2.
Pressure coefficient (a) and skin friction coefficient (b) on the surface of theGaussian bump. ∆n +max <
10, where ( s, n, z ) is the bump aligned coordinate system ( s and n are tangentand normal to the bump surface, respectively). Moreover, the () + superscript signifiesscaling by wall units using the friction velocity u τ and the viscous length scale l ν = ν/u τ .Note that since the DNS solution was not available at the time of the grid generation,the u τ profile from the preliminary RANS was used (see figure 2 for differences in wallfriction between RANS and DNS). The growth factor in the wall-normal direction waslimited to 5%. The streamwise spacing ∆s + = 15 was achieved everywhere by using thelocal value of the RANS friction velocity, however the growth or decay across adjacentelements was limited to 1% resulting in a smooth streamwise variation of ∆s . Closeto the separation region predicted by RANS, where u τ becomes ill-defined, an effectivefriction velocity was used instead to maintain adequate spacing in physical units. Thewall spacing ∆n +1 = 0 . ∆z + was fixed using the maximum friction velocityover the bump surface, which occurs near the bump peak. Since the peak skin frictionis overpredicted by RANS (see figure 2), the maximum spanwise spacing in the DNS isactually ∆z +max = 8. Similarly, the wall-normal spacing ∆n +max <
10 was driven by themaximum boundary layer thickness also computed by the preliminary RANS, which asshown in figure 1 is found at the outflow of the domain. As a result, spanwise and wall-normal spacing finer than 8 and 10 plus units, respectively, are present over the regionsdiscussed in this paper. For instance, ∆z + ≈ ∆n +max ≈ ∆n +max < et al. T long enough to satisfy T > t eddy everywherein the region of interest of this bump flow, where the eddy-turnover time is definedbased on the local edge velocity U e and boundary layer thickness δ as t eddy = U e /δ .Stationarity of the time- and span-averaged statistics was verified by comparing flowand turbulence quantities obtained from different time windows within the time interval T . The non-dimensional time step size was ∆t + = 0 .
11 based on the maximum averagefriction velocity of the DNS, which ensured a maximum Courant–Friedrichs–Lewy (CFL)number below one at each time step.
R. Balin and K. E. Jansen
All simulations presented in this work were performed with a stabilized finite elementmethod (Whiting & Jansen 1999) using tri-linear hexahedral elements and second orderaccurate, fully implicit time integration (Jansen et al. et al. (2009), in which tri-linear hexahedral elements were also used. Stabilization and time integration parameters(which affect numerical dissipation) chosen for this DNS follow the work of Trofimova et al. (2009) and were verified with a flat plate DNS to reproduce skin friction, velocity,and Reynolds stresses from a number of previous studies (Coles 1962; Schlatter et al. et al. et al. et al.
3. Results and Discussion
Wall and Integral Quantities
Time- and span-averaged pressure and skin friction coefficient profiles on the surfaceof the Gaussian bump obtained from the DNS are presented in figure 2. The coefficientsare defined in (3.1), where p ref is the reference wall pressure at location x/L = − . τ w is the wall shear stress rotated to the curvilinear coordinates ( s, n, z ). C p = p w − p ref12 ρ ∞ U ∞ C f = τ w ρ ∞ U ∞ (3.1)In this figure, the solution from the preliminary RANS simulation used to obtain theinflow profiles is also shown for comparison.Very good agreement is obtained for the pressure coefficient, indicating the successfulchoice of the DNS sub-domain and boundary conditions. The C p profile also outlinesthe series of pressure gradients experienced by the boundary layer. At the inflow to theDNS ( x/L = − . x/L = − .
29 where it switches to a strongFPG accelerating the flow over the upstream side of the bump. At the bump peak, arapid change from strong favorable to strong adverse occurs, with the latter persistinguntil about x/L = 0 .
40. Finally, a mild FPG helps the boundary layer recover afterincipient separation.Significant differences are observed for the skin friction coefficient predictions infigure 2. Fair agreement is only obtained in the initial mild APG and soon after the startof the strong FPG the curves deviate with RANS largely overpredicting C f over thebump. The DNS solution exhibits a much smaller peak and a local minimum-maximumimmediately downstream of the switch from FPG to APG. This feature of the wallshear stress is also present in the DNS of Uzun & Malik (2020) and is similar to the onesdocumented for other bump flows with strong FPG (Matai & Durbin 2019; Narasimha &Sreenivasan 1979; Warnack & Fernholz 1998). Of interest is also the streamwise positionof the skin friction maximum, which in both cases is located upstream of the bumppeak, but occurs further upstream in the DNS. Moreover, the DNS exhibits a muchsteeper decrease in wall shear between the location of C f maximum and x/L = 0 . C f approaching zero ( C f < × − over 0 . < x/L < .
27) but only becoming negativeover a short distance. Instantaneously, however, the flow at the wall does reverse directionand small confined separation bubbles are present.
NS of a Boundary Layer with Strong Pressure Gradients Figure 3.
Contours of instantaneous vorticity magnitude on the surface of the Gaussian bump.
Contours of the instantaneous vorticity magnitude on the surface of the bump infigure 3 elucidate some of the features of the DNS skin friction coefficient profile. Onthe upstream side of the bump, the footprint of typical near-wall structures is seen asstreamwise elongated streaks of high and low vorticity. As the boundary layer progressesthrough the strong FPG, the acceleration of the flow causes the wall vorticity to riseand the streaks to grow significantly both in width and length. Towards the end of theFPG, around x/L = − .
03 the wall vorticity drops and the streaks become weaker. Quietregions of low vorticity also form. This behavior explains the steep drop in C f between themaximum value of the curve and the bump peak. Starting slightly upstream of x/L = 0 . x/L = 0 .
05, spots of large vorticity arefound intermittent with the quiet regions, which visually resemble those characteristicof laminar-to-turbulent transition. Note that x/L = 0 .
05 is the location of the smalllocal maximum in the skin friction profile. These spots then appear to culminate in aregion of intense turbulent activity (0 . (cid:54) x/L (cid:54) .
10) with many small-scale and fairlyisotropic structures of large vorticity. Continuing further downstream, the intensity ofthe small-scale structures decreases as the skin friction also drops. Note that, due to thestrong APG which quickly brings the flow to incipient separation, the canonical streaksdo not reappear on the downstream side of the bump and by x/L = 0 .
15 the turbulentstructure at the wall changes dramatically once again.Due to the strong pressure gradients and geometry of the bump, the freestream(irrotational) flow is distorted and highly non-uniform ( ∂u s /∂n (cid:54) = 0). As a result, theclassical definitions of the boundary layer thickness δ and integral quantities suchas the displacement and momentum thicknesses, δ ∗ and θ respectively, are no longerapplicable. To resolve this issue, the definitions based on the generalized velocity (or“pseudo-velocity”), ˜ U , of Spalart & Watmuff (1993) were used instead. These are repeatedin (3.2)-(3.4) from Coleman et al. (2018) for convenience, where ω z is the mean spanwisevorticity. Note that the bump-aligned curvilinear coordinate system ( s, n, z ) is used forthe definitions instead of the freestream-aligned Cartesian system ( x, y, z ).˜ U ( s, n ) ≡ − (cid:90) n ω z ( s, n (cid:48) ) dn (cid:48) (3.2)˜ δ ∗ ( s ) ≡ − U e ( s ) (cid:90) ∞ nω z ( s, n ) dn (3.3)0 R. Balin and K. E. Jansen ˜ θ ( s ) ≡ −
2( ˜ U e ( s )) (cid:90) ∞ n ˜ U ( s, n ) ω z ( s, n ) dn − ˜ δ ∗ ( s ) (3.4)The edge velocity in the context of ˜ U is defined as˜ U e ( s ) ≡ ˜ U ( s, n → ∞ ) , (3.5)from which the 99 .
5% boundary layer thickness, ˜ δ , is computed as the height abovethe wall where ˜ U = 0 .
995 ˜ U e . This approach has been used successfully in other studies ofpressure gradient flows (Coleman et al. n → ∞ are evaluated correctly. Of course, the integralsmust be carried out over a finite height and it is thus important to ensure that the desiredquantities are sufficiently converged under changes to the truncation. Integration to 1 . . δ showed negligible difference in the integral quantities, and thus thelatter was used.Figure 4 shows the variation of the boundary layer integral quantities defined aboveover the bump surface. The streamwise extent of the domain is limited to − . (cid:54) x/L (cid:54) . L was used.The 99 .
5% boundary layer thickness generally follows the trends set by the pressuregradients, with a slight delay in its response. In the initial mild APG, ˜ δ grows almostby a factor of two until x/L = − .
25 which is downstream of the start of the FPG at x/L = − .
29. It then shrinks as it progresses through the rest of the FPG reaching alocal minimum just upstream of the bump peak. The strong APG on the downstreamside of the bump causes it to grow rapidly once again creating another local maximumat x/L = 0 .
30 downstream of incipient separation where ˜ δ is almost a factor of fourlarger than the inflow value. The displacement thickness follows similar trends, howeverit responds more quickly to changes in pressure gradient. The first local maximum islocated at x/L = − .
29, immediately downstream of the first change in sign of thepressure gradient. It is then reduced significantly in the strong FPG, dropping below theinflow value, and reaching a local minimum close to the bump peak at x/L = − . δ ∗ grows very rapidly to a local maximum which is almost 10 times largerthan the inflow value and is also located downstream of incipient separation. The thirdmeasure of thickness, ˜ θ , behaves similarly to the other two. The first local maximum islocated slightly downstream of the start of the FPG at x/L = − .
27, suggesting a slowerresponse to the pressure gradient relative to ˜ δ ∗ . The local minimum is immediatelyupstream of the bump peak with a value below that of the inflow. In the strong APG, ˜ θ also increases, although not as significantly as the displacement thickness.Figure 4(b) shows the shape factor H = ˜ δ ∗ / ˜ θ . After the inflow, H settles at a value of1 .
44 and then grows slightly to a local maximum of 1 .
48 as the boundary layer advectsthrough the mild APG. These are both typical values for turbulent boundary layers.The local maximum in this case is located at x/L = − .
35, which is well upstream ofthe switch from APG to FPG. The shape factor is then seen to drop through the startand middle sections of the FPG until a local minimum with value 1 .
27. Note that thestreamwise location of this extremum is at x/L = − .
11, which is significantly upstreamof the bump peak. During the remainder of the FPG, H grows and forms a plateaujust downstream of the bump peak where the oscillation in C f is also observed. Finallyit grows significantly in the strong APG on the downstream side of the bump untilanother local maximum located inside the region of incipient separation. After incipientseparation, H drops rapidly as the new boundary layer develops.The momentum thickness Reynolds number, also shown in figure 4(b), follows the NS of a Boundary Layer with Strong Pressure Gradients Figure 4.
Variation of different boundary layer thickness measures (a) and shape factor H and momentum thickness Reynolds number Re ˜ θ (b) over the bump. shape of the momentum thickness, however it is interesting to note the values of theReynolds number in the different regions of the flow. Initially Re ˜ θ grows from the inflowvalue of 1 ,
050 to 2 ,
200 at the same location of the local maximum in ˜ θ . Just upstreamof the bump peak, the minimum value of 600 is achieved, which is well above the lowestvalue identified for turbulence to exist. Downstream of incipient separation, the maximumvalue is around 5 , κ = κ ˜ δ reaches a maximum valueof 0 .
033 at x/L = − .
25. Note that κ in this context is the surface curvature and notthe streamline curvature since the latter varies with distance from the wall. At the peakof the bump, where the curvature is convex, a maximum of 0 .
055 is observed. Finally,in the concave region downstream of the bump ˆ κ = 0 .
065 at x/L = 0 .
26 due to theboundary layer thickness growing significantly. As discussed in Baskaran et al. (1991);Narasimha & Sreenivasan (1979); Schwarz & Plesniak (1996); So & Mellor (1973, 1975),streamline curvature effects are mainly observed in the outer region of the boundarylayer while the inner region is not directly affected unless the non-dimensional curvature κ + = κν/u τ is large. This parameter reaches an absolute maximum of 7 . × − at thepeak of the bump, indicating that direct curvature effects on the near-wall turbulencecan be considered negligible throughout the flow leaving pressure gradient effects to bedominant. 3.2. Mild Adverse Pressure Gradient
The initial disturbance experienced by the boundary layer is due to a weak to mildadverse pressure gradient. This force gradually increases from the inflow to a peak at x/L = − .
38 with strength K = − . × − (see figure 6), and then terminates at x/L = − .
29 where the pressure gradient changes sign. Concave streamline curvatureeffects are also present, gradually increasing from negligible values at the inflow to ˆ κ =0 .
028 at the end of the APG. Both APG and concave curvature effects are destabilizingto the boundary layer turbulence, increasing the intensity of the turbulent fluctuationsand the Reynolds shear stress.Changes to the streamwise velocity, turbulent kinetic energy (TKE), and Reynoldsshear stress due to the mild APG are shown in figure 5. The TKE is defined as k = (cid:16) u (cid:48) s + u (cid:48) n + u (cid:48) z (cid:17) . Further, note that at x/L = − .
38, ˆ κ = 0 .
01. The velocity exhibitsan increase in the wake region and a slight shortening of the logarithmic region which also2
R. Balin and K. E. Jansen
Figure 5.
Profiles of the streamwise velocity (a), turbulent kinetic energy (b), and Reynoldsshear stress (c) in the initial mild APG region. appears to have a slightly lower value of the intercept. Since concave curvature causes areduction in the wake (So & Mellor 1975), these affects are attributed to the APG and theincrease in Re ˜ θ . The TKE profile non-dimensionalized by the local friction velocity showsan increase mainly between 0 . (cid:54) n/ ˜ δ (cid:54) .
8, whereas the peak value increases onlymarginally. The non-dimensional Reynolds shear stress is also increased by the combinedAPG and curvature effects, but in this case the effects are evident everywhere below n/ ˜ δ < . Strong Favorable Pressure Gradient
This section focuses on the favorable pressure gradient region of the flow and describesthe effects of this force on the mean flow and turbulence. However, before presentingthe results of this study, it is insightful to review the full spectrum of FPG affects onturbulent boundary layers which have been studied in the literature. Not all the casesdiscussed are present in this flow, but such a discussion grounds the observations thatfollow in the available body of work. For mild streamwise gradients, these effects are
NS of a Boundary Layer with Strong Pressure Gradients et al. H = δ ∗ /θ . Overall, an FPG has a stabilizing effect on theboundary layer turbulence. Furthermore, the streamwise velocity may still show goodagreement with the standard (or a slightly displaced) logarithmic law (Smits et al. .
5% boundary layer thickness, anda reduction in the rate of wall-layer bursting.Under further sustained acceleration, the boundary layer enters a relaminarization(or reverse transitional) process. This is the second stage, and lasts until the processis complete resulting in a quasi-laminar boundary layer. It is important to highlightthat this is a gradual process, thus the onset does not imply that a quasi-laminar state isachieved. During relaminarization, the previously fully turbulent boundary layer developsa viscous dominated inner layer stabilized by the acceleration, while the turbulence inthe outer layer is distorted. Moreover, the skin friction is observed to drop while theflow still accelerates, the shape factor increases, and the relative (non-dimensionalizedby local u τ ) Reynolds stresses and TKE production drop significantly.The final stage occurs when the FPG is relaxed. The acceleration can no longer stabilizethe near-wall flow, and soon after the onset of instability a retransition process originatingnear the wall returns the boundary layer towards a fully turbulent state. This processis rapid and shares some similarities to laminar-to-turbulent transition, such as theformation and growth of turbulent spots (Blackwelder & Kovasznay 1972). Other notablefeatures of this process are a sudden rise in skin friction, a significant enhancement ofthe turbulent intensities, and a local maximum of the shape factor.Since relaminarization is a process that brings about significant changes to the bound-ary layer, it is of great interest to identify the onset and completion points. The lattercan be easily defined as the location where the effects of the Reynolds stresses on themean flow dynamics are negligible. Fluctuations are still present as a remnant from theupstream flow, but do not contribute to the development of the mean velocity. Thenear-wall bursting process has also ceased, therefore eliminating production of turbulentkinetic energy. The onset point, however, is more complicated to define and identify. Someprevious studies (Sreenivasan 1982; Warnack & Fernholz 1998), in fact, argued that it is4 R. Balin and K. E. Jansen an ill-defined exercise to determine a parameter and its critical value that uniquely definethe onset because the process is not catastrophic and sudden but rather it is gradual.Nevertheless, all the non-dimensional parameters and methods proposed in the literaturecan be used as guidelines to suggest whether or not relaminarization has ensued.Mean streamline curvature effects are also discussed in this section. The curvatureof the Gaussian bump changes sign from concave to convex at x/L = − .
14, which isapproximately in the middle of the FPG region. Moreover, the strongest concave effectsas measured by ˆ κ are also located within this region at x/L = − .
25. As a result, bothdestabilizing and stabilizing effects are present from this force and are interacting with theFPG. Concave curvature is known to increase the Reynolds shear stress and wall-normalturbulent transport of momentum, as well as reducing the wake of the streamwise velocity(So & Mellor 1975). By contrast, convex curvature can cause the Reynolds shear stress tovanish or change sign and the turbulent kinetic energy production to “turn off” furtherfrom the wall (So & Mellor 1973). Additionally, it reduces the extent of the logarithmiclaw and increases the wake. Decoupling of the inner and outer layers has been observedfor boundary layers under streamline curvature of both kind, with the effects mainlyfocused in the outer layer and negligible impact in the near wall region, although theskin friction is subject to change (convex curvature reduces C f and vice versa) (Baskaran et al. K defined in (1.2) is shown. Since this quantity is proportional to the gradientof the edge velocity ˜ U e , a negative value corresponds to a deceleration and thus an APG,while a positive value corresponds to an acceleration and thus an FPG. At the inlet, K isa small negative value, indicating the presence of a very weak APG at this location. Thisis followed by a mild APG ( K > − × − ) until x/L = − .
29 and then a strong FPGover the upstream side of the bump which approaches, but does not cross, the limit of K = 3 . × − marking the start of a reverse transitional process as outlined in previousstudies (Narasimha & Sreenivasan 1979; Sreenivasan 1982; Spalart 1986). The maximumvalue is in fact 2 . × − at x/L = − .
13. At the bump peak, the pressure gradientbecomes adverse again, reaching a maximum strength in terms of the accelerationparameter of − . × − . According to this quantity, therefore, onset of relaminarizationis not achieved, even considering the range 2 . × − (cid:54) K (cid:54) . × − identifiedwith DNS of a series of sink flows (Spalart 1986). However, one must point out that K has been criticized in previous studies (Patel & Head 1968; Sreenivasan 1982) for its lackof knowledge of the near-wall physics.The non-dimensional pressure gradient Λ = − ∂p w ∂s ˜ δ τ w (3.6)introduced by Narasimha & Sreenivasan (1973) is also shown in figure 6(a). This pa-rameter has been used to measure where pressure gradient forces dominate over the NS of a Boundary Layer with Strong Pressure Gradients Figure 6.
Pressure gradient parameters K and Λ (a) and pressure gradient ∆ p and shearstress gradient ∆ τ non-dimensionalized by inner units (b) over the bump. shear forces, thus useful in identifying the onset and completion of relaminarization(Narasimha & Sreenivasan 1973, 1979; Sreenivasan 1982). Values above 50 have beenproposed as a sign of the completion of the process and the achievement of a quasi-laminar state. Furthermore, based on the data presented in these studies, values of Λ between 10 and 25 can be correlated to the departure from fully turbulent flow and thusthe onset of relaminarization. The maximum value achieved in the FPG of the bump flowis around 19, indicating that this boundary layer is far from reaching the completion ofrelaminarization, but the process is likely underway.Other measures of the pressure gradient and its effects on the near-wall flow areshown in figure 6(b). These are ∆ p and ∆ τ defined in (1.1) and (1.3), respectively. Notsurprisingly, the two quantities are similar, with a small offset in the mild APG regionand then good agreement in the FPG. Note that due to the non-dimensionalizationby the friction velocity u τ , these two quantities become ill-defined in the vicinity ofincipient separation at x/L = 0 .
19. In the figure are also shown two horizontal lines,one at ∆ τ = − .
013 to mark the start of the deviation above the logarithmic law(discussed in detail later in this section) (Patel & Head 1968; Bradshaw 1969) andone at ∆ p = − .
025 to mark the onset of relaminarization (Narasimha & Sreenivasan1979; Spalart 1986). The non-dimensional shear stress gradient profile for the bumpcrosses the critical value of − .
013 at x/L = − .
22 relatively soon after the start of theFPG, but neither curve crosses the limit for relaminarization. The minimum value of ∆ p is − .
022 at x/L = − .
11. Similarly to the profile for the acceleration parameter K ,these parameters suggest that the favorable pressure gradient is indeed strong, enoughto expect the breakdown of the logarithmic law, but not enough for a relaminarizationprocess to start. It is worth mentioning, however, that Narayanan & Ramjee (1969)reported ∆ p = − .
020 as a critical value, which would imply that the bump boundarylayer does enter a relaminarization process at x/L = − .
17. The initial study by Patel(1965) even suggested ∆ p = − . − .
024 (Patel & Head 1968). This variation of the critical value for ∆ p in the literature highlights the uncertainty associated with marking the onset ofrelaminarization for general flows with parameters based on the pressure gradient.Integral quantities have also been investigated for strongly accelerated flows in con-nection to the onset of relaminarization (Narasimha & Sreenivasan 1973; Warnack &Fernholz 1998; Blackwelder & Kovasznay 1972; Cal & Castillo 2008). Combining theknown observations that the shape factor H decreases for turbulent boundary layers6 R. Balin and K. E. Jansen under FPG and that laminar boundary layers exhibit much larger values of H , theappearance of a local minimum in this quantity followed by a sharp rise while the flowis still being accelerated has been used to determine when the boundary layer departsfrom a fully turbulent state to approach a quasi-laminar one. Minimum values have beenobserved around 1 . H ∗ = ˜ δ ∗ / ˜ δ initially decreases in an FPG but then suddenly growsat the onset of relaminarization, therefore making a local minimum of H ∗ a marker ofreverse transition. For the Gaussian bump, the local minimum in H is co-located withthe peak of the FPG as measured by ∆ p and does not grow rapidly downstream of thispoint. Moreover, the minimum value of H is larger than what was reported in previousstudies. The local minimum of H ∗ is found just upstream of the bump peak where ˜ δ ∗ isalso a minimum and the FPG is very weak. Consequently, these integral parameters donot appear to be insightful in describing the onset of relaminarization for this flow. Theuncertainties associated with using H for this purpose were also discussed in Sreenivasan(1982).Another quantity not directly dependent on the pressure gradient that may be usefulis the skin friction coefficient. A sudden decrease in the wall shear stress producing alocal maximum while the flow is still undergoing acceleration has been identified as asign of the the relaminarization process taking effect in the near-wall region (Warnack &Fernholz 1998; Narasimha & Sreenivasan 1973). The reasoning being that as the near-wallregion progresses from a fully turbulent to a quasi-laminar state, the skin friction dropsconsistently with the known properties of laminar boundary layers. As shown in figure 2,the local maximum is located upstream of the bump peak at x/L = − .
055 where theflow is still experiencing a fairly strong acceleration ( K = 1 . × − ). This result is incontrast to what is expected of a fully turbulent boundary layer under the same pressuregradient, which resembles more closely the RANS prediction with a skin friction peakmuch closer to the end of the acceleration. Therefore, figure 2 suggests that the near-wallregion of the boundary layer over the bump is indeed experiencing relaminarization.Further credence to this claim comes by looking at instantaneous contours of vorticityat different heights within the boundary layer. Figure 7 shows slices of the flow at a fewlocations of constant n/ ˜ δ . Note that because ˜ δ varies with downstream position,these surfaces track local boundary layer height, not distance to wall. The vorticitymagnitude is normalized by the local (same streamwise location) time- and spanwise-averaged wall vorticity magnitude ω w in order to remove the increase of vorticity due toflow acceleration and highlight the fluctuations relative to their local mean wall value.Moreover, solid black lines across the width are placed to mark important landmarkswithin the FPG region. These are the start of the FPG at x/L = − .
29, the locationwhere the streamwise velocity departs above the logarithmic law at x/L = − . x/L = − .
11, andthe bump peak and end of the FPG at x/L = 0 .
00. Note that the same figures were alsoanalyzed at different time steps in order to confirm that the following features are alwayspresent and not just at the time instant shown here.Very close to the wall, the typical streamwise streaks of high and low vorticity arevisible at the beginning of the FPG region. These remain fairly unchanged until around x/L = − .
22, downstream of which the streaks appear to stretch in the streamwisedirection and reduce in relative intensity (recall the contours are normalized by thelocal mean wall vorticity). As the location of peak FPG strength is approached, thecharacter of the turbulence is significantly altered. The streaks become very elongatedin the streamwise direction while also thickening in the spanwise direction, and theirrelative intensity is significantly reduced. These changes continue past the peak FPG
NS of a Boundary Layer with Strong Pressure Gradients § n/ ˜ δ = 0 .
05 and 0 .
1, similar trends can be observed.Fully turbulent flow is present at the start of the FPG and is maintained until slightlydownstream of x/L = − .
22, after which the turbulent scales weaken relative to thewall vorticity, stretch significantly in the streamwise direction, and quiet regions of lowvorticity appear. At this height above the wall, the streaks of high vorticity fluctuationsturn into very thin and long ridges, while the quiet areas form “valleys” that are muchwider but equally as long. These are clear signs that the near-wall region of the boundarylayer is no longer fully turbulent, and instead becomes intermittent and approaches aquasi-laminar state. A relaminarization process is therefore taking place in this flow.The contours also show how this truly is a gradual process since no clear demarcationcan be identified from these instantaneous fields across which the turbulence suddenlychanges from fully turbulent to intermittent. Moreover, for this Gaussian bump the onsetof relaminarization occurs upstream of the peak strength of the FPG, where the pressuregradient parameters K and ∆ p are even further from the critical values reported in otherstudies (see figure 6).In the outer layer at n/ ˜ δ = 0 .
4, the turbulence is also very much affected by thestrong FPG. As the flow is being accelerated, the normalized vorticity fluctuations weakenwith the smaller scales decaying entirely while the larger scales remain, which is consistentwith the decrease in Reynolds number Re ˜ θ throughout the FPG. Contrarily to the near-wall region, however, the shape of the fluctuations does not appear to be significantlyimpacted.As was noted earlier, both concave and convex streamline curvature effects are actingduring the FPG, with the change from the former to the latter being located at x/L = − .
14. The vorticity contours in figure 7 clearly show that changes to the turbulence aretaking place where the destabilizing concave curvature is acting, well before the stabilizingconvex curvature takes effect. The pressure gradient effects are therefore dominating overthe curvature during most of FPG and the onset of relaminarization is attributed to theacceleration alone, rather than to the combination of the two forces.Finally, the onset of relaminarization can be implied by evidence of retransition (ora return to fully turbulent flow) where the strong FPG is relaxed and the APG starts.The details of this process are discussed later in § R. Balin and K. E. Jansen
Figure 7.
Instantaneous vorticity magnitude normalized by the local time- andspanwise-averaged wall vorticity at different locations within the boundary layer in the FPGregion of the bump. From top to bottom, the heights above the wall of the slices are n/ ˜ δ = 0 . , . , . , . n + ≈ , , , monotonically accelerated from a zero pressure gradient by variation of the top boundarycondition. The boundary layer over the Gaussian bump is deeply different; it experiencesan APG upstream of the acceleration which, while being mild, affects the velocity andstress profiles in a significant manner, there is a high degree of non-equilibrium inthe sense of the Clauser pressure gradient parameter, streamline curvature effects arepresent throughout the entire flow, and the Reynolds number is lower than what most NS of a Boundary Layer with Strong Pressure Gradients Figure 8.
Mean streamwise velocity profiles in the FPG region normalized by wall units. experimental studies were able to achieve. Furthermore, given the results discussed above,a definite marker or critical value marking the onset of relaminarization was not found.This is not a surprising conclusion, and in fact agrees with the interpretation offered bySreenivasan (1982) and seconded by Warnack & Fernholz (1998) that it is an ill-definedexercise to do so since relaminarization is a gradual process and not a catastrophic event.By extension, this study provides further evidence that critical parameters of any specificquantity should only be used as guidelines to suggest that reverse transition might beoccurring, but further analysis of the state of the turbulence is needed for a more definitestatement to be made. Moreover, a turbulence model based on these parameters is likelyto not be universal.Figure 8(a) shows the non-dimensional streamwise velocity profile at a number oflocations on the upstream side of the bump. At the start of the FPG ( x/L = − . κ = 0 .
41 and B = 5 .
0) but a logarithmic section remains and the wake is increased. Once again,the boundary layer entering the strong FPG is not the canonical ZPG flow. Soon afterthat at x/L = − .
22, the FPG effects become clear. The wake is significantly reducedand the profile lacks a linear region that would result from a logarithmic relationshipof the form u + s = ln ( n + ) /κ + B . There is, therefore, the breakdown of the logarithmiclaw with the velocity remaining for the most part below the law. As the flow progressesthrough the FPG, the effects become increasingly stronger. The velocity continues to riseabove the log law until the peak of the bump, with a thickening of the viscous sublayerand a continuous reduction of the wake in the outer layer. Note that convex curvatureeffects are also present downstream of x/L = − .
14, however the continued reduction ofthe wake downstream of this location indicate that the pressure gradient is dominatingover curvature on the velocity. This result is consistent with the experimental study ofSchwarz & Plesniak (1996).Additionally, figure 8(b) shows the velocity profiles at three locations in the proximityof the critical parameter ∆ τ = − .
013 where the initial departure above the logarithmiclaw has been identified (Bradshaw 1969). From figure 6, this critical parameter is observedat x/L = − .
22. Slightly upstream of this location, the velocity is entirely below thelog law in the range 10 (cid:54) n + (cid:54) (cid:54) n + (cid:54)
60. At the location of ∆ τ = − . n + = 35and is below otherwise. These results are in great agreement with Bradshaw (1969) andthe notion of Patel & Head (1968) that the non-dimensional shear stress gradient is0 R. Balin and K. E. Jansen a suitable quantity to measure the departure from the logarithmic law. Furthermore,following the description of strong FPG flows in Sreenivasan (1982), these effects arecharacteristic of the laminarescent boundary layer that may precede relaminarization.Regardless of the nomenclature used, when considered along with the non-dimensionalvorticity contours in figure 7, the velocity profiles provide further evidence that: 1) strongFPG cause the breakdown of the standard logarithmic law and a departure above it whilethe flow is still fully turbulent, 2) the breakdown of these standard laws comes with achange in the fundamental character of the turbulence (significant changes in the vorticalstructures in figure 7 can be seen to gradually occur across the x/L = − .
22 line at allvalues of n/ ˜ δ ), and 3) these changes to the velocity occur at much smaller valuesof the non-dimensional pressure gradients relative to the critical ones ( K = 1 . × − where ∆ p = ∆ τ = − . U ∞ . The same streamwise stationsas those in figure 8 are presented. The streamwise fluctuations develop a large innerpeak and an outer knee point, which progressively increases and decreases, respectively,as the turbulence advects through the pressure gradient. The TKE, which is heavilydependent on u (cid:48) s , follows very similar trends. The increase of the streamwise Reynoldsstress near the wall is consistent with the instantaneous contours of vorticity and the wallstreaks observed in figures 3 and 7. The wall-normal fluctuations u (cid:48) n show an increase inintensity in the initial part of the FPG before the deviation above the logarithmic law,however after this event a significant reduction in the peak value is observed. The peakvalue moves further from the wall, further reducing this component of the stress nearthe wall. The spanwise Reynolds stress behaves similarly, with only a slight change at x/L = − .
22, followed by a decrease further downstream. The shape does not appear tobe altered in a major way, rather the magnitude is scaled down. The fact that all threenormal components of the Reynolds stresses decrease in the outer layer, resulting in thereduction of the TKE, is consistent with the slice at n/ ˜ δ = 0 . − u (cid:48) s u (cid:48) n exhibits interesting features. At x/L = − .
22, the outer layer showsan increase in the stress, however, near the wall, the profile shows a sudden break (morevisible in figure 10 where the profiles are not as close together), suggesting a change in thenear-wall turbulence. Downstream of this location, the profiles develop a bi-modal shape,with an inner peak which appears to be relatively fixed in magnitude and distance fromthe wall ( n/ ˜ δ = 0 . x/L = − . et al. et al. x/L = − .
29. Moreover, in the presence of internal layers, the inner and outer layersbecome almost independent of each other, with the latter behaving similarly to a free-
NS of a Boundary Layer with Strong Pressure Gradients Figure 9.
Mean Reynolds stresses and turbulent kinetic energy in the FPG region of theGaussian bump normalized by U ∞ . shear layer (Baskaran et al. x/L = − .
29. Since the curvature changes direction approximately half way along theFPG, these two inner layer effects appear to be strongly related to the pressure gradient2
R. Balin and K. E. Jansen
Figure 10.
Mean turbulent kinetic energy (a), Reynolds shear stress (b), and total shearstress (c) in the FPG region of the Gaussian bump normalized by wall units. rather than the curvature. In the outer layer, the behavior of the Reynolds stresses isslightly more complex. The TKE is reduced from the start of the FPG and develops aknee point rapidly. On the contrary, the shear stress is observed to gradually increasein the outer layer relative to the profile at x/L = − .
29 until around x/L = − .
24 andthen decrease as shown in figure 9. A similar response is observed for the wall-normaland spanwise components as well. Interestingly, x/L = − .
24 is the location where theparameter ˆ κ indicates the concave curvature effects are a maximum upstream of thebump. This points to the fact that the outer layer shear stress is initially responding tothe curvature effects, which at this location dominate over the FPG effects. Eventually, asthe concavity is reduced and the strength of the FPG increases, pressure effects dominateand the Reynolds stresses in the outer layer are reduced. Downstream of x/L = − .
14, thestabilizing effects of the convex curvature are compounded with the FPG ones resultingin a rapid decrease of the shear in the outer layer.This analysis, particularly of the shear stress, strongly suggests the independence of theinner and outer layers. The former is dominated by favorable pressure gradient effects,while the latter responds to a combination of streamline curvature and pressure gradient.It is the inner layer physics, however, that are mostly responsible for the skin frictioncoefficient over the bump. Consequently, a turbulence model that hopes to be predictiveof a flow of this kind must be able to capture these near-wall effects. Note that theGaussian bump flow considered here is not an isolated test case of these effects, andsimilar examples exist in the literature (Wu & Squires 1998; Uzun & Malik 2018; Matai& Durbin 2019).Further interesting observations arise from the TKE and Reynolds shear stress profilesnormalized by wall units in figure 10. This non-dimensionalization is useful to show
NS of a Boundary Layer with Strong Pressure Gradients τ sn represents the total shear stress,which is the sum of the Reynolds and viscous stresses. As the turbulence advects throughthe FPG, the wall-shear normalized TKE is significantly reduced everywhere across theboundary layer, including at the peak. This trend of the near-wall TKE is in contrastto the one shown in figure 9, and indicates that while the streamwise fluctuations arestrengthening near the wall as the flow is being accelerated, they are not increasingfast enough relative to the increase in wall shear stress. Moreover, the drop in TKE inthe outer layer is magnified with this non-dimensionalization. These trends are in fullagreement with the vorticity contours in Fig 7 given the similar non-dimensionalizationused. The Reynolds shear stress exhibits similar trends to the TKE showing significantreduction due to the FPG. In this non-dimensional form, the inner peak is not constantand instead continues to diminish reaching values as low as 35% of the wall shear stress(a density of one was used in the simulation resulting in u τ = τ w ). The non-dimensionaltotal shear stress in figure 10 clearly shows the significant drop in turbulent stress relativeto the viscous stress both in the inner and outer layers. Considering both the Reynoldsand total shear stress profiles in the figure, at x/L = − .
05 the viscous stress dominatesover the turbulent one until about n + = 50, clearly pointing to a thickening of theviscous sublayer measured in wall units. This in turn leads to a very steep wall-normalgradient at a distance between 10 (cid:54) n + (cid:54)
60, which directly impacts the streamwisemomentum equation and thus the streamwise velocity as seen in figure 8. In other words,it is the steep gradient in the total shear stress that is responsible for the deviation abovethe logarithmic law. Once this gradient reaches a certain magnitude, measured in partby ∆ τ , the breakdown of the logarithmic law is to be expected. Note that while therelative size of the Reynolds shear stress is significantly reduced by the FPG, it is notnegligible and thus the relaminarization process does not reach completion and the flowis still turbulent, albeit only intermittently in the near-wall region. The outer layer, bycontrast, remains fully turbulent with the fluctuations decaying in intensity.Finally, the correlation coefficient defined as C τ = − u (cid:48) s u (cid:48) n (cid:113) u (cid:48) s (cid:113) u (cid:48) n (3.7)is plotted in figure 11 for the usual streamwise locations along the FPG. At x/L = − . − . < x/L (cid:54) − .
14, the streamwise and wall-normal fluctuations become slightly more negatively correlated (note the negative signin (3.7)) everywhere within the BL. Given the strength of the pressure gradient in thisregion, this change is attributed to the concave curvature of the mean streamlines whichis a maximum at x/L = − .
24. Furthermore, from the vorticity contours in figure 7 it isevident that the character of the turbulence is significantly altered by x/L = − .
14, yetthe profile of C τ at this location shows only a slight change. Only downstream of x/L = − .
14, where the curvature changes to convex, C τ decreases faster as the boundary layerapproaches the bump peak where ˆ κ is a local maximum. These results are in agreementwith Narasimha & Sreenivasan (1979); Spalart (1986); So & Mellor (1973), and supportother observations that the strong FPG does not change the correlation between thesetwo fluctuations of the turbulence even though the Reynolds stresses are significantlyreduced and the turbulence character is significantly altered during relaminarization.4 R. Balin and K. E. Jansen
Figure 11.
Correlation coefficient C τ in the FPG region of the Gaussian bump. Bump Peak and Strong Adverse Pressure Gradient
As shown by the skin friction coefficient in figure 2 and the contours of instantaneousvorticity in figures 3 and 7, in the vicinity of the bump peak where the FPG is relaxedand pressure gradient changes from favorable to adverse there is a sudden enhancementof the near-wall vorticity which leads to a rise in the wall shear. This feature is of courserelated to the significant weakening of the near-wall turbulence caused by the strongupstream FPG and the onset of relaminarization. It is, in fact, a partial retransitionto fully turbulent flow, where the word “partial” is used since the upstream flow didnot reach the quasi-laminar state associated with the completion of the relaminarizationprocess. In this section, the details of the flow and turbulence as they move through thissegment of the bump are discussed.Figure 12 shows slices of instantaneous vorticity magnitude normalized by the localtime- and spanwise-averaged vorticity at the wall over the bump peak. The slices aretaken at different heights within the boundary layer, in this case measured in wall unitswith n + . Note that streamwise changes in mean u τ are accounted for, thus the slice isat the same height above the wall in local wall units but not in physical units since u τ is not constant. Moreover, black lines across the domain are used to mark the locationof key events of the boundary layer flow. These are the peak strength of the FPG at x/L = − .
11, the bump peak and change in sign of the pressure gradient at x/L = 0 . C f at x/L = 0 .
05, and approximately half way betweenthe peak and incipient separation at x/L = 0 .
10. Finally, a logarithmic scale is used inorder to highlight turbulent structures with values of the normalized vorticity differingby orders of magnitude.Starting with the slice closest to the wall at n + = 5, turbulent spots, which maybe identified as clumps of small-scale turbulent structures of high intensity, are seento form as early as upstream of the bump peak where the FPG is relaxed. These areintermittent with regions of very quiet and weak vorticity fluctuations, and appear togrow in size and intensity as they progress downstream into the strong APG. As notedby other studies on relaminarizing boundary layers (Narasimha & Sreenivasan 1973;Blackwelder & Kovasznay 1972; Sreenivasan 1982), these are clear signs of a retransitionprocess taking place near the wall. As the FPG is relaxed, the stabilizing effect of theacceleration diminishes and instabilities are allowed to grow once again to produce anew fully turbulent internal layer. The presence of fluctuations of different scales inthe incoming flow makes this process often very sudden and energetic, as is evident in NS of a Boundary Layer with Strong Pressure Gradients Figure 12.
Instantaneous vorticity magnitude normalized by the local time- andspanwise-averaged wall vorticity at different locations within the boundary layer in thevicinity of the bump peak. From top to bottom, the heights above the wall of the slices are n + = 5 , , , figure 12. It is important to note that this phenomenon is a feature of the upstreamstrong FPG and relaminarization process, and is not due to the APG. In fact, previousstudies exhibited this phenomenon even with the strong FPG relaxing into a ZPG region.Nevertheless, the destabilizing effects of the APG certainly aid in the onset of instability6 R. Balin and K. E. Jansen and growth of the turbulent spots, accelerating the formation of fully turbulent flow. At x/L = 0 .
05, just four local boundary layer thicknesses downstream of the bump peak,the boundary layer at this height above the wall is effectively fully turbulent. The scalesat this location are small and fairly isotropic in shape, but grow in size and stretch inthe streamwise direction by x/L = 0 .
10. The structure of the canonical boundary layer(see figure 7 at the start of the FPG) is not recovered however.The other slices in figure 12 show the height of the turbulent spots and their growth ratein the wall-normal direction. In their initial stages, the turbulent spots extend well past n + = 30 but are barely visible at n + = 100. Signs of the new internal layer only becomevisible further downstream at this height. This indicates clearly that this phenomenonoriginates near the wall and propagates away from it as the internal layer grows into theboundary layer. The highest location above the wall, n + = 300, shows the flow in theouter wake of the boundary layer (see the velocity profiles in figure 13) as it flows fromthe FPG into the APG. In this region, no change in the size, structure, and strength ofthe turbulent scales is visible from the contours of vorticity. The only change is observeddownstream of x/L = 0 .
10 when the internal layer finally reaches this height above thewall. This behavior is consistent with the notion that the inner and outer layers of thisflow are almost independent of each other, with the latter behaving similarly to a free-shear layer subject to the strong convex curvature effects that are present around thepeak of the bump. By contrast, the near-wall physics are dominated by the pressuregradient and their associated relaminarization and retransition that are not directly feltin the outer layer.The mean streamwise velocity profiles are presented in figure 13 for a number oflocations around the bump peak and the initial part of the APG region. Note that x/L = 0 .
02 corresponds to the local minimum in the C f curve. At the end of theFPG and at the bump peak, significant deviation above the logarithmic law is stillvisible. In fact, the deviation continues to grow throughout the entire FPG with thelargest departure being at the peak and slightly downstream of it. At x/L = 0 .
02 thevelocity still exhibits the characteristics of the upstream FPG. This is consistent with thecontours in figure 12 since the flow is still heavily intermittent at this location. Furtherdownstream at x/L = 0 .
05, the effects of the partial retransition become visible. Thevelocity gradient in the standard logarithmic region is reduced along with the thicknessof the viscous sublayer due to the sudden surge in shear stress discussed later in thissection. At the last station plotted, the velocity approaches the standard log law evencloser and appears to have two distinct regions. The first is found below n + = 40 andresembles the canonical turbulent boundary layer shape with a buffer layer and a semi-linear section (although with a different slope and intercept from the standard values).This is caused by a new internal layer forming at the bump peak which has reachedheights around n + = 100 at this streamwise location. The second region appears to be avery pronounced wake above the underlying internal layer which maintains the shape ofthe upstream FPG profiles. This is once again in agreement with the contours in figure 12.Note that due to the strength of the APG and the imminent incipient separation of theflow at x/L = 0 .
19, the velocity does not reach agreement with the standard logarithmiclaw before separation. This result is consistent with the vorticity contours at the wallin figure 3 which show that the standard wall streaks do not develop within the APGregion.The Reynolds stresses for the same locations are presented in figure 14. Starting withthe turbulent kinetic energy, at the end of the FPG and at the bump peak the profileslook similar to the ones in the upstream FPG, with a larger peak than the x/L = − . NS of a Boundary Layer with Strong Pressure Gradients Figure 13.
Mean streamwise velocity profiles in the APG region normalized by wall units. component due to the acceleration of the flow. There is also a further reduction in theTKE of the outer layer due to the continued FPG and convex curvature effects. Enteringthe initial part of the APG, a significant surge the value of the peak is observed, almostdoubling at x/L = 0 .
05 relative to the bump peak and then decreasing again after thisstation. This inner peak is also increasing in thickness with downstream location. Theseare clear signs of the new internal layer produced by the partial retransition. The outerlayer TKE stays almost constant throughout the APG, providing further evidence of theseparation of these two sections of the boundary layer. As the underlying internal layergrows in thickness, the outer layer is simply displaced further from the wall, resulting ina significant growth of the overall boundary layer as seen in figure 4.The turbulent shear stress follows similar trends to the TKE, with some differencesmainly in the outer layer. The inner peak also rises in magnitude and thickens significantlydue to the new internal layer (more than a factor of three at x/L = 0 . x/L = 0 .
10. Sincethe pressure gradient effects are now destabilizing, these are clearly convex curvatureeffects that are still present until x/L = 0 .
14. Once again, the inner part of the boundarylayer is responding to the pressure gradient, while the outer layer is responding to themean streamline curvature.All three normal Reynolds stresses generally behave as expected, increasing nearthe wall and staying fairly constant or slightly decreasing in the outer layer. They doreact to retransition and the adverse pressure gradient slightly differently, however. Thestreamwise fluctuations respond the quickest, with the largest peak in figure 14 actuallyoccurring at x/L = 0 .
02. They also develop a tri-modal shape at x/L = 0 .
10 that is notseen in other components of the Reynolds stress tensor. The wall-normal and spanwisefluctuations respond later but more drastically (the wall-normal fluctuations increase bya factor of seven near the wall in the interval − . (cid:54) x/L (cid:54) . C τ in the initial part of the strong APG.8 R. Balin and K. E. Jansen
Figure 14.
Mean Reynolds stresses and turbulent kinetic energy in the APG region of theGaussian bump normalized by U ∞ . Relative to figure 11, the profiles at the end of the FPG and the bump peak show a slightreduction in correlation near the wall, but a more substantial reduction in the outerlayer due to the convex streamline curvature. As the boundary layer moves through theAPG, the correlation gradually increases near the wall and spreads away from it. This isanother indication of the second internal layer. In the outer layer, the C τ keeps decreasingthrough the APG, reaching negative values by x/L = 0 . NS of a Boundary Layer with Strong Pressure Gradients Figure 15.
Correlation coefficient C τ in the APG region of the Gaussian bump. that develops a second internal layer with significantly increased turbulent intensity andwall shear. This newly energized flow is more resilient to the strong APG, and thusonly exhibits incipient flow separation on the downstream side of the bump. Pressuregradient effects are focused in the near-wall region and are the main drivers of the skinfriction coefficient, whereas the outer layer is fairly independent and shows the continuedaction of convex streamline curvature. A more three-dimensional view of the internallayer and the turbulent spots is shown in figure 16 with an isosurface of the Q-criterion.The weakened near-wall turbulence and intermittent quiet regions at the bump peak arealso visible with the contours of vorticity magnitude at the wall. It must be stated onceagain that the formation of this internal layer is directly dependent on the upstreamFPG effects, and thus the two regions of the flow are highly connected. A turbulencemodel that does not predict the weakening of the near-wall turbulence due to the FPGand instead models the flow as being fully turbulent everywhere will also not predictthe partial retransition at the peak of the bump and consequently compute the wrongAPG response and separation. This is the case of the SA model used for the preliminarysimulations of the bump in figure 2. Alternatively, a model that does predict the FPGweakening effects but not the correct response at the bump peak will also not accuratelycompute the details of the separation. Adequate modeling of both of these physicalphenomena is required for correct prediction of the effects of strong pressure gradientsalternating in rapid succession as shown here, which is a common feature of aerodynamicflows. 3.5. Reynolds Stress Production Rates
Additional interesting and useful quantities for the analysis of the boundary layer overthe Gaussian bump are the production rates for the turbulent kinetic energy and theReynolds shear stress. These quantities are particularly useful to evaluate and improveRANS turbulence models.Figure 17 presents the production rates of the TKE and Reynolds shear stress overthe upstream side of the bump. These quantities are defined as P k = − u (cid:48) i u (cid:48) j ∂u i ∂x j P sn = − (cid:18) u (cid:48) s u (cid:48) k ∂u n ∂x k + u (cid:48) n u (cid:48) k ∂u s ∂x k (cid:19) , (3.8)respectively, where Einstein notation is used to represent the double contraction of twotensors and the indices i , j , and k take the three dimensions of the curvilinear coordinate0 R. Balin and K. E. Jansen
Figure 16.
Isosurface of the instantaneous Q-criterion, Q = 6 . × , and contours ofvorticity magnitude over the peak and downstream side of the Gaussian bump. Figure 17.
Production rate of the turbulent kinetic energy (a) and Reynolds shear stress (b)in the FPG region of the Gaussian bump non-dimensionalized by wall units. The vertical linesmark the boundary layer thickness ˜ δ for each streamwise location. system ( s, n, z ). Note the non-dimensionalization by wall units, which allows a comparisonof the production rates relative to the increase in wall shear due to the acceleration.Moreover, the distance to the wall is expressed in wall units to emphasize the near-wallregion where these quantities are the largest.At the start of the FPG, the profiles appear similar to the ones in the upstream mildAPG, although production of both quantities is notably increased in the outer later dueto the upstream APG and concave streamline curvature. Immediately downstream of thislocation, a significant reduction in the magnitude of the inner layer peaks is observed.Given the destabilizing concave curvature of the streamlines upstream of x/L = − . n + units, indicating once again the thickening of the viscoussublayer. The TKE production peak settles just above a non-dimensional value of 0.1(approximately a factor of 2.5 smaller relative to the start of the FPG) and undergoes onlya small change in the region − . (cid:54) x/L (cid:54) − .
05 between the peak favorable pressuregradient and the skin friction peak. The Reynolds shear stress production follows similartrends, but exhibits a larger reduction of the inner peak relative to the start of the FPG(approximately a factor of five smaller).Further from the wall in the outer layer ( n + (cid:38) x/L = − .
11. Once again, since convex curvature
NS of a Boundary Layer with Strong Pressure Gradients x/L = − .
14, this reduction is mainly attributed to thestrong FPG, although curvature effects are contributing past this location. The strongacceleration, therefore, effectively turns off production in the outer layer, suggesting onceagain how this region behaves as a free-shear layer with only a small dependence on thewall. Figure 17 therefore highlights the effects of a relaminarization process througha weakening of the near-wall turbulent production relative to the rising wall shear andacceleration. The profiles, however, do not become negligible in magnitude, giving furtherevidence of the fact that the relaminarization process does not complete and the flowover the bump is still partially turbulent.
4. Conclusions
Direct numerical simulation was performed of the turbulent boundary layer over aGaussian shaped bump. The smooth surface causes a series of alternating pressuregradients and mean streamline curvature effects which combine to form a multitudeof complex flow physics with significant deviation from standard turbulence behavior.The domain of interest for this study is the portion of the boundary layer from the inflowto the point of incipient separation, with particular emphasis on the strong favorablepressure gradient.Due to a strong acceleration, the boundary layer exhibits the stages of a relaminar-ization process. The standard logarithmic law breaks down soon after the start of theFPG where the pressure gradient and acceleration parameters are far below their criticalvalues. The non-dimensional shear stress parameter, however, was found to be predictiveof the start of the deviation above the logarithmic law. Instantaneous contours of vorticityindicate that the boundary layer is still fully turbulent where the velocity deviates abovethe logarithmic law, although the turbulent character and intensity are altered by theFPG.Under continued acceleration, the near-wall region gradually ceases to be fully tur-bulent and intermittent spots of quiet flow with low vorticity grow in size. This isconsidered to be the onset of relaminarization, the process by which fully turbulentflow progresses towards a quasi-laminar state. This is in spite of the pressure andacceleration parameters never reaching the critical values of ∆ p = − .
025 and K =3 × − , respectively, suggested in previous studies to mark the onset of this process.Relaminarization, however, does not complete and thus a quasi-laminar state is notachieved. In this region of the FPG, the Reynolds shear stress is reduced developing abi-modal shape and the skin friction coefficient drops. When normalized by the local wallstress, the significant reduction of the turbulent kinetic energy, the Reynolds stresses, andtheir production rates clearly describes the weakening of the turbulence relative to theacceleration of the flow.At the peak of the bump, as the pressure gradient changes to adverse, a partialretransition process is observed near the wall. The weakened flow from the upstreamacceleration experiences a sudden enhancement in turbulent intensity, producing spotsof large vorticity and a surge in turbulent kinetic energy and Reynolds shear stress. Theskin friction coefficient is also increased forming a local maximum. The newly energizedboundary layer is more resilient to the strong deceleration on the downstream side, onlyresulting in incipient separation.Two internal layers are also formed in the region of interest, causing the inner andouter regions of the boundary layer to be largely independent of each other during thestrong favorable and adverse pressure gradients. The flow near the wall is dominated bythe pressure gradients and is responsible for the skin friction coefficient, while further2 R. Balin and K. E. Jansen from the wall the turbulence behaves similarly to a free-shear layer subject to pressuregradients and streamline curvature. Curvature effects are found to be negligible in thenear-wall region.RANS computations with the Spalart-Allmaras model of the same Gaussian bumpflow resulted in a significant overprediction of the skin friction coefficient over mostof the favorable pressure gradient and over the bump peak, as well as a much largerseparation bubble due to the inability of the linear eddy viscosity closure to predict thecorrect pressure gradient effects. Similar results are observed with the SST k - ω modelin the context of wall-modeled LES (Balin et al. Declaration of Interests:
The authors report no conflict of interest.
REFERENCESAbe, H.
Journal of Fluid Mechanics , 563–598.
Abe, H., Mizobuchi, Y., Matsuo, Y. & Spalart, P. R.
Center for Turbulence Research Annual Research Briefs 2012 pp. 311–322.
Balin, R., Jansen, K. E. & Spalart, P. R.
Baskaran, V., Smits, A. J. & Joubert, P. N.
Journal of Fluid Mechanics , 47–83.
Baskaran, V., Smits, A. J. & Joubert, P. N.
Journal of FluidMechanics , 377–402.
Blackwelder, R. F. & Kovasznay, L. S. G.
Journal of Fluid Mechanics , 61–83. Bradshaw, P.
Journal of Fluid Mechanics , 387–390. Cal, R. B. & Castillo, L.
Physics of Fluids , 105106. Cavar, D. & Meyer, K. E.
Journal of Fluids Engineering , 111204.
Coleman, G. N., Rumsey, C. L. & Spalart, P. R.
Journal of FluidMechanics , 28–70.
NS of a Boundary Layer with Strong Pressure Gradients Coles, D. E.
Tech. Rep.
R-403-PR.The Rand Corporation.
Finnicum, D. S. & Hanratty, T. J.
AIChE Journal , 529–540. Greenblatt, D., Paschal, K. B., Yao, C-S., Harris, J., Schaeffler, N. W. & Washburn,A. E.
AIAA Journal , 2820–2830. Jansen, Kenneth E., Whiting, Christian H. & Hulbert, Gregory M. α method for integrating the filtered Navier-Stokes equations with a stabilized finite elementmethod. Computer Methods in Applied Mechanics and Engineering . Jimenez, J., Hoyas, S., Simens, M. P. & Mizuno, Y.
Journal of Fluid Mechanics , 335–360.
Kitsios, V., Atkinson, C., Sillero, J. A., Borrell, G., Gungor, A. G., Jimenez, J. &Soria, J.
International Journal of Heat and Fluid Flow , 129–136. Manhart, M. & Friedrich, R.
International Journal of Heat and Fluid Flow , 572–581. Matai, R. & Durbin, P.
Journal of Fluid Mechanics , 503–525.
Na, Y. & Moin, P.
Journal of Fluid Mechanics , 379–405.
Narasimha, R. & Sreenivasan, K. R.
Journal of Fluid Mechanics , 417–447. Narasimha, R. & Sreenivasan, K. R.
Advances inApplied Mechanics , 222–303. Narayanan, M. A. Badri & Ramjee, V.
Journal of Fluid Mechanics , 225–241. Patel, V. C.
Journal of Fluid Mechanics , 185–208. Patel, V. C. & Head, M. R.
Journal of FluidMechanics , 371–392. Schlatter, P., ¨Orl¨u, R., Li, Q., Brethouwer, G., Fransson, J. H. M., Johansson,A. V., Alfredsson, P. H. & Henningson, D. S. Re θ = 2500 studied through simulation and experiment. Physics of Fluids , 051702. Schwarz, A.C. & Plesniak, M. W.
Journal of Fluids Engineering , 787–794.
Shin, J. H. & Song, S. J.
Journal of Fluids Engineering ,011203.
Shur, M. L., Spalart, P. R., Strelets, M. Kh. & Travin, A. K.
Flow, Turbulence and Combustion , 63–92. Shur, M. L., Strelets, M. K., Travin, A. K. & Spalart, P. R.
AIAA Journal , 784–792. Skote, M. & Henningson, D. S.
Journal of Fluid Mechanics , 107–136.
Slotnik, J. P.
STO-MP-AVT-307 . Smits, A. J., Matheson, N. & Joubert, P. N.
Journal of Ship Research ,147. So, R. M. C. & Mellor, G. L.
Journal of Fluid Mechanics , 43–62. So, R. M. C. & Mellor, G. L.
Aeronautical Quarterly , 25–40. R. Balin and K. E. Jansen
Spalart, P.R.
Journal of Fluid Mechanics , 307–328.
Spalart, P. R. & Allmaras, S. R.
Recherche Aerospatiale , 5–21. Spalart, P. R., Belyaev, K. V., Garbaruk, A. V., Shur, M. L., Strelets, M. Kh. &Travin, A. K.
Flow, Turbulence and Combustion , 865–885. Spalart, P. R. & Coleman, G. N.
European Journal of Mechanics – B/Fluids , 169–189. Spalart, P. R. & Shur, M. L.
Aerospace Science and Technology , 297–302. Spalart, P. R. & Watmuff, J. H.
Journal of Fluid Mechanics , 337–371.
Sreenivasan, K. R.
ActaMechanica , 1–48. Trofimova, A. V., Tejada-Martinez, A. E. & Jansen, K.E.
Computers& Fluids , 924–938. Tsuji, Y. & Morikawa, Y.
The Aeronautical quarterly , 15–28. Uzun, A. & Malik, M. R.
AIAA Journal , 715–730. Uzun, A. & Malik, M. R.
AIAA AVIATION Forum . Virtual Event.
Warnack, D. & Fernholz, H. H.
Journal of Fluid Mechanics , 357–381.
Webster, D. R., DeGraaff, D. B. & Eaton, J. K.
Journal of Fluid Mechanics , 53–69.
Whiting, Christian H. & Jansen, Kenneth E.
Williams, Owen, Samuell, Madeline, Sarwas, E. Sage, Robbins, Matthew &Ferrante, Antonino
AIAA Scitech 2020 Forum . Orlando.
Wright, J., Balin, R., Patterson, J. W., Evans, J. A. & Jansen, K. E. arXiv:2010.3407543v1 [physics.flu-dyn] . Wu, X. & Squires, K. D.
Journal of Fluid Mechanics362