Fermi-Pasta-Ulam問題(FPU問題)で、振動子の動きをGIFアニメで見てみる。
非線形の再帰現象のさきがけとなったフェルミ・パスタ・ウラム問題。
http://oldweb.ct.infn.it/~rapis/corso-fsc2/FPU-recent-review-Ruffoetal05.pdf
一次元調和振動子を弱く非線形結合(2次や3次)して、最初に1つのモードにエネルギーを与えたら、全部のモードに分配されるかな?と思ってやってみたら最初のモードに戻ってしまったというもの。
例えば2次で結合する以下の方程式で、、、
d2xi/ dt2 = (xi+1 + xi-1 - 2xi) + α[(xi+1 - xi)2- (xi - xi-1)2]
α=1/4、N=64としてシンプレクティック8次でモード毎のエネルギーを計算(もちろんExcel VBAを使用)した結果がこれ。
最初のモードからいろんなモードにはエネルギーは行くものの、いつか戻ってくるという。
*FermiたちはMANIACというコンピュータでleap frog法で計算したそうです。
でこれはエネルギーの話。実際の振動子がそれぞれどう動いているか全くわからん。GIFアニメを作ってみた。
正弦波からだんだん気持ち悪くなっていく様子がよくわかる。ココログのアップロードファイルサイズ制限で長時間は表示できないので↓に静止画を。
各モードのピークになる時間でプロットしてる。
シンプレクティック8次のお話はこちら。
« せっかくの名曲が芸人さんの登場曲でだいなしになった例あれこれ、、、 | トップページ | 京都初詣シリーズ2013(その8 猿田彦神社) »
「学問・資格」カテゴリの記事
- 高周波(RF・マイクロ波・ミリ波・5G)関連ニュース2021年2月16日 IEEE Microwave Magazineの特集はオールデジタルのRFID、Microwave JournalはEバンド ミリ波通信に衛星や気球を使う話、アメリカの半導体企業がバイデンに投資を迫る、(2021.02.17)
- カオスを生じる電気回路、Chua’s circuitをLTspiceで回路シミュレーションしてみる。(2021.02.19)
- Labyrinth Chaos(迷宮カオス)を生むThomas-Rössler方程式のパラメータbを色々変えて、Python+Scipyでルンゲクッタ8次のDOP853(Dormand&Prince)を使って計算してGIFアニメ(2021.02.16)
- フィッツヒュー・南雲 (FitzHugh-Nagumo) 方程式をPython+Scipyでルンゲクッタ8次のDOP853(Dormand Prince)で計算。(2021.02.23)
- 「水晶振動子の等価回路計算」をカシオの高精度計算サイトkeisan.casio.jpの自作式としてUP! インピーダンスの大きさと位相がグラフ化できる。(2021.02.12)
« せっかくの名曲が芸人さんの登場曲でだいなしになった例あれこれ、、、 | トップページ | 京都初詣シリーズ2013(その8 猿田彦神社) »
コメント