REPORT DOCUMENTATION PAGE

Similar documents
Thermo-Kinetic Model of Burning for Polymeric Materials

REPORT DOCUMENTATION PAGE

A Realistic Thermal Mitigation Strategy for Solid Rocket Motors

REGENERATION OF SPENT ADSORBENTS USING ADVANCED OXIDATION (PREPRINT)

Report Documentation Page

REPORT DOCUMENTATION PAGE

REPORT DOCUMENTATION PAGE

REPORT DOCUMENTATION PAGE. Theoretical Study on Nano-Catalyst Burn Rate. Yoshiyuki Kawazoe (Tohoku Univ) N/A AOARD UNIT APO AP

Erik L. Swanberg 1 and Steven G. Hoffert 2. Veridian Systems 1, Autometric 2. Sponsored by Defense Threat Reduction Agency

REPORT DOCUMENTATION PAGE

System Reliability Simulation and Optimization by Component Reliability Allocation

Report Documentation Page

Effects of Constrictions on Blast-Hazard Areas

Analysis Comparison between CFD and FEA of an Idealized Concept V- Hull Floor Configuration in Two Dimensions. Dr. Bijan Khatib-Shahidi & Rob E.

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING ORGANIZATION REPORT NUMBER

AIR FORCE RESEARCH LABORATORY Directed Energy Directorate 3550 Aberdeen Ave SE AIR FORCE MATERIEL COMMAND KIRTLAND AIR FORCE BASE, NM

AN EMPIRICAL SHAPED CHARGE JET BREAKUP MODEL

Quantitation and Ratio Determination of Uranium Isotopes in Water and Soil Using Inductively Coupled Plasma Mass Spectrometry (ICP-MS)

P. Kestener and A. Arneodo. Laboratoire de Physique Ecole Normale Supérieure de Lyon 46, allée d Italie Lyon cedex 07, FRANCE

Report Documentation Page

Use of Wijsman's Theorem for the Ratio of Maximal Invariant Densities in Signal Detection Applications

Design for Lifecycle Cost using Time-Dependent Reliability

INFRARED SPECTRAL MEASUREMENTS OF SHUTTLE ENGINE FIRINGS

INFRARED SPECTROSCOPY OF HYDROGEN CYANIDE IN SOLID PARAHYDROGEN (BRIEFING CHARTS)

Attribution Concepts for Sub-meter Resolution Ground Physics Models

REPORT DOCUMENTATION PAGE

Crowd Behavior Modeling in COMBAT XXI

Comparative Analysis of Flood Routing Methods

Diagonal Representation of Certain Matrices

On Applying Point-Interval Logic to Criminal Forensics

Metrology Experiment for Engineering Students: Platinum Resistance Temperature Detector

Interim Research Performance Report (Monthly) 4. TITLE AND SUBTITLE 5a. CONTRACT NUMBER

SW06 Shallow Water Acoustics Experiment Data Analysis

PREDICTION OF HIGH TEMPERATURE GREASE LIFE USING A DECOMPOSITION KINETIC MODEL

CRS Report for Congress

STUDY OF DETECTION LIMITS AND QUANTITATION ACCURACY USING 300 MHZ NMR

Predicting Tropical Cyclone Formation and Structure Change

Parametric Models of NIR Transmission and Reflectivity Spectra for Dyed Fabrics

Coastal Engineering Technical Note

THE EULER FUNCTION OF FIBONACCI AND LUCAS NUMBERS AND FACTORIALS

Dynamics of Droplet-Droplet and Droplet-Film Collision. C. K. Law Princeton University

Optimizing Robotic Team Performance with Probabilistic Model Checking

SMA Bending. Cellular Shape Memory Structures: Experiments & Modeling N. Triantafyllidis (UM), J. Shaw (UM), D. Grummon (MSU)

Two-Dimensional Simulation of Truckee River Hydrodynamics

Broadband matched-field source localization in the East China Sea*

REPORT DOCUMENTATION PAGE

Development and Application of Acoustic Metamaterials with Locally Resonant Microstructures

A report (dated September 20, 2011) on. scientific research carried out under Grant: FA

USMC Enlisted Endstrength Model

Volume 6 Water Surface Profiles

REPORT DOCUMENTATION PAGE

High Resolution Surface Characterization from Marine Radar Measurements

USER S GUIDE. ESTCP Project ER

Real-Time Environmental Information Network and Analysis System (REINAS)

DIRECTIONAL WAVE SPECTRA USING NORMAL SPREADING FUNCTION

Estimation of Vertical Distributions of Water Vapor and Aerosols from Spaceborne Observations of Scattered Sunlight

uniform distribution theory

Ocean Acoustics Turbulence Study

Reverse Ion Acceleration by Laser-Matter Interaction

Expressions for the Total Yaw Angle

aas .C_ ( a u 00 E.cc CD ULL(h cucz 00)1 D7jw S C' cvcc cnl

Sensitivity of West Florida Shelf Simulations to Initial and Boundary Conditions Provided by HYCOM Data-Assimilative Ocean Hindcasts

Fleet Maintenance Simulation With Insufficient Data

LAGRANGIAN MEASUREMENTS OF EDDY CHARACTERISTICS IN THE CALIFORNIA CURRENT

IMPROVED SYNTHESIS AND REACTION CHEMISTRY OF FN 3

Predicting Deformation and Strain Behavior in Circuit Board Bend Testing

REPORT DOCUMENTATION PAGE

Acquisition Review Quarterly Spring 2003

REPORT DOCUMENTATION PAGE

Predictive Model for Archaeological Resources. Marine Corps Base Quantico, Virginia John Haynes Jesse Bellavance

Swash Zone Dynamics: Modeling and Data Analysis

REPORT DOCUMENTATION PAGE

Closed-form and Numerical Reverberation and Propagation: Inclusion of Convergence Effects

FLIGHT TEST VALIDATION OF AN ANALYTICAL METHOD FOR PREDICTING TRAILING CONE SYSTEM DRAG FORCE AND DROOP ANGLE REAGAN K. WOOLF

Effects of Closing Blocks on Hazard Areas

2008 Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies SPECTRAL ANALYSIS OF RADIOXENON

REPORT DOCUMENTATION PAGE

FRACTAL CONCEPTS AND THE ANALYSIS OF ATMOSPHERIC PROCESSES

PIPS 3.0. Pamela G. Posey NRL Code 7322 Stennis Space Center, MS Phone: Fax:

Direct Numerical Simulation of Aeolian Tones

Sums of the Thue Morse sequence over arithmetic progressions

Mass Transport by Second Mode Internal Solitary Waves

SWELL RATIO TESTER REPORT FY 2003

Flocculation, Optics and Turbulence in the Community Sediment Transport Model System: Application of OASIS Results

Molecular Characterization and Proposed Taxonomic Placement of the Biosimulant 'BG'

Z-scan Measurement of Upconversion in Er:YAG

An Invariance Property of the Generalized Likelihood Ratio Test

High-Fidelity Computational Simulation of Nonlinear Fluid- Structure Interaction Problems

Experimental and Theoretical Studies of Ice-Albedo Feedback Processes in the Arctic Basin

Modulation Instability of Spatially-Incoherent Light Beams and Pattern Formation in Incoherent Wave Systems

Determining the Stratification of Exchange Flows in Sea Straits

Sea Ice Model for Marginal Ice Zone

Modeling the Impact of Extreme Events on Margin Sedimentation

Estimation of Vertical Distributions of Water Vapor from Spaceborne Observations of Scattered Sunlight

NUMERICAL SOLUTIONS FOR OPTIMAL CONTROL PROBLEMS UNDER SPDE CONSTRAINTS

Improvements in Modeling Radiant Emission from the Interaction Between Spacecraft Emanations and the Residual Atmosphere in LEO

Wavelets and Affine Distributions A Time-Frequency Perspective

Mine Burial Studies with a Large Oscillating Water-Sediment Tunnel (LOWST)

Improved Parameterizations Of Nonlinear Four Wave Interactions For Application In Operational Wave Prediction Models

Super-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models

Transcription:

REPORT DOCUMENTATION PAGE Form Approved OMB No. 0704-0188 Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, gathering and maintaining the data needed, and completing and reviewing this collection of information. Send comments regarding this burden estimate or any other aspect of this collection of information, including suggestions for reducing this burden to Department of Defense, Washington Headquarters Services, Directorate for Information Operations and Reports (0704-0188), 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302. Respondents should be aware that notwithstanding any other provision of law, no person shall be subject to any penalty for failing to comply with a collection of information if it does not display a currently valid OMB control number. PLEASE DO NOT RETURN YOUR FORM TO THE ABOVE ADDRESS. 1. REPORT DATE (DD-MM-YYYY) 2. REPORT TYPE 3. DATES COVERED (From - To) 01-09-2010 Technical Paper & Briefing Charts 4. TITLE AND SUBTITLE 5a. CONTRACT NUMBER Mitigation of Fuel Fire Threat to Large Rocket Motors by Venting 5b. GRANT NUMBER 5c. PROGRAM ELEMENT NUMBER 6. AUTHOR(S) Kenneth J. Graham (Aerojet) 5d. PROJECT NUMBER R 5e. TASK NUMBER 5f. WORK UNIT NUMBER ARMY088S 7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING ORGANIZATION REPORT NUMBER Air Force Research Laboratory (AFMC) AFRL/RZSB AFRL-RZ-ED-TP-2010-368 4 Draco Drive Edwards AFB CA 93524-7160 9. SPONSORING / MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSOR/MONITOR S ACRONYM(S) Air Force Research Laboratory (AFMC) AFRL/RZS 11. SPONSOR/MONITOR S 5 Pollux Drive NUMBER(S) Edwards AFB CA 93524-7048 AFRL-RZ-ED-TP-2010-368 12. DISTRIBUTION / AVAILABILITY STATEMENT Approved for public release; distribution unlimited (PA #10444). 13. SUPPLEMENTARY NOTES For presentation at the 2010 Insensitive Munitions & Energetic Materials Technology Symposium, Munich, Germany, 11-14 October 2010. 14. ABSTRACT Venting of a container such as a rocket motor or a warhead case is a well-recognized method to potentially reduce the violent response of the system to a fuel fire threat. There have been many proposed rocket motor or warhead venting systems. The thermally-initiated venting system (TIVS) on the AMRAAM rocket motor has been shown to reduce violent response, by cutting the case with a linear shaped-charge. Graham has demonstrated the ARCAPS system in which a small insert of secondary propellant having a lower temperature than the main propellant grain reacts to perforate the rocket motor case, reducing the system response in both fast and slow cookoff. In studying the response of a 120mm mortar in fast cookoff, a manufactured vent was filled with ionomer plastic that melted at a particular temperature leading to a mild reaction. There are many other designs that include stress-risers, thermite plugs or inserts, slotted overwrapped designs and so on. 15. SUBJECT TERMS 16. SECURITY CLASSIFICATION OF: 17. LIMITATION OF ABSTRACT a. REPORT Unclassified b. ABSTRACT Unclassified c. THIS PAGE Unclassified 18. NUMBER OF PAGES SAR 33 19a. NAME OF RESPONSIBLE PERSON Dr. George C. Harting 19b. TELEPHONE NUMBER (include area code) N/A Standard Form 298 (Rev. 8-98) Prescribed by ANSI Std. 239.18

Mitigation of Fuel Fire Threat to Large Rocket Motors by Venting Kenneth J. Graham, Aerojet, Culpeper, VA Introduction: Venting of a container such as a rocket motor or a warhead case is a wellrecognized method to potentially reduce the violent response of the system to a fuel fire threat. There have been many proposed rocket motor or warhead venting systems. The thermally-initiated venting system (TIVS) on the AMRAAM rocket motor has been shown to reduce violent response, by cutting the case with a linear shaped-charge. Graham has demonstrated the ARCAPS system in which a small insert of secondary propellant having a lower temperature than the main propellant grain reacts to perforate the rocket motor case, reducing the system response in both fast and slow cookoff. In studying the response of a 120mm mortar in fast cookoff, a manufactured vent was filled with ionomer plastic that melted at a particular temperature leading to a mild reaction. There are many other designs that include stress-risers, thermite plugs or inserts, slotted overwrapped designs and so on. Problem: The question that is generally overlooked is what is the critical vent size to prevent overpressurization and how is it determined. The problem we are trying to solve is how to protect a large rocket motor while in the transportation mode typically truck and specialized trailer. This scenario provides the highest probability of a large rocket motor experiencing a fuel fire whether from a rupture and ignition of the truck s own fuel tanks in a crash, or running into some source of flammable fuel from another truck, a car, or even a service station gasoline pump. The basic solution to mitigation by venting is to understand the competition between pressure rise rate and pressure decay rate. Pressure Rise Rate: Kinney and Sewell [1] determined, from interior ballistics, the rate of pressure rise from combustion of an energetic material. The basic form is given in Equation (1) below: dp/dt = RT B /V * dn/dt (1) where dn/dt is the time rate of change of the number of moles of product gases. This equation may be replaced with one in which the variables are more easily measurable. Thus, dp/dt = RT B /V * /M * /(A-BT 0 ) * S B P (2) where: R = molar gas constant = 8.314 x 10-5 bar - m 3 /mol - K V = volume, m 3 T B = flame temperature, K M = formula mass product gas, kg/mol = density of explosive, kg/m 3 Distribution A: Approved for public release; distribution unlimited.

Re cipro ca l Burning Ra te, s/ m m T 0 = bulk temperature of explosive, K,A,B = energetic material constants (see below) S B = burn surface area, m 3 P = absolute pressure, bars The term [ /(A BT 0 )] represents the variation in burning rate with bulk explosive temperature. In experiments with Composition B explosive, it was found that the burning rate at ambient was 0.2 mm/s [2] and from thermal analysis via DSC, violent decomposition occurs about 513K [3]. This is considered the critical temperature as the burning rate is assumed to be infinite at this point. Utilizing the methodology of Andreev [2, 4, 5], we plot the reciprocal of burning rate against bulk explosive temperature. For Composition B, this gave the following energetic material constants: = 10-3 m/s-bar A = 12.04 B = 0.0235/K Thus: 1/burning rate = 12.04 0.0235T 0 14. 0 Reciprocal Burning Rate of Composition B vs. Bulk Explosive Temp erature 12. 0 10. 0 8. 0 6. 0 Ra te = 0.2 m m / s a t 300K 4. 0 2. 0 Critic a l T e m pe ra ture = 513K 0. 0 0 100 200 300 400 500 600 Bulk Ex plos iv e T e m pe ra ture, K Pressure Decay Rate: If the volume under consideration is vented, the flow through the vent tends to decrease the pressure. When the interior pressure exceeds the outside pressure by more than 0.8 bar, the flow velocity becomes sonic [6] and a very simple expression for the pressure-decrease results (equation 3). -dp/dt = (A v C D /V) a*p (3)

where: A = vent area, m 2 C D = discharge coefficient, 0.6 to 1.0 V = volume, m 3 a* = flow velocity, m/s P = absolute pressure, bars In the generic equation, the discharge coefficient C D was allowed to equal one, i.e., ideal flow. In actuality, flow through a square-edged orifice results in a coefficient of approximately 0.82 because of the vena contracta formed by the gases exiting the vent hole [7]. The sonic flow velocity of the gases through the vent hole, a*, is computed from the temperature of the products, and is also affected by compressible fluid flow. Thus: a* = (RT/M) 1/2 [k * (2k/k+1) 1/2 * (2/k+1) 1/k-1 ] (4) For a nominal combustion gas mixture with T = 2500K R= 8.31434 J/mol-K M=0.028 kg/mol k=1.27 a* is approximately 725 m/s. This estimate can be improved by knowing the actual product composition of the gases, the specific heat as a function of temperature, the actual flame temperature, which, of course, are different for each explosive material. Critical Vent Area: If the magnitudes of the pressure-decay and pressure-rise terms are equal, a critical condition results in which the pressure remains constant. This condition is met when the ratio of vent area to burning surface area is equal to a constant determined by the explosive constants and the initial temperature. The pressure-rise and pressure-decay equations can be combined. Thus: dp/dt = [(RT B * /M * /(A-BT 0 ) * S B ) (A v C D a*)] * (P/V) (5) If the vent-area to burn-surface-area ratio is less than the critical value, the pressure increases exponentially; if greater, the pressure decreases. Thus, the ratio is computed as: A v /S B = (RT B ) / [M C D a*(a-bt 0 )] (6) For the Composition B explosive cited above, and with an explosive density of 1700 kg/m 3, the predicted critical vent-area to burn-surface-area ratio as a function of bulk temperature is shown in Table 1.

Table 1. Critical Vent Area as a Function of Initial Explosive Temperature T 0 K Critical Ratio Av/S B 273 0.002161 288 0.002305 334 0.002896 Experiments: In NWC experiments with vented burning of Composition B explosive, it was found that using the discharge coefficient of 0.82 gave a conservative prediction of the demarcation between quiescent burning and violent reaction (Figure 2). All of the violent burns lie below the demarcation line while the quiescent ones essentially lie on the predicted line. 0.006 0.005 VEC T est Results Compared to Predicted Critical Vent A rea Critical Condition Comp B Contants C D = 0.82 A v /S B 0.004 0.003 0.002 0.001 Quiescent Burns Violent Pressurization 0.000 250 300 350 400 Initial Temperature, K Figure 2. NWC VEC Test Results Compared to Predicted Critical Vent Area for S B = 11.04 in 2. The US Air Force Weapons Laboratory at Kirtland AFB performed experiments on Composition B explosive similar to those done by NWC but with a somewhat larger initial burning surface area [7]. Their data is shown in Figure 3. Their experiments are in qualitative agreement with the NWC experiments. They found quiescent burns in the range of 0.003 to 0.121 A v /S B and violent pressurizations from 0.0028 to 0.0058 A v /S B.

A ir Force Venting T ests with Composition B 0.014 0.012 NWC Data at Same Temperature as AF Data NWC Data for Higher Initial Temperature 0.010 A v /S B 0.008 0.006 0.004 Air Force Safety Recommendation NWC Comp B Constants 0.002 0.000 Quiescent Burn Violent Reaction Figure 3. Air Force Vented Burning Studies with Composition B Explosive. Initial Bulk Temperature, T o = 300K; S B = 12.57 in 2. NWC Data (purple symbols) at T o = 288K and S B = 11.04 in 2. Summary of Experiments: Vent areas to prevent pressurization and violent reaction in these tests are significantly less than 1% of the burning surface area. Tests were conducted with end-burning test items. This formalism works well for items with bulk temperatures near ambient in particular, it works well in bullet impact of warheads where the bullet hole provides enough vent area to prevent overpressurization when the energetic material ignites and burns. Application to the fast cookoff scenario can be successful if the vent is created at a low enough energetic material bulk temperature. Ballistic Analysis Methodology: Another method for estimating the required venting area to prevent violent reaction relies on classical ballistic analysis. For this exercise, the Minuteman III first stage motor was chosen as an example large rocket motor [8]. Figure 4: Minuteman III First Stage Motor

The following typical aluminized propellant properties were assumed in the calculations: 70 o F Burning Rate: r b = 0.29 (P c /1000) 0.34 Temperature Coefficient: σ p = 0.001/ o F Characteristic Velocity: c* = 5172 ft/s Density: ρ = 0.065 lb/ft 3 Where: P c = chamber pressure in psia r b is burning rate in in/s For the initial analysis, the burning rate was adjusted to a temperature of 702 o F, and a single square-edged orifice was used as the vent. It was assumed that the whole exterior surface of the propellant grain ignited instantaneously between the case and the grain resulting in a burning surface area of 42,620 in 2 ; that all gases exited through the squareedged orifice; and that the motor surface was all at the same temperature. Analysis: The motor weight is 50,550 lb m and at 702 o F, the burning rate is calculated to be 0.546 (P c /1000) 0.34. First, compute the thrust using equation (7). F = P c A t C f η F (7) Where: F = Thrust, lb f A t = Throat area, in 2 (This is the vent size) C f = Thrust coefficient = 1.25 (exit cone with no expansion) η F = Thrust efficiency = 80% (square-edged orifice) Secondly, determine chamber pressure using equation (8). where: S B = the surface area, in 2 a = burning rate coefficient in the equation ap n, in/s g c = gravitational constant, 32.174 lb m -ft/lb f /s 2 n = burning rate exponent in the equation ap n P c = [(S B ρ c*a)/(a t g c )] (1/1-n) (8) We wish to keep thrust to < 80% of stage weight to prevent propulsion. Applying this to equation 7 we get equation (9): 40,202 = P c A t (1.25)(0.8) (9)

Outer Grain Pressure, psia Solving for P c through the use of equation (8) gives (10): P c = [(42629)(0.06519)(5172)(0.0521)/A t (32.174)] 1.515 (10) Which gives the solution: Outer grain pressure, P c = 4.99 psia and a required vent area of A t = 8053 sq. in. This analysis was applied over various temperatures to assess the required vent area. Figure 5 illustrates. It can be seen that early, lower temperature venting is definitely advantageous. 10000 Venting Analysis of Minuteman in Fuel Fire 20.0 Required Vent Area, in 2 8000 15.0 6000 10.0 4000 2000 5.0 Assumption: All gases exit through 1 rough-edged orifice. 0 0.0 200 300 400 500 600 700 800 Ignition Temperature, o F Figure 5. Effect of Ignition at Various Temperatures on Required Vent Area and Outer Grain Pressure. In terms of our original ratio of vent area to burning surface area ratio, A v /S B at various temperatures is shown in Table 2. Table 2. Vent Area to Burn Surface Ratio as a Function of Temperature Temperature, F Temperature, K A v /S B 300 422 0.059 400 477 0.077 500 533 0.106 600 589 0.141 700 644 0.189 702 645 0.190

Figure 6 illustrates the rectangular vent area for the MM III first stage at two extremes of surface temperature 700 o F (8000 in 2 area required) and 360 o F (3000 in 2 vent area). Figure 6. Vent area requirements as function of surface temperature Figure 7 illustrates the difference between the two cases studied in this paper. Vent A rea Ratios vs. T emperature 1 Bondline Surface Burning Higher Surface Temperature Log A v /S B 0.1 0.01 End Burning, Self-Heating Low Bulk Temperature 0.001 0 100 200 300 400 500 600 700 800 Temperature, o F Figure 7. Comparison of required vent area to burning surface area ratios Discussion: It can be seen that it is imperative to vent a cased energetic material subjected to fuel fire threat at as low a temperature as possible, consistent with its operational requirements and some margin of safety. Required vent areas are dramatically increased as temperature rises.

It should be noted that if the grain has a significant bore area and the flame reaches the bore, then increased vent area will be required. It is anticipated that the vent area should be on the side of the motor case rather than on the end to prevent launching the motor. Attempts at neutral thrust (vent in front, nozzle in the rear) have been successful but require an exceptionally uniform fuel fire. References 1. G. F. Kinney and R. G. S. Sewell, Venting of Explosions, NWC TM 2448, Naval Weapons Center. China Lake, Calif., NWC, July 1974. 2. Naval Weapons Center. Composition B Thermal Analysis. China Lake, CA, 3 March 1983. NWC Reg. Memo. 3242-27-83. 3. Johanason and Persson. Detonics of High Explosives. New York, Academic Press (1970), pp. 158-161. 4. K.J. Graham and R.G.S. Sewell, Burning-to-Violent-Reaction Transitions in Explosions, in 1980 JANNAF Propulsion Systems Hazards Subcommittee Meeting. Silver Spring, MD., Chemical Propulsion Information Agency, October 1980, p. 147. CPIA Pub 330. 5. Kenneth J. Graham, Explosive Response to Fragments: Venting Studies, NWC TP 6456, Naval Weapons Center, China Lake, CA, October 1986. AD-A176599. 6. B. W. Anderson. The Analysis and Design of Pneumatic Systems, Wiley, NY, 1967, p. 19. 7. Air Force Weapons Laboratory. Vulnerability of Nuclear Weapon Systems to Fire- Studies of Burning Explosives. Kirtland AFB, NM, AFWL, December 1963. Report RTD-TDR-63-3086 (DASA 1417). 8. Sutton, G.P., Rocket Propulsion Elements: An Introduction to the Engineering of Rockets, Sixth Edition, John Wiley & Sons, Inc., 1992.

A GenCorp Company Mitigation of Fuel Fire Threat to Large Rocket Motors by Venting Kenneth J. Graham, Technical Principal 540-854-2182 ken.graham@aerojet.com October 2010 1

A GenCorp Company The Problem Venting of a container such as a rocket motor or a warhead case is a well-recognized method to potentially reduce the violent response of the system to a fuel fire threat. AMRAAM TIVS ARCAPS 120mm Mortar with Ionomer-filled Vent Many others The problem we are trying to solve is how to protect a large rocket motor, perhaps the size of Minuteman or Peacekeeper, while in the transportation mode. What is the critical vent size to prevent overpressurization and how is it determined?

Large Motor Transport A GenCorp Company

Accidents Happen! A GenCorp Company

A GenCorp Company Solution The basic solution to mitigation by venting is to understand the competition between pressure rise rate and pressure decay rate. For Pressure Rise > Pressure Decay the system reacts violently For Pressure Rise = Pressure Decay the system is critically vented For Pressure Rise < Pressure Decay the system reacts mildly This is what we want!

A GenCorp Company Pressure Rise From interior ballistics, the rate of pressure rise from combustion of an energetic material is given by: dp/dt = RTB/V * dn/dt (1) where dn/dt is the time rate of change of the number of moles of product gases. This equation may be replaced with one in which the variables are more easily measurable. Thus, dp/dt = RT B /V * /M * /(A-BT 0 ) * S B P (2) R = molar gas constant = 8.314 x 10-5 bar - m 3 /mol - K V = volume, m 3 T B = flame temperature, K M = formula mass product gas, kg/mol = density of explosive, kg/m 3 T 0 = bulk temperature of explosive, K,A,B = energetic material constants (see below) S B = burn surface area, m 3 P = absolute pressure, bars

Re cipro ca l Burning Ra te, s/ m m A GenCorp Company Pressure Rise The term [ /(A BT 0 )] represents the variation in burning rate with bulk explosive temperature. Utilizing Andreev s method, we plot the reciprocal of burning rate against bulk explosive temperature. 14. 0 12. 0 10. 0 Reciprocal Burning Rate of Composition B vs. Bulk Explosive Temp erature For Composition B explosive = 10-3 m/s-bar A = 12.04 B = 0.0235/K 8. 0 6. 0 4. 0 2. 0 Rate = 0.2 m m / s at 300K Thus: 1/burning rate = 12.04 0.0235T 0 Critic a l T e m pe ra ture = 513K 0. 0 0 100 200 300 400 500 600 Bulk Ex plos iv e T e m pe ra ture, K

A GenCorp Company Pressure Decay When the interior pressure exceeds the outside pressure by more than 0.8 bar, the flow velocity becomes sonic [6] and a very simple expression for the pressure-decrease results (equation 3). -dp/dt = (A v C D /V) a*p (3) A = vent area, m 2 C D = discharge coefficient, 0.6 to 1.0 V = volume, m 3 a* = flow velocity, m/s P = absolute pressure, bars Flow through a square-edged orifice results in a coefficient of approximately 0.82 because of the vena contracta formed by the gases exiting the vent hole.

A GenCorp Company Pressure Decay The sonic flow velocity of the gases through the vent hole, a*, is computed from the temperature of the products, and is also affected by compressible fluid flow. Thus: a* = (RT/M) 1/2 [k * (2k/k+1) 1/2 * (2/k+1) 1/k-1 ] (4) a* is approximately 725 m/s for a nominal combustion gas mixture with:. T = 2500K R= 8.31434 J/mol-K M=0.028 kg/mol k=1.27

A GenCorp Company Critical Vent Area If the magnitudes of the pressure-decay and pressure-rise terms are equal, a critical condition results The pressure-rise and pressure-decay equations can be combined. dp/dt = [(RT B * /M * /(A-BT 0 ) * S B ) (A v C D a*)] * (P/V) (5) Rearrangement gives the relationship of vent area to burning surface area A v /S B = (RT B ) / [M C D a*(a-bt 0 )] (6) If A v /S B is greater than the critical value, pressure decreases. This is what we seek!

Critical Vent Area Ratio A GenCorp Company For the Composition B explosive cited previously, and with an explosive density of 1700 kg/m 3, the predicted critical vent-area to burn-surface-area ratio as a function of bulk temperature is: Table 1. Critical Vent Area as a Function of Initial Explosive Temperature T 0 K Critical Ratio Av/S B 273 0.002161 288 0.002305 334 0.002896 It doesn t take much vent area to prevent pressurization!

VEC Experiments A GenCorp Company NWC Composition B VEC T est Results Compared to Predicted Critical Vent A rea 0.006 0.005 Critical Condition Comp B Contants C D = 0.82 A v /S B 0.004 0.003 0.002 0.001 Quiescent Burns Violent Pressurization 0.000 250 300 350 400 Initial Temperature, K

VEC Experiments A GenCorp Company AFWL Kirtland Composition B A ir Force Venting T ests with Composition B 0.014 0.012 NWC Data at Same Temperature as AF Data NWC Data for Higher Initial Temperature 0.010 A v /S B 0.008 0.006 0.004 Air Force Safety Recommendation NWC Comp B Constants 0.002 0.000 Quiescent Burn Violent Reaction

A GenCorp Company Summary of Experiments Vent areas to prevent pressurization and violent reaction in these tests are significantly less than 1% of the burning surface area. Tests were conducted with end-burning test items. This formalism works well for items with bulk temperatures near ambient. Application to the fast cookoff scenario may be successful if the vent is created at a low enough energetic material bulk temperature.

A GenCorp Company Ballistic Analysis Minuteman III first stage motor was chosen as the example. The assumed propellant properties: Outer grain surface area: 42,629 sq. in. 70 o F Burning Rate: r b = 0.290 (P c /1000) 0.34 Temperature Coefficient: σ p = 0.001/ o F Characteristic Velocity: c* = 5172 ft/s Density: ρ = 0.0652 lb/ft 3 P c = chamber pressure in psia r b is burning rate in in/s Minuteman III First Stage Motor

A GenCorp Company Ballistic Analysis For the initial analysis, the burning rate was adjusted to a temperature of 702 o F, and a single square-edged orifice was used as the vent. Assumptions: The whole exterior surface of the propellant grain ignited instantaneously between the case and the grain All gases exited through the square-edged orifice The motor surface was all at the same temperature The Stage 1 weight is 50,550 lbf The 702 o F burning rate, r b = 0.546 (P c /1000) 0.34

A GenCorp Company MM III Ballistic Analysis First, compute the thrust using equation (7). F = P c A t C f η F (7) F = Thrust, lb f A t = Throat area, in 2 (NOTE: This is the vent size) C f = Thrust coefficient = 1.25 (exit cone with no expansion) η F = Thrust efficiency = 80% (square-edged orifice) Second, apply definition of the chamber pressure using equation (8) P c = [(S B ρ c*a)/(a t g c )]( 1/1-n) (8) S B = the surface area, in 2 a = burning rate coefficient in the equation ap n, in/s g c = gravitational constant, 32.174 lb m -ft/lb f /s 2 n = burning rate exponent in the equation ap n

MM III Ballistic Analysis A GenCorp Company We wish to keep thrust to < 80% of stage weight to prevent propulsion. Applying this to equation 7 we get equation (9): 40,202 = P c A t (1.25)(0.8) (9) Solving for Pc through the use of equation (8) gives (10): P c = [(42629)(0.06519)(5172)(0.0521)/A t (32.174)] 1.515 (10) The solution: Outer grain pressure, P c = 4.99 psia and a required vent area of A t = 8053 sq. in.

Outer Grain Pressure, psia MM III Ballistic Analysis A GenCorp Company This methodology was applied over a wide range of temperatures. Venting Analysis of Minuteman in Fuel Fire 10000 20.0 Required Vent Area, in 2 8000 15.0 6000 10.0 4000 2000 5.0 Assumption: All gases exit through 1 rough-edged orifice. 0 0.0 200 300 400 500 600 700 800 Ignition Temperature, o F Clearly, lower temperature venting is advantageous!

Vent Area Ratio for MM III Stage 1 A GenCorp Company Table 2. Vent Area to Burn Surface Ratio as a Function of Temperature for MM III Propellant in MM III Case. Temperature, F Temperature, K A v /S B 300 422 0.059 400 477 0.077 500 533 0.106 600 589 0.141 700 644 0.189 702 645 0.190

A GenCorp Company Effect of Surface Temperature at Time of Venting 700 o F -- 8000 in 2 vent area required 360 o F -- 3000 in 2 vent area required

Comparison of Methodologies A GenCorp Company Comparison of required vent area to burning surface area ratios for end burning and surface burning cased energetic grains 1 Vent A rea Ratios vs. T emperature Bondline Surface Burning Higher Surface Temperature Log A v /S B 0.1 0.01 End Burning, Self-Heating Low Bulk Temperature 0.001 0 100 200 300 400 500 600 700 800 Temperature, o F

Summary A GenCorp Company It is imperative to vent a cased energetic material subjected to a fuel fire threat at as low a temperature as possible, consistent with its operational requirements and a margin of safety. Required vent areas increase dramatically as the temperature rises If the grain has a significant bore area and the flame reaches the bore, then increased vent area will be required. Grains that burn cigarette fashion and slowly self-heat require less vent area than those exposed to an engulfing fuel fire where the whole outer surface area is heated. It is anticipated that the vent area should be on the side of the motor case rather than on the end to prevent launching the motor. A ballistics-based methodology has been presented to predict the critical vent area for a motor exposed to a fuel fire.