Physic Labs

流体力学

ナビエ–ストークス方程式

ナビエ–ストークス方程式は、慣性、圧力、体積力、粘性応力を組み合わせ、粘性流体の運動量収支を記述します。

ナビエ–ストークス方程式は流体要素にニュートンの第2法則を適用します。構成則と仮定により形は異なります。以下は一様・非圧縮・粘度一定のニュートン流体に対する微分形です。

ρ(∂u∂t+(u⋅∇)u)=−∇p+μ∇2u+ρf,∇⋅u=0\rho\left(\frac{\partial\mathbf u}{\partial t}+(\mathbf u\cdot\nabla)\mathbf u\right)=-\nabla p+\mu\nabla^2\mathbf u+\rho\mathbf f,\qquad\nabla\cdot\mathbf u=0

定義: 運動量方程式の項

ρ\rho は密度、u\mathbf u は速度、pp は圧力、μ\mu は粘度、f\mathbf f は単位質量当たりの体積力です。左辺は局所加速度と移流加速度、右辺は圧力勾配、粘性拡散、体積力を表します。

平板間クエット流の速度分布と粘度の影響を調べます。これは単純な層流解で、ナビエ–ストークス方程式の一般解ではありません。

運動量収支

物質加速度は局所変化 ∂tu\partial_t\mathbf u と移流 (u⋅∇)u(\mathbf u\cdot\nabla)\mathbf u を含みます。低 ReRe では粘性が速度勾配をならし、高 ReRe では移流と慣性の重要性が増して流れの予測が難しくなりがちです。

標準解:クエット流

ux(y)=UyH,0≤y≤Hu_x(y)=U\frac{y}{H},\qquad 0\le y\le H

間隔 HH の平行平板で下板を固定し、上板を速度 UU で動かします。流れ方向の圧力勾配がないニュートン流体の定常層流速度分布は線形です。すべりなし境界条件により、流体速度は各壁面の速度に一致します。

例: すき間のクエット流速

間隔 H=2.0 mmH=2.0\,\mathrm{mm} の平板で上板が U=0.10 m/sU=0.10\,\mathrm{m/s} で動きます。クエット流の中央の速度を求めましょう。

解答

u(H/2)=U(H/2)/H=U/2=0.050 m/su(H/2)=U(H/2)/H=U/2=0.050\,\mathrm{m/s} です。

クイックチェック

各項は物理的意味をもち、単位は体積当たりの力です。hoDmathbfu/Dt ho Dmathbf u/Dt は慣性、−ablap- abla p は圧力による力、muabla2mathbfumu abla^2mathbf u は粘性による運動量拡散、homathbff homathbf f は体積力です。非圧縮流体では ablacdotmathbfu=0 ablacdotmathbf u=0 と組み合わせます。壁面でのすべりなし条件などの境界条件がなければ解は定まりません。

示された非圧縮方程式はどの流体モデルに基づいていますか。

例の平板間クエット流で中央の速度はどれですか。

参考文献

  1. Frank M. White (2016). Fluid Mechanics
  2. G. K. Batchelor (1967). An Introduction to Fluid Dynamics