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. | トップページ | ローム、お前もか！ ソニー、パイオニアに続いてオカルトオーディオ、、、ワイヤボンドで音が変わる、、、半導体のバックグラインドで音が変わるといった人を思い出す。 »
「学問・資格」カテゴリの記事
- FIFAワールドカップカタール2022で今使われているボールはAdidasとKinexon共同開発のUWB＋慣性センサ内蔵で位置が高精度にわかる。ICはQorvo(Decawave)のDW1000を使ってる。タグはアンテナもリファレンスそのまま、アンカーは外部PA追加みたい。(2022.11.30)
- 「人の体の1日の水分の出入り（代謝回転量）の予測式」をカシオの高精度計算サイトkeisan.casio.jpに自作式としてUP!サイエンスに論文が出たということで。(2022.11.29)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース(11/21) IEEE Microwave Magazineにこれまでの無線通信についてのレビュー、Microwave Journalは6Gに向けたパワーアンプ、Bluetooth SIGが6GHz帯を検討、電波のバンドがめちゃくちゃややこしい件の解説。(2022.11.21)
- Raspberry PI 3を数年ぶりに触る。Mathematicaをインストールすると最新版の13.1が入った！新機能で分数階微分とその微分方程式が解けるようになってるので例題通りやってみた。(2022.11.10)
- 先延ばしの普遍的法則（課題提出、夏休みの宿題、進捗どうですか、などなど）N(t)=M/(D-t+C)-M/(D+C)を導出、カシオの高精度計算サイトkeisan.casio.jpにUP!締め切りギリギリじゃないとみんなやらない感覚にあってるな。(2022.11.05)
「日記・コラム・つぶやき」カテゴリの記事
- FIFAワールドカップカタール2022で今使われているボールはAdidasとKinexon共同開発のUWB＋慣性センサ内蔵で位置が高精度にわかる。ICはQorvo(Decawave)のDW1000を使ってる。タグはアンテナもリファレンスそのまま、アンカーは外部PA追加みたい。(2022.11.30)
- 「人の体の1日の水分の出入り（代謝回転量）の予測式」をカシオの高精度計算サイトkeisan.casio.jpに自作式としてUP!サイエンスに論文が出たということで。(2022.11.29)
- 新型コロナウイルス、日本の陽性者数＆ワクチン接種者数総計をプロット＆中国、韓国、アメリカ、ドイツ、フランス、イギリスの陽性者数もプロット（11/27更新) 日本も急増しているが、それ以上に中国の増え方がすごくなってる。韓国もじわじわ増えてきてる。(2022.11.27)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース(11/21) IEEE Microwave Magazineにこれまでの無線通信についてのレビュー、Microwave Journalは6Gに向けたパワーアンプ、Bluetooth SIGが6GHz帯を検討、電波のバンドがめちゃくちゃややこしい件の解説。(2022.11.21)
- 新型コロナウイルス、日本の陽性者数＆ワクチン接種者数総計をプロット＆中国、韓国、アメリカ、ドイツ、フランス、イギリスの陽性者数もプロット（11/13更新) 日本は完全に増加に転じた。人口比18.4%まで到達。中国と韓国はもう公式サイトで陽性者数発表しなくなった。。。(2022.11.13)
トラックバック
この記事へのトラックバック一覧です： FPU(Fermi-Pasta-Ulam)問題、非線形振動子でどうやってエネルギーがモードに分配されるか？をシンプレクティック8次で計算[解説編]:
« これが北朝鮮の毒針（ペンタイプ）か！2012年のCNNの記事より。金正男氏はこれで殺害されたのか、、、まさに007. | トップページ | ローム、お前もか！ ソニー、パイオニアに続いてオカルトオーディオ、、、ワイヤボンドで音が変わる、、、半導体のバックグラインドで音が変わるといった人を思い出す。 »
コメント