Stability criteria and critical runway conditions of propylene glycol manufacture in a continuous stirred tank reactor

Criterios de estabilidad y condiciones críticas de runaway de la producción de propilenglicol en un reactor continuo de tanque agitado

Here, a new method for the analysis of the steady state and the safety operational conditions of the hydrolysis of propylene oxide with excess of water, in a Continuous Stirred Tank Reactor (CSTR), was developed. For industrial operational typical values, at first, the generated and removed heat balances were examined. Next, the effect of coolant fluid temperature in the critical ignition and extinction temperatures (T_CI and T_CE, respectively) was analyzed. The influence of the heat exchange parameter (hS) on coolant and critical temperatures was also studied. Finally, the steady state operation areas were defined. The existence of multiple stable states was recognized when the heat exchange parameter was in the range 6.636 < hS kJ/(min.K) < 11.125. Unstable operation area was located between the T_CI and T_CE values, restricting the reactor operation area to the low stable temperatures.

En este trabajo se presenta un nuevo método para el análisis del estado estable y de las condiciones operativas seguras para la hidrólisis del óxido de propileno con exceso de agua en un reactor continuo de tanque agitado (CSTR). Inicialmente, para valores típicos de operación industrial, se examinaron los balances de calor generado y removido. Luego, se analizó el efecto de la temperatura del fluido de enfriamiento en las temperaturas de ignición y extinción (TCI y TCE, respectivamente). También se estudió el efecto del parámetro de transferencia de calor (hS) sobre las temperaturas críticas y del fluido de enfriamiento. Finalmente, se definieron las áreas de operación estables. Se reconoció la existencia de múltiples estados estables cuando el parámetro de transferencia de calor está en el intervalo de 6.636 < hS kJ/(min.K) < 11.125. El área de operación inestable está localizada entre los valores de TCI y TCE, restringiendo la operación del reactor al área de bajas temperaturas estables.

DOI: https://doi.org/10.15446/ing.investig.v35n2.46258

S. M. López-Zamora¹, I. Dobrosz-Gómez², and M. Á. Gómez-García³

---

¹ Sandra Milena López Zamora. Chemical Engineer, Universidad Nacional de Colombia, Manizales, Colombia. M.Sc. in Chemical Engineering. Affiliation: Universidad Nacional de Colombia, Manizales, Colombia. Grupo de Investigación en Procesos Reactivos Intensificados con Separación y Materiales Avanzados (PRISMA). E-mail: smlopezz@unal.edu.co.

² Isabel Dobrosz Gómez. Chemist, Lodz University of Technology, Lodz, Poland. M. Sc. Eng., Lodz University of Technology, Lodz, Poland. PhD., Lodz University of Technology, Lodz, Poland. Affiliation: Grupo de Investigación en Procesos Reactivos Intensificados con Separación y Materiales Avanzados (PRISMA). E-mail: idobrosz-gomez@unal.edu.co.

³ Miguel Ángel Gómez. Chemical Engineer, Universidad Nacional de Colombia, Manizales,Colombia. M.Sc. in Chemical Engineering, Universidad Nacional de Colombia, Bogotá, Colombia. PhD in Chemical Engineering, Université Louis Pasteur, Strasbourg I, France. Affiliation: Grupo de Investigación en Procesos Reactivos Intensificados con Separación y Materiales Avanzados (PRISMA).
E-mail: magomez@unal.edu.co.

How to cite: López-Zamora, S. M., Dobrosz-Gómez, I., & Gómez-García, M. Á. (2015). Stability criteria and critical runway conditions of propylene glycol manufacture in a continuous stirred tank reactor. Ingeniería e Investigación, 35(2), 56-60.

DOI: https://doi.org/10.15446/ing.investig.v35n2.46258

---

ABSTRACT

The present paper presents the development of a new method for the analysis of the steady state and the safety operational conditions of the propylene oxide hydrolysis with excess of water in a Continuous Stirred Tank Reactor (CSTR). At first, for the industrial operational typical values, the generated and removed heat balances were examined. Next, the effect of coolant fluid temperature on the critical ignition and extinction temperatures (TCI and TCE, respectively) was analyzed. The influence of the heat exchange parameter (hS) on coolant and critical temperatures was also studied. Finally, the steady state operation areas were defined. The existence of multiple stable states was recognized when the heat exchange parameter was in the range of 6.636 < hS kJ/(min.K) < 11.125. The Unstable operation area was located between the TCI and TCE values, restricting the reactor operation area to the low stable temperatures.

Keywords: Stability criteria, propylene oxide, propylene glycol, runaway conditions.

---

RESUMEN

Palabras clave: Criterios de estabilidad, óxido de propileno, propilen glicol, condiciones de runaway.

---

Received: October 15th 2014

Accepted: March 4th 2015

---

Introduction

Propylene glycol is an organic compound obtained by the hydrolysis of propylene oxide with excess of water. It is a colorless and odorless liquid. Nowadays it is produced on a large scale, and primarily used in the production of polymers (saturated polyester resins) and in the food processing (E-number E1520). Propylene glycol presents anti-freeze properties and it is used as fluid for heat transfer at low temperatures; as solvent in paints and cleaners; and as an intermediate for the synthesis of different chemical products such as alkyds (The Dow Chemical Company 2011).

The hydrolysis of propylene oxide with water is an exothermic reaction. Under certain conditions, it presents multiple steady states and/or runaway behavior, making the reactive operation potentially risky. For safety, it is very important to determine how stable the steady state is. Different methodologies for the safety operation of reactors and analysis of runaway conditions have been proposed. As a first approximation, thermal runaway prediction was performed by a lineal stability and overall energy balances analysis. It was established that critical conditions depend on the temperature curve. Thus, cooling in a chemical reactor can be controlled by thermal conduction in the reactive mixture (Semenov 1928; Kamenetskii 1969). More recently, the concept of parametric sensitivity analysis was proposed by Bilous & Amundson (1956). They suggested that when the system operates in the parametric sensible area, slight changes in the inlet parameters of the reactor generate excessive changes in the outlet parameters. Thus, different criteria based on the geometric conditions of temperature and conversion curves have been developed. For the Continuous Stirred Tank Reactor (CSTR) systems some examples can be pointed out: (i) Luo et al. (2000) studied the polymerization of phenol-formaldehyde reaction and (ii) Ojeda et al. (2014) studied the cyclo-trimethylene-triamine production process.

In the case of propylene oxide hydrolysis, Švandová et al. (2005) discussed the relation between software tools, primarily designed for mathematical modeling and simulation, and chemical reactors analysis. They identified the Hazard and Operability (HAZOP) zone for the hydrolysis reaction. Additionally, Vadapalli & Seader (2001) employed ASPEN Plus software for the study of multiple steady-states of this reaction, and Molnár et al. (2003) presented the sensitivity analysis of a CSTR with cooling jacket.

In this work, a new method for the analysis of steady state and safety operational conditions of the hydrolysis of propylene oxide with excess of water is proposed. Thus, at first, the generated and removed heat balances are analyzed. Next, the effect of coolant fluid in the critical ignition and extinction temperatures is presented. The influence of the heat exchange parameter (hS) on coolant and critical temperatures is also studied. Finally, the steady state operation areas are defined.

Mathematical model

For the propylene oxide hydrolysis, Fogler (2008) reports a first order kinetic (n=1) with respect to propylene oxide and zero order one with respect to water Equations (1-2).

(1)

(2)

The overall energy balance in a CSTR implies that the heat generation rate (QG), Equation (3), minus the heat removal rate (QR), Equation (4), equals the heat accumulation (Equation (5)).

(3)

(4)

(5)

If the reactive system operates adiabatically, conversion and reactive temperature, for a time t, can be correlated by the ratio between the transient heat release (ΔΗt) and the overall heat release (ΔΗtotal), as shown in Equation (6). Thus, the overall balance can be rewritten as follows (Equation (7)):

(6)

(7)

Uncontrolled reaction conditions (runaway) can occur when the heat generated in the process is higher than the removed one. In this case, the exothermic system goes through a quick increase in temperature (runaway condition). Semenov (1928) defined the sufficient, necessary, and critical runaway conditions as given by Equations (8) and (9).

(8)

(9)

Thus, considering steady state, from Equations (7) and (9), the critical temperature can be defined (Equation (10)), while the transition points, Equations (11)-(13), can be obtained by equaling the two roots of Equation (10). In fact, it presents two possible numerical solutions: the low one, which corresponds to the critical ignition temperature and the high one, which is the critical extinction temperature.

(10)

(11)

(12)

(13)

Thus, the analysis of the reactive system stability was performed based on the definitions presented above. Kinetic parameters and typical industrial operation conditions are given in Table 1.

Table 1. Kinetic parameters and operation conditions (Molnár et al. 2003).

Pre-exponential factor (k∞)	16.19 x 1012 1/h
Activation energy (Eα)	32400 Btu/mol
Propylen e oxide feed	19.523 kmol/h
Methanol feed	32.6 kmol/h
Water feed (containing 0.1%wt H2SO4)	364.15 kmol/h
Feed temperature	297 K
Maximum Temperature (Tmax)	348.23 K
Reactive volume	1.133 m3
Coolant feed temperature	288 K
Coolant rate	35 kmol/h

Results and discussion

The desirable operation conditions, defined from the transition point, the critical ignition and extinction temperatures, and the generated and removed heat curves are presented in Figure 1. The Figure also includes the overall energy balance expressed as the removed heat (straight lines) and generated heat (QG), at a coolant temperature of 288.15 K. As a reference, the feed temperature of the reactive mixture (Tfeed) and the vaporization temperature of propylene oxide (Tvap = 324.75 K) were indicated. It is important to notice that operation at temperatures higher than the vaporization one should be avoided.

Figure 1. Temperature effect on the generated and removed heat.

Table 2. Results for tangent and intersection points:

	Tc	TS	hS (kJ/(min.K))
Extinción	TCE = 331.39	TSE =291.65	11.125
Ignición	TCI =300.79	TSI =343.25	6.636

Two tangent and two intersection points (TCE, TCI, TSE, TSI) can be defined from Figure 1. They correspond to the critical hS slope values (Table 2). Thus, unstable operation conditions will correspond to the values 6.636 < hS kJ/(min.K) < 11.125. If the reactive system operates at these conditions, the existence of multiple steady states can be expected. In this region, three intersection points can be obtained: TS,L, TS,M and TS,H (Figure 1).

For an intermediate hS value of 8.8807, the intersection points were calculated (dotted line). These values are included in Table 3 and compared with the results reported by Molnár et al. (2003), obtained using a CSTR transient model. From the analysis of the TS,M value, one can see that a slight change in temperature can take the system into a higher (runaway) or lower operation temperature. Even if the higher operation state implies higher productivity, it is also related to a non-safety operation. Hence, operational values proposed by Molnár et al. (2003) can be considered as stable ones, taking into consideration that: firstly, operation temperature is lower than the ignition critical point (TS,L<TCI) and secondly it is near the intersection ignition point (TS,H~TSI).

Table 3. Results for the three intersection points of the unstable area

	hS (kJ/(min.K))	TS,L (K)	TS,M (K)	TS,H (K)
Unstable operation	8.8807	293.15	316.74	340.31
Molnár et al., 2003	6.9574	297.05	---	342.90

Figure 2 presents the effect of coolant temperature on the intersection and tangent points (TCI, TCE, TSI, TSE). Slight rise in coolant temperature results in an increase in TSE, TSL and TCI and in a decrease in TCE, TSH, and TSI. All lines converge to the transition point (Figure 2). Reactive system should operate in the stable state area without surpassing vaporization temperature.

Figure 2. Coolant temperature effect on Tc ad Ts values.

Thus, different operation areas can be defined:

(i) A stable high temperature area (TSH), limited by the final ignition temperature stable point (TSI) and extinction critical temperature (TCE),

(ii) A stable low temperature area (TSL), restricted by the final stable extinction temperature point and the critical ignition temperature (TCI);

(iii) An unstable area between critical extinction (TCE) and ignition temperatures (TCI); and

(iv) A stable area external to the final extinction (TSE) and ignition temperature points (TSI).

In order to avoid non-safety or non-controlled conditions, the reactive system should operate in the (iv) area.

Figure 3 displays the effect of coolant temperature on the slope of removed heat (hS). Maximum hS value corresponds to the transition point (hStr=16.74 kJ/(min.K)). This is a watershed of temperature: it increases as propylene oxide fraction increases
(see inserted figure in Figure 3). The value of hS increases as the Ta increases. When TC>TC,tr, the value of TC is TCE. If TC<TC,tr, TC is equal to TCI. When TS>TC,tr, the value of TS is TSI. When TS<TC,tr, TS is equal to TSE. Unstable operation area is located between the critical and stable ignition and extinction temperatures. High hS values are recommended in order to assure operation at stable area. Thus, for an analogous parameter, UA, Fogler (2008) (UA = 126.6 kJ/(min.K)) and Molnár et al., (2003)(UA = 240 kJ/(min.K)) guaranteed in their analysis that the removed heat will be higher than the generated one. However, no explanation related to the selection of these values was included. From the results of this study, it is possible to argue that in order to work in the stable low temperature area (TSL), the hS values should be higher than 9.38 kJ/(min.K).

Figure 3. Relation between hS parameter and coolant temperature (Ta) for the propylene oxide hydrolysis (propylene oxide
fraction = 0.0469).

Conclusions

A new method for the analysis of a steady state and safety operational conditions of the hydrolysis of propylene oxide to propylene glycol in a CSTR was developed. Thus, the existence of multiple stable states was recognized when the parameter hS, related to the reactor heat exchange, was in the range of 6.636 < hS kJ/(min.K) < 11.125.

The effect of feed temperature on the critical and stable ignition and extinction temperature locus was also considered and the transition point was calculated. This point separates the stable working area from the unstable one and allows defining different operational conditions. Unstable operation area was located between TCI and TCE, indicating that the operating point should be in the external area and under the vaporization temperature of propylene oxide. It restricts the reactor operation area to the stable low temperatures.

Nomenclature

C Concentration (mol/l)

Cp Heat capacity (J/(g.K))

Ea Activation energy (kJ/mol)

hS Heat exchange parameter (kJ/(min.K))

ΔΗt Transient heat release (J)

ΔΗtotal Overall released heat (J)

k Reaction rate constant (M1-n/min)

k∞ Pre-exponential factor

m Mass (g)

n Reaction order

Q Heat (kJ/min)

q Molar exothermic heat (kJ/mol)

R Gas universal constant (J/(gmol.K))

rA Reaction rate (mol/(l.min))

T Temperature (K)

t Reaction time (min)

UA Analogous to hS heat exchange parameter

V Volume of reactant (m3)

X Fractional conversion

Subscripts

A Species A

a Ambient

C Critical

CE Critical extinction

CI Critical ignition

feed Feed

G Generation

H High

L Low

M Medium

max Maximum

R Removal

S Stable

SE Stable extinction

SI Stable ignition

tr Transition

vap Vaporization

0 Initial

Greek letter

ρ Density (g/cm3)

References

Bilous, O. & Amundson, N. R. (1956). Chemical Reactor Stability. AIChE Journal, 2, pp.117–126.

DOI: 10.1002/aic.690020124.

Fogler, H.S. (2008). Elementos de ingeniería de las reacciones químicas (4ta Ed.) Pearson Education.

Kamenetskii, F. (1969). Diffusion and Heat Transfer in Chemical Kinetic (2da Ed.). P. Press, ed., New York.

Luo, K. M. et al. (2000). The critical runaway condition and stability criterion in the phenol–formaldehyde reaction. Journal of Loss Prevention in the Process Industries, 13(2), pp.91–108. DOI: 10.1016/S0950-4230(99)00071-6.

Molnár, A., Markos, J. & Jelemensky, L. (2003). Safety analysis of CSTR towards changes in operating conditions. Journal of Loss Prevention in the Process Industries, 16, pp.373–380. DOI: 10.1016/S0950-4230(03)00069-X.

Ojeda, J. C. et al. (2014). Análisis de la Sensibilidad Paramétrica del Proceso de Producción de Ciclo-Trimetileno-Triamina. Información tecnológica, 25(4), pp.153–162.

DOI: 10.4067/S0718-07642014000400017.

Semenov, N. N. (1928). Zur Theorie des Verbrennungsprozesses. Zeitschrift für Physik, 48(7-8), pp.571–582.

DOI: 10.1007/BF01340021.

Švandová, Z. et al. (2005). Steady States Analysis and Dynamic Simulation as a Complement in the Hazop Study of Chemical Reactors. Process Safety and Environmental Protection, 83(5), pp.463–471.

DOI: 10.1205/psep.04262.

The Dow Chemical Company. (2011). Hoja de datos técnicos,

Vadapalli, A. & Seader, J. D. (2001). A generalized framework for computing bifurcation diagrams using process simulation programs. Computers & Chemical Engineering, 25, pp.445–464. DOI: 10.1016/S0098-1354(01)00624-X.