CFD Analysis of Rocket Engine PDF
CFD Analysis of Rocket Engine PDF
J. Steelant
TOS-MPA
ESA - ESTEC
Acknowledgement: R. Schmehl
Noordwijk
The Netherlands
ESTEC
Division of Propulsion and Aerothermodynamics
Oxidizer dome
Injectors
Combustion chamber
Cooling channels
Aestus (Ariane 5)
ESTEC, 9-10th of Multi-Physics in Rocket Engines 2
October 2003 J. Steelant
ESTEC
Division of Propulsion and Aerothermodynamics
CFD-Analysis (1/3)
1. Geometrical aspects
- axisymmetric fuel dome with annular injector slits
- spatial : 20K nodes; temporal: 1 µsec
P_cent
P9D1
CFD-Analysis (2/3)
2. Numerical aspects
- Spatial: 3rd order Smart-scheme
- Time: backward Euler (1st order)
- Solver: Pressure-Correction method
3. Physical modelling:
- Transient laminar and isothermal Navier-Stokes equations
- Vapor transport equations:
MMH vaporization/condensation & presence of
non-condensable gas.
CFD-Analysis (3/3)
Vapor transport equation: f = vapour mass fraction
∂ρf r
+ ∇ • ( ρ uf ) = ∇ • (Γ∇f ) + Re − R f
∂t
è expressions for Re and Rc derived from Rayleigh-Plesset
equation with:
V 2 p sat − p V 2 p − p sat
Re = C e ch ρ v ρ l (1 − f ) Rc = C c ch ρ v ρ l f
σ 3 ρl σ 3 ρl
where Ce, Cs = phase change rate coefficients
σ = surface tension
psat = saturation pressure
Vch = characteristic velocity
ESTEC, 9-10th of Multi-Physics in Rocket Engines 7
October 2003 J. Steelant
ESTEC
Division of Propulsion and Aerothermodynamics
Conclusion