プログラミングの記事を書きます。偏微分方程式を数値計算で解いて、結果をアニメーションしたいと思ったときにPythonを使うと楽でしたので記事に残しておこうと思います(^^)/本記事の内容はこちらFortranで出した出力データをPythonでアニメーション作成する前回の記事で、アニメーション作成にはgnuplotやparaviewを使ったりしました。【Fortranコードあり】(前回の記事)Fortranで数値計算→paraviewでvtkファイルを読み込みアニメーション作成Fortranで解いた問題をParaviewを使ってvtkでアニメーション作成(前回の記事)Pythonで数値計算してPythonでアニメーション作成Pythonで2次元の熱伝導方程式をアニメーションが簡単にできた色々試してみて・・・・Fortranで数値計算してデータを出力しておいて、Pythonでサクッっとアニメーション作成作れるようにもしておこうかなって思いました。【本記事の内容】数値計算からアニメーション作成までのフローはこちら↓。※Fortranのコードとスクリプトは本記事にアップされています。※FortranのコードはFortran90/95の書き方にしているつもりですが正確に書けているかは不明です。※FortranはWSLを使ってインストールしています。Notebookを使っています。x}=0\tag{1}\end{align*}※\(u\):一定速度1次元の移流方程式は、下記の絵のように何かの物理量を形を変えずに一定速度\(u\)で運ぶ方程式を意味しています。1次元移流方程式を離散化する(1)式を離散化する場合は注意する必要があったんですよね。1次元移流方程式を陽解法で解くときには、数値的に安定な条件で離散化&条件設定をする必要がありました。【Fortranコードあり】カマキリこの記事は、なかなか教育的な記事が書けた気がする。陽解法で解くのあれば後退差分(風上差分)でないと数値的に安定には解いてくれないということがわかりました。↑Fortranで解くのは赤枠の離散化した式です。空間と時間の刻みにも制限がありました。thermal_explicitファイルの書き出し部分の解説↓Fortranを実行すると・・・(僕はgfortranを使っています)データがいっぱいできるので、これをPythonでアニメーション作成しようということになります。Nootebookで実行すると・・・アニメーション作成ができます(^^)/1次元移流方程式のアニメーション結果https://takun-physics.net/wp-content/uploads/2019/05/20190520_231631-2.mp4カマキリやったー。-【理想気体の比熱の温度依存性(8)】量子論で2原子分子の比熱を計算。低温極限で調和振動子の量子化エネルギーは凍結している。【CAE解析の目的とは?】CAEは実機の代替えなのか?精度を無視してパラメータスタディすることは意味がないのか?【Udemyがお勧め理由5つ】プログラミングを勉強するなら動画学習おすすめ【まとめ】比熱の温度依存性について。古典論と量子論から理解する。コメントを残すコメントをキャンセルプロフィールカマキリ(^^)大学以上の物理を勉強してブログに書いていきます。大学の専攻は物性理論で、Fortranを使って数値計算をしていました。CAE解析(流体解析)・毎日note更新しています。内容は、「現在勉強中の内容」「日々思ったこと」を書いています。1DCAEの操作方法などを動画にしています。youtube・Qiitaではプログラミング言語の基本的な内容をまとめています。