Reformulation of Density-Independent Smoothed Particle Hydrodynamics
with Riemann Solver: Godunov DISPH
౬ઙ1*ɾࢦಋڭ: ਖ਼1
1 ஜେֶ ཧֶ܈ ཧֶྨ, *[email protected]
1. SPH๏ͱ
• ཻࢠΛ༻͍ͯѹॖੑඇ೪ੑྲྀମ(ཧྲྀମ)ΛγϛϡϨʔγϣϯ͢
ΔͨΊͷܭࢉख๏.֨ࢠ๏ͱൺͨϝϦοτ͕͍͔ͭ͋͘Δ.
• ݩʑӉൃ.ۙۀCGͰ༻͍ΒΕ͍ͯΔ.
dv(r)
dt
= −
1
ρ(r)
∇P(r) + FAV
(α) (1)
du(r)
dt
= −
P(r)
ρ(r)
∇ · v(r) + GAV
(α) (2)
F(ri
) ≈
∫
F(r′
)W(|ri
− r′
|, h)d3r′
≈
∑
j
mj
ρj
Fj
W(|ri
− rj
|, h)
(3)
• ӡಈํఔࣜ(1),෦ΤωϧΪʔํఔࣜ(2)Λ,ࣜ(3)Λ༻͍ͯࢄԽ
͢Δ.W ΨγΞϯͷΑ͏ͳܗΛ͓ͯ͠Γ,पΓͷཻࢠͷӨڹΛ
ॏͶ߹ΘͤΔΠϝʔδ(i,j ཻࢠͷϥϕϦϯά). ཧؾମͷঢ়ଶ
ํఔࣜΛ༻͍Δͱ͕ࣜด͡Δ.
2.SPH๏ͷσϝϦοτ(৮ෆ࿈ଓ໘ͱিܸ)
৮ෆ࿈ଓ໘
ྲྀମͷڥք໘Ͱ͋Δ৮ෆ࿈ଓ໘, खͰนΛԡ͢ͱ͖ʹͰ͖Δෆ࿈
ଓ໘(खͱนͷ,ѹྗҰఆͰີ,ΤωϧΪʔෆ࿈ଓ)ͱࣅ͍ͯΔ.
SPHͰີࣜ(1)ͰٻΊΔ.
ρ =
∑
j
mj
Wij
(h) (4)
ͭ·Γ, ີͷۭؒ࿈ଓੑ͕Ծఆ͞Εͯ͠·͍ͬͯΔ. ·ͨ, ಋग़தͰ
ີͷۭؒඍΛߦ͍ͬͯΔ. ͜ΕΒͷཧ༝͔Β,SPH Ͱཧత
ʹີෆ࿈ଓͱͳΔ৮ෆ࿈ଓ໘͕͏·͘ѻ͑ͣʹඇཧతͳද໘ு
ྗΛൃੜͤͯ͞͠·͏͜ͱ͕ΒΕ͍ͯΔ.
ॲํᝦͱͯ͠,ີͷΘΓʹѹྗͷۭؒ࿈ଓੑ,ۭؒඍՄೳੑΛԾ
ఆͨ͠DISPH๏͕ଘࡏ. ඇཧతͳද໘ுྗΛ͑Δ͜ͱ͕Ͱ͖ΔҰ
ํ,িܸͰͷσϝϦοτSPHͱಉ༷ʹଘࡏ͢Δ.
িܸ
িܸԻͱಉ༷, ใͷൖʹΑΓى͜Δ. ཧతͳڥքͰͳ
͍. ѹྗ,ෆ࿈ଓ. ܭࢉͰিܸΛଊ͑ΔͨΊʹ,ਓ೪ੑ
߲((1),(2)ͷࣈͷ߲)ΛͬͯিܸΛ͋͑ͯΒ͔ʹ͢Δඞཁੑ͕
͋Δɻਓ೪ੑ߲ͷڧ͞Λௐઅ͢Δύϥϝʔλαਓ͕ؒखͰ,ܦݧ
ΛͱʹܾΊ͍ͯΔ.
ॲํᝦͱͯ͠,ཻࢠͷ૬ޓ࡞༻ܭࢉͷࡍʹRiemann solverͷղΛ༻͍Δ
͜ͱͰ,Ұͷύϥϝʔλͳ͠Ͱదͳ೪ੑ͕িܸྖҬʹೖΔ
Godunov SPH๏ʢGSPH๏ʣ͕ଘࡏ͢Δ. ͔͠͠,৮ෆ࿈ଓ໘ͷѻ͍
SPHͱಉ༷ʹෆશ.
3.Godunov DISPH๏ʹ͍ͭͯ
DISPH ๏ͱ GSPH ๏ΛΈ߹Θͤͨ Godunov DISPH๏Λ։ൃ. ৮ෆ
࿈ଓ໘Λਖ਼͘͠ѻ͑,దͳ೪ੑ͕ύϥϝʔλͳ͠ͰೖΔ.ҎԼͰͦ
ͷಋग़ʹ͍ͭͯܰ͘ड़Δ.
ཻࢠ i ͷඍখ෦ΤωϧΪʔมԽ, ཻࢠ͕֎෦͔Βड͚ΔମੵมԽ
Λ௨ͨ͡ࣄʹΑΔͷ(ॏ৺ͷҐஔมԽΛ௨ͨ͡ࣄӡಈΤωϧ
Ϊʔʹ)
dUi
= WV olume
i (5)
DISPHͰ(SPHͰ)
WV olume
i
= −Pi
dVi (6)
͜Εཻࢠi͕dtͷؒʹ֎෦͔Βड͚ΔѹྗɺPi
+ ϵͰ͋ΔͱԾఆ
͠ɺೋ࣍ͷඍখྔΛແࢹ͍ͯ͠Δ. (ྲྀମํఔࣜͷΤωϧΪʔͷࣜͰ
ಉ༷) ͜ͷԾఆͷԼͰ,ແݶݸͷཻࢠ,ແݶͷղ૾,ແݶখͷλΠ
ϜεςοϓΛ͑ਖ਼͍͠ղʹऩଋ͢Δͣ.
৽͍͠ߟ͑ํ
Pix
ΛɺཻࢠiΛத৺ͱͨ͋͠ΒΏΔํ,͔Βཻࢠi͕ड͚Δѹྗ
ͷ͋Δछͷۭ࣌ؒؒฏۉྔͱͯ͠
WV olume
i
= −Pix
dVi (7)
͕Γཱͭͱ͢Δɻͪ͜ΒͷԾఆͷԼͰ,ແݶݸͷཻࢠ,ແݶͷղ
૾,ແݶখͷλΠϜεςοϓΛ͑ਖ਼͍͠ղʹऩଋ͢Δͣ.
ۭ࣌ؒؒฏۉͷऔΓํʹҙੑ͕͋Δ. ࠓճ, ཻࢠ i ཻ͕ࢠ j ͔
Βड͚Δѹྗͷ࣌ؒฏۉΛ,ͦΕͧΕͷཻࢠͷཧྔΛೖྗͱͨ͠
Rieamnn solverͷղʹ͓͚Δ৮ෆ࿈ଓ໘্Ͱͷѹྗͱ͠,ۭؒฏۉ
దͳܗͰͱΔ. ͋ͱࣜ(7)͔ΒDISPH๏ͷಋग़ͱ΄΅ಉ͡Γํ
ͰGDISPH๏ͷӡಈํఔࣜ,ΤωϧΪʔํఔࣜΛಋग़Մೳ.
4. ςετܭࢉ(GDISPHͷ৮ෆ࿈ଓ໘ͱিܸʹର͢Δੑೳ)
ॳظ͕݅ѹྗฏߧͳͷͰ,
ॳظ݅ͷ··มಈ͠ͳ͍
͜ͱ͕ཧతͳղ.
SPH, GSPH৮ෆ࿈ଓ໘
͔Βͷද໘ுྗͷӨڹͰܗ
ؙ͕͘ͳͬͯ͠·͏.Ӎཻ͕
ද໘ுྗͰؙ͘ͳΔݮগͱ
ࣅ͍ͯΔ. DISPH, GDISPH
Ͱඇཧతͳද໘ுྗͷ
ޮՌ͕ͳ͘ͳ͍ͬͯΔ͜ͱ
͕Θ͔Δ.