system, agitation will have an effect on reaction rate;
however, it is considered unlikely that this alone explains
the extreme difference in behaviour observed in this instance.
Unfortunately, we were unable to recreate the runaway
reaction in the laboratory. The intention, had the reaction
been reproducible, was to study it using the techniques of
reaction calorimetry, adiabatic calorimetry, and in situ FTIR
analysis. The experimental data generated would have
allowed us to differentiate between the two proposed
mechanisms. However, this did not prove possible.
Reaction calorimetry is used to study a desired synthetic
reaction. It gives data not only on the overall heat liberated,
or absorbed, by a reaction but also on the rate of heat
production or uptake. In this case the overall heat released
would have been the same if the reaction had proceeded by
either separate dehydrogenation/hydrogenation or combined
hydrogen-transfer, but the profile of the heat flow would have
differed markedly. The separate dehydrogenation/hydrogena-
tion reaction would have had distinct endothermic and
exothermic portions, whilst the combined hydrogen-transfer
reaction would have been exothermic throughout.
Table 2. Estimated parameters used for vapour pressure
prediction
parameter
estimated valuea
critical temperature
critical pressure
acentric factor
386 °C
30.8 bar
0.39937
a
All parameters were estimated in Predict.
a thermal runaway, assuming that no further reactions occur
at elevated temperatures. If a further reaction does occur,
for example the exothermic decomposition of a product, the
maximum temperature and pressure reached may be greater
than that predicted from the heat of reaction alone.
Prediction of the maximum temperature reached can be
simplified by treating the reaction mass as an adiabatic
system. This is a valid assumption for large production scale
reactors. In this case the extremely rapid temperature rise
also supports this assumption, as the reactor heating/cooling
system would not have had the capacity to remove sufficient
heat in the time scale required.
A value of 2000 J/kg/K has been used for the specific
heat capacity of the reaction mixture (estimated value). The
predicted adiabatic temperature rises for the two hypotheses
can then be compared. In the first hypothesis the dehydro-
genation reaction occurs during the heating phase so that
the runaway is caused purely by hydrogenation of the double
bond. In this case the predicted adiabatic temperature rise is
approximately 400 K. This figure can be amended to 360 K
because the reaction went to 90% completion. Similarly, the
predicted adiabatic temperature rise for the combined
hydrogen-transfer scenario is approximately 164 K (at 90%
conversion).
Adiabatic calorimetry is applied to the study of runaway
reactions under heat-loss conditions similar to those on plant,
that is essentially adiabatic. It yields data on the rates of
temperature and pressure increase and also the final tem-
perature and pressure attained by a runaway. Such data would
have helped differentiate between the two proposed mech-
anisms.
In situ FTIR can be used to monitor the progress of a
reaction. In this instance it would have been used to
determine whether any unsaturated ketone was present in
the reaction mixture. The presence of such material would
suggest that the dehydrogenation and hydrogenation reactions
had proceeded separately.
The maximum temperature of the reaction mass after the
runaway is thus approximately 505 °C for the separate
dehydrogenation/hydrogenation case or approximately 309
Theoretical Considerations
°
C for the combined hydrogen-transfer case.
The vapour pressure of the system (effectively TCD
The fact that we could not recreate the incident practically
led us to consider a theoretical approach. Various methods
were used to derive an estimate of the heats of reaction of
the dehydrogenation and hydrogenation steps. These included
the use of heats of formation of structurally related com-
Ketone A at this point in time) can then be estimated using
a prediction package such as PREDICT. This package offers
6
a variety of methods for estimating vapour pressure, includ-
ing the Pitzer, Riedel, and GomezThdos methods. These
methods were evaluated for TCD Ketone A, and the Pitzer
method was chosen as the most representative, although it
is interesting to note that there was relatively little variation
between the data generated by each model. It should be stated
here that these figures are only approximations and that their
derivation required the estimation of other parameters such
as the critical temperature and pressure (Table 2). Therefore
they should be treated with caution.
2
pounds published in the literature and also the computer-
3
based predictive package, CHETAH (this package takes the
approach of using Benson groups4 to estimate the thermo-
,5
dynamic properties of a molecule). The heats of reaction
-
1
derived were considered, and values of +60 kJ mol for
-1
the dehydrogenation and -120 kJ mol for the hydrogena-
tion were taken as being representative. This means that a
-
1
value of -60 kJ mol was used for the heat of reaction of
the combined hydrogen-transfer reaction.
These estimated heats of reaction can be used to predict
the maximum temperature the reaction mass will achieve in
The data generated by the Pitzer model is presented in
Table 1 for selected temperatures. It can be seen that in the
region of 310 °C, the maximum temperature attainable for
the combined hydrogen-transfer hypothesis, the vapour
pressure of the system is approximately 11.6 bar. It was only
possible to predict the vapour pressure at temperatures up
(2) Pedley, J. B. Thermochemical Data and Structures of Organic Compounds;
Thermodynamics Research Centre: Texas, 1994; Vol. 1.
(
3) CHETAH, ver 4.4, The ASTM Chemical Thermodynamic and Energy
Release EValuation Program; American Society for Testing and Materi-
als: West Conshohocken, PA, 1990.
(
(
4) Benson, S. W.; Buss, J. H. J. Chem. Phys. 1958, 29, 546.
5) Benson, S. W.; Cruickshank, F. R.; Golden, D. M.; Haugen, G. R.; O’Neal,
H. E.; Rodgers, A. S.; Shaw, R.; Walsh, R. Chem. ReV. 1969, 69, 279.
(6) PREDICT, Prediction of Thermodynamics and Transport Properties;
Dragon Technology Inc., P.O. Box 16012, Golden, CO 80402, U.S.A.
306
•
Vol. 6, No. 3, 2002 / Organic Process Research & Development