🎓 Lesson 1 D1

Getting Started with Reaction Engineering and Kinetics

Reaction engineering and kinetics is the science of understanding how fast chemical reactions happen—and how to control them safely and efficiently in real-world systems like reactors or explosive detonations.

🎯 Learning Objectives

  • Explain the physical meaning of the rate law and distinguish between zero-, first-, and second-order reactions
  • Calculate reaction rate constants from experimental concentration–time data using linearized plots or numerical integration
  • Analyze and compare ideal reactor performance (CSTR vs. PFR) for a given reaction system using design equations
  • Apply the Arrhenius equation to predict how temperature changes affect reaction rate in mining-related processes (e.g., ANFO decomposition, acid leaching)
  • Design a simple batch reactor volume to achieve 90% conversion for a first-order exothermic reaction, accounting for thermal constraints

📖 Why This Matters

In mining and blasting engineering, reaction kinetics governs critical processes—from the detonation velocity and energy release of explosives like ANFO and emulsions, to heap leaching kinetics of gold ores with cyanide, and even the stability of reactive mine waste (e.g., acid rock drainage). Misjudging reaction rates can lead to incomplete blasting, hazardous gas buildup, inefficient metal recovery, or environmental failures. Mastering kinetics isn’t just theory—it’s the difference between predictable fragmentation and flyrock, or profitable extraction and costly reprocessing.

📘 Core Principles

Reaction kinetics begins with the rate law: rate = k [A]^m [B]^n, where k is the temperature-dependent rate constant and m,n are reaction orders determined experimentally—not from stoichiometry. The rate-determining step often controls overall behavior, especially in heterogeneous systems like ore leaching (solid–liquid) or explosive decomposition (solid–gas). We then integrate rate laws into reactor design equations: for a batch reactor, dN_A/dt = –r_A V; for a CSTR, V = F_{A0}(X)/–r_A; for a PFR, dX/dV = –r_A / F_{A0}. Coupling these with energy balances becomes essential when reactions are exothermic—common in blasting chemistry and sulfide oxidation.

📐 Arrhenius Equation

The Arrhenius equation quantifies how reaction rate constants change with temperature—a vital tool for predicting explosive sensitivity, leach tank residence time, or storage stability of reactive reagents.

Arrhenius Equation

k = A exp(–E_a / RT)

Relates rate constant k to absolute temperature T via activation energy E_a and pre-exponential factor A.

Variables:
SymbolNameUnitDescription
k Rate constant varies (e.g., s⁻¹, L·mol⁻¹·s⁻¹) Temperature-dependent kinetic parameter
A Pre-exponential factor same as k Frequency factor reflecting collision/orientation probability
E_a Activation energy J·mol⁻¹ Minimum energy barrier for reaction
R Universal gas constant J·mol⁻¹·K⁻¹ 8.314 J·mol⁻¹·K⁻¹
T Absolute temperature K Thermodynamic temperature of the system
Typical Ranges:
Explosive decomposition (e.g., ANFO): 120–180 kJ/mol
Acid leaching of oxides: 40–70 kJ/mol
Bio-oxidation of sulfides: 30–55 kJ/mol

💡 Worked Example

Problem: An ANFO decomposition reaction has a pre-exponential factor A = 1.2 × 10^13 s⁻¹ and activation energy E_a = 145 kJ/mol. Calculate k at 25°C (298 K) and 60°C (333 K). Use R = 8.314 J/mol·K.
1. Step 1: Convert E_a to J/mol → 145,000 J/mol
2. Step 2: Compute exponent term at 298 K: –E_a/(R·T) = –145000/(8.314 × 298) ≈ –58.57
3. Step 3: Calculate k_298 = A × exp(–58.57) ≈ 1.2×10¹³ × 1.8×10⁻²⁶ ≈ 2.2×10⁻¹³ s⁻¹
4. Step 4: Repeat for 333 K: –145000/(8.314 × 333) ≈ –52.43 → k_333 = 1.2×10¹³ × 1.5×10⁻²³ ≈ 1.8×10⁻¹⁰ s⁻¹
5. Step 5: Compare: k increases ~820× between 25°C and 60°C—highlighting strong thermal sensitivity relevant to summer storage safety.
Answer: k₂₉₈ = 2.2 × 10⁻¹³ s⁻¹; k₃₃₃ = 1.8 × 10⁻¹⁰ s⁻¹ — a >800-fold increase confirms why ANFO must be stored below 40°C to prevent autocatalytic decomposition.

🏗️ Real-World Application

At the Grasberg copper-gold mine (Indonesia), heap leaching kinetics of secondary sulfide ores were modeled using shrinking-core kinetics to optimize irrigation rate and acid addition. Field-measured Cu²⁺ concentration vs. time data revealed an apparent first-order dependence on acid concentration and half-order on oxygen partial pressure. By fitting the Arrhenius parameters (E_a = 62 kJ/mol, A = 3.8 × 10⁷ h⁻¹), engineers predicted leach cycle duration across seasonal temperature swings—reducing under-leaching by 22% and avoiding premature pad abandonment.

📋 Case Connection

📋 Bioethanol Fermentation Bioreactor Scale-Up with Inhibition Kinetics

Ethanol inhibition caused premature cessation at large scale despite matching nominal conditions

📋 Nitric Acid Absorption Tower Design for Tail-Gas Treatment

Incomplete absorption of NO and NO₂ due to slow liquid-phase oxidation kinetics and poor gas distribution

📚 References