îIZ_dgI`e'%úŶŒŀĒ5úŶÌĒ ō)ű&Ñġ:%9Ɖ:7+âĒ+ ŦĘ,ƈŕĀƁ±*ą½=ŗÀƈŖĸ=ŏ'&amp

(1)

3/&3 *794D!6)3B 8<

à ÈŎ

¹^Ƈ

fÇĢ Č

²

1m¼»¾ĞÊ¾ŭĭĵİŮÍğ¶¾Ãìƅ112-8551ÿsŮïsõò1-13-27Ɔ

2m¼»¾ĞÊ¾ŭŮÍğ¶¾İƅ112-8551ÿsŮïsõò1-13-27Ɔ Ƈ E-mail: [email protected]

\PI_VbgâĒ,ą½Ē'Ď/%ƈą½=ŗÀ9ÚŒ) 3ƈ»°#őżŕĀ Ɓ±ƈ9,Ĕīē5¶ģ+»¸Õ=|ēŪ*6Ûġ:%9Ɖ\PI_Vbg âĒ*9ĶŴħ)ŽíâĒ,ĴÕÔ'æÕÔq#ĳƃ9Ɖq#+ĳƃ+Ží

]Rc=ĎŢ9'ƈæÕÔ,ŖĸGKSĴÕÔ68ŰĿÒæÕÔ68~

3ƈØ7: Ńā,kŌĚ)îåœ7:9Ɖhñ&æÕÔ,aFaeJ_ħ)ëůñ ĲÔ*Ŭġ9+,žƉ

wi=Š0%ƈüĭĵ,æÕÔ+ŴŽíâĒ&aFaeJ_ħ)ëůñĲÔ*Ŭġ ƈ:*69ıēƀ+ĞâĒ=éĆ Ɖ

Key Words :Meshfree Method, Smoothed Particle Hydrodynamics, Moving Particle Semi-implicit

1.

îIZ_dgI`e'%úŶŒŀĒ5úŶÌĒ ō)ű&Ñġ:%9Ɖ:7+âĒ+

ŦĘ,ƈŕĀƁ±*ą½=ŗÀƈŖĸ=ŏ'&

8ƈą½Ē'¦-:%9Ɖą½Ē=ġ Ŗĸ+¿

ÀÝ5ĿÒ,ą½+ŗÀñĒ*6"%ƈ»ĥ)9 3ƈą½+ŗÀâĒ*ŵ9ĭĵº:%9Ɖ ƈ»°#őżŕĀƁ±ƈ9,Ĕīē5¶ģ+

»¸Õ=|ēŪ)(+ŕĀÂŝ*ą½Ē=Ŭġ 9+,ªž&9Ɖ

+6)ŋ÷+)ƈą½=ÚŒ')\PI_

VbgMeshfreeŖĸâĒ+ĭĵĨ>*ŏ<:%%

9Ɖ\PI_VbgĒ,ą½=ŗÀ) 3ƈŨń ŵî=Ő*,ƈ9Ɓ±+őî+ŖĸĘ+Þ³=

ġ%ƈ+Ɓ±+Ũńŵî=ĉĻ9ÚŒ9Ɖ ĉĻâĒ,ƈ»:9'ĴÕÔ'æÕÔqĳ ƃ9¹⁾Ɖ

q#+ĳƃ+Ží]Rc=ĎŢ9'ƈæÕÔ+

]Rc,æ=9 3ƈņÕñĲÔ=ŕÚŒ8ƈ ŖĸGKSĴÕÔ+]Rc68ŰƉ:*Âƈ ĴÕÔ+]Rc,+ñĲÔ* ý9ƅŬ£ÝƆ*

,Ŷþy=ĕ ):-)7)+&ƈĿÒ

æÕÔ+]Rc68~Ɖ

hñƈ\PI_VbgâĒ,Ŗĸm&Ũńŵî=ĉĻ 9 3ƈŖĸĘ+ı*Â%ƈą½Ē68¿Àħ*

Ŗĸ&9Ɖ

ē+ëůñĲÔ&9N-SñĲÔ,BAag'aF aeJ_q#+ÕÔ8ƈÖŊ,Ŋ'Ď/9'ƈı ēƀ=Ŗĸ9ÚŒ)Ę9Ɖt*ƈaFae J_ħ)ŖĸâĒ,ŖĸĘē+Ū''4*ı

%9 3ƈ\PI_VbgâĒ*Ŭ£9Ɖ a F a e J _ + \ P I _ V b g â Ē ' % ƈ SPH(Smoothed Particle Hydrodynamics)ƈMPS(Moving Particle Semi-implicit),ÏÑÛġ:%9ƉSPHƈMPS,Ĵ

ÕÔ+Ží]Rc=ġ9 3ƈØ7: ®,ƈ kŌĚ)å=9'Ĭ7:%9²⁾ƉhñƈĪç ħ*æÕÔ+]Rc=ġ%aFaeJ_ÕÔ+ëů ñĲÔ=Ží:-ƈŖĸĘ+ė´£&ƈ ŖĸĦí9'Ĭ7:%9Ɖ

wi=Š0%ƈüĭĵ&,ƈŖĸĘ=Àƈæ ÕÔ+]Rc=ġ ıēƀ=Ğ9âĒ+ųĦ=ĩ Ċ'9Ɖ

第23回地球環境シンポジウム講演集　2015年9月

42.

- 173 -

(2)

2.;"(' *('C,#

ĴÕÔ+ĶŴŽí]Rc,ƈ9´á+Ũńŵî

=Őö*ƈ+´á+¥ţ++Ű1uÐ¯=

Ɖhñƈ\PI_VbgâĒ*,ĴÕÔ*Â%ƈ

æÕÔ+ĶŴŽí]Rc9Ɖ:,ƈ9´á +ěĞŲ=Őö*ƈ9æŵî=ŗÀƈ+´

á+¥ţ+*ŵ9Þ³=ġ%æ9Ɖ

(1) ;"('C,#

ĴÕÔ+]Rc+²üĈÜ,Ô(2.1)*²$%9Ɖ Ô(2.1)=Ží9'ƈÔ(2.2)*)9Ɖ

φ(r )= φ(r'

)δ(r'

−r )dr'

∫

Ω ^(2.1)

φ(ri

)= 1 n0

φ(rj

)

∑

j≠i ^w(^r^j⁻^rⁱ⁾ ^(2.2)

,KCagŵîƈ,ÓĊXEScƈjuj,ÓĊ +Ť*9ļ½+Ĥ¢ƈ,Ű1ŵî&9ƉŤ

ļ½=ćŁ9ĺ«'Ð¯ļ½ŴşŽęŶ*Å )9'ƈÔ(2.2),Ô(2.1)* ý9Ɖ

Ô(2.2)*×"%ƈKCagŵî'ÙñĲÔ+ƀ=

:ø9'ƈÔ(2.3)ƈ(2.4)+ů]Rc.aWaI

?e]RcØ7:9Ɖ

∇φ _i= d n₀

φj−φi

r

j−r

i 2

∑

j≠i ^(r^j⁻^rⁱ⁾^w(^r^j⁻^rⁱ⁾ ^(2.3)

∇²φ

i= 2d

λn0

(φj−φi)

∑

j≠i ^w(^r^j⁻^rⁱ⁾ ^(2.4)

(2) *('C,#

üĭĵ,æÕÔ+Ží]Rcm+ıùÅqn Ē=ġ9ƉæÕÔ+Ží]Rc,æŵî+

î&ÙñĲÔ+ƀ=ŐƉ

φj−φi=a₁x_ij+a₂y_ij+a₃x_ij²+a₄x_ijy_ij+a₅y_ij² (2.5)

Ô(2.5)+6*ƈºƀÔ=ġ%ļ½ i+Ť&æ 9ƉŤļ½+=-ƈæŵî+î +ùÅqnèÀØ7:9Ɖ×"%ƈÔ(2.6)*į6

*QgagÆų=ġ:-ƈ7ÙñĲÔ=

Ží&9Ɖ

∇φ =(∂φ

∂x,∂φ

∂y)=(a1,a2)

∇²φ =∂²φ

∂x²+∂²φ

∂y² =2a3+2a5

(2.6)

(3) C,#2A

ĴÕÔ'æÕÔ+qĳƃ+Ží]Rc+Ží

ā=ĎŢ9 3ƈ9ĩĊŵî=ŗÀƈ qĳƃ+]Rc*69hŹÙ+ŽíŃā=ĎŢ 9Ɖ

%-1(a),ĴÕÔ+]Rc+ŖĸŃā=į%9Ɖ Ë¬,ļ½Ŕůň:9´£+ŖĸŃā=į%

8ƈÙ,Ğś'²üħ*¤'9Ɖ ƈ¡¬*įļ½kŔůň:9´£ƈÙ

,Ğś'»ĥ)9'Įř& Ɖ+

©,ƈļ½ą½Ĝ*ůň:%9+,ĴÕÔ +Ŭ£Ý=Ř9Ŷþy&9Ɖ

hñƈ%-1(b)*įæÕÔ+]Rc+Ńā,ƈļ½ Ŕůň.kŔůň+´£*%4ƈŖĸƁ

±+/%+ļ½&Ù,Ğś*ſÎ*Ť

=į%9Ɖ+Ńā,ƈļ½+ůň*Â9xÀ

=Ä)" '*9',ó7&8ƈ

æÕÔ+]Rczßļ½ůň+ůŖĸ=ŏ+*

úġ)âĒ&9'=į%9Ɖ

(a)ĴÕÔŽí]Rc*69ĶŴŽí

(b)æÕÔŽí]Rc*69ĶŴŽí

%-1ĴÕÔ'æÕÔŽí]Rc+ĎŢ

- 174 -

(3)

3..:'C,#

üĭĵ,ĶŴħ)Ží*ıùÅqnĒ=ġƈ öŴħ)Ží,ŷħ)ŕĒ=ġ9Ɖŷħ)ŕ Ē,ƈhNA[KQPW=ŸħčŹ'ŷħčŹq#*

%Ŗĸ=ŏƉŸħčŹ&,ľÝ'¹=Ŗĸƈ ļ½*ġƈļ½}ň=Ŗĸ9ƉŸħ)čŹ&Ø 7: ļ½+}ň,ŨńÔ=ĕ ) 3ƈŷħ)č Ź&®Y?MeñĲÔ=ŕƈØ7: ®=ġ%

ļ½+}ň=ŨńÔ=ĕ 6*Č9Ɖ

%-2MPSĒ+?cHbL[

%-2*į+,MPS+?cHbL[&8ƈĴÕÔ +Ží]Rc<:%9ƉÔ(3.1),MPS+Ží

®Y?MeñĲÔ=į%9Ɖ,Ÿħ)č ŹØ7: }ň7Ŗĸ ļ½îÁÒƅ9ļ½+

Ť+ļ½+łƆƈ,Ŗĸû+ļ½îÁÒ&

9Ɖ

∇²P=− ρ0

Δt²(n^*−n⁰

n⁰ ) ^(3.1)

æÕÔ+Ží]Rc,ƈŰ1ŵî=ġ%) 3ƈÔ(3.1)=ġ9'&)Ɖļ½Ť+

ěĞŲ+Ű1uÐ¯&ļ½+ěĞŲ=Ő%9+

&ƈļ½+}ňŨńÔ=ĕ )'ƈŰ1ŵî+

¸<9Ɖ+ 3ƈļ½îÁÒ,Ŗĸû'ĥ)9Ɖ ƈæÕÔ+Ží]Rc=ġ:-ƈŰ1ŵ î)Ɖ+ 3ƈļ½+}ň+¸,ē+ÁÒ+

¸=ô9'&)Ɖ+v<8*æÕÔ +Ží]Rc=ġ9´£&,ƈÔ*į6

*Ÿħ)čŹ&Ø7: ŧÒ+Ħí=®Y?Meñ ĲÔ+MgKƀ'%Ɖ

∇²P=− ρ

Δt²∇u^*_(3.2₎

4. *('C,#@=>0

(a)ĴÕÔŽí]Rc*69ďăÉ·Ŗĸ

(b)æÕÔŽí]Rc*69ďăÉ·Ŗĸ

%-3ĴÕÔ'æÕÔŽí]Rc*69ďăÉ·Ŗĸ

%-3,qċďăÉ·+ŖĸŃā=įƉŇš,®

=Ő%9ƉŖĸ+û,ļ½Ŕůň:%

9 3ƈĴÕÔ'æÕÔ,lñ'4¿Àħ*Ŗĸ

&9'9ƉöŴũ2'ƈļ½+ůňp :9 3ƈĴÕÔ+ŖĸŃā,»)®êpœ 7:9Ɖ:*Â%ƈæÕÔ,ŖĸûĴÕÔ '¤6*Ŗĸ=¿Àħ*ũŏ9Ɖ';ƈö Ŵ7*ũ2'ƈĴÕÔ+ŖĸŃā,9ĲÒ®

å%9ƈŖĸ,Ħí)Ɖhñƈæ ÕÔ+ŖĸŃā,ļ½+şŽŽ:ū&ƈŖĸĦ í%0'9Ɖ

- 175 -

(4)

5. /&3 *794D!6) 3

Ĺ*ť/ 6*ƈæÕÔ+Ží]Rc&Ŗĸ 9ö*ƈļ½îÁÒƅŖĸĘ+łƆhÀ*)9

Ř)Ɖ)<!ƈaFaeJ_ħ) ñĲÔ=

Ží9ź*ƈŖĸĘē''4*ı%9Ɖ +ŃāƈŖĸĘł)" 8Ľ)" 89Ɖ +¨ƂĘ+ŕđķ'%ƈh#ŉ7:9ñĒ,Ŗ ĸĘ=À9'&9Ɖ

MPSĒ+hNA[KQPW+ùÖ,®ů7ŧÒ

5}ň=ČƈØ7: Ńā,k+1NA[KQPW}

ň+ěĞŲƈ&9ƉŖĸ

Ę+}ň=À9 3*ƈ%-4+Ă+mį 6*

k+1NA[KQPW}ň+ěĞŲƈ

=Ŗĸ9ÚŒ9Ɖ

ħ)ñĒ,ƌƊƈ}ň+ěĞŲƈ

=ġ%ƈÔ(2.5)ƈ(2.6)7æŵî+î

=Đ39ƉƋƈŖĸû+ļ½}ň=æŵî*v

ƈ}ň+ěĞŲ=Đ39Ɖ

%-4üĭĵéĆ9ŖĸâĒ+?cHbL[

6.+1 )3@=>0

üĭĵéĆ âĒ&ƈdATcLî +D

^UQ@ē:+Ŗĸ=ŏ" Ɖ%-5,ŖĸŃā=į%

9Ɖ

(a)ļ½Ŕůň:9

(b)ļ½aeO[ůň:9

%-5éĆ9âĒ*69D^UQ@ē:+ŖĸŃā

7.

üĭĵ,\PI_VbgŖĸâĒ*9ĴÕÔ'

æÕÔ+ĶŴŽí]Rc=ĎŢ ƉĴÕÔ+Ķ ŴŽí]Rc+Ŭ£Ý*6"%kŌĚ)®êpĠ 9'=į Ɖ0 ƈæÕÔ+ĶŴŽí]Rc

=ö*ƈŖĸĘ+}ňė¸9´£ƈŖĸ Ħí%0'=Įř Ɖ

:*Â%ƈüĭĵ,ðâĒ=éĆ ƉMPS Ē*ùÅqnĒ=Äƈæ=ŏƉ:*6"%a FaeJ_ÕÔ+ëůñĲÔ=Ží%4ƈŖĸĘ=

À&9Ɖ

%-5+ŖĸŃā7œ9'ƈļ½+ůň*ŵ<7ƈ Ŗĸ=¿À*ŏ'&ƈüâĒ+úÝ=ćŘ Ɖ

$?-5

1)G.R.Liu Meshfree Methods, CRC Press, 2010, pp. 13-20.

2)ĝr, ŞµŚh, ĄĢ§o, MPS Ē*9ƄĿ

ųĦ, ŖĸÊ¾ŜĖ{śïŻ, 17, C-2-2, 2012

- 176 -

îIZ_dgI`e'%úŶŒŀĒ5úŶÌĒ ō)ű&amp;Ñġ:%9Ɖ:7+âĒ+ ŦĘ,ƈŕĀƁ±*ą½=ŗÀƈŖĸ=ŏ'&amp

3/&3 *794D!6)3B 8<

à ÈŎ

fÇĢ Č

42.

∫

∑

∑

∑

îIZ_dgI`e'%úŶŒŀĒ5úŶÌĒ ō)ű&Ñġ:%9Ɖ:7+âĒ+ ŦĘ,ƈŕĀƁ±*ą½=ŗÀƈŖĸ=ŏ'&amp