Phase-one development of a linear detonation tube test platform — a calculated stepping stone toward full RDE implementation — designed to experimentally validate a physics-based detonation model developed in-house. The program covers propulsion design, structural analysis, instrumentation, test stand construction, safety system architecture, and operational procedures from PDR through hot-fire. The numerical analysis phase is complete and published in AIAA proceedings.
AIAA Publication
Numerical Analysis and Precursor Detonation Tube Design for Continuously Rotating Detonation Engine Development
Sheelam, J.; Sharma, K.; Mooney, K. C.; Coulter, B. J.
AIAA Region II Student Conference 2026 · Embry-Riddle Aeronautical University, Prescott, AZ
Presents a coupled HLLC–Cantera–Lamé numerical framework that captures detonation-scale wave propagation within 2–3% of CJ predictions and confirms structural feasibility of the proposed hardware geometry — providing quantitative justification for fabrication and staged hot-fire testing.
View on AIAA ARC ↗ · DOI: 10.2514/6.2026-108393
Key Published Results
4,358 m/s
Detonation wave speed — within 2–3% of CJ prediction
33.8 MPa
Peak simulated pressure at t = 0.163 ms
4,923 K
Peak wave temperature
59.8 MPa
Max von Mises stress — 25% of yield limit (SS 304)
0 events
Structural violations — no yield, no UTS exceedance
<14 % Δ
ANSYS CFX vs. MATLAB solver agreement
The coupled Euler–finite-rate chemistry framework captured detonation-scale wave propagation within 2–3% of classical CJ velocity predictions. ANSYS CFX 3D simulations corroborated the reduced-order MATLAB solver's pressure-rise trends to within 14% after ignition kernel repositioning. Structural screening confirmed 17-4PH steel geometry supports the detonation environment with significant margin — providing quantitative justification for hardware fabrication and staged experimental testing.
Program Timeline
Physics-Based Numerical Model Complete
MATLAB 1D HLLC solver with Cantera finite-rate chemistry (GRI-Mech 3.0). CJ conditions, detonation wave propagation, compressible flow, transient thermal and structural loading. Wave speed within 2–3% of CJ theory.
ANSYS CFX Multidimensional Corroboration Complete
3D k-ω turbulence modeling for confinement behavior, shear-layer development, and DDT-relevant flow features. Pressure evolution consistent with MATLAB solver to within 14%.
Structural Feasibility — Lamé Analysis Complete
Transient pressure mapped to inner-wall von Mises stress via thick-cylinder Lamé relations. Max σ_vm = 59.8 MPa against SS 304 Sy = 240 MPa. No structural violations. Published in IOP proceedings.
AIAA Paper — Published Published
First and lead author. AIAA Region II Student Conference 2026. DOI: 10.2514/6.2026-108393. Full numerical framework documented and peer-reviewed.
DAQ & Instrumentation Design Complete
NI PXI at 10 MHz. Synchronized multi-channel sensor architecture. EMI mitigation across ignition and sensor lines. LabVIEW live acquisition + MATLAB post-processing pipeline.
Safety Architecture & SOP Pipeline Complete
PDR, CoDR, and safety board reviews completed. Full experimental SOP authored. Remote ignition with pre-ignition purge. Burst-diaphragm overpressure mitigation. EMI-shielded sensor/ignition lines.
Test Stand Fabrication Complete
Full hardware assembly. Regenerative cooling architecture, high-pressure ignition system, suppressor and wave-capture system for safe post-detonation depressurization without vacuum chamber infrastructure.
Hot-Fire Campaign Pending
C₂H₄/O₂ (1:3, φ = 1.2). Pre-ignition purge, remote initiation, 10 MHz synchronized data capture. Model validation against live pressure trace. Directly extends the published numerical framework.
System Specifications
| Oxidizer / Fuel | Oxygen / Ethylene (φ = 1.2, 1:3 ratio) |
| Initial Pressure | 1,114,575 Pa (~11 atm) · enables high-pressure detonation environment |
| Peak Pressure (DAF) | 115–117 MPa structural design target |
| Simulated Peak Pressure | 33.775 MPa at t = 0.163 ms, x = 1.75 ft |
| Detonation Wave Speed | 4,357.7 m/s (numerical) vs. 4,300–4,500 m/s CJ band — 2–3% deviation |
| Primary Material | 17-4 PH Stainless Steel, H1100 Heat Treatment |
| DAQ System | National Instruments PXI · 10 MHz acquisition |
| Software Stack | LabVIEW (acquisition) → MATLAB (post-processing & solver) |
| Computational Tools | MATLAB, Python, Cantera (GRI-Mech 3.0), ANSYS CFX, SOLIDWORKS, Fusion 360 |
| CFD Solver Type | HLLC finite-volume · operator-split chemistry · k-ω turbulence (CFX) · thick-cylinder Lamé structural |
| Program Status | Pre-hot fire — test stand complete · numerical framework published |
Engineering Approach
- Led test stand design as primary decision-maker — owned all hardware decisions across mechanical, electrical, and propulsion subsystems. Team of four engineers with supporting roles.
- Developed the custom MATLAB 1D HLLC solver from scratch: CJ/ZND theory, Euler equations, operator-split Cantera chemistry, shock/detonation front tracking, and Lamé structural screening — all integrated into a single coupled framework. Published.
- ANSYS CFX used for 3D multidimensional corroboration: k-ω turbulence, shear-layer development, and confinement behavior not resolvable in 1D. Identified and resolved secondary wave interference via ignition kernel repositioning.
- Designed regenerative cooling architecture, high-pressure ignition system, and suppressor/wave-capture system for post-detonation safe depressurization — no vacuum chamber infrastructure.
- Full operational pipeline: PDR, CoDR, safety board reviews, formal experimental SOP authored for every stage of the test campaign.
Media — Test Stand & System Architecture