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 Ͱඇཧతͳද໘ுྗͷ ޮՌ͕ͳ͘ͳ͍ͬͯΔ͜ͱ ͕Θ͔Δ.