Computational Simulation of Shock Oscillation Around A
Computational Simulation of Shock Oscillation Around A
Computational Simulation of Shock Oscillation Around A
and
Yasushi WATANABE‡
Japan Aerospace Exploration Agency, Tokyo, Japan
Nomenclature
h = throat height
hc = capture height
he = exit height
lFP = length from the flow plug head to the rear end of the intake
lr = ramp length
M = Mach number
P0 = stagnation pressure
PRth = theoretical total pressure recovery
Re = Reynolds number
T0 = stagnation temperature
I. Introduction
N EW supersonic transports (SST) are under development now. Number of environmental and economical
problems is remained to accomplish the development. From mid ‘90s, several research projects on the new
SST have been carried out in Europe, U.S.A., Russia, and Japan. Computational fluid dynamics (CFD) is a key
technique for designing the new SST.
*
Graduate Student, Department of Mechanical Systems Engineering, Koganei, Tokyo 184-8588, Japan.
†
Associate Professor, Department of Mechanical Systems Engineering, Koganei, Tokyo 184-8588, Japan.
‡
Researcher, Aviation Program Group, Mitaka, Tokyo 182-8522, Japan, Member AIAA.
1
American Institute of Aeronautics and Astronautics
Copyright © 2006 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.
We investigated the flow around supersonic air intakes by a CFD technique. The intake is an important part of
supersonic aircraft engines. The intake has a subsonic duct with a supersonic compression wedge. We solved both
the internal and external flow of a compressible viscous fluid, which is a challenging topic on current CFD
techniques.
In the present paper, we focused our attention to an intake ‘buzz’ flow.1 All the supersonic intakes exhibit a
marked instability of flow in subcritical operation below some value of the flow ratio, in the form of oscillation of
the shock waves, which is colloquially known as ‘buzz’, generated at the supersonic compression ramp.
Two remarkable buzz flows was caused by (1) inflow of shear layer (Ferri type2), (2) separation of the flow from
the ramp surface (Dailey type3). Dailey also observed a high frequency and small-amplitude oscillation under low
mass flow rate case, which is caused by acoustic resonance of the subsonic duct.4, 5 Recently, computational fluid
dynamics was applied to simulate the buzz flow.6, 7 However, since the flow is complicated, the mechanism has not
been fully understood yet.
We calculated the shock oscillation of a two-dimensional external-compression supersonic air intake. Especially,
we paid our attention to influence of a slit at the entrance of the subsonic duct. The slit was commonly placed for
bleed of the air, which is effective to reduce the shock oscillation. In addition, we examined the length of the duct
using two models having different lengths.
We found that high-frequency pressure oscillation was generated by the periodic separation on the wedge surface
due to the suction from the slit. Low-frequency component, on the other hand, was caused by the acoustic resonance
of the subsonic duct. The high-frequency component disappeared when the diffuser length was long.
A. Wedge
The angle of the wedge (ramp) is 10 degrees. The uniform stream (Mach number, M=1.64) passes an oblique
shock wave generated by the ramp. Theoretically, the flow is decelerated to M=1.295 and the recovery of total
pressure PRth=0.989. Subsequently, the flow is decelerated again to M=0.786 and PRth=0.979 by a normal shock
wave at the entrance of the subsonic diffuser (throat).
(a) Configuration
Figure 1. Two-dimensional model of an external-compression supersonic air-intake.
2
American Institute of Aeronautics and Astronautics
(b) Shock systems
Figure 1. Two-dimensional model of an external-compression supersonic air-intake.
B. Subsonic Diffuser
The characteristic length of subsonic diffuser is the height of the entrance (throat) h.
The length of subsonic diffuser is 10h. We designed the cross sectional area, which is corresponding to the
height in this two-dimensional model, along the streamline as follows.
The diffuser is separated into forward and rear parts, whose lengths are 3.5h and 6.5h, respectively. The forward
part is a bend. It has a constant deceleration rate; the Mach number gradient is dM/dh = -0.08. The turning angle is -
10 degree to align the direction of the flow parallel to the uniform external flow at the end of the forward part. The
rear part is a straight diffuser, in which the Mach number gradually decreases to a constant value (M=0.4).
3
American Institute of Aeronautics and Astronautics
C. Flow Plug
A triangle flow plug is used to control mass flow through the intake. The vertex angle of the flow plug is 20
degrees.
The cross sectional area of the diffuser exit can alter by insert length of the plug lFP. According to the insert
length, we can operate the intake as critical, subcritical and supercritical conditions, as shown in Fig. 2. The
supercritical condition is obtained when the insert length lFP is smaller than the length at the critical condition lFPC.
Oppositely, the subcritical condition is achieved when lFP is larger than lFPC. Remarkable shock oscillation is
observed in the subcritical condition.
A. Flow Solver
We used UPACS (Unified Platform for Aerospace Computational Simulation) developed by JAXA (Japan
Aerospace Exploration Agency) 9, 10 as the flow solver. The governing equations were Navier-Stokes equations. The
Roe scheme was used for spatial difference. The Euler implicit method was used for the time advancement. No
turbulent flow model was incorporated in the present calculation.
B. Numerical Grid
A numerical grid around the intake model is shown in Fig. 3. It was a structured grid and was divided into four
or five blocks. The initial grids were generated by an algebraic method. Then an elliptic method was applied to
improve the diagonality of the final grid. The total numbers of the grid points are from 295,458 to 442,809 points,
which depend on the geometry of the intake model. The minimum non-dimensional spacing of the grid is 2.0×10-5,
which is smaller than the thickness of viscous layer (laminar sublayer). A multi-grid method is used to connect the
blocks.
4
American Institute of Aeronautics and Astronautics
IV. Numerical Results and Discussion
5
American Institute of Aeronautics and Astronautics
(a) Numerical simulation (b) Experiment8
6
American Institute of Aeronautics and Astronautics
Table 1. Averaged mass flow ratio.
On the other hand, substantial differences are found between the standard and ‘without slit’. The averaged shock
position in the ‘without slit’ model is closer to the tip of the wedge than that of the standard one. This is primary
caused by the difference of total mass flow rate captured by the intake, which is the sum of the mass flows through
the subsonic duct and the slit. The total mass flow of ‘without slit’ model is much lower than that of the standard
model.
The dominant frequency of the shock oscillation was obtained from these sequential images: It is 4.2 kHz in Fig.
7(a), 4.6 kHz in Fig. 7(b), and 600 Hz in Fig. 7(c). We can find a large discrepancy between the standard and the
‘without slit’ models. The cause of this discrepancy is discussed in next subsection.
7
American Institute of Aeronautics and Astronautics
0 0
10 10
pressure data point: pressure data point:
4084[Hz] :b3 :z3
:b4 :z4
–2 –2
10 10
Normalized amplitude
Normalized amplitude
–4 –4
10 10
–6 –6
10 10
–8 –8
10 10
Standard type Standard type
Normalized amplitude
–4 –4
10 10
–6 –6
10 10
–8 –8
10 10
Long diffuser type Long diffuser type
Normalized amplitude
639[Hz]
–4 –4
10 10
–6 –6
10 10
–8 –8
10 10
Non–slit type Non–slit type
8
American Institute of Aeronautics and Astronautics
E. Discussion
Frequency spectra indicate that the high-frequency oscillation (4 kHz) is observed in the standard model as well as
the ‘long diffuser’ one. This oscillation is caused by the periodic suction through the slit of the air in the flow
separation on the wedge surface.
Figure 9 shows iso-Mach contour obtained from the numerical simulation. Figure 9(a) corresponds to the
standard model. In Fig. 9(a), a large separation of flow is observed on the surface of the wedge. The air in the
separation region is sucked into the slit due to the pressure difference between inside and outside of the intake. Just
after the suction, the mass flow into the subsonic diffuser becomes maximum value so that the normal shock is
placed its rear limit position. Then the separation is gradually developed behind the normal shock, which leads to
decrease of the mass flow into the subsonic diffuser. The position of normal shock is also gradually close to the tip
of the wedge due to the decrease of the mass flow into the subsonic diffuser. When the position reaches the front
limit, the air in the flow separation is sucked again.
Figure 9(b) is the iso-Mach contour of the ‘without slit’ model. The separation region of the ‘without slit’ model
is much larger than that of the standard model. In this case, periodic phenomena may be caused by the acoustic
resonance of the subsonic duct. The lowest resonant frequency is defined as
( )
f = 1 − M 2 c 4L (1)
where L, M, and c denote the length, averaged Mach number and speed of the sound in the diffuser, respectively.6
In the present calculation, f = 642 Hz with L = 76 mm, M = 0.7 and c = 380 m/s. This value is close to the dominant
frequency shown in Fig. 8(c).
V. Concluding remarks
The shock oscillation around supersonic air intakes was examined. Especially, we focused our attention to the
effect of the natural bleed through the slit. The principal conclusions are as follows.
1) Natural bleed causes high-frequency shock oscillation. This is due to the periodic suction of the air in flow
separation on the surface of the supersonic ramp (wedge).
2) The high-frequency component disappears in the subsonic diffuser in the case of ‘long diffuser’ model.
3) Natural bleed suppresses the pressure fluctuation at the end of the subsonic diffuser in the case of constant
mass flow ratio of the subsonic diffuser (0.7).
References
1
Seddon, J., and Goldsmith, E. L., Intake Aerodynamics, 2nd ed., Blackwell Science, 1999, pp. 245-265.
2
Ferri, A., and Nucci, L.M., “The Origin of Aerodynamic Instability of Supersonic Inlets at Subcritical Conditions, ” NACA
RM-L50K30, 1951.
3
Dailey, C. L., “Supersonic Diffuser Instability,” Journal of the Aeronautical Sciences, Vol.22, No.11, 1955, pp.733-749.
4
Sterbentz, W. H., and Evvard, J. C., “Criterions for Prediction and Control of Ram-jet Flow Pulsations,” NACA TN-3506,
1955.
9
American Institute of Aeronautics and Astronautics
5
Trimpi, R. L., “An Analysis of Buzzing in Supersonic Ram Jets by a Modified One-dimensional Nonstationary Wave
Theory,” NACA TN-3695, 1956.
6
Newsome, R. W., “Numerical Simulation of Near-critical and Unsteady Subcritical Inlet Flow Fields,” AIAA-Journal, Vol.
22, No. 10, 1984, pp. 1375-1379.
7
Lu, P. J., and Jain, L. T., “Numerical Investigation of Inlet Buzz Flow,” Journal of Propulsion and Power, Vol.14, 1998,
No.1, pp.90-100.
8
Nishizawa, U., Yamamoto, S., Kameda, M., and Watanabe, Y., “Shock Oscillation of an External-compression Supersonic
Air Intake,” Proceedings of 2005 JSASS-KSAS Joint International Symposium on Aerospace Engineering, Nagoya, Japan, 2005,
pp. 109-112.
9
Yamazaki, H., Enomoto, S., and Yamamoto, K., “A Common CFD Platform UPACS,” High Performance Computing: 3rd
International Symposium, ISHPC 2000, edited by M. Valero, K. Joe, M. Kitsuregawa, H. Tanaka, Lecture Notes in Computer
Science, Springer, Berlin, 2000, pp.182-190.
10
Takaki, R., Yamamoto, K., Yamane, T., Enomoto, S., and Mukai, J., “The Development of the UPACS CFD
Environment,” High Performance Computing: 5th International Symposium, ISHPC 2003, edited by A. Veidenbaum, K. Joe, H.
Amano, and H. Aiso, Lecture Notes in Computer Science, Springer, Berlin, 2003, pp.307-319.
10
American Institute of Aeronautics and Astronautics