Explicit and Approximated Solutions For Heat and Mass Transfer Problems With A Moving Interface

20
Explicit and Approximated Solutions

for Heat and Mass Transfer
Problems with a Moving Interface
Domingo Alberto Tarzia
CONICET and Universidad Austral
Argentina
1. Introduction
The goal of this chapter is firstly to give a survey of some explicit and approximated
solutions for heat and mass transfer problems in which a free or moving interface is
involved. Secondly, we show simultaneously some new recent problems for heat and mass
transfer, in which a free or moving interface is also involved. We will consider the following
problems:
1. Phase-change process (Lam-Clapeyron-Stefan problem) for a semi-infinite material:
i. The Lam-Clapeyron solution for the one-phase solidification problem (modeling the
solidification of the Earth with a square root law of time);
ii. The pseudo-steady-state approximation for the one-phase problem;
iii. The heat balance integral method (Goodman method) and the approximate solution for
the one-phase problem;
iv. The Stefan solution for the planar phase-change surface moving with constant speed;
v. The Solomon-Wilson-Alexiades model for the phase-change process with a mushy
region and its similarity solution for the one-phase case;
vi. The Cho-Sunderland solution for the one-phase problem with temperature-dependent
thermal conductivity;
vii. The Neumann solution for the two-phase problem for prescribed surface temperature at
the fixed face;
viii. The Neumann-type solution for the two-phase problem for a particular prescribed heat
flux at the fixed face, and the necessary and sufficient condition to have an
instantaneous phase-change process;
ix. The Neumann-type solution for the two-phase problem for a particular prescribed
convective condition (Newton law) at the fixed face, and the necessary and sufficient
condition to have an instantaneous phase-change process;
x. The similarity solution for the two-phase Lam-Clapeyron-Stefan problem with a
mushy region.
xi. The similarity solution for the phase-change problem by considering a density jump;
xii. The determination of one or two unknown thermal coefficients through an over-
specified condition at the fixed face for one or two-phase cases.
xiii. A similarity solution for the thawing in a saturated porous medium by considering a
density jump and the influence of the pressure on the melting temperature.
www.intechopen.com
Advanced Topics in Mass Transfer

440
2. Free boundary problems for the diffusion equation:
i. The oxygen diffusion-consumption problem and its relationship with the phase-change
problem;
ii. The Rubinstein solution for the binary alloy solidification problem;
iii. The Zeldovich-Kompaneets-Barenblatt solution for the gas flow through a porous
medium;
iv. Luikov coupled heat and mass transfer for a phase-change process;
v. A mixed saturated-unsaturated flow problem representing absorption of water by a soil
with a constant pond depth at the surface and an explicit solution for a particular
diffusivity;
vi. Estimation of the diffusion coefficient in a gas-solid system;
vii. The coupled heat and mass transfer during the freezing of the high-water content
materials with two free boundaries: the freezing and sublimation fronts.
2. Explicit solutions for phase-change process (Lam-Clapeyron-Stefan
problem) for a semi-infinite material
Heat transfer problems with a phase-change such as melting and freezing have been studied
in the last century due to their wide scientific and technological applications. A review of a
long bibliography on moving and free boundary problems for phase-change materials
(PCM) for the heat equation is shown in (Tarzia, 2000a). Some previous reviews on explicit
or approximated solutions were presented in (Garguichevich & Sanziel, 1984; Howison,
1988; Tarzia, 1991b & 1993). Some reviews, books or booklets in the subject are (Alexiades &
Solomon, 1993; Bankoff, 1964; Brillouin, 1930; Cannon, 1984; Carslaw & Jaeger, 1959; Crank,
1984; Duvaut, 1976; Elliott & Ockendon, 1982; Fasano, 1987 & 2005; Friedman, 1964; Gupta,
2003; Hill, 1987; Luikov, 1968; Lunardini, 1981 & 1991; Muehlbauer & Sunderland, 1965;
Primicerio, 1981; Rubinstein, 1971; Tarzia, 1984b & 2000b; Tayler, 1986).
2.1 The Lam-Clapeyron solution for the one-phase solidification problem (modeling
the solidification of the Earth with a square root law of time)
We consider the solidification of semi-infinite material, represented by x 0 > . We will find the
interface solid-liquid x s t ( ) = and the temperature T T x t ( , ) = of the solid phase defined by
( )
( ) ( )
( )
f
T x t if x s t t
T x t
T if s t x t
, 0 , 0
,
, 0
< < >
=

>
(1)
which satisfy the following free boundary problem:
( )
t xx
cT kT x s t t 0 , 0 , 0 = < < > (2)
( )
f
T t T T t
0
0, , 0 = < > (3)
( ) ( )
f
T s t t T t , , 0 = > (4)
( ) ( ) ( )
x
kT s t t s t t , , 0 = > (5)
www.intechopen.com
Explicit and Approximated Solutions for Heat and Mass Transfer

441
s(0) 0 = (6)
Eq. (2) represents the heat equation for the solid phase, k is the thermal conductivity, is
the mass density, c is the heat capacity, is the latent heat of fusion by unit of mass, T
0
is
the imposed temperature at the fixed face x 0 = , and the material is initially at the melting
temperature
f
T . The problem (2)-(6) is known in literature as the one-phase Stefan problem
(Lam-Clapeyron-Stefan problem) and the condition (5) as the Stefan condition. Free
boundary problems of this type were presented by the first time in (Lam & Clapeyron,
1831) in order to study the solidification of the Earth and was continued independently by
(Stefan, 1891a, b & 1990) in order to study the thickness of polar ice. We remark here that
Lam & Clapeyron found the important law for the phase-change interface with a square
root of time.
Theorem 1. (Lam-Clapeyron solution).
The explicit solution to the free boundary problem (2)-(6) is given by

f
T T
x
T x t T erf s t a t
f a t
0
0
( , ) ( ), ( ) 2
( ) 2

= + =
(7)
where
k
a
c
2
= = is the diffusion coefficient and 0 > is the unique solution to the

equation

Ste
E x x ( ) , 0
= > (8)
with

x
erf x u du E x x erf x x
2 2
0
2
( ) exp( ) , ( ) ( ) exp( ),
= =
(9)

f
c T T
Ste
0
( )
: Stefan number
, (10)
and the total heat flux at the fixed face = 0 x is given by
= =

2
0
( ) (0, ) ( )exp( )
t
x
Q t kT d s t . (11)
Proof.
We have the following properties:
E E E x x (0) 0, ( ) , ( ) 0, 0 = + = + > > . (12)
Remark 1.
From (4) we have
( ) ( ) ( ) ( ) ( )
x t
T s t t s t T s t t t , , 0, 0 + = > (13)
www.intechopen.com

442
and therefore the Stefan condition (5) is transformed in
( ) ( ) ( ) ( ) ( )
x t xx
k
kT s t t T s t t T s t t t
c
2
, , ( ), , 0 = = >
(14)
which implies that the problem (2)-(6) is always a nonlinear problem (Pekeris & Slichter,
1939).
Remark 2.
A generalization of the Lam-Clapeyron solution is given in (Menaldi & Tarzia, 2003) for a
particular source in the heat equation. A study of the behaviour of the Lam-Clapeyron
solution when the latent heat goes to zero is given in (Guzman, 1982; Sherman, 1971).
2.2 The pseudo-steady-state approximation for the one-phase problem
An approximated solution to problem (2)-(6) is given by the pseudo-steady-state
approximation which must satisfy the following conditions: (3)-(6) and the steady-state
equation
( )
xx
T x s t t 0 , 0 , 0 = < < > . (15)
Theorem 2 (Stefan, 1989a)
The solution to the problem (15), (3)-(6) is given by
( )
f
T T
T x t T x x s t t
s t
0
0
( , ) , 0 , 0
( )
= + < < > (16)

ap ap
Ste
s t a t ( ) 2 ,
2
= = (17)
Proof.
The solution to (15), (3) and (4) is given by (16). Therefore the condition (5) is transformed in
the ordinary differential equation

f
k T T s t s t
0
( ) / ( ) ( ) = (18)
with the initial condition (6), whose solution is given by

f
Ste
s t k T T t a t
2 2
0
( ) 2 ( ) /( ) 4
2
= = (19)
that is

f
k T T
s t t
0
2 ( )
( )
(20)
Remark 3.
If the Stefan number is very small, i.e.

f
c T T
Ste
0
( )
1
= <<
(21)
www.intechopen.com

443
then the solution to the equation (8) for the Lam-Clapeyron solution can be taken as
ap
,
given in (17). This can be obtained by using the following first approximation:
x f x x x
2
2
exp( ) 1, ( ) , 0 1
< << . (22)

Remark 4.
A study of sufficient conditions on data to estimate the occurrence of a phase-change
process is given in (Solomon et al., 1983; Tarzia & Turner, 1992 & 1999).
2.3 The heat balance integral method (Goodman method) and the approximate
solution for the one-phase problem
An approximated solution for the following fusion problem (similar to the solidification
problem (2)-(6))
( )
t xx
cT kT x s t t 0 , 0 , 0 = < < > (23)
( ) T t T t
0
0, 0 , 0 = > > (24)
( ) ( )
T s t t t , 0 , 0 = > (25)
( ) ( ) ( )
x
kT s t t s t t , , 0 = > (26)
s(0) 0 = (27)
is given by the heat balance integral method, known by the Goodman method
(Goodman,1958). This method consists of replacing the Stefan condition (26) by

x xx
T s t t T s t t t
c
2
( ( ), ) ( ( ), ), 0 = >
(28)
and the heat equation (23) by its integral on the domain s t (0, ( )) given by

s t s t s t s t
t t xx
x x x
d k
T x t dx T x t dx T s t t s t T x t dx T x t dx
dt c
k k
T s t t T t s t T t
c c k
( ) ( ) ( ) ( )
0 0 0 0
( , ) ( , ) ( ( ), ) ( ) ( , ) ( , )
[ ( ( ), ) (0, )] [ ( ) (0, )]

= + = =
= = +

(29)
that is

s t
x
d k
T x t dx s t T t
dt c k
( )
0
( , ) [ ( ) (0, )]
= +

. (30)
In order to solve (30), (28), (24), (25) and (27), we propose an approximated temperature
profile
T x t t s t x t s t x x s t t
2
( , ) ( )( ( ) ) ( )( ( ) ) , 0 ( ), 0 = + < < > (31)
www.intechopen.com

444
where t t ( ), ( ), = = and s s t ( ) = are real functions to be determined. Firstly, we can
obtain and as a function of s and, therefore, we solve the corresponding ordinary
differential equation for s s t ( ) = .
Theorem 3.
The Goodman approximated solution is given by:

+
=
1 2 1
( )
( )
Ste
t
c s t
,
+
=
0
2
( ) ( )
( )
( )
t s t T
t
s t
(32)

+ +
= = =
+ + +
0
1 2 1 2
( ) 2 , 3 ,
5 1 2
g g
Ste Ste cT
s t a t Ste
Ste Ste
(33)
Remark 5.
Other refinements of the Goodman method are given in (Bell, 1978; Lunardini, 1981;
Lunardini 1991). In (Reginato & Tarzia, 1993; Reginato et al, 1993; Reginato et al., 2000) the
heat balance method was applied to root growth of crops and the modelling nutrient
uptake. In (Tarzia, 1990a) the heat balance method was applied to obtain the exponentially
fast asympotic behaviour of the solutions in heat conduction problems with absorption.
2.4 The Stefan solution for the planar phase-change surface moving with constant
speed
When the phase-change interface is moving with constant speed we can consider the
following inverse Stefan problem: find the temperature T T x t ( , ) = and f t T t ( ) (0, ) = such
that:

xx t
k
T T x s t t
c
, 0 ( ), 0 ( )

= < < > = (34)
T s t t t ( ( ), ) 0, 0 = > (35)

x
kT s t t s t t ( ( ), ) ( ), 0 = > (36)
( ) s t m s s t mt ( ) 0, (0) 0 ( ) = > = = (37)
Theorem 4. (Stefan, 1989b & 1991)
The solution to (34)-(37) is given by

m
T x t mt x
c
( , ) [1 exp( ( ))]
(38)
and the temperature at the fixed face is variable in time given by the expression:

f
m t
f t T t T t
c
2
( ) (0, ) [exp( ) 1] 0 , 0
= = < = >
. (39)
Remark 6.
More details with respect to the inverse Stefan problem can be found in (Quilghini, 1967).
www.intechopen.com

445
2.5 The Solomon-Wilson-Alexiades model for the phase-change process with a
mushy region and its similarity solution for the one-phase case
We consider a semi-infinite material in the liquid phase at the melting temperature
f
T . We
impose a temperature
f
T T
0
< at the fixed face x 0 = , and the solidification process begins,
and three regions can be distinguished, as follows (Solomon et al., 1982):
i. the liquid phase, at temperature
f
T T = , occupying the region x r t t ( ), 0; > >
ii. the solid phase, at temperature
f
T x t T ( , ) < , occupying the region x s t t 0 ( ), 0 < < > ;
iii. the mushy zone, at temperature
f
T , occupying the region s t x r t t ( ) ( ), 0 < < > . We
make the following two assumptions on its structure:
a. the material in the mushy zone contains a fixed fraction (with constant 0 1 < < ) of
the total latent heat .
b. the width of the mushy zone is inversely proportional (with constant 0 > ) to the
temperature gradient at s t ( ) .
Therefore the problem consists of finding the free boundaries x s t ( ) = and x r t ( ) = , and the
temperature T T x t ( , ) = such that the following conditions are satisfied:
( )
t xx
cT kT x s t t 0 , 0 , 0 = < < > (40)
( )
f
T t T T t
0
0, , 0 = < > ; s r (0) (0) 0 = = (41)
( ) ( )
f
T s t t T t , , 0 = > (42)

x
kT s t t s t r t t ( ( ), ) [ ( ) (1 ) ( )], 0 = + > (43)

x
T s t t r t s t t ( ( ), )( ( ) ( )) , 0 = > . (44)
Theorem 5. (Solomon et al., 1982):
The explicit solution to problem (40)-(44) is given by:

f
T T
x
T x t T erf s t a t r t a t
erf a t
0
0
( )
( , ) ( ), ( ) 2 , ( ) 2
( ) 2

= + = = (45)
where

f
k
erf exp a
c
T T
2
0
( ) ( ),
2( )

= + =
(46)
and 0 > is the unique solution to the equation

0
( )
( ) , 0

= > =

f
c T T Ste
D x x Ste
(47)
with

f
D x xerf x x x erf x
T T
2 2 2
0
(1 )
( ) ( )exp( ) [exp( ) ( )]
2( )

= +
. (48)
www.intechopen.com

446
Remark 7.
The classical Lam-Clapeyron solution can be obtained for the particular case 1, 0 = = .
If the Stefan number is small, then an approximated solution for and is given by:

1
2
0
0
, [1 /( )]
2[1 (1 ) /( )]

= = +
+

f
f
Ste
T T
T T

. (49)
2.6 The Cho-Sunderland solution for the one-phase problem with temperature-
dependent thermal conductivity
We consider the following solidification problem for a semi-infinite material
cT x t k T T x t x s t t
x x
t
( , ) ( ( ) ( , )) , 0 ( ), 0 = < < > (50)
T t T T t
o
f
(0, ) , 0 = < > (51)
T s t t T t
f
( ( ), ) , 0 = > (52)
k T T s t t s t t
x
f
( ) ( ( ), ) ( ) , 0 = > (53)
where T(x,t) is the temperature of the solid phase, >0 is the density of mass, 0 > is the
latent heat of fusion by unity of mass, c >0 is the specific heat, x=s(t) is the phase-change
interface, T
f
is the phase-change temperature, T
o
is the temperature at the fixed face x=0. We
suppose that the thermal conductivity has the following expression:
k k T k T T T T
o o o
f
( ) [1 ( ) /( )] , = = + . (54)
Let
o
=k
o
/c be the diffusion coefficient at the temperature T
o
. We observe that if =0, the
problem (50)-(53) becomes the classical one-phase Lam-Clapeyron-Stefan problem.
Theorem 6. (Cho & Sunderland, 1974)
The solution to problem (50)-(54) is given by:

T T
o
f x
T x t T
o
t
o
( )
( , ) ( ) , , 0
( ) 2

= + = < <
(55)

o
s t t ( ) 2 = (56)

where x x ( ) ( )
= = is the modified error function, for > -1, the unique solution to the
following boundary value problem in variable x, i.e:

i x x x x x
ii
) [(1 ( )) ( )] 2 ( ) 0 , 0,
) (0 ) 0 , ( ) 1
+ + = >
+
= + =
(57)
www.intechopen.com

447
and the unknown thermal coefficients and must satisfy the following system of
equations:
( ) 0 = (58)

2 ( )
[1 ( )] 0
( ) ( ) c T T
f o

+ =

. (59)
Remark 8.
Explicit solutions are given in (Briozzo et al., 2007 & 2010; Briozzo & Tarzia, 2002; Natale &
Tarzia, 2006; Rogers & Broadbridge, 1988; Tirskii, 1959; Tritscher & Broadbridge, 1994)
where nonlinear thermal coefficients are considered and in (Natale & Tarzia, 2000; Rogers,
1986) for Storms materials.
2.7 The Neumann solution for the two-phase problem for prescribed surface
temperature at the fixed face
We consider a semi-infinite material with null melting temperature
f
T 0 = , with an initial
temperature C 0 < and having a temperature boundary condition B 0 > at the fixed
face x 0 = . The model for the two-phase Lam-Clapeyron-Stefan problem is given by: find
the free boundary x s t ( ) = , defined for t 0 > , and the temperature T T x t ( , ) = defined by

2
1
( , ) 0 ( ), 0
( , ) ( ), 0
( , ) ( ) , 0
> < < >
= = >
< < >
f
f
f
T x t T if x s t t
T x t T if x s t t
T x t T if s t x t
(60)
for 0 x > and 0 t > , such that (i=1: solid phase; i=2: liquid phase):

t xx
c T k T x s t t
2 2 2 2
0, 0 ( ), 0 = < < > , (61)

t xx
c T k T x s t t
1 1 1 1
0, ( ), 0 = > > , (62)

T x C x
1
( , 0) 0, 0, = < > (63)
T t B t
2
(0, ) 0, 0, = > > (64)

f
T s t t T t
1
( ( ), ) 0, 0 = = > , (65)

f
T s t t T t
2
( ( ), ) 0, 0 = = > , (66)
( ) ( )
x x
k T s t t k T s t t s t t
1 1 2 2
( ), ( ), ( ), 0 = > , (67)
s(0) 0 = . (68)
www.intechopen.com

448
Theorem 7. (Neumann solution (Webber, 1901))
The explicit solution to problem (61)-(68) is given by:

B x
T x t B erf x s t t
erf a a t
2
2 2
( , ) ( ), 0 ( ), 0
( / ) 2
= > (69)

B x
T x t C erfc s t x t
erfc a a t
1
1 1
( , ) ( ), ( ) , 0
( / ) 2
= + > (70)

k k
s t t a a
c c
2 2 2 1
2 1
2 1
( ) 2 ( , )

= = = (71)

where 0 > is the unique solution to the following equation:
( ) , 0 = > F x x x (72)
where

Bk x Ck x
F x F F
a a a a
2 1
2 1
2 1 2 1
( ) ( ) ( )

=

(73)

exp x exp x
F x F x erfc x erf x
erfc x erf x
2 2
1 2
( ) ( )
( ) , ( ) , ( ) 1 ( )
( ) ( )

= = = . (74)

Remark 9.
It is very interesting to answer the following question: When is the Neumann solution for a
semi-infinite material applicable to a finite material x
0
(0, ) ? (Solomon, 1979).
Taking into account that erf x for 2 x ( ) 1 , we have an affirmative answer for a short
period of time becauseT x t C
1 0
( , ) is equivalent to

x
erf
a t
0
1
( ) 1
2
(75)
that is

x
t
a
2
0
2
1
16
. (76)
Remark 10.
A generalization of Neumann solution is given in (Briozzo et al, 2004 & 2007b) for particular
sources in the heat equations for both phases. A study of the behaviour of the Neumann
solution when the latent heat goes to zero is given in (Tarzia & Villa, 1991). A generalization
of Neumann solution in multi-phase media is given in (Sanziel & Tarzia, 1989; Weiner, 1955;
Wilson, 1978 & 1982), and when we have shrinkage or expansion (Fi & Han, 2007; Natale et
al., 2010; Wilson & Solomon, 1986).
www.intechopen.com

449
2.8 The Neumann-type solution for the two-phase problem for a particular prescribed
heat flux at the fixed face, and the necessary and sufficient condition to have an
instantaneous phase-change process
If we consider the problem (61)-(68) by changing the boundary condition (64) at x 0 = by a
heat flux condition of the type

x
q
k T t
t
0
2 2
(0, ) = (77)
then we can obtain the following result:
Theorem 8. (Tarzia, 1981)
i. If q
0
verifies the inequality

Ck
q
a
1
0
1

> (78)
then we have an instantaneous change of phase and the corresponding explicit solution is
given by:

x
T x t A B erf x s t t
a t
2 2 2
2
( , ) ( ), 0 ( ), 0
2
= + > (79)

x
T x t A B erf s t x t
a t
1 1 1
1
( , ) ( ), ( ) , 0
2
= + > (80)

k k
s t w t a a
c c
2 2 2 1
2 1
2 1
( ) 2 ,

= = =

(81)
where

erf w a C
A w C B w
erfc w a erfc w a
1
1 1
1 1
( / )
( ) , ( )
( / ) ( / )
= = (82)

a q a q
A w erf w a B w
k k
2 0 2 0
2 2 2
2 2
( ) ( / ), ( )

= = (83)
and w 0 > is the unique solution to the equation

0
( ) , 0 = > F x x x , (84)
where

q Ck
F x x a F x a
a
2 2 0 1
0 2 1 1
1
( ) exp( / ) ( / )

= . (85)
ii. If q Ck a
0 1 1
/ the corresponding problem represents only a heat conduction
problem for the initial solid phase, and the temperature is given by
www.intechopen.com

450

q a x
T x t T x t C erfc x t
k a t
0 1
1
1 1
( , ) ( , ) ( ), 0, 0
2
= = + > > . (86)

Corollary 9 (Tarzia, 1981)
The coefficient that characterizes the free boundary s t t ( ) 2 = of Neumann solution
(69)-(74) must satisfy the following inequality:

B k c
erf
a C k c
2 2
2 1 1
( )
< . (87)
2.9 The Neumann-type solution for the two-phase problem for a particular prescribed
convective condition (Newton law) at the fixed face, and the necessary and sufficient
condition to have an instantaneous phase-change process
We consider the following free boundary problem: find the solid-liquid interface
x s t ( ) = and the temperature T x t ( , ) defined by

s
f
l
T x t if x s t t
T x t T if x s t t
T x t if x s t t
( , ) 0 ( ), 0,
( , ) ( ), 0,
( , ) ( ), 0,
< < >
= = >
> >
(88)
which satisfy the following equations and boundary conditions

t xx
s s s
T T x s t t , 0 ( ), 0 = < < > (89)

t xx
l l l
T T x s t t , ( ), 0 = > > (90)

s l f
T s t t T s t t T x s t t ( ( ), ) ( ( ), ) , ( ), 0 = = = > (91)

l l i
T x T t T x t ( , 0) ( , ) , 0, 0 = + = > > (92)

x
s s s
h
k T t T t T t
t
0
(0, ) ( (0, ) ), 0
= > (93)

x x
s s l l
k T s t t k T s t t s t t ( ( ), ) ( ( ), ) ( ), 0 = > (94)
s(0) 0 = (95)
where the subscripts s and l represent the solid and liquid phases respectively, is the
common density of mass and is the latent heat of fusion, and
f i
T T T
< < . We have the

following results:
Theorem 10. (Tarzia, 2004)
If the coefficient h
0
verifies the inequality

i f
l
i l
T T
k
h
T T
0
>
(96)
www.intechopen.com

451
there exists an instantaneous solidification process and then the free boundary problem (89)-
(95) has the explicit solution to a similarity type given by

l
s t t ( ) 2 = (97)

s
f
s s
s
s l
s s
h x
T T erf
k t
T x t T
h
erf
k
0
0
( )[1 ( )]
2
( , )
1 ( )
+
= +
+
(98)

l
l i i f
x
erfc
t
T x t T T T
erfc
( )
2
( , ) ( )
( )
= (99)
and the dimensionless parameter 0 > satisfies the following equation
F x x x ( ) , 0 = > (100)
where function F and the bs coefficients are given by

bx x
F x b b
erfc x b erf x b
2 2
1 3
2
exp( ) exp( )
( )
( ) 1 ( )

=
+
(101)

f
l
s l
h T T
b b
0
1
( )
0; 0
= > = >
(102)

l i f
s
s
c T T
h
b b
h
0
2 3
( )
0; 0
= > = >
(103)
Proof.
Function F has the following properties:

f l i f
l
h T T c T T
F b b
0
1 3
( ) ( )
(0 )

+

= =

(104)
( ) F F x x ( ) , 0, 0 + = < > (105)
Therefore, there exists a unique solution >0 of the Eq. (100) if and only if F(0 ) 0
+
> , that is
inequality (96) holds.
2.10 The similarity solution for the two-phase Lam-Clapeyron-Stefan problem with a
mushy region
We consider a semi-infinite material initially in the solid phase at the
temperature
f
C T 0 < = . We impose a temperature
f
B T 0 > = at the fixed face x 0 = , and the
fusion process begins, and three regions can be distinguished, as follows: (Tarzia, 1990b):
www.intechopen.com

452
i. the liquid phase, at temperature T T x t
2 2
( , ) 0 = > , occupying the region
x s t t 0 ( ), 0; < >
ii. the solid phase, at temperatureT T x t
1 1
( , ) 0 = < , occupying the region x r t t ( ), 0 > > ;
iii. the mushy zone, at temperature
f
T 0 = , occupying the region s t x r t t ( ) ( ), 0 < < > . We
make the following two assumptions on its structure:
a. the material in the mushy zone contains a fixed fraction (with constant 0 1 < < ) of
the total latent heat ;
b. the width of the mushy zone is inversely proportional (with constant 0 > ) to the
temperature gradient at s t ( ) .
Therefore, the problem consists of finding the free boundaries x s t x r t ( ), ( ) = = , and the
temperature:

T x t if x s t t
T x t if s t x r t t
T x t if r t x t
2
1
( , ) 0 0 ( ), 0
( , ) 0 ( ) ( ), 0
( , ) 0 ( ) , 0
> < < >
= >
< < >
(106)
defined for x 0 > andt 0 > , such that the following conditions are satisfied:

xx t
T T x s t t
2 2 2
, 0 ( ), 0 = < < > (107)

xx t
T T r t x t
1 1 1
, ( ) , 0 = < > (108)
s r (0) (0) 0, = = (109)
T s t t T r t t t
2 1
( ( ), ) ( ( ), ) 0, 0 = = > (110)

x x
k T r t t k T s t t s t r t
1 1 2 2
( ( ), ) ( ( ), ) [(1 ) ( ) ( )], = + (111)

x
T s t t r t s t t
2
( ( ), ) ( ( ) ( )) , 0 = > (112)
T x T t C x t
1 1
( , 0) ( , ) , 0, 0 = + = > > (113)
T t B t
2
(0, ) 0, 0 = > > (114)
Theorem 11. (Tarzia, 1990b)
i. The explicit solution to the problem (107)-(114) is given by

x x
T x t A B erf T x t A B erf
a t a t
1 1 1 2 2 2
1 2
( , ) ( ), ( , ) ( )
2 2
= + = + (115)

k k
s t t r t t a a
c c
2 2 2 1
2 1
2 1
( ) 2 , ( ) 2 ( , )

= = = = (116)
where
www.intechopen.com

453

Cerf
B C a
A B B A B
erf erfc erfc
a a a
1
2 2 1 1
2 1 1
( )
, , ,
( ) ( ) ( )

= = = = (117)

a
erf
B a a
2
2
2
2 2
( ) exp( ) ( )
2

= = + (118)
where 0 > is the unique solution to the equation

1 2
( ) ( ), 0 = > K x K x x (119)
with

k B x k B x x
K x F F F x
a a erfc x a a
a x x x
K x x erf F x
B a a erf x
2
2 1
1 2 1 1
2 1 2 1
2 2
2
2 2 2
2 2
( ) exp( )
( ) ( ) ( ), ( )
( )
exp( )
( ) [ exp( ) ( )], ( )
2 ( )
= =
= + =
(120)
Proof.
We have the following properties
K K K x
1 1 1
(0 ) , ( ) , 0, 0
+
= + + = < > , (121)
K K K x
2 2 2
(0 ) 0, ( ) , 0, 0
+
= + = + < > , (122)
and the thesis holds.
Remark 11
If the boundary condition (114) is replaced by a heat flux condition of the type (77) then we
will have an instantaneous change of phase if and only if the coefficient q
0
that characterizes
the heat flux (77) verifies an inequality (Tarzia, 1990b).
2.11 The similarity solution for the phase-change problem by considering a density
jump
We will consider the two-phase Lam-Clapeyron-Stefan problem for a semi-infinite material
taking into account the density jump under the change of phase. We will find the interface
s s t ( ) 0 = > (free boundary), defined for t 0 > , and the temperature

x t if x s t t
x t if x s t t
x t if x s t t
1
2
( , ) 0 0 ( ), 0,
( , ) 0 ( ), 0,
( , ) 0 ( ), 0,
< < < >
= = >
> > >
(123)
defined for x 0 > and t 0 > , such that they satisfy the following conditions:

xx t
x s t t
1 1 1
, 0 ( ), 0 = < < > (124)
www.intechopen.com

454

xx x t
s t x s t t
1 2
2 2 2 2
2
( ) , ( ), 0

+ = > > (125)

s t t t
s t t t
1
2
( ( ), ) 0, 0
( ( ), ) 0, 0
= >
= >
(126)

x x
k s t t k s t t s t t
1 1 2 2 1
( ( ), ) ( ( ), ) ( ), 0 = > (127)
x x
2 0
( , 0) 0, 0 = > > (128)
s(0) 0 = (129)
t d t
1
(0, ) 0, 0 = < > . (130)
Theorem 12 (Carslaw & Jaeger, 1959; Rubinstein, 1971)
The explicit solution to the free boundary problem (124)-(130) is given by

x
x t A B erf
a t
x
x t A B erf
a t
s t t
1 1 1
1
2 2 2 1
2
( , ) ( )
2
( , ) ( )
2
( ) 2 , 0
= +
= + +
= >
(131)
where

d
A d B
erf a
1 1
1
( ) , ( )
( )

= = , (132)

erf a
A B
erfc a erfc a
0 0 0
2 2
0 0
( )
( ) , ( )
( ) ( )

= = , (133)

a
a
a
1 2 2
1 0
2 2
, ,
1

= = =
+
, (134)
and is the unique solution to the following equation:
F x x x ( ) , 0 = > , (135)
with

k x k x
F x B x B x
a a a a
2 2
1 2
1 2 2 2
1 0 1 1 1 2
( ) ( )exp ( )exp

=

. (136)
Proof.
www.intechopen.com

455
F F F x x (0 ) , ( ) , ( ) 0, 0
+
= + + = < > . (137)
Theorem 13 (Bancora & Tarzia, 1985)
i. If we replace the boundary condition (130) by the following one given by:

x
q
k t t
t
0
1 1
(0, ) , 0 = > , (138)
then there exists an explicit solution corresponding to the free boundary problem (124)-(130)
and (138) if and only if the coefficient q
0
satisfies the inequality

k c
q
2 2 2
0 0

> . (139)
In this case the explicit solution is given by:

x
x t C D erf
a t
x
x t C D erf
a t
s t w t w
1 1 1
1
2 2 2 2
2
( , ) ( )
2
( , ) ( )
2
( ) 2 , 0
= +
= + +
= >
(140)
where

a q a q w
C w erf D w
k a k
1 0 1 0
1 1
1 1 1
( ) ( ), ( )

= = , (141)

erf w a
C w D
erfc w a erfc w a
0 0 0
2 2
0 0
( )
( ) , ( )
( ) ( )

= = , (142)

a
a
a
1 2 2
2 0
2 2
, ,
1

= = =
+
, (143)
and w is the unique solution to the following equation:
F x x x
0
( ) , 0 = > , (144)
with

q x k x
F x B x
a a a
2 2
0 2
0 2 2 2
1 1 0 1 2
( ) exp ( )exp

=

. (145)
Proof.

k
F q F F x x
a
2 0
0 0 0
1 2
1
(0 ) , ( ) , ( ) 0, 0

+

= + = < >

. (146)
www.intechopen.com

456
Remark 12 .
When the boundary condition at the fixed face x 0 = is given by (93) the explicit solution
was given in (Tarzia, 2007).
2.12 The determination of one or two unknown thermal coefficient through an over-
specified condition at the fixed face for one or two-phase cases
We consider the one-phase Lam-Clapeyron-Stefan problem with unknown thermal
coefficients. If we give an overspecified boundary condition at the fixed face x 0 = we can
determine one or two unknown coefficients following (Arderius et al., 1996; Cannon, 1963 &
1964; Garguichevich et al., 1985; Jones, 1962 & 1963; Tarzia, 1982,1983 & 1984)).
2.12.1 Determination of one unknown thermal coefficient through a one-phase case
The problem consists of finding the free boundary x s t ( ) = , the temperature T T x t ( , ) = , and
one unknown thermal coefficient chosen among { } k c , , , such that they must satisfy the
following conditions (we have a free boundary problem):

t xx
k
T T x s t t a
c
2
, 0 ( ), 0 ( )

= < < > = = (147)
T s t t t ( ( ), ) 0, 0 = > (148)

x
kT s t t s t t ( ( ), ) ( ), 0 = > , (149)
T t T
0
(0, ) 0 = > (150)

x
q
kT t t
t
0
(0, ) , 0 = > , (151)
s(0) 0 = , (152)
where q
0
is the coefficient that characterizes the heat flux at the fixed face x 0 = and it must
be obtained experimentally.
Theorem 14 (Tarzia, 1982)
Let T
0
and q
0
be determinated experimentally. The solution for the determination of one
thermal coefficient is given by:

T x
T x t T erf x s t t
a t
erf
a
0
0
( , ) ( ), 0 ( ), 0
2
( )
= < < > , (153)

s t t t ( ) 2 , 0 = > , (154)
where and the unknown thermal coefficient are computed in the summarized way in the
following Table 1:
www.intechopen.com

457
Case #
Formulae for
unknown
coefficients
Parameter is the
unique solution to the
equation
Restrictions on
data
1
k
c
q c
k
2
0
exp( )

=

T k c
erf x
q
x
0
0
( )
0
=
>

T k c
q
0
0
1
<
2
kT
erf q
q
erf
kcT
0
0
2
2 0
2
0
( )
( )

=
=

cT
E x
x
0
( )
0

=
>

-------
3
q
q
k erf
cT
2 0
2
2 0
2
0
exp( )
( )

=
=

As in Case 2 -------
4
q
q
c erf
kT
2 0
2
2 0
2
0
exp( )
( )

=
=

erf x k T
x
x q
x
2 0
2
0
( )
exp( )
0
=
>
k T
q
0
2
0
1
2
<

Table 1. Summary of the determination of one thermal coefficient through a one-phase
Lam-Clapeyron-Stefan problem (4 cases)
Remark 13.
The determination of one unknown thermal coefficient for phase-change problems with
temperature-dependent thermal conductivity of the type (54) was given in (Tarzia, 1998).
2.12.2 Determination of two unknown thermal coefficients for the one-phase case
If the interface solid-liquid is given by the law:
s t t t ( ) 2 , 0 = > (155)
where the coefficient 0 > was determined experimentally then the problem consists of
finding the temperature T T x t ( , ) = and two unknown thermal coefficients chosen among
{ } k c , , , such that they satisfy the conditions (147)-(152) (we have a moving boundary
problem).
Theorem 15 (Tarzia, 1983)
Let T
0
, and q
0
be determined experimentally. The solution for the determination of two
thermal coefficients is given by

T x
T x t T erf x s t t
a t
erf
a
0
0
( , ) ( ), 0 ( ), 0
2
( )
= < < > , (156)

www.intechopen.com

458
and the two coefficients are computed in the summarized way in the following Table 2.

Case #
Formulae for
unknown
coefficients
Parameter is the
unique solution to
the equation
Restrictions on data
1
q
c erf
T
q erf
k
T
0
0
0
0
( )
( )
=
=

q
0
log( )

=

q
0
1

>

2
q
c
k
2 0
2
2
exp( )
=
=

c T
xerf x
q
x
0
0
( )
0

=
>
---------
3
q
q c
k
2 0
2
0
2
exp( )
exp( )

cT
E x
x
0
( )
0

=
>

---------
4
k
c
q
2
2
2 0
exp( )
=
=

erf x kT
x q
x
0
0
( )
0

=
>

kT
q
0
0
1
2
<
5
k
c
c q
k
2
2
2
0
2
exp( )

=

As in case 4 As in case 4
6
q
k
c
q
2 0
2 2
0
exp( )
exp( )

=
=

As in case 4 As in case 4
Table 2. Summary of the simultaneous determination of two thermal coefficients through a
one-phase Lam-Clapeyron-Stefan problem (6 cases)
Remark 14.
The determination of thermal coefficients for the Solomon-Wilson-Alexiades mushy region
was obtained in (Tarzia, 1987). The simultaneous determination of two unknown thermal
coefficients for phase-change problems with temperature-dependent thermal conductivity
of the type (54) was given in (Salva & Tarzia, 2010) with a sensitivity analysis.
2.12.3 Determination of one or two unknown thermal coefficients through a two-phase
Lam-Clapeyron-Stefan problem
The determination of one or two unknown thermal coefficients for a two-phase
solidification or fusion problem was obtained in (Stampella & Tarzia, 1989). The
www.intechopen.com

459
determination of one or two unknown thermal coefficients for a mushy region was obtained
in (Gonzlez & Tarzia, 1996). Another variant for the simultaneous determination of the
thermal coefficients is given in (Tarzia, 1991a).
Remark 15.
Explicit solutions for the determination of unknown coefficients are given in (Briozzo et al.,
1999) for Storms type materials.
2.13 A similarity solution for the thawing in a saturated porous medium by
considering a density jump and the influence of the pressure on the melting
temperature
We consider the problem of thawing of a partialIy frozen porous medium, saturated with an
incompressible liquid. For a detailed exposition of the physical background we refer to
(Charach & Rubinstein, 1992; Fasano et al. 1993; Fasano & Primicerio, 1984; Nakano, 1990;
ONeill & Miller, 1985; Talamucci, 1997 & 1998). More specifically, we deal with the
following situations:
i. a sharp interface between the frozen part and the unfrozen part of the domain exists
(sharp, in the macroscopic sense);
ii. the frozen phase is at rest with respect to the porous skeleton, which will be considered
to be undeformable;
iii. due to the density jump between the liquid and solid phases, thawing can induce either
desaturation or water movement in the melting regon. We will consider the latter
situation, assuming that liquid is continuously supplied to keep the medium saturated.
The unknowns of the problem are the function x=s(t), representing the free boundary, and
the two functions u(x, t) and v(x, t) representing the temperature of the unfrozen and of the
frozen zone respectively which must satisfy the following conditions:

t xx x
u a u b s t u x s t t
1
( ) , 0 ( ), 0 = < < > (157)

t xx
v a v x s t t
2
, ( ), 0 = > > (158)
u s t t v s t t d s t s t t ( ( ), ) ( ( ), ) ( ) ( ), 0 = = > (159)

F x U x
k v s t t k u s t t s t s t s t t
2
( ( ), ) ( ( ), ) ( ) ( ) ( ( )) , 0 = + > (160)
u t B t (0, ) 0, 0 = > > . (161)
v x v t A x t ( , 0) ( , ) 0, 0, 0 = + = < > > (162)
s(0) 0 = (163)
with

U F W W
U U F F U U
W I I W I
I I W I
W
k k c
a a b d
c c c K
c c
d c c
K
2 2
1 1 2 2
2
, , ,
( )
, , ( )

= = = = = =

= = = =
(164)
www.intechopen.com

460
where
: porosity,
W I
and : density of water and ice,
c: specific heat at constant density,
U F
k k and : conductivity of the unfrozen and frozen zones,
u=v=0 : the melting point at atmospheric pressure,
: latent heat at u=0 ,
: coefficient in the CIausius-CIapeyron law,
0 > : viscosity of liquid,
K>0: hydraulic permeability,
B>0: boundary temperature at the fixed face x=0,
-A<0: initial temperature.
Theorem 16 (Fasano et al., 1999)
The free boundary problem (157) (163) has the similarity solution
s t t
1
( ) 2 , = (165)

x
t
m B
u x t B r pr dr
g p
1
( 2
2
2
0
( , ) exp( )
( , )
= + +
(166)

x x
m erfc A erf erf
t t
v x t
erfc
2
0
2 2
0
( )
2 2
( , )
( )

+

= (167)
if and only if the coefficient 0 > satisfies the following equation:
K B my H p y K F m y y y y
2 3
1 2
( ) ( , ) ( , ) , 0 = + > , (168)
where

y
p y
g p y r pyr dr H p y
g p y
2
2
0
exp(( 1) )
( , ) exp( ) , ( , )
2 ( , )

= + =
(169)

y
F m y A my
erfc y
2 2
2 0
0
exp( )
( , ) ( )
( )
= + (170)
and the constants K K p m and
1 2 0
, , , , , are defined as follows:

U F
k k
K K
1
1 2 0 1
2 1 2
0, 0, 0, 0

= > = > = > = > (171)
p b m d
2 3
1 1
2 , 2 0, 2 = = > = (172)
Moreover, the existence and uniqueness of the unknown coefficient 0 > depends on the
sign of the three parameters p m , and of the problem.
www.intechopen.com

461
If we replace the boundary condition (161) by the following one:

U x
q
k u t t
t
0
(0, ) , 0 = > , (173)
then we can consider the free boundary problem (157) (160), (173), (162)-(163) and we can
obtain the following result.
Theorem 17 (Lombardi & Tarzia, 2001)
The free boundary problem (157) (160), (173), (162)-(163) has the following similarity
solution:
s t t
1
( ) 2 ,
= (174)

x
t
U U
q q
u x t m g p r pr dr
K K
1
(2
2 2 0 1 0 1
0
2 2
( , ) ( ) ( , ) exp( )

= + +
(175)

m Aerf m x
v x t erf
erfc erfc t
2 2
0
0 0 1
( ) ( ) ( )
( , )
( ) ( ) 2

+
=

(176)
if and only if the coefficient 0
> satisfies the following equation:

q p y K F m y y y y
2 3
0 2
exp(( 1) ) ( , ) , 0 = + > , (177)
or its equivalent
Q y q y
0 0
( ) , 0 = > , (178)
where

K F m y y y
Q y y
p y
3
2
0 2
( , )
( ) , 0
exp(( 1) )
+ +
= >
. (179)
Moreover, the existence and uniqueness of the unknown coefficient 0
> depends on the

sign of the three parameters p m and , of the problem.
3. Explicit solutions for free boundary problems for the diffusion equation
Heat and mass transfer with phase change problems, taking place in a porous medium, such
as evaporation, condensation, freezing, melting, sublimation and desublimation, have wide
application in separation processes, food technology, heat and mixture migration in soils
and grounds, etc. Due to the non-linearity of the problem, solutions usually involve
mathematical difficulties and only a few exact solutions have been found. Mathematical
formulation of the heat and mass transfer in capillary porous bodies has been established by
(Luikov, 1964, 1966, 1975 & 1978). Some books or booklets on the subject for the diffusion
equation are (Crank, 1975; Duvaut, 1976; Fasano, 2005; Fasano & Primicerio, 1986; Froment
& Bischoff, 1979; Levenspiel, 1962; Primicerio & Gianni, 1989; Szekely et al., 1976).
www.intechopen.com

462
3.1 The oxygen diffusion-consumption problem and its relationship with the phase-
change problem
The diffusion-consumption of oxygen in absorbing tissue consists in finding the free
boundary x=s(t) and the concentration C(x,t) such that they satisfy the following conditions
(Crank & Gupta, 1972; Crank, 1984; Liapis et al., 1982):

t xx
C C x s t t T
0
1, 0 ( ), 0 = < < < < (180)
( ) ( ) C x H x x , 0 , 0 1 = , (181)
( )
x
C t G t C t F t t T
0
(0, ) ( ) ( 0, ( )), 0 = = < < , (182)
( ) C s t t t T
0
( ), 0, 0 = < < , (183)

x
C s t t t T
0
( ( ), ) 0, 0 = < < , (184)
s t T
0
(0) 1, 0 = < < . (185)
Remark 16.
We remark that s t ( ) does not appear in both conditions (183) and (184) on the free
boundary x=s(t) and for this reason the free boundary problem is of an implicit type; the
Lam-Clapeyron-Stefan problem is one of an explicit type.
Theorem 18. (Fasano, 1974; Schatz, 1966)
The free boundary problem of the diffusion-consumption of oxygen in absorbing tissue is
equivalent to the following Lam-Clapeyron-Stefan problem:

t xx
z z x s t t T
0
0, 0 ( ), 0 = < < < < (186)
z x h x x ( , 0) ( ), 0 1 = (187)

x
z t g t z t f t t T
0
(0, ) ( ) ( (0, ) ( )), 0 = = < < (188)
z s t t t T
0
( ( ), ) 0, 0 = < < (189)

x
z s t t s t t T
0
( ( ), ) ( ), 0 = < < (190)
s t T
0
(0) 1, 0 = < < . (191)
i. If s z T
0
( , , ) is a solution to the problem (186) (191) then s C T
0
( , , ) is solution to the
problem (180) (185) where we define:

s t s t
x
C x t d z y t dy
( ) ( )
( , ) [1 ( , )] ,
= +

(192)

t
x
H x d h y dy F t H f d
1 1
0
( ) [1 ( )] , ( ) (0) ( )
= + = +

(193)
www.intechopen.com

463

t
G t h x dx g d
1
0 0
( ) 1 ( ) ( ) , = + +

(194)
ii. If s C T
0
( , , ) is a solution to the problem (180) (185) then s z T
0
( , , ) is solution to the
problem (186) (191) where we define:

t
z x t C x t ( , ) ( , ) = (195)
h x H x g t G t f t F t ( ) ( ) 1, ( ) ( ), ( ) ( ) = = = (196)
Remark 17.
The oxygen diffusion-consumption free boundary problem was applied to the anaerobiosis
in saturated soil aggregates in (Gonzlez et al., 2008).
3.2 The Rubinstein solution for the binary alloy solidification problem
We consider a semi-infinite slab of a binary alloy consisting of two components A, B. Let C
be the concentration of A. We suppose that solidification of the alloy is governed by an
equilibrium phase diagram consisting of liquidus curve
L
T f C ( ) = , and a solidus curve
S
T f C C ( ), 0 1 = < < and we assume
L S
f f , to be monotonically increasing,
L S
f C f C ( ) ( ) >
and
A B
L S cr L S cr
f f T f f T (0) (0) , (1) (1) = = = = . Material is in its solid state if
S
T f C ( ) and liquid
if
L
f C T ( ) . If
S L
f C T f C ( ) ( ) < < then the material state is not well defined (it is known as
mushy region).
We consider that the semi-infinite alloy is initially liquid at constant temperature
in
T and
concentration
in
C , for which
L in in
f C T ( ) . Beginning at time t 0 = , a cold temperature
A
B cr
T T < is imposed at x 0 = . Freezing occurs with, in principle, a sharp phase change front
s s t ( ) = separating solid alloy x s t ( ( )) < from liquid alloy x s t ( ( )) > . The mathematical
formulation of the solidification process is given in (Rubinstein, 1971) as follows:
Find temperature T x t ( , ) , concentration C x t ( , ) and phase-change front x s t ( ) = , such that
the following conditions must be satisfied:

t xx
S S S
C D C x s t t , 0 ( ), 0 = < < > (197)

t xx
C DC x s t t , ( ), 0 = > >

(198)

t xx
s s s
T T x s t t , 0 ( ), 0 = < < > (199)

t xx
T T x s t t , ( ), 0 = > >

(200)

cr S S S
T T s t t f C s t t t ( ( ), ) [ ( ( ), )], 0 = = > (201)

cr L
T T s t t f C s t t t ( ( ), ) [ ( ( ), )], 0 = = >

(202)

x x
S S S
DC s t t D C s t t C s t t C s t t s t t ( ( ), ) ( ( ), ) [ ( ( ), ) ( ( ), )] ( ), 0 = >

(203)
www.intechopen.com

464

x x
S S
k T s t t k T s t t s t t ( ( ), ) ( ( ), ) ( ), 0 = >

(204)

A
S B cr
T t T T t (0, ) , 0 = < > (205)

in L in
T x T f C t ( , 0) ( ), 0 = > >
(206)

in
C x C x ( , 0) , 0 = >
(207)

x
S
C t t (0, ) 0, 0 = > (208)
s(0) 0 = (209)
where k D , , , , represent the mass density, the thermal conductivity, the thermal
diffusivity, the mass diffusion and the latent heat of fusion, being S and the subscripts
that denote the solid and liquid phase respectively.
Theorem 19 (Rubinstein, 1971; Solomon et al., 1983)
There exists a unique solution to the coupled free boundary problem (197)-(209), moreover,
the solidus C s t t
1
( ( ), ) and liquidus C s t t
2
( ( ), ) concentrations as well as the phase-change
temperature
Cr
T T s t t T s t t
1 2
( ( ), ) ( ( ), ) = = are constants in time. The explicit solution is given by
the following expressions:
s t t
1
( ) 2 = (210)

S S cr S S L L
C x t C T f C f C ( , ) , ( ) ( ), = = = (211)

in L in
S
erfc x D t
C x t C C C
erfc D
( 2 )
( , ) ( )
( )
= +

(212)

S
S B cr B
erf x t
T x t T T T
erf
( 2 )
( , ) ( )
( )
= + (213)

in cr in
S
erf x t
T x t T T T
erf
( 2 )
( , ) ( )
( )

= +

(214)
where and
cr
T (or equivalently
S
C and C
) must satisfy the following conditions:

in L
s cr
S L
C C
Q D T W
C C
1
( ), ( )
= =

(215)
where
Q x x x erfc x A k
2
( ) exp( ) ( ), / = =

(216)
www.intechopen.com

465

B in B S
x
W x T T T AQ x
x
1
( )
( ) [( ) ( / )]
1 ( )

= + +
+
(217)

S
S S
k E x
x
k Q x
( )
( )
( / )

=

(218)
Proof.
The system of equations (215) has a unique solution because of the following properties
Q Q Q (0 ) 0, ( ) 1, 0
+
= + = > (219)

B in
W T W T W
c
1 1 1
(0 ) , ( ) , 0
+
= + = + >
. (220)
Remark 18.
Some other references on the binary alloy solidification problems are (Alexandrov &
Malygin, 2006; Gupta et al., 1997; Tien & Geigen, 1967; Tien & Koump, 1970; Tsubaki &
Boley, 1977; Voller, 2006, 2008 a&b; White, 1985; Wilson et al., 1982). In (Cirelli & Tarzia,
2010) the binary alloy solidification problem (197)-(209) is solved by changing the boundary
condition (205) at the fixed face by a heat flux boundary condition of the type (77) or a
convective boundary condition of the type (93).
3.3 The Zeldovich-Kompaneets-Barenblatt solution for the gas flow through a porous
medium
The porous medium equation for a unidimensional material is given by

m
t xx
u u m ( ) , 1 = > (221)
which appears in a natural way, mainly to describe processes involving fluid flow, heat
transfer or diffusion, the flow of a gas through a porous medium and groundwater
infiltration (Vzquez, 2007). The diffusion coefficient of the equation (221) is

m
D u mu
1
( )

= (222)
assuming u 0 . Equation (221) is parabolic only at those points where u 0 and it is in
general a degenerate parabolic equation because degenerates wherever u 0 = .
Its main qualitative property with respect to the classical heat or diffusion equation is the
finite propagation which implies the appearance of a free boundary that separates the
regions where the solution u 0 > (where there is gas, according to the standard
interpretation of u as a gas density), from the empty region where u 0 = .
There exists an explicit solution to the Eq. (221) with a free boundary given by:
Theorem 20 (Zeldovich & Kompaneets, 1950; Barenblatt, 1952; Pattle, 1959)
The function
www.intechopen.com

466

m
x
for x s t t
u x t
s t s t
for x s t t
1
2
1
1
1 ( ), 0
( , )
( ) ( )
0 ( ), 0
>
(223)
with

m
m m
s t t
m
1
1
2 ( 1)
( ) ( 1)
1
+
+
= +

(224)
is a solution to the Cauchy problem for the equation (221) with initial data
u x u x x
0
( , 0) ( ), = given by

m
x
for x s
u x
s s
for x s
1
2
1
0
1
1 (0)
( )
(0) (0)
0 (0)
>
. (225)
3.4 The Luikov solution for the coupled heat and mass transfer for a phase-change
process
3.4.1 Drying with coupled phase-change in a porous medium
A semi-infinite porous medium is dried by maintaining a heat flux condition at x 0 = of the
type
q
t
0
, with q
0
0 > . Initially, the whole body is at uniform temperature t
0
and uniform
moisture potential u
0
. The moisture is assumed to evaporate completely at a constant
temperature, evaporation point
v
t . It is also assumed that the moisture potential in the first
region x s 0 ( ) < < , is constant at
v
u , where x s( ) = locates the evaporation front at time
0 > . It is further assumed that the moisture in vapour form does not take away any
appreciable amount of heat from the system. Neglecting mass diffusion due to temperature
variation, the problem can be expressed as (Cho, 1975; Gupta, 1974; Luikov, 1978; Santillan
Marcus & Tarzia, 2003):

t t
x x x s
x
2
1 1
1 2
( , ) ( , ), 0 ( ), 0

= < < >

(226)

v
u u x s
1
, 0 ( ), 0 = < < > (227)

m
t t Lc u
x x s
x c
2
2 2 2
2 2
2
( , ) , ( ), 0

= + > >

(228)

m
u u
x x x s
x
2
2 2
2
( , ) ( , ), ( ), 0

= > >

(229)
www.intechopen.com

467

q t
k x
x
0 1
1
, at 0, 0
= = >
(230)
t t x
2 0
in 0, 0 = > = (231)
u u x
2 0
in 0, 0 = > = (232)

v
t s t s t t x s
1 2 0
( ( ), ) ( ( ), ) on ( ) = = > = (233)

v
u s u s u u x s
1 2 0
( ( ), ) ( ( ), ) on ( ) = = < = (234)

m
t t ds
k s k s L x s
x x dt
1 2
1 2
( ( ), ) ( ( ), ) (1 ) on ( )

+ = =

(235)
where t
1
: temperature of the dried porous medium; t
2
: temperature of the humid porous
medium; u
2
:mass-transfer potential of the humid porous medium;
i
i ( 1, 2) = : thermal
diffusivity of the phase i;
12
: ratio of thermal diffusivities from phase 1 to phase 2;
m
:
moisture diffusivity;
m
c : specific mass capacity; c
2
: specific heat capacity;
i
k i ( 1, 2) = :
thermal conductivity of the phase i;
k
k
k
2
21
1
= ;
v
m
v
u u
K Lc
c t t
0
0
2 0
( )
( )
=

: Kossovitch
number; L : latent heat evaporation of liquid per unit mass-transfer potential; : coefficient
of internal evaporation;
m
: density of moisture;
u m
L
1
= : Luikov number, and
m v
L k t t
1 1 0
(1 ) ( ) 0. = >
Theorem 21 (Santillan Marcus & Tarzia, 2003)
i. If the Luikov number is equal to one, and the coefficient q
0
verifies the condition

v
k t t
q K
2 0
0 0
1
( )
( 2)
2

> + , (236)
then there exists one and only one solution 0 > to the following equation:

v
q k K
F x xF x K x K x x x
k t t
1 0 2 2 21 0
1 1 0 0
1 0
2 2
( )[ ( ) 2 2] exp( ) 2 , 0
( )

+ + = >
(237)

Furthermore, the solution to the problem (226)-(235) is given by:

v
x
u x u x s
1
1
( , ) , 0 ( ), 0,
2

= < < > = (238)

v
q
t x erf erf x s
k t t
0 1
1
1 0
( , ) 1 ( ( ) ( )), 0 ( ), 0
( )
= + < < >
(239)
www.intechopen.com

468

u
erfc
L
u x x s
erfc
2
( , ) , ( ), 0
( )

= > > (240)

erfc erfc K
t x s
erfc erfc erfc
2 2 0
2
( ) ( )
( ) [ exp( ) exp( )] , ( ), 0
( ) ( ) ( )

= + > > (241)
s
1
( ) 2 = . (242)
ii. If the Luikov number is different than one, (that is
1 m
) and the coefficient
0
q verifies the condition

u v
u
L K t t
q k
L
0 0
0 2
1
1
1

> +

+

(243)
then there exist one and only one solution 0 > to the equation:
x x x ( ) ( ), 0 = > , (244)
where

v
q
x x P x
t t
1 0 2
0
( ) exp( ) ( )
( )
= +
(245)
x k F x k x
2 1 1
( ) ( ) = + (246)

u
u u u
L K x
P x k F F x
L L L
0
2 1 1
1
( ) ( ) ( )
1

. (247)
Furthermore, the solution to the problem (226)-(235) is given by (238)-(240) , (242) and
x
1
2

=

u
u
u
u
x
erfc
L
erfc erfc K L
t x s
L erfc erfc
erfc
L
0
2
( ) ( )
( ) , ( ), 0.
1 ( ) ( )

= + + > >

(248)
iii. If the Luikov number
u
L verifies the condition

u
L
K
0
1
1
>
+
(249)
www.intechopen.com

469
then the temperature distribution t
2
reaches to a minimum value which is smaller than the
initial temperature or its limit value at +. The minimum value is attained when the
dimensionless variable
x
1
2

= takes the value

( )
u u u
u u
erfc
L K L L
L erfc K L
0
0
( 1) 1
log
1 ( )

+

=

. (250)
3.4.2 Other free boundary problem in a porous medium
There are some explicit solutions for the following free boundary problems for the diffusion
equation corresponding to evaporation, freezing, sublimation or desublimation processes in
(Lin, 1981, 1982a & 1982b; Mikhailov, 1975 & 1976; Santillan Marcus & Tarzia, 2000a & b).
The simultaneous determination of one or two unknown thermal coefficients of a semi-
infinite material through a desublimation process with coupled heat and moisture flows is
given in (Santillan Marcus & Tarzia, 2007; Santillan Marcus, et al., 2008).
3.5 A mixed saturated-unsaturated flow problem representing absorption of water by
a soil with a constant pond depth at the surface and an explicit solution for a
particular diffusivity
In wet soils, zones of saturation develop naturally in the vicinity of impermeable strata,
surface ponds and subterranean cavities. Hydrology must be concerned with transient flow
through coexisting unsaturated and saturated zones. Models of advancing saturated zones
necessarily involve a nonlinear free boundary problem (Broadbridge & White, 1990; Knigh
& Philip, 1974; Philip, 1957 & 1958; Warrick & Broadbridge, 1992).
We consider a homogeneous soil which initially has some uniform volumetric water
content
n
. At times t 0 > , water is supllied at the surface x 0 = under pressure head
0
.
Then, a mixed saturated-unsaturated flow problem representing absorption of water by a
soil with constant pond depth at the surface is presented. At any time t the zone of
saturation extends from x 0 = to x s t ( ) = . Assuming Darcys law and neglecting gravity, the
water flux is given by
( ) v K
x
(251)
where is the soil water matric potential and K is the hydraulic conductivity.
In the saturated zone we have

s
x t x s t
s t
0
0
( , ) , 0 ( )
( )

= < < (252)
and we obtain a free boundary problem for the unsaturated zone:

s
s t t t ( ( ) , ) , 0
+
= > (253)
www.intechopen.com

470
D x s t t
t x x
( ) , ( ), 0

= > >

(254)

s
s
D s t t K t
x s t
0
( ) ( ( ) , ) , 0
( )

+

= >
(255)

n
x t x s t t ( , 0) ( , ) , ( ), 0 = + = > > (256)
s(0) 0 = (257)
where
x : spatial coordinate,
t : time,
: volumetric water content,
n
: initial volumetric water content,
s
: volumetric water content at saturation,
: soil water matric potential,
0
: pond depth,
s
: soil water potential at x s t ( ) = ,
s 0
< <
K : hydraulic conductivity,
s
K : hydraulic conductivity at saturation,
D : soil water diffusivity (
d
D K
d
= ).
We consider the free boundary (253)-(257) where the position s t ( ) of the free boundary and
the water content field x t ( , ) must be determined; and we restrict our attention to the
special functional form of the soil water diffusivity

a
D
b
2
( )
( )
(258)
where a , and b are positive constants. With this form of diffusivity the nonlinear diffusion
equation (254) may be transformed to a linear diffusion equation. We consider the following
parameter:

n
s n
b
C 1
= >
. (259)
Remark 19.
In (Briozzo & Tarzia, 1998) a closed-form analytic solution can be obtained for a nonlinear
diffusion model under conditions of ponding surface. The explicit solution depends on a
parameter C (determined by the data of the problem ), according to two cases: C C
1
1 < < or
C C
1
, where C
1
is a constant which is obtained as the unique solution to an equation.
This results complements the study given in (Broadbridge, 1990) in order to established
when the explicit solution is available. The behaviour of the bifurcaton parameter C
1
as a
function of the driving potential is studied with the corresponding limits for small and large
www.intechopen.com

471
values. We also prove that the sorptivity is continuously differentiable as a function of
variable C .
3.6 Estimation of the diffusion coefficient in a gas-solid system
Looking for a competitive separation process like as the permeation, the development and
optimal choice of membrane materials become necessary. On this subject, equations
modelling the permeation process are required. The parameters contained in such a model
must be obtained from simple experiments. The knowledge of solubility and diffusivity are
very important to solve the separation problem.
We consider a polymeric membrane swelling for a hydrocarbon solution. The following
assumptions are considered: Once the gaseous component reaches a threshold concentration
on the gas-polymer interface, it diffuses through the membrane in the x direction being
immobilized by a quickly and irreversible transformation. Then a swelling front is generated
whose position is given by the free boundary x= s(t) , t >0 with the initial condition s(0)=0 .
Moreover, the hydrocarbon diffusion coefficient D in the saturated o swollen region of the
polymer is considered a constant for each experimental condition. A free boundary model
(Castro et al., 1987; Crank, 1975; Villa, 1987) with an overspecified condition for the one-
dimensional diffusion equation under the preceding assumptions is given:

t xx
c Dc x s t t , 0 ( ), 0 = < < > (260)
c s t t t ( ( ), ) 0, 0 = > (261)

x
Dc s t t s t t ( ( ), ) ( ), 0 = > , (262)
c t C t
0
(0, ) 0, 0 = > > (263)

t
x
A Dc d t t
0
(0, ) , 0 = >
, (264)
s(0) 0 = , (265)
where c=c(x,t) denotes the concentration profile of the hydrocarbon in the swollen area, s(t)
gives the position at time t of the free interface and separates the two regions in the
membrane, the saturated and unsaturated, D is the unknown diffusion coefficient in the
system, and
0
, and C are positive parameters and A is a positive constant which must be
obtained experimentally.
Theorem 22. (Destefanis et al., 1993)
The concentration profile and the free boundary position are given by:

C x
c x t C erf x s t t
Dt
erf
D
0
0
( , ) ( ), 0 ( ), 0
2
( )
= < < > , (266)

s t t t ( ) 2 , 0 = > , (267)
www.intechopen.com

472
and the unknown coefficients and D are obtained by the following expressions:

D rf
A C A
erf
AC A
2 2 2
2 2 2 2
0
2
0
exp( 2 )
e ( )
4 4
( ) exp( )
2 2

= =
= =
(268)
where is the unique solution to the equation:

C
E x x
0
( ) , 0

= > . (269)
Remark 20.
The methodology used in this determination of the unknown diffusion coefficient is a
variant of those developed in (Tarzia, 1982 & 1983) for the determination of thermal
coefficients for a semi-infinite material through a phase-change process.
3.7 The coupled heat and mass transfer during the freezing of high-water content
materials with two free boundaries: the freezing and sublimation fronts
Ice sublimation takes place from the surface of high water-content systems like moist soils,
aqueous solutions, vegetable or animal tissues and foods that freeze uncovered or without
an impervious and tight packaging material. The rate of both phenomena (solidification and
sublimation) is determined both by material characteristics (mainly composition, structure,
shape and size) and cooling conditions (temperature, humidity and rate of the media that
surrounds the phase change material). The sublimation process, in spite of its magnitude
being much less than that of freezing process, determines fundamental features of the final
quality for foods and influences on the structure and utility of frozen tissues. Modelling of
these simultaneous processes is very difficult owed to the coupling of the heat and mass
transfer balances, the existence of two moving phase change fronts that advance with very
different rates and to the involved physical properties which are, in most cases, variable
with temperature and water content.
When high water-content materials like foods, tissues, gels, soils or water solutions of
inorganic or organic substances, held in open, permeable or untightly-sealed containers are
refrigerated to below their initial solidification temperature, two simultaneous physical
phenomena take place:
Liquid water solidifies (freeze), and
Surface ice sublimates.
For the description of the freezing process, the material can be divided into three zones:
unfrozen, frozen and dehydrated. Freezing begins from the refrigerated surface/s, at a
temperature (T
if
) lower than that of pure water, due to the presence of dissolved materials,
and continues along an equilibrium line. Simultaneously, ice sublimation begins at the
frozen surface and a dehydration front penetrates the material, whose rate of advance is
again determined by all the abovementioned characteristics of the material and
environmental conditions. Normally this rate is much lower than that of the freezing front.
A complete mathematical model has to solve both, the heat transfer (freezing) and the mass
transfer (weight loss) simultaneously (Campaone et al., 2005a & b).
www.intechopen.com

473
Phase change is accounted for in the following way:
Solidification (freezing) as a freezing front (x = s
f
(t)) located in the point where material
temperature reaches the initial freezing temperature (T
if
), determined by material
composition. For temperatures lower than T
if
(the zone nearer to the refrigerated
surface) the amount of ice formed is determined by an equilibrium line (ice content vs
temperature and water content) specific to the material.
On the dehydration front (x = s
d
(t)) we impose Stefan-like conditions for temperature
distribution and vapor concentration.
We consider a semi-infinite material with characteristics similar to a very dilute gel (whose
properties can be supposed equal to those of pure water). The system has initial uniform
temperature equal to T
if
and uncovered flat surface which at time t=0 is exposed to the
surrounding medium (with constant temperature T
s
(lower than T
if
) and heat and mass
transfer coefficients h and K
m
). We assume that
s if
T T t T t
0
( ) , 0 < < > where T t
0
( ) is the
unknown sublimation temperature.
To calculate the evolution of temperature and water content in time, we will consider the
following free boundary problem: Find the temperatures ( )
d d
T T x t , = and ( )
f f
T T x t , = , the
concentrations ( )
va va
C C x t , = , the free boundaries ( )
d d
s s t = and ( )
f f
s s t = and the
temperature ( ) T T t
0 0
= at the sublimation front ( )
d
x s t = which must satisfy the following:
Differential equations at the dehydrated region:
( )
d
d d
d p d d
T T
C k x s t t
t x
2
2
, 0 , 0

= < < >

(270)
( )
va va
ef d
C C
D x s t t
t x
2
2
, 0 , 0

= < < >

(271)
Differential equations at the frozen region:
( ) ( )
f
f f
f p f d f
T T
C k s t x s t t
t x
2
2
, , 0

= < < >

(272)
Free boundary conditions at the sublimation front ( )
d
x s t = :

( ) ( ) ( ) ( ) ( )
d d f d
T s t t T s t t T t t
0
, , , 0 = = > (273)

( ) ( ) ( ) ( )
( )
f d d d
f d s s d
T s t t T s t t
k k L m s t t
x x
, ,
, 0

= >

(274)

( ) ( )
( )
va d
ef s d
C s t t
D m s t
x
,
=
(275)
( ) ( )
sat
va d
g g
c
b
T t MP T
C s t t Ma
R T t R T t
0
0 0
exp
( ) ( )
,
( ) ( )

= = (276)
www.intechopen.com

474
where ( ) ( )
va d
C s t t , is the equilibrium vapor concentration at ( ) T t
0
and the saturation
pressure
sat
P T ( ) is evaluated according to (Fennema & Berny, 1974).
Free boundary conditions at the freezing front ( )
f
x s t = :
( ) ( )
f f if
T s t t T t , , 0 = > (277)

( ) ( )
( )
f f
f f f f
T s t t
k m L s t t
x
,
, 0
= >
(278)
The convective boundary conditions at the fixed interphase x 0 = :

( )
( ) ( )
d
d d s
T t
k h T t T t
x
0,
0, , 0
= >
(279)

( )
( ) ( )
va
ef m va a
C t
D K C t C t
x
0,
0, , 0
= >
(280)
The initial conditions at t 0 = :

f d
s s (0) (0) 0 = = (281)

if
T T = for x 0 . (282)
We will solve the system (270) - (282) by using the quasi-steady method. In general, it is a
good approximation when the Stefan number tends to zero, i.e. when the latent heat of the
material is high with respect to the heat capacity of the solid material. This approximation
has often been used when modelling the freezing of high-water content materials.
Theorem 23. (Olguin et al., 2008)
The temperatures
f d
T T , and the concentration
va
C are given by the following expressions:

( ) ( ) ( ) ( )
d d
T x t A t B t x x s t t , , 0 , 0 = + < < > (283)
( ) ( ) ( ) ( )
a
v d
C x t D t E t x x s t t , , 0 , 0 = + < < > (284)
( ) ( ) ( ) ( ) ( )
f d f
T x t F t G t x s t x s t t , , , 0 = + < < > (285)

where ( ) A t B t D t , ( ), ( ) and ( ) E t as a function of ( ) T t
0
and ( )
d
s t , as well as ( ) F t and ( ) G t
as a function of ( ) T t
0
, ( )
d
s t and ( )
f
s t , given by the following expressions:
( )
( ) ( )
( )
( )
( )
( )
s d
s d
d
d d
d d
h
T t T s t
T t T k h
A t B t
h h
k
s t s t
k k
0
0
,
1 1
+
= =
+ +
(286)
www.intechopen.com

475
( ) ( )
m
a d a
ef g g
m
m m
ef
d d
ef ef
c c
b b
T t T t K
C s t Ma Ma C
D R T t R T t
K
D t E t
K K
D
s t s t
D D
0 0
0 0
exp exp
( ) ( )
( )
( ) ( )
,
1 ( ) 1 ( )

+
= =
+ +
(287)
( )
( ) ( ) ( )
( ) ( )
( )
( )
( ) ( )
f if d if
f d f d
T t s t T s t T T t
F t G t
s t s t s t s t
0 0
,

= =

(288)
and we obtain the following system of two ordinary differential equations and one algebraic
equation for ( )
d
s t , ( )
f
s t and ( ) T t
0
given by:

( ) ( )
( )
( ) ( )
( ) ( ) ( ) ( )
0
0
0
1 1 1
1 1
1
exp
( )
( )
1 ( )

+ + + + =

+

=

+

s s
f d f d
f if d if f if f d f
f f d
f
s m
a
m
if f g
d
ef
T t h T T h
s t hs t s t hs t
k T k T k T k k k
h
c
s t s t s t
b
k
T t L K
M a C
K
T k R T t
s t
D
(289)

( )
( ) ( )
( )
( ) ( )
( ) ( ) ( ) ( )
s s
f d f d
f if d f if if f d f
f if
d
s s
f f d
d
T t
h T T h
s t hs t s t hs t
k T k k T T k k k
k T
s t
m L h
s t s t s t
k
0
1 1 1
1 1
1

+ + + +

=

+

(290)

( )
( )
( ) ( )
o
f if if
f
f f f d
T t
k T T
s t
m L s t s t
1
=
(291)
( ) ( )
f d
s s 0 0 0 = = . (292)
Remark 21.
There exist some approximate or explicit solutions for some other free boundary problems
for the heat-diffusion equation, e.g.: model for a single nutrient uptake by a growing root
system by using a moving boundary approach; explicit estimate for the asymptotic
behavior of the solution of the porous media equation with absorption (reaction-diffusion
processes of a gas inside a chemical reactor); penetration of solvents in polymers; filtration
of water through oil in a porous medium; the Wen model for an isothermal mono-
catalytic diffusion-reaction process of a gas with a solid. The solid is chemically attacked
from its surface with a quick and irreversible reaction and, at the same time, a free
boundary begins, etc.
www.intechopen.com

476
4. Conclusion
We have given a review on explicit and approximated solutions for heat and mass transfer
problems in which a free or moving interface is involved. We have also showed some new
recent problems for heat and mass transfer in which a free or moving interface is also
involved.
5. Acknowledgements
This paper was partially sponsored by the project PIP No. 0460 of CONICET - UA (Rosario,
Argentina), and Grant FA9550-10-1-0023.
6. References
Alexandrov, D.V. & Malygin, A.P. (2006). Self-similar solidification of an alloy from a cooled
boundary. Int. J. Heat Mass Transfer, 49, 763-769.
Alexiades, V. & Solomon, A. D. (1993). Mathematical modelling of melting and freezing processes,
Hemisphere - Taylor & Francis, Washington.
Arderius, J.C.; Lara, M. & Tarzia, D.A. (1996). Experimental-numerical determination of
thermal coefficients through a phase-change process. Int. Comm. Heat Mass Transfer,
23, 745-754.
Bancora, A.B. & Tarzia, D.A. (1985). On the Neumann solution for the two-phase Stefan
problem including the density jump at the free boundary. Latin Amer. J. Heat Mass
Transfer, 9, 215-222.
Bankoff, S.G. (1964). Heat conduction of diffusion with change of phase. Advances in
Chemical Engineering, 5, 75-150.
Barenblatt, G.I. (1952). On some unsteady motions of a liquid or a gas in a porous medium.
Prink. Mat. Mekh, 16, 67-78 (in Russian).
Bell, G.E. (1978). A refinement of the heat balance integral method applied to a melting
problem. Int. J. Heat Mass Transfer, 21, 1357-1361.
Brillouin, M. (1930/31). Sur quelques problmes non rsolus de physique mathmatique
classique : Propagation de la fusion. Annales de lInst. H. Poincar, 1, 285-308.
Briozzo, A.C. & Tarzia, D.A. (1998). The explicit solution of a free boundary problem for a
nonlinear absorption model of mixed saturated-unsaturated flow. Advances in
Water Resources, 21, 713-721.
Briozzo, A.C. ; Natale, M.F. & Tarzia, D.A. (1999). Determination of unknown thermal
coefficients for Storm's type materials through a phase-change process.
International Journal of Non-Linear Mechanics, 34, 329-340.
Briozzo, A.C.; Natale, M.F. & Tarzia, D.A. (2004). An explicit solution for a two-phase Stefan
problem with a similarity exponential heat sources. MAT Serie A, 8, 11-19.
Briozzo, A.C.; Natale, M.F. & Tarzia, D.A. (2007a). Explicit solutions for a two-phase
unidimensional Lam-Clapeyron-Stefan problem with source terms in both phases.
J. Math. Anal. Appl., 329, 145-162.
Briozzo, A.C.; Natale, M.F. & Tarzia, D.A. (2007b). Existence of an exact solution for a one-
phase Stefan problem with nonlinear thermal coefficients from Tirskiis method.
Nonlinear Analysis, 67, 1989-1998.
www.intechopen.com

477
Briozzo, A.C.; Natale, M.F. & Tarzia, D.A. (2010). The Stefan problem with temperature-
dependent thermal conductivity and a convective term with a convective condition
at the fixed face. Comm. Pure and Applied Analysis, 9, 1209-1220.
Briozzo, A.C. & Tarzia, D.A. (2002). An explicit solution for the instantaneous two-phase
Stefan problem with nonlinear thermal coefficients. IMA J. Applied Math, 67, 249-
261.
Broadbridge, P. (1990). Solution of a nonlinear absorption model of mixed saturated-
unsaturated flow. Water Resour. Research, 26, 2435-2443.
Broadbridge, P. & White, I. (1988). Constant rate rainfall infiltration: A versatile nonlinear
model, 1, Analytic solutions. Water Resources Research, 24, 145-154.
Campaone, L.A., Salvadori, V.O. & Mascheroni, R.H. (2005a). Food freezing with
simultaneous surface dehydration. Approximate prediction of freezing time, Int. J.
Heat Mass Transfer, 48 (6), 1205-1213.
Campaone, L.A., Salvadori, V.O. & Mascheroni, R.H. (2005b). Food freezing with
simultaneous surface dehydration. Approximate prediction of weight loss during
freezing and storage, Int. J. Heat Mass Transfer, 48 (6), 1195-1204.
Cannon, J.R. (1963). Determination of an unknown coefficient in a parabolic differential
equation. Duke Math. J., 30, 313-323.
Cannon, J.R. (1964). Determination of certain parameters in heat conduction problems. J.
Math. Anal. Appl., 8, 188-201.
Cannon, J.R. (1984). The one-dimensional heat equation, Addison-Wesley, Menlo Park,
California.
Charach, Ch. & Rubisntein, I. (1992). Pressure-temperature effects in planar Stefan problems
with density change, J. Appl. Phys., 71, 1128-1137.
Carslaw, H.S. & Jaeger, J.C. (1959). Conduction of heat in solids, Clarendon Press, Oxford.
Castro, E.F., Gonzo, E.E. & Gottifredi, J.C. Gottifredi (1987). Thermodynamic of the
absorption of hydrocarbon vapor in poliethylene films, J. Mem. Sci., 31, 235-248.
Cho, S.H. (1975). An exact solution of the coupled phase change problem in a porous
medium. Int. J. Heat Mass Transfer, 18, 1139-1142.
Cho, S.H. & Sunderland, J.E. (1974). Phase change problems with temperature-dependent
thermal conductivity. J. Heat Transfer, 96C, 214-217.
Cirelli, M. B. & Tarzia, D.A. (2010), Explicit solutions for the Rubinstein binary-alloy
solidification problem with a heat flux or a convective condition at the fixed face,
In preparation.
Crank, J. (1975). The mathematics of diffusion, Clarendon Press, Oxford.
Crank, J. (1984). Free and moving boundary problems, Clarendon Press, Oxford.
Crank, J. & Gupta, R.S. (1972). A moving boundary problem arising from the diffusion of
oxygen in absorbing tissue.J. Inst. Math. Appl., 10, 19-33.
Destefanis, H.A.; Erdmann, E.; Tarzia, D.A. & Villa, L.T. (1993). A free boundary model
applied to the estimation of the diffusion coefficient in a gas-solid system.
International Communications Heat Mass Transfer, 20, 103-110.
Duvaut, G. (1976). Problmes frontire libre en thorie des milieux continus, Rapport de
Recherche No. 185, LABORIA IRIA, Rocquencourt.
Elliott, C.M. & Ockendon, J.R. (1982). Weak and variational methods for moving boundary
problems, Research Notes in Math., No. 59, Pitman, London.
www.intechopen.com

478
Fasano, A. (1974). Alcune osservazioni su una classe di problemi a contorno libero per
lequazione del calore, Le Matematiche, 29, 397-411.
Fasano, A. (1987). Las zonas pastosas en el problema de Stefan, Cuadern. Inst. Mat. Beppo Levi,
13, 1-60.
Fasano, A. (2005). Mathematical models of some diffusive processes with free boundaries, MAT
Serie A, 11, 1-128.
Fasano, A.; Guan, Y.; Primicerio, M. & Rubinstein, I. (1993). Thawing in saturated porous
media, Meccanica, 28, 103-109.
Fasano, A. & Primicerio, M. (1984). Freezing in porous media: a review of mathematical
model, In: Applications of Mathematics in Technology, V. Boffi & H. Neunzert (Eds.),
288-311, Teubner.
Fasano, A. & Primicerio, M. (1986). Nonlinear diffusion problems, Lecture Notes in Mathematics
No. 1224, Springer Verlag, Berlin.
Fasano, A.; Primicerio, M. & Tarzia, D.A. (1999). Similarity solutions in Stefan-like problems.
Math Models and Methods in Appl Sci, 9, 1-10.
Fennema, O, Berny, L.A. (1974). Equilibrium vapour pressure and water activity of food at
subfreezing temperature. In Proceedings of IV International Congress of Food Science
and Technology, 2, 27-35.
Fi, F. & Han, X. (2007). An one dimensional two-phase free boundary problem in an
angular domain. Nonlinear Anal. RWA, 8, 959-979
Friedman, A. (1964). Partial differential equations of parabolic type, Prentice Hall, Englewood
Cliffs, N. J.
Froment, G.F. & Bischoff, K.B. (1979). Chemical reactor analysis and design, J. Wiley, New York.
Garguichevich, G.G. & Sanziel, M.C. (1984). Una introduccin general a la resolucin
aproximada del problema de Stefan unidimensional. Cuadern. Inst. Mat Beppo Levi,
11, 167-177.
Garguichevich, G.G.; Sanziel, M.C. & Tarzia, D.A. (1985). Comparison of approximate
methods for the determination of thermal coefficients through a phase-change
problems. Int. Comm. Heat Mass Transfer, 12, 451-464.
Gonzlez, A.M.; Reginato, J.C. & Tarzia, D.A. (2008). A free-boundary model for
anaerobiosis in saturated soil aggregates. Soil Science, 173, 758-767.
Gonzlez, A.M. & Tarzia, D.A. (1996). Determination of unknown coefficient of a semi-
infinite material through a simple mushy zone model for the two-phase Stefan
problem. International Journal of Engineering Science, 34, 799-817.
Goodman, T.R. (1958). The heat balance integral and its application to problems involving a
change of phase. Trans. of the ASME, 80, 335-342.
Gupta, L.N. (1974). An approximate solution to the generalized Stefans problem in a porous
medium, Int. J. Heat Mass Transfer, 17, 313321.
Gupta, S.C. (2003). The classical Stefan problem. Basic concepts, modelling and analysis, Elsevier,
Amsterdam.
Gupta, S.C.; Sanziel, M.C. & Tarzia, D.A. (1997). A similarity solution for the binary alloy
solidification problem with a simple mushy zone model. Anales de la Academia
Nacional de Ciencias Exactas, Fsicas y Naturales, 49, 75-82.
www.intechopen.com

479
Guzman, J.A. (1982). Application of the heat transfer solution for the case of one-
dimensional solidification with negligible latent heat. Lat. Amer. J. Heat Mass
Transfer, 6, 41-51.
Hill, J.M. (1987). One-dimensional Stefan Problems : an introduction, Longman, Essex .
Howison, S.D. (1988). Similarity solutions to the Stefan problem and the binary alloy
problem. IMA J. Appl. Math., 40, 147-161.
Jones, Jr. B. F. (1962). The determination of a coefficient in a parabolic differential equation,
Part I : Existence and uniqueness. J. Math. Mech., 11, 907-918.
Jones, Jr. B. F. (1963). Various methods for finding unknown coefficients in parabolic
differential equation. Comm. Pure Appl. Math., 16, 33-44.
Knight, J. & Philip, J.R. (1974). Exact solutions in nonlinear diffusion. J. Eng. Math, 8, 219-227.
Lam, G. & Clapeyron, B.P. (1831). Memoire sur la solidification par refroidissement dun
globe liquide. Annales Chimie Physique, 47, 250-256.
Levenspiel, O. (1962). Chemical reaction engineering, J. Wiley, New York.
Liapis, A. I.; Lipscomb, G.G.; Crosser, O.K. & Liapis-Tsiroyianni, E. (1982). A model of
oxygen diffusion in absorbing tissue. Mathematical Modelling, 3, 83-92.
Lin, S. (1981). An exact solution of the sublimation problem in a porous medium. J. Heat
Transfer, 103, 165-168.
Lin, S. (1982a). An exact solution of the sublimation problem in a porous medium, Part II
With an unknow temperature and vapor concentration at the moving sublimation
front. J. Heat Transfer, 104, 808-811.
Lin, S. (1982b). An exact solution of the desublimation problem in a porous medium. Int. J.
Heat Mass Transfer, 25, 625-630.
Lombardi, A. & Tarzia, D.A. (2001). Similarity solutions for thawing processes with a heat
flux on the fixed boundary. Meccanica, 36, 251-264.
Luikov, A.V. (1964), Heat and mass transfer in capillary-porous bodies, Adv. HeatTransfer, 1,
123184.
Luikov, A.V. (1966). Heat and mass transfer in capillary-porous bodies, Pergamon Press, Oxford.
Luikov, A. V. (1968). Analytical heat diffusion theory, Academic Press, New York.
Luikov, A.V. (1975). Systems of differential equations of heat and mass transfer in capillary porous
bodies, Int. J. Heat Mass Transfer, 18, 114.
Luikov, A.V. (1978). Heat and mass transfer, MIR Publishers, Moscow.
Lunardini, V.J. (1981). Heat transfer in cold climates, Van Nostrand, New York.
Lunardini, V.J. (1991). Heat transfer with freezing and thawing, Elsevier, Amsterdam.
Menaldi, J.L. & Tarzia, D.A. (1993). Generalized Lam-Clapeyron solution for a one-phase
source Stefan problem. Comput. Appl. Math., 12, 123-142.
Mikhailov, M.D. (1975). Exact solution of temperature and mosture distributions in a porous
half-space with moving evaporation front. Int. J. Heat Mass Transfer, 18, 797-804.
Mikhailov, M.D. (1976). Exact solution for freezing of humid porous half-space. Int. J. Heat
Mass Transfer, 19, 651-655.
Muehlbauer, J.C. & Sunderland, J.E. (1965). Heat conduction with freezing or melting. Appl.
Mech. Reviews, 18, 951-959.
Nakano, Y. (1990). Quasi-steady problems in freezing soils, Cold Reg. Sci. Techn., 17, 207-226.
www.intechopen.com

480
Natale, M.F.; Santillan Marcus, E.A. & Tarzia, D.A. (2010). Explicit solutions for one
dimensional two-phase free boundary problems with either shrinkage or
expansion. Nonlinear Analysis Series B: Real World Applications, 11, 1946-1952.
Natale, M.F. & Tarzia, D.A. (2000). Explicit solutions to the two-phase Stefan problem for
Storm's type materials. Journal of Physics A: Mathematical and General, 33, 395-404.
Natale, M.F. & Tarzia, D.A. (2006). Explicit solutions for a one-phase Stefan problem with
temperature-dependent thermal conductivity. Bollettino dellUnione Matematica
Italiana, 9B, 79-99.
Olguin, M.C.; Salvadori, V.O.; Mascheroni, R.H. & Tarzia, D.A. (2008). An analytical
solution for the coupled heat and mass transfer during the freezing of high-water
content materials. International Journal of Heat and Mass Transfer, 51, 4379-4391.
O'Neill, K. & Miller, R.D. (1985). Explanation of a rigid ice model of frost heave, Water
Resour. Res., 21, 281-296.
Pattle, R.E. (1959). Diffusion from an instantaneous point source with concentration
dependent coefficient. Quart. J. Mech. Appl. Math, 12, 407-409.
Pekeris, C.L. & Slichter, L.B. (1939). Problem of ice formation. J. Appl. Physics, 10, 135-137.
Philip, J.R. (1957). The theory of infiltration, 4, Sorptivity and algebraic infiltration
equations. Soil Sci., 84, 257-264.
Philip, J.R. (1958). The theory of infiltration, 6, Effect of water depth over soil. Soil Sci., 85,
278-286.
Primicerio, M. (1981). Problemi di diffusione a frontiera libera. Bollettino Un. Mat. Italiana,
18A, 11-68.
Primicerio, M. & Gianni, R. (1989). La filtracin en medios porosos. Cuadern. Inst. Mat. Beppo
Levi, 18 1-92.
Quilghini, D. (1967). Sul problema inverso di quello di Stefan. Riv. Mat. Univ. Parma, 8, 131-
142.
Reginato, J.C.; Palumbo, M.C.; Moreno, I.S.; Bernardo, I.Ch. & Tarzia, D.A. (2000). Modeling
nutrient uptake using a moving boundary approach. Comparison with the Barber-
Cushman model. Soil Science Society of America Journal, 64, 1363-1367.
Reginato, J.C. & Tarzia, D.A. (1993). The balance integral method applied to root growth of
crops. International Journal of Engineering Science, 31, 61-70.
Reginato, J.C.; Tarzia, D.A. & Dzioba, M.A. (1993). A model for root growth with root
competition. Plant and Soil, 15, 185-196.
Rogers, C. & Broadbridge, P. (1988). On a nonlinear moving boundary problem with
heterogeneity: application of a reciprocal transformation. ZAMP, 39, 122-128.
Rogers, C. (1986). On a class of moving boundary problems in non-linear heat conduction:
Application of a Backlund transformation. Int. J. Nonlinear Mech., 21, 249-256.
Rubinstein, L.I. (1971). The Stefan problem, Translations of Mathematical Monographs, 27,
Amer. Math. Soc., Providence, R.I.
Salva, N.N. & Tarzia, D.A. (2010). Simultaneous determination of unknown coefficients
through a phase-change process with temperature-dependent thermal
conductivity, Submitted.
Santillan Marcus, E.A.; Natale, M.F. & Tarzia, D.A. (2008). Simultaneous determination of
two unknown thermal coefficients of a semi-infinite porous material through a
www.intechopen.com

481
desublimation moving boundary problem with coupled heat and moisture flows.
JP Journal of Heat and Mass Transfer, 2, 73-116
Santillan Marcus, E.A. & Tarzia, D.A. (2000a). Soluciones explcitas para el problema de
desublimacin en un semi-espacio hmedo poroso con una condicin de flujo de
calor. Informacin Tecnolgica, 11, 19-24.
Santillan Marcus, E.A. & Tarzia, D.A. (2000b). Explicit solution for freezing of humid porous
half-space with heat flux condition. International Journal of Engineering Science, 38,
1651-1665.
Santillan Marcus, E.A. & Tarzia, D.A. (2003). Exact solutions for drying with coupled phase
change in a porous medium with a heat flux condition on the surface. Computational
and Applied Mathematics, 22, 293-311.
Santillan Marcus, E.A. & Tarzia, D.A. (2007). Determination of one unknown thermal
coefficient of a semi-infinite porous material through a desublimation problem
with coupled heat and moisture flows. JP Journal of Heat and Mass Transfer, 1, 251-
270.
Sanziel, M.C. & Tarzia, D.A. (1989). Necessary and sufficient condition to obtain n phases in
a one-dimensional medium with a flux condition on the fixed face. Math. Notae, 33,
25-32.
Schatz, A. (1969). Free boundary problems of Stephan type with prescribed flux. J. Math.
Anal. Appl. 28, 569-580.
Sherman, B. (1971). Limiting behaviour in some Stefan problems as the latent goes to zero.
SIAM J. Appl. Math., 20, 319-327.
Solomon, A.D. (1979). An easily computable solution to a two-phase Stefan problem. Solar
Energy, 23, 525-528.
Solomon, A.D.; Wilson, D.G. & Alexiades, V. (1982). A mushy zone model with an exact
solution. Letters Heat Mass Transfer, 9, 319-324.
Solomon, A.D.; Wilson, D.G. & Alexiades, V. (1983). Explicit solutions to change problems.
Quart. Appl. Math., 41, 237-243.
Stampella, M.B. & Tarzia, D.A. (1989). Determination of one or two unknown thermal
coefficients of a semi-infinite material through a two-phase Stefan problem. Int. J.
Eng. Science, 27, 1407-1419.
Stefan, J. (1889a). Ueber einige Probleme der Theorie der Waermeleitung. Zitzungsberichte
der Kaiserlichen Akademie der Wissenschaften, Mathematisch-naturwissenschaftliche
classe, 98, 473-484.
Stefan, J. (1889b). Ueber die Theorie der Eisbildung, insbesondere ueber die Eisbildung im
Polarmeere. Zitzungsberichte der Kaiserlichen Akademie der Wissenschaften,
Mathematisch-naturwissenschaftliche Classe, 98, 965-983.
Stefan, J. (1890). Ueber die Theorie der Eisbildung. Monatshefte fur Mathematik und Physik, 1,
1-6.
Szekely, J.; Evans, J.W. & Sohn, H.Y. (1976). Gas-solid reactions, Academic Press, New York.
Talamucci, F. (1997). Analysis of the coupled heat-mass transport in freezing porous media,
Survey in Math. in Industry, 7, 93-139.
Talamucci, F. (1998). Freezing processes in saturated soils, Math. Models Meth. Appl. Sci., 8,
107-138.
www.intechopen.com

482
Tarzia, D.A. (1981). An inequality for the coefficient of the free boundary ( ) 2 s t t = of
the Neumann solution for the two-phase Stefan problem. Quart. Appl. Math., 39,
491-497.
Tarzia, D.A. (1982). Determination of the unknown coefficients in the Lame-Clapeyron
problem (or one-phase Stefan problem). Adv. Appl. Math., 3, 74-82.
Tarzia, D.A. (1983). Simultaneous determination of two unknown thermal coefficients
through an inverse one-phase Lam-Clapeyron (Stefan) problem with an
overspecified condition on the fixed face. Int. J. Heat Mass Transfer, 26, 1151-1158.
Tarzia, D.A. (1984a). A new variant for the simultaneous calculation of some thermal
coefficients of a semi-infinite material through a phase-change problem with an
over-condition on the fixed face. Latin Amer. J. Heat Mass Transfer, 8, 227-235.
Tarzia, D.A. (1984b). Soluciones exactas del problema de Stefan unidimensional. Cuadern.
Inst. Mat. Beppo Levi, 12, 5-36.
Tarzia, D.A. (1987). Determination of unknown thermal coefficients of a semi-infinite
material for the one-phase Lam-Clapeyron (Stefan) problem through the
Solomon-Wilson-Alexiades mushy zone model. Int. Comm. Heat. Mass Transfer,
14, 219-228.
Tarzia, D.A. (1990a). A variant of the heat balance integral method and a new proof of the
exponentially fast asymptotic behaviour of the solutions in heat conduction
problems with absorption. Int. J. Engineering Science, 28, 1253-1259.
Tarzia, D.A. (1990b). Neumann-like solution for the two-phase Stefan problem with a simple
mushy zone model. Computational and Applied Mathematics (ex Matemtica Aplicada
Computacional), 9, 201-211.
Tarzia, D.A. (1991a). Sobre una nueva variante para el clculo simultneo de coeficientes
trmicos de un material semi-infinito a travs de un problema directo o inverso de
Stefan a dos fases. Math. Notae, 35, 25-41.
Tarzia, D.A. (1991b). Six free boundary problems for the heat-diffusion equation. Revista Un.
Mat. Argentina, 37, 294-309.
Tarzia, D.A. (1993). Approximate and analytic methods to solve some parabolic free
boundary problems, In Free Boundary Problems Involving Solids (FBP90), J.M.
Chadam - H. Rasmussen (Eds.), Pitman Research Notes in Mathematics, Vol. 281,
pp. 190-196, Longman, Essex.
Tarzia, D.A. (1998). The determination of unknown thermal coefficients through phase
change process with temperature-dependent thermal conductivity. International
Communications in Heat and Mass Transfer, 25, 139-147.
Tarzia, D.A. (2000a). A bibliography on moving-free boundary problems for the heat-
diffusion equation. The Stefan and related problems, MAT - Serie A, 2, 1-297. See:
http://web.austral.edu.ar/descargas/facultad-cienciasEmpresariales/mat/Tarzia-
MAT-SerieA-2(2000).pdf
Tarzia, D.A. (2000b). Explicit solution to some free boundary problems for the heat equation,
In International Conference on Free Boundary Problems (FBP99), N. Kenmochi (Ed.),
GAKUTO International Series, Mathematical Sciences and Applications, Vol. 13,
pp. 385-398, Gottotosho, Tokyo.
Tarzia, D.A. (2004). An explicit solution for a two-phase unidimensional Stefan problem
with a convective boundary condition at the fixed face. MAT Serie A, 8, 21-27.
www.intechopen.com

483
Tarzia, D.A. (2007). Exact solution for a Stefan problem with convective boundary condition
and density jump. PAMM Proceedings on Applied Mathematics and Mechanics, 7,
1040307-1040308.
Tarzia, D.A. & Turner, C.V. (1992). A note on the existence of a waiting time for a two-phase
Stefan problem. Quart. Appl. Math., 50, 1-10.
Tarzia, D.A. & Turner, C.V. (1999). Estimation of the occurance of the phase-change process
in spherical coordinates. International Communications in Heat and Mass Transfer, 26,
559-568.
Tarzia, D.A. & Villa, L.T. (1991). Mathematical considerations on the heat transfer with
phase change with negligible latent heat. Latin American Applied Research, 21, 1-6.
Tayler, A.B. (1986). Mathematical models in applied mechanics, Clarendon Press, Oxford .
Tien, R.H. & Geiger, G.E. (1967). A heat-transfer analysis of the solidification of a binary
eutectic system. J. Heat Transfer, 89C, 230-234.
Tien, R.H. & Koump, V. (1970). Effect of density change on the solidification of alloys. J. Heat
Transfer, 92C, 11-16.
Tirskii, G.A. (1959). Two exact solutions of Stefans nonlinear problems. Soviet Physics Dokl.,
4, 288-292.
Tritscher, P. & Broadbridge, P. (1994). A similarity solution of a multiphase Stefan problem
incorporating general non-linear heat conduction. Int. J. Heat Mass Transfer, 37,
2113-2121.
Tsubaki, T. & Boley, B.A. (1977). One-dimensional solidification of binary mixtures. Mech.
Res. Comm., 4, 115-122.
Vzquez, J.L. (2007). The porous medium equation. Mathematical theory, Clarendon Press,
Oxford.
Villa, L.T. (1987). La ecuacin de la difusin y su aplicacin a problemas de frontera libre.
Cuadern. Inst. Mat. Beppo Levi, 14, 99-116.
Voller, V.R. (2006). A similarity solution for solidification of an under-cooled binary alloy.
Int. J. Heat Mass Transfer, 49, 1981-1985.
Voller, V.R. (2008a). A numerical method for the Rubinstein binary- alloy problem in the
presence of an under-cooled liquid. Int. J. Heat Mass Transfer, 51, 696-706.
Voller, V.R. (2008b). An enthalpy method for modeling dendritic growth in a binary alloy.
Int. J. Heat Mass Transfer, 51, 823-834.
Warrick, A.W. & Broadbridge, P. (1992). Sorptivity and macroscopic capillary length
relationships. Water. Resour. Researh, 28, 427-431.
Weber, H. (1901). Die partiellen Differential-Gleinchungen der Mathematischen Physik, nach
Riemann's Vorlesungen, t. II, Braunwschweig.
Weiner, J.H. (1955). Transient heat conduction in multi-phase media. British J. Appl. Physics,
6, 361-363.
White, R.E. (1985). The binary alloy problem: existence, uniqueness, and numerical
approximation. SIAM J. Numer. Anal., 22, 205-244.
Wilson, D.G. (1978). Existence and uniqueness for similarity solutions of one dimensional
multi-phase Stefan problems. SIAM J. Appl. Math., 35, 135-147.
Wilson, D.G. (1982). Lagrangian coordinates for moving boundary problems. SIAM J. Appl.
Math, 42, 1195-1201.
www.intechopen.com

484
Wilson, D.G. & Solomon, A.D. (1986). A Stefan-type problem with void formation and its
explicit solution. IMA J. Appl. Math., 37, 67-76.
Wilson, D.G. & Solomon, A.D & Alexiades, V. (1982). A shortcoming of the explicit solution
for the binary alloy solidification problem. Letters Heat Mass Transfer, 9, 421-428.
Zeldovich, Y.B. & Kompaneets, A.S. (1950). Towards a theory of heat conduction with thermal
conductivity depending on the temperature. In collection of papers dedicated to 70
th

anniversary of A.F. Ioffe, Izd. Akad. Nauk. SSSR, Moscow, 61-72 (in Russian).

www.intechopen.com
Edited by Prof. Mohamed El-Amin
ISBN 978-953-307-333-0
Hard cover, 626 pages
Publisher InTech
Published online 21, February, 2011
Published in print edition February, 2011
InTech Europe
University Campus STeP Ri
Slavka Krautzeka 83/A
51000 Rijeka, Croatia
Phone: +385 (51) 770 447
Fax: +385 (51) 686 166
www.intechopen.com
InTech China
Unit 405, Office Block, Hotel Equatorial Shanghai
No.65, Yan An Road (West), Shanghai, 200040, China
Phone: +86-21-62489820
Fax: +86-21-62489821
This book introduces a number of selected advanced topics in mass transfer phenomenon and covers its
theoretical, numerical, modeling and experimental aspects. The 26 chapters of this book are divided into five
parts. The first is devoted to the study of some problems of mass transfer in microchannels, turbulence, waves
and plasma, while chapters regarding mass transfer with hydro-, magnetohydro- and electro- dynamics are
collected in the second part. The third part deals with mass transfer in food, such as rice, cheese, fruits and
vegetables, and the fourth focuses on mass transfer in some large-scale applications such as geomorphologic
studies. The last part introduces several issues of combined heat and mass transfer phenomena. The book
can be considered as a rich reference for researchers and engineers working in the field of mass transfer and
its related topics.
How to reference
In order to correctly reference this scholarly work, feel free to copy and paste the following:
Domingo Alberto Tarzia (2011). Explicit and Approximated Solutions for Heat and Mass Transfer Problems
with a Moving Interface, Advanced Topics in Mass Transfer, Prof. Mohamed El-Amin (Ed.), ISBN: 978-953-
307-333-0, InTech, Available from: http://www.intechopen.com/books/advanced-topics-in-mass-
transfer/explicit-and-approximated-solutions-for-heat-and-mass-transfer-problems-with-a-moving-interface

Explicit and Approximated Solutions For Heat and Mass Transfer Problems With A Moving Interface

Uploaded by

Document Informationclick to expand document information

Copyright:

Available Formats

Explicit and Approximated Solutions For Heat and Mass Transfer Problems With A Moving Interface

Uploaded by

Document Information

Original Description:

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Explicit and Approximated Solutions For Heat and Mass Transfer Problems With A Moving Interface

Uploaded by

Copyright:

Available Formats

20

Explicit and Approximated Solutions

= = is the diffusion coefficient and 0 > is the unique solution to the

= + < < > (16)

< << . (22)

< < >

= = + > > . (86)

< < . We have the

< < >

< < < >

> > >

+ = > > (125)

= < < > , (153)

= < < > , (156)

> satisfies the following equation:

> depends on the

) must satisfy the following conditions:

= + < < >

= < < > , (266)

You might also like