Euler方程式

Dr. SSS 2022/09/17 - 20:35:24 流体力学
はじめに

理想流体の運動を記述するEuler方程式の導出を行う。


keywords: Lagrange微分, 理想流体, 対流微分, 粘性

内容

理想流体

流体は,その各要素が,境界面を通して様々な方向に力を及ぼし,結果として要素間での摩擦(粘性)や熱交換が生じる。 境界面を通した面積当たりの力は,ストレス(応力)テンソル(stress tensor)と呼ばれるテンソル量$\bm{\sigma}$によって表現できる。 成分で$\sigma_{ij}$と表したとき,添え字${ij}$は,$i$方向に垂直な面を通して$j$方向に作用する成分であることを指している。

ここでは,ストレステンソルが等方成分のみを持ち,その他の成分に対応する粘性は無視できる,理想流体(ideal fluid)あるいは完全流体(perfect fluid)と呼ばれる流体を考える。 この場合,仮定より,ストレステンソルはスカラー圧力$p$を用いて

\begin{align} \bm{\sigma}=p\bm{I} \end{align}

のように表せる。 ここで$\bm{I}$は3次元の単位テンソル

\begin{align} \bm{I} = \left( \begin{array}{ccc} 1&0&0 \\ 0&1&0 \\ 0&0&1 \end{array} \right) \end{align}

である。

すると,圧力$p$を持つある体積に働く力は,その境界上の積分により

\begin{align} -\oint p\bm{I}\cdot d\bm{S} \end{align}

で与えられる。 ここで,面積要素$dS$に対して外向きの法線方向を正に取っていることに注意しよう。 これにGaussの定理を適用して体積積分に変えると

\begin{align} -\oint \nabla \cdot (p\bm{I})dV = -\int \nabla pdV \end{align}

が得られる。 つまり,流体中の単位体積には$-\nabla p$という力が働くことが分かる。

この事実を用いて流体の運動方程式を構築してみる。 古典的な流体を考えているため,流体は古典的な質点の集まりとみなせる。 よって,各流体要素もNewtonの運動方程式と同様の形の運動方程式に従うことが予期される。 古典的粒子に働く力に相当する量は今得られたから,後は単位体積当たりの質量,すなわち質量密度$\rho_m$に,流体要素のフロー速度$\bm{u}$の時間微分をかけてイコールの関係で結ぶことで

\begin{align} \rho_m \frac{d\bm{u}}{dt} = -\nabla p \end{align}

が得られる。



対流微分

大きさや形のない質点の場合と異なる流体の運動の性質として,要素ごとに速度が異なりうることから,速度が一般に位置$\bm{x}$に依存する関数であるということが挙げられる。 つまり,例え各点ごとに流体の流れが定常であったとしても,流れが一様であるとは言えないから,異なる流体要素間で速度を比較した場合,その差は$0$とは限らない。 こうした速度の位置依存性や,それを通した時間依存性は$\bm{u}=\bm{u}(\bm{x}(t),t)$と表せる。 よってこの時間微分は

\begin{align} \label{eq:fluid_derivative} \frac{d\bm{u}}{dt} = \left(\frac{\pd\bm{u}}{\pd t}\right)_{\bm{x}} + \frac{d\bm{x}}{dt}\cdot \left(\frac{\pd}{\pd \bm{x}}\right)_{t}\bm{u} \end{align}

となる。 陽に書いてある通り,右辺1項目が位置を固定して測定した速度変化率であり,2項目が時間変化を無視し,流れに沿った位置の変化に伴って測定される速度変化率を表している。 流体の流れに沿った位置の時間変化とは,すなわち$\bm{u}=d\bm{x}/dt$であるから,結局(\ref{eq:fluid_derivative})は

\begin{align} \label{eq:convective_derivative} \frac{d\bm{u}}{dt} = \frac{\pd\bm{u}}{\pd t} + (\bm{u}\cdot\nabla)\bm{u} \end{align}

と表せる。 ここに現れる,流体の流れに沿った時間変化率を与える微分演算子は,Lagrange微分(Lagrange derivative)物質微分(material derivative)対流微分(convective derivative)など,様々な呼び方をされるもので

\begin{align} \frac{D}{Dt} \equiv \frac{\pd}{\pd t} + (\bm{u}\cdot\nabla) \end{align}

の左辺のように,特別な記号を用いて表させることも多い。

Euler方程式

表現(\ref{eq:convective_derivative})を用いれば,理想流体の運動方程式は

\begin{align} \frac{\pd\bm{u}}{\pd t} + (\bm{u}\cdot\nabla)\bm{u} = -\frac{1}{\rho_m}\nabla p \end{align}

と書き換えられる。 その他,重力などの外力が加わる場合は,それによる加速度$\bm{g}$を右辺に加え

\begin{align} \frac{\pd\bm{u}}{\pd t} + (\bm{u}\cdot\nabla)\bm{u} = -\frac{1}{\rho_m}\nabla p +\bm{g} \end{align}

となる。 理想流体の運動を記述するこの方程式を,Euler方程式(Euler equation)という。


参考文献


自然科学に関する質問やノート作りは『AfterSchool』で。