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

球面衝突

渦点の数（VN）＝４
シートの数（M）＝２
シートの分割数（DN）＝２０４８
時間刻み幅（h）＝０．０１
δ＝０.１

渦の強さ（Γ）

g[0]=0.3
g[1]=-0.1
g[2]=0.2
g[3]=0.2

渦点の位置

vz[0]=-cos(base); vx[0]=sqrt(1.0-vz[0]*vz[0]) * cos(2.0*pi); vy[0]=sqrt(1.0-vz[0]*vz[0]) * sin(2.0*pi);

vz[1]=cos(base); vx[1]=sqrt(1.0-vz[1]*vz[1]) * cos(2.0*pi); vy[1]=sqrt(1.0-vz[1]*vz[1]) * sin(2.0*pi);

vz[2]=cos(base); vx[2]=sqrt(1.0-vz[2]*vz[2]) * cos(2.0*pi*1.0/3.0+ shift ); vy[2]=sqrt(1.0-vz[2]*vz[2]) * sin(2.0*pi*1.0/3.0+ shift );

vz[3]=cos(base); vx[3]=sqrt(1.0-vz[3]*vz[3]) * cos(2.0*pi*2.0/3.0+ shift ); vy[3]=sqrt(1.0-vz[3]*vz[3]) * sin(2.0*pi*2.0/3.0+ shift );

base=pi/4.0

T=0.0

T=1.0

T=2.0

T=3.0

T=6.0

T=7.0

T=8.0

T=9.0

T=10.0

T=11.0

T=12.0