Environmental Engineering Reference
In-Depth Information
5.7.3 Method of Continuation for Spherically-Symmetric Problems
We now apply the method of continuation to solve
r 2
u t
τ 0 +
A 2 1
r 2
u
u tt
=
+
f
(
r
,
t
) ,
0
<
r
< + ,
0
<
t
,
r
r
(5.93)
u r
(
0
,
t
)=
0
,
u
(
r
,
0
)= ϕ (
r
) ,
u t (
r
,
0
)= ψ (
r
) .
By the solution structure theorem, we can first seek the solution of
u t τ 0 +
A 2
u tt =
Δ
u
,
0
<
r
< + ,
0
<
t
,
u r (
0
,
t
)=
0
,
(5.94)
u
(
r
,
0
)=
0
,
u t (
r
,
0
)= ψ (
r
) .
should be an even
function of r . To obtain the solution of PDS (5.94), consider an auxiliary problem
Based on the given boundary condition u r (
0
,
t
)=
0, the u
(
r
,
t
)
u t
τ 0 +
A 2
u tt =
Δ
u
, − <
r
< + ,
0
<
t
,
(5.95)
ψ (
r
) ,
0
r
,
u
(
r
,
0
)=
0
,
u t (
r
,
0
)= Ψ (
r
)=
ψ (
r
) ,
r
<
0
,
where
Ψ (
r
)
comes from an even prolongation of
ψ (
r
)
.
By function transformation v
(
r
,
t
)=
ru
(
r
,
t
)
, PDS (5.95) is transformed to
v t
τ
A 2 v rr , − <
0 +
v tt =
r
< + ,
0
<
t
,
r
(5.96)
ψ (
r
) ,
0
r
,
v
(
r
,
0
)=
0
,
v t (
r
,
0
)=
ψ (
) ,
<
.
r
r
r
0
Its solution is available in Eq. (5.30). Thus the solution of PDS (5.94) is,
u
(
r
,
t
)=
W ψ (
r
,
t
)
t
I 0 b
2
τ 0 r + At
r
2 Ar e
1
2
2
=
(
At
)
(
r
ρ )
ρψ ( ρ )
d
ρ ,
(5.97)
At
where
0.
Finally, the solution of PDS (5.93) reads, by the solution structure theorem,
ψ ( ρ )
should be
ψ ( ρ )
if
ρ <
1
W ϕ (
t
τ 0 +
u
=
r
,
t
)+
W ψ (
r
,
t
)+
W f τ (
r
,
t
τ )
d
τ ,
(5.98)
t
0
1 2 a τ 0 ,
where f τ =
f
(
r
, τ )
, b
=
ϕ (
r
)= ϕ (
r
)
, f
(
r
,
t
)=
f
(
r
,
t
)
.
 
Search WWH ::




Custom Search