ストークス流れ（低Reynolds数）
ストークス流れ（低Reynolds数）の理論基礎
概要
先生、ストークス流れって普通のNS方程式とどう違うんですか？
Reynolds数が非常に小さい（$Re \ll 1$）流れでは、慣性力が粘性力に比べて無視できる。NS方程式から移流項（非線形項）を落とした線形方程式がストークス方程式だ。微小粒子の沈降、マイクロ流体デバイス、生物流体力学などで重要になる。
ストークス方程式
定常・非圧縮ストークス方程式は次の通り。
運動量方程式:
連続の式:
NS方程式から $\rho(\partial\mathbf{u}/\partial t + \mathbf{u}\cdot\nabla\mathbf{u})$ を除いた形だ。圧力勾配と粘性力が完全にバランスしている。
どういう条件で移流項を無視できるんですか？
Re数で判定する。
| 対象 | 典型的なRe | ストークス近似 |
|---|---|---|
| 細菌の遊泳 | $10^{-5}$〜$10^{-4}$ | 非常に良好 |
| 精子の遊泳 | $10^{-2}$ | 良好 |
| 粒子沈降（$d=10\mu$m, 水中） | $10^{-3}$ | 良好 |
| 高粘度ポリマー流 | $10^{-2}$〜$10^{-1}$ | おおよそ有効 |
| マイクロ流路（$L=100\mu$m） | $0.1$〜$10$ | 場合による |
ストークスの抵抗法則
球（半径 $R$）が粘性流体中を一様速度 $U$ で運動するときの抵抗力は、ストークスが1851年に導いた。
これがストークスの抵抗法則だ。抵抗係数で書くと、
$Re = \rho U (2R)/\mu$ で定義。$Re < 1$ でこの関係が成り立つ。
ストークス流れの特性
ストークス流れには独特の数学的性質がある。
- 可逆性: 時間反転しても同じ流れ場（動画を逆再生しても区別できない）
- 瞬時応答: 慣性がないので、境界条件の変化に即座に応答
- 線形性: 解の重ね合わせが成立。複数の粒子の相互作用を解析的に扱える
- 一意性: 境界条件が与えられれば解は一意
可逆性って面白いですね。普通の流れではあり得ない。
Taylor-Couette装置で粘性流体中のインクの帯を「攪拌」してから逆回転すると元に戻る、というデモンストレーションが有名だ。これはストークス流の可逆性の直感的な実証だ。
細菌は「後ろ向きに泳いでも前に進む」世界に住んでいる
細菌（大腸菌など）が鞭毛を回転させて泳ぐ仕組みは、ストークス流れ（Re≪1）の世界でしか成立しません。「ホタテ貝の定理（スカロップ定理）」という有名な結果があり、「可逆的な動き（往復運動）ではRe≪1の流れで推進力はゼロになる」と証明されています。これが人間のような泳ぎ方を細菌ができない理由です。細菌は鞭毛を螺旋状に回転させることで非可逆な動きを実現し、この逆説を回避しています。マイクロロボットや医療ナノマシンの設計でも、この原理は重要な制約条件になっています。
ストークス流れ（低Reynolds数）の数値計算手法
ストークス方程式の数値解法
ストークス方程式を数値的に解く方法を教えてください。
ストークス方程式は線形なので、NS方程式より解きやすい。ただし速度-圧力の連成（鞍点問題）は残る。
鞍点問題としての定式化
ストークス方程式をFEMで離散化すると、以下の鞍点型の連立方程式になる。
ここで $A$ は粘性項の剛性マトリクス、$B$ は発散演算子の離散版。
ゼロ対角ブロックがあると解きにくそうですね。
その通り。inf-sup条件（LBB条件）を満たす要素の組み合わせが必要だ。
| 速度要素 | 圧力要素 | LBB安定 | 名称 |
|---|---|---|---|
| P2（二次三角形） | P1（一次三角形） | 安定 | Taylor-Hood要素 |
| P1+bubble | P1 | 安定 | MINI要素 |
| Q2（二次四角形） | Q1（一次四角形） | 安定 | 標準的 |
| P1 | P1 | 不安定 | 安定化が必要 |
| P1 | P0 | 不安定 | 使用不可 |
FVMでの解法
有限体積法ではSIMPLEアルゴリズムがそのまま使えるが、Re数が極めて小さいため特有の配慮が必要。
- 圧力の緩和係数: 通常より大きく（0.5〜0.8）設定可能。移流項がないため安定しやすい
- 移流スキーム: 不要（移流項がゼロだから）
- 収束: 線形問題なので通常数十反復で収束
境界積分法（BEM）
ストークス流れの線形性を活用した強力な手法が境界積分法だ。グリーン関数（ストークスレット）を使い、領域全体でなく境界面上だけに未知数を配置する。
ここで $G_{ij}$ はOseen-Burgers テンソル（自由空間のグリーン関数）。3D問題では体積メッシュが不要になり、計算量が大幅に削減される。
ストークス流れだと境界だけで解けるんですね。NS方程式ではこれは無理ですよね？
そう。NS方程式の非線形性のためにグリーン関数が存在しない。ストークス流れの線形性が可能にする特権的な手法だ。
ストークス流れの数値解法——境界要素法(BEM)が有限体積法より有利な理由
ストークス方程式（Re→0の低Re極限）は線形偏微分方程式であるため、境界要素法（BEM）が特に有効だ。BEMは流体領域内部をメッシュ化せず、境界面のみを離散化するため、3次元問題でも2次元のメッシュしか必要としない。特に外部流れ（開領域）では無限遠の境界を自然に扱え、有限体積法が人工的に設定しなければならない「遠方境界」が不要になる。ただしBEMの行列（N²の密行列）はメッシュ数Nに対して格納容量O(N²)、求解コストO(N³)が必要で、大規模問題には高速多重極法（FMM: Fast Multipole Method）を組み合わせてO(N log N)に圧縮する手法が使われる。流体中の粒子群・気泡群の多体問題にも有効な手法だ。
ストークス流れ（低Reynolds数）の実務適用
実践ガイド
ストークス流れの解析を実務で行うケースを教えてください。
主な応用分野を挙げよう。
粒子沈降速度の計算
球形粒子の終端沈降速度はストークスの抵抗法則から直接求まる。重力と抵抗力が釣り合う条件で、
| 粒子 | 直径 | 流体 | 沈降速度 | Re |
|---|---|---|---|---|
| 花粉 | 30 $\mu$m | 空気 | 2.7 cm/s | 0.05 |
| 砂粒 | 100 $\mu$m | 水 | 0.8 cm/s | 0.08 |
| 血球 | 8 $\mu$m | 血漿 | 0.4 $\mu$m/s | $3\times10^{-6}$ |
| 霧粒 | 10 $\mu$m | 空気 | 0.3 cm/s | 0.02 |
$Re > 1$ では修正が必要ですか？
その通り。Oseen補正やSchiller-Naumann式を使う。
CFDの粒子追跡（DPMモデル）ではこの補正付き抗力モデルがデフォルトで使われている。
マイクロ流体デバイスの設計
マイクロ流路（幅 10〜100 $\mu$m）では $Re$ が $O(1)$〜$O(10^{-2})$ であり、ストークス近似が有効。
| 設計パラメータ | 計算方法 |
|---|---|
| 圧力損失 | $\Delta p = \frac{12\mu L Q}{wh^3}$（矩形断面の近似） |
| 混合長 | $L_{\text{mix}} \sim Pe \cdot w$（Pe数 $= Uw/D$） |
| Dean流二次流れ | $De = Re\sqrt{d/R}$ で強度を評価 |
マイクロ流路だと混合が難しいって聞いたんですけど？
ストークス流の可逆性のためだ。乱流混合が起きないので拡散に頼るしかない。そこでジグザグ流路やヘリングボーン構造でカオス的移流を誘起する設計が使われる。
CFDでの設定のコツ
低Re数の解析特有の注意点を整理しよう。
| 項目 | 設定 | 理由 |
|---|---|---|
| 乱流モデル | Laminar（層流） | Re < 1 では乱流は発生しない |
| 移流スキーム | 中心差分でOK | Pe数が小さい |
| 密度の扱い | 非圧縮 | Ma数は極めて小さい |
| 浮力 | Boussinesq or 密度差 | 沈降問題では重力必須 |
| 収束 | 容易（数十反復） | 線形問題に近い |
低Re数の解析は収束が楽でいいですね。
そう。ただし、物理現象がスケール依存（表面張力、電気浸透流、ファンデルワールス力等）になるので、支配的な力を正しく見極める必要がある。
MEMS流路のストークス流れ——マイクロバルブ設計でのCFD適用事例
MEMSデバイスの微細流路（幅1〜100μm）ではReynolds数がRe<<1のストークス流れ（クリープ流れ）が支配的で、慣性力は無視できる。医療用薬液送達ポンプのマイクロバルブ設計では、弁体と弁座の隙間（サブμmオーダー）でのストークス流れ解析が重要だ。実際の設計事例では、スリット幅0.5μmのシリコンマイクロバルブのCFD解析（ストークスソルバ使用）で、弁前後の圧力差0.1kPaに対する流量係数Cdを予測し、製造公差±0.1μmが流量に与える感度を事前評価することで初期プロトタイプの試作回数を3回から1回に削減した事例が学術誌Microfluidics and Nanofluidicsで報告されている。
ストークス流れ（低Reynolds数）のソフトウェア比較
低Re数流れの商用ツール対応
ストークス流れを解くのに適したソルバーはどれですか？
汎用CFDソルバーはすべて層流モードでストークス流れを解ける。ただし、特殊な応用ではFEMやBEMが有利な場合がある。
ソルバー比較
| ソルバー | 手法 | 低Re流れの特徴 | 適用分野 |
|---|---|---|---|
| Ansys Fluent | FVM | Laminar設定で対応。DPMで粒子追跡 | 粒子沈降、フィルタ |
| STAR-CCM+ | FVM | Laminar + DEM連成が強力 | 粉体、スラリー |
| COMSOL | FEM | Creeping Flow物理モデルあり | マイクロ流体、MEMS |
| OpenFOAM | FVM | icoFoam（層流専用） | 学術、マイクロ流体 |
| Stokesflow (BEM) | BEM | 境界要素法。多粒子問題に最適 | 生物流体、コロイド |
COMSOLにCreeping Flowっていう専用モデルがあるんですね。
そう。COMSOLのCreeping Flow物理インターフェイスでは移流項が最初から除かれたストークス方程式を解く。これにより安定性が向上し、極端に低いRe数でも安定して計算できる。
粒子追跡モデルの比較
| 手法 | Fluent | STAR-CCM+ | COMSOL |
|---|---|---|---|
| DPM (Discrete Phase Model) | 対応 | 対応 | Particle Tracing |
| DEM (Discrete Element Method) | DEM連成あり | DEM標準搭載 | 制限あり |
| 粒子間相互作用 | 簡易モデル | 完全DEM対応 | 簡易 |
| ブラウン運動 | Brownian forceオプション | 対応 | 対応 |
| 電気泳動力 | UDFで実装 | Field force | Electrophoresis |
マイクロ流体デバイス設計に適したツール
マイクロ流体の設計ではCOMSOLが特に強い。マルチフィジックスの統合が容易で、以下のモジュールを組み合わせられる。
- Creeping Flow + Transport of Diluted Species（拡散混合）
- Creeping Flow + AC/DC Module（電気浸透流）
- Two-Phase Flow + Surface Tension（液滴操作）
マルチフィジックスが絡むとCOMSOLの強みが出るんですね。
逆に、単純な粒子沈降や管内流の圧損計算なら、FluentやOpenFOAMのLaminarモードで十分だ。ツール選定は解きたい物理の複雑さで決まる。
ストークス流れ対応CFDツール——COMSOL Multiphysicsのクリープ流れ解析機能
低Reynolds数（Stokes流れ）解析に特に優れたCFDツールとして、COMSOL Multiphysicsは多くの研究者の支持を得ている。COMSOLは「Creeping Flow（クリープ流れ）インターフェース」を持ち、慣性項を除いたストークス方程式を直接解く設定が容易だ。また生物流体力学（血液・粘液）や電気浸透流（Electroosmotic Flow）との多物理連成が直感的なGUIで設定でき、マイクロ流体研究での採用率が高い。一方OpenFOAMにはStokesFoamは存在せず、simpleFoamの緩和係数を極端に絞って低Re解を得る対処が必要。汎用CFD（Fluent）は収束は遅いものの低Re解も得られるが、収束に要するイテレーション数がStokes専用ソルバの10〜100倍になることが多い。
ストークス流れ（低Reynolds数）の先端研究
先端トピック
ストークス流れの研究って今でも進んでるんですか？
生物流体力学、マイクロロボティクス、コロイド科学など、低Re数の世界は非常に活発な研究分野だ。
マイクロスイマーの流体力学
ストークス流れの可逆性から、「帆立貝定理」（Scallop theorem, Purcell 1977）が導かれる。時間可逆な運動（1自由度の開閉のみ）では推進できない。
微生物は以下の戦略で推進する。
- 回転鞭毛（大腸菌）: ヘリカル波の伝播、非相反運動
- 繊毛打（ゾウリムシ）: 有効打と回復打の非対称性
- 波動運動（精子）: 正弦波の伝播
帆立貝みたいに開閉するだけでは泳げないんですね。
ストークス流の可逆性の直接的帰結だ。マイクロロボットの設計でも、非相反運動の実現が核心的課題になっている。
ストークス方程式の高度な解析解
球以外の形状に対するストークス抵抗の解析解も重要だ。
| 形状 | 抵抗力 | 備考 |
|---|---|---|
| 球 | $F = 6\pi\mu RU$ | Stokes (1851) |
| 楕円体（長軸方向） | $F = \frac{16\pi\mu aU}{2\ln(2a/b) - 1}$ | Oberbeck (1876), 細長体近似 |
| 円板（法線方向） | $F = 16\mu RU$ | Lamb |
| 円柱（無限長、2D） | 解なし（ストークスのパラドックス） | Oseen補正が必要 |
2D円柱で解がないんですか？
有名な「ストークスのパラドックス」だ。2次元では粘性流のストークス方程式が境界条件を満たす解を持たない。Oseen方程式（移流項の線形化版）で解消される。3次元球では問題なく解が存在する。
非ニュートン流体のストークス流
高粘度のポリマーや血液はニュートン流体ではない。せん断減粘性（shear-thinning）を示す。
Cross模型やCarreau模型がCFDでよく使われる。CFDソルバーでは材料モデルとして設定する。
アクティブマター（活物質）の集団運動
多数のマイクロスイマーの集団運動は「アクティブマター」として物理学の最先端テーマだ。各スイマーが生成するストークスレット的な流れ場の相互作用で、秩序的パターンや乱流的挙動が創発する。
低Reの世界は単純に見えて、実はものすごく豊かな物理があるんですね。
ストークスの時代から170年以上経った今も、新しい発見が続いている。Re数が小さいからといって物理が単純なわけではないのだ。
ストークス流れのフロンティア——マイクロスイマーの推進メカニズムとCFD
細菌や精子などの微生物（マイクロスイマー）は Re<<1 のストークス流れ域で自己推進する。この流れ域では「スカロップ定理（Scallop Theorem）」が成立し、往復対称な形状変化では前進できない——細菌が螺旋状鞭毛を回転させて進む理由はここにある。最先端の研究では、このストークス流れの物理的制約を利用したナノロボット（薬剤輸送用マイクロマシン）の推進機構設計にCFDが使われている。境界要素法（BEM: Boundary Element Method）は流体-構造連成が少ない計算コストで実現でき、ストークス流れ解析に最適な手法だ。MITやハーバード大でのマイクロスイマーCFD研究は医療応用への実用化を目指す活発なフィールドとなっている。
ストークス流れ（低Reynolds数）のトラブル対応
トラブルシューティング
低Re数の解析で起きやすいトラブルってありますか？
高Re数とは違う種類の問題が起きる。精度と物理モデルの問題が中心だ。
1. 圧力場が安定しない（FEMの場合）
症状: 圧力にチェッカーボード状の振動、または圧力がゼロに近い異常値。
原因: 速度-圧力の要素の組み合わせがinf-sup条件（LBB条件）を満たしていない。
対策:
- Taylor-Hood要素（P2/P1）を使用
- 安定化手法（PSPG, GLS）を適用
- COMSOLならCreeping Flow物理で自動的に安定な定式化が選ばれる
2. 粒子の抵抗力がストークス法則と合わない
CFDで球の抗力を計算したら、$6\pi\mu RU$ と10%以上ずれるんですが…
| 原因 | 対策 |
|---|---|
| 計算領域が小さい | 球径の50倍以上の外部領域を確保 |
| メッシュが粗い | 球表面に30要素以上（周方向） |
| 壁面の影響 | 無限遠条件（遠方境界で一様流） |
| Re > 0.1 | Oseen補正を考慮: $C_D = (24/Re)(1 + 3Re/16)$ |
ストークス解は無限領域を仮定しているので、有限の計算領域では壁面効果が不可避。外部境界を十分遠くに置くことが重要だ。
3. マイクロ流路の流量が理論値と合わない
症状: ハーゲン-ポアズイユの理論流量と大きくずれる。
確認事項:
- 流路断面が円管でない場合、断面形状に応じた修正係数が必要
- 入口の助走区間（$L_e \approx 0.06 Re D$、ストークス流ではほぼゼロ）
- 壁面のスリップ条件: マイクロスケールではKnudsen数 $Kn = \lambda/L$ が有限になり、no-slipが成り立たない場合がある
4. 非ニュートン流体での収束困難
症状: べき乗則流体で粘度が局所的に非常に大きく（または小さく）なり、収束しない。
対策:
- 粘度の上下限を設定（$\mu_{\min} < \mu_{\text{eff}} < \mu_{\max}$）
- Carreau模型を使用（ゼロせん断粘度と無限せん断粘度で上下を規定）
- 緩和係数を下げる
- 初期値としてニュートン流体（平均粘度）の解を使用
5. 多粒子問題での相互作用の過小評価
粒子がたくさんある場合、一つずつ計算するのはダメですか？
孤立粒子の仮定は粒子間距離が十分大きいときだけ有効。体積分率が0.1を超えると粒子間の流体力学的相互作用が無視できない。
対策:
ストークス流でも実務上のトラブルは色々あるんですね。
方程式が線形だからといって解析が簡単とは限らない。スケールの問題、壁面効果、非ニュートン性など、低Re特有の課題を理解しておくことが大切だ。
ストークス流れの計算が「合わない」—— Re数を疑え
ストークス流れ（Re≪1）の解析で最もよくある失敗は、「低速だからストークス近似でいいだろう」という思い込みです。Re数は流速だけでなく代表長さにも依存する。マイクロチャネル内の 1μm 粒子は Re≈0.001 で完全にストークス領域ですが、同じ流速でも 1mm の気泡になると Re≈1 に近づき慣性項が効いてくる。解析が実験値とずれ始めたら、まず問題の代表長さを見直す——これがストークス流れトラブルの鉄則です。
関連トピック
なった
詳しく
報告