ExcelでFermi-Pasta-Ulam問題をシンプレクティック法を使って計算してみる。
この前、シンプレクティック法（8次）で3体問題を計算したけど、次にフェルミ・パスタ・ウラムの問題（いわゆるFPU問題）を解いてみた。
元論文がwebに転がっていた。
http://www.physics.utah.edu/~detar/phycs6720/handouts/fpu/FermiCollectedPapers1965.pdf
使った式はこれ。
α=1/4、N=32としてシンプレクティック8次でモード毎のエネルギーを計算した結果がこれ。
うん、元論文の結果を再現してるな。最初のモードに戻って行ってる。
しかし結構これ時間かかるぞ。上の計算で10分くらい。FPUの時代はMANIAC Iってコンピュータでleap frog法を使ったみたいなんで、一体どれくらいの時間かかったんだろう。
MANIACについては以下に。
http://electronic-journal.seesaa.net/article/6436566.html
FPUについての発展はこんなのが出てた。
« ランダウ・リフシッツ方程式（教科書じゃない） | トップページ | ガストフロント/ダウンバーストって？（福井の突風） »
「学問・資格」カテゴリの記事
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2020年9月17日 IEEE microwave magazine10月号はマイクロ波顕微鏡特集、 Microwave Journalはミリ波パッケージ。Vayyarの4Dイメージングレーダの仕組み、MIPIのA-PHY、など。(2020.09.17)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2020年8月17日 Microwave Journal 8月号はドローンで衛星通信をテスト、220GHzの半導体テストシステム。DesignConで導体の粗さモデルと因果律について解説。Qualcommとワシントンポストの5Gクラス、Samsung Galaxy Watch3分解(2020.08.17)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2020年8月10日 Microwave Magazine 9月号はシリコンフォトニクスをマイクロ波・ミリ波通信に使う、コロナと5G、IEEE RFICでimecが低消費電力60GHzレーダを発表、Qorvoの5G向けフィルタ解説、など。(2020.08.11)
- Python＋Scipyでルンゲクッタ8次のDOP853(Dormand＆Prince)を使う（その16) 時間依存のシュレーディンガー方程式のトンネル効果を計算してGIFアニメに。(2020.07.20)
- Python＋Scipyでルンゲクッタ8次のDOP853(Dormand＆Prince)を使う（その15) 蔵本モデルを計算してGIFアニメに。 (2020.07.18)
« ランダウ・リフシッツ方程式（教科書じゃない） | トップページ | ガストフロント/ダウンバーストって？（福井の突風） »
コメント