研究結果(12月22日)ver.2

回転系(シート二枚、δ=0.0)の振る舞い

渦点の数(VN)=4
時間刻み幅(h)=0.01

渦の強さ(Γ)

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/15.0
shift=pi/4.0
上記の条件をもとに、初期配置をz軸に対してpi/2傾けた.

T=0.0
T=0.5
T=1.0
T=1.5

ハミルトニアン

点と点の距離

t=0~60
t=0~10
t=0~1.5

結果

衝突しない。理由はver1と同様. 違うところは、渦点に対して、シートが明らかに渦点の周りを巻くような形を形成した点。

最終更新:2008年12月23日 09:17