Aerodynamics Investigation for Prospective Spacecraft Vehicle in the Transitional Regime
Aerodynamics Investigation for Prospective Spacecraft Vehicle in the Transitional Regime
International Journal of Aeronautical and Space Sciences. 2013. Sep, 14(3): 215-221
Copyright ©2013, The Korean Society for Aeronautical Space Science
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License ( which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : July 07, 2013
  • Accepted : September 06, 2013
  • Published : September 30, 2013
Export by style
Cited by
About the Authors
Khlopkov Yuri Ivanovich
Zay Yar Myo Myint
Khlopkov Anton Yurievich

The basic quantitative tool for the study of hypersonic rarefied flows is the direct simulation Monte Carlo method (DSMC). The DSMC method requires a large amount of computer memory and performance and is unreasonably expensive at the first stage of spacecraft design and trajectory analysis. A possible solution to this problem is approximate engineering methods. However, the Monte Carlo method remains the most reliable approach to compare to the engineering methods that provide good results for the global aerodynamic coefficients of various geometry designs. This paper presents the calculation results of aerodynamic characteristics for spacecraft vehicles in free molecular, the transitional and continuum regimes using the local engineering method. Results and methods would be useful to calculate aerodynamics for new-generation hypersonic vehicle designs.
1. Introduction
Its aerospace system is the main indicator of the international power of any country concerned with guaranteeing security, the standard of scientific research and technical development. When designing the launch system of cargo and people into orbit, we used high-altitude aviation and rocket technology. In the future, almost all the projects are to design new space systems for multiple use - to launch loads into orbit and return from it. First, this is related to the use of winged spacecraft for descent and landing on usual aerodromes. To estimate the aerodynamic loads and to predict the landing area, one should know at the stage of initial design, the coefficients of aerodynamic forces and moments of spacecraft for varying temperatures, densities and flight velocities for numerous possible descent trajectories.
The development of spacecraft and rocket technologies requires reliable data on the aerodynamic and aerothermodynamic characteristics for hypersonic speed in an entire range of flow regimes, i.e., from the continuum flow regime to the free-molecular regime. During de-orbiting, the spacecraft passes through the free molecular, then through the transitional regime and the final flight is in the continuum flow.
It is well known that for flight in the upper atmosphere, where it is necessary to take into account the molecular 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 [1] . While aircraft are moving in a low atmosphere, the problems are reduced 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. It is natural to create engineering methods, justified by cumulative data of experimental, theoretical and numerical results, enabling the prediction of aerodynamics characteristics of complex bodies in the transitional regime [2] .
The difficulty in the development of a hypersonic vehicle is caused by quite a number of problems of modeling fullscale flight conditions in wind tunnels. The analysis of the aerodynamics and aerothermodynamic characteristics of a hypersonic vehicle at high-altitude requires numerous numerical calculations [3 - 5] .
To correctly simulate hypersonic flows, the flows must be understood and modeled correctly and this truer than in the numerical simulation of hypersonic flows. Hypersonics must be dominated by an increased understanding of fluid mechanics reality and an appreciation between reality and the modeling of that reality [6] . 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 determination of thermal and mechanical load [7] . Multi-parametric calculations can be performed only by using an approximation engineering approach. Computer modeling allows rapid analysis of the aerodynamic characteristics of hypersonic vehicles using theoretical and experimental research in aerodynamics of hypersonic flows.
The basic quantitative tool for the study of hypersonic rarefied flows is the direct simulation Monte Carlo method (DSMC) [8 , 2] . The DSMC method requires a large amount of computer memory and performance and is unreasonably expensive at the initial stage of spacecraft design and trajectory analysis. Approximate engineering methods pose a possible solution to this problem. These methods allow us to calculate aerodynamic characteristics of vehicles for multiple variants of free-stream parameters within a reasonable time and to refine the results at the most important segments of the flight trajectory by the DSMC method. The early work of [2] indicated that the local engineering method could have a significant effect on aerodynamic characteristics of various hypersonic vehicles.
The purpose of this work is to create an engineering program for aerodynamic characteristics calculation of the perspective hypersonic vehicles. This method is suited to calculate for hypersonic vehicles, taking into account the influence of Reynolds number, and provide good results for various vehicle shapes.
2. Calculation Methods for Hypersonic Aerodynamics
Difficulties of solving aerodynamic problems of overflow on bodies of rarefied gas stimulated the development of engineering semi-empirical methods based on the accumulated experimental data. The most suitable method to compute aerodynamic forces of a spacecraft relies on bridging formulae. When modeling the natural conditions, it is necessary to consider the basic similarity criteria. The hypersonic aerodynamic coefficients most commonly used parameters are:
Mach number M
PPT Slide
Lager Image
Knudsen number Kn
PPT Slide
Lager Image
Reynolds number Re
PPT Slide
Lager Image
where, λ - the mean free path, L ref - reference length of vehicle, μ - the viscosity coefficient. In hypersonic flow calculation more efficiently accepted to use Reynolds number than Knudsen number.
The famous Newtonian sine squared law indicated that the force varies as the square of the sine of the deflection angle. More than half century later it was indicated that Newton’s sine-squared law was not very accurate and, indeed, the preponderance of fluid dynamic experience up to the present day confirms this finding, the exception to this is the modern world of hypersonic aerodynamics [9] . The pressure coefficient at all forward facing points of the body is
  • Cp= 2sin2θ
where, θ is the local inclination of the surface to the free stream.
An alternative representation of the pressure coefficient for hypersonic flow, replacing the factor 2 by the appropriate value of the stagnation-point pressure coefficient at the flow condition of interest, which is known as modified Newtonian flow theory [10]
  • Cp=Cp,mod2sin2θ
where Cp ,mod could be found experimentally or by using the Rayleigh’s pitot formula [11 , 12] .
PPT Slide
Lager Image
Where, p stag - the stagnation pressure, γ - specific heat ratio.
For the last 50 years, the introduction of modified
Newtonian theory continues to be used often to generate approximation for the pressure acting on configurations in hypersonic flows [13 - 15] .
In this work, the expressions for the elementary pressure forces and friction forces are applied in the form described in [2 , 16 , 17] . The local formulae take directly into account the geometry of the vehicle and calculate pressure and skin friction distribution on the body surface.
  • p=p0sin2θ+p1sinθ
  • τ = τ0sinθ cosθ
here, p 0 , p 1 , τ 0 (coefficients of the flow regime) are dependent on the Reynolds number Re 0 = ρ V L 0 ( T 0 ), in which the viscosity coefficient μ 0 is calculated at stagnation temperature T 0 . The most important parameter is the temperature factor Tw / T 0 . T w is the surface temperature.
The dependency of the coefficients of the regime in the hypersonic case must ensure the transition to the freemolecular values at Re 0 →0, and the values corresponding to the Newton theory, methods of thin tangent wedges and cones at Re 0 . By the analysis of computational and experimental data, the empirical formulas are proposed, as follow
  • p0=p∞+ [p∞(2-αn)-p∞]p1/z,
  • p1=zexp[-(0.125+0.078tw)Re0϶фф],
PPT Slide
Lager Image
PPT Slide
Lager Image
  • Re0϶фф= 10-mRe0,m= 1.8(1-h)3
Where, h is a relative lateral dimension of the body, which is equal to the ratio of its height to its length.
PPT Slide
Lager Image
Calculation diagram of prospective spacecraft vehicle “Clipper”
The technique was proved to be good for the calculation of hypersonic flow for convex and not very thin bodies.
The calculation fully reflects a qualitative behavior of drag force coefficient CD as a function of the medium rarefaction within the whole range of the angles of attack, and provides a quantitative agreement with the experiment and calculation through the Boltzmann equation with an accuracy of 5%. Considering the accuracy of the relation of the locality method, it 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 θ˃˃1 [2] .
Thus, the locality method in the calculation of aerodynamic characteristics of the bodies in the hypersonic flow of rarefied gas in the transitional regime gives a good result for CD for a wide range of bodies, and a qualitatively sufficient result for lift force coefficient CL . A comparison of the experimental and calculated results (dashed curve) of drag coefficient CD for the sphere using the local engineering method (M = 5.15-10, Tw / T 0 = 1) is performed in fig. 2 [16] . It can be seen that the influence is small. Figure 3 presents the results of drag coefficients CD (α) for the space shuttle with the use of this method.
PPT Slide
Lager Image
Drag coefficient CD for a sphere at various Reynolds number [16]
PPT Slide
Lager Image
Drag coefficient of Space Shuttle at the various stages of flight trajectories [2]
3. Computational Results and Discussion
The coefficients of drag force CD , lift force CL and pitching moment MZ ware calculated according to equations, as below
PPT Slide
Lager Image
L ref , S ref – references length and surface; F i M – resultant
PPT Slide
Lager Image
Drag coefficients CD for “Clipper” (tw = 0.01)
PPT Slide
Lager Image
Lift coefficients CL for “Clipper” (tw = 0.01)
force acting on the vehicle and moment, respectively.
The calculation has been performed through the method described in the previous 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 t w = Tw / T 0 = 0.01; velocity ratio M =15, Reynolds number Re 0 = 0, 10, 10 2 , 10 4 .
The dependencies of the drag force coefficient CD , lift force coefficient CL and pitching moment coefficient MZ on angle of attack α from 0 to 90 deg for the prospective spacecraft vehicle “Clipper, TsAGI model” ( fig. 1 ) [17 , 18] are presented in figs 4 - 9 . Comparisons of the numerical values of CD , CL , MZ for the spacecraft vehicle using the DSMC method, Newton’s method and local engineering method are also presented in figs 4 - 6 . Thus, the dependencies of the CD , CL , MZ on the Reynolds number are significant in all ranges of angle of attack. This comparison led to the corrections of the computation of the aerodynamic characteristics coefficients.
From these results, it can be seen that the influence of results by DSMC, Newton’s method and local engineering method at Re 0 →0 (in the free molecular regime) is not significant. The growth of a temperature ratio for the free molecular regime, leads to a considerable growth of the absolute value of CL and to the loss of symmetry of CL (α)
PPT Slide
Lager Image
Pitching moment coefficients MZ for “Clipper” (tw = 0.01)
PPT Slide
Lager Image
Dependencies of drag coefficients CD (α) for “Clipper” on various Reynolds number
in respect to the coordinate origin, and these tendencies were noted to by analysis of the results obtained by the Monte Carlo method. Thus, it can see that the local engineering method gives good results in the calculation of aerodynamic characteristics of hypersonic vehicles in rarefied gas flow. The methods considered do not take into account the influence of the interaction of a boundary layer with the inviscid hypersonic flow at the large numbers of Re 0 .
It has to be pointed out that when the Reynolds number increased, the drag coefficients CD of vehicle diminished, which can be explained by the decrease of normal and tangent stresses ( fig. 7 ). At high Reynolds number Re 0 ≥ 106, characteristics are almost unchanged. The dependency CL (α) is increased at a high Reynolds number, which can be explained by the decrease of normal and tangent stresses ( fig. 8 ). The values of MZ (in respect to the center of mass of a uniform body) are quite sensitive to the variation of Re 0 . MZ changes its sign to less than zero at
PPT Slide
Lager Image
Dependencies of lift coefficients CL (α) for “Clipper” on various Reynolds number
PPT Slide
Lager Image
Dependencies of pitching moment coefficients MZ (α) for “Clipper” on various Reynolds number
Re 0 ~ 10 2 . At Re 0 ~ 10 4 , the value of MZ = - 0.03 at the angle of attack is reached at α ≈ 40 deg ( fig. 9 ). As it follows from these results, the coefficient of moment of tangent MZ is less sensitive to the variations of the temperature ratio, but its value means that this factor should be necessarily taken into account by the analysis of the variations of the body’s orientation, when this body is subject to the action of a flow of strongly rarefied gas.
4. Conclusions
The analysis of different approaches to calculate aerodynamic characteristics of perspective spacecraft vehicles in rarefied gas flow were performed. The engineering method to calculate hypersonic aerodynamics in the transitional regime is described. The calculation results of aerodynamic characteristics for a spacecraft vehicle by the engineering method in transitional regime with various Reynolds numbers are presented. Comparison to Newton’s method and the DSMC method are described. We can reduce calculation errors in vehicle design projects by multiple calculations using this method. The above engineering method in transitional regimes gives a good result for a wide range of bodies. The obtained data can be applied in the future for hypersonic vehicle design.
The authors thank all colleagues Lipatov I.I., Tsipenko V.G., Galkin V.S., Erofeev A.I., Tolstykh A.I., Gorelov S.L., Zharov V.A., Voronich I.V. for participation in useful discussions. The support by RFBR (research project No. 11-07-00300-a) is cordially appreciated by the authors.
Kogan N., M. 1969 Rarefied Gas Dynamic Plenum New York
Belotserkovskii O., M. , Khlopkov Y., I. 2010 Monte Carlo Methods in Mechanics of Fluid and Gas World Scientific Publishing Co. N-Y, London, Singapore, Beijing, Hong Kong
Moss J., N. , Glass C., G. , Hollis B., R. , Van, Norman J., W. 2006 “Low Density Aerothermodynamics of the Inflatable Reentry Vehicle Experiment” AIAA Paper 2006-1189
Horvath T., J. , Berry S., A. , Merski N., R. 2004 “Hypersonic Boundary/Shear Layer Transition for Blunt to Slender Configurations – A NASA Langley Experimental Perspective” Enhancement of NATO Military Flight Vehicle Performance Czech Republic 4-7 October RTO-MPAVT-111.
Ivanov M., S. 2011 “Statistical Simulation of Reentry capsule Aerodynamics in Hypersonic Near-Continuum Flows”, RTO lecture series, Models and Computational Methods for Rarefied Flows EN-AVT-194-18, (AVT-194)
Neumann R., D. 1988 “Missions and requirements. Special Course Aerothermodynamics of Hypersonic Vehicles” Neuilly sur Seine France AGARD Report 761
Hirschel E., H. 2005 Basics of Aerothermodynamics Springer-Verlag Berlin, Heidelberg, New York
Bird G., A. 1994 Molecular Gas Dynamics and the Direct Simulation of Gas Flows Oxford University Press New York
Anderson J., D. 1988 Hypersonic and High-Temperature Gas Dynamics McGraw-Hill New York
Lees L. 1955 “Hypersonic Flow” Processing of 5th International Aeronautical Conference Los Angeles, New York 241 - 275
Bertin J., J. , Cummings R., M. 2006 “Critical Hypersonic Aerothermodynamic Phenomena” Annual Review of Fluid Mechanics 38 129 - 157    DOI : 10.1146/annurev.fluid.38.050304.092041
Shapiro A., H. 1953 The Dynamics and Thermodynamics of Compressible Fluid Wiley New York
Glass C., E. , Holden M., S. , Wieting A., R. 1989 “Effect of Leading Edge Sweep on Shock/Shock Interference at Mach 8” Reno, NV AIAA Paper 89-0271
Reubush D., E. , Omar M., E. 1989 “Pressure and Heat Transfer Investigation of a Modified NASP Baseline Configuration at M=6” Reno, NV AIAA Paper 89-0246
Fujimoto K. , Fujii K. , Tsuboi N. 2003 “CFD Prediction of the Aerodynamic Characteristics of Capsulelike Configurations for the Future SSTO Development” 41st AIAA Aerospace Science Meeting Reno, NV AIAA Paper 03-0912
Galkin V., S. , Erofeev A., I. , Tolstykh A., I. 1977 “Approximate method of calculation of the aerodynamic characteristics of bodies in a hypersonic rarefied gas” Proceedings of TsAGI (in Russian) (1833) 6 - 10
Myint Zay Yar Myo , 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 (5) 551 - 566    DOI : 10.1615/TsAGISciJ.v41.i5.50
Vaganov A., V. , Drozdov S. , Kosykh A., P. , Nersesov G., G. , Chelysheva I., F. , Yumashev V., L. 2009 “Numerical Simulation of Aerodynamics of Winged Reentry Space Vehicle” TsAGI Science Journal 40 (2) 131 - 149    DOI : 10.1615/TsAGISciJ.v40.i2.10