Next Article in Journal
Contributions to the Teaching and Learning of Fluid Mechanics
Next Article in Special Issue
Experiments in Shock-Vortex Interactions
Previous Article in Journal
On the Use of Probability-Based Methods for Estimating the Aerodynamic Boundary-Layer Thickness
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Simulation of Pressure Fluctuation near an Expansion Corner in a Supersonic Flow of M = 3.01

Department of Engineering Mechanics, Tsinghua University, Beijing 100084, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 1 July 2021 / Revised: 23 July 2021 / Accepted: 23 July 2021 / Published: 28 July 2021
(This article belongs to the Special Issue High Speed Flows)

Abstract

:
The influence of the expansion corner on pressure fluctuation is an important subject in supersonic flow around high-speed vehicles. Past studies have clarified how the expansion corner alters the root-mean-square of the fluctuating pressure coefficient ( C p r m s ) and the power spectral density (PSD) without considering how these fluctuating properties are related to compressible waves. In this paper, we use characteristics to determine the direction of wave propagation and identified three zones—U-zone, M-zone and D-zone—within which both C p r m s and PSD are likely to display different behaviors across the boundary layer. The U-zone is upstream of the characteristic line of the second family and passing through the corner. The D-zone is downstream of the characteristic line of the first family and passing through the corner. The middle zone lies between the U-zone and D-zone. The results of C p r m s and PSD at different layers within the boundary layer are obtained using numerical computation through a Detached Eddy Simulation (DES). It is found that in the U-zone and D-zone, both C p r m s and PSD are the same in different layers within the boundary layer. In the M-zone, however, both C p r m s and PSD may vary in different layers and this variation occurs in the high-frequency band upstream of the corner and mid-frequency band downstream of the corner. A feedback mechanism is tentatively used to explain the difference of spatial distribution of fluctuation properties inside the M-zone.

1. Introduction

The fluctuating pressure acting on the surface of the vehicle induces structural vibration, which may cause damage to the aircraft structure and a strong noise environment that affects the normal operation of airborne instruments, including the reliability and safety of weapons and equipment [1]. The intensity of the pressure fluctuations is commonly characterized by the root-mean-square fluctuating pressure coefficient ( C p r m s ), Power Spectral Density (PSD) and correlation coefficient [2]. Pressure fluctuation beneath a supersonic turbulent boundary layer (c.f. S. Beresh and J. Henfling and R. Spillers and B. Pruett [2]) may be amplified by shock–boundary layer interaction (c.f. M. Holden [3], H. Babinsky and J. Harvey [4]) and altered by geometry, such as forward step (V.Bibko and B. Efimtsov and V. Kuznetsov [5]) and expansion corners (c.f. [6]). Here in this paper, we consider pressure fluctuations due to expansion corners, which are typical geometric configurations that change the local flow properties, leading to the formation of a complex fluctuating environment [6]. This topic has received a number of studies.
Fluctuating pressure for subsonic flow around expansion corners has been studied experimentally by J. Robertson [7], who measured the surface fluctuating pressure in the vicinity of the expansion corner and found that the peak of the fluctuating pressure exists at the reattachment point, and this point moves backwards when the Mach number increases. Moreover, the fluctuating pressure of the expansion corner has been studied analytically by K. Plotkin and J. Roberson [6] based on the experimental data on the wall, and they found the relationship between the C p r m s and the local Mach number of the separated flow caused by several expansion corners (cone-cylinder, back step, etc.). They found that the C p r m s is closely related to the local Mach number and that the fluctuating pressure is stronger at low Mach numbers and decreases at high Mach numbers. Later, X. Ligong and L. Zhenhuan [8] also used the same formula to predict the C p r m s at the expansion corner, he took the characteristic length as the displacement thickness of the local boundary layer and obtained results that are more consistent with the previous experiment. Furthermore, the surface pressure fluctuations of subsonic turbulent flow downstream of small expansion corners were found to be normally distributed through the expansion process but were severely attenuated [9]. In addition, the fluctuating pressure on an expansion corner has been investigated numerically using the Detached Eddy Simulation (DES) method by S. Deck and P. Thorigny [10]. They focused on the surface fluctuating pressure statistical properties and found that the spectral analysis of the pressure fluctuations has shown different frequency contributions depending on the location considered in the recirculation bubble. Close to the expansion corner, the spectrum displays a peak near a normalized frequency, which has been shown to represent the footprint of the vertical motion of the bubble.
Fluctuating pressure for transonic flow around expansion corners has been studied experimentally by D. Depres and P. Reijasse and J. Dussauge [11]. The experiments have been carried out to investigate the unsteady fluctuating surface-pressure field on afterbodies of revolution at transonic speeds with a freestream Mach number 0.85, and it was found that the spectra of fluctuating pressure in the whole region exhibit a well-defined periodicity that corresponds to the formation of large-scale structures in the wake, and when the shear layer reattaches near the end of the protruding wall, the rms pressure fluctuations reach their maximum value. Moreover, Z. Rui, R. Jili and R. Fang [12] used the large eddy simulation (LES) method to simulate the structure of a typical expansion corner numerically, and they recorded data on how pressure changes over time on the surface. Additionally, they found that the fluctuating pressure in the turbulent boundary layer, separation zone and other flow structures changes significantly with space; the fluctuating pressure in the separation zone is generally higher than the turbulent boundary layer; and the separation reaction on the shock wave will cause a stronger fluctuating pressure environment and shock waves to self-oscillate [13]. Then, based on previous experimental data, they provided an empirical formula for the P S D of the fluctuating pressure of the expansion-reflection separation flow. Combining the calculation results of LES, they introduced spatial-related information based on the empirical formula published by K. Plotkin and J. Roberson [6], which improved the prediction accuracy of engineering algorithms for fluctuating pressure [14].
Fluctuating pressure for supersonic flow around expansion corners has been studied experimentally by C. Kungming and F. Lu [9]. In this work, surface pressure fluctuations of Mach 8 turbulent flow past a 2.5- and a 4.25-deg expansion corner maintained a Gaussian distribution but were severely attenuated by the expansion process. The pressure fluctuations did not recover to those of an equilibrium turbulent flow even though the mean pressure reached downstream inviscid values in four to six boundary-layer thicknesses. The fluctuations were convected with a velocity comparable to that on a flat plate, and they maintained their identities longer for the stronger expansion. The damping of pressure fluctuations at hypersonic Mach numbers, even by small corner angles, may be exploited in fatigue design. Moreover, fluctuating pressure near expansion corners has been studied experimentally by J. Dawson and M. Samimy [15], wherein multipoint wall pressure measurements were used to investigate the response of a Mach 3.01, fully developed, compressible, turbulent boundary layer to centered and gradual expansions of both 7- and 14-deg deflection. Furthermore, they found that although rms fluctuation levels decrease across the expansions, the rms normalized by the local static pressure remains nominally constant. Just downstream of the expansions, normalized power spectra are more concentrated at low frequencies than upstream, suggesting small-scale turbulence is quenched. This spectra alteration is more prominent for centered expansions and larger deflections. The spectra evolve very quickly after the centered expansions and very slowly after the gradual expansions. In the area near the expansion corner, there are separation points and reattachment points, a shock wave will appear at each location, and the fluctuating pressure environment will be more complicated [16].
K. Plotkin and J. Roberson [6] analyzed the typical subsonic, transonic and supersonic flow field and its fluctuating pressure of the reentry body with two shapes: cone-cylinder and cone-cylinder-skirt. They pointed out that at subsonic speeds, fluctuating pressure is caused by the appearance of the separation zone, and at transonic speeds, the airflow behind the cone-cylinder shoulder may reach supersonic speed, and shock waves will induce the separation of the boundary layer and cause strong fluctuating pressure. At supersonic speeds, there is a certain degree of separation after the expansion corner, and the fluctuating pressure environment will be more complicated.
Past studies have focused on fluctuating pressure on the surface. Here in this paper, we consider the spatial distribution of fluctuation properties in the vicinity of the expansion corner when the flow is supersonic. The problem we will use is the experimental model of J. Dawson and M. Samimy [15]. In Section 2, we will give the details of this model and present the method we use to obtain the pressure fluctuation properties. We will use characteristic lines to split the region near the corner into three zones—the U-zone, M-zone and D-zone—and study how Cprms and PSD change across the boundary layer in different zones. The description of characteristics and definition of the U-zone, M-zone and D-zone will be given in Section 3. The computed spatial distribution of Cprms and PSD in three zones around the expansion corner is presented in Section 4. The behavior of spatial distribution of fluctuating properties in the three zones will be summarized in Section 5, where we tentatively provide a possible feedback mechanism to explain the observed behavior. The conclusion will be summarized in Section 6.

2. The Expansion Corner Problem and Numerical Method for Simulation

In this paper, the experimental model of J. Dawson and M. Samimy [15], with an expansion of 14 deg, as shown in Figure 1, is used. The upstream flow is supersonic, with a freestream Mach number M a = 3.01 . The Reynolds number based on the momentum thickness ( δ 2 = 0.37 mm) is R e δ 2 = 24,700. The incoming pressure is p = 22 KPa. J. Dawson and M. Samimy [15] measured the normalized rms pressure fluctuations and power spectra on the body surface.
Since we are interested in spatial distribution of the pressure fluctuation properties, we need numerical simulation to obtain these properties. Note that pressure fluctuation for flow around an expansion corner has been considered numerically using Reynolds Averaged Navier–Stokes (e.g., R. Soni and N. Arya and A. De [17]), Detached Eddy Simulation (DES) (e.g., S. Deck and P. Thorigny [10]), large eddy simulation (LES) (e.g., M. Grilli and S. Hickel and N. Adams [18]), and direct numerical simulation (DNS) (e.g., M. Kopera and R. Kerr and H. Blackburn and D. Barkley [19]).
To balance the accuracy and time requirement, we use the well-established DES method, which combines the features of the classical RANS formulations [20] with elements of LES method.
The domain and boundary conditions are shown in Figure 1. The left side is supersonic inlet, and the right boundary is a supersonic exit. The wall has a no-slip condition. A small segment of the free-sliding wall is set at the front edge of the non-slip wall. The calculation results are compared with J. Dawson and M. Samimy [15] in the paper to verify the accuracy and reliability of the CFD calculation method we chose.
The flow equation is solved by means of the second-order-accurate scheme of implicit Roe using finite difference simulation. The time step is fixed to Δ t = 5 e 7 s, which corresponds to 10 4 time steps for one passage of the model at freestream velocity. The freestream turbulent in tensity is T u = 3 % , and there are three grids containing 7.5, 10.36 and 23.6 million nodes for computation, respectively. Several grid spaces have been considered in successive refinement to study the influence of the mesh density on the fluctuating pressure of the expansion corner. Based on the results of analysis, a grid of 10.36 million nodes with refinement inside the boundary and near the expansion corner is finally used. The grid has a wall normal resolution y + ~ 0.5 and x + ~ 0.5 in the densified area (similarly as in the work of [21]), and the increasing rate is 1.1. A study on the time independence is performed and the results of fluctuating pressure at four different times are given in Figure 2. In the early stage of the calculation, the PSD keeps changing with time. Until 0.37 s, the PSD does not vary in time. Thus, we will extract the data of 0.45 s.
As displayed in Figure 3, the numerical results of static pressure are comparable with the experimental results of J. Dawson and M. Samimy [15].
As displayed in Figure 4, the numerical results of the normalized fluctuating pressure on the wall are comparable with the experimental results of J. Dawson and M. Samimy [15]. In Figure 4, x is the position of flow direction, p ¯ 2 is the root-mean-square fluctuating pressure, p is the static pressure, x / δ 0 = 0 is the location of the expansion corner, and the thickness of the boundary layer δ 0 at the starting point of the expansion is 9.2 mm.
Here, we use a Detached Eddy Simulation (DES) to obtain the instantaneous flow field. The instantaneous pressure p ( t ) at a given point in the flow field (not just on the body surface) is recorded, and C p r m s is computed as C p r m s = p ¯ 2 0.5 ρ U 2 , p ¯ 2 = 1 T 0 T ( p ( t ) p a v g ) d t , where T is the time interval for sampling, p a v g is the average pressure during the entire sampling process, and p ¯ 2 is the root-mean-square (rms) of the pressure fluctuation, which represents the total intensity of pressure. The power spectral density (PSD) (named ϕ ( ω ) ) is obtained by the Fourier transform of the pressure change in the time domain ϕ ( ω ) = 1 T 0 T p ( t ) e i ω t d t , where ω is the frequency (Hz).

3. Characteristics and Definition of U-Zone, M-Zone and D-Zone

Characteristics are well defined for inviscid and isentropic flow. Outside the boundary layer, the flow may be regarded as inviscid and isentropic. For inviscid and isentropic flow, the characteristic line is defined by
d y d x = 1 λ
where λ is the eigenvalue. For the characteristic line of the first family, λ = λ 1 , and for characteristic line of the second family, λ = λ 2 , where
λ 1 = u v a 2 M a 2 1 v 2 a 2 λ 2 = u v + a 2 M a 2 1 v 2 a 2
Here, u and v are the local flow velocity components; M a is the local Mach number; a is the local sound speed.
One way to obtain the characteristic lines is to use the local flow parameters obtained by inviscid CFD and then use Equations (1) and (2) to integrate the characteristic lines. Here, we use the local flow parameters from CFD results and apply Equations (1) and (2) to find the characteristic lines as if the flow is inviscid. Note that inside the boundary layer and very close to the wall, M a < 1 so that we have no real values of λ 1 and λ 2 . In such regions, we simply extend the characteristic lines to the wall using straight lines.
The characteristic lines and Mach contours in the vicinity of the expansion corner are displayed in Figure 5. Not surprisingly, the Mach number is constant along the characteristic line of the second family and connecting the expansion corner.
Now we define three zones using the two characteristic lines of the first and second families connecting the expansion corner O. These two characteristic lines, “ C U ”and “ C D ”, are shown in Figure 6. The intersection of C U and the edge of the boundary layer is denoted by “U”; the intersection of C D and the edge of the boundary layer are denoted by “D”.
The region near the expansion corner is divided into three zones:
(1)
U-zone is the region upstream of the characteristic line C U .
(2)
M-zone is the region between the characteristic lines C U and C D .
(3)
D-zone is the region downstream of the characteristic line C D .
In the next section, we will display how C p r m s and PSD vary across the boundary layer inside these three zones.
The distribution of the measuring points is shown in Figure 7a,b. Measuring points are evenly distributed in U-zone, D-zone, M-zone, and along the vertical direction of the wall surface. The measuring points are arranged at different layers within the boundary layer in these three zones:
(a)
wall surface, with y + = 0 and labelled point “1”.
(b)
buffer layer, with y + 30 and labelled point “2”.
(c)
log-law layer, with y + 300 and labelled point “3”.
(d)
outer layer, with y + > 1000 and labelled point “4”.
For each layer, position “o” is at the expansion corner, positions “o−”, “a”, “b”, “c”, “d”, “e”, “f”, “g”, “h”, “i”, “j” and “k” are upstream of the expansion corner, and positions “o+”, “l”, “m”, “n”, “o”, “p”, “q”, “r”, “s”, “t”, “u”, “v”, “w”, “x”, “y” and “z” are downstream of the expansion corner. The subscripts 1 , 2 , 3 and 4 for these points correspond to the four layers.
For positions “j” and “k”, all the four layer points ( 1 , 2 , 3 , 4 ) are upstream of C U and are inside the U-zone. For positions “y” and “z”, the four layer points ( 1 , 2 , 3 , 4 ) are downstream of C D and are inside the D-zone.
Usually, the pressure fluctuation properties on the wall surfaces were considered. Here, we consider these properties at the four different layers. They have the same values at these four layers at a fixed horizontal location in the case of a pure boundary layer. We wonder if they may have different values in the vicinity of the expansion corner, and this will become clear in the next section.

4. The Cp rms and PSD in the U-Zone, M-Zone and D-Zone

Now we display the pressure fluctuation properties C p r m s and P S D at the four different layers (wall, buffer layer, log-law layer and outer layer) in the U-zone, M-zone and D-zone.

4.1. The Distribution of C p r m s in Various Zones

The distribution of C p r m s at the four layers are displayed in Figure 8. The abscissa in Figure 8 corresponds to the different positions in Figure 7a,b. The ordinate corresponds to the value of C p r m s . The four curves correspond to the four layers.
Looking at the positions of the points, as shown in Figure 7a,b, we observe that for points below the characteristic line C U , i.e., inside the U-zone, the C p r m s is almost constant in the vertical direction, i.e., they are the same at different layers of the boundary layers. For points below the characteristic line C D , i.e., inside the D-zone, the C p r m s is almost constant in the vertical direction, i.e., they are the same at different layers of the boundary layers.
In the M-zone, the C p r m s changes in different layers of the boundary layer. Inside the M-zone, the C p r m s is the highest in the buffer zone (point “2”) and is the lowermost in the outer layer (point “4”).
In approaching the expansion corner from the upstream along the flow direction, starting from position i for the measuring points that lie between characteristic line “ C U ” and characteristic line “ C D ”, the values of C p r m s began to differ across the vertical direction of the boundary layer, and the difference gradually increased along the flow direction. After the expansion corner, the difference of C p r m s gradually decreases at different layers in the direction normal to the surface within boundary layer. Until position x, this difference almost shrinks to zero and continues to develop downstream along the flow direction. The values of C p r m s across the vertical direction of the boundary layer do not change.

4.2. The P S D in Various Zones

The distribution of the P S D at position “d~k” (sufficiently far upstream from the expansion corner) is displayed in Figure 9. For positions k and j, where all the four layer points are inside the U-zone, the curves of the P S D at y + = 0 (wall surface), y + = 30 (buffer layer), y + = 300 (log-law layer) and y + > 1000 (outer layer) almost overlap; that is, the frequency properties of fluctuating pressure do not change much in the vertical direction within the U-zone. For positions i and h, the values of P S D at y + > 1000 (outer layer) is lower than other layers in the high-frequency band, while for positions g~d, the values of P S D in the high-frequency band gradually decrease when the distance from the wall surface increases.
The distribution of P S D at positions “c~m” (in the vicinity of the expansion corner) is displayed in Figure 10. For positions c~o−, the P S D in the high-frequency band gradually decreases with the increasing distance from the wall surface. In approaching the expansion corner along the flow direction, the difference of PSD in the high-frequency band gradually increases, and at position o−, this difference reaches the maximum value. In return, along the streamwise direction starting from position o, the P S D becomes to show obvious differences in the mid-frequency band, and with increasing distance from the wall surface, the P S D in the mid-frequency band gradually decreases.
The distribution of P S D at positions “n~v” (sufficiently far downstream from the expansion corner) is displayed in Figure 11 and, at position “w~z”, is displayed in Figure 12. For positions n~z, the P S D in the mid-frequency band gradually decreases with the increasing distance from the wall surface, and in the direction away from the expansion corner, the difference of PSD in the mid-frequency band gradually decreases. For positions y and z, where all the four layer points are inside the D-zone, the curves of P S D at y + = 0 (wall surface), y + = 30 (buffer layer), y + = 300 (log-law layer) and y + > 1000 (outer layer) almost overlap along the vertical direction of the wall surface within the boundary layer; that is, the frequency properties of fluctuating pressure do not change much in “D-zone”.

5. Summary of Fluctuating Properties and Feedback Mechanism

Based on the observation presented in Section 4, we use Figure 13 to summarize the distribution of the fluctuating properties in the vicinity of the expansion corner. There are three zones, U-zone, M-zone and D-zone, as defined in Section 3, using the characteristic lines C U and C D .
In the U-zone and D-zone, the two key parameters of fluctuating pressure— C p r m s and P S D —are almost constant along the vertical direction of the wall surface within the boundary layer.
The M-zone can be divided into two subzones: M 1 zone is upstream of the expansion corner, and M 2 zone is downstream. The C p r m s changes significantly along the vertical direction of the wall surface within boundary layer. In the M 1 zone, when approaching the expansion corner, the difference of the C p r m s gradually increases along the vertical direction of the wall surface, and the PSD varies in the high-frequency band. In the M 2 zone, the difference of C p r m s in each layer gradually decreases to zero when the distance from the expansion corner becomes large, and the variation of the PSD occurs in the mid-frequency band downstream of the corner. Moreover, in the upstream of the expansion corner, the difference of P S D in the high-frequency band gradually increases when approaching expansion corner and the difference of P S D in mid-frequency band gradually decreases to disappear along the flow direction.
It is desirable to understand why the fluctuation properties have different spatial distributions in the M-zone compared to the U-zone and D-zone. Here, we tentatively provide an explanation. For this purpose, we design a triangle (with one curved side) in Figure 13. Point O of this triangle is near the expansion corner. Point P is upstream so that the side OP lies inside the subsonic region within the boundary layer. The point Q is above and near the characteristic line C U , such that PQ is perpendicular to the wall and the curved side QO is close and parallel to the characteristic line C U .
The expansion corner is supposed to yield a perturbation of the pressure fluctuation, and this perturbation can propagate upstream along Q P to the point P since O P is in the subsonic region within the boundary layer (information can propagate upstream in subsonic flow). This perturbation can then propagate to point Q along P Q . Since Q O is along the propagation direction of the characteristic line C U , the above-mentioned perturbation (originally from the expansion corner) can propagate back to O. This defines a feedback loop, which may explain why the pressure fluctuation properties change inside the M 1 -zone compared to that inside the U-zone. The feedback mechanism can be similarly defined to explain the change of fluctuation properties in the M 2 -zone compared to that inside the D-zone.
It is also interesting to see the correlation of instantaneous pressure at two different points once these points cross the characteristic line C U or C D .
The correlation coefficient r A B of the instantaneous pressure at two different points A and B is
r A B = c o v ( p A ( t ) , p B ( t ) ) V a r [ p A ( t ) ] · V a r [ p B ( t ) ]
where p A ( t ) is the instantaneous pressure at point A; p B ( t ) is the instantaneous pressure at point B; c o v ( p A ( t ) , p B ( t ) ) is the covariance of p A ( t ) and p B ( t ) ; V a r [ p A ( t ) ] is the variance of p A ( t ) ; and V a r [ p B ( t ) ] is the variance of p B ( t ) .
The points A and B will be both inside zones U, M or D or in different zones among U , M and D. The correlation coefficients computed by Equation (3) are displayed in Figure 14, where three curves are displayed. For the wall surface and buffer layer curve, point A is in the wall layer and point B is in the buffer layer, for the buffer layer and log-law layer curve, point A is in the buffer layer, and point B is in the log-law layer. For the log-law layer and outer layer curve, point A is in the log-law layer, and point B is in the outer layer. The abscissa corresponds to different positions, and the ordinate corresponds to the correlation coefficient.
We observe that, inside the U-zone and D-zone, the correlation coefficient between any two layers is near 1.
If point A is inside the U-zone and point B is inside the M-zone, the correlation coefficient is less than 1. Consider, for instance, position c. At this position, the wall surface is in the “U-zone”, and the buffer layer is in the “M-zone”, so the correlation coefficient between wall surface and buffer layer is 0.96.
If both points are inside the M-zone, the correlation coefficient is much lower than 1. Consider, for instance, position c, the log-law layer and outer layer are in the “M-zone”, so the correlation coefficient between buffer layer and log-law layer is 0.71, and the correlation coefficient between log-law layer and outer layer is 0.62.
If point A is inside the D-zone and point B is inside the M-zone, the correlation coefficient is less than 1. Consider, for instance, position n. At this position, the wall surface is in the “D-zone”, and the buffer layer is in the “M-zone”, so the correlation coefficient between wall surface and buffer layer is 0.94. Moreover, both the log-law layer and outer layer are in the “M-zone”, so the correlation coefficient between buffer layer and log-law layer is 0.45, and the correlation coefficient between log-law layer and outer layer is 0.24.

6. Conclusions

In this work, the fluctuating pressure in the vicinity of the expansion corner in supersonic flow is studied using numerical results from DES simulation. We use characteristics (as if the flow is inviscid and isentropic) to determine the direction of wave propagation and identify three different zones — U-zone, M-zone and D-zone. These zones are bound by the characteristic lines of the first and second families and start from the expansion corner (in an approximative way).
In the U-zone, which is upstream of the upstream characteristic line connecting the expansion corner, both C p r m s and PSD are almost constant along the vertical direction within the boundary layer, and the correlation coefficient between two vertical points within this zone is close to 1. In the D-zone, which is downstream of the downstream characteristic line connecting the expansion corner, both C p r m s and PSD are almost constant along the vertical direction within the boundary layer, and the correlation coefficient between two vertical points within this zone is close to 1.
In the “M-zone”, however, the values of C p r m s differ across the vertical direction of the boundary layer, and the difference gradually increases along the flow direction in regions upstream of the expansion corner. Downstream of the expansion corner, the difference of C p r m s gradually decreases and almost shrinks to zero finally. Moreover, the variation of PSD occurs in the high-frequency band upstream of the corner and mid-frequency band downstream of the corner. The correlation of C p r m s between two points with at least one point inside the M-zone is less than 1.
A feedback mechanism is tentatively used to explain the difference of spatial distribution of fluctuation properties inside the M-zone: the change of fluctuation properties by the expansion corner can propagate upstream from the subsonic region inside the boundary layer and then propagate vertically to the M-zone, where this perturbation can be brought back to the expansion corner.

Author Contributions

L.Z. and Z.-N.W. have equally contributed to the development of the present paper, including the conceptualization, methodology, data acquisition and analysis, interpretation of the results, and manuscript preparation. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported partly by the National Key Project (Grant No. GJXM92579) and the National Science and Technology Major Project 2017-II-003-0015.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We are grateful to all reviewers who provided valuable suggestions that improve the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shiwei, M.; Huabing, J. Discussion on prediction methods of fluctuating pressure environments of flow fields surrounding the aircraft. Equip. Environ. Eng. 2021, 18, 14–22. [Google Scholar]
  2. Beresh, S.J.; Henfling, J.F.; Spillers, R.W.; Pruett, B.O.M. Fluctuating wall pressures measured beneath a supersonic turbulent boundary layer. Phys. Fluids 2011, 23, 75–110. [Google Scholar] [CrossRef] [Green Version]
  3. Holden, M.S. Shock Wave-Turbulent Boundary Layer Interaction in Hypersonic Flow. AIAA J. 1972, 1–30. [Google Scholar]
  4. Babinsky, H.; Harvey, J.K. Shock Wave-Boundary-Layer Interactions; Cambridge University Press: Cambridge, UK, 2011; pp. 1–461. [Google Scholar]
  5. Bibko, V.N.; Efimtsov, B.M.; Kuznetsov, V.B. Pressure fluctuations on a plate ahead of a step. Uch. Zap. TsAGI 1986, 17, 170–173. [Google Scholar]
  6. Plotkin, K.; Roberson, J. Prediction of Space Shuttle Fluctuating Pressure Environments, Including Rocket Plume Effects; NASA: San Jose, CA, USA, 1973; pp. 1–163.
  7. Robertson, J. Wind Tunnel Investigation of the Effects of Reynolds Number and Model Size on the Steady and Fluctuating Pressures Experienced by Cone Cylinder Missile Configurations Attransonic Speeds; AEDC: Ascension Parish, LA, USA, 1967; pp. 152–161.
  8. Ligong, X.; Zhenhuan, L. Analysis and prediction of fluctuating pressure environment of reentry aircraft. Acta Aerodyn. 1991, 4, 457–464. [Google Scholar]
  9. Kungming, C.; Lu, F.K. Damping of surface pressure fluctuations in hypersonic turbulent flow past expansion corners. AIAA J. 1993, 31, 1229–1234. [Google Scholar]
  10. Deck, S.; Thorigny, P. Unsteadiness of an axisymmetric separating-reattaching flow: Numerical investigation. Phys. Fluids 2007, 19, 65–103. [Google Scholar]
  11. Depres, D.; Reijasse, P.; Dussauge, J.P. Analysis of Unsteadiness in Afterbody Transonic Flows. AIAA J. 2004, 42, 2541–2550. [Google Scholar] [CrossRef]
  12. Rui, Z.; Jili, R.; Fang, R. An improved empirical prediction formula for fluctuating pressure on the wall of transonic rotating body. J. Astron. 2016, 37, 1179–1184. [Google Scholar]
  13. Qiang, L.; Xiaolong, C.; Yi, L. Numerical Simulation Research on Transonic BuffetingLoad of Supercritical Airfoil. Tactical Missile Technol. 2018, 192, 50–57. [Google Scholar]
  14. Shujie, R.; Shouyun, Z.; Guirong, Y. Forecast of rocket transonic pressure environmentbased on RANS/NLAS. Solid Rocket Technol. 2011, 34, 418–422. [Google Scholar]
  15. Dawson, J.; Samimy, M. Effects of Expansions on a Supersonic Boundary Layer: Surface Pressure Measurements. AIAA J. 1994, 32, 2169–2177. [Google Scholar] [CrossRef]
  16. Tsutsumi, S.; Takaki, R.; Takama, Y. Hybrid LES /RANS simulations of transonic flowfieldaround a rocket fairing. In Proceedings of the 30th AIAA Applied Aerodynamics Conference, New Orleans, LA, USA, 25–28 June 2012; pp. 25–41. [Google Scholar]
  17. Soni, R.K.; Arya, N.; De, A. Numerical simulation of supersonic separating-reattaching flow through RANS. J. Phys. Conf. Ser. 2017, 822, 12–37. [Google Scholar]
  18. Grilli, M.; Hickel, S.; Adams, N.A. Large-eddy simulation of a supersonic turbulent boundary layer over a compression-expansion ramp. Int. J. Heat Fluid Flow 2013, 42, 79–93. [Google Scholar] [CrossRef]
  19. Kopera, M.A.; Kerr, R.M.; Blackburn, H.M.; Barkley, D. Direct numerical simulation of turbulentflow over a backward-facing step. J. Fluid Mech. 2014, 24–47. [Google Scholar] [CrossRef] [Green Version]
  20. Spalart, P.R.; Jou, W.; Strelets, M.; Allmaras, S. Comments on the feasibility of LES for wings andon a hybrid RANS/LES approach. In Proceedings of the 1st AFOSR International Conferenceon DNS/LES, Ruston, LA, USA, 4–8 August 1997; pp. 181–195. [Google Scholar]
  21. Guilmineau, E.; Deng, G.; Wackers, J. Numerical simulation with a DES approach for automotive flows. J. Fluids Struct. 2011, 27, 807–816. [Google Scholar] [CrossRef]
Figure 1. A model of the expansion probelm. The Mach angles μ 1 = a r c sin 1 M a 1 = 19.4 and μ 2 = a r c sin 1 M a 2 = 15.3 are marked in the figure.
Figure 1. A model of the expansion probelm. The Mach angles μ 1 = a r c sin 1 M a 1 = 19.4 and μ 2 = a r c sin 1 M a 2 = 15.3 are marked in the figure.
Fluids 06 00268 g001
Figure 2. Time-independent check.
Figure 2. Time-independent check.
Fluids 06 00268 g002
Figure 3. Normalized pressure p along the flow direction on the wall ( p 1 is the pressure in the freestream).
Figure 3. Normalized pressure p along the flow direction on the wall ( p 1 is the pressure in the freestream).
Fluids 06 00268 g003
Figure 4. The distribution of normalized fluctuating pressure along the flow direction.
Figure 4. The distribution of normalized fluctuating pressure along the flow direction.
Fluids 06 00268 g004
Figure 5. Characteristics in the Mach number contour map.
Figure 5. Characteristics in the Mach number contour map.
Fluids 06 00268 g005
Figure 6. The U-zone, M-zone and D-zone.
Figure 6. The U-zone, M-zone and D-zone.
Fluids 06 00268 g006
Figure 7. The distribution of measuring points.
Figure 7. The distribution of measuring points.
Fluids 06 00268 g007
Figure 8. The distribution of C p r m s in different layers within boundary layer near expansion corner.
Figure 8. The distribution of C p r m s in different layers within boundary layer near expansion corner.
Fluids 06 00268 g008
Figure 9. PSD and turbulent intensity at positions “d~k”.
Figure 9. PSD and turbulent intensity at positions “d~k”.
Fluids 06 00268 g009
Figure 10. PSD and turbulent intensity at positions “c~m”.
Figure 10. PSD and turbulent intensity at positions “c~m”.
Fluids 06 00268 g010
Figure 11. PSD and turbulent intensity at positions “n~v”.
Figure 11. PSD and turbulent intensity at positions “n~v”.
Fluids 06 00268 g011
Figure 12. PSD and turbulent intensity at positions “w~z”.
Figure 12. PSD and turbulent intensity at positions “w~z”.
Fluids 06 00268 g012
Figure 13. The pressure fluctuation in the U-zone, M-zone and D-zone.
Figure 13. The pressure fluctuation in the U-zone, M-zone and D-zone.
Fluids 06 00268 g013
Figure 14. The distribution of the correlation coefficient between different layers near expansion corner.
Figure 14. The distribution of the correlation coefficient between different layers near expansion corner.
Fluids 06 00268 g014
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, L.; Wu, Z.-N. Numerical Simulation of Pressure Fluctuation near an Expansion Corner in a Supersonic Flow of M = 3.01. Fluids 2021, 6, 268. https://0-doi-org.brum.beds.ac.uk/10.3390/fluids6080268

AMA Style

Zhang L, Wu Z-N. Numerical Simulation of Pressure Fluctuation near an Expansion Corner in a Supersonic Flow of M = 3.01. Fluids. 2021; 6(8):268. https://0-doi-org.brum.beds.ac.uk/10.3390/fluids6080268

Chicago/Turabian Style

Zhang, Lei, and Zi-Niu Wu. 2021. "Numerical Simulation of Pressure Fluctuation near an Expansion Corner in a Supersonic Flow of M = 3.01" Fluids 6, no. 8: 268. https://0-doi-org.brum.beds.ac.uk/10.3390/fluids6080268

Article Metrics

Back to TopTop