物質微分

物質微分（ぶっしつびぶん、英: material derivative）とは流れに乗って移動する流体粒子の物理量 (温度や運動量)の時間変化率のことで、連続体力学の概念の一つである。固定された場所での物理量の時間変化でなく、流れに乗って動く仮想的な「観測者」が観た物理量の時間変化を記述する。 物質微分はラグランジュ描像に基づく時間変化をオイラー描像に基づく時間変化で記述したものである。物体固有の時間変化を記述するものなので物質微分 D / D t {\displaystyle \mathrm {D} /\mathrm {D} t} は偏微分 ∂ / ∂ t {\displaystyle \partial /\partial t} と違いガリレイ不変（英語版）である。 名称としては他に、物質時間微分、流れに乗って移動するときの微分、実質微分、ラグランジュ微分などとも呼ばれる。 速度場 v {\displaystyle {\boldsymbol {v}}} の流れにおける、スカラー場 φ ( x , t ) {\displaystyle \varphi ({\boldsymbol {x}},t)} およびベクトル場 A ( x , t ) {\displaystyle {\boldsymbol {A}}({\boldsymbol {x}},t)} の物質微分は以下のように表される。 D φ D t = ∂ φ ∂ t + v ⋅ ∇ φ {\displaystyle {\frac {\mathrm {D} \varphi }{\mathrm {D} t}}={\frac {\partial \varphi }{\partial t}}+{\boldsymbol {v}}\cdot \nabla \varphi } D A D t = ∂ A ∂ t + v ⋅ ∇ A {\displaystyle {\frac {\mathrm {D} {\boldsymbol {A}}}{\mathrm {D} t}}={\frac {\partial {\boldsymbol {A}}}{\partial t}}+{\boldsymbol {v}}\cdot \nabla {\boldsymbol {A}}} ここで、それぞれの式の右辺第2項を移流項、対流項と呼び、非一様な物理量の分布の中を移動したことで観測される物理量の変化率を表す。 スカラー場 φ ( x , t ) {\displaystyle \varphi ({\boldsymbol {x}},t)} の物質微分は直観的には流れに乗って動く物体から見た場合における φ ( x , t ) {\displaystyle \varphi ({\boldsymbol {x}},t)} の変化率を表す。（ベクトル場の場合も同様）。 実際、位置のグラフ u ( t ) {\displaystyle {\boldsymbol {u}}(t)} で記述される質点の軌道は速度場 v ( x , t ) {\displaystyle {\boldsymbol {v}}({\boldsymbol {x}},t)} にそっているので、 d u d t = v ( u , t ) {\displaystyle {\frac {\mathrm {d} {\boldsymbol {u}}}{\mathrm {d} t}}={\boldsymbol {v}}({\boldsymbol {u}},t)} となるから(なお、この性質を満たす u ( t ) {\displaystyle {\boldsymbol {u}}(t)} を流跡線という)、ライプニッツ則から d φ ( u , t ) d t = ( ∇ φ ) ⋅ d u d t + ∂ φ ∂ t = v ⋅ ∇ φ + ∂ φ ∂ t = D φ D t {\displaystyle {\frac {\mathrm {d} \varphi ({\boldsymbol {u}},t)}{\mathrm {d} t}}=(\nabla \varphi )\cdot {\frac {\mathrm {d} {\boldsymbol {u}}}{\mathrm {d} t}}+{\frac {\mathrm {\partial } \varphi }{\mathrm {\partial } t}}={\boldsymbol {v}}\cdot \nabla \varphi +{\frac {\mathrm {\partial } \varphi }{\mathrm {\partial } t}}={\frac {\mathrm {D} \varphi }{\mathrm {D} t}}} が成立する。 上では一粒子しかない場合を想定したが、初期時刻における位置 X {\displaystyle {\boldsymbol {X}}} でパラメトライズされた粒子の族 u ( X , t ) {\displaystyle {\boldsymbol {u}}({\boldsymbol {X}},t)} を考えた場合も、 X {\displaystyle {\boldsymbol {X}}} を固定して同様の証明を行う事で、同様の式が導ける。更に粒子の族を連続体にまで拡張したものが物質微分である。 以上の説明から分かるように、物質微分 D φ D t {\displaystyle {\frac {\mathrm {D} \varphi }{\mathrm {D} t}}} は物質に固定して観測する座標系（物質表示）における時間微分を表すが、それに対し通常の偏微分 ∂ φ ∂ t {\displaystyle {\frac {\partial \varphi }{\partial t}}} は空間上に固定された座標系（空間表示）における時間微分であるといえる。 定常流はすべての物理量のオイラー描像的時間変化率が ∂ / ∂ t = 0 {\displaystyle {\partial /\partial t}=0} となる流れであるが、ラグランジュ描像的時間変化率が D / D t = 0 {\displaystyle {\mathrm {D} /\mathrm {D} t}=0} となるとは限らないことに注意すべきである。 一つの流線に着目する。流線上のある点からの道のりを s {\displaystyle s} 、流線の単位接ベクトルを e s {\displaystyle {\boldsymbol {e}}_{s}} と表す。速度ベクトルは流線に接しているので、定常流における物質微分は D D t = ∂ ∂ t + v ⋅ ∇ = 0 + ( v e s ) ⋅ ∇ = v ∂ ∂ s ( ∵   e s ⋅ ∇ = ∂ ∂ s ) {\displaystyle {\begin{aligned}{\mathrm {D} \over \mathrm {D} t}&={\partial \over \partial t}+{\boldsymbol {v}}\cdot \nabla \\&=0+(v{\boldsymbol {e}}_{s})\cdot \nabla \\&=v{\partial \over \partial s}&(\because ~{\boldsymbol {e}}_{s}\cdot \nabla ={\partial \over \partial s})\end{aligned}}} となり、流線方向の変化率に速さをかけたものに等しいことが導かれる。 これから、定常流( ∂ / ∂ t = 0 {\displaystyle {\partial /\partial t}=0} )でも、流線に沿って物理量が変化するなら D / D t ≠ 0 {\displaystyle {\mathrm {D} /\mathrm {D} t}\neq 0} であることがわかる。 応用で重要なのは速度の物質微分すなわち加速度である。定常流、つまり、速度の時間変化がない流れでも、流体粒子の加速度は0とは限らない。定常流でも、流線に沿って速度の大きさは変化しうるし、流線に沿って速度の方向が変わる（流線が曲がる）こともありうる。これを式に表すと、 D v D t = ∂ ∂ s ( v 2 2 ) e s − v 2 R e r {\displaystyle {\mathrm {D} {\boldsymbol {v}} \over \mathrm {D} t}={\partial \over \partial s}\left({v^{2} \over 2}\right){\boldsymbol {e}}_{s}-{v^{2} \over R}{\boldsymbol {e}}_{r}} ただし、 s {\displaystyle s} は流線上のある点からの道のり、 r {\displaystyle r} は瞬間的な曲率中心からの距離、 R {\displaystyle R} は流線の曲率半径、 e s {\displaystyle {\boldsymbol {e}}_{s}} は接線方向の単位ベクトル、 e r {\displaystyle {\boldsymbol {e}}_{r}} は半径方向の単位ベクトルを表す。 加速度の流線方向の成分は流線にそった速さの変化率に対応し、加速度の法線方向の成分は流線が曲がることによる向心加速度に対応する。 外力のない非粘性のバロトロピック流体の定常な流れを考える。非粘性流体の流れを記述するオイラー方程式 D v D t = − 1 ρ ∇ p + f {\displaystyle {\mathrm {D} {\boldsymbol {v}} \over \mathrm {D} t}=-{1 \over \rho }\nabla p+{\boldsymbol {f}}} は定常、外力がない、バロトロピックという条件では ∂ ∂ s ( v 2 2 ) e s − v 2 R e r = − ∇ ∫ d p ρ {\displaystyle {\partial \over \partial s}\left({v^{2} \over 2}\right){\boldsymbol {e}}_{s}-{v^{2} \over R}{\boldsymbol {e}}_{r}=-\nabla \int {\mathrm {d} p \over \rho }} と変形できる。 方程式の両辺にそれぞれ e s , e r {\displaystyle {\boldsymbol {e}}_{s},\,{\boldsymbol {e}}_{r}} を内積でかけることで、流線方向(接線)成分、半径方向（主法線）成分は、 ∂ ∂ s ( v 2 2 ) = − ∂ ∂ s ∫ d p ρ ∴   ∂ ∂ s ( v 2 2 + ∫ d p ρ ) = 0 − v 2 R = − ∂ ∂ r ∫ d p ρ | r = R ∴   ∂ p ∂ r = ρ v 2 r ( o r   ∂ ∂ r ∫ d p ρ = v 2 r ) {\displaystyle {\begin{aligned}{\partial \over \partial s}\left({v^{2} \over 2}\right)&=-{\partial \over \partial s}\int {\mathrm {d} p \over \rho }\qquad \quad \therefore ~{\partial \over \partial s}\left({v^{2} \over 2}+\int {\mathrm {d} p \over \rho }\right)=0\\-{v^{2} \over R}&=-\left.{\partial \over \partial r}\int {\mathrm {d} p \over \rho }\right|_{r=R}\quad \therefore ~{\partial p \over \partial r}=\rho {v^{2} \over r}\quad \left(\mathrm {or~} {\partial \over \partial r}\int {\mathrm {d} p \over \rho }={v^{2} \over r}\right)\end{aligned}}} と表せる。ただし、方向微分の性質： e s ⋅ ∇ = ∂ ∂ s , e r ⋅ ∇ = ∂ ∂ r {\displaystyle {\boldsymbol {e}}_{s}\cdot \nabla ={\partial \over \partial s},\,{\boldsymbol {e}}_{r}\cdot \nabla ={\partial \over \partial r}} を使った。 第1式がベルヌーイの定理、第2式が流線曲率の定理に対応する。 移流項における ∇ φ {\displaystyle \nabla \varphi } は スカラー量の勾配であるが、対流項における ∇ A {\displaystyle \nabla {\boldsymbol {A}}} はベクトル量の共変微分である。ベクトル量の対流項 v ⋅ ∇ A {\displaystyle {\boldsymbol {v}}\cdot \nabla {\boldsymbol {A}}} を ( v ⋅ grad ) A {\displaystyle ({\boldsymbol {v}}\cdot \operatorname {grad} ){\boldsymbol {A}}} と記述することがあるが、この表示はデカルト座標系でしか等価でないことに注意すべきである（スカラー量の対流項 v ⋅ ∇ φ {\displaystyle {\boldsymbol {v}}\cdot \nabla \varphi } については ( v ⋅ grad ) φ {\displaystyle ({\boldsymbol {v}}\cdot \operatorname {grad} )\varphi } と等価である）。 共変微分を使わずに一般の座標系で成り立つ表現としては D A D t = ∂ A ∂ t + 1 2 { grad ⁡ ( v ⋅ A ) + rot ⁡ v × A + rot ⁡ A × v − rot ⁡ ( v × A ) + v div ⁡ A − A div ⁡ v } {\displaystyle {\frac {\mathrm {D} {\boldsymbol {A}}}{\mathrm {D} t}}={\frac {\partial {\boldsymbol {A}}}{\partial t}}+{\frac {1}{2}}\left\{\operatorname {grad} ({\boldsymbol {v}}\cdot {\boldsymbol {A}})+\operatorname {rot} {\boldsymbol {v}}\times {\boldsymbol {A}}+\operatorname {rot} {\boldsymbol {A}}\times {\boldsymbol {v}}-\operatorname {rot} ({\boldsymbol {v}}\times {\boldsymbol {A}})+{\boldsymbol {v}}\,\operatorname {div} {\boldsymbol {A}}-{\boldsymbol {A}}\,\operatorname {div} {\boldsymbol {v}}\right\}} がある。特に加速度の回転形表示 D v D t = ∂ v ∂ t + grad ⁡ ( | v | 2 2 ) − v × rot ⁡ v {\displaystyle {\frac {\mathrm {D} {\boldsymbol {v}}}{\mathrm {D} t}}={\frac {\partial {\boldsymbol {v}}}{\partial t}}+\operatorname {grad} \left({\frac {|{\boldsymbol {v}}|^{2}}{2}}\right)-{\boldsymbol {v}}\times \operatorname {rot} {\boldsymbol {v}}} は重要である。 曲線直交座標系 r = r ( q 1 , q 2 , q 3 ) {\displaystyle {\boldsymbol {r}}={\boldsymbol {r}}(q^{1},q^{2},q^{3})} における対流項 v ⋅ ∇ A {\displaystyle {\boldsymbol {v}}\cdot \nabla {\boldsymbol {A}}} の j {\displaystyle j} 成分は以下のように与えられる。 [ v ⋅ ∇ A ] j = ∑ k { v k h k ∂ A j ∂ q k + A k h k h j ( v j ∂ h j ∂ q k − v k ∂ h k ∂ q j ) } {\displaystyle [{\boldsymbol {v}}\cdot \nabla {\boldsymbol {A}}]_{j}=\sum _{k}\left\{{\frac {v_{k}}{h_{k}}}{\frac {\partial A_{j}}{\partial q^{k}}}+{\frac {A_{k}}{h_{k}h_{j}}}\left(v_{j}{\frac {\partial h_{j}}{\partial q^{k}}}-v_{k}{\frac {\partial h_{k}}{\partial q^{j}}}\right)\right\}} ただし、 h k = | ∂ r ∂ q k | = g k k {\displaystyle h_{k}=\left|{\partial {\boldsymbol {r}} \over \partial q^{k}}\right|={\sqrt {g_{kk}}}} （ g i j {\displaystyle g_{ij}} は計量テンソル）である。 先で述べたように [ ( v ⋅ grad ) A ] j = ∑ k v k ( 1 h k ∂ ∂ q k ) A j {\displaystyle [({\boldsymbol {v}}\cdot \operatorname {grad} ){\boldsymbol {A}}]_{j}=\sum _{k}v_{k}\left({\frac {1}{h_{k}}}{\frac {\partial }{\partial q^{k}}}\right)A_{j}} とはデカルト座標系 ( h 1 = h 2 = h 3 = 1 ) {\displaystyle (h_{1}=h_{2}=h_{3}=1)} においてのみ等しい。 A = v {\displaystyle {\boldsymbol {A}}={\boldsymbol {v}}} とした時の物質微分(＝加速度)の対流項に現れる第2項 ∑ k { v k h k h j ( v j ∂ h j ∂ q k − v k ∂ h k ∂ q j ) } {\displaystyle \sum _{k}\left\{{\frac {v_{k}}{h_{k}h_{j}}}\left(v_{j}{\frac {\partial h_{j}}{\partial q^{k}}}-v_{k}{\frac {\partial h_{k}}{\partial q^{j}}}\right)\right\}} は曲線直交座標系で現れる見かけの力に対応する。 実際、 L = m 2 ∑ k ( h k q ˙ k ) 2 {\displaystyle L={\frac {m}{2}}\sum _{k}(h_{k}{\dot {q}}^{k})^{2}} に対して ( d / d t ) ( ∂ L / ∂ q ˙ j ) − ( ∂ L / ∂ q j ) = 0 {\displaystyle ({\mathrm {d} /\mathrm {d} t})({\partial L/\partial {\dot {q}}^{j}})-({\partial L/\partial q^{j}})=0} を計算すると、 d v j d t − ∑ k { v k h k h j ( v j ∂ h j ∂ q k − v k ∂ h k ∂ q j ) } = 0 {\displaystyle {\mathrm {d} v_{j} \over \mathrm {d} t}-\sum _{k}\left\{{\frac {v_{k}}{h_{k}h_{j}}}\left(v_{j}{\partial h_{j} \over \partial q^{k}}-v_{k}{\partial h_{k} \over \partial q^{j}}\right)\right\}=0} が得られる。ただし、 v k = h k q ˙ k {\displaystyle v_{k}=h_{k}{\dot {q}}^{k}} であり、 h ˙ k = ∑ i ∂ h k ∂ q i q ˙ i {\displaystyle {\dot {h}}_{k}=\sum _{i}{\partial h_{k} \over \partial q^{i}}{\dot {q}}^{i}} を使う。 時間 t {\displaystyle t} の代わりに固有時間 τ {\displaystyle \tau } による物質微分を構成することもできる．具体的に，次式で与えられる[要出典]． D A μ ( x ) D τ = v ν ∇ ν A μ . {\displaystyle {\frac {\mathrm {D} A^{\mu }(x)}{\mathrm {D} \tau }}=v^{\nu }\nabla _{\nu }A^{\mu }.} v ν {\displaystyle v^{\nu }} は流体の四元速度， A μ ( x ) {\displaystyle A^{\mu }(x)} は四元時空を変数とするベクトルである．とくに， A μ ( x ) {\displaystyle A^{\mu }(x)} を速度とすると，測地線方程式となり，零になる． D v μ ( x ) D τ = v ν ∇ ν v μ = v ν ( ∂ ν v μ + Γ ν λ μ v λ ) = ∂ x ν ∂ τ ∂ v μ ∂ x ν + Γ ν λ μ v λ v ν = ∂ 2 x μ ∂ τ 2 + Γ ν λ μ ∂ x λ ∂ τ ∂ x λ ∂ τ = 0. {\displaystyle {\begin{aligned}{\frac {\mathrm {D} v^{\mu }(x)}{\mathrm {D} \tau }}&=v^{\nu }\nabla _{\nu }v^{\mu }=v^{\nu }\left(\partial _{\nu }v^{\mu }+\Gamma _{\nu \lambda }^{\mu }v^{\lambda }\right)\\&={\frac {\partial x^{\nu }}{\partial \tau }}{\frac {\partial v^{\mu }}{\partial x^{\nu }}}+\Gamma _{\nu \lambda }^{\mu }v^{\lambda }v^{\nu }\\&={\frac {\partial ^{2}x^{\mu }}{\partial \tau ^{2}}}+\Gamma _{\nu \lambda }^{\mu }{\frac {\partial x^{\lambda }}{\partial \tau }}{\frac {\partial x^{\lambda }}{\partial \tau }}=0.\\\end{aligned}}}