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年1月16日 5G、レーダで使われるミリ波フェイズドアレイアンテナ各社まとめ、HuaweiはSAWまで取り込む、iPhone 11 ProがFBIにロック解除される、など。(2020.01.16)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2020年1月15日 iPhone 5G mmWのAiP、TSMCとともにASEにも発注と報道、Microwave Journal 1月号の特集はレーダ＆アンテナ、ResonantのWhitepaperなど。(2020.01.15)
- ”楕円の円周と一方の径からもう一方の径を求める。”をカシオの高精度計算サイトkeisan.casio.jpに自作式としてUP.(2020.01.13)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2020年1月13日 IMECがbeyond 5GのRF半導体技術（III-V/CMOS)について解説、ミンチークオはiPhone 5GモデルはSub6GHzもミリ波も同時に出ると語る、など。(2020.01.13)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2020年1月4日 Qorvoが6GHz対応のWiFiを支持＆5G RF for dummies 2nd edを公開、TSMCがiPhone mmW向けにInFO AiPの注文を受ける、Akouticsが軍用レーダ向けBAW出荷など。(2020.01.04)
« ランダウ・リフシッツ方程式（教科書じゃない） | トップページ | ガストフロント/ダウンバーストって？（福井の突風） »
コメント