7

Computational Flowfield Analysis of a Planetary Entry Vehicle

Antonio Viviani1 and Giuseppe Pezzella2
1Seconda 2Centro

Universit degli Studi di Napoli, via Roma 29, 81031 Aversa, Italiano Ricerche Aerospaziali, via Maiorise, 81043 Capua Italy

1. Introduction
Computational Fluid Dynamics (CFD) analysis represents a key technology within planetary entry vehicle design. Safe landing of vehicles re-entering from space requires, in fact, an accurate understanding of all physical phenomena that take place in the flowfield past the hypersonic vehicle to assess its aerodynamics and aerothermodynamics performance. CFD allows to significantly reduce the number of in-flight and plasma windtunnel (PWT) experimental test campaigns and to account for real-gas flow features, which are difficult to reproduce in ground-test facilities. Flight measurements collected during reentry have demonstrated that real gas effects strongly influence both aerodynamics and aerothermal loads of hypervelocity vehicles. On the other hand, trajectory calculation for atmospheric re-entry involves determination of vehicle aerodynamics and aerothermodynamics. As a consequence, accurate modeling of flow physics, in particular flow chemistry is fundamental to reliably design re-entry vehicles. In this chapter, we stress this point with an application to a capsule-type crew return vehicle (CRV) for the International Space Station (ISS) support servicing. However, high accuracy in modeling flow and chemistry coupling may produce only a small increase in the numerical results accuracy, despite the high modeling efforts and the increased computational cost. So, one must balance the theoretical and computer time effort needed to use a more general and sophisticated model against the expected accuracy of results. The question then arises as to what extent the number of reactions, coefficients, reaction mechanism, etc. influence the flow. To answer this question, a step-by-step numerical investigation has been carried out to examine the influence of the chemical reactions, its mechanisms and kinetics, and of thermal non-equilibrium on the air flows past the CRV, in the framework of a low Earth orbit (LEO) scenario. Two-dimensional axisymmetric and three-dimensional Navier-Stokes computations are performed, for perfect gas and reacting gas mixture in thermal and chemical non-equilibrium, and for several chemical reaction mechanisms. In particular, simulations are computed with different wallsurface boundary conditions: non-catalytic wall (NCW), partially catalytic wall (PCW), fully catalytic wall (FCW) to underline the effect of the heat shield catalyticity on the vehicle aerodynamic heating. The work confirms that high-temperature transport phenomena markedly influence the vehicle flowfield and, in turn, the vehicle aerodynamics and aerothermodynamics, but it also stresses that, with an acceptable loss of results accuracy, we

www.intechopen.com

124

Numerical Simulations - Examples and Applications in Computational Fluid Dynamics

do not need to use models of such high complexity, and therefore considerable computing time can be saved.

2. Real gas effects and re-entry hypersonic flight

During atmospheric descent, re-entry vehicles encounter several flow regimes and thermochemical phenomena: they fly from free molecular to fully continuum phases and, when in continuum, from laminar to fully turbulent flows. When freestream enthalpy is large enough the flow passing through the bow shock dissociates resulting in a several species reacting mixture flow around the vehicle. The thermal and chemical characteristics of the gas in the shock layer are altered depending on the atomic and molecular structure of the air species (Sarma, 1995). For instance, when flow velocity is low, energy is absorbed only into particles vibration and rotation degrees of freedom (dof). But as velocity sufficiently increases, the thermal energy of the gas becomes comparable with the energy associated with a whole range of gas phase chemical processes, such as the excitation of molecular modes of vibration; the dissociation of oxygen and nitrogen; the formation of other chemical species through recombination reactions; the ionisation of both molecular and atomic species. As a consequence, the flowfield chemical composition around the re-entry vehicle varies spatially and temporally and, because shock layer molecules continuously exchange their energy between the translational and internal dof, the air can result in a thermal-and/or chemical non-equilibrium mixture. Then, the microscopic structure of the mixture species, affecting the ways in which energy may be redistributed, influences the specific heat ratio (), the chemical reaction rates, and the transport properties. These quantities, in turn, affect the dynamics of the flow as well as shock and expansion waves (i.e. pressure, temperature, and velocity distributions), the chemical energy diffused to the surface (i.e. the chemical contribution to the heat flux at the wall), the boundary layer structure (i.e. the heat flux and shear stress). In particular, the flow chemical dissociation results in a large density ratio () across the strong bow shock, which markedly influences the capsules aerodynamics. In fact, influences the shock shape, the stand-off distance, and the wall-surface pressure that, at the stagnation point (e.g., Cpmax), reads:

C p max = C pt 2 =
where , in the hypersonic limit, is:

2 Pt 2 P Pt 2 = 1 2 2 q P M

(1)

= lim

1 1 = 2 + 1

(2)

High temperature effects also modify the hypersonic capsule-vehicle aerodynamics and aerothermodynamics by means of a very abrupt change in the CRV trim angle of attack (trim). This is due to the shift of the sonic line position at the vehicle leeside because of the change in , thus affecting the CRV pitching moment (CMY) (Hassan et al., 1993) and, hence, the capsule static stability that is a critical requirement for a re-entry vehicle, because static instability could lead to catastrophic failure if the thermal shield is not protecting the vehicle anymore. Real gas effects also influence vehicle aeroheating since thermal protection material (TPM) could promote the chemical recombination at wall of flowfield atomic

www.intechopen.com

Computational Flowfield Analysis of a Planetary Entry Vehicle

125

species thus increasing the overall heat flux up to two times or more than the value of a noncatalytic wall (Scott, 1997). For instance, the reactions considered above, taking place only in the gas phase, are termed as homogeneous chemical reactions and differ from the heterogeneous ones that, instead, occur near the vehicle wall involving gas and solid species. They can be catalyzed by the TPM and, being exothermic, contribute to the vehicle aeroheating. Thus, the TPM, promoting or preventing species recombination at wall, depending on its catalyticity, plays an important role in the aerodynamic heating. Neglecting conduction into the heatshield and radiation from the gas, the energy balance at the vehicle surface reads:
N Tv,i T 4 $ r = Tw q = tr + ,i n n w i =1
NS Yi + C Pi dT + h o fi D i n w i =1 To w T

(3)

The first term, on the right-hand side, is the conductive heat-flux, the second one is the vibrational contribution, and the last one is the species diffusion contribution that strongly depends on the catalytic properties of thermal protection system (TPS). Therefore, the heatshield should be a poor catalyst (Anderson, 1973). Of course the entire above mentioned scenario depends on the kind of re-entry (i.e. orbital or superorbital one). For example, flowfield computation involving ionized species, as for superorbital reentries, requires at least 11 chemical species with 20 reactions, whereas for lower velocity reentries, 5 non-ionized species and 17 reactions are sufficient (Sarma, 1995). Therefore, a reliable numerical simulation of re-entry flows can be very challenging, depending on the more or less correct and accurate modelling of the flowfield thermochemical processes. In this framework simulation problems may arise as the coupling of flow and chemistry leads to a stiff problem due to differences in reaction rate characteristic times (Anderson, 1989); dissociation rate coefficients can differ by orders of magnitude and, since reaction rates are very difficult to measure, different values may exist for the same coefficient. As a result, the appropriate set of reactions to be used represents a very relevant choice, especially if one considers that, in general, an increased model complexity does not entail a greater accuracy of numerical results, despite the higher computational cost needed for increased reactions set. Moreover, when one increases the number of chemical reactions, numerical results can be more influenced by the effect of the uncertainty in input data, such as species transport coefficients, relaxation times for thermal and chemical non-equilibrium. So, it could be important to simplify the reaction mechanisms, by reducing as much as possible the number of chemical reactions, without loss in accuracy but greatly reducing computing time.

3. CRV concept and re-entry flight scenario

The re-entry system is an Apollo-like capsule measuring about 5 m in diameter (D), with a nose radius (RN) of 6.05 m; the sidewall angle () of 33 deg and the overall vehicle height of 3.8 m (see Fig. 1). The offset centre of gravity (cg) is located at x/D=0.26 and y/D=-0.0353. This vehicle concept represents a scaled-up version of the ARD capsule, which is a flying test bed successfully experimented by ESA in October 1998 (Walpot, 2001). The reference mission scenario considered for the CRV is the re-entry from the ISS orbit performed by a vehicle weighting about 9 ton, starting from the atmospheric entry interface (hE=120 km) with VE=8 km/s inertial, and E=-2 deg. The re-entry flight scenario is given in

www.intechopen.com

126

Numerical Simulations - Examples and Applications in Computational Fluid Dynamics

both altitude-velocity map and Mach-Reynolds plane in Fig. 2. These re-entry trajectories have been computed by means of the ENTRY (ENtry TRajectrY) code developed at SUN (Viviani, 2006). The blue curve is a ballistic descent trajectory, while the red one refers to a lifting return since the capsule, flying trimmed at =20 deg constant over the critical heating regime, is employing aerodynamic lift to sustain the descent flight path. As shown, the capsule, moving from a very rarefied atmosphere to a denser one, shifts from the free molecular flow (FMF) regime, where Kn10 and individual molecular collisions are important, to the transition one, where 10-3<Kn<10 and slip effects are important, and then to the continuum regime where Kn10-3. For instance, the similarity parameter that governs these different flow regimes is the Knudsen number (Bertin, 1994):

Kn L

ref

M = 1.25 Re L L ref

(4)
ref

where Lref is the characteristic length of the body (e.g., the capsule diameter). Therefore, vehicle aerodynamics must also consider the effect of the Kundsen number. Re-entry trajectories of Fig.2 result in an aerothermal environment that must be accurately predicted for a reliable TPS design (Pezzella et al., 2007). To this end, trajectory freestream flight conditions have been used to perform numerical computations, as reported hereinafter.

3.8 (m)

= 33 (deg) RN= 6.05 (m)

HRSI RCC
5.0 (m)
Fig. 1. CRV Vehicle configuration with quotes and TPS layout with high reusable surface insulator (HIRS) and reinforced carbon-carbon (RCC) tiles

4. Design approaches and numerical tools

The preliminary AErodynamic DataBase (AEDB) of the CRV has been provided according to the space-based design approach (Prabhu, 2004), which dictates the generation of data set as function of a number of independent parameters (i.e. M, Re, , ). On the other hand, the preliminary AeroThermodynamic DataBase (ATDB) has been computed following the trajectory-based design approach, which consists in performing aerothermal computations at a finite number of critical points of the nominal re-entry trajectory (Olynick, 1998).

www.intechopen.com

Computational Flowfield Analysis of a Planetary Entry Vehicle

10 x 10
4

127

lifting trajectory, AoA=20 Ballistic trajectory 9

8 ReD=1e5 7 ReD=5e5 6 Altitude [m] ReD=1e6 5 ReD=5e6 4 ReD=1e7

1 M=4 0 0 1000 M=8 M=12 3000 M=16

M=20

M=24

2000

4000 5000 Velocity [m/s]

6000

7000

8000

9000

a)
50

40 Mach number

Free Molecular

Transition

Continuum

30

M/Re=10 M/Re=0.1 M/Re=0.001 CRV Ballistic CRV Lifting H=87 km H=120 km

20

10

0 1.E-02

1.E+00

1.E+02 1.E+04 Reynolds Number, ReLref

1.E+06

1.E+08

b) Fig. 2. Re-entry flight scenario; a) altitude-velocity map; b) Mach-Reynolds map with constant Knudsen numbers Computational analysis of the CRV flowfield is performed by means of the solver code Fluent together with several user defined functions (UDF), developed by the authors, in order to manage vibrational relaxation, several catalyticity models, radiative equilibrium at the wall and other boundary conditions. Computations have been carried out on multiblock structured grids. A close-up view of both 2-D and 3-D mesh can be seen in Fig. 3. The grid used for 3D calculations consists of 32 blocks with about 900.000 cells. For each computational case, a new grid has been created to properly accommodate for the detached bow shock location. See (Viviani et al., 2008) for further details on the computational grids.

www.intechopen.com

128

Numerical Simulations - Examples and Applications in Computational Fluid Dynamics

Fig. 3. Close-up view of 3-D (on top) and 2-D axi-symmetric computational mesh domains

5. The model and the numerical technique

The mathematical model describing the flowfield around a hypervelocity vehicle deals with balance equations for a multi-species chemically reacting gas mixture supplemented with an appropriate set of chemical reactions (i.e. the reactions mechanism) and with equations modelling species vibrations relaxation (i.e. thermal non-equilibrium). 5.1 Flowfield governing equations The governing equations are made up of equations for mass conservation, total momentum balance, total energy (without the vibrational one) conservation, individual species balance and vibrational energy conservation. The full set of equations for a laminar viscous compressible continuum flow in thermal and chemical non-equilibrium, assuming the air as a mixture of NS perfect gases and Nv vibrating species, written in the integral conservation is (Gnoffo et al., 1989):

f where W = , u, v, w,et , 1 ,....., Ns-1 , ev1 ,....., evNv is the unknown state vector of the conserved quantities, in which is the fluid density, u, v and w are the momentum densities, et is the total internal energy per unit mass, i and evi are, respectively, the density and the vibrational energy of the ith species while evi takes into account for vibrational

f* f f f f f f W d + Finv + Fvis n dS + A* inv + A vis dV = dV t V rV S V

(5)

www.intechopen.com

Computational Flowfield Analysis of a Planetary Entry Vehicle

129

f f energy conservation. F and A* , splitted into an inviscid and a viscous part, are the flux vector and the axis-symmetric terms matrixes, respectively. is equal to 1 for axisf T symmetric flows and 0 for 2D and 3D flows. = [0 ,0 ,0 ,0 ,0 , 1 ,....., Ns 1 , 1 ,......, N ] is the source terms vector. It defines the mass and energy exchange among the species as a result of the chemical reaction rate and the energy transfer due to the internal energy excitation processes (Bertin, 1994). Finally, V is the arbitrary control volume cell, S is its closed f boundary control surface and n is the outward normal unit vector (Anderson, 1989). Eq.(5) can be written in differential form as follows: Continuity:

f f + V = 0 t

( )
f

(6)

Momentum:

f V t

f ( ) +

( V ) ( VV ) + p = 2
ff f f f

o s

(7)

Energy:

( et ) t

f f f f f = + 2 T V + + e p V ( ) t
( Yi ) t

( )

f f $i e $vj V + hi J i hi o i j i
s

(8)

Species:

f f f f $i + VYi + J i =

(9)

Vibrational energy:

evj t

f $i In these equations, V is the velocity vector, Yi is the mass fraction of the ith species and, is the rate of change of i due to chemical reactions, Ji is the diffusive flux of ith species,
which arises due to concentration gradients, Mi and hi are, respectively, the molecular weight and enthalpy of ith species, p is the pressure, is the viscosity and is the thermal conductivity. For each species the perfect gas model applies and the Daltons law holds:

f ) +

( Ve ) = e$
f
vj

(10)

vj

p = pi
i

(11)

where pi is the partial pressure of the ith species of the mixture. As a consequence, the following relation for density reads:

p R0T Yi / M i
i

(12)

where R0 is the universal gas constant. The internal energy of the mixture is defined as:

www.intechopen.com

130

Numerical Simulations - Examples and Applications in Computational Fluid Dynamics

e = (Yi ei )
i

(13)

where ei, the internal energy of the single component gas, is the sum of the energies representing the different degrees of freedom of the molecules. Finally, the enthalpy is:
h = (Yi hi )
i

(14)

Computation of the diffusive fluxes requires knowledge of the transport coefficient. 5.2 Transport properties For pure species from kinetic theory of gases (Anderson, 1989) follows that: Viscosity:

i =
Thermal conductivity:

2.6693 106 M iT

i2 i

(15)

i =
Mass diffusivity:
Dij =

15 i R0 4 Mi
3

4 c pi M i 1 + 3 15 R0

(16)

0.0188 T 2

2 Dij p ij

( Mi + M j ) / MiM j

(17)

Global transport properties of the gas mixture rely on semi-empirical rules such as Wilkes mixing rule for viscosity and thermal conductivity :
a=
i i a i 1 1 1 2 1 Mi 2 a i 2 Mi 4 + 1+ 1 jj 8 Mj aj Mj

a = ,

(18)

where i is the mole fraction of species i and ai (equal to i or i) is obtained by kinetic theory of gases. For the diffusion coefficient of the ith species in the mixture the multicomponent diffusion coefficient is applied:

Di =

(1 i )

(19)

Di , j

with Di,j evaluated by kinetic theory. Finally, vibrational relaxation is modelled using a Landau-Teller formulation, where relaxation times are obtained from Millikan and White, assuming simple harmonic oscillators (Bertin, 1994).

www.intechopen.com

Computational Flowfield Analysis of a Planetary Entry Vehicle

131

5.3 Chemical species and reactions mechanism Within a LEO re-entry scenario (e.g., no flowfield ionization occurs), the gas is approximated as a finite-rate chemistry mixture of N2, O2, N, O, NO species (Gnoffo et al., 1999). The elementary reactions mechanism, governing the species in high-temperature air, deals with three dissociation reactions and two exchange reactions, as reported in Tables 1a, 1b; there M, namely reacting partner or third body, can be any of the five reacting species, thus providing or removing collision energy. Efficiencies of the third body are also reported in Tables 1a, 1b, since they are employed in computations to increase CPU time efficiency. The reactions mechanism results in a system of 17 chemical reactions, with 17 forward and backward reactions rate coefficients. No 1 2 3 4 5 6 Reaction
O2 + M = 2O + M
N2 + M = 2N + M
N2 + N = 3N

Af,r T (m3/kgmole s) (k) 3.60x1015 1.90x1014 4.085x1019 3.90x1017 3.20x106 7.00x1010 T T T T T T

f,r
-1.0 -0.5 -1.5 -1.5 1.0 0.0

Ea f,r (J/kgmole) 4.947x108 9.395x108 9.395x108 6.277x108 1.638x108 3.159x108

Third body efficiency O2=9, N2=2, O=25, N=NO=1 O2=1, N2=2.5, O=N=NO=1 O2=N2=1, O=N=NO=20 -

NO + M = N + O + M NO + O = O2 + N N 2 + O = NO + N

Table 1a. Reaction rate parameters in Eq. (21), Dunn & Kang model (Gnoffo et al., 1989) Af,r (m3/kgmole s) 1.00x1019 3.00x1019 1.10x1014 2.40x106 1.80x1011

No 1 2 3 4 5

O2 + M = 2O + M

Reaction

T
(k) Ta Ta Ta T T

f,r
-1.5 -1.6 0.0 1.0 0.0

Ea f,r (J/kgmole)

Third body efficiency

N2 + M = 2N + M

NO + M = N + O + M

N 2 + O = NO + N

NO + O = O2 + N

4.947x108 O2=N2=NO=0.2, O=N=1 O2=N2=NO=0.233, 9.412x108 O=N=1 O =N2=0.05, 2 6.277x108 O=N=NO=1 1.598x108 3.193x108 -

Table 1b. Reaction rate parameters in Eq. (21), Park model (Park et al., 1993) Of the simpler sets of homogeneous reactions, it is standard to use the following three chemical reactions, known as Zeldovich process (Sarma, 1995):

O 2 + N 2 = 2O + N 2 N 2 + O = NO + N NO + O = O 2 + N

(20)

It considers only Oxygen dissociation, due to collisions with molecular Nitrogen, and two exchange reactions. This model can be explained considering that the gas is so hot that the Oxygen completely dissociates, while the Nitrogen does neither dissociate completely nor as fast as the Oxygen. Finally, exchange reactions are important because they determine the

www.intechopen.com

132

Numerical Simulations - Examples and Applications in Computational Fluid Dynamics

speed of Nitrogen dissociation. Chemical reactions proceed with forward rates, kf,r that $ i ) of the species transport equation, Eq. (9). They are expressed appear in the source terms ( in the Arrhenius form as:

k f ,r = k f ,r (T ) = k f T a Tvb = Af ,rT

f ,r

Ea exp f ,r RT o

(21)

where T and the constants depend on the model kinetics (see Tables 1a, 1b). The Dunn-Kang model uses one temperature to describe all the energy modes (e.g. T = T , a b thermal equilibrium), while the Park model assumes that T can be T, Tv, or T Tv (namely rate controlling temperature) depending on the reaction (see Table 1b). So Parks twotemperature model provides more accurate results because it uses T to describe translational and rotational energy modes and Tv for vibrational and electron-translational modes. 5.4 Boundary conditions Eq. (3) states that the properties of a surface are emissivity () and wall catalyticity (i.e. kwi). As atoms of dissociated flow strike vehicle surface, the catalyticity property of the wall is $ w i 0) in the boundary layer problem to implemented by means of a production term (i.e. solve. Then, steady-state mass atomic conservation at the wall states that the production of ith species, due to the catalytic recombination rate, must be balanced by the rate of diffusion to the surface:
$ a is given by Goulards relationship: The source term

$ a ) w = ( a va ) w (
$ a = kwa ( wYiw

(22)

(23)

where p is the reaction order and kwi is the catalytic reaction rate (Anderson, 1989). The diffusive flux a va is expressed by means of Ficks law and then:

$ a ) = kwa ( Yi ) ( w w

$ w i = 0), the TPS surface is When the TPM does not promote any particular reaction (i.e., called NCW (i.e. kwi=0); on the opposite situation, when the TPM can activate any reactions, the TPS surface is called FCW (i.e. kwi). Between these two limit cases (i.e. 0<kwi<), the vehicle surface is considered as PCW and the heat fluxes to the vehicle greatly differ depending on the value of kwi. Furthermore, when a low conductive TPS protects the vehicle, the radiative equilibrium condition holds at vehicle surface. To account for this condition, during numerical simulations, the wall temperature is calculated by StephanBoltzman law and is implemented by means of a Newton-Raphson approach.

Y = Da i n w

(24)

5.5 Numerical technique The governing equations, together with the proper boundary conditions, are discretized using a cell-centered finite volume formulation on a structured multiblock grid (see Fig. 3). For the single-mesh cell, the discretized equation reads:

www.intechopen.com

Computational Flowfield Analysis of a Planetary Entry Vehicle

f dWc 1 + dt Vc

(F
6 f =1

f
c inv

f + Fc vis

f f f f* n + Ac inv + Ac* vis = H c r

133

(25)

where f is the index of the cell face and c refers to the single cell of the computational domain. The inviscid fluxes at cell interfaces are calculated by using a flux difference splitting (FDS) Riemann solver, since upwind methods are particularly suitable for high speed flows. However, second order accuracy is not automatically reached. For this reason, a second order essentially non oscillatory (ENO) technique for the reconstruction of cell interface values is employed (Anderson, 1995). The viscous fluxes are calculated by central differencing, i.e. computing the gradients of flow variables at cell interfaces by means of Gauss theorem. The method is second order accurate in space. Time integration is performed by employing both an explicit single-stage (Euler forward) algorithm and an explicit five stage Runge-Kutta scheme, coupled with an implicit evaluation of the chemical and vibrational source terms, under the hypothesis of time marching approach to reach the steady solution for the flow.

6. Reliability of numerical study

A reliable flow simulation dictates the validation of the theoretical models describing the high temperature effects in a hypersonic flow by means of PWT and free flight experimental data. To this end, in order to assess the reliability of present results a numerical rebuilding of ELECTRE (Muylaert et al., 1999) and ARD (Walpot, 2001) test campaigns were performed. Experimental and numerical comparisons of results for both the test models were reported. 6.1 ELECTRE test article in HEG wind-tunnel ELECTRE test article (see Fig. 4) consists of a blunt conical surface 0.4 m long, semiaperture cone angle of 4.6 deg, and hemispherical nose radius of 0.035 m. It was tested in flight and in PWT, becoming a standard reference model to study non-equilibrium hypersonic flow past blunt-body configurations (Muylaert et al., 1999). The CFD computational domain, shown in Fig.4, consists of 60x120 cells with a minimum normal wall spacing of 10-5m which a grid sensitivity analysis has shown to be necessary to obtain a sufficient resolution of the flowfield features.
Y (m)

0.2

0.1

0.1

0.2

0.3

0.4

X (m)

Fig. 4. ELECTRE test article geometry and axi-symmetric mesh domain

www.intechopen.com

134

Numerical Simulations - Examples and Applications in Computational Fluid Dynamics

Test conditions are summarized in Tab.2. They correspond to operating conditions of the HEG PWT located at DLR Gottingen at which two different test cases were run considering alternatively the specimen wall as NC, and FC for N and O species. T (K) 790 P (Pa) 430

(kg/m3) 1640x10-6

V (m/s) 5919

M (-) 9.7

Re/m (1/m) 270x103

YO (-) 0.179

YN (-) 1.0x10-6

YNO (-) 3.3x10-2

YO2 (-) 3.6x10-2

Table 2. Reference calibration point conditions of the HEG nozzle (Muylaert et al., 1999) As an early assessment of work several engineering evaluations have been accomplished by means of ENTRY as, for example, pressure coefficient and convective heat transfer distributions on the test article wall. ENTRY quickly evaluates the pressure distribution by using the surface impact method, typical of hypersonics:
C p = C pt 2 sen 2

(26)

where Cpt2 is the stagnation-point pressure coefficient that depends on the flow theory one considers, while is the local slope body angle (Bertin, 1994). For example, Newtonian flow theory states that Cpt2 is equal to 2 while in the case of Modified Newtonian theory (see Eq.(1)) it reads:

P 2 C pt 2 = t 2 1 2 P M

(27)

For the aeroheating, ENTRY evaluates the convective heat transfer around both hemispherical nose and spherically capped cone, according to the Lees theory (Lees, 1956). Lees approach shows that, for a blunt cone with nose radius RN and semiaperture vertex angle c, at any point on the cone surface, the ratio of heat transfer to the stagnation value

$ co reads: q

Nose

$ w ( ) 2 sin cos 2 q = $ co q D ( )
3 $ w (s ') s' q s' B ( c ) + = A ( c ) $ co q RN RN 1 2

(28)

Eq.(28) is valid for s/RN cotc, where s is the curve length measured along the cone surface of the effective sharp-nosed cone and, for high flight Mach number:
A(c ) 3 sin c 2 3 1 B(c ) 16 sin 4 c c 2 D() cot 3 c = c
2

Cone skirt

(29)

1 1 D() 2 sin 4 + (1 co s 4 ) 2 8

www.intechopen.com

Computational Flowfield Analysis of a Planetary Entry Vehicle

135

CFD aerothermodynamic computations have been performed with different wall catalytic boundary conditions for the test article. The computations refer to fully laminar nonequilibrium flow conditions with model temperature fixed to Tw=300 K. The flowfield past the test bed is shown in Fig.5, where the Mach number contour field is plotted, comparing the results for perfect gas (upper side) and real gas model (lower side).

Perfect gas model Real gas model

Mach Number: 0.3 0.8 1.3 1.7 2.2 2.7 3.2 3.7 4.1 4.6 5.1 5.6 6.0 6.5 7.0 7.5 8.0 8.4 8.9 9.4

Fig. 5. Mach number contours field. Comparison between perfect gas (top) and real gas model. Detail on the nose region of the test bed (right) Pressure coefficient and wall heat flux distributions for different wall catalytic conditions are reported respectively in Fig. 6 and Fig. 7, where CFD results were compared with engineering and available experimental results (i.e. shot157, shot159, shot164) (Muylaert et al., 1999). The Cp comparison of numerical, experimental and engineering results shows good agreement in the first part of the test specimen, with differences only at the model end. In particular, Fig. 6 highlights that no differences exist between NC and FC wall boundary conditions, as expected. Moreover, ENTRY results (see Fig. 6-b) compare well with the Modified Newtonian theory, where the Cpt2 has been evaluated using the inflow specific heat and Mach number as provided by CFD analyses. As shown, the pressure on the probe nose decreases with curve length from the stagnation point, becoming constant on the conical skirt (e.g., =c=const in Eq.(26)). Modified Newtonian theory, however, loses in accuracy as highlighted by CFD since the surface pressure depends on many factors such as the interaction of compression and of expansion waves which originate from body curvature, reflection from the bow shock and slip lines due to the rotationality introduced by the curved bow shock wave. In particular, as the cone is very slender, the surface pressure in the expanding flow decreases so slowly that the asymptotic (sharp cone) value is not reached, resulting in an underexpanded flow. The heat flux distribution (see Fig. 7) shows an agreement with the numerical FC solution (the red curve) on the nose of the test article while on the rear part of cone there is a mismatch between experimental data and CFD results, as already seen in the case of Cp. Both these mismatches could be probably caused by flowfield perturbations due to the support arm, located at the end of the test bed. Anyway Fig.7-b shows that CFD results compare well with Lees theory.

www.intechopen.com

136

Numerical Simulations - Examples and Applications in Computational Fluid Dynamics

10

Fluent cat Fluent noncat shot 157 shot 159 shot 164 off axis
0

10 Pressure coefficient, cp

10

-1

10

-2

0.1

0.2

0.3

0.4 0.5 0.6 0.7 Nondimensional length along axis, x/L

0.8

0.9

a)

10

Newtonina theory Modified newtonian theory

Fluent 10 Pressure coefficient, Cp

0

10

-1

10

-2

3 4 Nondimensional curve length, s/Rn

b) Fig. 6. a) Cp comparison between numerical, engineering, and experimental data, b) detail on the nose probe (Muylaert et al., 1999)

www.intechopen.com

Computational Flowfield Analysis of a Planetary Entry Vehicle

137

10

Fluent cat Fluent noncat shot 157 shot 159 shot 164 off axis 10 Heat flux density, W/m2
7

10

10

0.2

0.4 0.6 Nondimensional length along axis, x/L

0.8

a)
10
0

Lees Theory Fluent

Nondimensional heat flux

10

-1

10

-2

0.5

1 1.5 2 Nondimensional curve length, s/Rn

2.5

b) Fig. 7. Heat flux; a) comparison between numerical, engineering, and experimental data; b) detail on the nose probe (Muylaert et al., 1999)

www.intechopen.com

138

Numerical Simulations - Examples and Applications in Computational Fluid Dynamics

6.2 ARD model in S4 wind-tunnel S4 test campaigns provide pressure profiles on the forebody centerline of the atmospheric re-entry demonstrator (ARD) capsule (see Fig.1) at the freestream conditions of Tab. 3 (Walpot, 2001). P0 (bar) T0 (k) ReD M P (Pa) T (k) Twall (k) Xcg/D Ycg/D AoA (deg) 85 1151 967237.3 9.92 211.3 55.7 300 0.26 0.0353 -20 25 1108 319208 9.72 71.17 55.7 300 0.26 0.0353 -20

Table 3. ONERA S4 exit conditions (Walpot, 2001)

2.0 1.8 1.6 Pressure Coefficient 1.4 1.2 1.0 0.8 0.6 0.4 0.2 0.0 -1.2 -0.8 -0.4 0.0 s/Rb 0.4 0.8 1.2
CFD W alpot Exp. data Po=25/85 bar CFD present computation

Fig. 8. Pressure coefficient: comparison between present results and data in (Walpot, 2001) Results of present computations in terms of Cp on forebody centerline are summarized and compared with experimental data in Fig. 8, where CFD data provided by Walpot are also reported as further benchmark (Walpot, 2001). As one can see, experimental and numerical data compare well, thus confirming the reliability of numerical results.

7. CRV flowfield computational analysis

Axi-symmetric and fully three-dimensional simulations were performed at the freestream conditions listed in Table 4, with the far field composed of 79% N2 and 21% O2. In order to appreciate how vehicle flowfield depends on real gas effects, numerical computations have been performed in a step-by-step approach. For instance, starting with the perfect gas model, flowfield computation accounts for the chemistry considering first the flow in equilibrium conditions and then in non-equilibrium ones. In the latter case, the effects of both reaction mechanism and chemical kinetics are considered, for example, by means of Zeldovich reaction mechanism, and Dunn-Kang and Park kinetic models. Finally,

www.intechopen.com

Computational Flowfield Analysis of a Planetary Entry Vehicle

139

thermal non-equilibrium and wall catalyticity are accounted for, providing a complete overview of high temperature effects on CRV aerodynamics and aerothermodynamics. Altitude (Km) 50 57 57 57 57 57 57 Mach (-) 10 19 19 19 19 16 12 Pressure (Pa) 79.78 32.78 32.78 32.78 32.78 32.78 32.78 Temperature (k) 270.65 255.27 255.27 255.27 255.27 255.27 255.27 AoA (deg) 0 10 17.5 21 28 28 28

Table 4. Freestream conditions of CFD computations All the computations considered in this work have been performed by means of an Intel Core Duo E7300 at 2.66 GHz. 7.1 Axi-symmetric computations A general overview of the flowfield past the vehicle is shown in Fig. 9, where contours of pressure and translational temperature in the forebody flowfield are plotted. They refer to a perfect gas computation performed at M=19, H=57 km and =0 deg (e.g., peak heating conditions of the ballistic trajectory of Fig. 2). Figure 10 shows a comparison of the non-dimensional temperature profile between perfect gas (PG), equilibrium gas (EG), chemical non-equilibrium gas with NCW and chemical nonequilibrium gas with FCW, as evaluated along with the stagnation line. As clearly shown, the temperature is large enough to cause the complete Oxygen dissociation in the shock layer while Nitrogen partially dissociates. Moreover, differences can be found in temperature peak, stand-off distance and also in the equilibration trend for translational temperature. In particular, in the case of chemical non-equilibrium computation, the temperature profile on the stagnation line exhibits a sharp discontinuity at the shock wave and a large overshooting value due to the finite rate dissociation of molecules.

Fig. 9. Contours of static pressure (left) and temperature for =0 deg, M=19, and H=57 km

www.intechopen.com

140

Numerical Simulations - Examples and Applications in Computational Fluid Dynamics

body
Fig. 10. Comparison of static temperature for =0 deg, M=19, and H=57 km, along with the stagnation line for different chemical models Figure 10 also shows that the shock layer becomes thinner because high temperature phenomena absorb heat, thus decreasing the effective specific heat ratio . As a result, the gas compressibility changes, which lead to changes in the shock wave shape around the vehicle. In particular, as the green curve highlights, this phenomenon occurs more prominently in the equilibrium flow than in a non-equilibrium one. When we account for the influence of chemical kinetics, we have provided in Fig. 11 the comparison of non-dimensional temperature profiles, along the stagnation line, for the results of Dunn-Kang and Park models.

Fig. 11. Comparison of translational temperature, along the stagnation line, for different chemical models and reaction mechanisms, for =0 deg, M=19, and H=57 km. As shown, chemical kinetics slightly changes both stand-off distance and the peak temperature. Differences are also in the shape of temperature profiles in the shock layer,

www.intechopen.com

Computational Flowfield Analysis of a Planetary Entry Vehicle

141

because for the Dunn-Kang model the flow equilibrates before reaching the boundary layer, differently from the case of Park model. As to the effects of the reaction mechanism, Fig. 11 (right side) also reports the effect of the Zeldovichs model on translational temperature along the stagnation line. As one can see, this reaction mechanism over-estimates both temperature peak and stand-off distance. In particular, the temperature profile of complete reaction mechanism (the blue curve), shows that flowfield in the shock layer tends toward equilibrium faster than Park and Zeldovich results. The effects of vibrational relaxation can be appreciated in Fig. 12, where the comparison between translational temperature distributions, along the stagnation line, is reported for Park kinetics. Results shows only slight differences between the case of thermal equilibrium and nonequilibrium computation. When vibrational equilibrium holds (green curve) there is no incubation time for vibration to relax, as highlighted by the temperature rise of the curve. Therefore, we can conclude that at the peak heating the flowfield around the capsule is characterized by an almost thermal equilibrium conditions (remember that capsule features a large forebody radius, RN=6.05 m). Hence, chemistry is active just behind the shock and energy goes only in chemistry thus yielding lower translational temperature in the shock layer. In fact, even if the exact functional dependence is unknown, it is agreed that a reaction mechanism depending only on T over-predicts the amount of dissociation. Thermal equilibrium conditions are also confirmed by results summarized in Fig. 12 (right side), where the comparison between the translational and vibrational non-dimensional temperatures is reported. As we can see, temperature profiles are almost overlapping except across the shock, with the vibrational temperature of molecular nitrogen that slightly lags behind the others (e.g., the energy transfer takes a certain number of collisions to proceed). In particular, O2, NO and N2 quickly equilibrate.

Fig. 12. Translational and vibrational temperatures along with the stagnation line for =0 deg, M=19, and H=57 km for NCW. Comparison between thermal equilibrium (TEQ) and non-equilibrium (TNEQ) conditions for Park kinetics

www.intechopen.com

142

Numerical Simulations - Examples and Applications in Computational Fluid Dynamics

7.2 Three-dimensional computations In the framework of 3-D CFD computations, several AoAs have been considered (see Tab.4). Figure 13 shows the static temperature contours on capsule symmetry plane and two flowfield cross sections and the static pressure contours on the CRV surface at M=19 and =20 deg, considering the flow as a reacting gas mixture. As shown, the CRV bow shock structure around the descent vehicle can be appreciated as well. Figure 14 reports the flowfield contours of Oxygen (O) and Nitrogen (N) mass fraction (on the right side) for M=19, H=57 km, and =20 deg. As expected, the oxygen is fully dissociated. Concerning capsule static stability, the transition in sonic line location is shown in Fig. 15 to Fig. 17 for different AoA and Mach numbers. Figure 15 reports sonic line location for four cases, involving two different Mach numbers (e.g., 10 and 19) and all the AoA considered in the computations. As shown, the sonic line location and its shape markedly depend on freestream Mach number, altitude and AoA, thus highlighting that capsule attitude conditions are highly influenced during descent. In order to appreciate the effect of finite rate chemistry, the sonic line comparison between PG and RG computations both for M=16 and M=19 is summarized in Fig. 16. In both cases for the PG solution the flowfield around almost all the capsule heat shield is entirely subsonic with consequent high pressure distribution. Therefore, as the sonic line shifts due to the chemical reactions, the flow becomes entirely supersonic and the pressure decreases, thus confirming that pitching moment (CMy) and trim are affected by real gas effects. In particular, differences between the green and blue lines on the capsule afterbody, underline that at M=19 CMy is affected by the chemical kinetics while no influences are expected when CRV is flying at M=16 (see left side of Fig. 16).

Fig. 13. Static temperature on CRV symmetry plane and two cross sections at M=19, H=57 km and =20 deg. Static Pressure contours on capsule forebody

www.intechopen.com

Computational Flowfield Analysis of a Planetary Entry Vehicle

143

Fig. 14. Flowfield contours of Oxygen and Nitrogen (right) mass fraction for M=19, H=57 km, and =20 deg on CRV symmetry plane and two cross sections. Static Pressure contour on capsule forebody

21

Fig. 15. Sonic line location in the capsule pitch plane for different AoA and Mach number

28 deg

M=19 AoA=28 deg

Fig. 16. Sonic line location in the capsule pitch plane for M=16, M=19, H=57 km and =28 deg. Comparison between PG and RG computations

www.intechopen.com

144

Numerical Simulations - Examples and Applications in Computational Fluid Dynamics

Fig. 17. Sonic line location in the capsule pitch plane at =28 deg and H=57 km. Comparison among M =12, 16 and 19 Finally, Fig. 17 shows the sonic line for M=12, 16, and 19 at the same , thus highlighting the effect of the Mach number. Therefore, region of vehicle static instability could be expected during re-entry, depending on the capsules c.g. location.

8. CRV aerodynamic analysis

Aerodynamic features of primary interest for an axi-symmetric vehicle are lift (CL), drag (CD), and pitching moment (CMy) coefficients, which are calculated according to Eq. (30).
Ci = Fi 1 2 S ref v 2 i = L, D
CM j = Mj 1 2 v Lref S ref 2 j =Y

(30)

where Lref=5.0 m (i.e., CRV diameterlongitudinal reference length); Sref=(D2)/4=19.6 m2 (i.e. CRV maximum cross-section area). Pole coordinates for the CMy calculation are (1.3,0.1765,0) m (i.e., vehicle cg). Based on the reentry scenario of Fig. 2 the AEDB has been generated for FMF, transitional regime and continuum flow. No lateral directional analysis has been taken into account in this work. CRV aerodynamic appraisal within FMF and transitional regime entails Direct Simulation Monte Carlo Method (DSMC) computations and a very simple relationship to bridge the transitional flow regime from continuum to FMF one (Pezzella et al., 2009). CRV continuum aerodynamics refers to both engineeringbased analysis and CFD-based analysis. 8.1 Engineering-based aerodynamics of CRV As an early assessment of the CRV continuum aerodynamics several engineering evaluations have been accomplished by ENTRY considering that when the capsule travels at hypersonic velocities the aerodynamic forces are dominated by pressure effects. Viscous and other effects, such as base drag, represent only about ten percent of the total. So the flow may be approximated as inviscid and surface inclination methods (SIM), like Modified Newtonian theory (MN), can be used (Bertin, 1994). To this end let us consider Fig.18 where the capsules outer mold line (OML) is shown. Assuming that the freestream particles impact only on the frontal area of the body and cannot curl around it, combining Eq. (26) and Eq. (1), MN theory suggests that the pressure coefficient (Cp) becomes:

www.intechopen.com

Computational Flowfield Analysis of a Planetary Entry Vehicle

145

CP = (2 ) sen 2

(31)

By integrating Eq. (31) over the whole vehicle surface one is able to evaluate aerodynamic forces acting on the vehicle, both for zero lift and AoA cases.
y

+CMy CY CD

+ u

CP x CL

CX

Fig. 18. Capsule OML used for surface impact method aerodynamics 8.1.1 Zero lift aerodynamics The zero lift drag coefficient of the CRV is given by:

CDo =

where is the capsule frustum angle (i.e. 33 deg).

2 2 (2 ) 2 (1 + cos ) CP r sin cos d = 2 Sref 0

(32)

8.1.2 AoA aerodynamics When the capsule is flying at an AoA the pressure coefficient on vehicles wall facing flow, reads:

n ) = ( 2 )( co s co s + sin sin ) CP = ( 2 ) ( u
2

(33)

So, when < (i.e. the sidewall surface is shadowed) the axial and normal force coefficients are, respectively:

CX =

(2 ) cos 2 (1 + cos 2 ) + 1 sen 2 sen 2

2 2

CY =

(2 ) cos sen sen 2

2

(34)

(35)

Therefore, the capsule lift to drag ratio (E) is:

www.intechopen.com

146

Numerical Simulations - Examples and Applications in Computational Fluid Dynamics

E=

C L C Y cos C X sen = C D C Y sen + C X cos

(36)

By imposing the pitch moment balance around the vehicle heat shield nose we have:

CM Y =

x cp L ref

CY =

(2 ) x
2L ref

cp

cos sen sen 2

(37)

where xcp is the abscissa of the vehicle center-of-pressure (cp) that, as said before, is expected to be influenced by the real gas effects, thus affecting capsule trim angle and in turn its descent flight due to the lower aerodynamic efficiency of vehicle (Park et al., 1992). Then, for very high Mach numbers, say larger than five, SIM states that the aerodynamic characteristic of the capsule does not depend on vehicle velocity, but rather on the and the geometric angle of the body shape . Further expressions to compute lift and drag coefficients are:

CD = CDo + 12(1 CDo ) sen 2

6(6 5C Do ) sen 4

+ 4(6 5CDo ) sen 6

(38)

5 CL = 2(1 CDo ) (3 CDo ) sen 2 sen 2

(39)

8.2 CFD-based aerodynamics of CRV Curves of aerodynamic efficiency and pitching moment coefficient are shown in Fig. 19. Note that, CMY, shows that the pitching moment derivative CM is negative in the AoA range of 150-180 deg, signifying that the CRV is statically stable for this range of (provided that the capsule cg is close to the moment reference point). Other numerical data, provided in (Crowder et al., 1969), are also reported to highlight the accuracy of results. As to the effect of chemical kinetics on CRV aerodynamics, Fig. 20 shows that at M=16 no differences are expected for aerodynamic forces, passing from Dunn-Kang to Park kinetics. On the contrary, at M=19 both CL and CD increase ranging from PG through DK and Park kinetics. Differences of about 2% and 5% are observed concerning the value of CL and CMY, respectively. These conclusions for M=19 are confirmed by the displacement of the abscissa of the vehicle centre-of-pressure (xcp) nondimensionalized with respect to its value for MN (Viviani et al., 2010). The increase in Cp at the stagnation region and its decrease over the remaining region lead to a forward (toward nose) shift of centre-of-pressure, or equivalently, positive (nose-up) pitching moment, as the flow decreases. Figure 20 also shows that at M=16, the xcp in the case of PG solution decreases of about 10% with respect to the MN estimation and there are no differences between the values provided by both the chemical models. At M=19, instead, even if the value of xcp in the case of PG solution decreases again of about 10% with respect to the MN estimation, the two chemical models provide values that differ from each other of about 2%. Note that the latter difference may be dangerous if neglected when designing vehicle thermal shield layout. Analyses of the effect of chemical kinetics on capsule aerodynamics at =28 deg versus Mach number, can be found in (Viviani et al., 2010).

www.intechopen.com

Computational Flowfield Analysis of a Planetary Entry Vehicle

147

L/D 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0.0 130 140 150 AoA (deg) 160 170 180
CFD RG M =12 CFD RG (Tab.1a) M =16 CFD RG (Tab.1b) M =16 CFD RG (Tab.1a) M =19 CFD RG (Tab.1b) M =19 CFD PG M =19 CFD PG M =19 (AoA=21 ) CFD RG (Tab.1a) M =19 (AoA=17.5 ) CFD PG M =19 (AoA=10 ) CFD RG (Tab.1a) M =19 (AoA=10 ) CFD PG M =10 CFD RG (Tab.1a) M =10 Crow de r & M oote

0,08

CM Y(pole @cg)

0,06

0,04

0,02
C rowder Moote C F DR GM= 12 C F DR G(Tab.1a) M= 16 C F DR G(Tab.1b) M= 16 C F DR G(Tab.1a) M= 19 C F DR G(Tab.1b) M= 19 C F DP GM= 19 C F DP GM= 19 (A oA = 21 ) C F DR G(Tab.1a) M= 19 (A oA = 17.5 ) C F DP GM= 19 (A oA = 10 ) C F DR G(Tab.1a) M= 19 (A oA = 10 ) C F DP GM= 10 C F DR G (Tab.1a) M= 10

0,00

-0,02

-0,04

-0,06 130 140 150

A oA(deg)

160

170

180

Fig. 19. Capsule L/D and CMy versus . Comparison among present and literature data (Crowder et al., 1969) Finally, Fig. 21 shows, in a step-by-step approach, the effects of reaction mechanism on CRV aerodynamics when the capsule is flying at =28 deg, M=19 and H=57 km. As shown, the results for Zeldovich model compare globally well with those of the complete reaction mechanism, i.e. 17 reactions. Such a result is very interesting considering that Zeldovich results differ from those of full reaction mechanism of 1% only for all the aerodynamic coefficients. In this case, in fact, the CPU time efficiency increases of about 40 %. Therefore, the gain in terms of solution speed up is very high if compared to the loss in

www.intechopen.com

148

Numerical Simulations - Examples and Applications in Computational Fluid Dynamics

accuracy of 1% only. The same consideration can be made in the case of aerodynamic results obtained in the case of O2, N2 and NO dissociation only, even if here the CMy differs of about 5 %. When results refer to O2 and N2 dissociation only, CPU time efficiency increases of about 45% but the CMy accuracy evaluation is in the range of 10%. Finally, when aerodynamic coefficients account for only O2 dissociation the speed up efficiency reaches about 65%, whereas when flow dissociation is neglected (e.g., PG solution) the gain is nearly 70%. In this case, however, CMy differs of about 15% and 25%, respectively. As a result, Fig. 21 highlights that CPU speed up depends on the accuracy expected in vehicle pitching moment assessment. Therefore, this confirms that CMy is one of the most critical parameters within vehicle aerodynamic design.

AoA=28 deg 1.02 1.00 0.98 0.96 C-/C-(MN) 0.94 0.92 0.90 0.88 0.86 0.84 0.82 MN PG DK
AoA=28 deg 1.02 1.00 0.98 Xcp/Xcp(MN) 0.96 0.94 0.92 0.90 0.88 0.86 0.84 MN PG DK Park
M=16 M=19

Lift @ M=16 Lift @ M=19 Drag @ M=16 Drag @ M=19

Park

Fig. 20. Effect of chemical reactions on CRV aerodynamics at = 28 deg for M=16 and M=19. Comparison among MN, PG, Dunn-Kang and Park computations

www.intechopen.com

Computational Flowfield Analysis of a Planetary Entry Vehicle

1.010 1.005 1.000 0.995 0.990 CL, CD 0.985 0.80 0.980 0.975 0.970 0.965 0.960 PG .O2 diss .O2+N2 diss .O2+N2+NO Zel'dovich diss 17 react.s 0.75 Lift Coeff. Drag Coeff. Pitching mom. Coeff. 0.70 0.65 0.60 CMy 0.85 1.05 1.00 0.95 0.90

149

Fig. 21. Effects of reaction mechanism on CRV aerodynamics at =28 deg, M=19 and H=57 km 8.3 Free molecular flow and transitional aerodynamics of CRV Drag coefficients, evaluated by means of DSMC simulations, have been reported and compared in Fig. 22, with those assessed through a quick engineering estimation provided by a bridging relationship between FMF results and continuum ones. For instance, a very simple relationship to bridge the transitional flow regime from continuum regime to FMF reads:
C DTransitional = C DContinuum + ( C DFM C DContinuum ) C D

(40)

2.5 2.3 2.1

D

Drag Coefficient, C

1.9 1.7 1.5 1.3 1.1 0.9 0.7 0.5 1.0E-04 1.0E-03
85 km 105 km 115 km

125 km

200 km

Bridging formula @ AoA=0 deg DSMC @ AoA=0 deg

1.0E-02 1.0E-01 1.0E+00 Knudsen number, KnLref

1.0E+01

1.0E+02

Fig. 22. Drag coefficients vs Knudsen number for =0 deg (Pezzella et al., 2009) where the normalized coefficient C i uses Knudsen number as the independent parameter:

www.intechopen.com

150

Numerical Simulations - Examples and Applications in Computational Fluid Dynamics

CD =

10-3<Kn<10 and C DContimuum and C DFM are the aerodynamic drag coefficient in continuum and FM flow regimes, respectively. The effect of rarefaction on the aerodynamic drag is clearly shown considering that CD at =0 deg increases of about 24% passing from 85 km to 125 km, whereas the drag at H=200 km is 25% higher than the one at 85 km.

C D C DContinuum = F(Kn ) = sen 2 ( 3 + Log 10 Kn ) C DFM C DContinuum 8

(41)

9. CRV aerothermodynamic features

The CRV aeroheating analysis reveals two critical regions on the vehicle surface (Viviani et al., 2007). They are the stagnation point (capsule flying at =0 deg) and the vehicle side corner (capsule flying at 0 deg). Several Navier-Stokes computations have been performed assuming the chemically reacting gas model, considering alternatively the heat shield surface as NC, PC and FC wall. The computations refer to fully laminar flow conditions with capsule surface temperature fixed at 300 K or in radiative equilibrium conditions. Note that, due to relatively low entry velocities (e.g., no radiation heat flux applies) only convective heat flux are taken into account. Moreover no heat shield ablation and recession were assumed for simplicity. In order to illustrate the high temperature real gas effects in air, Fig. 23 shows the comparison of heat flux on the forebody centerline, in the cases of PG, EG and chemical nonequilibrium gas; for the latter case, the results for NCW and FCW are also reported. All these heat flux profiles are evaluated for cold wall conditions (e.g., Tw=300 K).

Fig. 23. Real gas effect on CRV forebody heat flux at = 0 deg, M=19 and H=57 km. Tw=300 K

www.intechopen.com

Computational Flowfield Analysis of a Planetary Entry Vehicle

151

As shown, the wall heat flux is higher in the case of chemical equilibrium condition and it is closer to the heat flux for PG simulation. This can be explained considering that as the wall is cold, the chemical equilibrium leads to the recombination of dissociated atoms. When the boundary layer (BL) is in chemical non-equilibrium, CFD simulations confirm that the higher heat flux is attained for FCW. As shown, this value is lower than that for EG, but it is very large compared with the case of NCW, as expected. Therefore, we can conclude that if the gas in the BL is in equilibrium (e.g., fast recombination) then a surface catalyst will not have any effect on the formation of molecules. In this case, in fact, atoms recombine and liberate their energy of dissociation to the gas in the BL. This added heat tends to increase the heat flux to the surface via thermal conduction, thus indicating that the recombination of atoms is more important than in the case of a FCW. For this reason, we regard the equilibrium condition as the reference condition in much of the state of the art TPS design activities. The equilibrium condition hypothesis, however, may lead to an excessively conservative assessment of vehicle aeroheating. Therefore, realistic heat flux assessment over vehicle surface demands non-equilibrium flowfield computations with a full reaction mechanism, as wall catalyticity plays a significant role when assessing vehicle aeroheating.

Fig. 24. Effects of chemical kinetics (left side) and of vibrational relaxation (right side) on capsule forebody total heat flux at = 0 deg , M=19 and H=57 km for Tw=300 K As far as effect of chemical kinetics is concerned, Fig. 24 shows that in the case of a FCW a $ c is due to the energy released by the recombination of atoms so that the large part of q influence of chemical kinetics is small. On the contrary the differences in the heat flux profile, as evaluated for NCW, underline that the role of reaction rate may be very important in order to assess vehicle aerodynamic heating. Therefore, we can conclude that the chemical model is negligible for a catalytic wall and is significant in the non-catalytic case. Of course the latter conclusion is expected depending on freestream conditions. For instance, at high altitude the shock layer has not dissociated since pressure and density are low. Therefore, the number of particle collisions is low so that the chemical reactions are not significantly activated (e.g., the flow is nearly frozen): there is not much energy involved in dissociation. Most of the energy is in translational modes. As flying altitude decreases the

www.intechopen.com

152

Numerical Simulations - Examples and Applications in Computational Fluid Dynamics

density suddenly increases and the chemical reactions are activated. At very low altitude, the shock layer is so dense that the BL is close to equilibrium (e.g., atoms recombine in the BL before they have a chance to strike the wall). Hence, the BL is heated and there is more potential for heat transfer by conduction (i.e., reactions rate comeback to have a little effect). Figure 24 also displays the effect of vibrational relaxation on the wall heat flux. As shown, the heat flux profile in the case of FCW conditions depends on vibrational relaxation more than in the case of NCW (Viviani et al., 2007). In order to illustrate the magnitude of the effects of catalytic activity on TPS, Fig. 25 shows the comparison of the heat flux along the forebody centreline at wall radiative equilibrium conditions, between the cases of FCW, PCW and NCW, in the case of lifting re-entry (3D computations). As one can see, the overheating caused by the catalytic action is potentially very large compared with the case of NCW. In particular, the largest difference occurs at the sphere-cone junction (corner fillet) where large changes in the flow gradients along the surface occur. Therefore, the corner radius is the dominant geometric feature for the convective heating (instead of heat shield radius of curvature), and it is confirmed that significant reduction in convective heat flux occurs if the thermal shield is built with a noncatalytic TPM.
700

600

FCW PCW NCW

500 Heat flux density, KW/m2

400

300

200

100

0 -1.5

-1

-0.5 0 0.5 Nondimensional forebody length, s/Rb

1.5

Fig. 25. Heat flux on the forebody OML. Nondimensional surface distance measured from centerline. Lifting trajectory. Vehicle trimmed at =20 deg

10. Conclusions
Real gas effects on the design of a Crew Return Vehicle for the International Space Station have been taken into account. CFD analysis has been performed to simulate the re-entry of an Apollo-shaped capsule in order to determine the effect of thermal and chemical nonequilibrium on flowfield past the vehicle. Numerical results confirm that high temperature air conditions strongly affects both the aerodynamics and the aerothermodynamics of capsule. In particular, a comparison is made between CFD computations for perfect gas and reacting gas mixture, in order to bring into evidence the effects of finite rate chemistry, the

www.intechopen.com

Computational Flowfield Analysis of a Planetary Entry Vehicle

153

role of reaction mechanism with the related chemical kinetics and, finally, the influence of wall catalyticity. Model simplification is also investigated by considering a restricted set of predominant chemical reactions such as those of the Zeldovich model. Results comparisons for CRV aerodynamics confirm the strong effect of the real gas behavior on the capsule static stability. In particular, aerodynamic coefficients obtained with the Zeldovich model are within 1% of that of a solution with complete reaction mechanism, but the simulation speed up efficiency is of about 40 %. Therefore, the gain in terms of solution speed up is very high if compared with the loss in accuracy of 1% only. Anyway, in the framework of vehicle aerodynamics, the CPU speed up depends on the accuracy expected in vehicle pitching moment assessment thus confirming that the CMy is one of the most critical vehicle design parameters. On the contrary as far as vehicle aerothermodynamic is concerned a reliable heat flux assessment over vehicle surface demands flowfield computations with a full reaction mechanism as wall catalyticity plays a significant role in the assessment of vehicle aeroheating. Finally, the work underlines that the exact prediction of the heat transfer and chemical environment is crucial for the design of the vehicle TPS. In fact, the possibility of reducing the heat loads on the surface of space vehicles has been highlighted.

11. References
Anderson, J. D. (1989). Hypersonic and High Temperature Gas Dynamics, McGraw-Hill Book Company, ISBN 0-07-001671-2, New York. Anderson, J. D. (1995). Computational Fluid Dynamics: The Basics with Applications, McGrawHill Book Company, ISBN 0-07-001685-2, New York. Anderson, L. A. (1973). Effects of Surface Catalytic Activity on Stagnation Heat Transfer Rates. AIAA Journal, Vol. 11, No.5, page numbers (649-656), doi: 10.2514/3.6806. Bertin, J. J. (1994). Hypersonic Aerothermodynamics, AIAA Education series, ISBN 1-56347-0365, Washington. Crowder R. S.; Moote J. D., (1969). Apollo Entry Aerodynamics. Journal of Spacecraft and Rockets, Vol. 6, No. 3, page numbers (302-307), doi: 10.2514/3.29589. Gnoffo, P. A.; Weilmuenster K. J.; Hamilton, H. H.; Olynick, D. R.; Venkatapathy, E. (1999). Computational Aerothermodynamic Design Issues for Hypersonic Vehicles. Journal of Spacecraft and Rockets, Vo. 36, No. 1, page numbers (21-43). doi: 10.2514/2.3430. Gnoffo, P. A.; Gupta, R. N.; Shinn, J. (1989). Conservation Equations and Physical Models for Hypersonic Air Flows in Thermal and Chemical Non-equilibrium. NASA TP 2867. Hassan, B.; Candler, G. V.; Olynick, D. (1993). Thermo-Chemical Non-equilibrium Effects on the Aerothermodynamics of Aerobraking Vehicles. Journal of Spacecraft and Rockets, Vol. 30, No. 6, page numbers (647-655). doi: 10.2514/3.26369. Lees, L., (1956). Laminar Heat Transfer over Blunt-Nosed Bodies at Hypersonic Flight Speeds, Jet Propulsion, Vol. 26, page numbers (259-269). Muylaert, J.; Walpot, L.; Wennemann, D. (1999). A Review of European Code-Validation studies in High-Enthalpy Flow. Philosophical Transactions. The Royal Society. London. A 357, page numbers (2249-2278). Olynick, D. (1998). Trajectory-Based Thermal Protection System Sizing for an X-33 Winged Vehicle Concept. Journal of Spacecraft and Rockets, Vol. 35, No. 3, page numbers (249257). doi: 10.2514/2.3338.

www.intechopen.com

154

Numerical Simulations - Examples and Applications in Computational Fluid Dynamics

Park, C.; Lee Seung-Ho. (1993) Validation of Multi-Temperature Nozzle Flow Code Noznt. Proceedings of the 28th AIAA Thermophysics Conferences, Orlando (FL), USA. paper AIAA-93-2862. Park, C.; Yoont S. (1992). Calculation of Real-Gas Effects on Blunt body Trim Angles. AIAA Journal, Vol. 30, No. 4, page numbers (1146-1146). doi: 10.2514/3.48931. Pezzella, G.; Battista, F.; Schettino, A.; Marini, M.; De Matteis, P. (2007). Hypersonic Aerothermal Environment Preliminary Definition of the CIRA FTB-X Re-entry Vehicle. Proceedings of the West-East High Speed Flow Field Conferences; Moscow, RUSSIA. Nov. (19-22). Pezzella, G.; Votta, R. (2009). Finite Rate Chemistry Effects on the High Altitude Aerodynamics of an Apollo-Shaped Reentry Capsule, 16th AIAA/DLR/DGLR International Space Planes and Hypersonic Systems and Technologies Conference, Bremen, Germany, Oct. 19-22. AIAA-2009-7306. Prabhu, D. K. (2004). System Design Constraints-Trajectory Aerothermal Environments. RTO AVT/VKI Lecture Series in Critical Technologies for Hypersonic Vehicle Development. Rhode St. Gense, Belgium. May (10-14). Sarma, G. S. R. (1995). Aerothermochemistry for Hypersonic Technology. VKI Lecture series 1995-04. Rhode St. Gense, Belgium. Scott, C. D. (1987). The Effects of Thermochemistry, Non-equilibrium, and Surface Catalysis in the Design of Hypersonic Vehicles. 1th Joint Europe-US Short Course on Hypersonic. GAMNI-SMAI. Paris, France Dec. (7-11). Viviani, A.; Pezzella, G.; Cinquegrana, D. (2006). Aerothermodynamic Analysis of an Apollo-like Re-entry Vehicle. Proceedings of the 14th AIAA/AHI Space Planes and Hypersonic Systems and Technologies Conferences. 6-9 November, Canberra (AU); paper AIAA-2006-8082. Viviani, A.; Pezzella, G.; Borrelli, S. (2008). Effect of Finite Rate Chemical Models on the Aerothermodynamics of Re-entry Capsules. Proceedings of the 15th AIAA Space Planes and Hypersonic Systems and Technologies Conferences. 28 April-1 May, Dayton (OH); paper AIAA-2008-2668. Viviani, A.; Pezzella, G. (2007). Catalytic Effects on Non-Equilibrium Aerothermodynamics of a Re-entry Vehicle. Proceedings of the 45th AIAA Aerospace Sciences Meeting and Exhibit. 8-11 January, Reno, NE (USA); paper AIAA-2007-1211. Viviani, A., Pezzella, G. (2007). Influence of Surface Catalyticity on Re-entry Aerothermodynamics and Heat Shield. Proceedings of the 39th AIAA Thermophysics Conference. 25-28 June, Miami, FL (USA); paper AIAA-2007-4047. Viviani, A., Pezzella G. (2010). Computational Flowfield Analysis over a Blunt-Body Reentry Vehicle. Journal of Spacecraft and Rockets. Vol.47, No. 2, page numbers (258-270). ISSN 0022-4650. doi: 0.2514/1.40876. Walpot, L. (2001). Numerical Analysis of the ARD Capsule in S4 Wind Tunnel. ESA SP.487.

www.intechopen.com

Numerical Simulations - Examples and Applications in Computational Fluid Dynamics Edited by Prof. Lutz Angermann

ISBN 978-953-307-153-4 Hard cover, 440 pages Publisher InTech

Published online 30, November, 2010

Published in print edition November, 2010 This book will interest researchers, scientists, engineers and graduate students in many disciplines, who make use of mathematical modeling and computer simulation. Although it represents only a small sample of the research activity on numerical simulations, the book will certainly serve as a valuable tool for researchers interested in getting involved in this multidisciplinary eld. It will be useful to encourage further experimental and theoretical researches in the above mentioned areas of numerical simulation.

How to reference

In order to correctly reference this scholarly work, feel free to copy and paste the following: Antonio Viviani and Giuseppe Pezzella (2010). Computational Flowfield Analysis of a Planetary Entry Vehicle, Numerical Simulations - Examples and Applications in Computational Fluid Dynamics, Prof. Lutz Angermann (Ed.), ISBN: 978-953-307-153-4, InTech, Available from: http://www.intechopen.com/books/numericalsimulations-examples-and-applications-in-computational-fluid-dynamics/computational-flowfield-analysis-of-aplanetary-entry-vehicle

InTech Europe

University Campus STeP Ri Slavka Krautzeka 83/A 51000 Rijeka, Croatia Phone: +385 (51) 770 447 Fax: +385 (51) 686 166 www.intechopen.com

InTech China

Unit 405, Office Block, Hotel Equatorial Shanghai No.65, Yan An Road (West), Shanghai, 200040, China Phone: +86-21-62489820 Fax: +86-21-62489821