An LX-10 kinetic model calibrated using simulations of multiple small-scale thermal safety tests

A new chemical kinetic model for the beta-delta transition and decomposition of LX-10 (95% octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine, 5% Viton A binder) is presented here. This model implements aspects of previous kinetic models but calibrates the model parameters to data sets of three experiments: differential scanning calorimetry (DSC), thermogravimetric analysis (TGA), and one-dimensional time to explosion (ODTX). The calibration procedure contains three stages: one stage uses open-pan DSC and TGA to develop a base reaction for formation of heavy gases, a second stage features closed-pan DSC to ascertain the autocatalytic behavior of reactant gases attacking the solid explosive, and a final stage adjusts the rate for the breakdown of heavy reactant gases using ODTX experimental data. The resultant model presents a large improvement in the agreement between simulated DSC and TGA results and their respective experiments while maintaining the same level of agreement with ODTX, scaled thermal explosion, and laser heating explosion times when compared to previous models.     

A comparison of direct methods to determine n-th order kinetic parameters of solid thermal decomposition for use in fire models

Pyrolysis models for burning solids in fire simulations are sensitive to the values of the activation energy, frequency factor, and reaction order that characterize the thermal decomposition of the solid to gaseous fuel, so direct measurement of these kinetic parameters is recommended, and simple methods are preferred. Three direct methods were evaluated with regard to the ability of their kinetic parameters to reproduce the thermal decomposition of five polymers measured by differential thermogravimetric analysis using the reaction order model. It was found that the two multiple heating rate methods produced identical, physically based kinetic parameters, while the peak property method produced non-physical kinetic parameters. However, all of these kinetic parameters in a single-step reaction order model gave reasonably good conversion histories for non-charring and moderately charring polymers. For a highly charring polymer, the conversion histories were poorly described without a multiple step reaction. The temperature at the maximum rate of conversion was found to be essentially independent of the reaction order, which decouples the frequency factor from the reaction order in the direct kinetic methods. Any of the direct methods are sufficiently accurate to obtain kinetic parameters for pyrolysis models because of the inherent spatial and temporal averaging of reaction rates at the burning surface of a thick solid and the uncertainty in the heat transfer mechanisms and thermo-physical parameters used in the models.     

Application of global kinetic models to HMX beta-delta transition and cookoff processes

The reduction of the number of reactions in kinetic models for both the HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine) beta-delta phase transition and thermal cookoff provides an attractive alternative to traditional multi-stage kinetic models due to reduced calibration effort requirements. In this study, we use the LLNL code ALE3D to provide calibrated kinetic parameters for a two-reaction bidirectional beta-delta HMX phase transition model based on Sandia instrumented thermal ignition (SITI) and scaled thermal explosion (STEX) temperature history curves, and a Prout-Tompkins cookoff model based on one-dimensional time to explosion (ODTX) data. Results show that the two-reaction bidirectional beta-delta transition model presented here agrees as well with STEX and SITI temperature history curves as a reversible four-reaction Arrhenius model yet requires an order of magnitude less computational effort. In addition, a single-reaction Prout-Tompkins model calibrated to ODTX data provides better agreement with ODTX data than a traditional multistep Arrhenius model and can contain up to 90% fewer chemistry-limited time steps for low-temperature ODTX simulations. Manual calibration methods for the Prout-Tompkins kinetics provide much better agreement with ODTX experimental data than parameters derived from differential scanning calorimetry (DSC) measurements at atmospheric pressure. The predicted surface temperature at explosion for STEX cookoff simulations is a weak function of the cookoff model used, and a reduction of up to 15% of chemistry-limited time steps can be achieved by neglecting the beta-delta transition for this type of simulation. Finally, the inclusion of the beta-delta transition model in the overall kinetics model can affect the predicted time to explosion by 1% for the traditional multistep Arrhenius approach, and up to 11% using a Prout-Tompkins cookoff model.     

Thermal decomposition of monomethylhydrazine: Shock tube experiments and kinetic modeling

The thermal decomposition of gaseous monomethylhydrazine (MMH) was studied by recording MMH absorption at 220 nm of the reacting gas behind a reflected shock wave at temperatures of 900–1370 K, pressures of 140–450 kPa, and in mixtures containing 97.5–99 mol% argon. Based on previous work (Sun and Law; J Phys Chem A 2007, 111(19), 3748–3760), a kinetic mechanism was developed over extended temperature and pressure ranges to model these experimental data. Specifically, the temperature and pressure dependence of the unimolecular rate coefficients on the dissociation of MMH and the associated radicals were calculated by the QRRK/Master equation analysis at temperatures of 300–2000 K and pressures of 1–100 atm based on published thermochemical and kinetic parameters. They were then fitted using the Troe formalism and incorporated in the kinetic model. This unadjusted model was then used to predict the MMH decomposition profiles at different temperatures and pressures for seven groups of MMH/Ar mixtures and the half‐life decomposition times from shock tube experiments. Good agreement was observed below 940 K and above 1150 K for the diluted MMH/Ar mixtures. The model predictions further show that the overall MMH decomposition rate follows first‐order kinetics, and that the N–N bond scission is the most sensitive reaction path for the modeling of the homogeneous decomposition of MMH at elevated pressures. However, the model predictions deviate from the experimental data with the incubation period of ca. 100 μs observed in the 1030–1090 K temperature range, and it also predicts longer ignition delays for highly concentrated MMH/Ar mixtures. The discrepancy between the model predictions and experimental data at these special conditions of MMH decomposition was analyzed. © 2008 Wiley Periodicals, Inc. Int J Chem Kinet 41: 176–186, 2009     

Ab initio kinetics of gas phase decomposition reactions

The thermal and kinetic aspects of gas phase decomposition reactions can be extremely complex due to a large number of parameters, a variety of possible intermediates, and an overlap in thermal decomposition traces. The experimental determination of the activation energies is particularly difficult when several possible reaction pathways coexist in the thermal decomposition. Ab initio calculations intended to provide an interpretation of the experiment are often of little help if they produce only the activation barriers and ignore the kinetics of the decomposition process. To overcome this ambiguity, a theoretical study of a complete picture of gas phase thermo-decomposition, including reaction energies, activation barriers, and reaction rates, is illustrated with the example of the β-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) molecule by means of quantum-chemical calculations. We study three types of major decomposition reactions characteristic of nitramines: the HONO elimination, the NONO rearrangement, and the N-NO(2) homolysis. The reaction rates were determined using the conventional transition state theory for the HONO and NONO decompositions and the variational transition state theory for the N-NO(2) homolysis. Our calculations show that the HMX decomposition process is more complex than it was previously believed to be and is defined by a combination of reactions at any given temperature. At all temperatures, the direct N-NO(2) homolysis prevails with the activation barrier at 38.1 kcal/mol. The nitro-nitrite isomerization and the HONO elimination, with the activation barriers at 46.3 and 39.4 kcal/mol, respectively, are slow reactions at all temperatures. The obtained conclusions provide a consistent interpretation for the reported experimental data.     

Oxygen concentration and modeling thermal decomposition of a high‐performance material: A case study of polyimide (Cirlex)

Kinetic decomposition models for the thermal decomposition of a high‐performance polymeric material (Polyimide, PI) were determined from specific techniques. Experimental data from thermogravimetric analysis (TGA) and previously elucidated decomposition mechanism were combined with numerical simulating tool to establish a comprehensive kinetic model for the decomposition of PI under three atmospheres: nitrogen, 2% oxygen, and synthetic air. Multistaged kinetic models with subsequent and competitive reactions were established by taking into consideration the different types of reactions that may be occurring during the thermal decomposition of the material (chain scission, thermo‐oxidation, char formation). The decomposition products and decomposition mechanism of PI which was established in our previous report allowed for the elucidation of the kinetic decomposition models. A three‐staged kinetic thermal decomposition pathway was a good fit to model the thermal decomposition of PI under nitrogen. The kinetic model involved an autocatalytic type of reaction followed by successive nth order reactions. Such types of models were set up for the evaluation of the kinetics of the thermal decomposition of PI under 2% oxygen and in air, leading to models with satisfactory fidelity.     

Influence of particle size and kinetic parameters on tire pyrolysis

One of the most important parameters that can significantly affect the cost of the tire pyrolysis process is the time needed for thermal decomposition of the tire material. In this work, the influence of particle size and kinetics of thermal decomposition on the pyrolysis time was studied. The apparent kinetic parameters of tire thermal decomposition were estimated using three different approaches based on thermogravimetry measurements. In separate experiments, tire particles with different sizes were pyrolyzed under isothermal conditions in a laboratory flow reactor at different residence times of the particles in the reactor and the data recorded were employed in the determination of tire conversion during the thermal decomposition. A mathematical model of tire pyrolysis considering heat conduction in the tire particles was developed. All three types of estimated kinetic parameters were used to determine the conversion behavior at isothermal conditions and the results were compared with experimental data obtained. The model was used also to calculate the pyrolysis time in an industrial scale reactor at different temperatures and particle sizes.     

A kinetic study of the thermal degradation of 3-methylaminopropylamine inside AlPO<Subscript>4</Subscript>-21

Kinetics of the thermal decomposition of 3-methylaminopropylamine which was used as a structure-directing agent in the synthesis of AlPO₄-21 has been studied under isothermal and non-isothermal conditions. The decomposition is a single-step reaction occurring in the 573–663 K range. It is a phase-boundary-controlled process, described by the ‘F2/3, R3’ kinetic model. The activation energy values obtained under the non-isothermal and isothermal conditions lie in the 173–151 kJ mol^–1 range.     

Comparative kinetic study of decomposition of some diazepine derivatives under isothermal and non-isothermal conditions

Thermal analysis is one of the most widely used methods for studying the solid state of pharmaceutical substances. TG/DTG and DSC curves provide important information regarding the physical properties of the pharmaceutical compounds (stability, compatibility, polymorphism, kinetic analysis, phase transitions etc.). The purpose of a kinetic investigation is to calculate the kinetic parameters and the kinetic model for the studied process. The results are further used to predict the system’s behaviour in various circumstances. A kinetic study regarding the diazepam, nitrazepam and oxazepam thermal decomposition was performed, under non-isothermal and isothermal conditions and in a nitrogen atmosphere, for the temperature steps: 483, 498, 523, 538 and 553 K. The TG/DTG data were processed by three methods: isothermal model-fitting, Friedman’s isothermal-isoconversional and Nomen-Sempere non-parametric kinetics. In the model-fitting methods the kinetic triplets (f(α), A and E _a) that defines a single reaction step resulted in being at variance with the multi-step nature of diazepines decomposition. The model-free approach represented by isothermal and non-isothermal isoconversional methods, gave dependences of the activation energies on the extent of conversion. It is very difficult to obtain an accord with the similar data which resulted under non-isothermal conditions from a previous work. The careful treatment of the kinetic parameters obtained in different thermal conditions was confirmed to be necessary, as well as a different strategy of experimental data processing.     

Thermokinetic parameters and thermal hazard evaluation for three organic peroxides by DSC and TAM III

The thermokinetic parameters were investigated for cumene hydroperoxide (CHP), di-tert-butyl peroxide (DTBP), and tert-butyl peroxybenzoate (TBPB) by non-isothermal kinetic model and isothermal kinetic model by differential scanning calorimetry (DSC) and thermal activity monitor III (TAM III), respectively. The objective was to investigate the activation energy (E _a) of CHP, DTBP, and TBPB applied non-isothermal well-known kinetic equation to evaluate the thermokinetic parameters by DSC. We employed TAM III to assess the thermokinetic parameters of three liquid organic peroxides, obtained thermal runaway data, and then used the Arrhenius plot to obtain the E _a of liquid organic peroxides at various isothermal temperatures. In contrast, the results of non-isothermal kinetic algorithm and isothermal kinetic algorithm were acquired from a highly accurate procedure for receiving information on thermal decomposition characteristics and reaction hazard.     

迪开石热分解及重结晶的非等温动力学

采用热分析(TG-DTG-DSC)和X射线衍射(XRD)技术研究了福州迪开石在动态空气气氛条件下的热分解及重结晶过程, 利用TG和DSC数据分别对迪开石的热分解和重结晶过程进行动力学分析. 由Friedman法求得热分解过程的表观活化能, 以此为初始值, 通过非线性回归法拟合得到了热分解过程最可能的反应机理和动力学参数; 将迭代的等转化率法和最小偏差法相结合计算得到了高温下迪开石重结晶过程最可能的反应机理和动力学参数. 研究结果表明, 迪开石在450~750 ℃内发生热分解, 脱去2个水分子, 生成无定形的准迪开石, 此过程为两步连串反应, 第一步为随机成核和随后生长机理An, 第二步为自催化反应机理CnC; 在975~1050 ℃内, 准迪开石重结晶转化为莫来石, 此过程对应的反应机理为An.     

Thermal and kinetic evaluation of biodiesel derived from soybean oil and higuereta oil

The purpose of this work is to evaluate the thermal and kinetic behavior of biodiesel derived from soybean and higuereta oils. The thermogravimetric profiles of biodiesels indicated one step associated with decomposition of the ethyl esters. The thermal profiles were compared with reference diesel. Kinetic parameters were obtained by thermal analysis, estimating reaction order, pre-exponential factor and activation energy. For estimation of the kinetic parameters the Coats & Redfern method was used. The selection criterion of best fit was based on the correlation coefficient of the linear regression and the compensation effect of the kinetic parameters.     

Kinetics of thermal decomposition of aseptic packages

Kinetics of thermal decomposition of aseptic packages (e.g. Tetrapak cartons) and pyrolysis of this waste in a laboratory flow reactor was studied. Three different models for the calculation of the reaction rate and the determination of apparent kinetic parameters of thermal decomposition were used. The first method assumes a two stage thermal decomposition and the kinetic parameters were determined by fitting a derivative thermogravimetric (DTG) curve to experimentally determined thermogravimetric data of whole aseptic cartons. The second method uses kinetic parameters determined by fitting DTG curves to thermogravimetric data of individual components of aseptic packages. The last method was a multi-curve isoconversion method assuming a change of kinetic parameters with the increasing conversion. All types of the determined kinetic parameters were used in a mathematical model for thermal decomposition of mini briquettes made from aseptic packages at the temperature of 650°C. The model calculated also the heat conduction in the particles and it was verified by an independent set of experiments conducted in a laboratory screw type flow reactor.     

Toluene and benzyl decomposition mechanisms: elementary reactions and kinetic simulations

The high temperature decomposition kinetics of toluene and benzyl were investigated by combining a kinetic analysis with the ab initio/master equation study of new reaction channels. It was found that similarly to toluene, which decomposes to benzyl and phenyl losing atomic hydrogen and methyl, also benzyl decomposition proceeds through two channels with similar products. The first leads to the formation of fulvenallene and hydrogen and has already been investigated in detail in recent publications. In this work it is proposed that benzyl can decompose also through a second decomposition channel to form benzyne and methyl. The channel specific kinetic constants of benzyl decomposition were determined by integrating the RRKM/master equation over the C(7)H(7) potential energy surface. The energies of wells and saddle points were determined at the CCSD(T) level on B3LYP/6-31+G(d,p) structures. A kinetic mechanism was then formulated, which comprises the benzyl and toluene decomposition reactions together with a recently proposed fulvenallene decomposition mechanism, the decomposition kinetics of the fulvenallenyl radical, and some reactions describing the secondary chemistry originated by the decomposition products. The kinetic mechanism so obtained was used to simulate the production of H atoms measured in a wide pressure and temperature range using different experimental setups. The calculated and experimental data are in good agreement. Kinetic constants of the new reaction channels here examined are reported as a function of temperature at different pressures. The mechanism here proposed is not compatible with the assumption often used in literature kinetic mechanisms that benzyl decomposition can be effectively described through a lumped reaction whose products are the cyclopentadienyl radical and acetylene.     

Detailed experimental and kinetic modeling study of 3-carene pyrolysis

3-Carene is an important potential biofuel with properties similar to the jet-propellant JP-10. Its thermal decomposition and combustion behavior is to date unknown, which is essential to assess its quality as a fuel. A combined experimental and kinetic modeling study has been conducted to understand the initial decomposition of 3-carene. The pyrolysis of 3-carene was investigated in a jet-stirred quartz reactor at atmospheric pressure, at temperatures varying from 650 to 1050 K, covering the complete conversion range. The decomposition of 3-carene was observed to start around 800 K, and it is almost complete at 970 K. Online gas chromatography shows that primarily aromatics are generated which suggests that 3-carene is not a good fuel candidate. The potential energy surface for the initial decomposition pathways determined by KinBot shows that a hydrogen elimination reaction dominates, giving primarily cara-2,4-diene. Next to this molecular pathway, radical pathways lead to aromatics via ring opening. The kinetic model was automatically generated with Genesys and consists of 2565 species and 9331 reactions. New quantum chemical calculations at the CBS-QB3 level of theory were needed to calculate rate coefficients and thermodynamic properties relevant for the primary decomposition of 3-carene. Both the conversion of 3-carene and the yields of the primary products (ie, benzene and hydrogen gas) are well predicted with this kinetic model. Rate of production analyses shows that the dominant pathways to convert 3-carene are hydrogen elimination reaction and radical chemistry.     

A kinetic expression for the pyrolytic decomposition of polytetrafluoroethylene

Despite the fact that the thermal decomposition of polytetrafluoroethylene has been extensively studied over the past six decades, some inconsistencies regarding the kinetic parameters, e.g. the order of the reaction, remain. Representative kinetic data are essential for practical purposes such as reactor design and scaling. In general the literature data refer to homogeneous bulk heating, whereas the case of the non-homogeneous heating of a single particle has not received attention. Data (reaction rate and pre-exponential factor) applicable to this latter case were experimentally determined from isothermal thermogravimetric analyses of the depolymerisation reaction of PTFE. The kinetic data obtained on coarse granules (800-1000 μm) are reported here. The rate law is consistent with a shrinking particle kinetic model, with chemical kinetics controlling phase-boundary movement. The mass loss rate is directly proportional to surface area. A rate law applicable to this case, and useable for geometries of arbitrary shape, is derived.     

Studies on multistep thermal decomposition behavior of polytriazole polyethylene oxide‐tetrahydrofuran elastomer

Polytriazole polyethylene oxide‐tetrahydrofuran (PTPET) is an energetic propellant elastomer that is prepared using glycidyl azide polymer and trifunctional alkynyl‐terminated polyethylene oxide‐tetrahydrofuran. Its thermal decomposition, determined using thermogravimetic analysis, showed two mass‐loss peaks largely related to the decomposition of azide groups and the main chain. Flynn‐Wall‐Ozawa and Kissinger‐Akahira‐Sunose methods were deployed to obtain kinetic triplet parameters of PTPET thermal decomposition by the traditional model‐free method; the Coats‐Redfern approach was used as the model‐fitting method. Kinetics analysis indicated that the mechanism of the two‐step reactions were the primary‐reaction of first order and the power‐law phase reaction of the 2/3 order. The first decomposition stage of PTPET had an activation energy (E_a) of 113 to 116 kJ/mol while the second was 196 to 210 kJ/mol. The thermal decomposition of PTPET with different heating rates and mechanisms showed good kinetic compensation effects, the gas products being further studied with TG‐FTIR.     

Thermal behavior and decomposition kinetics of Formex-bonded explosives containing different cyclic nitramines

Thermal behavior and decomposition kinetics of Formex-bonded PBXs based on some attractive cyclic nitramines, such as 1,3,5-trinitro-1,3,5-triazinane (RDX) and 1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX). Actually, cis-1,3,4,6-tetranitrooctahy droimidazo-[4,5-d]imidazole (BCHMX) and 2,4,6,8,10,12-hexanitro-2,4,6,8,10, 12-hexaazaisowurtzitane (CL-20), was investigated by means of nonisothermal thermogravimetry (TG) and differential scanning calorimetry (DSC). It was found that the mass loss rate of PBXs involved in this research depends greatly on heating rate and the residue of the decomposition of these PBXs decreases with the heating rate. The onset of the exotherms was noticed at 215.4, 278.7, 231.2 and 233.7 °C with the peak maximum at 235.1, 279.0, 231.2 and 233.7 °C for RDX-Formex, HMX-Formex, CL-20-Formex, and BCHMX-Formex, respectively. Their corresponding exothermic changes were 1788, 1237, 691, and 1583 J g^?1. It was also observed that the dependence on the heating rate for onset temperatures of HMX- and BCHMX-based PBXs was almost the same due to their similar molecular structure. In addition, based on nonisothermal TG data, the kinetic parameters for thermal decomposition of these PBXs were calculated by isoconversional methods. It was shown that the Formex base has great effects on the activation energy distribution of nitramines. It was further found that the kinetic compensation effects occurred during the thermal decomposition of nitramine-based PBXs, and they almost have the same compensation effects due to similar decomposition mechanism.