Appendix H

Solution of the Problem (Chapter 9)

Using the dimensionless expressions in Eqs. (9.28)(9.38) and the following approximate expressions at interfaces z=0 and z=h3,

p1zz=h3=p1p1z=h3h1/2andp2zz=0=p21z=0p2h2/2

si1_e  (H.1)

Eqs. (9.8)(9.21) becomes

ω1pD1tDw1ΔDpD1+2kvD1pD1pD1zD=1=0

si2_e  (H.2)

ω2pD2tDw2ΔDpD2+2kvD2pD2pD2zD=0=0

si3_e  (H.3)

pDj0whentD0,j=1,2

si4_e  (H.4)

pDj0whenrD,j=1,2

si5_e  (H.5)

wjrDpDjrDrD=1=qDj,j=1,2

si6_e  (H.6)

pwDj=pDjsjrDpDjrDrD=1,j=1,2

si7_e  (H.7)

qLD=fpwD1pwD2

si8_e  (H.8)

or

qLD=αDpwD1pwD2

si9_e  (H.9)

qD1=1qLDCD1dpwD1dtD

si10_e  (H.10)

qD2=qLDCD2dpwD2dtD

si11_e  (H.11)

kvD1hD1pD1zDzD=1=kvD3pD3zDzD=1

si12_e  (H.12)

kvD2hD2pD2zDzD=0=kvD3pD3zDzD=0

si13_e  (H.13)

where ΔD=1rDrDrDsi14_e is a dimensionless differential operator. Eqs. (9.23)(9.26) for the middle layer becomes

ω3pD3tDkvD32pD3zD2=0

si15_e  (H.14)

pD30fortD0

si16_e  (H.15)

pD1zD=1=pD3zD=1

si17_e  (H.16)

pD2zD=0=pD3zD=0

si18_e  (H.17)

After Laplace transformation, Eqs. (H.14), (H.15) become

d2p¯D3dzD2υ2p¯D3=0

si19_e  (H.18)

where

υ=ω3σ/kvD3

si20_e  (9.58)

The solution of Eq. (H.18) under boundary conditions Eqs. (H.16), (H.17) is

p¯D3=p¯D1zD=1sinhυzD+p¯D2zD=0sinhυ1zD/sinhυ

si21_e  (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¯D1zD=1υkvD3sinhυp¯D2zD=0=2kvD1p¯D1

si22_e  (H.20)

υkvD3sinhυp¯D1zD=1+2kvD2+υkvD3cothυp¯D2zD=0=2kvD2p¯D1

si23_e  (H.21)

Solve Eqs. (H.20), (H.21) to obtain the following:

p¯D1zD=1=1G01+kvD3υcothυ2kvD2p¯D1+kvD3υcschυ2kvD1p¯D2

si24_e  (H.22)

p¯D2zD=0=1G0kvD3υcschυ2kvD2p¯D1+1+kvD3υcothυ2kvD1p¯D2

si25_e  (H.23)

where

G0=1+kvD3υcothυ21kvD1+1kvD2+ω3σkvD34kvD1kvD2

si26_e  (9.57)

The Laplace transform of Eqs. (H.2), (H.3) is

ω1σp¯D1w1ΔDp¯D1+2kvD1p¯D1p¯D1zD=1=0

si27_e  (H.24)

ω2σp¯D2w2ΔDp¯D2+2kvD2p¯D2p¯D2zD=0=0

si28_e  (H.25)

Substituting Eqs. (H.22), (H.23) into Eqs. (H.24), (H.25) results in the following:

w1ΔDp¯D1+G1p¯D1Gp¯D2=0

si29_e  (H.26)

w2ΔDp¯D2Gp¯D1+G2p¯D2=0

si30_e  (H.27)

where

G1=ω1σ+υcothυ+ω3σ2kvD2kvD3G0

si31_e  (9.54)

G2=ω2σ+υcothυ+ω3σ2kvD1kvD3G0

si32_e  (9.55)

G=υkvD3G0sinhυ

si33_e  (9.56)

Considering Eq. (H.5), the solution of Eqs. (H.26), (H.27) has the form

p¯D1=AK0λrD

si34_e  (H.28)

p¯D2=DK0λrD

si35_e  (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λ2w1AGD=0

si36_e  (H.30)

GA+G2λ2w2D=0

si37_e  (H.31)

The eigenequation is

λ4G1w1+G2w2λ2+G1G2G2w1w2=0

si38_e  (H.32)

The two eigenvalues are

λ1,2=12G1w1+G2w2±14G1w1+G2w22+G2G1G2w1w212

si39_e  (9.51)

From Eqs. (H.31), (H.30), we find

D=β2Aforeigenvalueλ1

si40_e

and

A=β1Dforeigenvalueλ2

si41_e

where

β1=G/G1λ22w1

si42_e  (9.52)

β2=G/G2λ12w2

si43_e  (9.53)

Using

λ12+λ22=G1w1+G2w2

si44_e

we have

β2=β1w1/w2

si45_e  (H.33)

So the solution is

p¯D1=AK0λ1rDλ1K1λ1+Dβ1K0λ2rDλ2K1λ2

si46_e  (H.34)

and

p¯D2=Aβ2K0λ1rDλ1K1λ1+DK0λ2rDλ2K1λ2

si47_e  (H.35)

where K1 is the modified Bessel function of the first order. The Laplace transforms of Eqs. (H.6), (H.7) are

wjrDdp¯DjdrDrD=1=q¯Dj,j=1,2

si48_e  (H.36)

p¯wDj=p¯DjsjrDdp¯DjdrDrD=1,j=1,2

si49_e  (H.37)

Substituting Eqs. (H.34), (H.35) into Eqs. (H.36), (H.37) results in the following:

w1Aσ+Dσβ1=q¯D1

si50_e  (9.44)

w2Aσβ2+Dσ=q¯D2

si51_e  (9.45)

AσKλ1+s1+Dσβ1Kλ2+s1=p¯wD1

si52_e  (9.46)

Aσβ2Kλ1+s2+DσKλ2+s2=p¯wD2

si53_e  (9.47)

where function K(x) is defined by

Kx=K0xxK1xln2γx

si54_e  (9.59)

The Laplace transform of Eqs. (H.8)(H.11) gives

q¯D1=1σq¯LDσCD1p¯wD1

si55_e  (9.48)

q¯D2=q¯LDσCD2p¯wD2

si56_e  (9.49)

or

q¯LD=f¯pwD1pwD2

si57_e  (9.50a)

q¯LD=αDp¯wD1p¯wD2

si58_e  (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+β2f¯σCD2p¯wD2/1β1β2w1

si59_e  (H.38)

Dσ=f¯σCD2p¯wD2+β11σf¯σCD1p¯wD1/1β1β2w2

si60_e  (H.39)

Substituting Eqs. (H.38), (H.39) into Eqs. (9.46), (9.47) and solving p¯wD1si61_e for p¯wD2si62_e and, we have:

p¯wD1=M4M5M2M6M1M4M2M3

si63_e  (9.60)

and

p¯wD2=M1M6M3M5M1M4M2M3

si64_e  (9.61)

where

M1=w1w21β1β2+σCD1w2Kλ1+s1β1β2Kλ2+s1

si65_e  (9.62)

M2=w1β1σCD2Kλ2Kλ1

si66_e  (9.63)

M3=w1β1σCD1Kλ2Kλ1

si67_e  (9.64)

M4=w1w21β1β2+σCD2w1Kλ2+s2β1β2Kλ1+s2

si68_e  (9.65)

M5=w2Kλ1+s11σ+β21f¯+w1β1Kλ2+s1f¯+β11σf¯

si69_e  (9.66)

M6=w2β2Kλ1+s21σ+β21f¯+w1Kλ2+s2f¯+β11σf¯

si70_e  (9.67)

H.2 When Leakage Function f=αDpwD1pwD2si71_e

Substituting Eqs. (9.50b), (H.38), (H.39) into Eqs. (9.46), (9.47) and solving for p¯wD1si61_e and p¯wD2si62_e, we have

p¯wD1=L5L4L2L6σL1L4L2L3

si74_e  (9.68)

and

p¯wD2=L1L6L3L5σL1L4L2L3

si75_e  (9.69)

where

L1=w1w21β1β2+w2Kλ1+s1αD1β2+σCD1β1w1Kλ2+s1αD1β1σβ1CD1

si76_e  (9.70)

L2=w2Kλ1+s1αD1β2σβ2CD2+β1w1Kλ2+s1αD1β1+σCD2

si77_e  (9.71)

L3=w2β2Kλ1+s2αD1β2+σCD1w1Kλ2+s2αD1β1σβ1CD1

si78_e  (9.72)

L4=w1w21β1β2w2β2Kλ1+s2αD1β2σβ2CD2+w1Kλ2+s2αD1β1+σCD2

si79_e  (9.73)

L5=w2Kλ1+s1β1β2Kλ2+s1

si80_e  (9.74)

L6=w1β1Kλ2Kλ1

si81_e  (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

w11β1Aσ+w21β2Dσ=1σσCD1p¯wD1+CD2p¯wD2

si82_e  (9.81)

Solving Eqs. (9.47), (9.81), we have

A=p¯wD2w21β21σσCD1p¯wD1+CD2p¯wD2Kλ2+s2/L7

si83_e  (H.40)

D=β21σσCD1p¯wD1+CD2p¯wD2Kλ1+s2p¯wD21β1w1/L7

si84_e  (H.41)

where

L7=1β2w2β2Kλ1+s2w11β1Kλ2+s2

si85_e  (H.42)

Substituting Eqs. (H.40), (H.41) into Eq. (9.46) and solving for p¯wD1si61_e, we have

p¯wD1=Kλ1+s11β2w2p¯wD21σσCD2p¯wD2Kλ2+s2/L8+β1Kλ2+s1β2Kλ1+s21σσCD2p¯wD2w11β1p¯wD2/L8

si87_e  (9.82)

where

L8=1β2w2β2Kλ1+s2w11β1Kλ2+s2σCD1Kλ1+s1Kλ2+s2β1β2Kλ2+s1Kλ1+s2

si88_e  (9.83)

From Eqs. (9.45), (9.49), (9.50a), we have

f¯=w2Aσβ2+Dσ+σCD2p¯wD2

si89_e  (9.84)

where A(σ) and D(σ) are determined by Eqs. (H.40), (H.41).

..................Content has been hidden....................

You can't read the all page of ebook, please click here login for view all page.
Reset
3.15.25.32