流体力学：ナビエ–ストークス方程式、レイノルズ数、乱流

流体力学は、液体や気体がどのように動くかを説明する分野です。入門レベルの問題では、主に3つの考え方が中心になります。ナビエ–ストークス方程式は力が運動をどう変えるかを示し、レイノルズ数はどの効果が重要かを判断する助けになり、乱流は擾乱が成長したときに現れる不規則で混合の強い流れを表します。

流体力学とは何か

流体はせん断を受けても形を一定に保てないため、流れながら変形し続けます。流体力学では、速度、圧力、密度、場合によっては温度といった量が、空間や時間とともにどう変化するかを追います。

この同じ枠組みは、管内流れ、血流、天気、空気力学、海流などに現れます。細部は違っても、繰り返し出てくる問いは同じです。何が流れを生み出しているのか、どの力が支配的なのか、そして流れは滑らかか不安定か、という点です。

ナビエ–ストークス方程式が教えてくれること

密度と粘性が一定の非圧縮ニュートン流体では、よく使われる形の1つは次のとおりです。

$$\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}$$

これに加えて、非圧縮条件

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

を満たします。

ここで $\mathbf{u}$ は速度、$p$ は圧力、$\rho$ は密度、$\mu$ は粘性係数、$\mathbf{f}$ は重力のような単位質量あたりの体積力です。

記号は少し重たく見えますが、考え方は単純です。左辺は、動いている流体要素がどのように加速するかを表します。右辺は、その加速が圧力差、隣り合う層による粘性抵抗、外力によって生じることを示しています。

ただし、この形が常に使えるわけではありません。流体が圧縮性をもつ場合、非ニュートン流体である場合、あるいは温度によって性質が大きく変わる場合には、モデルも変える必要があります。

レイノルズ数は最初の目安になる

レイノルズ数は、慣性の効果と粘性の効果を比べる無次元量です。

$$Re = \frac{\rho U L}{\mu} = \frac{U L}{\nu}$$

ここで $U$ は代表速度、$L$ は代表長さ、$\nu$ は動粘性係数です。

$Re$ が小さいと、粘性のなめらかにする効果が強く、流れは層流になりやすくなります。$Re$ が大きいと、慣性の影響が強くなり、擾乱は減衰するよりも成長しやすくなります。

ただし、これはあくまで目安であって、万能な切り替え値ではありません。滑らかな円管内流れでは、一般に $Re \lesssim 2300$ で層流的なふるまいが見られますが、遷移は形状、表面粗さ、流入時点ですでにどれだけ乱れているかにも左右されます。

計算例：管の中を流れる水

室温付近の水が、直径 $D = 0.02\ \mathrm{m}$ の管を平均速度 $U = 1.0\ \mathrm{m/s}$ で流れているとします。動粘性係数を

$$\nu \approx 1.0 \times 10^{-6}\ \mathrm{m^2/s}$$

とします。

$L = D$ とすると、レイノルズ数は

$$Re = \frac{U D}{\nu} = \frac{(1.0)(0.02)}{1.0 \times 10^{-6}} = 2.0 \times 10^4$$

となります。

したがって $Re \approx 20{,}000$ です。滑らかな円管内の内部流れとしては、これは通常の層流域を大きく超えています。したがって、層流モデルよりも乱流モデルから考え始めるほうがずっと安全です。

これがレイノルズ数の大きな価値です。速度場全体を与えてくれるわけではありませんが、単純な層流の見方が早い段階で破綻しそうかどうかを教えてくれます。

乱流とは何か

乱流は、単に「ぐちゃぐちゃな流れ」という意味ではありません。強く不規則な速度変動があり、多くの長さスケールにわたって混合が起こる運動です。

多くの乱流では、エネルギーは大きなスケールで流れに入り、より小さなスケールへと受け渡され、最終的には粘性によって熱として散逸します。この多重スケール構造こそが、乱流を細部まで計算するのが難しい大きな理由の1つです。

運動そのものは依然としてナビエ–ストークス方程式に従いますが、現実的な乱流に対して厳密な解析解が得られることはまれです。実際には、工学では実験、数値シミュレーション、簡略化したモデルに頼ることが多くなります。

流体力学でよくある間違い

レイノルズ数を魔法のしきい値だと思う

$Re$ は流れを分類する助けになりますが、乱流への遷移があらゆる状況で1つの値で決まるわけではありません。

代表長さを間違える

$Re$ の値は $L$ に依存します。管内流れでは通常 $L$ は直径ですが、球のまわりの流れや平板上の流れでは、別の長さスケールを使うのが適切です。

ナビエ–ストークス方程式の形を誤って使う

上に書いた形は、密度と粘性が一定の非圧縮ニュートン流体を仮定しています。圧縮性流れ、非ニュートン流体、温度依存性の強い性質をもつ流体では、別のモデル化が必要です。

粘性と密度を混同する

密度は、ある体積の中にどれだけ質量が詰まっているかを表します。粘性は、変形やせん断に対する抵抗の大きさを表します。どちらも流体力学に現れますが、役割は異なります。

乱流には構造がないと思い込む

乱流は不規則に見えますが、それでも渦、境界層、大規模なコヒーレント構造のような組織的な特徴をもっています。

流体力学はどこで使われるか

流体力学は、動く液体や気体が関わるあらゆる場面で使われます。たとえば、航空機や自動車の空気力学、ポンプやパイプライン、天気予報、心血管内の流れ、化学反応器、環境中の輸送現象などです。

方程式全体が複雑でも、実際の問いは同じです。何が流れを生み出しているのか。どの力が最も重要なのか。流れは滑らかに保たれそうか、それとも遷移や混合を見込むべきか、ということです。

似た問題に挑戦してみよう

ストローの中の流れ、シャワーの配管の流れ、あるいは車の窓から手を出したときの手のまわりを流れる空気について、レイノルズ数を見積もってみましょう。速度、長さスケール、流体の種類を変えるだけで、整った流れのままでいる場合と乱流になる場合がある理由がよくわかります。