Page 113 - 《水资源与水工程学报》2023年第2期
P. 113
0
! 2# òóô, =: SPHvõö÷ø¶Õùúûüýþ5./ 1 9
2
/
c ρ 0 γ8½。 Ó5; >VU÷f 5 k Wendland VUwax;£
$æ åV(|iæÔÕ8½,Colagrossi ()f; ê7Þß÷fñÏFÕê7ÞßÐ
[13]
, ®'E3Ù(3) iAghÔÕ χ ²*ÅÂ £x;·O, Ük|Íx, ¿ê7Þß'¿ß[
~Ù=, Ú±¯»@23Õ) [14] ²«aþ 0//F|i<mu; GÀEF#Û
V(¦7, ÐÔ?±oE3G: Éǰ, ÷fÀɬÕ/{ [17] $#ÛÉ
ρ
]
2
P =B [ ( ) γ -1 +χ -βρ (4) -{; @23¸]fÐ÷l$êUè,
ρ 0 '#ÛÉ/F@3|i, @23¸G<Ç
Ùi: U βÉɬÕ34。 |u [16] , Ú±È8½Ï, ÷f
2.2 N-S W'h SPH£ @23¸²Apvn£]fw/®。
Ük SPH)vRz, SPHÎÙ N-S E
3 5u©Úñ
3ÖE3ÓFÖE3Éï«\。
ÖE3: 3.1 ©°ªA÷
dρ i ®x¼zûi8W×x®Ý9 Chanson
= ∑ mv W (5)
dt j jij i ij , [18] 'ÈB5\&H
ûÙ¿ß[×x, ?@
FÖE3: |ÇÈ_®×xÇÈ$ú³¿²x¼U
dv P +P j è®u。¿ûÙ¿ß[®®G 2.70
i
i
=- ∑ m i W +g (6)
dt j j ρ i ρ j ij m, G 1.00m, û0G 21.8°, ?@.Â
Ùi:mGx;;v =v-v G i x;Ó j x; ® 0.60m`D0ûX, i û a=
ij
i
j
Kµ,m/s;WG SPH`a>VU; i W G i 0.25m, + b=0.10m, ûÕ0oÕp, ®
ij
x; j x;>VUS。
«¬ 1 。' \ z û ® × x i, 6 C À
GÀ"#|ifÙ, Ð*Å/ 3
0.046~0.182m /s /©ªËû/oÄn
ëiÜDÔÕEF,Antuono , [15] 'ÖE3
Ó0/¹¸, JÝ9ÀÃiÚ 2l Q222
i»pv), ,»À δ-SPH»Ù。 2
1, Ñê/G 0.058m /s, ëGf21\^_a
dρ i
0 ∑
- )·
= ∑ mv W +2 δ hc· [( ρ j ρ i þlÓa0/Á0/ÎÏ, г
dt j jij i ij j
x¼Uè®'ÙÚ/o\§f¸Ó|u。
r W · m j ] (7)
ij i ij
2
r
ij ρ j
Ùi: δGpv)U, 'J|ièG 0.1;
h G`a®,m;r =r-r G i x;_ j x;(
j
ij
i
),m。
'FÖE3i, |÷f Monaghan [16] ®
@23¸ÐApv, $ê|@3iÜD
8W, ] SPH|á=。 H 1 ÕÖ[MêA÷Z[H( ¯ê:m)
dv
i =- ∑ ( P +P j + ) i W +g (8) 3.2 S{÷
i
m
dt j j ∏ ij ij
ρ i ρ j
J1f DualSPHysics_v4.4 Uè|, U
-αμ ij 珋
c
{ ij (v·r <0) è®^4_×x® 1∶1 úût$, G¥^
ij
ij
∏ = ρ ij (9)
ij {ôûx;´É$|u34, JUè
0 (v·r >0) |÷fæ åV( SPH®, U讫¬ 2
ij
ij
Ùi: ij 。?@»/ê7Þßó5²wa0o=,
∏ G@23¸;αG@23¸U, #=À
2
2
2
2
=hv·r/(r+η), η =0.01h Ük|Íx, R>x;£(ó5G 0.005m±,
3¸pvh; μ ij ij ij ij
f£ÛVx;((ÑûnSj³z@9; ª¼§%8½ú®\«m+|u, f
珋 =0.5(c+c) G¬LK,m/s。
c i j ®s 453000 x;, |±£G 20s, 1f Gefore
ij
6CiÃõ4ÇUó5Ó(|E7 RTX2080GPU@ 1.8GHz àÊ GPU¾8
G: ±£s³÷f Verlet s³¶·á*x;K ½, 8½p±G 34~35h。