🎓 Lesson 11
D5
CrunchFlow for Multi-Redox Systems: Modeling Mn, Fe, and S Cycling
CrunchFlow is a computer tool that simulates how water, chemicals, and minerals interact underground—especially how manganese, iron, and sulfur change forms and move through mine waste over time.
🎯 Learning Objectives
- ✓ Explain the role of redox coupling in controlling Mn, Fe, and S speciation in mine drainage
- ✓ Design a CrunchFlow input file for a multi-redox system including kinetic mineral reactions and electron acceptor/donor constraints
- ✓ Analyze simulated concentration breakthrough curves to diagnose redox zonation and reaction rate limitations
- ✓ Calculate electron balance residuals to verify redox consistency in CrunchFlow output
- ✓ Apply CrunchFlow results to predict long-term acid rock drainage (ARD) or neutral mine drainage (NMD) behavior in sulfidic tailings
📖 Why This Matters
Mine waste often contains pyrite, sphalerite, and Mn-oxides—minerals that react over decades to generate acidic, metal-laden drainage. But real systems don’t behave as isolated redox pairs: Fe oxidation consumes O₂, enabling sulfate reduction; Mn oxide precipitation scavenges dissolved Fe²⁺; and H₂S can reduce MnO₂, restarting cycles. Ignoring these couplings leads to wildly inaccurate predictions of metal release and treatment longevity. CrunchFlow is one of only a few tools capable of resolving these feedbacks—making it essential for robust closure planning and regulatory compliance.
📘 Core Principles
Reactive transport in multi-redox systems requires simultaneous solution of mass conservation (advection + dispersion), chemical equilibrium (log K constraints), and kinetic rate laws (e.g., surface-controlled oxidation). Redox couples are linked via shared electrons: e.g., oxidation of Fe²⁺ produces electrons used to reduce SO₄²⁻ → HS⁻. CrunchFlow tracks total electron equivalents (‘e⁻-equivalents’) and enforces stoichiometric coupling (e.g., 8e⁻ per SO₄²⁻ reduced; 1e⁻ per Fe²⁺ oxidized). It uses an implicit finite-difference solver with automatic time stepping and handles mineral surface area evolution, pH buffering, and secondary mineral saturation limits. Crucially, it distinguishes between *thermodynamically favored* and *kinetically accessible* reactions—e.g., pyrite oxidation may dominate over slower MnO₂ reduction even if the latter is more thermodynamically favorable.
📐 Electron Balance Residual
The electron balance residual quantifies numerical consistency in redox-coupled simulations. A large residual (>10⁻⁶ mol e⁻/L) indicates incomplete convergence or missing redox pathways. CrunchFlow reports this automatically—but users must validate it to trust predictions.
💡 Worked Example
Problem: In a 1D column simulation (1 m length, 0.25 cm/s flow), CrunchFlow outputs at t = 1 yr: ∑oxidation e⁻ produced = 1.37 × 10⁻³ mol e⁻/L; ∑reduction e⁻ consumed = 1.3698 × 10⁻³ mol e⁻/L.
1.
Step 1: Compute difference: EBR = |1.37e−3 − 1.3698e−3| = 2.0 × 10⁻⁷ mol e⁻/L
2.
Step 2: Compare to tolerance threshold (1 × 10⁻⁶ mol e⁻/L)
3.
Step 3: Since EBR < tolerance, redox accounting is consistent; model output is numerically reliable.
Answer:
The result is 2.0 × 10⁻⁷ mol e⁻/L, which falls within the safe range of <1 × 10⁻⁶ mol e⁻/L.
🏗️ Real-World Application
At the Duck Pond Mine (Newfoundland, Canada), CrunchFlow modeled 30-year evolution of a Mn-rich tailings impoundment where initial porewater contained Fe²⁺, Mn²⁺, and SO₄²⁻ under mildly reducing conditions. The model included kinetic oxidation of pyrite (rate law from Williamson & Rimstidt, 1994), microbial sulfate reduction (Monod kinetics), and MnO₂(s) reductive dissolution by HS⁻. Simulations predicted a transient Mn-oxide 'scavenging front' migrating downward at ~0.15 m/yr—verified by sequential extraction profiles from core samples collected at 5- and 10-year intervals. This validated the model’s ability to capture coupled Fe-Mn-S redox buffering, directly informing the design of a shallow oxygen barrier for long-term Mn immobilization.
🔧 Interactive Calculator
🔧 Open Mine Waste Characterization & Geochemical Modeling Calculator📋 Case Connection
📋 Copper Mine Waste Rock Stockpile ARD Mitigation at Escondida Extension
High-pyrite waste rock (>3.2% S) stockpiled without cover; predicted ARD onset within 5 years
📋 Gold Tailings Geochemical Stabilization at Granny Smith Mine (WA)
Arsenic-rich tailings (up to 120 mg/kg As) exhibiting elevated As leaching under oxidizing conditions
📋 Limestone Mine Neutral Drainage Management at Mount Read Complex (Tasmania)
Historic waste dumps containing carbonate-hosted Pb-Zn mineralization generating neutral metal leachate (Zn >15 mg/L, Cd...
📋 Iron Ore Mine Waste Rock Long-Term Stability at Brockman 4 (Pilbara)
Massive hematite-goethite waste rock (low sulfide but high Mn/Al) showing delayed acidity and Al leaching post-construct...
📋 Coal Mine Spoil Geochemical Capping at Hunter Valley Reclamation Project
Spoil with pyritic shale interbeds generating ARD despite initial alkaline overburden; inconsistent capping led to local...