テスラノートテスラノート物理シミュレーションとAIとその応用。非圧縮性流体2階微分の差分近似まず、前回一階微分に対してしかやっていなかったので、2階微分に対しての差分近似から始めます。\qquad\end{align}$$拡散項は2階微分でした。いうことで、2階微分に対する差分近似を考えてみます。仮想的に上のようなメッシュを考えてみます。メッシュの中心に値が入っているという、いわゆるスタガード格子と呼ばれるものです。i+1/2、i-1/2の点があると考えます。すると、以下のように考えて、差分近似をすることができます。\label{eq:scale-factor-3}\end{align}$$よしいけた。1/2の部分を考えるなんて怪しげなんて思うでしょうか。こんなんじゃ納得できないぃ。いう人のために、Taylor展開でやっておきます。考え方は前回と同じです。1階微分を消して、2階微分を残せばいい。\label{eq:scale-factor-7}\end{align}$$お、最初から2次精度です。最初から2次精度であり、次は4次精度なので、まあほとんどこの中心差分近似が使われるが、一応他にもあります。導き方は前回の記事の3次精度や4次精度と同じなので、暇だったらやってみるといいかもしれません。\label{eq:scale-factor-10}\end{align}$$2階微分と拡散2階の偏微分というのは、しょっちゅう目にすると思います。2階偏微分を目にしない日はないと言っていいくらいです。\label{eq:scale-factor-11}\end{align}$$これはみんな大好きポアソン方程式だー。圧力のとこで出できます。電磁気でもひたすら見かけます。右辺のfが0になっているのがラプラス方程式です。これらは基本的に値が広がっていくような現象を表しています。ポアソン方程式は電磁気のイメージがあると思いますが、右辺の値が湧き出して広がっていくイメージで、ラプラス方程式は右辺が0なので、均一に広がって馴染んでいくイメージです。他にも、熱の拡散方程式なんてがあります。どれも拡散っぽいイメージがあります。?波動方程式なんかはちょっと違うイメージですが。熱の拡散方程式を考えてみましょう。こやつが一番わかりやすいです。物理現象としては、個体の熱が隣に伝わっていって、均一になっていくというものです。\label{eq:scale-factor-12}\end{align}$$右辺はベクトル的に書けば、divであり、以下のように縮約を取ることに注意しなければならなりません。\label{eq:scale-factor-13}\end{align}$$これは、まずフーリエの法則から導き出せます。フーリエの法則とは、次のようなものです。\label{eq:scale-factor-14}\end{align}$$Jは熱流速密度[W/m2]、αは熱伝導率[W/(m・K)]です。この法則は何を意味しているか分かるでしょうか?意味しているのは非常に簡単なことです。簡単の為に1次元で考えてみます。\label{eq:scale-factor-15}\end{align}$$この式が意味していることは簡単です。位置によって熱の勾配があると、勾配の逆方向に熱が流れ込みます。熱が高い方から低い方へ移動していくということです。x}\)は正の値になり、この勾配の逆方向に熱伝導率をかけた割合で熱が流入するということになります。次に、次図のような微小領域が蓄える熱量を考えてみます。この微小空間の熱収支を考えてみます。\label{eq:scale-factor-16}\end{align}$$のようになります。何が言いたいかというと、熱伝導において空間の2階微分が示すものはどれだけ熱が残るのかということを示しています。熱の収支こそがその場所の熱の変化量で、それを個々に考えていったものが熱拡散方程式というわけです。熱はフーリエの法則により伝わるんだけど、一瞬で伝わるわけでもないし、じゃあどういう風に伝わっていくの?というのを示すのが熱拡散方程式というわけです。流体の拡散項(=粘性項)も同じで、流体は速度差があると分子粘性によって速度が渡されていき、つまり運動量が熱のように高い方から低い方へ伝わっていくような現象が起こり、拡散していくのです。それを表すのが拡散項であり、拡散項は別名で粘性項と呼ばれます。熱の場合は係数が熱伝導率でしたが、流体の場合は動粘性係数になります。余談ではありますが、熱伝導率\(\alpha\)は必ず正だということを強調しておきます。これは物理現象を考えれば分かることで、もし熱伝導率が負だと、熱が低い方から高い方へと流れていくことになります。当たり前のことですが、認識だけしておきましょう。流体力学の重要な定理にベルヌーイの定理というものがあります。プログラムは、以下の書籍の内容やプログラムを参考にして作ってます。コメントをキャンセルメールアドレスが公開されることはありません。