「研究結果(12月22日)ver.1」の編集履歴(バックアップ)一覧はこちら
「研究結果(12月22日)ver.1」(2008/12/23 (火) 09:13:02) の最新版変更点
追加された行は緑色になります。
削除された行は赤色になります。
回転系(δ=0.0、シート二枚)の振る舞い
渦点の数(VN)=4 &br()
時間刻み幅(h)=0.01 &br()
&br()
***渦の強さ(Γ)&br
g[0]=3.0&br
g[1]=-1.0&br
g[2]=2.0&br
g[3]=2.0&br
&br
***渦点の位置&br
vz[0]=-cos(base);
vx[0]=sqrt(1.0-vz[0]*vz[0]) * cos(2.0*pi+ shift);
vy[0]=sqrt(1.0-vz[0]*vz[0]) * sin(2.0*pi+ shift);
vz[1]=cos(base);
vx[1]=sqrt(1.0-vz[1]*vz[1]) * cos(2.0*pi+ shift);
vy[1]=sqrt(1.0-vz[1]*vz[1]) * sin(2.0*pi+ shift);
vz[2]=cos(base);
vx[2]=sqrt(1.0-vz[2]*vz[2]) * cos(pi*2.0/3.0+ shift );
vy[2]=sqrt(1.0-vz[2]*vz[2]) * sin(pi*2.0/3.0+ shift );
vz[3]=cos(base);
vx[3]=sqrt(1.0-vz[3]*vz[3]) * cos(pi*4.0/3.0+ shift );
vy[3]=sqrt(1.0-vz[3]*vz[3]) * sin(pi*4.0/3.0+ shift );
base=pi/15.0
shift=pi/4.0
|T=0.0|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=27&file=%E5%88%9D%E6%9C%9F%E5%80%A4%EF%BC%91.gif,width=250,height=250)
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=27&file=%E5%88%9D%E6%9C%9F%E5%80%A4%EF%BC%92.gif,width=250,height=250)
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=27&file=%E5%88%9D%E6%9C%9F%E5%80%A4.gif,width=250,height=250)
|T=1.0|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=27&file=100.gif,width=250,height=250)
|T=2.0|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=27&file=200.gif,width=250,height=250)
|T=3.0|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=27&file=300.gif,width=250,height=250)
|T=4.0|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=27&file=300.gif,width=250,height=250)
|T=5.0|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=27&file=500.gif,width=250,height=250)
|T=6.0|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=27&file=600.gif,width=250,height=250)
|T=10.0|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=27&file=1000.gif,width=250,height=250)
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=27&file=1000z.gif,width=250,height=250)
|T=12.0|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=27&file=1200.gif,width=250,height=250)
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=27&file=1200z.gif,width=250,height=250)
|T=14.0|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=27&file=1400.gif,width=250,height=250)
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=27&file=1400z.gif,width=250,height=250)
***ハミルトニアン
|t=0~60|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=27&file=h.gif,width=300,height=300)
|t=0~10|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=27&file=h0-1.gif,width=300,height=300)
***点と点の距離
|渦1と渦2|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=27&file=kyori0-6.gif,width=300,height=300)
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=27&file=kyori0-1.gif,width=300,height=300)
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=27&file=kyori1-1.8.gif,width=300,height=300)
***結果
点渦は衝突しない。
なぜなら、非回転系と比較した際、このコンピュータでは10^{-5}程度接近しないと衝突するとは言えなかったが、
回転系で上記の初期値では、点渦の距離が10^{-2}程度にまでしか接近しないため。