「研究結果（１２月２日）ver.2」の編集履歴（バックアップ）一覧に戻る

研究結果（１２月２日）ver.2 - (2008/12/02 (火) 20:00:42) のソース

*シートのない場合
渦点の数（VN）＝４ &br()
時間刻み幅（h）＝０．０１ &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

&br()
*δ＝０．０で回転する球の場合　
T=12.0（秒）