ストークス流れ（低Reynolds数） — トラブルシューティングガイド
この記事は統合版に移行しました
より充実した内容を stokes-flow.html でご覧いただけます。
より充実した内容を stokes-flow.html でご覧いただけます。
ストークス流れ（低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特有の課題を理解しておくことが大切だ。
Coffee Break よもやま話
ストークス流れの計算が「合わない」—— Re数を疑え
ストークス流れ（Re≪1）の解析で最もよくある失敗は、「低速だからストークス近似でいいだろう」という思い込みです。Re数は流速だけでなく代表長さにも依存する。マイクロチャネル内の 1μm 粒子は Re≈0.001 で完全にストークス領域ですが、同じ流速でも 1mm の気泡になると Re≈1 に近づき慣性項が効いてくる。解析が実験値とずれ始めたら、まず問題の代表長さを見直す——これがストークス流れトラブルの鉄則です。
関連トピック
この記事の評価
ご回答ありがとうございます！
参考に
なった
なった
もっと
詳しく
詳しく
誤りを
報告
報告