研究結果（１２月２日）ver.1 - 数値計算 ver.秋山

シートのない場合

渦点の数（VN）＝４
時間刻み幅（h）＝０．０１

渦の強さ（Γ）

g[0]=3.0
g[1]=-1.0
g[2]=2.0
g[3]=2.0

渦点の位置

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

δ＝０．０で非回転球の場合　

T=12.0（秒）

点と点の距離

vx[0]:vx[1]

vx[0]:vx[2]

vx[0]:vx[3]

vx[1]:vx[2]

vx[1]:vx[3]

vx[2]:vx[3]

vy[0]:vy[1]

vy[0]:vy[2]

vy[0]:vy[3]

vy[1]:vy[2]

vy[1]:vy[3]

vy[2]:vy[3]

vz[0]:vz[1]

vz[0]:vz[2]

vz[0]:vz[3]

vz[1]:vz[2]

vz[1]:vz[3]

vz[2]:vz[3]

Jベクトル

横軸（時間（秒）） 縦軸(値）
※Jx,Jy,Jzが全て一定値(0.0)で保存。

ハミルトニアン

結果

点と点の比較からわかるように時間が６０秒過ぎると(vx[0],vy[0],vz[0])以外の渦点は衝突する。

δ＝0.1で非回転球の場合

T=12.0（秒）

点と点の距離

vx[0]:vx[1]

vx[0]:vx[2]

vx[0]:vx[3]

vx[1]:vx[2]

vx[1]:vx[3]

vx[2]:vx[3]

vy[0]:vy[1]

vy[0]:vy[2]

vy[0]:vy[3]

vy[1]:vy[2]

vy[1]:vy[3]

vy[2]:vy[3]

vz[0]:vz[1]

vz[0]:vz[2]

vz[0]:vz[3]

vz[1]:vz[2]

vz[1]:vz[3]

vz[2]:vz[3]

Jベクトル

横軸（時間（秒）） 縦軸(値）
※Jx,Jy,Jzが全て一定値(0)で保存された。

ハミルトニアン

結果

周期的に各渦点は接近するが、δ＝0.1の効果により衝突はせず、周期的な接近を繰り返す。