1. Introduction
The pressure fluctuation of the flow near the surface of aircraft may induce damage-causing structural vibration and bothering noise environment that affects the normal operation of airborne instruments [
1]. For this reason the problem of pressure fluctuation for an airfoil has received considerable interests.
For subsonic flow, the wall fluctuating pressure generated around a NACA0012 airfoil with the attack angle of 10
has been studied numerically by Wang and Tian [
2], the frequency property shows that the fluctuating pressure is dominated by the vortex shedding for stationary situation and the flapping for the pitching situation, respectively. For the stationary and low frequency pitching airfoil, the fluctuation of the high frequency band induced by the vortex scattering are crucial, which is decreased by the high frequency pitching. In addition, Awasthi et al. [
3] used the large eddy simulation (LES) method to simulate the characterization of the wall fluctuating pressure, the computations of pressure fluctuation around the airfoil and near the trailing-edge are consistent with the measurements in the low-to-mid frequency band, and the histogram distribution of surface pressure fluctuation within the separation bubble upstream of the airfoil shows serious streamwise difference reflecting the instantaneous nature of the flow in this region.
For transonic flow, Hillenherms and Limberg [
4] has experimentally measured the fluctuating pressure on a pitching airfoil. They found that, on the upper side of the pitching airfoil, maximum fluctuation levels occur for Mach numbers of 0.72 to 0.77 which is not only in the region of the shock but overall on the upper side, since this is likely due to a higher fluctuation level of the incoming flow. Hasan and Alam [
5] used numerical simulation to investigate the fluctuating pressure over a supercritical airfoil for a fixed free-stream Mach number of 0.77 and at angle of attack varying from 2
to 7
. They found that the values of peak pressure fluctuation is increased with an increase of angle of attack, and the location of peak pressure fluctuation is shifted toward the leading edge of the airfoil with increasing angle of attack. Chen, Xu and Lu [
6] used detached-eddy simulation (DES) to study pressure fluctuation at Mach number
and discovered that the velocity of downstream-propagating waves in the separated boundary layer is close to the convection speed of the coherent vortical structures, the power spectral density function (
PSD) properties show that there are various spectral scalings with the frequency in different flow regions.
For supersonic flow, Messiter [
7] combined the method of multiple scales and matched asymptotic expansions to derive the disturbances caused by an oscillating airfoil. They found that the oscillatory part of the pressure at a leading-edge shock wave diminishes fleetly with increasing distance from the edge. Later, fluctuating pressure around airfoil has been studied experimentally by Fleeter [
8] for six different oscillating airfoil surfaces in a supersonic flowfield, and the unsteady pressure fluctuation properties is correlated with Mach number.
Past studies of airfoil pressure fluctuation have focused on the wall surface. In generally, pressure fluctuation propagates at the speed of sound, so it is imaginable that the pressure fluctuation does not vary significantly across the boundary layer. However, at least for academic purpose, it is interesting to investigate whether the variation of pressure fluctuation is indeed negligible across the boundary layer.
For a pure turbulent boundary layer, Panton, Lee and Moser [
9] used direct numerical simulation (DNS) method to investigate the correlation of spatial fluctuating properties in turbulent layers, and the root-mean-square (RMS) pressure fluctuations in the inner and outer layers are correlated for different Reynolds numbers. They discovered that there are small variations of fluctuating pressure across the boundary layer and the outer layer profiles of the fluctuation properties correlate very well for various Reynolds numbers, however, the correlation in inner profiles is not excellent, and the overlap matching area has logarithmic behavior.
Zhang and Wu [
10] studied the fluctuating pressure near the expansion corner in supersonic flow and found that near the corner the fluctuation pressure varies across the boundary layer. Moreover, there exists a region bounded by the characteristics lines that the variation of the pressure fluctuation is pronounced. They identified three zones: U-zone, M-zone and D-zone inside which the variation of the pressure fluctuation exhibits various behaviors, due to that disturbance in supersonic flow propagate along inside the Mach cone.
In this paper, we consider the spatial pressure fluctuation properties inside the supersonic flow region of a NACA0012 airfoil at sonic inflow condition, that is, the upstream flow Mach number is one. The reason to consider this case is that there is a limiting characteristics line in the supersonic zone, upstream and downstream of which the pressure waves propagate in different directions. It is expected that the spatial variation of pressure perturbation may exhibit similar interesting phenomena like in the supersonic expansion corner case, and more undesired properties due to the existence of limiting characteristics line.
Here, in this work, we relate the spatial distribution of the fluctuating pressure to the propagation of the characteristics around NACA0012 airfoil, at sonic inflow. In
Section 2, we will introduce the sonic flow problem and numerical problem. The division of the flow domain into U-zone and D-zone according to the limiting characteristic line will be described in
Section 3. The results of the spatial distribution of fluctuation properties in the two different zones are given in
Section 4, where the pressure distributions properties will be measured by the root-mean-square fluctuating pressure coefficient (
) and power spectral density function (
) [
11]. The special behaviors of spatial distribution of
and
in U-zone and D-zone will be presented in
Section 5. Finally, conclusion will be given in
Section 6.
5. Summary of Fluctuating Pressure Properties
Based on the fluctuation properties displayed in
Section 4, we use
Figure 11 to summarize the special behaviors of the fluctuating pressure near the limiting characteristics in a sonic flow around NACA0012 airfoil. There are two zones, U-zone and D-zone, as introduced in
Section 3.
In the U-zone which is upstream of the limiting characteristics, there is spatial variation of across the boundary and this variation becomes important when approaching the limiting characteristics. In the D-zone which is downstream of the limiting characteristics, the spatial variation of gradually increases when approaching the limiting characteristic line. The spatial variation of is the largest along the limiting characteristic line, showing an interesting property of the limiting characteristic line.
Upstream of the limiting characteristic line, PSD shows spatial variation mainly in the middle-frequency band, downstream of the limiting characteritic line the spatial variation of PSA occurs mainly in the high-frequency band, and along the limiting characteristic line this variation occurs over the entire-frequency band.
It is funny to see the correlation of unsteady pressure at two different layers near the limiting characteristic line. The correlation coefficient (named
) of the transient pressure data at two positions
A and
B is defined by
where
is the transient pressure at a certain point
x;
is the covariance of
and
;
is the variance of
.
The distribution of
near the limiting characteristics calculated by Equation (
6) is shown in
Figure 12, where we displayed
between the wall surface and the buffer layer,
between the buffer layer and log layer, and
between the log layer and the outer layer. It is strange that the correlation between different points is the minimal along the limiting characteristic line. Thus, not only the limiting characteristic line has no interaction with the adjacent characteristic lines, the interaction between various locations along this limiting characteristic line is also the smallest.
We observe that a phenomenon occurs both in the U-zone and D-zone. If point A and point B are far from the limiting characteristics(such as position “a” and “g”), the correlation coefficient is a little bit less than 1. If point A and point B are along the limiting characteristics (position “d”), the correlation coefficient is much lower than 1.
6. Conclusions
In this work, the pressure fluctuation properties near the limiting characteristics in a sonic flow (the mach number of inflow is ) around NACA0012 airfoil is studied using numerical method based on DES. The characteristic lines are obtained using the formula for characteristics and using numerical results of flow as input data. The limiting characteristic line is used to divide the supersonic flow zone into two zones: U-zone and D-zone. U-zone is upstream of the limiting characteristic line and D-zone is downstream of this line.
We studied the pressure fluctuation properties inside the U-zone and D-zone, and along the limiting characteristic line. In the U-zone which lies between the sonic line and the limiting characteristics, the values of differ along the characteristic lines within the boundary layer, and the difference gradually increases approaching the limiting characteristics, the variation of PSD occurs in the mid-frequency band in this region, and the correlation coefficient between two adjacent layers is less than 1.
Along the limiting characteristics, both and PSD change significantly in different layers within the boundary layer and the variation of PSD occurs in the entire-frequency band. The correlation coefficient between two adjacent layers is the smallest.
In the D-zone which lies between the limiting characteristics and shock wave, in the process of leaving the limiting characteristics, the variation of gradually decreases in different layers along the characteristic lines, the difference of PSD occurs in the high-frequency band in this region, and the correlation coefficient between two adjacent layers is less than 1.
Thus, the fluctuation pressure and its spatial variation across the boundary layer is strongly determined by the limiting characteristic line. The fluctuation pressure is the largest along the limiting characteristic line, while the correlation coefficient between two adjacent points is the smallest along the limiting characteristic line. Since the waves upstream and downstream of the limiting characteristic line propagate away from the limiting characteristic line, the present finding raises an unanswered question: why the pressure fluctuation is the largest along this line. Away from the limiting characteristic line, the fluctuation pressure decays. There is a spatial variation of the pressure fluctuation across the boundary layer and this spatial variation is in the mid-frequency band in the U-zone, in the high-frequency band in the D-zone, and in the entire-frequency band along the limiting characteristics line. This study revealed a special behavior of the pressure fluctuation along the limiting characteristic line, which not only enriches our knowledge about transonic flow but also raises the question why the pressure fluctuation is the largest along the limiting characteristic line.