非圧縮性流体

前回、の続きでRunge-Kutta法という天才的なアルゴリズムをみていきます。前回同様、関数を

$$
\begin{align}
&\frac{d u}{d t} =  f(t,u) \tag{1 ...

非圧縮性流体

前回、MAC法で方程式を導きました。ここで格子の考え方を理解しておかないと面倒なので、様々な格子について理解していきます。

MAC法での議論をするので、前回導いたMAC法の方程式を書いておきます。

$$

非圧縮性流体

非圧縮性流体力学の方程式はテンソル表記で下のように表せます。

$$
\begin{align}
&質量保存則 :\hspace{ 30pt } \frac{ \partial u_i }{ \partia ...

非圧縮性流体

次は陰解法で

プログラムは、以下の書籍の内容やプログラムを参考にして作ってます。プログラムは、一部改修・デバッグして使用しています。超いい本なので、ぜひ読んでください。

プログラムは朝倉書店のホームベーじからダ ...

非圧縮性流体

環境について

開発環境を整えます。

僕はWindowsはあまり触ったことないのと、そもそもパソコン5台(Macbook, MacbookPro, 自作Ubuntu, ラズパイ2台)も持ってるのにWindowsを ...

非圧縮性流体

まずは陽解法で

プログラムは、以下の書籍の内容やプログラムを参考にして作ってます。プログラムは、一部改修・デバッグして使用しています。超いい本なので、ぜひ読んでください。

プログラムは朝倉書店のホームベーじから ...

非圧縮性流体

まずは簡単なやつから

プログラムの計算では、以下の書籍の内容やプログラムを参考にして作ってます。プログラムは、一部改修・デバッグして使用しています。超いい本なので、ぜひ読んでください。

プログラムは朝倉書店のホ ...

非圧縮性流体

熱エネルギー保存則

前回で「非圧縮非等温流体」の運動量保存則に修正しました。一応書いておきます。

$$
\begin{align}
&\frac{ \partial u_i }{ \pa ...

非圧縮性流体

熱流体の扱い方

前回でやっと、「非圧縮等温流体」の運動量保存則が完成しました。一応書いておきます。

美しいですね。流体の挙動を表す式が、こんなシンプルな方程式にまとまるなんて。これを作ったNavierさんとSt ...

非圧縮性流体

重力項の消去

引き続き、流体の運動量保存則とはどのようなものになるかを考えていきます。やっとです。やっと今回で完結です。

前回、運動量保存則はここまで進んだのでした。

$$
\begin{a ...