粗密波の伝播速度と圧縮性の効果

問題の所在

流体力学において、運動量とインパルスの関係は一見すると奇妙な性質を示す。特に非圧縮性流体を無限領域で考えると、流体全体の運動量積分

$$
\int \B{u}\,\diff V
$$

は必ずしも通常の意味でよく定義されない。速度場が無限遠で十分速く減衰しないため、この積分は絶対収束せず、積分領域の取り方に依存することがある。

一方で、渦度分布から定義される流体力学的インパルス

$$
\B{I} = \frac{1}{2}\int \B{x}\times \boldsymbol{\omega}\,\diff V
$$

は、無限領域中の局在した渦運動に対して意味を持つ量である。非保存的外力が存在しない場合、このインパルスは保存量となる。流体力学的インパルスは、「静止状態からその流れを瞬間的に生成するために必要な全衝撃力」と等価に解釈できる。

しかし、ここで一つの疑問が生じる。非圧縮性流体では、圧力変化が無限速度で全空間に伝わると仮定される。そのため、ある場所に衝撃力を加えた瞬間、その影響が無限遠にまで及ぶように扱われる。この仮定が、運動量に関するいくつかの逆説的な結果を生んでいるのである。

この逆説が「非圧縮性」という理想化に由来することを示す点にある。わずかな圧縮性を導入すれば、圧力は音速で伝播し、運動量積分の発散的性質は解消される。

非圧縮性流体における逆説

非圧縮性流体では、圧力は拘束条件として働く。すなわち、速度場は常に

$$
\nabla\cdot \B{u}=0
$$

を満たさなければならず、そのために圧力場が瞬時に調整される。

この「瞬時に調整される」という点が重要である。現実の流体では圧力変化は音波として有限速度で伝わる。しかし非圧縮性流体では音速を無限大とみなすため、局所的な力の作用が全空間へ即座に伝達される。

その結果、無限領域における運動量の解釈が曖昧になる。流体運動量は流体力学的インパルスに表面積分項を加えた形で表され、この表面積分項は無限遠へ持っていっても消えるとは限らず、領域の形状に依存する。

つまり、非圧縮性流体においては、単純に「流体全体の運動量」を考えることが危うい。運動量そのものよりも、渦度分布に基づくインパルスの方が、無限領域の渦運動を記述するうえで自然な量となる。

わずかな圧縮性を導入する

圧力が無限の速度で伝播する問題を解決するため、わずかに圧縮性をもつ流体を考える。扱うのは、静止していた流体に原点で集中衝撃力を加えた場合である。

流体は静止状態から出発し、時刻 $t=0$ に集中した衝撃力を受ける。このとき、速度、圧力、密度の擾乱を小さいものとして線形化する。基礎方程式は次の形で与えられる。

$$
\left\{
\begin{split}
\, &\rho_0\frac{\partial \B{u}}{\partial t} = -\nabla p + \rho_0\,\B{I}\delta(\B{x})\delta(t) \EE
\, &\frac{\partial \rho}{\partial t} = -\rho_0\nabla\cdot\B{u} \EE
\, & p-p_0=c^2(\rho-\rho_0)
\end{split}
\right.
$$

ここで、$\rho_0$ は静止状態の密度、$p_0$ は静止状態の圧力、$c$ は音速、$\delta$ はデルタ関数である。また、$\rho_0\B{I}$ が加えられた流体力学的インパルスを表す。

この設定の意味は明快である。非圧縮性流体では無限遠まで瞬時に伝わっていた圧力応答を、圧縮性流体では有限速度 $c$ で伝わる波として扱うのである。

渦度は原点に集中する

運動方程式の回転を取ると、圧力項は消える。したがって、渦度は集中衝撃力の回転によって生成される。渦度は

$$
\boldsymbol{\omega} = \nabla \delta(\B{x})\times \B{I}\,H(t)
$$

の形で原点に集中する。ここで $H(t)$ はヘヴィサイド関数である。

これは、渦度そのものが空間全体に広がるわけではなく、原点に特異的に存在することを意味する。原点以外では渦度がゼロであるため、速度場はポテンシャル流れとして表すことができる。

$$
\B{u} = \nabla \phi \qquad (\B{x}\neq 0)
$$

すなわち、渦度は原点に集中し、その外側には音波として伝播する圧力・速度場が形成される。

圧力擾乱は波動方程式に従う

速度と密度を消去すると、圧力擾乱は波動方程式に従う。次の形の方程式が得られる。

$$
\frac{1}{c^2}\frac{\partial^2 p}{\partial t^2}- \nabla^2 p = -\B{I}\cdot \nabla \delta(\B{x})\delta(t)
$$

この式は、原点で加えられた衝撃力が圧力波を発生させることを表している。右辺は単なる点音源ではなく、$\B{I}$ の方向性をもつ双極子的な音源である。この波動方程式の解は、半径

$$
r=ct
$$

の球面上に特異的な波面をもつ。すなわち、時刻 $t$ において、衝撃の影響は半径 $ct$ の球の内側にだけ到達している。球の外側、すなわち

$$
r>ct
$$

では、圧力擾乱も速度場もまだ存在しない。これは非圧縮性流体との決定的な違いである。非圧縮性流体では力の効果が瞬時に無限遠まで及ぶが、圧縮性流体では効果は音速で有限に広がる。

球面内部の流れは双極子流れになる

波面の内側、すなわち

$$
r<ct
$$

では、速度場は双極子によるポテンシャル流れと同じ形を持つ。速度ポテンシャルは

$$
\phi = -\frac{1}{4\pi\rho_0}\B{I}\cdot\nabla \frac{H(t-r/c)}{r}
$$

の形で表される。球面内部では $H(t-r/c)=1$ であるため、速度場は実質的に

$$
\phi \sim -\frac{1}{4\pi\rho_0} \B{I}\cdot\nabla \left( \frac{1}{r} \right)
$$

という双極子ポテンシャルになる。

これは、非圧縮性流体でインパルスをもつ局在渦が遠方に作る双極子的な速度場と対応している。つまり、圧縮性を導入しても、波面の内側では非圧縮性理論に現れる双極子場が再現される。ただし、その場は最初から全空間に存在するのではなく、音速で拡大する球の内側に限って存在するのである。

「失われた運動量」はどこへ行くのか？

ここで最も重要な点が現れる。

非圧縮性理論では、インパルスと実際の運動量の間にずれが生じる。たとえば、無限大の球領域で流体運動量を評価すると、流体力学的インパルスそのものとは一致せず、表面項の寄与が現れる。これは「運動量がどこかへ消えた」ように見える。

しかし、圧縮性を導入すると、この解釈は明確になる。

球面 $r=ct$ の内側には双極子流れが存在し、その運動量は全インパルスの一部に対応する。一方、残りの「失われた」運動量は半径 (ct) の球面状の波面によって運ばれる。

すなわち、運動量は消えたのではない。音波として外側へ伝播する球面波の前線が担っているのである。

これは、非圧縮性極限で見えていた逆説の正体を示している。非圧縮性流体では音速 $c\to\infty$ とみなすため、球面波の前線は瞬時に無限遠へ飛ぶ。そのため、運動量の一部が無限遠の表面項に隠れているように見える。圧縮性流体として扱えば、その運動量は有限速度で伝播する波面上に存在することが明らかになる。

非圧縮性近似の意味

この議論は、非圧縮性流体力学の限界を示している。

非圧縮性近似は、低マッハ数の流れや、音波の伝播を問題にしない流れでは極めて有効である。しかし、運動量が空間的にどのように伝達されるかを厳密に問う場合、非圧縮性近似は物理的な伝播過程を消してしまう。

非圧縮性理論では、圧力は楕円型方程式によって決まり、境界条件の影響が瞬時に全領域へ及ぶ。一方、圧縮性理論では、圧力擾乱は波動方程式に従い、音速で伝わる。この違いが、運動量の解釈に直接反映される。

したがって、流体力学的インパルスを理解するうえでは、非圧縮性理論だけでなく、その背後にある圧縮性極限を意識することが重要である。インパルスは非圧縮性流体における便利な数学的量であるだけでなく、圧縮性流体における運動量伝達の極限としても解釈できる。

まとめ

今回の要点は、次のように整理できる。

第一に、非圧縮性流体における運動量の逆説は、圧力が無限速度で伝わるという仮定に由来する。

第二に、わずかな圧縮性を導入すると、圧力擾乱は音速 $c$ で伝わる球面波として記述される。

第三に、衝撃力によって生成された流れは、球面 $r=ct$ の内側では双極子流れとして振る舞い、外側ではまだ静止状態のままである。

第四に、非圧縮性理論で「失われた」ように見える運動量は、実際には球面状の音波前線によって運ばれている。

第五に、非圧縮性極限とは、この音波前線が無限遠へ瞬時に移動する極限である。そのため、運動量の一部が無限遠の表面項として現れる。

このように、圧縮性の効果は単なる補正ではない。流体インパルスと運動量の関係を物理的に理解するための鍵である。非圧縮性流体力学に現れる数学的な不思議さは、圧縮性流体の有限伝播速度を考慮することで、より自然な姿を取り戻すのである。