= *G*.

(2b)

The initial temperature at the beginning of freeze is

(2c)

The maximum depth at any time to which the temperature disturbance will be felt is *X *+ δ. Then

(2d)

The energy balance at the phase change interface for the freeze process is

.

(2e)

The energy balance at the freezing front can also be written as two equations (Lunardini 1981b)

*T *(X , *t*)

2

2

+ *k*2

-*k*1 1

= ρ1lα1

(3)

2

*T *( X, *t*)

2

= ρ2lα 2

-*k*1

+ *k*2 2

(4)

.

Because of the initial temperature distribution, during freeze the heat flow to the interface from the

thawed region will exceed the geothermal heat flow until equilibrium is established. Likewise, during a

thaw period the heat flow from the thawed zone will be less than the deep geothermal heat flow.

An approximate solution to this problem will be obtained using the heat balance integral technique (see

Lunardini 1991). In this method, the differential equations are solved on average over a finite volume of

material rather than at each point of the region. The integration of the energy equations over the regions

where temperature changes are occurring, 0 ≤ *x *≤ *X *+ δ, detailed by Lunardini (1981b) is

∫

∫ T2 ( x, *t*)dx - ρ1l*X *+ (ρ2c2 - ρ1c1)Tf X

ρ1c1 T1( x, *t*)dx + ρ2c2

0

- ρ2c2 ( X + δ)*T*o + ( X + δ) = -*k*1 1

+ *k*2G .

(5)

2

Quadratic temperature profiles in regions 1 and 2 that satisfy the boundary conditions are chosen as

2

+ (a1X - ∆*T*1)

(6)

* X*

* X*

(x - *X*) 2

(

)

[

]

- *GX *+ ∆*T*

(7)

δ

δ2

where

8