FPU(Fermi-Pasta-Ulam)問題、非線形振動子でどうやってエネルギーがモードに分配されるか?をシンプレクティック8次で計算[解説編]
先ほど 「FPU(Fermi-Pasta-Ulam)問題、非線形振動子でどうやってエネルギーがモードに分配されるか?をシンプレクティック8次で計算。エネルギーではなくて、実際の変位をGIFアニメにしてみた。あまり見たことないのでは?詳細はのちほどブログにて。」 というのをTweetした。
FPU(Fermi-Pasta-Ulam)問題、非線形振動子でどうやってエネルギーがモードに分配されるか?をシンプレクティック8次で計算。エネルギーではなくて、実際の変位をGIFアニメにしてみた。あまり見たことないのでは?詳細はのちほどブログにて。 pic.twitter.com/POQSixHgEE
— tomo (@tonagai) 2017年2月15日
ではその詳細を。この手の計算は昔いっぱいしてこのブログでも取り上げていたのですが、さっき
昔は計算担当の女性職員として名も残らないことが普通だったと思いますが、時代を経て取り上げられるようになったこともあります。
Fermi, Pasta, Ulam, and a mysterious ladyhttps://t.co/hcTrw1qpFT https://t.co/n8L7AEXAMP— Paul Painlevé@JPN (@Paul_Painleve) 2017年2月15日
というPainleve先生のTweetを見て再掲しようと思った。Nasaの女性”コンピュータ”を描いた話題の映画Hidden Figuresも作られたことだし。 では過去の記事の再掲:
非線形の再帰現象のさきがけとなったフェルミ・パスタ・ウラム問題。
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アニメにしたのが最初に書いたTweetでした。
シンプレクティック8次のお話はこちら。
« これが北朝鮮の毒針(ペンタイプ)か!2012年のCNNの記事より。金正男氏はこれで殺害されたのか、、、まさに007. | トップページ | ローム、お前もか! ソニー、パイオニアに続いてオカルトオーディオ、、、ワイヤボンドで音が変わる、、、半導体のバックグラインドで音が変わるといった人を思い出す。 »
「学問・資格」カテゴリの記事
- 高周波(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)
「日記・コラム・つぶやき」カテゴリの記事
- ExcelのLET関数+SEQUENCE関数で数値計算シリーズ(その5)ワンライナーでラマヌジャンの円周率公式を使って計算する。n=2で収束した。(2021.03.01)
- mRNAワクチンやCRISPR cas9のことが全然わからんので昔買ったブルーバックスの「アメリカ版 大学生物学の教科書」を3巻まで読む、、、と思ったら11年ぶりに完全改訂されたものが出てた!(2021.03.02)
- 新型コロナウイルス、中国、日本、韓国、アメリカ、ドイツ、フランス、イギリスでの感染者数を指数関数&ロジスティック関数&Log-Logプロットでべき関数フィッティングした(2/28更新)さすがにどの国も伸びは鈍ったか。(2021.03.01)
トラックバック
この記事へのトラックバック一覧です: FPU(Fermi-Pasta-Ulam)問題、非線形振動子でどうやってエネルギーがモードに分配されるか?をシンプレクティック8次で計算[解説編]:
« これが北朝鮮の毒針(ペンタイプ)か!2012年のCNNの記事より。金正男氏はこれで殺害されたのか、、、まさに007. | トップページ | ローム、お前もか! ソニー、パイオニアに続いてオカルトオーディオ、、、ワイヤボンドで音が変わる、、、半導体のバックグラインドで音が変わるといった人を思い出す。 »
コメント