A general criterion to define runaway limits in chemical reactors

A general runaway criterion valid for single as well as for multiple reaction types, i.e. consecutive, parallel, equilibrium, and mixed kinetics reactions, and for several types of reactors, i.e. batch reactor (BR), semibatch reactor (SBR) and continuous stirred tank reactor (CSTR) has been developed. Furthermore, different types of operating conditions, i.e. isoperibolic and isothermal (control system), have been analysed. The criterion says that we are in a runaway situation when the divergence of the system becomes positive (div>0) on a segment of the reaction path. The results show that this is a general runaway criterion than can be used to calculate the runaway limits for chemical reactors. The runaway limits have been compared with previous criteria. A considerable advantage, over existing criteria, is that it can be calculated on-line using only temperature measurements and, hence, it constitutes the core of an early warning runaway detection system we are developing.     

On the divergence criterion for runaway detection: Application to complex controlled systems

It is well-known that, for certain values of the operative parameters influencing the dynamic behavior of a chemical reactor, a phenomenon known as thermal runaway (that is, a loss of the reactor temperature control) may arise. Such a situation can be really dangerous because above a certain threshold temperature value unwanted side reactions or, worse, decompositions of the reacting mixture may be triggered evolving high amounts of flammable or toxic gases that can cause reactor pressurization and, eventually, its explosion. For this reason, since the beginning of the previous century a number of studies concerning the prediction of the so called runaway boundaries has been carried out. In this work, a modified version of the divergence criterion for runaway detection, originally developed by Zaldívar and co-workers, is presented. Such a modified divergence criterion is capable of treating whatever type of complex controlled reacting system (taking into account not only temperature control but also dosing strategies) and its reliability has been demonstrated for isoperibolic semibatch reactors using literature experimental data concerning the nitration of 4-Chlorobenzotrifluoride in mixed acids and the nitric acid oxidation of 2-octanol to 2-octanone and further carboxylic acids.     

Thermal risk in batch reactors: Theoretical framework for runaway and accident

Thermal safety and risk of accidents are still challenging topics in the case of batch reactors carrying exothermic reactions. In the present paper, the authors develop an integrated framework focusing on defining the governing parameters for the thermal runaway and evaluating the subsequent risk of accident. A relevant set of criteria are identified in order to find the prior conditions for a thermal runaway: failure of the cooling system, critical temperature threshold, successive derivatives of the temperature (first and second namely) vs. time and no detection in due time (reaction time) of the runaway initiation. For illustrative purposes, the synthesis of peracetic acid (PAA) with hydrogen peroxide (HP) and acetic acid (AA) is considered as case study. The critical and threshold values for the runaway accident are identified for selected sets of input data. Under the conditional probability of prior cooling system failure, Monte Carlo simulations are performed in order to estimate the risk of thermal runaway accident in batch reactors. It becomes then possible to predict the ratio of reactors, within an industrial plant, potentially subject to thermal runaway accident.     

Probe into effect of the initiator on the thermal runaway hazards of methyl methacrylate bulk polymerization

The bulk polymerization of methyl methacrylate (MMA) is of great importance in chemical industry, but the polymerization process is highly hazardous, and few reports have focused on the effect of initiators on its thermal hazards. In this work, to thoroughly explore the thermal hazard characteristics, the runaway behavior of MMA bulk polymerization is investigated by a combination of thermodynamics experimental and kinetics theoretical methods. The results indicate that the presence of initiator exhibits an undesirable thermal hazard to the MMA bulk polymerization, and its exothermic behavior is also greatly influenced by the type and concentration of initiator. For azobisisoheptanenitrile (ABVN), azodiisobutyronitrile (AIBN) and dibenzoyl peroxide (BPO) initiators as examples, the AIBN-initiated reaction has the shortest adiabatic induction period (39.51 min), whereas the BPO-initiated polymerization exhibits the strongest maximum temperature-rising rate and maximum pressure-rising rate. Under adiabatic runaway, the temperature and pressure change significantly with increasing AIBN concentration, revealing a great potential risk of thermal runaway. Kinetic parameters are calculated to further understand the thermal runaway mechanisms, showing a strong agreement with the adiabatic experimental data. Finally, based on the cooling failure scenario, severity grading is determined by the evaluation criteria. The current work provides extensive data as a reference and guidance for the process design and optimization of MMA bulk polymerization from the perspective of safety.     

Gas-phase thermal explosions in catalytic direct oxidation of alkenes

The metal-based catalytic oxidation of alkenes to the corresponding epoxides is playing a significant role in the modern chemical industry. Nevertheless, these key processes are still lacking proper understanding with respect to the gas-phase runaway behaviour (thermal explosion) and to the hot spot formation on the catalytic surface, under the typical process conditions.This work aims to enlighten these aspects by considering either the catalytic or the gas-phase chemistry for the development of reactor operative diagrams, in order to define the best-operating conditions with respect to the selectivity, the productivity, and the process safety aspects.The proposed methodology has been applied to the oxidation of ethylene and propylene for the direct oxidation process by pure oxygen, considering a detailed kinetic model accounting for the homogeneous reactions, coupled with the heterogeneous catalytic mechanisms.Sensitivity and reaction path analyses were performed to individuate the ruling species and reactions determining the transition to runaway conditions.     

The critical runaway condition and stability criterion in the phenol–formaldehyde reaction

The procedure of phenol–formaldehyde polymerization is a rather important and complicated reaction in the chemical industry. This exothermic polymeric reaction releases a huge amount of heat. The high amount of energy accumulated and increasing temperature in this reaction process always lead to runaway reaction and a hazard situation owing to the high released heat and improper operation. In this investigation, we used sodium hydroxide as alkali–catalyst in the phenol–formaldehyde polymerization and estimated the reaction kinetics parameters to evaluate the thermal hazard conditions. The critical temperatures and stable criteria of the runaway reaction in this exothermic polymerization were evaluated. This technique is important and useful for safe operation in the phenol–formaldehyde polymerization process.     

Dynamic modeling and bifurcation analysis for the methyl isocyanate hydrolysis reaction

The aim of this work was to characterize the methyl isocyanate hydrolysis reaction and to identify its operational criteria. The parametric sensitivity and dynamic stability methodologies were performed at the Bhopal disaster circumstances, over the relevant operating parameter space. Stable and unstable conditions, bifurcations points, turning points and oscillatory behavior were determined. The combined methodology give useful guidance on the operative conditions selection and the appropriate strategy to overcome hazardous situations. The obtained results demonstrated high sensitivity to small perturbations (thermal runaway) and prevalent oscillatory behavior. Moreover, the following critical parameters for the studied dynamic system were defined: the inverse residence time of 1.5700103 and the heat transfer coefficient of 752.394.     

Thermal risk in batch reactors: Case of peracetic acid synthesis

The present paper deals with accidents risk in batch reactors. It identifies the conditions for the occurrence of a thermal runaway and develops a probabilistic approach to assess the relevant risk. It investigates also the conditions for optimal synthesis of peracetic acid (PAA) with hydrogen peroxide (HP) and acetic acid (AA). The kinetic model of reversible reaction and side reaction of PAA synthesis is used to predict reactor temperature and molar ratio of PAA by ASPEN PLUS software. A sensitivity analysis is performed under different conditions such as constant temperature or adiabatic process with different concentrations of sulfuric acid. Assuming a prior cooling system failure, the conditions for reaction runaway triggering a thermal accident are identified in the case of PAA synthesis. Monte Carlo simulations are used in order to calculate the conditional probability of accident and optimize the synthesis of PAA.     

Inherent thermal runaway hazard evaluation method of chemical process based on fire and explosion index

Thermal runaway hazard assessment provides the basis for comparing the hazard levels of different chemical processes. To make an overall evaluation, hazard of materials and reactions should be considered. However, most existing methods didn't take the both into account simultaneously, which may lead the assessment to a deviation from the actual hazard. Therefore, an integrated approach called Inherent Thermal-runaway Hazard Index (ITHI) was developed in this paper. Similar to Dow Fire and Explosion Index(F&EI) function, thermal runaway hazard of chemical process in ITHI was the product of material factor (MF) and risk index (RI) of reaction. MF was an indicator of material thermal hazards, which can be determined by initial reaction temperature and maximum power density. RI, which was the product of probability and severity, indicated the risk of thermal runaway during the reaction stage. Time to maximum rate under adiabatic conditions and criticality classes of scenario were used to indicate the runaway probability of the chemical process. Adiabatic temperature rise and heat of the desired reaction and secondary reaction were used to determine the severity of runaway reaction. Finally, predefined hazard classification criteria was used to classify and interpret the results obtained by this method. Moreover, the method was validated by case studies.     

化学放热系统热爆炸潜在危险性评价

讨论化学放热系统的热稳定性和临界条件,用化学反应物无消耗的假设推导化学放热系统热失控(热爆炸)时的动力学参数临界值,得到热失控的判据、临界点火温度和熄火温度。提出用系统安全指数概念来评价放热反应系统发生热爆炸的潜在危险性,分析化学放热系统的平衡域。用硝酸甲酯分解爆炸实例,说明如何利用安全指数对具有热爆炸可能性的系统的潜在危险性进行定量评价,其预测结果与实验结果一致。     

Influence of calorimetric approximations on the design of emergency relief systems

The design of an emergency relief system (that is, a pressure safety valve or a rupture disk) for vessels, which may involve runaway reactions, requires knowledge of the chemical kinetics of the reactions involved. When safety-related problems are considered this is usually achieved using calorimetric tests, coupled with some suitable approximations on the kinetics of the reacting system. In this work we have analysed the extent to which the precise knowledge of the chemical kinetics influences the size of the relief system device for different reaction conditions. Decision criteria are proposed to identify situations where approximations in the kinetic mechanism lead to underestimation in the venting area.     

某恒温间歇反应的热失控研究

为得到某恒温间歇反应冷却系统临界温度Tc并评估该反应体系的热危险性,基于间歇反应体系热参数的敏感性,探讨数值计算方法和选用4种反应热失控临界判据得到的不同Tc值。结果表明,数值计算方法得到了反应体系的Tc为31.7℃,Semenov判据、Sliding判据、Da/St判据得到的Tc值分别为10℃,17℃和27℃;无量纲绝热温升B给出的结论是,在允许的工艺温度范围内,该反应对任意的冷却系统温度都处于难以控制状态。数值计算方法及各判据得到的Tc的较大差异性说明各方法均有一定的局限性。在运用参数敏感性分析的基础上,结合风险矩阵图方法对系统冷却失效的热危险性进行评估,得到不同Tc下风险可接受、有条件接受及不可接受对应的参数范围的物料累积度Xac与最大反应速率到达时间θ间的关系。     

Diphenylmethane diisocyanate self-polymerization: Thermal hazard evaluation and proof of runaway reaction in gram scale

Thermal analysis by differential scanning calorimetry and thermogravimetric/differential thermal analysis mass spectrometry, adiabatic calorimetry, a gram-scale heating test, and infrared spectroscopy were performed to evaluate the thermal hazards of diphenylmethane diisocyanate (MDI) and prove the occurrence of a runaway reaction. The self-polymerization of MDI was found to occur at about 340 °C under rapid heating conditions. Carbon dioxide was eliminated and heat was generated to allow polymerization. Under adiabatic and closed conditions, the runaway reaction of MDI can begin at least from 220 °C. Besides it is highly probable that the runaway reaction of MDI can begin from a lower temperature in an actual process scale. More heat was generated than in the previous case and the pressure rose rapidly. A closed 2-mm-thick glass vessel exploded because of the runaway reaction of MDI even if the temperature was lower than 300 °C. Therefore, MDI could cause fatal runaway reactions below 300 °C, where MDI had been assumed to self-polymerize by eliminating carbon dioxide previously.     

Using thermal analysis with kinetic calculation method to assess the thermal stability of azo compound on construction materials

The application of construction polymers in engineering and alternative materials has always occupied a place in the market. In the production process of polymer resins, initiators can be used to lower the polymerization reaction energy threshold, which can improve reaction efficiency and reduce energy loss. However, as a commonly used energetic substance in the polymerization process, azos have caused related process hazards due to their exothermic characteristics. Because of this, it is essential to examine and analyze the thermal hazard characteristics of emerging azo substances, such as 2-cyanopropan-2-imemicarbazide (CABN). Although previous literature performs the calculation on related thermal hazard parameters of CABN, there is still exists a void for discussion in estimating the reaction model to avoid analogous hazards and enhance the existing thermal analysis. Based on the past literature, the reaction model is improved with thermogravimetric analysis as evidence. The revised thermal hazard parameters are calculated as the basis of control and mitigation measures, the kinetic model is used to estimate the modified safety parameters, and in the judgment of the runaway reaction, the critical temperature of the runaway is found by analyzing the influence of slight changes in ambient temperature on the reaction temperature. The results show that the critical temperature that causes CABN to enter the runaway reaction is delayed, and the hazard is lower than in the storage situation. Therefore, the thermal hazard to CABN mainly focuses on the safety environment and measures during storage.     

Parametric sensitivity analysis for thermal runaway in semi-batch reactors: Application to cyclohexanone peroxide reactions

The semi-batch reactors (SBRs) system, which is widely used in industrial processes, possesses an intrinsic parametric sensitivity, in which infinitesimal disturbances of input parameters can result in large variations in output variables. In this work, local parametric sensitivity analysis (PSA) was used to understand parameter variations and global PSA was conducted to examine the interaction of input parameters. The effects of these parameters on the output of the system model were analyzed based on the Monte Carlo method with Latin hypercube sampling and the extended Fourier amplitude sensitivity test model. The results showed that the evolution of thermal behaviors in SBRs were observed: marginal ignition; thermal runaway; and the quick onset, fair conversion, and smooth temperature profile. The threshold point of transition from marginal ignition to thermal runaway was at the maximal value of local sensitivity, for which the slope with respect to cooling temperature equaled zero. Moreover, the sequence of the global sensitivity of six common input parameters was computed and evaluated. The reliability of the numerical models was verified by using our previous experimental results of cyclohexanone peroxide reaction. This comprehensive sensitivity analysis could provide valuable operating information to improve chemical process safety.     

Using thermal analysis with kinetic calculation method to assess the thermal stability of 2-cyanopropan-2-yliminourea

Azo compounds, which are commonly used as initiators and blowing agents, are also typical self-reactive materials capable of undergoing runaway reaction during storage and transportation, which can cause severe fires and accidents. To ensure the thermal safety of azo compounds in the process, transportation, and applications, this study investigated 2-cyanopropan-2-yliminourea, which can also be called V-30. First, thermal decomposition characteristics under the non-isothermal conditions were obtained using differential scanning calorimetry. Second, the collected data were combined with a mathematical model to evaluate the primary thermal hazard during the process for V-30. Then, based on a heat-transfer model, the self-accelerating decomposition temperature (SADT) was extrapolated for consideration and non-consideration of consumption of chemicals. The results showed that SADT of V-30 was less than 80 °C. Therefore, it is essential to avoid a temperature beyond SADT or the cooling system will fail. The influence of consumption was also considered for SADT in this study.     

The thermal hazard evaluation of 1,1-di (tert-butylperoxy) cyclohexane by DSC using non-isothermal and isothermal-kinetic simulations

1,1-Di (tert-butylperoxy) cyclohexane (DTBPH) has been widely employed in the chemical industry. Unfortunately, organic peroxides have been involved in many serious fires and explosions in manufacturing processes, storage, and transportation. This study investigated the thermokinetic parameters by isothermal kinetic and non-isothermal-kinetic simulation, using differential scanning calorimetry (DSC) tests. DSC was applied to assess the kinetic parameters, such as kinetic model, frequency factor (ln k₀), activation energy (E_a), reaction order, and heat of reaction (ΔH_d). Comparisons of non-isothermal and isothermal-kinetic model simulation led to a beneficial kinetic model of thermal decomposition to predict the thermal hazard of DTBPH. Simulations of a 0.5 L Dewar vessel and 25 kg barrel commercial package in liquid thermal explosion models were performed and compared to the results in the literature. From the results, the optimal conditions for use of DTBPH to avoid violent runaway reactions during the storage and transportation were determined. This study established the features of thermal decomposition that could be executed as a reduction of energy potential and storage conditions in view of loss prevention.     

基于VSP2的引发剂质量分数对醋酸乙烯聚合反应失控特性及安全泄放设计影响研究

为了解决醋酸乙烯聚合反应失控所引起的超压问题,通过VSP2绝热量热仪研究了醋酸乙烯聚合反应的失控特性,并通过Leung's法对某醋酸乙烯聚合反应器的安全泄放面积进行了计算;然后,在其他条件不变的情况下,研究引发剂质量分数对失控特性和泄放面积的影响,结果表明,引发剂质量分数对反应总放热量的影响不大,体系绝热温升为105~115℃;但引发剂质量分数越大,失控反应的最大温升速率和最大压升速率越大。这是因为引发剂质量分数越大,在相同泄放压力和最大累积压力下,单位质量反应物的放热速率就越大,也就需要更大的泄放面积;最后,引入无量纲数W~*、G~*和A~*,拟合出它们与引发剂质量分数X*的关系式,结果表明,在研究范围内所需安全泄放面积随引发剂质量分数线性增大。     

Causes and consequences of thermal runaway incidents—Will they ever be avoided?

A study of runaway incidents involving thermal chemical reactions in the UK over the past 25 years (1988–2013) has been carried out. The objective of this study is to determine possible causes of thermal runaway incidents. A statistical analysis of the underlying problems that led to thermal runaway incidents has been provided. A comparison of the current study on thermal runaway incidents with those identified prior to 1988 has been carried out. This study clearly shows that lessons have not been learnt from thermal runaway incidents caused by operator errors, management failures and lack of organised operating procedures. These factors have been the possible causes of about 77% of all the thermal runaway incidents analysed in this study. The number of fatalities and injuries as a result of thermal runaway incidents has increased by ∼325% and ∼279%, respectively, in the last 25 years even though the number of incidents was significantly less. On the basis of this analysis, several recommendations have been proposed that could help to minimise the risks associated with any thermal runaway incidents in the future.