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についての発展はこんなのが出てた。
« ランダウ・リフシッツ方程式（教科書じゃない） | トップページ | ガストフロント/ダウンバーストって？（福井の突風） »
「学問・資格」カテゴリの記事
- ミリ波パッチアンテナは比誘電率が高い方が効率が上がる？、、、ということを最近何人かに聞かれたのでちょっと説明。(2021.11.30)
- 購入回数の確率を求めるNBDモデル（負の二項分布）の確率分布を与えるパラメータKの計算をカシオの高精度計算サイトkeisan.casio.jpにUP!(林先生の森岡毅さんインタビュー見て数学マーケティングというのに興味持って確率思考の戦略論を買って読んだので）(2021.11.26)
- 購入回数の確率を求めるNBDモデル（負の二項分布）の計算をカシオの高精度計算サイトkeisan.casio.jpにUP!(林先生の森岡毅さんインタビュー見て数学マーケティングというのに興味持って確率思考の戦略論を買って読んだので）(2021.11.25)
- 「円の弧長・弦長・矢高・半径のどれか2つを与えて残りを計算」するWindowsアプリ(.exe)を作ってみた。keisan.casio.jpとExcel VBA版のVB.net移植版。(2021.11.24)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2021年11月16日 IEEE Microwave Magazineで金星からの通信、Microwave JournalではRF GaN、DuPontがRogers（低損失PCB）を買収、など。(2021.11.16)
« ランダウ・リフシッツ方程式（教科書じゃない） | トップページ | ガストフロント/ダウンバーストって？（福井の突風） »
コメント