Skip to content

A1.ナビエ・ストークス方程式

物体の運動方程式はma=Fであった。 流体の運動方程式は、ナビエ・ストークス方程式と呼ばれ、以下の形を持つ。

vt+(v)v= 1ρp+μρ2v+F

ここでは、この方程式のかんたんな導出と、円筒座標系への変換について書いておく。

流体素片の運動方程式の導出

いま、とある流体素片があったとして、この流体素片の運動方程式を導出してみよう。 流体素片の体積をδVとしておく。すなわち、この流体素片の質量は

()=ρδV

となっている。この流体素片の加速度はラグランジュ微分を用いてdvdt と与えられそうなので、運動方程式ma=Fは、

ρδVdvdt=δFbody+δFsurface

と書ける。ここで、右辺のちからを2つに分けた。右辺第一項は、流体素片の全体に掛かる力である(体積力)。例えば重力など。そのため、この力はかんたんに

δFbody=ρδVF

と書けるであろう。これは質点の運動方程式と変わらない。 問題になるのは第二項、表面力である(連続体の力学では面積力とか言うかも)。この力は、流体素片の表面にかかる力だとしよう。

ここで、流体素片の面積要素をベクトルdS(太字)とおくと、おそらく表面力は面積要素の大きさに比例すると思われるので、二階のテンソルを用いて

(δFsurface)i=PijdSij

とかけそうな気がしてくる。これを流体素片の表面全体で積分するとき、体積積分に直すときは一回微分しておけばよかったので

(δFsurface)i=PijxjdV

これを元の式に戻すと流体素片の運動方程式は、i成分を取り出して書くと

ρdvidt=ρFiPijxj

という形に書ける。

表面力の分解 - 圧力勾配力と粘性項

ここで、この右辺第二項を、圧力項と粘性項に分けて書いてみよう。クロネッカーのデルタを思い出して、対角成分だけ別に書くと

Pij=pδij+πij

と分けることができる。

もし、第二項πijが0であれば、この場合は、表面力は全てその法線方向に働く。これを理想流体とよんでいる。 だが我々は、日常的な経験の中で、流体の抵抗力を知っている。例えば「ねばりけ」というのは、とある流体素片が動いているときに、そのそばにいるが、速度ベクトルの先にはいない別の流体素片に対してなにかしら力が作用している、ということである。

…ちょっと例を出そう。流れ場にシアがあった場合、πijはどんな形を持つだろうか。

(途中。P79あたり。5Viscous Flows)

圧力項

どこでもドアを月につないだらどうなるだろう?たぶん地球の空気が凄まじい勢いで月に噴出するであろう。これが圧力勾配力である。 何も難しいことは言っておらず、もし圧力に空間分布があり、勾配があれば、すなわち圧力pのグラジエントに比例して力がかかる効果である。 圧力は

p=ρcs2

だったことを思い出すと、原始惑星系円盤では基本的に圧力は動径方向には内側ほど高く、鉛直方向には中心面ほど高い。 すなわち、圧力勾配力は外側に押し出すように働いている。

粘性項

(まだ書いてない)

円筒座標系への変換

まず、円筒座標系における発散等をおさらいしよう。 ψをスカラー場、A=Are^r+Aθe^θ+Aze^zをベクトル場とすると、下記が成り立つ。

ψ=ψre^r+1rψθe^θ+ψze^z,A=1rr(rAr)+1rAθθ+Azz,×A=(1rAzθAθz)e^r+(ArzAzr)e^θ+1r[r(rAθ)Arθ]e^z,2ψ=1rr(rψr)+1r22ψθ2+2ψz2.

(人生で一度くらいは導いてみてもいいかもしれない)。

Chouduri, Appendix C1, P393

これをナビエ・ストークス方程式に当てはめると、円筒座標系におけるナビエ・ストークス方程式が得られる。

vrt+vrvrr+vθrvrθ+vzvrzvθ2r=1ρpr+ν(2vrr2+1r22vrθ2+2vrz2+1rvrr2r2vθθvrr2)+Fr,vθt+vrvθr+vθrvθθ+vzvθz+vrvθr=1ρrpθ+ν(2vθr2+1r22vθθ2+2vθz2+1rvθr+2r2vrθvθr2)+Fθ,vzt+vrvzr+vθrvzθ+vzvzz=1ρpz+ν(2vzr2+1r22vzθ2+2vzz2+1rvzr)+Fz.

Chouduri, Appendix C2, P393

円筒座標系での連続の式

ついでだし、円筒座標系における連続の式もメモしておこう。 ベクトルの形で書いた連続の式

ρt+(ρv)=0

に対して、先程の円筒座標系におけるをただ代入すると

ρt+1rr(rρvr)+1rθ(ρvθ)+z(ρvz)=0

となる