# Reaction Calorimetry (RC1) in Process Safety & Scale-Up: Heat Flow, Heat Balance, Accumulation & MTSR Calculation

While Differential Scanning Calorimetry (DSC) screens raw material thermal stability, Reaction Calorimetry (e.g., METTLER TOLEDO RC1e / RC1 mx) measures the dynamic heat generated by the desired synthesis reaction under real-world recipe conditions (dosing profiles, temperature setpoints, agitation, mass transfer, and stoichiometry).

In semi-batch stirred tank reactors (STR)—the workhorse of Active Pharmaceutical Ingredient (API) and fine chemical plants—heat release is controlled primarily by reagent addition rates. Reaction calorimetry provides quantitative chemical engineering data to prevent thermal runaways, determine plant cooling jacket requirements, calculate unreacted reagent accumulation ($X_{accum}$), derive the Maximum Temperature of the Synthesis Reaction ($MTSR$), and assign the Stössel Criticality Class (Class 1 to 5).

# 1. Operating Principles: Heat Flow vs. Heat Balance Calorimetry

Reaction calorimeters operate inside automated glass or metallic benchtop reactors ($0.5 - 2.0\text{ L}$) equipped with precision temperature control loops, calibrated mass dosing pumps, torque sensors, and calibration heaters.

# A. Heat Flow Calorimetry (HFC)

In Heat Flow Calorimetry, the reactor temperature $T_r$ is held constant (isothermal mode) or ramped according to a set profile by dynamically adjusting the jacket fluid temperature $T_j$.

The heat flow through the reactor wall $q_{flow}$ is described by Fourier's law of heat conduction:

where:

• $U(t)$ is the overall heat transfer coefficient ($\text{W/(m}^2\cdot\text{K)}$).
• $A(t)$ is the wetted heat transfer area ($\text{m}^2$), which increases continuously during reagent dosing.
• $T_r - T_j$ is the temperature difference between the reaction mass and the cooling jacket.

# B. Heat Balance Calorimetry (HBC)

In Heat Balance Calorimetry, the jacket flow rate $\dot{m}_j$ is held constant at a high velocity, and the inlet jacket temperature $T_{j,in}$ and outlet jacket temperature $T_{j,out}$ are continuously monitored.

The heat removed by the jacket is calculated directly from the enthalpy gain of the heat transfer fluid (HTF):

This method is independent of heat transfer area $A(t)$ and wall fouling, making it ideal for crystallization, polymerization, or multiphase reaction masses where $U$ varies unpredictably.

# 2. Governing Equations & Mathematical Data Integration

The primary objective of a reaction calorimetry run is solving the instantaneous reactor energy balance equation to isolate the chemical heat generation rate $q_{rxn}(t)$ ($\text{W}$ or $\text{W/kg}$).

# A. Comprehensive Energy Balance Equation

The dynamic energy balance around a batch or semi-batch reactor is expressed as:

where:

• $q_{flow} = U A (T_r - T_j)$ is the jacket heat exchange rate.
• $q_{accum} = \left( m_r C_{p,r} + m_{vessel} C_{p,vessel} \right) \frac{dT_r}{dt}$ is the thermal accumulation due to reactor mass temperature change.
• $q_{dos} = \dot{m}_{feed} \cdot C_{p,feed} \cdot (T_r - T_{dos})$ is the sensible heat required to heat incoming feed from dosing temperature $T_{dos}$ to reactor temperature $T_r$.
• $q_{reflux} = \dot{m}_{evap} \cdot \Delta H_{vap}$ is the heat removed by reflux condensation.
• $q_{loss}$ is ambient heat loss through reactor head space (calibrated prior to run).
• $q_{ins} = \omega \cdot \tau$ is mechanical agitator power dissipation (usually negligible except in high-viscosity masses).

# B. In-Situ Calibration of Thermal Transmittance ($U \cdot A$)

Because fluid viscosity, mass density, and liquid level change during dosing, $U \cdot A$ must be evaluated before, during, and after the reaction using an internal electrical calibration heater:

where $q_{cal}$ is a known precise electrical power input (e.g., $15.0\text{ W}$).

+---------------------------------------------------------------------------------------+
|                    STEP-BY-STEP RC1 MATHEMATICAL INTEGRATION                          |
+---------------------------------------------------------------------------------------+
| 1. Baseline Calibration: Fire qcal pulse --> Compute initial (U·A)1                   |
| 2. Reagent Dosing: Log Tr(t), Tj(t), m_feed(t) --> Compute qflow(t) and qdos(t)       |
| 3. Isolate qrxn(t): qrxn(t) = U·A(t)·(Tr - Tj) + (m·Cp)·(dTr/dt) + qdos              |
| 4. Total Enthalpy: ΔHtotal = ∫[t0 to tend] qrxn(t) dt                                  |
| 5. Specific Enthalpy: ΔHrxn = ΔHtotal / n_limiting (kJ/mol)                           |
+---------------------------------------------------------------------------------------+

# C. Unreacted Reagent Accumulation ($X_{accum}$)

In semi-batch processes, if the chemical reaction rate $r_{rxn}$ is slower than the feed dosing rate $\dot{m}_{feed}$, unreacted reagent accumulates in the vessel.

The instantaneous fractional conversion $Y(t)$ and unreacted accumulation $X_{accum}(t)$ are calculated by integrating the heat flow rate $q_{rxn}(t)$:

If a reaction is purely dosing-controlled (instantaneous conversion), $X_{accum}(t) \approx 0$. If a reaction is kinetically controlled, $X_{accum}(t) \to 1.0$ ($100\%$ accumulation).

# D. Calculation of $MTSR$ (Maximum Temperature of Synthesis Reaction)

If cooling is completely lost at the moment of maximum reagent accumulation $X_{accum,max}$, the accumulated unreacted feed will react adiabatically, driving the reactor temperature up to the Maximum Temperature of the Synthesis Reaction ($MTSR$):

where:

• $T_p$ is the desired process operating temperature ($^\circ\text{C}$).
• $\Delta T_{ad,rxn} = \frac{-\Delta H_{total}}{m_{batch} \cdot C_{p,batch}}$ is the adiabatic temperature rise of the total synthesis reaction.

# 3. Stössel Criticality Classification Hierarchy

Developed by Dr. Francis Stoessel, this framework assigns process safety risk based on the thermodynamic hierarchy of four key temperatures:

1. $T_p$: Process Operating Temperature.
2. $MTSR$: Maximum Temperature of the Synthesis Reaction.
3. $T_D$ (or $T_{D,24}$): Decomposition Onset Temperature (where Time-to-Maximum-Rate $\text{TMR}_{ad} = 24\text{ hours}$, obtained from ARC testing).
4. $T_{max}$: Maximum Technical Temperature (boiling point $T_{bp}$ or emergency relief setpoint $T_{relief}$).

# Detailed Breakdown of Stössel Criticality Classes

# 4. Sample Analysis Illustration & Step-by-Step Result Derivation

To illustrate how RC1 experimental data evaluates synthesis safety and assigns Stössel criticality, examine a real-world case study: Semi-batch Nitration of an Aromatic Intermediate ($1.0\text{ L}$ RC1 reactor).

# A. Experimental Setup & Recipe Parameters

• Batch Mass ($m_{batch}$): $1.25\text{ kg}$ ($1250\text{ g}$).
• Reagent Dosing Mass ($m_{feed}$): $250\text{ g}$ Nitric acid solution fed over $90\text{ minutes}$.
• Process Setpoint ($T_p$): $35.0\text{ }^\circ\text{C}$.
• Specific Heat Capacity ($C_{p,batch}$): $1.95\text{ J/(g}\cdot\text{K)}$.
• Atmospheric Boiling Point ($T_{max}$): $110.0\text{ }^\circ\text{C}$.
• ARC Secondary Decomposition Onset ($T_{D,24}$): $95.0\text{ }^\circ\text{C}$.

+---------------------------------------------------------------------------------------+
|                       RC1 EXPERIMENTAL INTEGRATION WORKFLOW                           |
+---------------------------------------------------------------------------------------+
|  Time Period (min)   |  Dosing Rate  | Avg qflow (W) |  Accumulated Heat ∫ qrxn dt   |
+----------------------+---------------+---------------+-------------------------------+
|  t = 0 to 15 min     |   2.78 g/min  |    12.4 W     |         11.16 kJ              |
|  t = 15 to 90 min    |   2.78 g/min  |    38.5 W     |        173.25 kJ              |
|  t = 90 to 120 min   |   0.00 g/min  |    17.0 W     |         30.60 kJ (Post-rxn)   |
+---------------------------------------------------------------------------------------+
                                           |
                                           v
+---------------------------------------------------------------------------------------+
| 1. Total Reaction Enthalpy:  ΔHtotal = -215.01 kJ                                     |
| 2. Total Adiabatic Rise:     ΔTad,rxn = 215,010 J / (1250 g × 1.95 J/g·K) = 88.2 K   |
| 3. Post-Dosing Heat Area:    Unreacted Heat = 30.60 kJ                                |
| 4. Max Accumulation Ratio:   Xaccum = 30.60 kJ / 215.01 kJ = 0.142 (14.2% Accumulation) |
+---------------------------------------------------------------------------------------+

# B. Step-by-Step MTSR Calculation & Criticality Decision

# 5. Scale-Up Application & Heat Transfer Calculations

Reaction calorimetry data provides the exact parameters required for scaling up reactions from bench scale ($1\text{ L}$) to production reactors ($10,000\text{ L}$).

# A. Industrial Cooling Jacket Duty Check

The maximum required cooling duty $Q_{cool,max}$ ($\text{kW}$) for a commercial vessel of volume $V$ is given by:

where $S.F. \approx 1.25$ is a $25\%$ engineering safety factor.

The minimum required heat transfer area $A_{min}$ for the plant reactor is checked against available jacket area:

# 6. Where to Use Reaction Calorimetry in API Scale-Up

1. Process Optimization (Kilolab Stage): Optimize dosing time $t_{feed}$ to minimize cycle time while keeping heat generation within plant cooling limits ($q_{rxn} < U A \Delta T_{max}$).
2. Thermal Safety Assessment: Calculate $MTSR$ and establish Stössel Criticality Class to select SIL safety instrumented systems.
3. Green Chemistry & Yield Engineering: Monitor heat release kinetics to detect side-reactions, slow mass transfer limitations, or catalyst decay.