Geothermal Convection in Porous Media

Download as pdf or txt
Download as pdf or txt
You are on page 1of 2

1

Notes on
1.63 Advanced Environmental Fluid Mechanics
Instructor: C. C. Mei, 2002
[email protected], 1 617 253 2994

November 6, 2002

6.4 Geothermal convection in porous media


Motivations: Geothermal Energy, Pipe line in permafrost, thermal stress around buried
nuclear waste.
S. Chandrasekhar: Hydrodynamic and Hydromagnetic Stabililty, Oxford, 1961.
D. A. Nield and A. Bejan, Convection in porous media , Springer-Verlag, 1992.
Let us begin with the effective equations over macroscales.
Recall the definitions of seepage velocity,
1 ZZZ
u= q dV (6.4.1)
Ω Ωf

and of porosity:
Ωf
n= (6.4.2)

Mass conservation for an incompressible fluid:

∇·u=0 (6.4.3)

Momentum conservation (Darcy ’s law)


µ
0 = −∇p − u + ρf g (6.4.4)
k
Inview of the small size of the solid grains heat diffuses by conduction almost instantly in
time scales of geothermal interest. We therefore assume the temperature in both phases to
be equal,
Ts = Tf = T (6.4.5)
The energy conservation law in each phase is :
Solid
∂T
(1 − n)(ρC)s = (1 − n)Ks ∇2 T (6.4.6)
∂t
fluid:
∂T
n(ρC)f + (ρC)f u · ∇T = nKf ∇2 T (6.4.7)
∂t
2

Adding the two we get the energy equation for the mixture
∂T
(ρC)m + (ρC)f u · ∇T = Km ∇2 T (6.4.8)
∂t
where
(ρC)m = (1 − n)(ρC)s + n(ρC)f , Km = (1 − n)Ks + nKf (6.4.9)
Defining
(ρC)m
σ= (6.4.10)
(ρC)f
and
Km
κm = (6.4.11)
(ρC)f
as the mixture diffusivity, we can rewrite (6.4.8) as

∂T
σ + u · ∇T = κm ∇2 T (6.4.12)
∂t
Note: Strictly speaking, the above equation ought to be derived from the microscale
considerations, i.e., by homogenization analysis.
Equation of state:
ρf = ρo [1 − β(T − To )] (6.4.13)
In summary the governing eqautions are (6.4.3), (6.4.4), (6.4.12) and 6.4.13).

You might also like