Scramjet Internal Flow
Scramjet Internal Flow: Theoretical Foundations
Overview
Professor, a scramjet is an engine that combusts fuel while keeping the air supersonic without slowing it down, right? Why do we do that?
In a ramjet, air at M>5 is decelerated to subsonic speeds, but during this process, the temperature reaches over 4000K, causing the air to dissociate and lose energy. A scramjet (Supersonic Combustion Ramjet) avoids this total temperature rise by combusting while the air remains supersonic.
Mixing and combusting fuel in a supersonic flow must be extremely difficult, right?
Yes. In a flow moving at M=2-3, fuel must be injected, mixed, ignited, and combustion completed within milliseconds. The residence time is on the order of just 1 ms.
Governing Equations
What equations govern the flow in a scramjet?
The Navier-Stokes equations including chemical reactions (RANS or LES), plus transport equations for each chemical species. For hydrogen fuel, a 9-species (H₂, O₂, H₂O, OH, H, O, HO₂, H₂O₂, N₂) chemical reaction model is standard.
In one dimension, Rayleigh flow (duct flow with heating) is fundamental, and the enthalpy relationship between inlet and outlet is:
The combustion efficiency $\eta_{comb}$ is typically 0.7-0.9 and is key to engine performance.
So combustion efficiency is the key. The 1D model doesn't seem that complex...
The 1D model is used for conceptual design, but actual scramjet flowpaths are complex and three-dimensional. Oblique shock trains at the inlet, shock-vortex interactions around injectors, thermal choking in the combustor, recombination reactions in the nozzle... all need to be predicted by CFD.
Combustion Stability and Flameholding
How do you hold a flame in a supersonic flow?
There are three main methods.
- Cavity Flameholder: Creates a recess (cavity) on the wall to form a recirculation zone that holds hot gases.
- Strut Injection: Inserts a strut (thin plate) into the center of the flowpath and injects fuel into its wake.
- Oblique Shock-Induced Combustion: Uses a shock wave to compress and heat the fuel-air mixture for ignition.
The cavity method, in particular, has been adopted and flight-proven in HIFiRE and the X-51A Waverider. A cavity with a length/depth ratio (L/D) greater than 5 is called an open cavity, and less than 5 is a closed cavity, each with different recirculation patterns.
The X-51A actually flew, right?
In 2013, it successfully achieved 240 seconds of supersonic combustion flight at M=5.1. It used JP-7 hydrocarbon fuel and was designed to utilize endothermic decomposition (breaking down into lighter hydrocarbons via endothermic reactions for simultaneous cooling).
The Historical Dilemma of the Scramjet as a "Dream Engine"
The scramjet concept has been researched since the 1950s, so why isn't it practical yet? One fundamental reason is the chicken-and-egg problem: "To fly it, you first need another propulsion system to accelerate it to supersonic speeds." A scramjet itself only functions above Mach 4-5, requiring a combined system that uses rockets or turbojets for takeoff before switching over. The X-43A used a rocket booster to accelerate before firing its scramjet for just 10 seconds. The governing equations look elegant, but practical application is a continuous series of gritty engineering trade-offs—that's also what makes scramjets fascinating.
Computational Methods for Scramjet Internal Flow
Selection of Turbulent Combustion Models
What kind of turbulent combustion models are used for scramjet combustion CFD?
A challenge specific to supersonic combustion is that the turbulent time scale and chemical reaction time scale are close (Damköhler number near 1). This means laminar flame models that ignore turbulence effects cannot be used.
| Model | Features | Application |
|---|---|---|
| Finite-Rate/Eddy-Dissipation (FR/ED) | Controlled by the slower of reaction rate and turbulent mixing rate. | RANS rough estimation |
| EDC (Eddy Dissipation Concept) | Fine structure reactor model. Compatible with detailed chemistry. | RANS detailed analysis |
| Flamelet/Progress Variable (FPV) | References pre-computed flamelet tables. | Recommended for LES |
| Transported PDF Method | Solves transport equations for probability density functions of chemical species. | High accuracy but high computational cost |
Which is more appropriate, LES or RANS?
For the design stage, RANS+EDC is practical, but for detailed prediction of combustion efficiency or flame structure, LES+FPV or LES+Finite-Rate Chemistry is needed. For LES of mixing-controlled combustion, balancing the detail of the chemical reaction mechanism with computational cost is crucial; the 9-species 19-reaction mechanism for hydrogen (Jachimowski mechanism) is standard.
Numerical Analysis of Supersonic Mixing
How is fuel-air mixing calculated?
Transverse injection into a supersonic flow is a typical problem setup. A bow shock forms where the fuel jet intersects the main flow, and a pair of counter-rotating vortices (CVP) is generated in the jet wake. This vortex structure promotes mixing.
This momentum flux ratio $J$ is the parameter governing jet penetration depth. $J=1-5$ is typical; larger $J$ means the jet penetrates deeper into the main flow.
Reproducing this interaction in CFD must require quite high resolution.
A mesh size of 1/20th the jet diameter or smaller is needed. RANS (SST k-omega) tends to underestimate CVP strength, so DES or higher is desirable for mixing prediction.
Chemical Reaction Mechanism Reduction
Detailed chemical reaction mechanisms are computationally expensive, right?
The 9-species 19-reaction mechanism for hydrogen-air is still manageable, but hydrocarbon fuels (JP-7, ethylene, etc.) require detailed mechanisms containing hundreds of species and thousands of reactions. Directly incorporating this into CFD is unrealistic, so the following reduction techniques are used.
- Skeletal reduction: Removes low-importance species/reactions (reduces ~200 reactions to ~30).
- QSSA (Quasi-Steady State Approximation): Steady-state approximation for short-lived radicals.
- ISAT (In-Situ Adaptive Tabulation): Stores and reuses calculation results of chemical source terms in a table.
- FGM (Flamelet Generated Manifold): Compresses to a 2D table of mixture fraction and progress variable.
ISAT is implemented in Fluent, right?
Yes. Fluent's ISAT feature can accelerate detailed chemical mechanism calculations by 10-100 times. The first timestep is slow, but speed increases dramatically as the table accumulates.
Boundary Conditions and Inlet Conditions
How are scramjet inlet conditions set?
The combustor inlet is actually the state after compression by the inlet shock waves. Typically, this is M=2-3, static temperature 800-1500 K, static pressure 50-200 kPa. This is given as the inlet boundary condition using total temperature, total pressure, and Mach number. Fuel injection is set at the injection location using a mass-flow-inlet. The outlet is a supersonic outflow condition.
Burning Out in Under 1 Millisecond—The Insane Design Conditions of Scramjets
In a scramjet, fuel (hydrogen) has a residence time inside the engine of only about 0.5 to 1 millisecond. Within that time, fuel and air must be mixed, ignited, and combustion completed. For comparison, gasoline engine combustion takes about 10 milliseconds, and jet engines take several milliseconds. To solve this combustion in CFD, numerical methods for "chemically reacting flows" that handle both chemical reaction time scales and flow time scales simultaneously are needed, and a single timestep setting can cause the calculation to blow up. In practice, a staged approach of "first checking the flow field without combustion (cold flow), then adding reactions" is standard procedure.