層流管内流れ（Hagen-Poiseuille） — 先端技術と研究動向

その通り。パイプ流の遷移は流体力学の未解決問題の一つだ。Hagen-Poiseuille 流は全てのReで線形安定（線形安定性解析で不安定モードが見つからない）だが、実際にはRe $\approx 2000\text{--}2300$ で乱流に遷移する。これは「亜臨界遷移」と呼ばれる。

有限振幅の擾乱が必要なんだ。線形安定性解析は微小擾乱しか扱わない。有限振幅の3D擾乱が加わると、非線形効果で「遷移的成長（transient growth）」を経て乱流に遷移する。

Re $\approx 2000$ 付近では「パフ（puff）」と呼ばれる局在化した乱流構造が現れる。

• パフ (Re $\approx 1700\text{--}2300$): 乱流の塊が上流に鋭い前面、下流になだらかな後面を持つ。流れ方向に移流されつつ、分裂や消滅を繰り返す
• スラグ (Re > $2300$ 程度): 乱流領域が上流・下流の両方向に拡大。最終的に管全体が乱流化

Avila et al. (2011, Science) はDNSと実験で、パフの分裂・消滅が統計的過程であることを示し、臨界Re $\approx 2040$ を決定した。これは統計力学の directed percolation と同じ普遍性クラスに属する。

スペクトル法が標準的だ。Willis & Kerswell のコード「openpipeflow」がオープンソースで公開されている。円柱座標のFourier（軸方向・方位角方向）+ Chebyshev/有限差分（径方向）の組み合わせだ。

有限体積法のCFDコード（OpenFOAM等）でもDNS級の計算は可能だが、スペクトル法に比べて同じ精度を達成するのに $5\text{--}10$ 倍のメッシュ点数が必要になる。

基本的には使える。マイクロ流路（$D \sim 10\text{--}100 \, \mu$m）では Re が低い（$Re \ll 1$ が典型的）ので、慣性項は完全に無視でき、Stokes 流れの世界だ。

ただし以下の効果が重要になる場合がある。

• スリップ流れ: Knudsen 数 $Kn = \lambda / D > 0.01$ で壁面スリップが顕著に。連続体仮定の限界
• 電気浸透流: 壁面の電荷による Debye 層の効果。Hagen-Poiseuille 流にプラグ流が重畳
• 非ニュートン効果: 血液などの非ニュートン流体では速度プロファイルが放物線から逸脱

そうだ。COMSOL の Microfluidics モジュールでは、電気浸透流、電気泳動、拡散-反応を含むマルチフィジックス解析が可能だ。