「研究結果(12月10日)」の編集履歴(バックアップ)一覧はこちら
「研究結果(12月10日)」(2008/12/10 (水) 17:48:24) の最新版変更点
追加された行は緑色になります。
削除された行は赤色になります。
*ハミルトニアンの信頼性について
初期位置や渦の強さは、前回の結果を引き継ぐ。
**δ=0.0で非回転の場合
***|h=0.001 && t=12000|
各時間におけるハミルトニアンの値。&br()
&br()
・0<=t<=10948 ⇒ -0.497523 &br()
・10948<=t<=10998 ⇒ -0.497522~-0.497520 &br()
・10999<=t<=11011 ⇒ -0.497520~-0.497510 &br()
・11012<=t<=11019 ⇒ -0.497510~-0.497500 &br()
・t=11020 ⇒ -0.497523 &br()
・t=11021 ⇒ -0.497563 &br()
・t=11022 ⇒ -0.497628 &br()
・t=11023 ⇒ -0.497751 &br()
&br()
という結果になったため、少数第4桁目が5から6に変わる時刻(t=11022)以降は
信頼性の低いものとして切り捨てて考えることとする。
この時、δ=0.0で非回転という条件下での、t=11021での渦点と渦点の距離は、&br()
|渦1と渦2|
0.000950
|渦1と渦3|
0.000759
|渦2と渦3|
0.000858
&br()
この距離をもって数値計算では渦点は衝突すると定義する。
**δ=0.0で回転の場合
***|h=0.001 && t=12000 && Ω=2.0π/12.0|
各時間におけるハミルトニアンの値。&br()
&br()
・0<=t<=10943 ⇒ -0.497523 &br()
・10943<=t<=10997 ⇒ -0.497522~-0.497520 &br()
・10998<=t<=11011 ⇒ -0.497520~-0.497510 &br()
・11012<=t<=11019 ⇒ -0.497510~-0.497500 &br()
・t=11020 ⇒ -0.497523 &br()
・t=11021 ⇒ -0.497555 &br()
・t=11022 ⇒ -0.497619 &br()
・t=11023 ⇒ -0.497739 &br()
&br()
という結果になったため、こちらも少数第4桁目が5から6に変わる時刻(t=11022)以降は
信頼性の低いものとして切り捨てて考えることとする。
この時、δ=0.0で回転する球という条件下での、t=11021での渦点と渦点の距離は、&br()
|渦1と渦2|
0.000949
|渦1と渦3|
0.000760
|渦2と渦3|
0.000859
&br()
***結果
Ω=2.0π/12.0の場合は、非回転の時と同様で、
同時刻において衝突すると言っても良いであろう結果が得られたことになる。
&br()
&br()
***|h=0.001 && t=12000 && Ω=2.0π/6.0|
各時間におけるハミルトニアンの値。&br()
&br()
・0<=t<=10937 ⇒ -0.497523 &br()
・10943<=t<=10995 ⇒ -0.497522~-0.497520 &br()
・10995<=t<=11010 ⇒ -0.497520~-0.497510 &br()
・11010<=t<=11014 ⇒ -0.497510~-0.497500 &br()
・11014<=t<=11019 ⇒ -0.497500~-0.497490 &br()
・t=11020 ⇒ -0.497508 &br()
・t=11021 ⇒ -0.497537 &br()
・t=11022 ⇒ -0.497596 &br()
・t=11023 ⇒ -0.497711 &br()
&br()
という結果になったため、こちらは、特に大きく少数第4桁目が特に大きく変化した時刻(t=11022)以降を信頼性の低いものとして切り捨てて考えることとする。
この時(t=11021で)の渦点と渦点の距離は、&br()
|渦1と渦2|
0.000949
|渦1と渦3|
0.000761
|渦2と渦3|
0.000859
&br()
***結果
Ω=2.0π/6.0の場合も、同様な結果が得られ、上のような条件下では、衝突すると言っても良い結果が得られたことになる。
&br()
&br()
***|h=0.001 && t=12000 && Ω=2.0π/2.0|
各時間におけるハミルトニアンの値。&br()
&br()
・0<=t<=10861 ⇒ -0.497523 &br()
・10862<=t<=10984 ⇒ -0.497522~-0.497520 &br()
・10985<=t<=10999 ⇒ -0.497520~-0.497510 &br()
・11000<=t<=11007 ⇒ -0.497510~-0.497500 &br()
・11008<=t<=11011 ⇒ -0.497500~-0.497490 &br()
・t=11012 ⇒ -0.497485 &br()
・t=11013 ⇒ -0.497477 &br()
・t=11014 ⇒ -0.497466 &br()
・t=11015 ⇒ -0.497452 &br()
・t=11016 ⇒ -0.497435 &br()
・t=11017 ⇒ -0.497416 &br()
・t=11018 ⇒ -0.497399 &br()
・t=11019 ⇒ -0.497385 &br()
・t=11020 ⇒ -0.497374 &br()
・t=11021 ⇒ -0.497372 &br()
&br()
という結果になったため、これは、少数第4桁が桁落ちした時刻(t=11008)以降を信頼性の低いものとして切り捨てて考えることとする。
この時(t=11008で)の渦点と渦点の距離は、&br()
|渦1と渦2|
0.001888
|渦1と渦3|
0.001496
|渦2と渦3|
0.001692
&br()
***結果
Ω=2.0π/2.0の場合は、非回転で衝突すると定義した渦点の距離とは、単純に約2倍の差がある。
これは、少々衝突したとは言いがたい結果となった。&br()
しかし、ハミルトニアンによる信頼性は低下するが、仮に(t=11021)以下を切捨てるとすると、
|渦1と渦2|
0.000956
|渦1と渦3|
0.000774
|渦2と渦3|
0.000870
&br()
となり、衝突したといえる距離となる。