🎓 Lesson 20
D5
Co-Simulation of BEME Thermal Load, Ventilation Flow, and Electrical Grid Response
Co-simulation means running models of heat from battery-electric machines, airflow in mine tunnels, and power demand on the electrical grid all at once to see how they affect each other.
🎯 Learning Objectives
- ✓ Analyze thermal–ventilation–grid coupling using co-simulation output data
- ✓ Design BEME duty cycles that satisfy thermal limits, airflow requirements, and grid voltage drop constraints
- ✓ Calculate transient heat generation rates for battery packs under variable load profiles
- ✓ Explain how ventilation airflow rate modulates BEME surface temperature and subsequent grid reactive power demand
- ✓ Apply IEC 61850-10 and ISO 50001 principles to validate co-simulation boundary conditions
📖 Why This Matters
Underground mines deploying Battery-Electric Mobile Equipment (BEME) face a hidden systems challenge: heat from batteries raises tunnel air temperature, forcing ventilation fans to work harder—which draws more power, stressing an already constrained grid. If modeled separately, these effects are underestimated. Co-simulation reveals cascading failures before commissioning—like a 2023 incident at a Canadian nickel mine where uncoordinated BEME rollout caused 12% grid voltage sag, triggering ventilation shutdown and 4-hour production stoppage. This lesson equips you to prevent such failures by integrating physics across domains.
📘 Core Principles
Co-simulation rests on three pillars: (1) Domain decomposition—thermal (battery & motor loss), fluid (CFD-based airflow), and electrical (grid impedance + harmonic load models) are solved in separate but synchronized solvers; (2) Interface coupling—data exchange occurs at defined time steps via standardized middleware (e.g., FMI 3.0); and (3) Feedback fidelity—e.g., rising air temperature increases motor winding resistance → higher current draw → greater grid voltage drop → reduced inverter efficiency → more heat. Real-time thermal inertia (τ ~ 4–12 min for battery packs) and ventilation lag (τ ~ 30–90 s for main fans) dictate minimum co-simulation step size (typically 1–5 s). Stability requires convergence checks across domains at each step—failure here indicates design incompatibility, not solver error.
📐 Transient Thermal Load Coupling
This formula links BEME heat generation to ventilation airflow and grid voltage response. It quantifies how convective cooling efficiency affects downstream electrical demand.
Coupled Thermal-Ventilation-Grid Power Index (CTGPI)
CTGPI = (Q_waste / Q_conv) × (Z_grid × S_sc / P_rated)Dimensionless index quantifying integrated system stress: ratio of unremoved thermal load to convective cooling capacity, scaled by grid impedance impact.
Variables:
| Symbol | Name | Unit | Description |
|---|---|---|---|
| Q_waste | Waste heat generation | kW | Total thermal loss from battery, motor, and inverter during duty cycle |
| Q_conv | Convective heat removal | kW | Heat transferred from equipment surfaces to airflow via convection |
| Z_grid | Transformer impedance | pu (per unit) | Percent impedance rating of primary distribution transformer |
| S_sc | Grid short-circuit capacity | kVA | Available fault level at point of common coupling (PCC) |
| P_rated | Rated power of BEME unit | kW | Nameplate mechanical output power |
Typical Ranges:
Well-integrated BEME fleet: 1.2 – 2.8
Marginal ventilation + aging grid: 3.5 – 7.1
💡 Worked Example
Problem: A 35-ton BEME hauler operates at 85 kW average power over 10-min cycle. Battery efficiency = 92%, motor/inverter efficiency = 94%. Ventilation airflow = 45 m³/s at inlet temp = 22°C; tunnel cross-section = 28 m²; measured air temp rise = 3.7°C. Grid short-circuit MVA = 120; transformer impedance = 6.2%. Calculate CTGPI and interpret.
1.
Step 1: Compute total waste heat: Q_waste = P_in × (1 − η_battery × η_motor) = 85 kW × (1 − 0.92×0.94) = 85 × (1 − 0.8648) = 85 × 0.1352 ≈ 11.5 kW
2.
Step 2: Compute convective heat removal: Q_conv = ṁ × c_p × ΔT = (ρ × V̇) × c_p × ΔT = (1.2 kg/m³ × 45 m³/s) × 1005 J/kg·K × 3.7 K ≈ 54 kg/s × 1005 × 3.7 ≈ 200.3 kW
3.
Step 3: Compute CTGPI = (Q_waste / Q_conv) × (Z_grid × S_sc / P_rated) = (11.5 / 200.3) × (0.062 × 120,000 kVA / 85 kW) ≈ 0.0574 × (7440 / 85) ≈ 0.0574 × 87.53 ≈ 5.02
Answer:
CTGPI = 5.02 — exceeds the safe threshold of 3.0, indicating insufficient ventilation cooling relative to grid capacity. Remediation required: increase airflow ≥55 m³/s or reduce simultaneous BEME duty.
🏗️ Real-World Application
At the Kevitsa Mine (Finland), LKAB and ABB co-simulated 12 BEME units (loaders, trucks) with 3D CFD ventilation model (ANSYS Fluent) and real-time grid model (ETAP v22.1) over 24-hr shift cycles. Key findings: (1) Peak ventilation demand occurred 4.2 min after BEME startup—not at peak power—due to thermal inertia; (2) Grid voltage dip exceeded IEEE 1547-2018 limit (±5%) when >7 units started within 15 s; (3) Installing 200-kW battery-buffered UPS at substation reduced dip by 68%. The validated co-simulation informed revised dispatch logic and ventilation ramp-up protocols—now codified in LKAB’s BEME Integration Standard v3.1.
🔧 Interactive Calculator
🔧 Open Battery-Electric Mobile Equipment (BEME) Deployment Calculator📋 Case Connection
📋 Polish Hard Coal Mine BEME Ventilation Integration
Legacy ventilation system designed for diesel exhaust dilution; inability to reduce airflow without violating methane di...