Using the dimensionless expressions in Eqs. (9.28)–(9.38) and the following approximate expressions at interfaces z=0 and z=h3,
∂p′1∂z∣∣z=h3=p1−p'1|z=h3h1/2and∂p′2∂z∣∣z=0=p12∣∣z=0−p2h2/2
(H.1)
Eqs. (9.8)–(9.21) becomes
ω1∂pD1∂tD−w1ΔDpD1+2kvD1(pD1−p′D1∣∣zD=1)=0
(H.2)
ω2∂pD2∂tD−w2ΔDpD2+2kvD2(pD2−p′D2∣∣zD=0)=0
(H.3)
pDj→0whentD→0,j=1,2
(H.4)
pDj→0whenrD→∞,j=1,2
(H.5)
wj(rD∂pDj∂rD)rD=1=−qDj,j=1,2
(H.6)
pwDj=(pDj−sjrD∂pDj∂rD)rD=1,j=1,2
(H.7)
qLD=f(pwD1,pwD2)
(H.8)
or
qLD=αD(pwD1−pwD2)
(H.9)
qD1=1−qLD−CD1dpwD1dtD
(H.10)
qD2=qLD−CD2dpwD2dtD
(H.11)
kvD1hD1∂p′D1∂zD∣∣zD=1=kvD3∂p′D3∂zD∣∣zD=1
(H.12)
kvD2hD2∂p′D2∂zD∣∣zD=0=kvD3∂p′D3∂zD∣∣zD=0
(H.13)
where ΔD=1rD(rD∂∂rD) is a dimensionless differential operator. Eqs. (9.23)–(9.26) for the middle layer becomes
ω3∂p′D3∂tD−kvD3∂2p′D3∂z2D=0
(H.14)
p′D3→0fortD→0
(H.15)
p′D1∣∣zD=1=p′D3∣∣zD=1
(H.16)
p′D2∣∣zD=0=p′D3∣∣zD=0
(H.17)
After Laplace transformation, Eqs. (H.14), (H.15) become
d2p¯′D3dz2D−υ2p¯′D3=0
(H.18)
where
υ=ω3σ/kvD3−−−−−−−−√
(9.58)
The solution of Eq. (H.18) under boundary conditions Eqs. (H.16), (H.17) is
p¯′D3=[p¯′D1|zD=1sinhυzD+p¯′D2|zD=0sinhυ(1−zD)]/sinhυ
(H.19)
Substituting Eq. (H.19) into the Laplace transforms of conditions Eqs. (H.12), (H.13) and noticing Eq. (H.1), we have
[2kvD1+υkvD3cothυ]p′¯¯¯D1∣∣zD=1−υkvD3sinhυp′¯¯¯D2∣∣zD=0=2kvD1p¯D1
(H.20)
−υkvD3sinhυp¯′D1|zD=1+(2kvD2+υkvD3cothυ)p¯′D2|zD=0=2kvD2p¯D1
(H.21)
Solve Eqs. (H.20), (H.21) to obtain the following:
p¯′D1|zD=1=1G0[(1+kvD3υcothυ2kvD2)p¯D1+kvD3υcschυ2kvD1p¯D2]
(H.22)
p¯′D2|zD=0=1G0[kvD3υcschυ2kvD2p¯D1+(1+kvD3υcothυ2kvD1)p¯D2]
(H.23)
where
G0=1+kvD3υcothυ2(1kvD1+1kvD2)+ω3σkvD34kvD1kvD2
(9.57)
The Laplace transform of Eqs. (H.2), (H.3) is
ω1σp¯D1−w1ΔDp¯D1+2kvD1(p¯D1−p¯′D1|zD=1)=0
(H.24)
ω2σp¯D2−w2ΔDp¯D2+2kvD2(p¯D2−p¯′D2|zD=0)=0
(H.25)
Substituting Eqs. (H.22), (H.23) into Eqs. (H.24), (H.25) results in the following:
−w1ΔDp¯D1+G1p¯D1−Gp¯D2=0
(H.26)
−w2ΔDp¯D2−Gp¯D1+G2p¯D2=0
(H.27)
where
G1=ω1σ+(υcothυ+ω3σ2kvD2)kvD3G0
(9.54)
G2=ω2σ+(υcothυ+ω3σ2kvD1)kvD3G0
(9.55)
G=υkvD3G0sinhυ
(9.56)
Considering Eq. (H.5), the solution of Eqs. (H.26), (H.27) has the form
p¯D1=AK0(λrD)
(H.28)
p¯D2=DK0(λrD)
(H.29)
where K0 is the modified Bessel function of zero order, and A and D are arbitrary functions of Laplace variable σ.
Substituting Eqs. (H.28), (H.29) into Eqs. (H.26), (H.27), we have
(G1−λ2w1)A−GD=0
(H.30)
−GA+(G2−λ2w2)D=0
(H.31)
The eigenequation is
λ4−(G1w1+G2w2)λ2+G1G2−G2w1w2=0
(H.32)
The two eigenvalues are
λ1,2=[12(G1w1+G2w2)±14(G1w1+G2w2)2+G2−G1G2w1w2−−−−−−−−−−−−−−−−−−−−√]12
(9.51)
From Eqs. (H.31), (H.30), we find
D=β2Aforeigenvalueλ1
and
A=β1Dforeigenvalueλ2
where
β1=G/(G1−λ22w1)
(9.52)
β2=G/(G2−λ21w2)
(9.53)
Using
λ21+λ22=G1w1+G2w2
we have
β2=−β1w1/w2
(H.33)
So the solution is
p¯D1=AK0(λ1rD)λ1K1(λ1)+Dβ1K0(λ2rD)λ2K1(λ2)
(H.34)
and
p¯D2=Aβ2K0(λ1rD)λ1K1(λ1)+DK0(λ2rD)λ2K1(λ2)
(H.35)
where K1 is the modified Bessel function of the first order. The Laplace transforms of Eqs. (H.6), (H.7) are
wj(rDdp¯DjdrD)rD=1=−q¯Dj,j=1,2
(H.36)
p¯wDj=(p¯Dj−sjrDdp¯DjdrD)rD=1,j=1,2
(H.37)
Substituting Eqs. (H.34), (H.35) into Eqs. (H.36), (H.37) results in the following:
w1[A(σ)+D(σ)β1]=q¯D1
(9.44)
w2[A(σ)β2+D(σ)]=q¯D2
(9.45)
A(σ)[K(λ1)+s1]+D(σ)β1[K(λ2)+s1]=p¯wD1
(9.46)
A(σ)β2[K(λ1)+s2]+D(σ)[K(λ2)+s2]=p¯wD2
(9.47)
where function K(x) is defined by
K(x)=K0(x)xK1(x)≈ln2γx
(9.59)
The Laplace transform of Eqs. (H.8)–(H.11) gives
q¯D1=1σ−q¯LD−σCD1p¯wD1
(9.48)
q¯D2=q¯LD−σCD2p¯wD2
(9.49)
or
q¯LD=f¯(pwD1,pwD2)
(9.50a)
q¯LD=αD(p¯wD1−p¯wD2)
(9.50b)
Eqs. (9.44)–(9.50) provides the solution of the problem.
H.1 When Leakage Function f Is Known
Solving A(σ) and D(σ) from Eqs. (9.44), (9.45) and using Eqs. (9.48)–(9.50a), (H.33), we have
A(σ)=[1σ−f¯−σCD1p¯wD1+β2(f¯−σCD2p¯wD2)]/(1−β1β2)w1
(H.38)
D(σ)=[f¯−σCD2p¯wD2+β1(1σ−f¯−σCD1p¯wD1)]/(1−β1β2)w2
(H.39)
Substituting Eqs. (H.38), (H.39) into Eqs. (9.46), (9.47) and solving p¯wD1 for p¯wD2 and, we have:
p¯wD1=M4M5−M2M6M1M4−M2M3
(9.60)
and
p¯wD2=M1M6−M3M5M1M4−M2M3
(9.61)
where
M1=w1w2(1−β1β2)+σCD1w2{K(λ1)+s1−β1β2[K(λ2)+s1]}
(9.62)
M2=w1β1σCD2[K(λ2)−K(λ1)]
(9.63)
M3=w1β1σCD1[K(λ2)−K(λ1)]
(9.64)
M4=w1w2(1−β1β2)+σCD2w1{K(λ2)+s2−β1β2[K(λ1)+s2]}
(9.65)
M5=w2[K(λ1)+s1][1σ+(β2−1)f¯]+w1β1[K(λ2)+s1][f¯+β1(1σ−f¯)]
(9.66)
M6=w2β2[K(λ1)+s2][1σ+(β2−1)f¯]+w1[K(λ2)+s2][f¯+β1(1σ−f¯)]
(9.67)
H.2 When Leakage Function f=αD(pwD1,pwD2)
Substituting Eqs. (9.50b), (H.38), (H.39) into Eqs. (9.46), (9.47) and solving for p¯wD1 and p¯wD2, we have
p¯wD1=L5L4−L2L6σ(L1L4−L2L3)
(9.68)
and
p¯wD2=L1L6−L3L5σ(L1L4−L2L3)
(9.69)
where
L1=w1w2(1−β1β2)+w2[K(λ1)+s1][αD(1−β2)+σCD1]−β1w1[K(λ2)+s1][αD(1−β1)−σβ1CD1]
(9.70)
L2=−w2[K(λ1)+s1][αD(1−β2)−σβ2CD2]+β1w1[K(λ2)+s1][αD(1−β1)+σCD2]
(9.71)
L3=w2β2[K(λ1)+s2][αD(1−β2)+σCD1]−w1[K(λ2)+s2][αD(1−β1)−σβ1CD1]
(9.72)
L4=w1w2(1−β1β2)−w2β2[K(λ1)+s2][αD(1−β2)−σβ2CD2]+w1[K(λ2)+s2][αD(1−β1)+σCD2]
(9.73)
L5=w2{K(λ1)+s1−β1β2[K(λ2)+s1]}
(9.74)
L6=w1β1[K(λ2)−K(λ1)]
(9.75)
H.3 When Leakage Function f Is Unknown
Adding Eqs. (9.44), (9.45) and using Eqs. (9.48), (9.49), (H.33), we have
w1(1−β1)A(σ)+w2(1−β2)D(σ)=1σ−σ(CD1p¯wD1+CD2p¯wD2)
(9.81)
Solving Eqs. (9.47), (9.81), we have
A={p¯wD2w2(1−β2)−[1σ−σ(CD1p¯wD1+CD2p¯wD2)][K(λ2)+s2]}/L7
(H.40)
D={β2[1σ−σ(CD1p¯wD1+CD2p¯wD2)][K(λ1)+s2]−p¯wD2(1−β1)w1}/L7
(H.41)
where
L7=(1−β2)w2β2[K(λ1)+s2]−w1(1−β1)[K(λ2)+s2]
(H.42)
Substituting Eqs. (H.40), (H.41) into Eq. (9.46) and solving for p¯wD1, we have
p¯wD1=[K(λ1)+s1]{(1−β2)w2p¯wD2−(1σ−σCD2p¯wD2)[K(λ2)+s2]}/L8+β1[K(λ2)+s1]{β2[K(λ1)+s2](1σ−σCD2p¯wD2)−w1(1−β1)p¯wD2}/L8
(9.82)
where
L8=(1−β2)w2β2[K(λ1)+s2]−w1(1−β1)[K(λ2)+s2]−σCD1{[K(λ1)+s1][K(λ2)+s2]−β1β2[K(λ2)+s1][K(λ1)+s2]}
(9.83)
From Eqs. (9.45), (9.49), (9.50a), we have
f¯=w2[A(σ)β2+D(σ)]+σCD2p¯wD2
(9.84)
where A(σ) and D(σ) are determined by Eqs. (H.40), (H.41).