S. Q. Wang and S. T. Thynell

Similar documents
Decomposition and Combustion of Monomethylhydrazinium Nitrates

First Principles Study of the Ignition Mechanism for Hypergolic Bipropellants: N,N,N 0,N 0 -Tetramethylethylenediamine

Reactive Dynamics Study of Hypergolic Bipropellants: Monomethylhydrazine and Dinitrogen Tetroxide

The Pennsylvania State University. The Graduate School. College of Engineering STUDY OF CONDENSED PHASE REACTIONS BETWEEN

Rocket Propulsion Overview

DISTRIBUTION A: Distribution approved for public release.

Ionic Liquids as Hypergolic Fluids

Combustion. Indian Institute of Science Bangalore

Globally Harmonized System of Classification and Labelling of Chemicals (GHS) Classification criteria for substances and mixtures Physical hazards

The Concept of Equilibrium

Energetic Ionic Liquids

g of CO 2 gas is at a temperature of 45 o C and a pressure of 125 kpa. What is the volume of the container? 11 L

Flame-Spread Characteristics of n-decane Droplet Arrays at Different Ambient Pressures in Microgravity

Rocket Propulsion Prof. K. Ramamurthi Department of Mechanical Engineering Indian Institute of Technology, Madras

Rocket Propulsion. Combustion chamber Throat Nozzle

CHAPTER 1 INTRODUCTION TO HIGH ENERGY MATERIALS. Among exothermic reactions, some are extensively used in dayto-day

1.23 Gas Calculations

Gas Volumes and the Ideal Gas Law

GHS: Physical Hazards

Technology of Rocket

Laminar flame speed (burning velocity) reactants. propagating flame front. products. reactants

A Gas Uniformly fills any container. Easily compressed. Mixes completely with any other gas. Exerts pressure on its surroundings.

DISTRIBUTED BY: National Technical Information Service U. S. DEPARTMENT OF COMMERCE 5285 Port Royal Road. Sprioigfield Va

A).5 atm B) 1 atm C) 1.5 atm D) 2 atm E) it is impossible to tell

3. Write a detailed note on the following thrust vector control methods:

Lecture 9 Hybrid Rocket Propulsion Liquefying Fuels

Energy Changes in Chemical Reactions

Chemical reactions. C2- Topic 5

Stoichiometry SUPPLEMENTAL PROBLEMS CHAPTER 12. 3Si(s) 2N 2 N 4. (g) 0 Si 3. (s) PO 4. the reaction. Cr(s) H 3. (aq) 0.

Globally Harmonized System of Classification and Labelling of Chemicals (GHS)

2012 AP CHEMISTRY FREE-RESPONSE QUESTIONS

Chapter 11 Gases 1 Copyright McGraw-Hill 2009

Chapter 5 Test. Directions: Write the correct letter on the blank before each question.

5.) One mole of an ideal gas expands isothermally against a constant pressure of 1 atmosphere. Which of the following inequalities is true?

10 States of Matter. Aubrey High School AP Chemistry. Period Date / / 10.2 Problems - Liquids and Gases

Chapter 6. Thermochemistry

B 2, C 2, N 2. O 2, F 2, Ne 2. Energy order of the p 2p and s 2p orbitals changes across the period.

Rocket Dynamics. Forces on the Rocket

Defense Technical Information Center Compilation Part Notice

EFFECT OF KOH-ACTIVATED RICE HUSK ON THERMAL DECOMPOSITION OF HYDROXYLAMINE NITRATE MONOPROPELLANT

Characteristics of some monopropellants (Reprinted from H. Koelle, Handbook of Astronautical Engineering, McGraw-Hill, 1961.)

relatively narrow range of temperature and pressure.

Investigating the Phase of High Entropy. Pressure

Chemistry Chapter 16. Reaction Energy

Write a balanced equation for the thermal decomposition of calcium nitrate. ... (2)

Chemically-Augmented Pulsed Laser-Ramjet

U N I T T E S T P R A C T I C E

Presentation Start. Zero Carbon Energy Solutions 4/06/06 10/3/2013:; 1

Combustion and Flame

Name AP Chemistry / / Chapter 5 Collected AP Exam Free Response Questions Answers

Chapter 3 Test Bank. d. The decomposition of magnesium oxide produces 2.4 g of magnesium metal and 3.2 g of oxygen gas.

Kc is calculated for homogeneous reactions using the concentrations of the reactants and products at equilibrium:

AP Chapter 5: Gases Name

Unit 7: Stoichiometry Homework Packet (85 points)

Additional Calculations: 10. How many joules are required to change the temperature of 80.0 g of water from 23.3 C to 38.8 C?

Chapter 5 Thermochemistry

PLEASE USE THE INFORMATION TECHNOLOGY DIVISION MCQ ANSWER SHEET FOR QUESTION 1

Chemistry 40S Chemical Kinetics (This unit has been adapted from

Experimental and QM Investigation of Decomposition and Early Chemical Kinetics of Ammonia Borane

Droplet Evaporation Behavior of Kerosene/Nano Aluminum Fuels at High Pressure Environment

Useful Information to be provided on the exam: 1 atm = 760 mm Hg = 760 torr = lb/in 2 = 101,325 Pa = kpa. q = m C T. w = -P V.

Gas UNIT - READINESS ASSESSMENT QUIZ (RAQ) THIS QUIZ WILL BE PACED WITH CLICKER QUESTIONS

Test Booklet. Subject: SC, Grade: HS CST High School Chemistry Part 2. Student name:

AP* Chapter 6. Thermochemistry

11.1 Survey of Spacecraft Propulsion Systems

SUPPLEMENTARY TOPIC 3 ENERGY AND CHEMICAL REACTIONS

HIGH PRESSURE METHANE-OXYGEN COMBUSTION KINETIC ANALYSIS

Chapter 10 Practice. Name: Multiple Choice Identify the choice that best completes the statement or answers the question.

Chemistry Lab Fairfax High School Invitational January 7, Team Number: High School: Team Members Names:

13.1 States of Matter: A Review 13.2 Properties of Liquids A. Evaporation B. Vapor Pressure C. Surface Tension 13.3 Boiling Point and Melting Point

Chemistry Entrance Material for Grade 10 to

REPORT DOCUMENTATION PAGE

10/16/2018. Why study gases? An understanding of real world phenomena. An understanding of how science works.

Name Quarterly Practice # 1 Period

Chemistry 104 Final Exam Content Evaluation and Preparation for General Chemistry I Material

Section 1 - Thermochemistry

Q1. (a) State what is meant by the term activation energy of a reaction. (1)

Unit 6 Solids, Liquids and Solutions

Chapter 9. Chemical Equilibrium

Lecture 9 Laminar Diffusion Flame Configurations

Unit 10 Stoichiometry Chapter 11 of your textbook

Gas Volumes and the Ideal Gas Law

Energy in Chemical Reaction Reaction Rates Chemical Equilibrium. Chapter Outline. Energy 6/29/2013

Chapter 11. Preview. Lesson Starter Objectives Pressure and Force Dalton s Law of Partial Pressures

Name SCS- Date. Experiment 5: Hydrogen Formation and Reaction with Oxygen

17.2 Thermochemical Equations

10 Enthalpy changes Answers to Activity and Practice questions

ENERGETICS. Energy changes which take place during a chemical reaction can be represented on an energy (enthalpy) profile diagram.

Chapter 15 Equilibrium

12.2. The Ideal Gas Law. Density and Molar Mass of Gases SECTION. Key Terms

AP CHEMISTRY 2009 SCORING GUIDELINES

Review Article Advances in Hypergolic Propellants: Ignition, Hydrazine, and Hydrogen Peroxide Research

Atomization. In Flame Emission

Multiple Choice Identify the letter of the choice that best completes the statement or answers the question.

1. Balance the following chemical equations: a. C 8 H 18 + O 2 à CO 2 + H 2 O. b. B 5 H 9 + O 2 à B 2 O 3 + H 2 O. c. S 8 + Cl 2 à S 2 Cl 2

CHEMISTRY. Section II (Total time 95 minutes) Part A Time 55 minutes YOU MAY USE YOUR CALCULATOR FOR PART A.

Lecture 1: Vapour Growth Techniques

Chemistry 102b First Exam. CIRCLE the section in which you are officially registered:

CHAPTER 13. States of Matter. Kinetic = motion. Polar vs. Nonpolar. Gases. Hon Chem 13.notebook

(E) half as fast as methane.

Transcription:

Paper # 070HE-0247 Topic: Heterogeneous Combustion 8 th U. S. National Combustion Meeting Organized by the Western States Section of the Combustion Institute and hosted by the University of Utah May 19 22, 2013 Experimental Investigation of Pressure Effect on Ignition Delay of Monomethylhydrazine, 1,1 Dimethylhydrazine, Tetramethylethylenediamine and 2 Dimethylaminoethylazide with Nitric Acid S. Q. Wang and S. T. Thynell Department of Mechanical and Nuclear Engineering, The Pennsylvania State University, University Park, Pennsylvania 16802, USA A drop-on-pool impingement setup, coupled with a high-speed camera, was used to investigate the ignition delay of several hydrazine-, amine- and azide-based fuels, including monomethylhydrazine (MMH), 1,1-dimethylhydrazine (UDMH), tetramethylethylenediamine (TMEDA) and 2-dimethylaminoethylazide (DMAZ), with white fuming nitric acid (WFNA) as the oxidizer. The drop-on-pool impingement experiments were conducted in an N 2 purged chamber with pressures ranging from -40 kpa to 600 kpa (gauge). In general, the ignition delay of hydrazines (MMH and UDMH) decreases with increasing pressure, whereas the ignition delay of TMEDA and DMAZ increases with increasing pressure. In addition to the pressure effect, the ignition delay of MMH/WFNA is also significantly affected by the type of droplet-pool interaction. The droplet stays within the pool after impingement at sub-atmospheric pressures (type I), while the droplet is ejected from the pool at elevated pressures (type II). The critical pressure is about 0 kpa (g) at which either interaction type was observed. The ignition delay of type I interaction is much longer than that of type II at that pressure. The ignition delay of a mixture of 33.3% DMAZ and 66.7% TMEDA was also studied. This mixture has a shorter ignition delay than both pure TMEDA and DMAZ. In addition, the pressure effect on the ignition delay of this mixture is not significant. 1. Introduction Hypergolic propellants are pairs of liquid fuels and oxidizers in which ignition occurs spontaneously upon contact between the two liquids, thereby eliminating the need for a complex ignition system [1]. The reliable restart capability of these types of engines makes them ideal for spacecraft maneuvering. An important criterion to evaluate the performance of hypergolic bipropellants is the ignition delay, which is defined as the time from the physical contact of liquid fuel and oxidizer to the achievement of ignition. A long ignition delay can cause a hard-start problem in which unburnt fuel and oxidizer in the combustion chamber accumulate, then suddenly ignite to generate damaging pressure spikes [2]. One major factor that determines the ignition delay is the pre-ignition chemical reactions of the two reactants. For example, monomethylhydrazine (MMH, CH 3 NHNH 2 ) has a shorter ignition delay than tetramethylethylenediamine [TMEDA, (CH 3 ) 2 NCH 2 CH 2 (CH 3 ) 2 ] when red fuming nitric acid (RFNA) is used as the oxidizer, because the hydrogen abstraction from an amino group by NO 2 has a lower barrier than that from the methyl group [3-5]. Another major factor that can significantly affect the ignition delay is the extent of initial contact and mixing between liquid fuel and oxidizer. For the same fuel and oxidizer pair, the ignition delay in an engine test is shorter than that obtained from a drop test [6], because the engine test has a better mixing than drop test. In addition, ignition delay obtained by the same method (i.e., drop test) may also vary from one to another because it may be sensitive to many other factors such as ambient pressure and temperature, size of droplet and pool, and even the relative speed of impingement. Therefore, the measurement of ignition delay is meaningful only when the test method and conditions are clarified. In an earlier work, a drop-on-pool impingement setup was developed to study the hypergolic interaction and measure the ignition delay of various hypergolic pairs in an open environment and at ambient pressure [7, 8]. In this work, however, a

modified drop test setup is developed to study the pressure effect on the ignition delay of four fuels with white fuming nitric acid (WFNA). The molecular structure of the four fuels, monomethylhydrazine (MMH), 1,1-dimethylhydrazine (UDMH), tetramethylethylenediamine (TMEDA) and 2-dimethylaminoethylazide (DMAZ), are shown in Fig. 1. The two hydrazines are well-known fuels which have been deployed for decades in rocket engines and spacecraft [1]. TMEDA and DMAZ are less toxic than hydrazines and are considered as potential alternative fuels for MMH [8-10]. In addition, a mixture of DMAZ and TMEDA, with a mass ratio of 2:1, was claimed to have both a density impulse and ignition delay that are comparable to MMH [6]. Therefore, the pressure effect on the ignition delay of this mixture is also examined in this work. 2. Experimental Setup NH NH 2 (MMH) N NH 2 (UDMH) N CH 2 CH 3 CH 2 N CH 3 (TMEDA) Figure 1. Molecular structures of the considered fuels N CH 2 CH 2 N N + N - (DMAZ) A schematic diagram of the drop-on-pool impingement setup is shown in Fig. 2. About 100 µl WFNA is deposited in a glass vial (inner diameter is 0.5 inch). The glass vial is then placed in a bottom cylindrical chamber with an inner diameter of 2 inches and length of 5 inches. A syringe dispenses the fuel, and it is placed above the glass vial and points to the center of the oxidizer pool. The syringe is kept in the top cylindrical chamber, and the drop is released by the motion from the piston of a pneumatic actuator. The chamber is first purged by N 2 and is then either pressurized by the same purge gas or the pressure is reduced by a vacuum pump. In each test, about 10 µl fuel is loaded in the syringe so that only one droplet (about 7 µl) is produced. The drop-on-pool impingement interaction as well as the subsequent gas-phase ignition and combustion processes were recorded by a Phantom V710 high-speed camera, through an optical access glass window with a thickness of 1 inch. A He-Ne laser and a corresponding photodiode, placed underneath the syringe but above the glass vial, triggers the camera through a reduction in the photodiode signal as the drop of fuel intersects the laser beam. The ignition delay is defined as the time delay between the contact of two liquids and the occurrence of luminosity. The high speed videos were acquired at a frame rate of 5000 fps. Therefore the temporal resolution of ignition delay measurement is 0.2 ms. He-Ne laser syringe purge flow distributor purge outlet piston (driven by pneumatic actuator) photodiode fuel droplet WFNA N 2 gas cylinder vacuum pump Figure 2. Schematic diagram of strand burner 3. Results and Discussion 3.1 MMH and UDMH The drop-on-pool impingement interactions and ignition delays of MMH and UDMH with WFNA were studied from sub-atmospheric to elevated pressures. For MMH/WFNA, two different types of drop-pool interactions were observed. At pressures below 0 kpa (g), the droplet stays in the nitric acid pool after the impact (type I), while at elevated pressures the droplet is ejected from the oxidizer pool and then suspended above the pool surface (type II). Typical physical observations from the two types of interaction are illustrated in Fig. 3a and 3b, which show selected images from a test at reduced pressure (-20 kpa (g)) and a test at elevated pressure (50 kpa), respectively.

Type I: as shown in Fig. 3a, the droplet impinges on the nitric acid pool at t = 0. An aerosol cloud, which is believed to be MMH nitrates [7], is formed along the path of the falling droplet through the reactions between the fuel droplet and the nitric acid vapor that is confined in the glass vial. Upon impact, liquidphase reactions will occur at the contact surface which forms nitrates, generate abundant heat, and produce plenty of gases. It is generally believed that a gas layer (or vapor layer) will be formed between the droplet and pool surface [11]. In this case, the gas flow is not strong enough to eject the droplet from the pool. Meanwhile, the vapors of the two reactants react above the pool to form a nitrate cloud as shown at t = 80 ms. When the temperature rises to a certain level, the accumulated nitrate cloud decomposes rapidly followed by a gas-phase ignition which occurs at t = 87 ms. A stable cone-shape diffusion flame sustained for about 350 ms. Type II: as shown in Fig. 3b, the reactions between the droplet and pool surface generate a gas flow that is strong Ignition dealy, ms 100 10 1 0.1 0 100 200 300 400 Gauge Pressure, kpa (a) t = 0 t = 80 ms t = 87 ms t = 300 ms (b) t = 0 t = 8 ms t = 20 ms t = 200 ms (c) t = 0 t = 4.6 ms t =10 ms t = 20 ms Figure 3. Droplet-on-pool impingement: a) MMH/WFNA at -20 kpa; b) MMH/WFNA at 50 kpa; 3) UDMH/WFNA at 0 kpa (gauge pressure) enough to eject the fuel droplet from the nitric acid pool at t = 8 ms. Luminosity is observed at the moment of ejection, and a diffusion flame is formed between the droplet and the pool. With the establishment of a force balance between the droplet gravity and the repelling force from gas flow, the MMH droplet is floating on and rolling along the WFNA pool surface smoothly for about 200 ms followed by rapid disintegration into many smaller droplets. The relation between the gravity (G) of droplet and the repelling force (F) from the gas flow determines whether the droplet will be ejected or not. If F > G, the droplet will be ejected. Gravity (G) is a constant in all tests if one assumes that applied chamber pressure has negligible effect on size of droplet, which is reasonable because surface tension is usually not sensitive to pressure. Therefore, the gas flow, which generates the repelling force, should increase with UDMH / WFNA MMH / WFNA Figure 4. Ignition delay of MMH/WFNA and UDMH/WFNA pressure to explain the above observations. The gas flow may come from two paths: 1) the gases generated by the liquid-phase reactions, which should be independent of pressure because the density change is negligible; and 2) the gases formed by the gas-phase reactions in the vapor layer, which should increase with pressure because the vapor concentrations should increases with pressure. The existence of a very thin vapor layer between the droplet and pool surface is claimed by Daimon et al. [11] who visualized the layer using a special technique. Theoretically, there should be a critical pressure P c, at which the repelling force (F c ) is equal to gravity (G). A small perturbation may lead to totally different droplet-on-pool interactions and cause huge difference on ignition delays. It is interesting to note that this critical pressure is

around 0 kpa (g) for both interaction types I and II, with an ignition delay of 80 and 14 ms, respectively. All tests at subatmospheric pressures of -20 and -40 kpa (g) follows type I, whereas all tests at elevated pressures of 50, 100, 200, 300 and 400 kpa (g) follows type II interaction. The ignition delays of MMH/WFNA at these pressures are plotted in Fig. 4. It should be noted the ignition delays measured in this work have a resolution of 0.2 ms since a frame rate of 5000 fps is used. The data in the dashed frame in Fig. 4 indicate that the ignition delay is shorter than 0.2 ms. A gas-phase ignition is usually achieved through the gasification of reactants and subsequent gas-phase reactions, therefore the chamber pressure can affect ignition delay through two competing mechanisms. On one hand, gasification of reactants requires more heat at higher chamber pressures since boiling point increases with pressure. The delay in gasification with increasing pressure may thus cause a delay in ignition. One the other hand, increased pressure may have a positive effect on gas-phase reactions through increasing the concentrations of gaseous reactants. In the case of MMH/WFNA, the second mechanism controls the ignition delay, and therefore it decreases with increasing pressure. As mentioned earlier, interaction type I has a much longer ignition delay than type II even at the same pressure; a plausible explanation is not available. The speculation is made that when the droplet stays in the pool, a layer of nitrate salts may form between the two liquid reactants, which prevents or reduces the rates of the reactions and lead a longer ignition delay. However, small variations in the experimental conditions, such as slightly different droplet size or concentration of reactants in the vial, may affect the thickness of the nitrate layer and contribute to the uncertainty in the type of the observed interaction. Ignition delay of UDMH/WFNA at pressures ranging from -40 kpa to 400 kpa (g) is also studied. Figure 3c shows selected images from a drop-on-pool impingement test at 1 atm. UDMH reacts with WFNA more violently and have a shorter ignition delay compared to MMH/WFNA. The droplet is ejected from the pool and breaks into many smaller droplets. The repelling force from the gas flow is larger than the gravity of droplet (F > G) even at a reduced pressure of -40 kpa (g), therefore only one type of interaction (droplet ejection) was observed in all tests. It is not clear whether a critical pressure will also exist at even further reduced pressures. Similar to MMH/WFNA, the ignition delay of UDMH/WFNA decreases with increasing pressure, as shown in Fig. 4. 3.2 TMEDA, DMAZ and their mixture The drop-on-pool impingement interactions and the ignition delays of a tertiary amine TMEDA, and an azide DMAZ and a mixture of 66.7% TMEDA with 33.3% DMAZ were studied at various reduced and elevated pressures. Figures 5a, b and c show selected Ignition dealy, ms 100 DMAZ / WFNA TMEDA / WFNA (a) t = 0 100 ms 218 ms 230 ms 10 Mixture / WFNA (b) t = 0 33 ms 50 ms 100 ms (c) t = 0 5 ms 11 ms 20 ms Figure 5. Droplet-on-pool impingement at 500 kpa: a) DMAZ/WFNA, b) TMEDA/WFNA, and c) Mixture/WFNA 0 100 200 300 400 500 600 Gauge Pressure, kpa Figure 6. Ignition delay of TMEDA/WFNA, DMAZ/WFNA, and their mixture/wfna images from drop tests of DMAZ/WFNA, TMEDA/WFNA, and mixture/wfna at 500 kpa, respectively. H-abstraction from C-H has a much higher barrier (13.4 kcal/mol) [3] than that from N-H (8.3 kcal/mol) [4], therefore a gas flow from the rapid gas-phase reactions, which is strong enough to eject the fuel droplet, was not formed. Instead, the early reactions between TMEDA (or DMAZ) vapor and nitric acid vapor form a dense nitrate cloud, as shown Figure 5a and b. The accumulation of heat from nitrates formation finally leads to

the decomposition of nitrates and a gas-phase ignition. The ignition delays of TMEDA, DMAZ, and their mixture at pressures from -40 to 600 kpa (g) were measured and plotted in Fig. 6. In general, the ignition delay increases with increasing pressure. As discussed earlier, the pressure shows an opposite effect on ignition delays of hydrazines, for which the controlling factor on ignition delays is gas-phase reactions, such as H abstraction from amino groups. However, in the case of TMEDA and DMAZ, H abstraction from methyl groups is relatively slow at pre-ignition temperatures. Instead, the controlling factor on ignition delay is the temperature rise from nitrate-salt formation reactions in the liquid phase, which takes a longer time at higher pressures due to higher boiling points. This may also explain why TMEDA has a shorter ignition delay than DMAZ. More heat is released by TMEDA nitrate formation compared to DMAZ nitrate formation, because one mole of TMEDA can react with two moles of HNO 3 to form a dinitrate [8, 12]. It should be noted that the evolved gas-phase reactants, including HNO 3 and the fuels, are rapidly consumed to form a particulate cloud of ionic nitrates, and these reactions are also exothermic. The mixture of TMEDA and DMAZ has a shorter ignition delay than pure TMEDA and DMAZ, as shown in Fig. 6. Compared to pure DMAZ, the addition of TMEDA to DMAZ increases the heat release from salt formation reactions, thus reduces the ignition delay. Compared to pure TMEDA, the addition of DMAZ brings in an azide group which may react at relatively lower temperatures than the temperature required to decompose TMEDA salts or other secondary reactions such as H abstraction from a methyl group. A combination of these two factors leads to a shorter ignition delay. In addition, the ignition delay of this mixture shows little dependence on pressure. The ignition delay is around 10 ms which is close to the value (9 ms) obtained in Stevenson s work [6]. It is should be noted that for DMAZ/WFNA, the ignition delay of DMAZ/WFNA increases with decreasing pressure at reduced pressures, as shown in Fig. 6. A minimum ignition delay is observed at about 0 kpa (g). 4. Conclusions The pressure effect on ignition delays of several fuels with nitric acid was examined. Pressure may affect the ignition delay through two different mechanisms. 1) Increasing pressure can reduce ignition delay if secondary gas-phase reactions, such as H abstraction reactions, are important and are the controlling factors of ignition delay. Typical examples are MMH and UDMH. 2) Increasing pressure can delay the gasification of reactants due to an increased boiling point, and therefore, it can cause a longer ignition delay if the ignition is mainly determined by the heat release from salt formation reactions in the liquid phase. Typical examples are TMEDA and DMAZ, since their secondary reactions (i.e., H abstraction from C-H) are very slow at pre-ignition temperatures. In addition to the pressure effect, a change of physical interaction type between droplet and pool can cause a huge difference on ignition delay. An example is MMH/WFNA at 0 kpa (g), which has an ignition delay of about 14 ms if the droplet is ejected and 80 ms if the droplet stays in the pool. Another interesting finding in this work is that a mixture of TMEDA and DMAZ has a shorter ignition delay than both pure TMEDA and DMAZ. The ignition delay of this mixture shows little dependence on pressure, and thus controlled by liquid-phase reactions whose rates are unaffected by the applied pressure. Acknowledgements This material is based upon work supported by, or in part by, the U.S. Army Research Laboratory and the U.S. Army Research Office under grant number W911NF-08-1-0124. References [1] G. P. Sutton, History of Liquid Propellant Rocket Engines, American Institute of Aeronautics and Astronautics, Alexandria, VA, 2006, 5-21. [2] R. T. Holtzmann, Chemical Rockets, Marcel Dekker: New York, 1969, 123 124. [3] W.-G. Liu, S. Dasgupta, S. V. Zybin, W. A. Goddard, Journal of Physical Chemistry A 115(2011) 5221-5229. [4] M. J. McQuaid, Y. Ishikawa, Journal of Physical Chemistry A, 110(2006) 6129-6138. [5] W.-G. Liu, S. Q. Wang, S. Dasgupta, S. T. Thynell, W. A. Goddard, S. Zybin, R. A. Yetter, Combustion and Flame 160(2013) 970-981. [6] W. H. Stevenson, Hypergolic Liquid or Gel Fuel Mixtures. U.S. Patent No. 0127551, 2008. [7] S. Q. Wang and S. T. Thynell, Combustion and Flame 159(2012) 438-447. [8] S. Q. Wang and S. T. Thynell, A. Chowdhury, Energy & Fuels 24(2010) 5320-5330. [9] D. M. Thompson, Tertiary Amine Azides in Hypergolic Liquid or Gel Fuels Propellant Systems. U.S. Patent No. 6013143, 2000.

[10] M. J. McQuaid, W. H. Stevenson, D. M. Thompson, Computationally Based Design and Screening of Hypergolic Multiamines Technical Report, NTIS No. ADA433347; U.S. Army Research Laboratory: Aberdeen Proving Ground, MD, 2004. [11] W. Daimon, Y. Gotoh, I. Kimura, Journal of Propulsion and Power 7(1991) 946-952. [12] S. Q. Wang, S. T. Thynell, ACS Symposium Series 1117(2012) 51-66