Journal of Applied Mathematics and Physics, 2014, 2, 115-123
Published Online April 2014 in SciRes.
How to cite this paper: Khlopkov, Yu.I., et al. (2014) Computer Modelling of Aerothermodynamic Characteristics for Hy-
personic Vehicles. Journal of Applied Mathematics and Physics, 2, 115-123.
Computer Modelling of Aerothermodynamic
Characteristics for Hypersonic Vehicles
Yuri Ivanovich Khlopkov, Anton Yurievich Khlopkov, Zay Yar Myo Myint
Department of Aeromechanics and Flight Engineering, Moscow Institute of Physics and Technology, Zhukovsky,
Received Dec emb er 2013
The purpose of this work is to describe the suitable methods for aerodynamic characteristics cal-
culation of hypersonic vehicles in free molecular flow and the transitional regimes. Moving of the
hypersonic vehicles at high altitude, it is necessary to know the behavior of its aerodynamic cha-
racteristics for all flow regimes. Nowadays, various engineering approaches have been developed
for modelling of aerodynamics of aircraft vehicle designs at initial state. The engineering method
that described in this paper provides good results for the aerodynamic characteristics of various
geometry designs of hypersonic vehicles in the transitional regime. In this paper present the cal-
culation results of aerodynamic characteristics of various hypersonic vehicles in all range of re-
gimes by using engineering method.
Aerodynamic Characteristics, Computational Aerodynamics, Hypersonic Technology, Rarefied Gas
Dynamics; Engineering Method, Aerodynamics in Transitional Regime
1. Introduction
Theoretical studies of hypersonic flows associated with the creation of “Space Shuttle” to transport people and
cargo into Earth orbit began in the last century. Research had been focused mainly on the department of TsAGI
named after N.Y. Zhukovsky (Central Aerohydrodynamic Institute). Practical work on the creation of aerospace
systems had been instructed engineering centre of the experimental design bureau named after A.I. Mikoyan.
Air Force Research Institute has developed an original concept of space system, which efficiently integrated the
ideas of the aircraft, rocket plane and space object in 1960. The project was called “Spiral” and represented as a
complex system. Powerful hypersonic aircraft (weight 52 tons, length 38 m, wingspan 16.5 m), which was dis-
persed to six times the speed of sound (Mach = 6), then at height of 28-30 km from its back, manned orbital
plane 10 ton (8 m long and 7.4 m span) was supposed to start.
The “Spiral” (Figure 1) was a response to the U.S. space program to create an interceptor reconnaissance
bomber, the X-20 “Dyna Soar” (Figure 2). As shown, the implementation of project “Dyna Soar” is not suc-
cessful as a “Spiral”. In the end both projects have been folded, although at different stages of development [1].
Since 1980 aerospace vehicle programs are developed in many developed countries as the U.S. and the Soviet
Yu. I. Khlopkov et al.
Figure 1. Russian project “Spiral”.
Figure 2. USA project “Dyna Soar”.
Union. For example, in England “HOTOL” (Figure 3), Germany “Zenger” (Figure 4), France “Hermes” (Fig-
ure 5), Japan “Hope” (Figu re 6), China “Shenlong” (Figur e 7) and India “AVATAR” (Figure 8). Some of
them have been folded.
Russia is developing new generation reusable spacecraft programs “Clipper”, “RUS” to deliver crews and
cargo to low Earth orbit and the space station since 2000. The first flight is expected in 2015. United States is
developing a new spacecraft “Orion” to deliver crew and cargo to low Earth orbit and back. First full-size test
flight is in 2014, the first flight to the moon planned in 2020.
The DARPA (Defense Advanced Research Projects Agency) is developed “FALCON” (Force Application
and Launch from CONtinental United States) since 2003. The Falcon Hypersonic Technology Vehicle (HTV-2)
is a multiyear research and development effort to increase the technical knowledge base and advance critical
technologies to make long-duration hypersonic flight a reality. Falcon HTV-2 is an unmanned, rocket-launched,
maneuverable aircraft that glides through the Earth’s atmosphere at incredibly fast speedsMach 20.
With the development of space and rocket technologies, it is required the reliable methods on the aerodyna-
mic and aerothermodynamic characteristics modelling of hypersonic vehicles in the whole range of flow regi-
mes, i. e., from the continuum flow regime up to the free-molecular regime. During de-orbiting, the spacecraft
passes through the free molecular, then through the transitional regime and the finalized flight is in the contin-
uum flow.
As we have known that for flight in the upper atmosphere, where it is necessary to take into account the mole-
cular structure of a gas, kinematics models are applied, in particular, the Boltzmann equation and corresponding
numerical methods of simulation. In the extreme case of free-molecular flow, the integral of collisions in the
Boltzmann equation becomes zero, and its general solution is a boundary function of distribution, which remains
constant along the paths of particles [2]. While aircraft are moving in a low atmosphere, the problems are re-
duced to the problems that can be solved in the frame of continuum theory or, to be more precise, by application
of the Navier-Stokes equations and Euler equations [2,3]. On the transition interval between the free molecular
and continuum regimes numerical methods of solving the Boltzmann equation and its model equations are being
used with success [4].
Yu. I. Khlopkov et al.
Figure 3. HOTOL.
Figure 4. Zenger.
Figure 5. Hermes.
Figure 6. Hope.
Figure 7. Shenlong.
Figure 8. AV A T AR .
Yu. I. Khlopkov et al.
To correctly simulate hypersonic flows, the flows must be understood and modeled correctly and this is more
true than in the numerical simulation of hypersonic flows. Hypersonic must be dominated by an increased un-
derstanding of fluid mechanics reality and an appreciation between reality and the modeling of that reality [5].
The benefits of numerical simulation for flight vehicle design are enormous: much improved aerodynamic shape
definition and optimization, provision of accurate and reliable aerodynamic data and highly accurate determina-
tion of thermal and mechanical load [6]. Multi-parametric calculations can be performed only by using an ap-
proximation engineering approach. Computer modeling allows to quickly analyze the aerodynamic characteris-
tics of hypersonic vehicles by using theoretical and experimental research in aerodynamics of hypersonic flows.
Nowadays, the basic quantitative tool for study of hypersonic rarefied flows is direct simulation Monte Carlo
method (DSMC) [3] [7]. DSMC method required large amount of computer memory and unreasonable expen-
sive at the initial stage of vehicle design and trajectory analysis. The solution for this problem is the approximate
local engineering methods [8-11]. The Monte Carlo method remains the most reliable approach, together with
the local engineering methods that provide good results for the global aerodynamic coefficients. The early work
of [3] indicated that local engineering methods could have significant effect on aerodynamic characteristics of
various hypersonic vehicles. It is natural to create engineering methods, justified by cumulative data of experi-
mental, theoretical and numerical results, enabling the prediction of aerodynamics characteristics of complex
bodies in the transitional regime [3].
The purpose of this work is to calculate aerodynamic characteristics of perspective space vehicle “Clipper”
and hypersonic technology vehicle “Falcon HTV-2” by using engineering method. This engineering method is
suitable to calculate with taking into account the various Reynolds number, and provide good results for various
hypersonic vehicle designs.
2. Calculation Methods
2.1. Method for Hypersonic Aerodynamics in the Free Molecular Flow Regime
The problem of determination of the aerodynamic characteristics in free molecular flow is set in a usual way. In
the case of hypersonic flow, the values of pressure p and tangential τ forces and heat flux q for the element of a
surface look as below
( )
ργ 1
2σ2sin σπ sin
τσ2sin cos
ργ 1
σ1 2sin
= −
where στaccommodation coefficient, γspecific heat ratio, Tw, Tsurface temperature and flow tempera-
ture respectively . The experiment provides the most reliable results for the formulation of a connection between
the flow as like the oncoming particles characteristics and the reflected particles through the coefficients of ac-
commodation. The results of aerodynamic characteristics of semi-sphere, cone with the various semi-angles by a
vertex, blunted semi-cones, and blunted semi-cones with wings are presented in [12].
2.2. Method for Hypersonic Aerodynamics in the Transitional Regime
In this work we use the expressions for the elementary pressure forces and friction forces are applied in the form
described in [3,11].
where, coefficients p0, p1, τ0 (coefficients of the flow regime) are dependent on the Reynolds number Re0 =
ρVL/µ0, in which the viscosity coefficient µ0 is calculated at stagnation temperature T0. Except Reynolds
number the most important parameter is the temperature factor Tw/T0, where T0, Tw are the stagnation tempera-
sin θ +sinθpp p=
τ = τ sinθ cosθ
Yu. I. Khlopkov et al.
ture and surface temperature respectively.
The dependency of the coefficients of the regime in the hypersonic case must ensure the transition to the
free -molecular values at Re0→0, and to the values corresponding to the Newton theory, methods of thin tangent
wedges and cones, at Re0→∞. On the basis of the analysis of computational and experimental data, the empirical
formulas are proposed [15-17]
Re10 Re
1.8(1 )
= −
whe r e h is a relative lateral dimension of the apparatus, which is equal to the ratio of its height to its length.
The technique proved to be good for the calculation of hypersonic flow of convex, not very thin, and spatial
bodies. The calculation fully reflects a qualitative behavior of drag force coefficient CD as a function of the me-
dium rarefaction within the whole range of the angles of attack, and provides a quantitative agreement with ex-
periment and calculation through the Boltzmann equation with an accuracy of 5%. On the accuracy of the rela-
tion of the locality method can be said that they are applied with the smallest error in the case of the bodies that
are close to being spherical, and are not applied in the case of very thin bodies, when the condition is M sin θ >>
This locality method for calculation of aerodynamic characteristics of the bodies in the hypersonic flow of ra-
refied gas in the transitional regime gives a good result for CD for a wide range of bodies, and a qualitatively
right result for lift force coefficient CL. In this case, it is necessary to involve more complete models that take
into account the presence of the boundary layer [3]. In early papers [13,15-19] described the results of aerody-
namic characteristics of various hypersonic vehicles by using this method.
3. Results
The results of the calculation of the coefficients of drag force CD, lift force CL, pitching moments MZ with value
of angle of attack α from 0 to 90 deg for Russian perspective space vehicle “Clipper, TsAGI model” (Figure 9)
[14] and USA perspective hypersonic technology vehicle “Falcon HTV-2” (Figure 10) are presented.
The calculation has been carried out through the method described on the above section within the range of
angles of attack α from 0 deg up to 90 deg with a step of 5 deg. The parameters of the problem are the following:
ratio of heat capacities
= 1.4; temperature factor Tw/T0 = 0.01; Reynolds number Rе0 = 0, 10, 100, 10000; ve-
locity ratio = 15.
The dependencies of CD(
), CL(
) and MZ(
) are presented in Figures 11-13. It can be seen from these results
that when the Reynolds number increased, the drag coefficients CD of vehicle diminished which can be ex-
Figure 9. Geometry view of
space vehicle “Clipper”.
[(2 α)] /
p ppppz
∞∞ ∞
=+ −−
2 1/2
τ3.7 2[6.88exp(0.00720.000016)]R RR
Re 44
= +
Yu. I. Khlopkov et al.
Figure 10. Geometry view of hyper-
sonic vehicle “Falcon HTV-2”.
Figure 11. Drag coefficients CD for space vehicle “Clipper”.
Figure 12. Lift coefficients CL for space vehicle “Clipper”.
Yu. I. Khlopkov et al.
plained by the decrease of normal and tangent stresses. At high Reynolds number Re0 ≥ 106, characteristics al-
most not changed. The dependency CL(
) is increased at high Reynolds number which can be explained by the
decrease of normal and tangent stresses.
The values of MZ are quite sensitive to the variation of Re0. MZ changes its sign less than zero at Re0 ~ 102. At
Re0 ~ 104, the value of Mz = 0.03 at the angle of attack is reached at α 40 deg. Results by using local engi-
neering method are compared with the results obtained by DSMC and Newtonian methods.
The dependencies of CD(α), CL(α) and MZ(α) of hypersonic vehicle “Falcon HTV-2” are presented in Figures
14-16. In Figure 14, it can be seen with the increasing of Reynolds number, the drag coefficients CD of vehicle
decreased. It can be explained that by the decrease of normal and tangent stresses. Drag coefficients CD of Fal-
con more than Clipper. The dependency CL(α) is increased, and the value is reached to 0.54 at Re0 ~ 104. The
values of MZ are quite sensitive to the variation of Re0, changes its sign at α ~ 5 deg.
Figure 13. Pitching moment coefficients MZ for space
vehicle “Clipper”.
Figure 14. Drag coefficients CD for hypersonic vehicle
“Falcon HTV-2”.
Yu. I. Khlopkov et al.
Figure 15. Lift coefficients CL for hypersonic vehicle “Falcon
HTV -2”.
Figure 16. Pitching moment coefficients MZ for hypersonic
vehicle “Falcon HTV-2”.
4. Conclusions
Many analytical and numerical methods to estimate the aerothermodynamics of hypersonic vehicles were ap-
peared with the development of space launch technologies. This paper presents different methods to calculate
aerodynamic characteristics of various perspective hypersonic vehicles in rarefied gas flows. With the use of this
engineering method, it had been calculated aerodynamic characteristics for various hypersonic vehicle design in
rarefied gas flow. Results are compared with the DSMC and traditional Newtonian method. Methods which de-
scribed above give good results and suitable to calculate aerodynamics for various hypersonic vehicle designs.
The reported study was partially supported by RFBR (research project No. 11-07-00300-а).
Yu. I. Khlopkov et al.
The support by RFBR (research project No. 11-07-00300-а) is cordially appreciated by the authors.
[1] Khlopkov, Yu . I., Chernyshev, S.L., Khlopkov, A.Yu. and Myint , Z.Y.M. (2013) Introduction to Specialt y. High-Speed
Aerial Vehicles. MIPT, Moscow.
[2] Koga n, N.M. (1969) Rarefied Gas Dynamic. P len um, New York. -6381-9
[3] Belo t serko vskii , O.M. and Khlopko v, Yu .I. (2010) Monte Carlo Methods in Mechanics of Fluid and Gas. World Scien-
tific Publishing Co. N-Y, London, Singapore, Beijing, Hong Kong.
[4] Gu s e v , V.N. (1993) High-Altitude Aerothermodynamics. Journal of Fluid Dynamics, 28, 269-276.
[5] Neuman n, R.D. (19 88 ) Missions and Requi rements. Special Course Aerothermodynamics of Hypersonic Vehicles.
AGARD Report 761, Neuilly sur Seine, France.
[6] Hirsch el, E.H. (2005) Basics of Aerothermodynamics. Progress in Astronautics and Aeronautics, AIAA, Springer-
Verlag, Berlin/Heidelberg/New York.
[7] Bird, G. A. (1994) Molecular Gas Dynamics and the Direct Simulation of Gas Flows. Oxford University Press.
[8] Kotov, V., Lych ki n, E., Resheti n, A. and Shelkonogov, A. (1982) An Approximate Method of Aerodynamics Calcula-
tion of Complex Shape Bodies in a Transition Region. Proc. of 13th International Conference on Rarefied Gas Dy-
namics, Plenum Press, New York, 1, 487-494.
[9] Morsa, L., Zuppardi, G., Schettino, A. and Votta, R. (2010) Analysis of Bridging Formulae in Transitional Regime.
Proc. of 27th International Symposium on Rarefied Gas Dynamics, Pacific Grove, California, 10-15 July.
[10] Vashchenkov, P.V., Ivano v, M.S. and Kr ylov, A.N. (2010) Numerical Simulations of High-Altitude Aerothermody-
namics of a Promising Spacecraft Model. Proc. of 27th International Symposium on Rarefied Gas Dynamics, Pacific
Grove, California, 10-15 July, 1337-1342.
[11] Galki n, V.S., Er ofeev, A.I. and Tolstykh, A.I. (1977) Approximate Method of Calculation of the Aerodynamic Char-
acteristics of Bodies in a Hypersonic Rarefied Gas. Proceedings TsAG I , 18 33 , 6-10.
[12] Khlopkov, Yu. I. (2006) Statistical Modeling in CFD. MIPT, Moscow.
[13] M yin t, Z.Y.M. and Khlopkov, A.Yu. (2010 ) Aerodynamic Characteristics of an Aircraft with a Complex Shape Taking
into Account the Potential of Molecular Flow Interaction with a Surface. TsAGI Science Journal, 41, 551-566.
[14] Vagan ov, A.V., Drozdov, S., Kosykh, A.P., Nersesov, G.G., Chel yshe va, I.F. and Yumashev, V.L. (2009) Numerical
Simulation of Aerodynamics of Winged Reentry Space Vehicle. TsAGI Science Journal, 40, 131-149.
[15] Khlopkov, Yu.I., Myint, Z.Y.M. and Khlopkov, A.Yu. (2013) Aerodynamic Investigation for Prospective Aerospace
Vehicle in the Transitional Regime. International Journal of Aeronautical and Space Sciences, 14, 215-221.
[16] M yin t, Z.Y.M., Khlopkov, Yu.I. and Khl opko v, A.Yu. (2013 ) Aerothermodynamics Investigation for Future Hyper-
sonic Aerospace Sys tem. Proc. of 4th International Conference on Science and Engineering, Yangon, Myanmar, 9-10
[17] Khlopkov, Yu.I., Zharov, V.A., Myin t, Z.Y.M. and Khlopkov, A.Yu. (20 13 ) Aerodynamic Characteristics Calculation
for New Generation Space Vehicle in Rarefied Gas Flow. Universal Journal of Physics and Application, 1, 286-289.
[18] M yin t, Z.Y.M., Khlopkov, Yu.I., Khlopkov, A.Yu. and Polyakov , M.S. (2013) Computational Analysis of Aerod y-
namic Characteristics for Hypersonic Vehicles. International Journal of Applied and Fundamental Research.
[19] Khlopkov, Yu.I., Chernyshev, S.L., Zharov, V.A., Myint, Z.Y.M., Khlopkov, A.Yu., Polyakov, M.S. and Zin, K. (2014)
Modern Trends in the Development of Reusable Aerospace System. Asian Journal of Applied Sciences, 2.