Geoscience Reference
In-Depth Information
It is convenient to introduce new dimensionless functions g 1 D .B 1 C 2B 2 /=B 0
and g 2 D .B 1 B 2 /=B 0 . Summing and subtracting Eq. ( 7.85 ) and ( 7.86 ), and
rearranging Eq. ( 7.87 ), leads to
@ t g 1 D m @ r g 1 C
2 @ r V C 2 V
r
;
2@ r g 1
r
(7.89)
@ t g 2 D m @ r g 2 C
r 2
2@ r g 2
r
6g 2
V
r ;
C @ r V
(7.90)
6g 2
r D 0:
@ r .g 1 C 2g 2 / C
(7.91)
Such a form of the equations are preferable because either of Eqs. ( 7.89 )or( 7.90 )
includes only a single unknown function. Note that Eq. ( 7.89 ) can be transformed
to a classical 1D diffusion equation via changing the unknown function g 1 .r;t/ by
the function g 0 1 .r;t/ D g 1 .r;t/=r. The solution of this equation at zero initial and
boundary conditions are known (Surkov 1989a )
Z
Z
t
dt 0
r 0
g 3 r;r 0 ;r 0 r 0 @ r 0 V C 2V dr 0 ;
2
1=2 r
g 1 .r;t/ D
(7.92)
0
0
where V D V .r 0 ;t 0 / is the mass velocity of the medium. Here we made use of the
following abbreviations:
g 3 r;r 0 ;r 0 D exp
!
exp
! ;
.r r 0 / 2
r 0
.r C r 0 / 2
r 0
(7.93)
and
r 0 D 2 m t t 0 1=2 :
(7.94)
Substituting Eq. ( 7.92 )forg 1 into Eq. ( 7.91 ) we come to a differential equation
of the first order with respect to the function g 2 . The solution of this equation takes
the form
Z
r
r 03 @ r 0 g 1 r 0 ;t dr 0 :
1
2r 3
g 2 .r;t/ D
(7.95)
0
So the functions g 1 and g 2 given by Eq. ( 7.92 ) and Eq. ( 7.95 ) are the solutions of
Eq. ( 7.89 ) and ( 7.91 ).
 
Search WWH ::




Custom Search