Mathcad 11 CASA_06-27.mcd Appendix 1 γ s 76972.857 kg m 2 s 2 = pcf 157.087464 kg m 2 s 2 = ρ s γ s g := ρ s 7849.047053 kg m 3 = A 8 := A 8 = m 1 ρ 0 a ⋅ A := m 1 38.09371 kg m 2 = h a 62.752 := h 0.004857 m = m 2 ρ s h ⋅ := m 2 38.124515 kg m 2 = m 2 m 1 - m 1 0.000809 = m 3 ρ s π a h + ( ) 2 ⋅ π a 2 ⋅ - ⋅ 2 π ⋅ a 1 2 h ⋅ + ⋅ := m 3 38.124515 kg m 2 = m 3 m 1 - m 1 0.000809 = R 1 12.478 := E K ρ 0 ρ s 1 ν 2 - ( ) ⋅ ⋅ 12.468436 = Skalak, R. (1956), Ref. 4 "An extension of the theory of water hammer." Transactions of the ASME 78, 105-116. (PhD Thesis, Columbia University, New York, USA, 1954, Ref. 2) Table 1 g 32.174049 ft s 2 = g 9.80665 m s 2 = a 1 ft ⋅ := a 0.3048 m = c 5000 ft s ⋅ := c 1524 m s = ρ 0 1.94 slug ft 3 ⋅ := ρ 0 999.834908 kg m 3 = K ρ 0 c 2 ⋅ := K 2.322193 10 9 × Pa = E 30 10 6 ⋅ psi ⋅ := E 2.068427 10 11 × Pa = ν 0.3 := ν 0.3 = pcf lbf ft 3 := γ s 490 pcf ⋅ := Mathcad 11 CASA_06-27.mcd Appendix 2 There is something wrong with Skalak's equations [74] and [76]: the -c 0 2 term in the last factor must be deleted in [74] and [76] misses a square root. c 2 c 0 0.981061 = c 2 c 3.464192 = c 1 c e 0.999004 = c 1 c 0.643647 = c e c 0.644288 = c 0 c 3.531067 = Check with Table 1 Eq. (76) c 2 5279.429013 m s = c 2 c 2A ⋅ R ⋅ R + R 2 1 ν 2 - ( ) ⋅ + 2A ⋅ R ⋅ R + R 2 1 ν 2 - ( ) ⋅ + 2 4R 2 ⋅ 1 ν 2 - ( ) ⋅ 2A ⋅ R + ( ) ⋅ - + 2 2A ⋅ R + ( ) ⋅ 0.5 ⋅ := Eq. (76) c 1 980.917501 m s = c 1 c 2A ⋅ R ⋅ R + R 2 1 ν 2 - ( ) ⋅ + 2A ⋅ R ⋅ R + R 2 1 ν 2 - ( ) ⋅ + 2 4R 2 ⋅ 1 ν 2 - ( ) ⋅ 2A ⋅ R + ( ) ⋅ - - 2 2A ⋅ R + ( ) ⋅ 0.5 ⋅ := Eq. (1) c e 981.895414 m s = c e c 1 2K ⋅ a ⋅ Eh ⋅ + := Korteweg - Joukowsky A 7.993536 = A ρ 0 a ⋅ m 2 := R R 2 := R 2 12.468436 = R 2 c 0 2 c 2 := Nomenclature c 0 5381.346495 m s = c 0 Eh ⋅ m 2 1 ν 2 - ( ) ⋅ := Wave speeds Mathcad 11 CASA_06-27.mcd Appendix 3 c t 5133.477378 m s = c t c 0 0.953939 = With FSI: γ2 c2 F c2 t + ν 2 ρ f ρ t ⋅ D ee ⋅ c2 F ⋅ + := γ2 2.883457 10 7 × m 2 s 2 = λ2 1 1 2 γ2 γ2 2 4 c2 F ⋅ c2 t ⋅ - - ⋅ := λ2 1 9.622 10 5 × m 2 s 2 = λ2 3 1 2 γ2 γ2 2 4 c2 F ⋅ c2 t ⋅ - + ⋅ := λ2 3 2.787 10 7 × m 2 s 2 = λ 1 λ2 1 := λ 1 980.917501 m s = c F c t ⋅ 5.17868431352064 10 6 × m 2 s 2 = λ 3 λ2 3 := λ 3 5279.429013 m s = λ 1 λ 3 ⋅ 5.17868431352065 10 6 × m 2 s 2 = λ 3 c t 1.028431 = Check with Table 1 λ 1 c 0.643647 = λ 1 c e 0.999004 = λ 3 c 3.464192 = λ 3 c 0 0.981061 = ρ f ρ 0 := ρ f 999.834908 kg m 3 = Tijsseling D 2a ⋅ := D 0.6096 m = ee h := ee 0.004857 m = ρ t ρ s := ρ t 7849.047053 kg m 3 = Longitudinal wave speeds Classical: c2 F K ρ f 1 KD ⋅ E ee ⋅ 1 ν 2 - ( ) ⋅ + 1 - ⋅ := c F c2 F := c F 1008.806299 m s = c F c e 1.027407 = c2 f K ρ f 1 KD ⋅ E ee ⋅ + 1 - ⋅ := c f c2 f := c f 981.895414 m s = c f c e 1 = c2 t E ρ t := c t c2 t := Mathcad 11 CASA_06-27.mcd Appendix 4 d 2 11.477473 m 3 s = d 1 2.925721 m 3 s = d 2 ca 2 ⋅ A 4 + ( ) c 2 c 5 c 2 c 3 1 R + ( ) ⋅ - c 2 c R ⋅ + 16 - c 2 c 2 ⋅ 2A ⋅ R + ( ) ⋅ 8R ⋅ 2A ⋅ 1 + ( ) ⋅ + 8R 2 ⋅ 1 ν 2 - ( ) + ⋅ := Eq. (78) d 1 ca 2 ⋅ A 4 + ( ) c 1 c 5 c 1 c 3 1 R + ( ) ⋅ - c 1 c R ⋅ + 16 - c 1 c 2 ⋅ 2A ⋅ R + ( ) ⋅ 8R ⋅ 2A ⋅ 1 + ( ) ⋅ + 8R 2 ⋅ 1 ν 2 - ( ) + ⋅ := Eq. (78) Nominal length of the wave front according to Tijsseling d 2old 11.477844 m 3 s = d 1old 2.925652 m 3 s = d 2old ca 2 ⋅ A 4 + ( ) c 2 c 5 c 2 c 3 1 R + ( ) ⋅ - c 2 c R ⋅ + c 2 c h 2 12 a 2 ⋅ ⋅ 1 2R 2 ⋅ ν ⋅ + ( ) ⋅ - 16 - c 2 c 2 ⋅ 2A ⋅ R + ( ) ⋅ 8R ⋅ 2A ⋅ 1 + ( ) ⋅ + 8R 2 ⋅ 1 ν 2 - ( ) + ⋅ := Eq. (78) d 1old ca 2 ⋅ A 4 + ( ) c 1 c 5 c 1 c 3 1 R + ( ) ⋅ - c 1 c R ⋅ + c 1 c h 2 12 a 2 ⋅ ⋅ 1 2R 2 ⋅ ν ⋅ + ( ) ⋅ - 16 - c 1 c 2 ⋅ 2A ⋅ R + ( ) ⋅ 8R ⋅ 2A ⋅ 1 + ( ) ⋅ + 8R 2 ⋅ 1 ν 2 - ( ) + ⋅ := Eq. (78) Nominal length of the wave front according to Skalak λ 3 c 2 1 = λ 1 c 1 1 = Check Skalak - Tijsseling
13
Embed
Skalak, R. (1956), Ref. 4 Wave speeds ( ) Nomenclature 11 CASA_06-27.mcd Appendix 1 γs 76972.857 kg m2s2 = pcf 157.087464 kg m2s2 ρs γs g:= ρs 7849.047053 kg m3 A := 8 A 8= …
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Mathcad 11 CASA_06-27.mcd Appendix 1
γs 76972.857kg
m2
s2
= pcf 157.087464kg
m2
s2
=
ρs
γs
g:= ρs 7849.047053
kg
m3
=
A 8:= A 8=
m1
ρ0 a⋅
A:= m1 38.09371
kg
m2
=
ha
62.752:= h 0.004857 m=
m2 ρs h⋅:= m2 38.124515kg
m2
=m2 m1−
m1
0.000809=
m3
ρs π a h+( )2
⋅ π a2
⋅− ⋅
2 π⋅ a1
2h⋅+
⋅
:= m3 38.124515kg
m2
=m3 m1−
m1
0.000809=
R1 12.478:=E
K
ρ0
ρs 1 ν2
−( )⋅
⋅ 12.468436=
Skalak, R. (1956), Ref. 4"An extension of the theory of water hammer."Transactions of the ASME 78, 105-116.(PhD Thesis, Columbia University, New York, USA, 1954, Ref. 2)
Table 1
g 32.174049ft
s2
= g 9.80665m
s2
=
a 1 ft⋅:= a 0.3048m=
c 5000ft
s⋅:= c 1524
m
s=
ρ0 1.94slug
ft3
⋅:= ρ0 999.834908kg
m3
=
K ρ0 c2
⋅:= K 2.322193 109
× Pa=
E 30 106
⋅ psi⋅:= E 2.068427 1011
× Pa=
ν 0.3:= ν 0.3= pcflbf
ft3
:=
γs 490 pcf⋅:=
Mathcad 11 CASA_06-27.mcd Appendix 2
There is something wrong with Skalak's equations [74] and [76]: the -c02 term in
the last factor must be deleted in [74] and [76] misses a square root.
c2
c0
0.981061=c2
c3.464192=
c1
ce
0.999004=c1
c0.643647=
ce
c0.644288=
c0
c3.531067=
Check with Table 1
Eq. (76)c2 5279.429013m
s=
c2 c2 A⋅ R⋅ R+ R
21 ν
2−( )⋅+ 2 A⋅ R⋅ R+ R
21 ν
2−( )⋅+
2
4 R2
⋅ 1 ν2
−( )⋅ 2 A⋅ R+( )⋅−+
2 2 A⋅ R+( )⋅
0.5
⋅:=
Eq. (76)c1 980.917501m
s=
c1 c2 A⋅ R⋅ R+ R
21 ν
2−( )⋅+ 2 A⋅ R⋅ R+ R
21 ν
2−( )⋅+
2
4 R2
⋅ 1 ν2
−( )⋅ 2 A⋅ R+( )⋅−−
2 2 A⋅ R+( )⋅
0.5
⋅:=
Eq. (1)ce 981.895414m
s=ce
c
12 K⋅ a⋅
E h⋅+
:=
Korteweg - Joukowsky
A 7.993536=Aρ0 a⋅
m2
:=
R R2:=R2 12.468436=R2
c02
c2
:=
Nomenclaturec0 5381.346495m
s=c0
E h⋅
m2 1 ν2
−( )⋅
:=
Wave speeds
Mathcad 11 CASA_06-27.mcd Appendix 3
ct 5133.477378m
s=
ct
c0
0.953939=
With FSI:
γ2 c2F c2t+ ν2
ρf
ρt
⋅D
ee⋅ c2F⋅+:= γ2 2.883457 10
7×
m2
s2
=
λ211
2γ2 γ2
24 c2F⋅ c2t⋅−−
⋅:= λ21 9.622 10
5×
m2
s2
=
λ231
2γ2 γ2
24 c2F⋅ c2t⋅−+
⋅:= λ23 2.787 10
7×
m2
s2
=
λ1 λ21:= λ1 980.917501m
s= cF ct⋅ 5.17868431352064 10
6×
m2
s2
=
λ3 λ23:= λ3 5279.429013m
s= λ1 λ3⋅ 5.17868431352065 10
6×
m2
s2
=
λ3
ct
1.028431=Check with Table 1
λ1
c0.643647=
λ1
ce
0.999004=λ3
c3.464192=
λ3
c0
0.981061=
ρf ρ0:= ρf 999.834908kg
m3
= Tijsseling
D 2 a⋅:= D 0.6096 m=
ee h:= ee 0.004857m=
ρt ρs:= ρt 7849.047053kg
m3
=
Longitudinal wave speeds
Classical:
c2FK
ρf
1K D⋅
E ee⋅1 ν
2−( )⋅+
1−
⋅:= cF c2F:= cF 1008.806299m
s=
cF
ce
1.027407=
c2fK
ρf
1K D⋅
E ee⋅+
1−
⋅:= cf c2f:= cf 981.895414m
s=
cf
ce
1=
c2tE
ρt
:= ct c2t:=
Mathcad 11 CASA_06-27.mcd Appendix 4
d2 11.477473m
3
s=d1 2.925721
m3
s=
d2 c a2
⋅
A 4+( )c2
c
5c2
c
3
1 R+( )⋅−c2
c
R⋅+
16−c2
c
2
⋅ 2 A⋅ R+( )⋅ 8 R⋅ 2 A⋅ 1+( )⋅+ 8 R2
⋅ 1 ν2
−( )+
⋅:=Eq. (78)
d1 c a2
⋅
A 4+( )c1
c
5c1
c
3
1 R+( )⋅−c1
c
R⋅+
16−c1
c
2
⋅ 2 A⋅ R+( )⋅ 8 R⋅ 2 A⋅ 1+( )⋅+ 8 R2
⋅ 1 ν2
−( )+
⋅:=Eq. (78)
Nominal length of the wave front according to Tijsseling
d2old 11.477844m
3
s=d1old 2.925652
m3
s=
d2old c a2
⋅
A 4+( )c2
c
5c2
c
3
1 R+( )⋅−c2
c
R⋅+
c2
c
h2
12 a2
⋅
⋅ 1 2 R2
⋅ ν⋅+( )⋅−
16−c2
c
2
⋅ 2 A⋅ R+( )⋅ 8 R⋅ 2 A⋅ 1+( )⋅+ 8 R2
⋅ 1 ν2
−( )+
⋅:=Eq. (78)
d1old c a2
⋅
A 4+( )c1
c
5c1
c
3
1 R+( )⋅−c1
c
R⋅+
c1
c
h2
12 a2
⋅
⋅ 1 2 R2
⋅ ν⋅+( )⋅−
16−c1
c
2
⋅ 2 A⋅ R+( )⋅ 8 R⋅ 2 A⋅ 1+( )⋅+ 8 R2
⋅ 1 ν2
−( )+
⋅:=Eq. (78)
Nominal length of the wave front according to Skalak
λ3
c2
1=λ1
c1
1=
Check Skalak - Tijsseling
Mathcad 11 CASA_06-27.mcd Appendix 5
L2 5 s⋅( ) 15.669309 m=L2 5 s⋅( ) 51.408494 ft=
L2 1 s⋅( ) 9.163467 m=L2 1 s⋅( ) 30.06387 ft=
L1 5 s⋅( ) 9.935248 m=L1 5 s⋅( ) 32.595957 ft=
L1 1 s⋅( ) 5.810168 m=L1 1 s⋅( ) 19.062231 ft=
Check with Table 1
Eq. (64)L2 t( )3 π⋅
3d2 t⋅⋅
Γ1
3
sinπ
3
⋅
:=L1 t( )3 π⋅
3d1 t⋅⋅
Γ1
3
sinπ
3
⋅
:=
3 π⋅
Γ1
3
sinπ
3
⋅
4.062354=Γ1
3
2
3π
3
1
2
Γ2
3
⋅⋅→sinπ
3
1
23
1
2⋅→
c2
c
h2
12 a2
⋅
⋅ 1 2 R2
⋅ ν⋅+( )⋅
A 4+( )c2
c
5c2
c
3
1 R+( )⋅−c2
c
R⋅+
0.000032−=
c1
c
h2
12 a2
⋅
⋅ 1 2 R2
⋅ ν⋅+( )⋅
A 4+( )c1
c
5c1
c
3
1 R+( )⋅−c1
c
R⋅+
0.000024=
Error because of "wrong" (h/a)2 term in Eq. (78)
h
a
2
0.000254=R 12.468436=A 7.993536=
Mathcad 11 CASA_06-27.mcd Appendix 6
0 1 2 3 4 5 6 7 8 9 100
5
10
15
20
L1 t( )
L2 t( )
t
TOL 1012−
:=even function to be integrated
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 51
0
1
2
sin η η3
+( )η
η
Mathcad 11 CASA_06-27.mcd Appendix 7
0
9.3
ηsin η β η
3⋅+( )
η
⌠⌡
d 0.840919869378856=β 10:=
0
20
ηsin η β η
3⋅+( )
η
⌠⌡
d 1.11015525927976=β 1:=
0
43
ηsin η β η
3⋅+( )
η
⌠⌡
d 1.42315940650369=β1
10:=
0
93
ηsin η β η
3⋅+( )
η
⌠⌡
d 1.56376138708481=β1
100:=
0
201
ηsin η β η
3⋅+( )
η
⌠⌡
d 1.57083144591109=β1
1000:=
0
430
ηsin η β η
3⋅+( )
η
⌠⌡
d 1.57077652953821=β1
10000:=
0
907
ηsin η β η
3⋅+( )
η
⌠⌡
d 1.57076569930658=β1
100000:=
0
1843
ηsin η β η
3⋅+( )
η
⌠⌡
d 1.57082809089641=β1
1000000:=
π
21.570796=
Eq. (49)
0
6600
ηsin η β η
3⋅+( )
η
⌠⌡
d 1.57093027532399=β 0:=
Determination of the maximum possible upper boundaries of integrals:
Mathcad 11 CASA_06-27.mcd Appendix 8
0
20
ηsin η β η
3⋅+( )
η
⌠⌡
d 0.432399540407618=β 1−:=
0
43
ηsin η β η
3⋅+( )
η
⌠⌡
d 1.73656497765972=β1−
10:=
0
94
ηsin η β η
3⋅+( )
η
⌠⌡
d 1.71751324086010=β1−
100:=
0
202
ηsin η β η
3⋅+( )
η
⌠⌡
d 1.18218831662671=β1−
1000:=
0
441
ηsin η β η
3⋅+( )
η
⌠⌡
d 1.56865597358518=β1−
10000:=
0
962
ηsin η β η
3⋅+( )
η
⌠⌡
d 1.70193923364200=β1−
100000:=
0
2115
ηsin η β η
3⋅+( )
η
⌠⌡
d 1.62570169385300=β1−
1000000:=
π
21.570796=
0
6600
ηsin η β η
3⋅+( )
η
⌠⌡
d 1.57093027532399=β 0:=
0
2.0
ηsin η β η
3⋅+( )
η
⌠⌡
d 0.599018563818040=β 1000:=
0
4.3
ηsin η β η
3⋅+( )
η
⌠⌡
d 0.681123567409668=β 100:=
Mathcad 11 CASA_06-27.mcd Appendix 9
β 10−:=
0
9.3
ηsin η β η
3⋅+( )
η
⌠⌡
d 0.123095673983503−=
β 100−:=
0
4.3
ηsin η β η
3⋅+( )
η
⌠⌡
d 0.347949976012880−=
β 1000−:=
0
2.0
ηsin η β η
3⋅+( )
η
⌠⌡
d 0.444274408294637−=
Maximum possible upper boundaries b of integrals as function of β:
b β( ) 1843 β 0.000001≤if
907 0.000001 β< 0.00001≤if
430 0.00001 β< 0.0001≤if
201 0.0001 β< 0.001≤if
93 0.001 β< 0.01≤if
43 0.01 β< 0.1≤if
20 0.1 β< 1≤if
9.3 1 β< 10≤if
4.3 10 β< 100≤if
2.0 100 β< 1000≤if
0 otherwise
:=
Mathcad 11 CASA_06-27.mcd Appendix 10
Define range of n values for β:
n 1000:=
i 0 n..:= j i( )10
ni⋅ 7−:= β i( ) 10
j i( ):= β i( )
110 -7
1.02310 -7
1.04710 -7
1.07210 -7
1.09610 -7
1.12210 -7
1.14810 -7
1.17510 -7
1.20210 -7
1.2310 -7
1.25910 -7
1.28810 -7
1.31810 -7
1.34910 -7
1.3810 -7
1.41310 -7
= b β i( )( )1843
1843
1843
1843
1843
1843
1843
1843
1843
1843
1843
1843
1843
1843
1843
1843
=
i
0
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
= j i( )
-7
-6.99
-6.98
-6.97
-6.96
-6.95
-6.94
-6.93
-6.92
-6.91
-6.9
-6.89
-6.88
-6.87
-6.86
-6.85
=
Maximum possible upper boundaries b of integrals as function of β:
1 .107
1 .106
1 .105
1 .104
1 .103
0.01 0.1 1 10 100 1 .103
0
500
1000
1500
2000
b β i( )( )
β i( )
Mathcad 11 CASA_06-27.mcd Appendix 11
Eq. (49) for negative (f=wn) and positive (g=wp) values of β:
f(x) is the integral from 0 to minus b(x) OR f(x) is minus the integral from 0 to b(x)
f x( )
0
b x( )
ηsin η x η
3⋅−( )
η
⌠⌡
d−:= g x( )
0
b x( )
ηsin η x η
3⋅+( )
η
⌠⌡
d:=
wni f β i( )( ):= wpi g β i( )( ):=
wp
0
0
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
1.571
1.571
1.571
1.571
1.571
1.571
1.571
1.571
1.571
1.571
1.571
1.571
1.571
1.571
1.571
1.571
=wn
0
0
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
-1.550
-1.591
-1.616
-1.570
-1.544
-1.525
-1.537
-1.547
-1.556
-1.564
-1.565
-1.559
-1.547
-1.534
-1.526
-1.536
=
ββi β i( ):= zni1−
3β i( )
:= zpi1
3β i( )
:=
ββ
0
0
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
110 -7
1.02310 -7
1.04710 -7
1.07210 -7
1.09610 -7
1.12210 -7
1.14810 -7
1.17510 -7
1.20210 -7
1.2310 -7
1.25910 -7
1.28810 -7
1.31810 -7
1.34910 -7
1.3810 -7
1.41310 -7
= zn
0
0
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
-215.443
-213.796
-212.162
-210.539
-208.93
-207.332
-205.747
-204.174
-202.613
-201.064
-199.526
-198.001
-196.487
-194.984
-193.494
-192.014
= zp
0
0
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
215.443
213.796
212.162
210.539
208.93
207.332
205.747
204.174
202.613
201.064
199.526
198.001
196.487
194.984
193.494
192.014
=
Mathcad 11 CASA_06-27.mcd Appendix 12
Skalak, R. (1956), Ref. 4"An extension of the theory of water hammer."Transactions of the ASME 78, 105-116.