物理を解説 ♪
ツイッター用のシェアボタン フェイスブック用のシェアボタン はてなブックマーク用のシェアボタン ライン用のシェアボタン
ツイッター用のシェアボタン フェイスブック用のシェアボタン はてなブックマーク用のシェアボタン ライン用のシェアボタン

方程式の一覧

参照しやすいように一か所にまとめておこうと思った。
作成:2020/9/29


あれ? 第 1 部まだ終わらないの?

基礎方程式をひと揃い集めるのにかなりの話数を費やしたので,どの式をどこで説明したのだったかが分かりにくくなってしまっていることだろう.探し回らなくてすむように,完成した方程式ばかりをここにまとめておくことにしよう.

ここまで書いてきて思うのは,初学者にとっての困難の一つは,同じ意味の数式に複数の表現方法があるというところだろう.ベクトルで表したり,成分で表したり,略記号で表したりする.それは著者の単なる好みの問題ではなくて,その場その場で利点や欠点があるものだから,どれかに統一して書き続けるというわけにもいかない.とにかく読者に慣れてもらうしかないのである.


連続の方程式

流体の質量保存則である.流れている流体が途中で消えたり現れたりしないことを表している.
(1)div(ρv) = ρt
この式は次のような意味である.
(2)(ρvx)x+(ρvy)y+(ρvz)z = ρt
これはアインシュタインの縮約記法で次のように表しても同じである.
(3)(ρvi)xi = ρt
この偏微分を展開して次のように書いてもいい.
(4)ρvixi + viρxi = ρt
右辺を左辺に持ってきて左辺第 2 項とまとめればラグランジュ微分の形になっているので,次のようにも書ける.
(5)ρvixi + DρDt = 0
左辺第 1 項を略記号で置き換えればまた少し印象が変わったりもする.
(6)ρdivv + DρDt = 0
読者の皆さんに期待するのは,これらの式が全て同一のものであるということに慣れてもらって,どの形で飛んできても即座に「連続の方程式だな」と把握してもらえるようになることである.難しい要求をしていることは分かっている.

急いで無理に暗記するのではなく,何度も式をいじっているうちに「こんなの当たり前じゃん何なら自力で変形するのも今すぐできる」と思えるくらいになってもらうのが良いと思う.


連続の方程式(密度が一定の場合)

密度ρが変化しない場合には連続の方程式は次のような簡単な形になる.
(7)divv = 0
これは次のような意味である.
(8)vxx+vyy+vzz = 0
縮約記法を使えば次のようにも書ける.
(9)vixi = 0



ラグランジュ微分

これはわざわざ書く予定ではなかったが,今後もよく使うので載せておこう.ラグランジュ微分の定義は次のようなものである.
(10)DADt  At + vxAx+vyAy+vzAz
これは流体と一緒に流れている人から見た,その人のいる地点での物理量Aの時間的な変化率を意味している.物理量Aを省いた次のような表記もよく使われる.
(11)DDt  t + vxx+vyy+vzz
右辺の最初の 3 つの項は速度ベクトルvとベクトル演算子 = (x,y,z)との内積の形になっているので,次のように表記することもある.
(12)DDt  t + v


運動量保存則

流体の運動量保存則の式は縮約記法を使って次のように書ける.
(13)(ρvi)t + (ρvivj)xj = fi + Tijxj
fiは流体の単位体積あたりに働く外力のi方向成分を意味している.Tijは応力テンソルである.

これを「運動量方程式」と呼ぶこともあるし,「コーシーの運動方程式」と呼ぶこともある.Tijに具体的な形を代入したものが「ナヴィエ・ストークス方程式」であるが,この段階でも「ナヴィエ・ストークス方程式」と呼ぶことがある.

(13) 式の左辺は「連続の方程式」を前提にすることでラグランジュ微分に書き換えることができるので,(13) 式は次のように表すこともできる.
(14)ρDviDt = fi + Tijxj
流体力学では「連続の方程式」は当たり前のように成り立っていることが前提となるので,(13) 式であろうが (14) 式であろうが,どちらを使っても構わない.

(12) 式のような表記を使えば次のように表すことも可能だろう.
(15)ρ(vit + vvi) = fi + Tijxj
先ほどから成分表示と縮約記法によって一つの式で簡潔に表しているが,実は次のような 3 成分の式であることは忘れないでいてもらいたい.
ρ(vxt + vvx) = fx + Txxxx + Txyxy + Txzxz(16)ρ(vyt + vvy) = fy + Tyxxx + Tyyxy + Tyzxzρ(vzt + vvz) = fz + Tzxxx + Tzyxy + Tzzxz
ここで出てきた (13) ~ (16) 式まで,全て同一の内容である.


ナヴィエ・ストークス方程式

ニュートン流体を仮定したときの応力テンソルは次のように表すことができる.
(17)Tij = pδij + λΘδij + μeij
μは粘性率,λは第二粘性率,Θdivvの略記,eijというのは具体的には次のような式である.
(18)eij = vixj + vjxi
この応力テンソルTijを運動量の保存則の式 (13) に代入したものが「ナヴィエ・ストークス方程式」である.普通は粘性係数は定数だと仮定して整理した次のような式を使う.
(19)(ρvi)t + (ρvivj)xj = fi  pxi + (λ+μ)Θxi + μvi
この式の左辺は (13) 式 ~ (15) 式の左辺のように,どのような形に表してもいい.

例えば (14) 式の左辺のような形を採用してやれば,ベクトル表記にも都合が良くて,次のように書くこともできる.
(20)ρDvDt = f  p + (λ+μ)Θ + μv


バロトロピー流体の仮定

流体の密度ρが圧力pだけで決まるという仮定ができる流体をバロトロピー流体と呼ぶのだった.
(21)ρ = ρ(p)
液体の密度は圧力によってほとんど変化しないが,そのようなものもバロトロピー流体の一種と考えていい.通常の気体の密度は温度Tにも依存するのでρ = ρ(p,T)であり,これは状態方程式である.(21) 式のような仮定は,特殊な制限化での状態方程式であると言える.

例えば等温条件下の理想気体では次のような関係が成り立つと考えてよい.Cは定数である.
(22)ρ = Cp
断熱条件下の理想気体では比熱比γを使って次のような関係が成り立つと考えてよい.Cは定数である.
(23)ρ = Cp1/γ


エネルギー保存則

通常は「連続の方程式」と「ナヴィエ・ストークス方程式」と「バロトロピー流体の仮定」の合計 5 つの式があれば未知関数の数と方程式の数が一致するので,それを基礎方程式として議論を進めればよい.常に成り立っているはずのエネルギー保存則のことは無視して構わない.

エネルギー保存則は流体の内部エネルギーなどを含んだ熱力学的な式として表されることになる.
(24)DEDt + pD(1/ρ)Dt = 1ρ(Q  divq + Φ)
Eは単位質量あたりの内部エネルギーである.Qは「加熱率」で,単位体積当たり,単位時間あたりに外部から与える熱エネルギーである.qは「熱流量」であり,熱伝導によって単位時間に単位面積を通過する熱エネルギーである.「フーリエの法則」によるとqは温度勾配に比例するので次のような式で表される.
q = kT(25)= (kTx, kTy, kTz)
比例係数kは「熱伝導率」である.

Φは「散逸関数」と呼ばれ,粘性が原因となって「単位体積」「単位時間」あたりに発生する摩擦熱である.それは次のような式で表される.
(26)Φ  λΘ2 + μ2eijeij
このeijは (18) 式の定義と同じである.

このように,エネルギー保存則を基礎方程式に取り入れようとすると未知関数の数が増えて,かえってややこしくなってしまう.しかしバロトロピー流体の仮定が使えない場合などには (21) 式のような条件に縛られない普通の状態方程式と一緒にして基礎方程式の一つとして使われることもある.




趣味の物理学書店

趣味で量子力学2の広告バナー