研究結果(11月14日)

球面衝突
渦点の数(VN)=4
シートの数(M)=2
シートの分割数(DN)=2048
時間刻み幅(h)=0.01
δ=0.1

渦の強さ(Γ)

g[0]=1.5
g[1]=-0.5
g[2]=1.0
g[3]=1.0

渦点の位置

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/8.0

T=0.0
T=1.0
T=2.0
T=3.0
T=4.0
T=5.0
T=6.0
T=7.0
T=8.0
T=9.0
T=10.0
T=11.0
T=12.0
最終更新:2008年11月15日 02:34