Scramjet Internal Flow
Theory and Physics
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.
Physical Meaning of Each Term
- Temporal Term $\partial(\rho\phi)/\partial t$: Think of the moment you turn on a faucet. At first, the water comes out spluttering and unstable, but after a while, the flow becomes steady, right? This term describes that "period of change." The pulsing of blood flow from a heartbeat, or the flow fluctuation each time an engine valve opens and closes—all are unsteady phenomena. So what is steady-state analysis? It's looking only at "after sufficient time has passed and the flow has settled down"—meaning setting this term to zero. This drastically reduces computational cost, so solving steady-state first is a basic CFD strategy.
- Convection Term $\nabla \cdot (\rho \mathbf{u} \phi)$: If you drop a leaf into a river, what happens? It gets carried downstream by the current. This is "convection"—the effect where fluid motion transports things. Warm air from a heater reaching the far side of a room is also convection, where the air "carrier" transports heat. Here's the interesting part—this term contains "velocity × velocity," making it nonlinear. That is, as the flow speeds up, this term rapidly strengthens, making control difficult. This is a root cause of turbulence. A common misconception: "Convection and conduction are similar" → They're completely different! Convection is carried by flow, conduction is transmitted by molecules. There's an order-of-magnitude difference in efficiency.
- Diffusion Term $\nabla \cdot (\Gamma \nabla \phi)$: Have you ever put milk in coffee and left it? Even without stirring, it naturally mixes after a while. That's molecular diffusion. Next question—honey or water, which flows more easily? Obviously water, right? Honey has high viscosity ($\mu$), so it flows poorly. Higher viscosity strengthens the diffusion term, making the fluid move "sluggishly." In low Reynolds number flows (slow, viscous), diffusion dominates. Conversely, in high Re flows, convection overwhelms and diffusion plays a supporting role.
- Pressure Term $-\nabla p$: When you push a syringe plunger, liquid shoots out forcefully from the needle, right? Why? Because the plunger side is high pressure, the needle tip is low pressure—this pressure difference provides the force that pushes the fluid. Dam discharge works on the same principle. On a weather map, where isobars are tightly packed? That's right, strong winds blow. "Flow is generated where there is a pressure difference"—this is the physical meaning of the pressure term in the Navier-Stokes equations. A point of confusion here: "Pressure" in CFD is often gauge pressure, not absolute pressure. If results go haywire when switching to compressible analysis, mixing up absolute/gauge pressure might be the cause.
- Source Term $S_\phi$: Heated air rises—why? Because it becomes lighter (less dense) than its surroundings and is pushed up by buoyancy. This buoyancy is added to the equation as a source term. Other examples: chemical reaction heat from a gas stove flame, Lorentz force acting on molten metal in a factory's electromagnetic pump... These are all actions that "inject energy or force into the fluid from the outside," expressed by source terms. What happens if you forget a source term? In a natural convection analysis, forgetting buoyancy means the fluid doesn't move at all—a physically impossible result where warm air from a heater doesn't rise in a winter room.
Assumptions and Applicability Limits
- Continuum Assumption: Valid for Knudsen number Kn < 0.01 (mean free path ≪ characteristic length).
- Newtonian Fluid Assumption: Shear stress and strain rate have a linear relationship (non-Newtonian fluids require viscosity models).
- Incompressibility Assumption (for Ma < 0.3): Treat density as constant. For Mach numbers above 0.3, compressibility effects must be considered.
- Boussinesq Approximation (Natural Convection): Density variation is considered only in the buoyancy term; constant density is used in other terms.
- Non-applicable Cases: Rarefied gases (Kn > 0.1), supersonic/hypersonic flows (requires shock capturing), free surface flows (requires VOF/Level Set, etc.).
Dimensional Analysis and Unit Systems
| Variable | SI Unit | Notes / Conversion Memo |
|---|---|---|
| Velocity $u$ | m/s | When converting from volumetric flow rate for inlet conditions, pay attention to cross-sectional area units. |
| Pressure $p$ | Pa | Distinguish between gauge and absolute pressure. Use absolute pressure for compressible analysis. |
| Density $\rho$ | kg/m³ | Air: ~1.225 kg/m³ @20°C, Water: ~998 kg/m³ @20°C |
| Viscosity Coefficient $\mu$ | Pa·s | Be careful not to confuse with kinematic viscosity $\nu = \mu/\rho$ [m²/s]. |
| Reynolds Number $Re$ | Dimensionless | $Re = \rho u L / \mu$. A criterion for Laminar/turbulent transition. |
| CFL Number | Dimensionless | $CFL = u \Delta t / \Delta x$. Directly related to time-step stability. |
Numerical Methods and Implementation
Turbulent Combustion Model Selection
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.
Upwind Differencing (Upwind)
1st-order Upwind: Large numerical diffusion but stable. 2nd-order Upwind: Improved accuracy but risk of oscillations. Essential for high Reynolds number flows.
Central Differencing
2nd-order accurate, but numerical oscillations occur for Pe > 2. Suitable for low Reynolds number diffusion-dominated flows.
TVD Schemes (MUSCL, QUICK, etc.)
Maintain high accuracy while suppressing numerical oscillations via limiter functions. Effective for capturing shocks and steep gradients.
Finite Volume Method vs Finite Element Method
FVM: Naturally satisfies conservation laws. Mainstream in CFD. FEM: Advantageous for complex shapes and multiphysics. Mesh-free methods like SPH are also developing.
CFL Condition (Courant Number)
Explicit methods: CFL ≤ 1 is the stability condition. Implicit methods: Stable even for CFL > 1, but affects accuracy and iteration count. LES: CFL ≈ 1 recommended. Physical meaning: Information should not travel more than one cell per timestep.
Residual Monitoring
Convergence is judged when residuals for the Continuity Equation, momentum, and energy drop by 3-4 orders of magnitude. The mass conservation residual is particularly important.
Relaxation Factors
Pressure: 0.2-0.3, Velocity: 0.5-0.7 are typical initial values. Reduce the factor if diverging. Increase after convergence to accelerate.
Internal Iterations for Unsteady Calculations
Iterate within each timestep until a steady solution converges. Internal iteration count: 5-20 iterations is a guideline. Residuals should...
Related Topics
なった
詳しく
報告