🎓 Lesson 17
D5
Void Fraction Prediction Using Drift-Flux Model
Void fraction is the percentage of space in a pipe or blast-induced fracture network that is filled with gas (like air or explosive gases) instead of liquid or solid material.
🎯 Learning Objectives
- ✓ Calculate void fraction using the drift-flux model given phase velocities, densities, and flow rates
- ✓ Analyze the impact of drift velocity and distribution parameter on void fraction prediction accuracy
- ✓ Explain the physical significance of the distribution parameter C₀ and drift velocity V₉ in multiphase flow through fractured rock or blast-generated void networks
- ✓ Apply the drift-flux correlation to estimate gas holdup in post-blast venting scenarios during underground mine de-gassing
📖 Why This Matters
In underground mining and post-blast ventilation, predicting how much gas occupies fractures and cavities—i.e., the void fraction—is essential for safe degassing, preventing methane accumulation, and optimizing dewatering systems. Underestimating void fraction can lead to inadequate ventilation design, risking explosive atmospheres; overestimating it may cause oversized, costly infrastructure. The drift-flux model provides a physically grounded, yet computationally efficient way to predict this key parameter where traditional homogeneous models fail—especially in inclined, low-velocity, or density-stratified flows common in blast-damaged rock masses.
📘 Core Principles
The drift-flux model treats two-phase flow as a mixture with a weighted average velocity (j = jₗ + j_g), then introduces corrections for local phase slip. It defines void fraction α via the fundamental relation: j_g = C₀·j·α + V₉·α(1−α), where C₀ is the distribution parameter (reflecting flow geometry and turbulence) and V₉ is the drift velocity (representing buoyancy-driven relative motion). For vertical upward flow in fractured media, C₀ ≈ 1.1–1.2 and V₉ ≈ 0.3–0.8 m/s are empirically validated. In blast contexts, the effective 'channel' is not a smooth pipe but a stochastic fracture network—so C₀ must be adjusted for tortuosity and aperture distribution, often using digital rock physics calibrations.
📐 Drift-Flux Void Fraction Formula
The drift-flux model rearranges to solve explicitly for void fraction α when total superficial velocity j and gas superficial velocity j_g are known. This quadratic form enables direct calculation without iteration under steady-state assumptions.
💡 Worked Example
Problem: Given: total superficial velocity j = 0.45 m/s, gas superficial velocity j_g = 0.12 m/s, distribution parameter C₀ = 1.15, drift velocity V₉ = 0.48 m/s — calculate void fraction α in a blast-fractured limestone horizon.
1.
Step 1: Write the drift-flux equation: j_g = C₀·j·α + V₉·α(1−α)
2.
Step 2: Rearrange into standard quadratic: V₉·α² + (C₀·j − V₉)·α − j_g = 0 → 0.48α² + (1.15×0.45 − 0.48)α − 0.12 = 0 → 0.48α² + 0.0375α − 0.12 = 0
3.
Step 3: Solve using quadratic formula: α = [−0.0375 + √(0.0375² + 4×0.48×0.12)] / (2×0.48) = [−0.0375 + √0.232] / 0.96 ≈ (−0.0375 + 0.482)/0.96 ≈ 0.444/0.96 ≈ 0.463
Answer:
The void fraction is α ≈ 0.463 (or 46.3%), which falls within the typical range of 0.35–0.55 for moderately fractured, upward-venting blast zones in carbonate rock.
🏗️ Real-World Application
At the Stillwater Platinum Mine (Montana, USA), engineers used the drift-flux model to size post-blast ventilation ducts in stopes with high methane potential. Field measurements of gas flow rate (2.1 m³/min) and water cut (68% liquid by volume) in a 300-mm-diameter fractured raise were combined with core-derived fracture aperture distributions to calibrate C₀ = 1.18 and V₉ = 0.52 m/s. Predicted α = 0.41 matched neutron porosity log estimates within ±3%, enabling accurate fan selection and reducing degassing time by 22% versus homogeneous-model designs.
✏️ Student Exercise
A student is analyzing venting from a recently blasted stope in porphyritic granite. Measured total volumetric flow is 0.85 m³/s in a 450-mm-diameter drift. Gas flow rate is 0.31 m³/s. Assume C₀ = 1.12 (based on image analysis of blast fracture maps) and V₉ = 0.39 m/s (from lab-scale air-water flow in replica fracture networks). Calculate the void fraction α. Then, determine whether the result suggests stratified (α < 0.25), slug (0.25 ≤ α ≤ 0.55), or annular (α > 0.55) flow regime per the Hewitt–Robertson classification.