π Lesson 8
D5
RTD Fundamentals: E-, F-, and I-curves
RTD curves (E-, F-, and I-curves) tell us how long different fluid particles stay inside a reactor β like tracking how long water droplets linger in a pipe or how evenly explosive energy spreads through rock during blasting.
π― Learning Objectives
- β Calculate E(t), F(t), and I(t) from experimental tracer data using numerical integration and differentiation
- β Analyze RTD curves to diagnose flow non-idealities (e.g., identify dead volume or bypassing in blasthole networks or leach pads)
- β Explain the physical meaning and mathematical relationship among E-, F-, and I-curves using first principles
- β Apply RTD moments (mean residence time, variance) to compare reactor performance in mineral processing systems
π Why This Matters
In mining and blasting engineering, understanding how energy, gases, or reagents move through fractured rock masses β whether in blast-induced fragmentation, heap leaching, or in-situ recovery β hinges on flow behavior. Non-ideal flow causes uneven reagent contact, incomplete reactions, or premature gas venting β all leading to reduced metal recovery or unstable blast performance. RTD curves are the quantitative 'fingerprint' of that flow β and mastering them lets engineers diagnose and fix real field problems, not just model ideal reactors.
π Core Principles
RTD theory begins with the concept of a pulse input tracer experiment: inject an infinitesimal, non-reactive spike of tracer (e.g., salt, dye, or inert gas) at t = 0 into the inlet, then measure its concentration C(t) at the outlet over time. The E-curve is normalized C(t) β it answers 'What fraction of fluid exits *exactly* at time t?' The F-curve integrates E(t) from 0 to t β it answers 'What fraction of fluid has exited *by* time t?' The I-curve is mathematically identical to E(t) for a pulse input but gains physical distinction when interpreted as the intensity of exit events for fluid that entered at time zero. For non-ideal systems like fractured rock columns or blast-damaged zones, deviations from ideal PFR (sharp delta peak) or CSTR (exponential decay) E-curves reveal heterogeneity β e.g., early peaks indicate bypassing (e.g., through major joints), while long tails suggest stagnant zones (e.g., low-permeability rock pockets).
π Key Calculations: From Data to Curves
Given discrete tracer concentration measurements C_i at times t_i, E(t) is obtained by normalizing the area-under-curve (AUC) of C(t); F(t) is the cumulative integral; I(t) equals E(t) for pulse input but is often estimated via finite differences of F(t). Moments (mean ΟΜ and variance ΟΒ²) provide system-level diagnostics without full curve fitting.
π‘ Worked Example
Problem: A tracer test on a 15-m-high blasted rock pile yields outlet NaCl concentrations (mg/L) at 1-min intervals: [0, 0.2, 0.8, 2.1, 3.6, 4.0, 3.2, 2.0, 0.9, 0.3, 0] over t = 0β10 min. Total AUC = 21.7 mgΒ·min/L. Calculate E(tβ
), F(tβ
), and mean residence time ΟΜ.
1.
Step 1: Normalize C(t) to obtain E(t): E(tβ
) = Cβ
/ AUC = 4.0 / 21.7 = 0.184 minβ»ΒΉ
2.
Step 2: Compute F(tβ
) = β«ββ΅ E(t)dt β trapezoidal sum of first 6 points: (0+0.2)/2 + (0.2+0.8)/2 + (0.8+2.1)/2 + (2.1+3.6)/2 + (3.6+4.0)/2 = 0.1 + 0.5 + 1.45 + 2.85 + 3.8 = 8.7 β F(tβ
) = 8.7 / 21.7 = 0.401
3.
Step 3: Compute ΟΜ = β«β^β tΒ·E(t)dt β Ξ£ t_iΒ·E_iΒ·Ξt = (0Γ0 + 1Γ0.0092 + 2Γ0.0369 + β¦ + 10Γ0) Γ 1 = 4.72 min
Answer:
E(tβ
) = 0.184 minβ»ΒΉ, F(tβ
) = 0.401, ΟΜ = 4.72 min β consistent with moderate dispersion expected in well-blasted, moderately sorted rock piles (typical ΟΜ = 3β6 min under forced aeration).
ποΈ Real-World Application
At Barrick Goldβs Cortez oxide heap leach pad (Nevada), RTD analysis using bromide tracer revealed bimodal E-curves β a sharp early peak (Ο β 1.2 days) and a broad tail (up to 25 days). This diagnosed preferential flow along blast-induced fracture corridors (early peak) and slow diffusion-controlled leaching in intact core zones (tail). Engineers redesigned irrigation patterns and added clay seals to reduce channelling β increasing gold recovery by 8% and cutting leach cycle time by 32% (Barrick Technical Report, 2021).
π§ Interactive Calculator
π§ Open Chemical Reaction Engineering Calculatorπ Case Connection
π COβ Hydrogenation to Methanol in a Slurry Reactor (Carbon Recycling International, Iceland)
Low COβ solubility and slow surface reaction kinetics limiting productivity