ストークス流れ（低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の流れで推進力はゼロになる」と証明されています。これが人間のような泳ぎ方を細菌ができない理由です。細菌は鞭毛を螺旋状に回転させることで非可逆な動きを実現し、この逆説を回避しています。マイクロロボットや医療ナノマシンの設計でも、この原理は重要な制約条件になっています。
各項の物理的意味
- 時間項 $\partial(\rho\phi)/\partial t$：蛇口をひねった瞬間を思い浮かべてください。最初は水がバタバタと不安定に出て、しばらくすると安定した流れになりますよね？ この「変化している最中」を記述するのが時間項です。心臓の拍動で血流が脈打つのも、エンジンのバルブが開閉するたびに流れが変動するのも、すべて非定常現象。では定常解析とは？ 「十分時間が経って流れが落ち着いた後」だけを見る——つまりこの項をゼロにする。計算コストが大幅に下がるため、まず定常で解いてみるのがCFDの基本戦略です。
- 対流項 $\nabla \cdot (\rho \mathbf{u} \phi)$：川に落ち葉を落としたらどうなりますか？ 流れに乗って下流に運ばれますよね。これが「対流」——流体の動きが物を運ぶ効果です。暖房の温風が部屋の端まで届くのも、空気という「運び屋」が熱を対流で輸送しているから。ここが面白いところ——この項は「速度×速度」を含むため非線形です。つまり、流れが速くなるとこの項が急激に強くなり、制御が難しくなる。これが乱流の根本原因です。よくある勘違い：「対流と伝導は同じようなもの」→ 全然違います！ 対流は流れが運ぶ、伝導は分子が伝える。桁違いの効率差があります。
- 拡散項 $\nabla \cdot (\Gamma \nabla \phi)$：コーヒーにミルクを入れて放置したことはありますか？ かき混ぜなくても、しばらく経つと自然に混ざりますよね。あれが分子拡散です。では次の質問——ハチミツとお水、どちらが流しやすいですか？ 当然お水ですよね。ハチミツは粘性（$\mu$）が高いから流れにくい。粘性が大きいと拡散項が強くなり、流体は「もったりした」動きになります。レイノルズ数が小さい流れ（ゆっくり、ドロドロ）では拡散が支配的。逆にRe数が大きい流れでは対流が圧倒し、拡散は脇役になります。
- 圧力項 $-\nabla p$：注射器のピストンを押すと、液体が針先から勢いよく出ますよね？ なぜでしょう？ ピストン側が高圧、針先が低圧——この圧力差が流体を押す力になるからです。ダムの放水も同じ原理。天気図で等圧線がギュッと密になっている場所では？ そう、強風が吹きます。「圧力差があるところに流れが生まれる」——これがナビエ-ストークス方程式の圧力項の物理的意味。ここでの勘違いポイント：CFDの「圧力」は絶対圧ではなくゲージ圧のことが多い。圧縮性解析に切り替えたとたんに結果がおかしくなる場合、絶対圧/ゲージ圧の混同が原因かもしれません。
- ソース項 $S_\phi$：暖められた空気が上に昇る——なぜでしょう？ 周囲より軽く（密度が低く）なったから、浮力で押し上げられるのです。この浮力はソース項として方程式に追加されます。他にも、ガスコンロの炎で化学反応熱が発生する、工場の電磁ポンプで金属溶湯にローレンツ力がかかる…これらはすべて「外部から流体にエネルギーや力を注入する」作用であり、ソース項で表現します。ソース項を忘れるとどうなるか？ 自然対流の解析で浮力を入れ忘れると、流体は一切動かない——冬の部屋で暖房をつけたのに暖かい空気が上に行かない、という物理的にありえない結果になります。
仮定条件と適用限界
- 連続体仮定：クヌッセン数 Kn < 0.01（分子平均自由行程 ≪ 代表長さ）で成立
- ニュートン流体仮定：せん断応力と歪み速度が線形関係（非ニュートン流体では粘度モデルが必要）
- 非圧縮性仮定（Ma < 0.3の場合）：密度を一定として扱う。マッハ数0.3以上では圧縮性効果を考慮
- ブシネスク近似（自然対流）：密度変化を浮力項のみで考慮し、他の項では一定密度を使用
- 適用外ケース：希薄気体（Kn > 0.1）、超音速・極超音速流れ（衝撃波捕捉が必要）、自由表面流れ（VOF/Level Set等が必要）
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 速度 $u$ | m/s | 入口条件で体積流量から換算する際、断面積の単位に注意 |
| 圧力 $p$ | Pa | ゲージ圧と絶対圧の区別。圧縮性解析では絶対圧を使用 |
| 密度 $\rho$ | kg/m³ | 空気: 約1.225 kg/m³@20°C、水: 約998 kg/m³@20°C |
| 粘性係数 $\mu$ | Pa·s | 動粘性係数 $\nu = \mu/\rho$ [m²/s] との混同に注意 |
| レイノルズ数 $Re$ | 無次元 | $Re = \rho u L / \mu$。層流/乱流遷移の判定指標 |
| CFL数 | 無次元 | $CFL = u \Delta t / \Delta x$。時間刻みの安定性に直結 |
数値解法と実装
ストークス方程式の数値解法
ストークス方程式を数値的に解く方法を教えてください。
ストークス方程式は線形なので、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）
ストークス流れの線形性を活用した強力な手法が境界積分法だ。グリーン関数（ストークスレット）を使い、領域全体でなく境界面上だけに未知数を配置する。
$$ u_j(\mathbf{x}) = -\frac{1}{8\pi\mu}\oint_S G_{ij}(\mathbf{x}, \mathbf{y})f_i(\mathbf{y})\,dS(\mathbf{y}) + \frac{1}{8\pi}\oint_S T_{ijk}(\mathbf{x}, \mathbf{y})u_i(\mathbf{y})n_k(\mathbf{y})\,dS(\mathbf{y}) $$
ストークス流れの線形性を活用した強力な手法が境界積分法だ。グリーン関数（ストークスレット）を使い、領域全体でなく境界面上だけに未知数を配置する。
ここで $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)に圧縮する手法が使われる。流体中の粒子群・気泡群の多体問題にも有効な手法だ。
風上差分（Upwind）
1次風上: 数値拡散が大きいが安定。2次風上: 精度向上するが振動のリスク。高レイノルズ数流れでは必須。
中心差分（Central Differencing）
2次精度だが、Pe数 > 2で数値振動が発生。低レイノルズ数の拡散支配流れに適する。
TVDスキーム（MUSCL、QUICK等）
リミッタ関数により数値振動を抑制しつつ高精度を維持。衝撃波や急勾配の捕捉に有効。
有限体積法 vs 有限要素法
FVM: 保存則を自然に満足。CFDの主流。FEM: 複雑形状・マルチフィジックスに有利。SPH等のメッシュフリー法も発展中。
CFL条件（クーラン数）
陽解法: CFL ≤ 1が安定条件。陰解法: CFL > 1でも安定だが、精度と反復回数に影響。LES: CFL ≈ 1を推奨。物理的意味: 1タイムステップで情報が1セル以上進まないこと。
残差モニタリング
連続の式・運動量・エネルギーの各残差が3〜4桁低下で収束と判断。質量保存の残差は特に重要。
緩和係数
圧力: 0.2〜0.3、速度: 0.5〜0.7が一般的な初期値。発散する場合は緩和係数を下げる。収束後は上げて加速。
非定常計算の内部反復
各タイムステップ内で定常解に収束するまで反復。内部反復数: 5〜20回が目安。残差がタイムステップ間で変動する場合は時間刻みを見直す。
SIMPLE法のたとえ
SIMPLE法は「交互に調整する」手法。まず速度を仮に求め（予測ステップ）、その速度で質量保存が満たされるよう圧力を補正し（補正ステップ）、補正された圧力で速度を修正する——このキャッチボールを繰り返して正解に近づく。2人で棚を水平にする作業に似ている：片方が高さを合わせ、もう片方がバランスを取り、これを交互に繰り返す。
風上差分のたとえ
風上差分は「川の流れに立って上流の情報を重視する」手法。川の中にいる人が下流を見ても水の出所は分からない——上流の情報が下流を決めるという物理を反映した離散化手法。精度は1次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
実践ガイド
ストークス流れの解析を実務で行うケースを教えてください。
主な応用分野を挙げよう。
粒子沈降速度の計算
球形粒子の終端沈降速度はストークスの抵抗法則から直接求まる。重力と抵抗力が釣り合う条件で、
| 粒子 | 直径 | 流体 | 沈降速度 | 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で報告されている。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
低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倍になることが多い。
選定で最も重要な3つの問い
- 「何を解くか」：ストークス流れ（低Reynolds数）に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端トピック
ストークス流れの研究って今でも進んでるんですか？
生物流体力学、マイクロロボティクス、コロイド科学など、低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研究は医療応用への実用化を目指す活発なフィールドとなっている。
トラブルシューティング
トラブルシューティング
低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 に近づき慣性項が効いてくる。解析が実験値とずれ始めたら、まず問題の代表長さを見直す——これがストークス流れトラブルの鉄則です。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——ストークス流れ（低Reynolds数）の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告