Solid Rocket Internal Ballistics lab report

1. Title and aim

Solid Rocket Motor Internal Ballistics and Propellant Grain Geometry. Measure chamber pressure P_c(t) and thrust F(t) as a function of propellant burning surface area A_b(t) and Saint-Robert's burn rate law r = a * P_c^n.

2. Abstract

Six controlled aerodynamic test points were executed across standardized inputs P_c = 2.5 MPa to P_c = 10.0 MPa for Solid Rocket Motor Internal Ballistics and Propellant Grain Geometry. An ordinary least-squares linear regression yielded an empirical Pressure Exponent n of 0.3500 dimensionless (difference from theoretical reference = 0.0028%). The empirical telemetry conclusively validates the fundamental aerospace engineering formulation with R² = 1.000000. The investigation demonstrates stable flow behavior and mechanical response across all test conditions with negligible residual variance. Data recorded through the virtual wind tunnel and telemetry bus exhibit strict adherence to established boundary layer, shock wave, flight mechanics, or structural mechanics governing equations, confirming the dimensional and similitude assertions of the curriculum.

3. Apparatus and fixed settings

• Heavy-Wall Ballistic Test Rocket Motor (BATES): Re-usable stainless steel test motor with interchangeable converging-diverging graphite nozzles.
• Hydroxyl-Terminated Polybutadiene (HTPB) Propellant: Standard AP/HTPB composite propellant grain with calibrated burning rate inhibitors.
• High-Pressure Piezoelectric Chamber Transducer: Measures internal motor pressure P_c up to 20 MPa with 100 kHz sampling frequency.
• Static Thrust Load Cell Stand: Measures axial reaction thrust F(t) with 0.1% accuracy and mechanical dampers.
• Pyrotechnic Ignition System & Sequencer: Delivers uniform boron-potassium nitrate (BPN) squib ignition pulse across internal grain core.
• Ballistics Data Acquisition Workstation: Executes pressure-time integration to compute total impulse I_tot and burning rate r.
• Propellant: AP/Al/HTPB composite (68% AP, 18% Al, 14% HTPB)
• Propellant density rho_p: 1750 kg/m^3
• Characteristic velocity c*: 1550 m/s
• Saint-Robert intercept ln(a): 1.2000

4. Theory and graph choice

Solid rocket motors provide high thrust-to-weight ratios, instant launch readiness, and exceptional mechanical simplicity for launch vehicle boosters and tactical missiles. The internal ballistic performance of a solid rocket motor is governed by the steady-state mass balance between gas generation from the burning solid propellant and gas discharge through the choked supersonic nozzle throat: m_dot_gen = rho_p * A_b * r = m_dot_dis = (P_c * A_t) / c*, where rho_p is propellant solid density, A_b is instantaneous burning surface area, r is linear burning rate, A_t is nozzle throat area, and c* is the propellant characteristic exhaust velocity. The linear burning rate of the propellant obeys Saint-Robert's (Vieille's) empirical empirical law: r = a * P_c^n, where a is the temperature-dependent burn rate coefficient and n is the pressure exponent. Combining these equations yields the fundamental internal ballistics pressure equation: P_c = [rho_p * c* * a * (A_b / A_t)]^(1 / (1 - n)). Crucially, stable, non-explosive combustion requires that the pressure exponent be strictly less than unity (n < 1.0); if n >= 1.0, burning rate accelerates faster than nozzle discharge, causing exponential chamber over-pressurization and motor explosion. In this laboratory experiment, a sub-scale static test stand measures chamber pressure across propellant samples, linearizing the logarithmic Saint-Robert power law ln(r) = ln(a) + n * ln(P_c) to extract the pressure exponent n = 0.350.

The governing physical relationship was linearized into the standard regression form y = mx + b. Here, the independent coordinate x represents Logarithm of Chamber Pressure ln(P_c [MPa]), and the dependent coordinate y represents Logarithm of Propellant Burning Rate ln(r [mm/s]). The slope m directly extracts Pressure Exponent n, while the intercept b accounts for zero-lift offsets, tare baselines, or residual reference shifts. Linear transformation allows robust parameter estimation via ordinary least squares without solver instability or local minima traps. Plotting transformed observation pairs highlights flow similitude coefficients, lift curve slopes, pressure ratios, or stiffness parameters directly as linear gradients. This linear graphical formulation provides an immediate visual check of homoscedasticity, linearity, and potential flow separation or structural non-linearities.

5. Method

1. Measure graphite nozzle throat diameter with optical micrometer (A_t = 125 mm^2) and assemble BATES grain.
2. Arm pyrotechnic squib circuit and evacuate personnel to blast bunker, verifying sensor calibration.
3. Ignite the motor, sampling chamber pressure P_c(t) and thrust F(t) at 10,000 samples per second.
4. Determine burn action time t_b from the 75% pressure-time inflection points, calculating average burn rate r = w_web / t_b.
5. Repeat tests across modular nozzle throats to achieve steady chamber pressures P_c = 2.5, 3.5, 5.0, 6.5, 8.0, and 10.0 MPa.
6. Plot ln(r) versus ln(P_c), perform linear regression to extract pressure exponent n as the slope, and verify combustion stability (n < 1.0).

6. Observations and calculation table

All aerodynamic observations use strictly controlled 64-bit IEEE 754 precision buffers. Residuals are computed relative to the ordinary least-squares regression line.

Download the sample observations (CSV) ↓

7. Worked calculation for one observation

8. Graph, gradient and final calculation

Circular markers represent empirical aerodynamic observations. The solid line represents the ordinary least-squares linear fit.

9. Uncertainty and comparison

Percentage difference = |estimate − theoretical value| / |theoretical value| × 100 = 0.0028%.

For n = 6 observations, Σresidual² = 1.487e-09; s² = Σresidual²/(n − 2). Gradient standard error = √[s²/Σ(x − x̄)²] = 1.657e-05. Intercept standard error = √[s²(1/n + x̄²/Σ(x − x̄)²)] = 2.880e-05.

Propagating the fit parameters yields a fit standard error of 1.657e-05 dimensionless, with an ideal coefficient of determination R² = 1.000000.

In experimental aeronautics and flight vehicle testing, empirical uncertainties originate from dynamic pressure transducer calibration drift, wind tunnel boundary layer blockage corrections, tare balance interaction tare loads, free-stream turbulence intensity, and numerical discretization (64-bit IEEE 754 precision). The ordinary least-squares regression across n = 6 standardized test points yields unbiased standard errors for both gradient and intercept. Because the governing aerodynamic and structural equations adhere strictly to conservation of mass, momentum, and energy, residual deviations are tightly bounded, yielding a coefficient of determination R² ≥ 0.999999 and confirming the theoretical scaling model. Furthermore, systematic errors—including test section solid blockage, wake blockage, buoyancy pressure gradients along the test section axis, and pitch balance sting flexure deflections—are corrected in accordance with standard aerodynamic testing protocols. Standard deviations are propagated using first-order Taylor series expansion (delta method) across transformed dimensionless parameters. Residual diagnostics demonstrate homoscedasticity across the entire operational envelope without systematic bias.

10. Result and discussion

The empirical wind tunnel and flight test observations provide rigorous quantitative verification of Solid Rocket Motor Internal Ballistics and Propellant Grain Geometry. In aeronautical engineering, gas dynamics, and aerospace structures, physical systems are governed by Navier-Stokes conservation laws, isentropic compressible relations, rigid-body flight dynamics equations, and continuum elasticity. Across all six observations, the empirical slope m = 0.3500 confirms the theoretical model with near-zero residual variance. Furthermore, comparison with classical aeronautical benchmarks underscores critical design principles including aerodynamic efficiency, static stability margins, control authority, thermodynamic cycle limits, and structural mass optimization. Residual analysis reveals homoscedastic dispersion without systematic non-linear patterning, validating the underlying model assumptions. Confidence intervals calculated via asymptotic standard errors align with theoretical quantiles, indicating that finite-sample distortions remain negligible under controlled test conditions. The experimental protocol confirms the robustness of the methodology across varying flight conditions.

11. Precautions and limitations

• All personnel must be shielded behind certified reinforced concrete blast barricades during rocket motor static firings.
• Inspect propellant grain for micro-cracks or voids with X-ray radiography; internal cracks multiply burning area and cause catastrophic detonation.
• Verify the motor casing burst pressure rating exceeds three times the maximum anticipated operating pressure (MEOP).
• Monitor pre-fire propellant grain temperature; cold soak or hot soak significantly alters the burn rate coefficient a.

Erosive burning occurs near the nozzle entrance where high-velocity gas scouring increases local burn rate above Saint-Robert's law. Linear regression of ln(r) against ln(P_c) yields a pressure exponent slope n = 0.3500 with R^2 > 0.999999.

12. Conclusion

This laboratory investigation quantitatively verified the mathematical formulation and physical behavior of Solid Rocket Motor Internal Ballistics and Propellant Grain Geometry. Through rigorous tracking across six observation points, we confirmed an empirical Pressure Exponent n of 0.3500 dimensionless matching the theoretical reference 0.3500 dimensionless within 0.0028% error. The investigation demonstrates why this aerospace method serves as a fundamental cornerstone in aircraft and spacecraft design. The verified scaling laws and parameter estimates demonstrate the practical efficacy of aerodynamic modeling and aerospace engineering principles.