Appl. Math. Mech. -Engl. Ed. DOI 10.1007/s10483-017-2248-6 c
Shanghai University and Springer-Verlag Berlin Heidelberg 2017
Applied Mathematics and Mechanics (English Edition)
Near-wall behaviors of oblique-shock-wave/turbulent-boundary-layer interactions∗ Mingsheng YE1 , Ming DONG1,2,† 1. Department of Mechanics, Tianjin University, Tianjin 300072, China; 2. Tianjin Key Laboratory of Modern Engineering Mechanics, Tianjin 300072, China
Abstract A direct numerical simulation (DNS) on an oblique shock wave with an incident angle of 33.2◦ impinging on a Mach 2.25 supersonic turbulent boundary layer is performed. The numerical results are confirmed to be of high accuracy by comparison with the reference data. Particular efforts have been made on the investigation of the near-wall behaviors in the interaction region, where the pressure gradient is so significant that a certain separation zone emerges. It is found that, the traditional linear and logarithmic laws, which describe the mean-velocity profiles in the viscous and meso sublayers, respectively, cease to be valid in the neighborhood of the interaction region, and two new laws of the wall are proposed by elevating the pressure gradient to the leading order. The new laws are inspired by the analysis on the incompressible separation flows, while the compressibility is additionally taken into account. It is verified by the DNS results that the new laws are adequate to reproduce the mean-velocity profiles both inside and outside the interaction region. Moreover, the normalization adopted in the new laws is able to regularize the Reynolds stress into an almost universal distribution even with a salient adverse pressure gradient (APG). Key words shock wave, turbulent boundary layer, direct numerical simulation (DNS), adverse pressure gradient (APG), separation Chinese Library Classification O354.3, O354.5, O357.4 2010 Mathematics Subject Classification 76F40, 76J20, 76L05
1
Introduction
Shock wave-turbulent boundary layer interaction (SWTBLI) is one of the most ubiquitous yet perplexing phenomena occurring in a wide variety of internal and external aerodynamics applications. This phenomenon could lead to significantly negative effects on the performance of high-speed vehicles and turbomachines, e.g., increasing the surface drag and thermodynamical loading and inducing the flow separation. Simultaneously, the unsteady motion of both the shock system and the separation bubble accompanies the SWTBLI, which operates in an extremely low frequency in comparison with the coherent structure of the turbulence and hence causes an increment in the complexity[1] . Through the experiments on a variety of models associated with the boundary layers with laminar, transitional, and turbulent separations, Chapman et al.[2] proposed an early remarkable ∗ Received Jan. 13, 2017 / Revised Mar. 4, 2017 Project supported by the National Natural Science Foundation of China (Nos. 11472189 and 11332007) † Corresponding author, E-mail:
[email protected]
2
Mingsheng YE and Ming DONG
concept, i.e., free interaction. It was found that the pressure rise to separation is independent of the inducing separation mode for either laminar or turbulent separation in supersonic flows. The plateau pressure rise is also independent for laminar separation, but it significantly depends on the model geometry for turbulent separation. The first asymptotic description of a supersonic laminar boundary layer through a free interaction was provided independently by Stewartson and Williams[3] and Neiland[4] . According to their arguments, in the close neighborhood O(Re3/8 L) of the shock/boundary-layer interaction region, three layers, which are expressed in terms of the magnitude O(Re3/8 L) for the upper deck, the magnitude O(Re1/2 L) for the main deck, and the magnitude O(Re5/8 L) for the lower deck, emerge in the wall-normal direction when Re → ∞, where Re is the Reynolds number based on the distance from the local interaction location to the leading edge L. A series of subsequent work, as reviewed by Adamson and Messiter[5] , has depicted the whole scenario of the interactions of the shock waves with laminar boundary layers. It has been revealed that, under the laminar conditions, the heat transfer rates, as well as other dynamical values, even with the high temperature effects, can be predicted by the computational fluid dynamics (CFD) accurately with adequate meshes if the freestream conditions are properly characterized (see Ref. [1]). However, the investigations on the interactions with the turbulent boundary layers are insufficient. This may be due to that the numerical predictions from the Reynolds averaging Navier-Stokes (RANS) equations with turbulence modelling are unsatisfactory and sometimes disastrous, especially on predicting the heat flux and surface drag[1,5–6] . Extensive experimental and numerical investigations[7–17] have been performed to provide comprehension on the relative physical features and serve as the databases for CFD validation. One factor that affects the accuracy of the CFD prediction attributes to the uncertainty of the temperature boundary condition at the wall. As pointed out by Edwards[18], in some cases, a pre-test ambient temperature is designated as the wall temperature of the CFD calculation, but the run time of the wind tunnel is long enough to render this assumption questionable. Also, the SWTBLIs often result in an order of magnitude enhancement on the heat transfer rates, which leads to locally high-temperature phenomena. Therefore, it is rational to take the real-gas effect and the possible chemical reaction into account. However, as illustrated by Knight et al.[19] , the applied chemical models would be crucial to the numerical results, which reinforces the uncertainty of the numerical prediction. Even when the comparisons are made between the RANS solutions and those from the direct numerical simulations (DNSs) or large eddy simulations (LESs), salient discrepancies are always found both qualitatively and quantitatively[17] . The rational improvement on the accuracy of the turbulence models on predicting SWTBLIs requires further understanding on the features of the flow motion, especially the near-wall behaviors. Loginov et al.[20] examined the development of the mean-velocity profiles of a supersonic flow past a compression ramp from the LES solution. A turbulent-laminar reverse transition or relaminarization in the reverse-flow region is observed. It is well-known that for a turbulent boundary layer with zero pressure gradient (ZPG), the mean-velocity profiles in the near-wall region satisfy the linear law in the viscous sublayer and the logarithmic law in the mesolayer. Such behaviors are found to be independent of the flow conditions such as the Mach number, the Reynolds number, and the temperature conditions at the wall (see Wilcox[21] and Duan et al.[22] ). However, with the presence of a sufficiently strong adverse pressure gradient (APG), the traditional laws of the wall cease to be valid. As proposed by Skote and Henningson[23] from their study on the near-wall behavior of the incompressible turbulent boundary layer with separation, the mean-velocity is described by a pair of new laws corresponding to two sublayers. The new relations require the balance among the convection, the shear stress (including both the viscous and Reynolds stresses), and the pressure gradient, and are proven to be adequate via the DNS data. The normalization adopted in the new relations are also able to reproduce the Reynolds stresses in the incompressible turbulent boundary layer with separation (see Gungor
Near-wall behaviors of oblique-shock-wave/turbulent-boundary-layer interactions
3
et al.[24] ). The objective of the present paper is to shed some additional light on the near-wall behavior of the interactions in the SWTBLI. The laws of the wall as proposed by Skote and Henningson[23] for the incompressible separation are followed, but the compressible effect is considered additionally due to the supersonic nature.
2
Numerical methods
2.1 Physical model In order to demonstrate the physical model to be studied in this paper, we plot the contours of the instantaneous temperature on a spanwise slice (see Fig. 1). We focus on the interaction between an impinging oblique shock wave and a turbulent boundary layer on a flat plate. The Cartesian coordinates (x, y, z), which, respectively, denote the streamwise, wall-normal, and spanwise directions, are employed, in which the reference length δin is selected as the momentum thickness of the boundary layer at the inlet of the computational domain. It is clear from Fig. 1 that, the flow remains laminar in the very early stage. Since the artificial perturbation with the form of blowing and suction at the wall is introduced near the inlet, the laminar-turbulent transition emerges at about x = 220. The turbulence develops further downstream until it reaches the outlet of the computational domain. Moveover, an oblique shock wave is introduced from the inviscid region, and impinges on the turbulent boundary layer at about x = 420, which generates even higher level fluctuations in the turbulent boundary layers. Simultaneously, a reflect shock wave is observed downstream of the interactive region. This interaction is often accompanied with flow separation and low-frequency unsteadiness, which challenges the accuracy of turbulence-modeling prediction and consequently is of our interest. /P-SFGMFDUJOHDPOEJUJPO
*ODJEFOUTIPDLXBWF
/POQIZTJDBM[POF
-BNJOBSJOGMPXZ
3FHJPOPGCMPXJOH BOETVDUJPO
Fig. 1
5SBOTJUJPO 'VMM-EFWFMPQFE *OUFSBDUJPO 3FDPWFSZ [POF UVSCVMFOU[POF [POF [POF Y
Contour of the instantaneous temperature on a spanwise slice
2.2 Numerical methods The governing equations are the three-dimensional Navier-Stokes equations for a perfect gas, which can be found in Qin and Dong[25] . The DNS code is kindly provided by Prof. Xinliang LI of Institute of Mechanics, Chinese Academy of Science. The code has been extensively used for a variety of simulations involving turbulent boundary layers and shock waves[26–28] . All the physical quantities are normalized by the freestream quantities, i.e., ρ=
ρ∗ , ρ∞
ui =
u∗i , u∞
p=
p∗ , ρ∞ u2∞
T =
T∗ , T∞
µ=
µ∗ , µ∞
4
Mingsheng YE and Ming DONG
where ρ, ui (i = 1, 2, 3), T , p, and µ denote the density, the velocities, the temperature, the pressure, and the dynamic viscosity, respectively, the superscript ∗ represents the dimensional quantities, and the subscript ∞ represents the dimensional quantities at the freestream. The Mach number is defined by M a∞ = u∞ /c∞ , and the Reynolds number is defined by Re = ρ∞ u∞ δin /µ∞ , where c∞ is the sound speed. The seventh-order weighted essentially nonoscillatory (WENO) scheme[9] is used for the nonlinear convection terms, and the eighth-order central scheme is used for the viscous terms. The time integration is performed by use of an explicit forth-order Runge-Kutta scheme. 2.3 Computational setup The computational parameters for the simulation of the SWTBLI conducted in this paper are listed in Table 1. A blowing and suction slot between x = 20 and x = 40 (see Fig. 1) is introduced as a perturbation generator, which subsequently leads to the emergence of the instability modes of the Tollmien-Schlichting type. Due to their amplification, the TollmienSchlichting waves eventually trigger the laminar-turbulent transition. This approach is also used by Pirozzoli and Grasso[9]. The turbulence (see Fig. 1) appears downstream of the transitional region when x is approximately 220, and there is a sufficiently long distance for the development before interacting with the incident shock wave. Table 1
Numerical configurations of the simulation
M a∞
Re
β/(◦ )
T∞ /K
Lx × Ly × Lz (δin )
Nx × Ny × Nz
2.25
3 700
33.2
169.44
520 × 60 × 27
2 530 × 180 × 251
Simultaneously, we introduce an oblique shock with the angle β = 33.2◦ from the upper boundary of the computational domain. The physical quantities ahead of the shock are related to those behind the shock by the Rankine-Hugoniot jump condition. The incident shock wave impinges on the turbulent boundary layer, and generates a reflected shock wave. Through resolution study, the computational domain and the grid point number are finally chosen as those listed in Table 1. The grid points in the streamwise direction are uniformly located except in the shock-turbulent-interacting region (x ∈ [390, 450]), where the mesh is clustered due to the high resolution requirement. In the wall-normal direction, the grid points are clustered in the near-wall region according to a hyperbolic sine transformation, and the uniform grid points are located in the spanwise direction. In the full-developed turbulent region, the grid spacing, expressed in terms of the viscous wall unit, is ∆x+ × ∆y + × ∆z + = 7.73 × 0.4 × 6.0, where (x+ , y + , z + ) =
uτ (x, y, z). νw
In the above equation, νw is the kinematic viscosity at the wall, and uτ is the friction velocity defined by s ∂e u uτ = νw . (1) ∂y y=0 Throughout the paper, the subscript “w” represents the mean quantities at the wall, and the bar and the tilde over the quantities represent the Reynolds and Favre averaging, respectively. The minimum grid spacing in the streamwise direction occurs in the shock-turbulent-interacting region, where ∆x+ = 3.03.
Near-wall behaviors of oblique-shock-wave/turbulent-boundary-layer interactions
5
Because the simulation starts from the laminar boundary layer, the Blasius similarity solution is introduced as the inflow profile. The no-slip and non-penetration boundary conditions are imposed at the wall, except at the blowing and suction slot. The wall is assumed to be isothermal, and the wall temperature is approximately equal to the adiabatic wall temperature of the turbulent boundary layer, i.e., Tw = 332 K. The non-reflecting boundary condition is used at the upper boundary. In order to diminish the perturbation-reflection from the outflow boundary, an additional sponge region is used from x = 520 to x = 575.6 with 120 grids stretching towards the downstream direction. The periodic boundary condition is used in the spanwise direction.
3
Baseline validation
In order to compare our simulation results with the reference data, we first compare the main properties of the turbulent profile at x = 380 with those from Pirozzoli and Bernardini[14] and Morgan et al.[17] (see Table 2). Among the symbols in the table, φ represents the flow deflection angle, i.e., 2 cot β(M a2 sin2 β − 1) ∞ , (2) φ = arctan M a2∞ (γ + cos 2β) + 2 Reδ , Reθ , and Reδ∗ are the local Reynolds numbers whose reference lengths are the the nominal thickness (δ0 ), the momentum thickness (θ), and the displacement thickness (δ ∗ ) of the boundary layer, respectively, H = δ ∗ /θ denotes the shape factor, and the skin-friction coefficient is defined by τw Cf = , ρ∞ u2∞ /2 where τw is the shear stress at the wall defined by τw = µw (
∂e u )y=0 . ∂y
The three oncoming Mach numbers in Table 2 are fairly close, but the Reynolds numbers of our simulation are apparently larger than those of the reference papers. It is satisfactory that quite good agreement on the sensitive quantities such as the shape factor and friction velocity is observed, and the slightly smaller value of Cf in our simulation is anticipated due to the larger local Reynolds number. Table 2
Comparison of the simulation parameters and mesh scale
Parameter
Present DNS
Ref. [14]
Ref. [17]
M a∞
2.25
2.28
2.28
φ
8.1
8
8 26 067
Reδ
43 031
29 000
Reθ
3 563
2 330
2 295
Reδ∗
12 543
8 320
7 817
H
3.53
3.6
3.53
Cf
2.34 ×
103
2.51 ×
103
2.61 × 103
uτ /u∞
0.046 7
0.049 3
0.045 1
+ + ∆x+ min × ∆yw × ∆z
3.03 × 0.4 × 6.0
5.6 × 0.93 × 6.6
−
The mean streamwise velocity at the reference location x = 380 is plotted in Fig. 2. Two empirical curves are shown in Fig. 2, which, respectively, represent the linear law + u+ c = y
(3)
6
Mingsheng YE and Ming DONG
%/4 V DZ V DMOZ
V D
Fig. 2
Z
Mean velocity at the reference station xref = 380
in the viscous sublayer and the logarithmic law 1 ln y + + B κ
u+ c =
(4)
with κ = 0.4 and B = 5.5 in the mesolayer. The velocity and the wall-normal coordinate in Eq. (4) are normalized by wall units, i.e., u+ c =
uc , uτ
y+ =
yuτ , νw
(5)
where uc =
Z
0
u e
ρ 2 de u ρw 1
(6)
is the van Driest transformation for compressible flows. The numerical mean-velocity agrees well with the two empirical curves in the regions y + 6 5 and 306y + 6 100, respectively. The implication is that the fully-developed turbulent boundary layer has already been formed at this location. It has to be noticed that an obvious discontinuity is observed at y + ≈ 1 000, which is induced by the incident shock wave. The distribution of the skin-friction coefficient Cf and the wall pressure pw in the vicinity of the interaction region are plotted in Fig. 3, where the streamwise coordinate is rescaled by (x − x0 )/Lsep as suggested by Morgan et al.[17] , in which x0 = 436.69 denotes the nominal shock impingement point, and the separation length Lsep is defined as the distance between the separation point (xsep = 414.16) and the reattachment point (xrea = 428.7). From Fig. 3, we can see that our results agree well with those of Refs. [14] and [17]. The only difference is that the Cf curve in our simulation is lower than the reference curves outside the interacting zone, which is attributed to the relatively larger Reynolds number that we have selected. Under the normalization suggested by Pirozzoli and Bernardini[14] , Fig. 4 exhibits the contours of the averaging and fluctuation velocities near the interaction region, where the subscript “rms” denotes the root-mean-square. In the sub-figures, the streamwise and wall-normal coordinates are normalized as follows: ξ=
x − x0 , L
η=
y , L
(7)
respectively, where the interaction length L is defined as the distance between x0 and the origin point of the reflected shock (xon = 402). From the figures, we can see that our results agree with those of Fig. 3 in Ref. [14] and Fig. 12 in Ref. [17].
Near-wall behaviors of oblique-shock-wave/turbulent-boundary-layer interactions
QX
$G
%/4 3FG<> 3FG<>
YY -TFQ
%/4 3FG<> 3FG<>
B $G
Fig. 3
Z-
Z-
Skin friction coefficient Cf and wall pressure pw
B V
C W
Z-
Z-
YY -TFQ C QX
D VSNT
Fig. 4
7
E VĄ WĄ
Contours of the averaging and the fluctuation velocities near the interaction region
Figure 5 shows the distribution of root-mean-squire of the wall pressure, the streamwise velocity, and the Reynolds stress. Although the overall agreement between our results and those of Ref. [14] is satisfactory, yet not surprisingly, there remain some nuances in the detail. (i) A slightly higher wall pressure fluctuation arises at both the upstream and the downstream limits, which is understandable since the Reynolds number chosen in our simulation is larger. (ii) The curves from our simulation oscillates more apparently, which is believed to be related to the insufficiency of the statistical time interval. Because the frequency of the shock oscillation is extremely low, it is too costly to include the shock motion in our sampling series. However, the time interval is much larger than the large-scale coherent structures in the turbulent flow. As pointed out by Piponniau et al.[12] , this could at most lead to an error of 1% on the first- and second-order statistics. This is why our results agree with those of Pirozzoli and Bernardini[14] quantitatively.
8
Mingsheng YE and Ming DONG
VSNT
QSNT
%/4 3FG<>
%/4 3FG<>
B QSNT
C VSNTBUZ-
VĄ WĄ
%/4 3FG<>
D VĄ WĄ BUZ-
Fig. 5
4
Comparison of the fluctuations near the interaction region
Statistical properties
4.1 Mean flow Figure 6 plots the distribution of the skin-friction coefficient (Cf ). At the sufficient upstream locations, the values of Cf stay almost unchanged, drop dramatically in the early stage of the interaction region, and become negative after about x = 414.16, which indicates the separation point. After attaining the minimum value at about x = 416, the values of Cf vary gradually, until they reach the close neighborhood at the reattachment point at about x = 429. Then, the Cf curve increases with x monotonically. Simultaneously, as shown in Fig. 7, an APG is produced in the separation zone. The pressure starts to rise at the same dropping location of Cf . After peaking at about x = 417.9, the pressure gradient gradually decreases when x approaches further downstream. However, its value has not reached zero up to the outflow boundary of the computational domain, implying that the flow field has not yet recovered to the upstream turbulent state. The distribution of the heat flux qw is plotted in Fig. 8, where ∂ Te qw = − k¯ , ∂y y=0
and k is the thermal conductivity. In the upstream undisturbed zone, because the wall temperature Tw is close to the turbulent recovery temperature, the wall heat flux is quite small. However, in the vicinity of the interaction zone, a large amount of heat is transferred to the
Near-wall behaviors of oblique-shock-wave/turbulent-boundary-layer interactions
QX
$G
QX Q Y
Y
Fig. 6
Distribution of the skin-friction coefficient
Q Y
9
Y
Fig. 7
Distributions of the wall pressure and its gradient
h
RX
Y
Fig. 8
Distribution of the heat transfer rate
wall and the qw peaks inside the separation zone. In order to monitor the detailed flow feature, 14 sampling points are selected, whose streamwise coordinates and their particular significance are summarized in Table 3. Table 3
Streamwise locations of the sampling points
Section
x
Comment
s1
390.0
Upstream of the interaction region (ZPG)
s2
406.0
Near xon
s3
409.6
s4
413.7
s5
416.0
s6
418.0
s7
420.0
s8
422.6
Near xsep Peak location of the pressure gradient
s9
424.0
s10
426.0
s11
429.3
Near xrea
s12
436.7
Nominal shock impingement point x0
s13
445.0
Downstream of the interaction region
s14
490.0
10
Mingsheng YE and Ming DONG
Figures 9(a) and 9(b) plot the profiles of the mean streamwise velocity and temperature, respectively. The incident and reflected shocks in the inviscid region are apparent from the profile saltation. In the interacting region, where remarkable APG exists, some certain deficiencies appear on both the velocity profiles and the temperature profiles. With the decrease in the APG, the profiles gradually recover to the fully-developed turbulent profile. In the simulation, the separation zone is so tiny that it is inconspicuous in Fig. 9(a). We therefore zoom in the the mean-velocity profiles in the close neighborhood of the wall, i.e., 0 6 y 6 0.1 (see Fig. 9(c)). Two separation bubbles with different sizes are clearly observed from the reverse-velocity zones, i.e., the larger one appears from s4 to s8, and the smaller one appears in the vicinity of s10. The two bubbles corresponds to the two valleys of the Cf curve (see Fig. 6).
T
T
T
T
T
T
T
T
T
T
T
T
T
T
T
T
T
T
T
T
Z
_ V
T
_ B 6 T
T
T
T
T
T
T
Z
_ 5 _ C 5
T
T
T
T
T
T
T
Z
Fig. 9
_ V _ D ;PPN-JOPG6JOŏZŏ
Wall-normal profiles at different stremwise locations, where the scales of s2, s3, · · · , s14 are the same as those of s1 in (a) and (b), and the scales of s5, s6, · · · , s10 are the same as that of s4 in (c)
4.2 Reverse velocity in the separation zone A clearer observation on the separation bubbles are shown in Fig. 10(a), where the contours of the mean-velocity and the streamlines in the near-wall region are plotted.
Near-wall behaviors of oblique-shock-wave/turbulent-boundary-layer interactions
11
As pointed out by Zheltovodov[29–30] , the reverse flow ahead of the forward-facing step can be considered as a near-wall jet, and the laminar and the turbulent profiles satisfy different distributions. Based on the suggestion, Loginov et al.[20] demonstrated that there was a relaminarization effect in the separation bubble induced by a supersonic flow past a compression ramp. For the present simulation, we perform the same analysis on the mean-velocity profiles in the separation zone (see Fig. 10(b)). In the figure, the mean-velocity is normalized by the maximum value of the reverse flow umax , and the wall-normal coordinate is normalized by the location ζ1/2 at which the mean-velocity reaches half of its maximum value. The red solid and dashed lines in Fig. 10(b), respectively, represent the laminar and turbulent near-wall-jet profiles suggested by Vulis and Kashkarov[31]. The velocity profiles from our simulation scatter around the laminar-jet profile, and keep far away from the turbulent profile, which implies that the flow in the separation zone is very close to the laminar state. No reverse transition from turbulent to laminar is observed for the oblique shock wave impingement on the turbulent boundary layers, which is different from the cases of the supersonic flow past compression ramps in Ref. [20].
-BNJOBS 5VSCVMFOU
Z
Z
Y _ B $POUPVSTPG VBOETUSFBNMJOFT
Fig. 10
_ V V
NBY
C .FBOWFMPDJUZQSPGJMFT
Mean-velocity in the vicinity of the separation zone
Furthermore, in the vicinity of the separation zone, the APG induced by the oblique shock wave balances the total shear stress to the leading-order accuracy, i.e., ∂p ∂τ ⋍ , ∂x ∂y
(8)
where the total shear stress τ is the sum of the viscous stress and the Reynolds shear stress, i.e., ∂e u ′′ v ′′ . τ ⋍µ − ρu] (9) ∂y Figure 11 plots the wall-normal distributions of both the viscous stress and the Reynolds stress in the near-wall region of the separated zone. The conspicuously common feature in the separation zone is that the viscous stress is dominant whereas the Reynolds stress is almost 0. In comparison with the behavior of the turbulent boundary layer with ZPG, the shear rate of the viscous stress is no more zero, but becomes nearly a constant for each x. As determined by Eq. (8), the shear rate balances the APG. Therefore, we compare the two quantities at different x (see Table 4). Both the shear rates of the viscous stress and the pressure gradient peak are at the same location (s6), and quantitative agreement is observed for all the locations.
12
Mingsheng YE and Ming DONG
h
7JTDPVT 5VSCVMFOU
h 7JTDPVT 5VSCVMFOU
Z
B T
Z
C T
h
7JTDPVT 5VSCVMFOU
h 7JTDPVT 5VSCVMFOU
Z
D T
Z
E T
Fig. 11
Wall-normal distributions of the viscous stress (solid line) and the Reynolds stress (dashed line) in the reverse-flow region, where the circles represent the upper boundary of the separation zone
Table 4
Wall shear rates of the viscous stress and the pressure gradient in the reverse-flow region (×10−3 )
Section “ ” ∂ u ˜ µ ¯ ∂∂y ∂y ∂p ∂x
s5
s6
s7
s8
s9
s10
8.00
9.12
8.21
5.78
5.06
4.82
8.68
10.27
9.45
6.84
5.94
5.51
4.3 Near-wall behavior of the mean-velocity in the interaction region 4.3.1 Introducing the pressure-velocity scale For turbulent boundary layers with zero or small pressure gradient, the wall friction velocity uτ , which is defined in Eq. (1), is selected as the reference velocity to demonstrate the near-wall behavior of the mean-velocity. Such a behavior is also described in Fig. 2 of the current paper. It is found that the linear law (3) and the logarithmic law (4) emerge in the viscous and meso sublayers, respectively. However, the presence of the APG in the SWTBLIs breaks the balance of the viscous stress and the Reynolds stress in the turbulent boundary layer. Therefore, the traditional laws of the wall cease to be valid. Furthermore, ( ∂∂yu˜ )y=0 in the reverse flow becomes negative, which leads to the invalidity of the definition (1).
Near-wall behaviors of oblique-shock-wave/turbulent-boundary-layer interactions
13
For the incompressible turbulent boundary layer with the APG, Skote and Henningson[23] suggested to choose the pressure velocity scale up =
ν d¯ 1/3 w p ρw dx
(10)
as the reference velocity for normalization, and redefined the dimensionless velocity up = and the dimensionless wall-normal coordinate yp = y
u e up
(11)
up . νw
(12)
Figure 12(a) shows the distributions of uτ and up in the neighbourhood of the interaction region, where the definition of uτ is changed to s ∂e u uτ = νw . ∂y y=0
(13)
γ = uτ /up
(14)
In the upstream of the separation zone, the value of uτ is about 0.047, which is much larger than the value of up (about 0.003). In the vicinity of the separation point, uτ falls rapidly, and is overwhelmed by up . In the separation zone, uτ and up become comparable until they are in the close neighborhood of the reattachment point. After that uτ recovers to the dominant factor. In order to better compare the magnitudes of uτ and up , we introduce
to denote their ratio. Figure 12(b) exhibits the distribution of γ near the interacting region. The ratio is much larger than 1 outside the separation zone. However, the ratio is of the order 1 inside the separation zone, and is close to zero near the separation and reattachment points.
V VQ
VVQ
V VQ
Y B V BOEVQ
Fig. 12
Y C V VQ
Distributions of uτ , up , and γ = uτ /up
14
Mingsheng YE and Ming DONG
4.3.2 Mean velocity in the near-wall region Since the leading-order terms in the x-momentum equation are the pressure gradient, the viscous shear stress, and the Reynolds shear stress, Eqs. (8) and (9) can be recast to d ] d de u dp − (ρu′′ v ′′ ) = . µ dy dy dy dx
(15)
′′ v ′′ , µ, and ρ are functions of y since the Here, we assume that, p is a function of x, and u e, u] quantities vary slowly along the other direction. Multiplying Eq. (15) by νw /(ρw u3p ), applying Eqs. (10)–(12), and introducing
D(y) =
ρ , ρw
W (y) =
ν , νw
′′ v ′′ p u] ′′ v ′′ = u] u2p
give
p
′′ v ′′ ) dup d(Du] d DW − = 1. p p dy dy dy p
(16)
Inserting Eqs. (11) and (12) into Eq. (13) yields
dup = γ2. dy p y=0
(17)
Integrate Eq. (16) with respect to y p , apply the condition (17), and notice D(0) = W (0) = 1. Then, we obtain p dup ′′ v ′′ = ±γ 2 + y p , DW p − Du] (18) dy where signs + and − represent the profiles without and with separation, respectively. It is noticed that the left-hand side of Eq. (18) is essentially the total shear stress. 4.3.3 Viscous sublayer In the region very close to the wall, i.e., y p = O(1), the Reynolds stress is negligible, i.e., ′′ v ′′ ≈ 0. −ρu]
Therefore, the dominant shear stress is the viscous stress, and this region is referred to as the viscous sublayer. Integrate Eq. (18) with respect to y p , and notice (D, W, up ) → (1, 1, 0) when y → 0. Then, we can obtain an explicit relation between up and y p , i.e., up = ±γ 2 y p +
(y p )2 . 2
(19)
Sufficiently far from the separation zone, where γ ≫ 1, Eq. (19) recovers to the linear law of the turbulent boundary layer with ZPG, i.e., up = γ 2 y p .
(20)
Moreover, in the close neighborhood of the separation and the reattachment points, where γ ≪ 1, Eq. (19) becomes a parabolic law, i.e., up = (y p )2 /2, in which the pressure gradient plays the leading role.
(21)
Near-wall behaviors of oblique-shock-wave/turbulent-boundary-layer interactions
15
4.3.4 Mesolayer The Reynolds stress comes to the leading order in the mesolayer, where y p ≫ 1. We herein introduce a new characteristic velocity scale u∗ , which is defined by u2∗ = u2p (±γ 2 + y p ),
(22)
where the signs + and − have the same definitions as those in Eq. (18). Multiplying Eq. (18) by u2p /u2∗ gives p Du2p dup ′′ v ′′ ] W − u = 1. (23) u2∗ dy p Based on the Boussinesq hypothesis, the Reynolds stress can be expressed as follows: ′′ v ′′ = ν −u] T
de u , dy
(24)
where νT is the eddy viscosity. According to the concept of the mixing length, νT ∼ κ2 y 2
de u , dy
(25)
where κ = 0.4 is the Karman constant (see Ref. [21]). Applying Eqs. (24) and (25) to Eq. (23) and noticing y p ≫ 1, we have √ dup u∗ yp D p = . (26) dy κup It can be seen that the viscous stress plays a minor role in the mesolayer. For compressible flows, the van Driest transformation is introduced, i.e., upc =
Z
up
√ Ddup .
(27)
0
Then, the solution of Eqs. (26) is obtained to be upc =
p p 1 (γ ln y p + 2 γ 2 + y p − 2γ ln( γ 2 + y p + γ)) + C κ
(28)
for the profiles without separation and upc
1 p p 2 y − γ 2 − 2γ arctan = κ
r
yp − 1 +C γ2
(29)
for the profiles with separation, where C is constant. If γ ≫ 1, where the pressure gradient is very small, Eq. (28) recovers to the logarithmic law upc =
γ ln y p γ ln γ + + Bγ, κ κ
(30)
where B = 5.5 is the constant from the classical turbulent theory. In the vicinity of the separation and the reattachment points, where γ ≪ 1, both Eqs. (28) and (29) satisfy a squireroot law 2√ p y + C, (31) upc = κ implying a natural communication when the profiles cross the zero-shear points.
16
Mingsheng YE and Ming DONG
4.3.5 Verification by the DNS data Figure 13 shows the comparisons of the near-wall mean-velocity between the DNS data and the above analysis, where the dashed lines are from Eq. (19) for the viscous sublayer and Eq. (28) or Eq. (29) for the mesolayer, while the dot-dashed lines are from Eq. (20) for the viscous sublayer and Eq. (30) for the mesolayer. The latter are essentially the same as the linear and logarithmic laws (3) and (4) of the classical turbulent theory for the ZPG. At all the streamwise locations, the present laws of the wall can predicate the mean-velocity profiles to high accuracy in both the viscous sublayer and the mesolayer. However, the traditional near-wall laws show poor agreement with the DNS data, especially inside the separation zone, which indicates the significant role played by the pressure gradient. Admittedly, the constant C in Eqs. (28) and (29) are derived from curve fitting instead of from analytical approach. The reason is that the matching condition from the outer region of the boundary layer is not determined due to the effect of the impinging shock wave. Nevertheless, the importance of the new law of the wall is clear because the dominant mechanism is obtained. It is also noted that the buffer layer connecting the viscous sublayer and mesolayer thickens inside the separation zone. We however have not derived the explicit relation for this layer due to its complexity. For the turbulent boundary layer with the ZPG, the second-order statistics in the near-wall region, such as the Reynolds stress, also exhibit universal relations (local Reynolds number ′′ v ′′ /(ρ u2 ) against y + . However, for the SWTindependent) when plotting −hρu′′ v ′′ i+ = −ρu] w τ BLI problem, such curves at different streamwise locations differ from each other remarkably (see Fig. 14(a)). Gungor et al.[24] suggested that the characteristic stress scale should be ρw u2p instead of ρw u2τ for the turbulent boundary layer with the APG. Noticing that the APG only concentrates in the close neighborhood of the separation zone for the current problem, we therefore choose ρw u2∗ as the reference stress scale. Thus, we introduce − hρu′′ v ′′ i∗ = − y∗ =
u∗ y, νw
′′ v ′′ ρu] , ρw u2∗
(32) (33)
and plot −hρu′′ v ′′ i∗ against y ∗ for different streamwise locations (see Fig. 14(b)). All the curves exhibit the same trend both qualitatively and quantitatively.
5
Conclusions
The interaction of an impinging oblique shock wave with a supersonic turbulent boundary layer at the Mach number 2.25 is investigated by the DNS. The main purpose of this paper is to provide a comprehensive analysis on the near-wall behaviors of the interaction, such that further works on improving the accuracy of the turbulent modelling could be inspired. Our numerical results are compared with the reference data from Pirozzoli and Bernardini[14] and Morgan et al.[17] with similar configurations, and are shown to be of high accuracy and reliable. Additional efforts have been made on providing the detailed turbulent statistics such as the distribution of the surface friction coefficient, the pressure gradient, the heat flux, and the mean flow profiles. These results are adequate to serve as a standard data-set on further turbulence-modelling validations. The reverse flow inside the separation zone is found to be laminar through the analysis of the local mean profiles, which implies that the Reynolds stress here is negligible. Simultaneously, the viscous stress linearly increases with the wall-normal coordinate y, instead of being a constant
Near-wall behaviors of oblique-shock-wave/turbulent-boundary-layer interactions
%/4 1SFTFOU 0SJHJOBM
VQD
ZQ B T $
%/4 1SFTFOU 0SJHJOBM
VQD
VQD
%/4 1SFTFOU 0SJHJOBM
VQD
17
ZQ D T $
%/4 1SFTFOU 0SJHJOBM
%/4 1SFTFOU 0SJHJOBM
VQD
VQD
ZQ E T $
%/4 1SFTFOU 0SJHJOBM
ZQ C T $
ZQ F T $
%/4 1SFTFOU 0SJHJOBM
ZQ G T $
%/4 1SFTFOU 0SJHJOBM
VQD
VQD
ZQ H T $
Fig. 13
ZQ I T $
Validation of the law of the wall in the interaction region
18
Mingsheng YE and Ming DONG
Z B LJ ρVąWąLj
Fig. 14
T T T T T T
LJ ρVąWąLj
LJ ρVąWąLj
T T T T T T
Z
C LJ ρVąWąLj
Distributions of the Reynolds stress scaled by ρw u2τ and ρw u2∗ in the separation zone
for the ZPG cases. The shear rate of the viscous stress agrees perfectly with the local pressure gradient, indicating the balance between the viscosity and the pressure gradient in this region. Therefore, a new u e-y relation in the viscous sublayer is formed due to the additionally dominate role of the pressure gradient, which, being different from the traditional linear law for the ZPG, exhibits a parabolic dependence (19) (with a negative sign in front of γ) of the mean-velocity on y. Right above this sublayer, the Reynolds stress comes to the leading order as well. Therefore, the leading balance is among the viscosity, the pressure gradient, and the Reynolds stress. This region is referred to as the buffer layer, and no explicit u e-y relation has been obtained so far. As y further increases to the mesolayer, the viscosity becomes negligible, and the balance between the pressure gradient and the Reynolds stress determines the relation (29). In the regions upstream and downstream of the separation zone, the parabolic law (19) (with a positive sign in front of γ) of the viscous sublayer still holds, and it gradually recovers to the linear law (20) when the pressure gradient reduces successively. Moreover, the relation in the mesolayer is expressed as Eq. (28) in the attached region, and also recovers to the logarithmic law (30) in the small pressure-gradient limit. The two relations (29) and (28) in the mesolayer communicate with each other at the separation and reattachment points for the common squireroot dependent of u e on y. The new laws of the wall proposed in this paper are proved to be effective on reproducing the mean-velocity profiles from the DNS results both inside and outside the separation zone. The distribution of the Reynolds stress in the near-wall region is remarkably regularized under the normalization approach utilized in the new law of the wall. Acknowledgements
The authors are grateful to Professor Xinliang LI of Institute of Mechanics, Chinese Academy of Sciences for kindly providing the DNS code. Professors Heng ZHOU and Jisheng LUO of Tianjin University are also acknowledged for their valuable discussion and suggestions. This work is carried out at the National Supercomputer Center in Tianjin, and the calculations are performed on TianHe-1 (A).
References [1] Gaitonde, D. V. Progress in shock wave/boundary layer interactions. Progress in Aerospace Sciences, 72, 80–99 (2015) [2] Chapman, D. R., Larson, H. K., and Kuehn, D. M. Investigation of separated flows in supersonic and subsonic streams with emphasis on the effect of transition. NASA Technical Reports, NASA, United States (1957)
Near-wall behaviors of oblique-shock-wave/turbulent-boundary-layer interactions
19
[3] Stewartson, K. and Williams, P. G. Self-induced separation. Proceedings of the Royal Society of London, A312, 181–206 (1969) [4] Neiland, V. Y. Theory of laminar boundary layer separation in supersonic flow. Fluid Dynamics, 4, 33–35 (1972) [5] Adamson, T. C. and Messiter, A. F. Analysis of two-dimensional interactions between shock waves and boundary layers. Annual Review of Fluid Mechanics, 12, 103–138 (1980) [6] Roy, J. C. and Blottne, G. F. Review and assessment of turbulence models for hypersonic flows. Progress in Aerospace Sciences, 42, 469–530 (2006) [7] Settles, G. S. and Dodson, L. J. Hypersonic shock/boundary layer interaction database. NASA Technical Reports, NASA, United States (1991) [8] Knight, D., Yan, H., Panaras, A. G., and Zheltovodov, A. Advances in CFD prediction of shock wave turbulent boundary layer interactions. Progress in Aerospace Sciences, 39, 121–184 (2003) [9] Pirozzoli, S. and Grasso, F. Direct numerical simulation of impinging shock wave/turbulent boundary layer interaction at M = 2.25. Physics of Fluids, 18, 1–17 (2006) [10] Ganapathisubrama, B., Clemens, N. T., and Dolling, D. S. Low-frequency dynamics of shockinduced separation in a compression ramp interaction. Journal of Fluid Mechanics, 636, 397–436 (2009) [11] Humble, R. A., Scarano, F., and van Oudheusden, B. W. Unsteady flow organization of a shock wave/turbulent boundary layer interaction. Journal of Fluid Mechanics, 635, 47–74 (2009) [12] Piponniau, S., Dussauge, J. P., Debieve, J. F., and Dupont, P. A simple model for low-frequency unsteadiness in shock-induced separation. Journal of Fluid Mechanics, 629, 87–108 (2009) [13] Touber, E. and Sandham, N. D. Large-eddy simulation of low-frequency unsteadiness in a turbulent shock-induced separation bubble. Theoretical and Computational Fluid Dynamics, 23, 79–107 (2009) [14] Pirozzoli, S. and Bernardini, M. Direct numerical simulation database for impinging shock wave/turbulent boundary-layer interaction. AIAA Journal, 49, 1307–1312 (2011) [15] Priebe, S. and Martin, M. P. Low-frequency unsteadiness in shock wave-turbulent boundary layer interaction. Journal of Fluid Mechanics, 699, 1–49 (2012) [16] Priebe, S., Tu, J. H., Rowley, C. W., and Martin, M. P. Low-frequency dynamics in a shock-induced separated flow. Journal of Fluid Mechanics, 807, 441–477 (2016) [17] Morgan, B., Duraisamy, K., Nguyen, N. S., Kawai, S., and Lele, K. Flow physics and RANS modelling of oblique shock/turbulent boundary layer interaction. Journal of Fluid Mechanics, 729, 231–284 (2013) [18] Edwards, J. R. Numerical simulations of shock/boundary layer interactions using time-dependent modeling techniques: a survey of recent results. Progress in Aerospace Sciences, 44, 447–465 (2008) [19] Knight, D., Longo, J., Drikakis, D., Gaitonde, D., Lani, A., Nompelis I., Reimann, B., and Walpot, L. Assessment of CFD capability for prediction of hypersonic shock interactions. Progress in Aerospace Sciences, 48-49, 8–26 (2012) [20] Loginov, M. S., Adams, N. A., and Zheltovodov, A. A. Large-eddy simulation of shockwave/turbulent-boundary-layer interaction. Journal of Fluid Mechanics, 565, 135–170 (2006) [21] Wilcox, D. C. Turbulence Modeling for CFD, DCW Industries, United States (1998) [22] Duan, L., Beekman, I., and Martin, M. P. Direct numerical simulation of hypersonic turbulent boundary layers: part 3, effect of Mach number. Journal of Fluid Mechanics, 672, 245–267 (2011) [23] Skote, M. and Henningson, D. S. Direct numerical simulation of a separated turbulent boundary layer. Journal of Fluid Mechanics, 471, 107–136 (2002) [24] Gungor, A. G., Maciel, Y., Simens, M. P., and Soria, J. Scaling and statistics of large-defect adverse pressure gradient turbulent boundary layers. International Journal of Heat and Fluid Flow, 59, 109–124 (2016)
20
Mingsheng YE and Ming DONG
[25] Qin, H. and Dong, M. Boundary-layer disturbances subjected to free-stream turbulence and simulation on bypass transition. Applied Mathematics and Mechanics (English Edition), 37(8), 967–986 (2016) DOI 10.1007/s10483-016-2111-8 [26] Li, X. L., Fu, D. X., Ma, Y. W., and Liang, X. DNS of Shock/Boundary Layer Interaction Flow in a Supersonic Compression Ramp, Springer-Verlag, Heidelberg, 729–737 (2011) [27] Li, X. L., Leng, Y., and He, Z. W. Optimized sixth-order monotonicity-preserving scheme by nonlinear spectral analysis. International Journal for Numerical Methods in Fluids, 73(6), 560– 577 (2013) [28] Li, X. L. Direct numerical simulation techniques for hypersonic turbulent flows. Hangkong Xuebao/Acta Aeronautica et Astronautica Sinica, 36, 147–158 (2015) [29] Zheltovodov, A. A. Analysis of properties of two-dimensional separated flows at supersonic speeds conditions (in Russian). Investigations of Near-wall Flows of Viscous Gas, USSR Academy of Sciences, Russia, 59–94 (1979) [30] Zheltovodov, A. A. Shock waves/turbulent boundary-layer interactions — fundamental studies and applications. Journal of Chemical Physics, 67, 677–698 (2013) [31] Vulis, L. A. and Kashkarov, W. P. Theory of the Viscous Liquid Jets, Nauka, Moscow (1965)