「研究結果(12月2日)ver.1」の編集履歴(バックアップ)一覧はこちら
「研究結果(12月2日)ver.1」(2008/12/02 (火) 20:01:28) の最新版変更点
追加された行は緑色になります。
削除された行は赤色になります。
*シートのない場合&br()
渦点の数(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/8.0
shift=pi/4.0
&br()
*δ=0.0で非回転球の場合
T=12.0(秒)
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=322-1-8pi-non-0.0.gif,width=250,height=250)
***点と点の距離
|vx[0]:vx[1]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-x-01-0.0.gif,width=350,height=350)
|vx[0]:vx[2]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-x-02-0.0.gif,width=350,height=350)
|vx[0]:vx[3]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-x-03-0.0.gif,width=350,height=350)
|vx[1]:vx[2]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-x-04-0.0.gif,width=350,height=350)
|vx[1]:vx[3]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-x-05-0.0.gif,width=350,height=350)
|vx[2]:vx[3]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-x-06-0.0.gif,width=350,height=350)
|vy[0]:vy[1]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-y-01-0.0.gif,width=350,height=350)
|vy[0]:vy[2]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-y-02-0.0.gif,width=350,height=350)
|vy[0]:vy[3]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-y-03-0.0.gif,width=350,height=350)
|vy[1]:vy[2]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-y-04-0.0.gif,width=350,height=350)
|vy[1]:vy[3]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-y-05-0.0.gif,width=350,height=350)
|vy[2]:vy[3]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-y-06-0.0.gif,width=350,height=350)
|vz[0]:vz[1]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-z-01-0.0.gif,width=350,height=350)
|vz[0]:vz[2]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-z-02-0.0.gif,width=350,height=350)
|vz[0]:vz[3]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-z-03-0.0.gif,width=350,height=350)
|vz[1]:vz[2]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-z-04-0.0.gif,width=350,height=350)
|vz[1]:vz[3]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-z-05-0.0.gif,width=350,height=350)
|vz[2]:vz[3]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-z-06-0.0.gif,width=350,height=350)
***Jベクトル
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=j-non-0.0.gif,width=350,height=350)
横軸(時間(秒))
縦軸(値)&br()
※Jx,Jy,Jzが全て一定値(0.0)で保存。
&br()
***ハミルトニアン
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=hami-non-0.0.gif,width=350,height=350)
***結果
点と点の比較からわかるように時間が60秒過ぎると(vx[0],vy[0],vz[0])以外の渦点は衝突する。
&br()
*δ=0.1で非回転球の場合
T=12.0(秒)
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=322-1-8pi-non-0.1.gif,width=250,height=250)
***点と点の距離
|vx[0]:vx[1]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-x-01-0.1.gif,width=350,height=350)
|vx[0]:vx[2]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-x-02-0.1.gif,width=350,height=350)
|vx[0]:vx[3]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-x-03-0.1.gif,width=350,height=350)
|vx[1]:vx[2]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-x-04-0.1.gif,width=350,height=350)
|vx[1]:vx[3]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-x-05-0.1.gif,width=350,height=350)
|vx[2]:vx[3]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-x-06-0.1.gif,width=350,height=350)
|vy[0]:vy[1]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-y-01-0.1.gif,width=350,height=350)
|vy[0]:vy[2]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-y-02-0.1.gif,width=350,height=350)
|vy[0]:vy[3]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-y-03-0.1.gif,width=350,height=350)
|vy[1]:vy[2]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-y-04-0.1.gif,width=350,height=350)
|vy[1]:vy[3]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-y-05-0.1.gif,width=350,height=350)
|vy[2]:vy[3]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-y-06-0.1.gif,width=350,height=350)
|vz[0]:vz[1]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-z-01-0.1.gif,width=350,height=350)
|vz[0]:vz[2]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-z-02-0.1.gif,width=350,height=350)
|vz[0]:vz[3]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-z-03-0.1.gif,width=350,height=350)
|vz[1]:vz[2]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-z-04-0.1.gif,width=350,height=350)
|vz[1]:vz[3]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-z-05-0.1.gif,width=350,height=350)
|vz[2]:vz[3]|
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=non-z-06-0.1.gif,width=350,height=350)
***Jベクトル
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=j-non-0.0.gif,width=300,height=300)
横軸(時間(秒))
縦軸(値)&br()
※Jx,Jy,Jzが全て一定値(0)で保存された。
&br()
***ハミルトニアン
#ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=23&file=hami-non-0.1.gif,width=350,height=350)
***結果
周期的に各渦点は接近するが、δ=0.1の効果により衝突はせず、周期的な接近を繰り返す。
&br()