*シートのない場合&br() 渦点の数(VN)=4 &br() 時間刻み幅(h)=0.01 &br() δ=0.1 &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 ***non-rotating sphere |T=12.0| #ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=22&file=322-1-non-8pi-12y.gif,width=250,height=250) #ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=22&file=322-1-non-8pi-12.gif,width=250,height=250) ***結果 渦点は衝突する。 ***rotating sphere 角速度(Ω)=2.0*pi/12.0 の場合$br() |T=12.0| #ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=22&file=322-1-ro-8pi-12.gif,width=250,height=250) #ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=22&file=322-1-ro-8pi-x.gif,width=250,height=250) #ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=22&file=322-1-ro-8pi-z.gif,width=250,height=250) 角速度(Ω)=2.0*pi/24.0 の場合$br() |T=12.0| #ref(http://www32.atwiki.jp/20073001?cmd=upload&act=open&pageid=22&file=322-1-ro-8pi-24.gif,width=250,height=250)