**[９月２週目] テスト **[９月３週] シートの分割数を増やし、時間発展させると解が爆発してしまう。 [[研究結果（９月１６日）]] **[９月４週] 月曜日にM２の皆で、研究内容の発表会を行いました。 それに伴い、論文をもう一度読み返し、Γの値も論文通りに定め、 計算しなおしています。 しかし、まだ、正しい結果が出ていないのが現状です。 **[１０月１週] まだ、論文通りの絵は描けていないが、 球面からはずれることはなくなった。 **[１０月２週] 研究結果を添付します。 [[研究結果（１０月７日）]] 前回のUPから、先輩方とかなり詰めた確認をしているのですが、 論文通りには描写できていません。 ちょっとお手上げ状態になりつつあります。 **[１０月３週] 研究結果を添付します。 [[研究結果（１０月１４日）]] メインプログラムを渦点と渦シートの位置情報を出力するプログラムと、 datファイルから読み込んだデータをＧＬＳＣで描画するプログラムに、 改良した。 **[１０月４週] 1.サブルーチンuv_sを改正した。 2.0.1ずつ時間発展させたものを出力するようにした。 3.GLSCで描画する時に行っていた軸のずらしを、 メインの計算プログラムで行うようにした。 研究結果の絵は今週の中頃には載せる予定。 **[１０月５週] 実行速度を上げるためサブルーチンus_pとus_sを一つにまとめ us_psとした。 現在各速度場の見直し中。 現時点でのプログラムアップ（vor-data5.c）。 **[１０月３０日] 完成した絵の添付。 [[研究結果（１０月３０日）]] ***[１１月３日] 球面衝突について（近藤君のを参照） Γ：3.2.2.-1 で計算したタイプ。 [[研究結果（１１月３日）]] ***[１１月５日] １１月３日の条件下で、 Γの値だけ1/10に変更したものです。 図は、効率よく表と裏に分けるプログラムが まだ完成していないため、見にくい図のままです。 今後、改善でき次第、表と裏に分けてUPしていきます。 [[研究結果（１１月５日）]] **[１１月１４日] 渦点を北極近くに配置しました。 図を改良しましたので、添付します。 gammaの値（3,2,2,-1）を1/2にして計算しています。 [[研究結果（１１月１４日）]] **[１１月１８日] １１月１４日のデータにおいて、円の三等分を改正しました。 また、渦点の位置をより北極に近づけました。 点渦の距離は、数値だけだと見にくいため、グラフにして アップデートしました。 [[研究結果（１１月１８日）]] **[１１月２５日] １．非回転シートなし球面 ２．回転シートなし球面 ３．回転シートあり球面 [[研究結果（１１月２２日）]] 上記の３パターンの渦点の位置 を定量的に比較しグラフで表した。 [結果] 条件１では、渦点は衝突したが、 条件２、３では、衝突せず。 条件２の場合にのみ角速度を小さくすれば 衝突する位置がみえてくる。 **[１２月３日] シートがない条件のもと、 δ＝０．１とδ＝０．０での 非回転球で、渦点と渦点の距離をはかり、 比較を行った。 [[研究結果（１２月２日）ver.1]] 上の条件に回転する球という部分だけを 置き換えて比較を行った。 [[研究結果（１２月２日）ver.2]] **[１２月７日] １２月２日の研究結果の欠陥を改正。 [[研究結果（１２月７日）]] **[１２月１０日] ・ハミルトニアンの信頼性について ・渦点が衝突するといってよいライン引き [[研究結果（１２月１０日）]] **[１２月２２日] 回転系においてbase=pi/15とした 点渦の振る舞いと衝突について [[研究結果（１２月２２日）ver.1]] 初期値をz軸に対してpi/2傾けた。 [[研究結果（１２月２２日）ver.2]] **[１２月２８日] 回転系においてbase=pi/20として、 点渦の挙動だけを図示しました。 シートの強さを極端に小さくすると、 衝突したようにみえる挙動を示したが、 数値的には衝突したといえる数値には至らなかった。 [[研究結果（１２月２８日）]] **[一言] - テストコメント -- 秋山 (2008-09-10 10:56:26) - テストです。 -- 秋山ログインなし (2008-09-10 10:57:00) - 足跡 -- 近藤 (2008-09-10 16:36:21) - メールみました。内容確認しました。計算頑張ってください。 -- 坂上です (2008-09-10 21:29:32) - 結果見ました。なんか変ですね。考えられるのは点渦とシートの離散点の距離がほとんど０になっているかもしれないということですが、距離を毎ステップ確認しながら計算してみてはどうでしょうか？ それからとりあえず現在のプログラムを私のところに送ってください。こちらでも確認してみますから。 -- 坂上です (2008-10-08 16:06:35) - 結果みました。数値計算そのものはうまくいってるみたいですね。ただ点渦の強さが強すぎるみたいですね。大事なのはこの点渦の比なので、全体を１／１０ぐらいにして計算してみたらどうなりますか？シートの数は二枚でもいいかもしれません。もちろん渦の強さの比は3:2:2:-1ですよ。それから図ですが、透けて見える方法だととても見にくいので、表と裏を別々に表示するようにしてください。修士論文にのせるにしてもそのほうがいいと思います。 -- 坂上です (2008-11-04 08:10:47) - 次は、Γの値すべてを１/１０にして計算してみます。そして、図を表と裏で表示するようにします。 -- 秋山 (2008-11-04 19:51:05) - 点渦の初期配置ですが，あまり緯度を低くせずもっと北極近くにとって計算して見てください．11月３日の結果そのものはとてもおもしろいですが，渦の影響でシートがかなり引き延ばされていて近似の点数が足りなくなっているようです．これは点渦とシートの距離が近いことが原因ですので，点渦をもっと北極の近くにとるといいのではないでしょうか？ 近藤君がやっている問題ではそのパラメータはΘoで与えられています． -- 坂上です (2008-11-05 09:09:36) - 現在、点渦の位置を北極に近づけたものを計算中です。 -- 秋山 (2008-11-05 16:42:16) - 結果みました。点渦の初期配置が近藤君のものと同じになっていないような...。 円を三等分したら角度は0, 2/3*pi, 4/3*piですよ。やりなおしてください。また図はぐっと見やすくなりましたね。点渦が背景シートの影響で衝突するかしないかが問題になってくるので、この数値計算結果を利用して、各点渦の距離を測って、その様子をアップデートしてください。 -- 坂上です (2008-11-16 15:38:24) - 距離の関数，おかしくないですか？ 私の論文で出てきた曲線とは全く違うように思います．定義を確認してください．衝突するかどうかの判断は慎重にしましょう．あと不思議なのは衝突後に数値計算ができていることです．衝突したら速度が無限大になるので，それ以上の計算はできないと思いますがどうなのでしょうか？それから，回転のない時にハミルトニアンの値が変化している結果があるようですが，これは既に数値計算の精度が悪くなっていると思います．もう少し詳しくは後でメールします． -- 坂上です (2008-12-03 17:13:38) - 確認しました．精度の問題だったのですね．距離が限りなく０に近づくという時間ですが，どれくらい近くなっていますか？もちろん数値的には距離が０になるとは主張できませんから，ある程度「小さい」ことをもって衝突することを示すしかなさそうです．δが０で非回転の場合はかなり信頼度が高く衝突するので，あなたの計算でどれくらいまで距離が近くなっているか確認してください．また非回転の場合はハミルトニアンが保存するので，ハミルトニアンの値が大きく変化（たとえば何桁目の値がずれるなどの条件をつけてみて）してしまったら，それは「精度」が落ちていると見なして，以後の計算はしないようにしましょう．もはやそれは信頼性がありません．そうすれば，おそらくδが０の場合は最初に点が近づいた時点で計算はストップしなければならないでしょう．結果はとても興味深いので，回転していていも衝突する可能性はありますから，慎重にハミルトニアンの保存などを考慮しながら計算してください．回転していたらハミルトニアンは保存しませんが，その変化は「ゆるやか」で急激にあたいを変えたりすることはないはずですから，そのあたりの値の変化も見ながら「精度よく」「慎重に」計算を薦めてください． -- 坂上です (2008-12-08 07:17:50) - 結果みました．回転（２シート）の場合の距離のグラフを見せてください．それから，回転系の場合はハミルトニアンは変化するので，この結果で述べているようなハミルトニアンの値の変化による議論は正しくありません．もともとシートとの相互作用でハミルトニアンが変化するわけですから．それよりはたとえば点渦の距離が１０＾-5wo -- 坂上です (2008-12-17 18:59:59) - 途中でボタンを押してしまいました．続きです．たとえば衝突しそうな点渦の距離が10^-5を下回るかどうかを調べてください．さらにそれを5x10^-6とか, -- 坂上です (2008-12-17 19:01:48) - このWikiページは特定のキー操作に反応してしまうようです．続きです．さらにそれを5x10^-6とか1x10^-6とかにしたらどうなるかを調べてみてはどうでしょうか？数値計算なので本当に０になるかどうかは判定できませんが，どれくらい「ちいさい」値になるかは判定できます． -- 坂上です (2008-12-17 19:03:07) - わかりました。計算に約３日ほどかかってしまうので、処理できしだいUPします。 -- 秋山 (2008-12-18 00:20:21) - 計算結果みました。点渦の初期配置ですが、（例２のように）傾きを変える必要はありません。傾きは０のままで、baseの値を少しずつかえてみてください。また、点渦の距離のグラフを見ていてきになったのですが、距離が10^5とかそういうオーダーになっているのはもはや数値計算が破綻している証拠です。そもそも点渦の距離は０から４の値しかとりません。（気がついていましたか？）ですから４を越えるような値が出ている事自体計算がダメであることを物語っています。数値計算はt=60までもする必要はないでしょう。だいたい衝突が予測される時刻（１０程度から１５程度まで）で計算してください。そうすれば計算時間も１／４で済むはずです。 -- 坂上です (2008-12-23 17:48:34) - 結果みました。ずいぶんよくなりましたね。結論としては回転があると衝突はしないということですね。baseの値を変えても -- 坂上です (2008-12-29 02:02:13) #comment(title_name=名前,title_msg=コメント,nsize=15,size=100,vsize=1)
