京都旅行3日目、最終日。 最終日は、奈良方面へ行ってみることにした。奈良といえば、思いつくメジャーな観光スポットがいくつかあるが、まず法隆寺に行った。聖徳太子の時代に建立されたこともあって、ある意味、京都の有名どころの神社仏閣よりも地味であ…

4/29京都旅行2日目。2日目は太秦の東映映画村から。実際の映画のセットや楽屋も見学できるかなと期待していたが、あまり実際に使われているものは少なく、アトラクション的なものが中心であった。 東映映画村の後は、亀岡までJR嵯峨野線で上がっていき、保津…

このGW前半（4/27〜4/29）に2泊３日で京都に旅行に行ってきた。 まず、初日（4/27）から。 京都駅に13:30頃に到着して、河原町のホテルに移動し14:00 頃にホテルにチェックイン。 チェックイン後は、四条大橋や祇園周辺を散策しながら、三十三間堂まで歩いて…

gmshインストールメモ これは、gmshをインストールした際のメモである。 インストール環境 OS： macOS Montery（12.0.1） CPU： Apple M1 インストーラーの取得&インストール gsmhの公式ホームページ中のmacOS(ARM)をダウンロード。 gmsh公式ホームページ ht…

概要・背景 物体周りの流れの解析を定常解析で実施する場合、特に形状が複雑な場合、目標残差に収束せず計算が何ステップにわたり継続し、計算がなかなか終了できない状況も起こりうる。 そのような場合、残差だけではなく、評価したい物理量の履歴をプロッ…

概要 本記事は、OpenFOAMのsimpleFomaソルバー（定常乱流場解析）の残差をプロットするPythonスクリプトを作成したメモである。 使用環境 使用環境は以下になる。 - OpenFOAM ver2312 - Python 3.10 OpenFOAMにおける設定 simpleFOAMの計算において、logファ…

概要 OpenFOAMでは、化学反応や燃焼を計算するソルバーとして、reactingFoamが用意されている。本記事は、reactingFoamの基本的な機能を確認した際のメモである。 実施した例題 reactingFoam内にあるチュートリアルのうち、以下の例題を実施した。 - counter…

概要 昨年の5月にOpenFOAM v2112をインストールをしていたが、その後更新されたバージョンのインストールが滞っていた。 これは、最新バージョン（2024年3月30日現在）をインストールした際のメモである。 前回（2022年5月11日）同様に、以下のGabriel Gerle…

前回、ロジスティック方程式について、数値解法を自作の関数を作成して解いた。 akh1116.hatenablog.com 今回は、scipyライブラリのintegrateモジュールを用いてコーディングをし、解いてみる。 自作で関数を作成するより、モジュールに渡す引数に注意すれば…

今回は、Pythonによる常微分方程式を解く際によく使われるscipyライブラリのメモである。 Scipyによる常微分方程式 Pythonを用いて常微分方程式を解く際は、ライブラリを用いるのが最も手軽で、かつ高精度な結果が得られる。 常微分方程式を解くライブラリと…

常微分方程式の数値計算例題-ロジスティック方程式 ここでは、常微分方程式の数値計算の例題として、ロジスティック方程式を扱ってみる。 ロジスティック方程式 ロジスティック方程式は、生物の個体数の増加を表すモデルであり、以下のような形の方程式とし…

本稿から具体的な数値計算の手法についての話に入る。 まずは、古典力学の基本方程式であるNewtonの運動方程式を数値的に解く手法から。 古典力学と常微分方程式 古典力学は、基本的には質点の運動はNewtonの運動方程式で記述される。 Newtonの運動方程は、…

Matplotlib基本事項メモ ここでは、数値計算結果をグラフなどの可視化の際に便利なMatplotlibライブラリの基本的事項を記しておく。 Matplotlib Matplotlibはグラフ描画用のライブラリである。数値計算の結果を簡易的に確認するだけではなく、論文やプレゼン…

現在、以下の書籍で数値計算（特に物理に関する）におけるPythonの活用を学習している。 Pythonによる計算物理作者:大槻純也森北出版Amazon まずその１として、Pythonで数値計算コードを書く際に威力を発揮するNumpyライブラリごく基本的な事項のメモである…

目的 この記事は、MacでC言語の開発環境を整備した際のメモである。 コンパイラのインストール MacではXcodeをインストールすることで、Mac専用のC言語開発環境であるclangが自動的にインストールされる。 xcodeのインストールについては、ターミナルを開い…

函館旅行から大晦日に帰ってきて、年越しは昨年同様自宅でゆっくり過ごした。 旅行先での年越しも非日常感があり、良いものであるが、やっぱり自宅が一番落ち着く。 天気が良かったので、湘南海岸まで散歩に出かけて、帰り途中に近くの龍口寺で初詣をしてき…

2023年の仕事納めは28日で、翌日の29日〜31日の2泊3日の旅程で函館へ行ってきた。 29日は妻が仕事納めで、就業後直で東京駅19:20の最終の新幹線で函館へ向かい、函館の宿泊ホテルにチェックインしたのは日付が回っていたこともあり、初日は移動のみ。29日の…