層流管内流れ（Hagen-Poiseuille） — トラブルシューティングガイド
より充実した内容を laminar-pipe-flow.html でご覧いただけます。
よくある問題と対策
こんな簡単な問題で失敗することあるんですか？
意外とあるんだ。Hagen-Poiseuille は「うまくいくはず」の問題だからこそ、合わないときに問題の切り分けがしやすい。
1. 速度プロファイルが放物線にならない
確認ポイント:
- パイプが十分長いか: $L > 0.06 Re \cdot D$ でなければ、出口断面で速度は完全発達していない。Re=100 なら $L > 6D$ が必要
- 入口条件: 一様流を入口に与えた場合、助走区間を経て放物線に発達する。入口断面で放物線を期待してはいけない
- メッシュの軸対称性: 非構造メッシュの場合、断面内の非対称性がプロファイルを歪ませることがある
2. 圧力降下が理論と合わない
圧力降下が理論値の $128\mu LQ / (\pi D^4)$ とずれるんですが。
確認ポイント:
| 原因 | 対策 |
|---|---|
| 助走区間の付加圧力降下 | 完全発達領域内の2点間で $\Delta p$ を計算 |
| 単位系の不整合 | $\mu$ と $\nu$ の単位を確認。$\mu = \rho \nu$ |
| 出口境界条件 | Pressure Outlet ($p = 0$) が出口にあるか確認 |
| 数値拡散 | 1次精度スキームでは圧力勾配の精度も低下 |
3. 収束が遅い / 収束しない
Hagen-Poiseuille は線形問題なので、適切に設定すれば数百反復以内に収束するはずだ。収束が遅い場合：
- 緩和係数が低すぎる: SIMPLE法の場合、圧力緩和 $\alpha_p = 0.3$、速度 $\alpha_U = 0.7$ が標準的。不必要に下げていないか確認
- メッシュの歪み: 非直交性が高いと収束が悪化。
checkMeshでnon-orthogonality を確認 - 境界条件の不整合: 入口で質量流入があるのに出口が壁面になっている等の致命的ミス
4. 壁面せん断応力がゼロ
wallShearStress がゼロと出力されるんですが。
壁面条件が slip になっていないか確認すること。noSlip（OpenFOAM では fixedValue uniform (0 0 0)）が正しい。
また、OpenFOAM の wallShearStress post-processing ユーティリティは、乱流モデルが有効でないと正しく計算されない場合がある。層流の場合は wallShearStress ではなく、wallGradU から手動で $\tau_w = \mu (\partial u / \partial n)_{wall}$ を計算するほうが確実だ。
5. 周期境界条件が機能しない
OpenFOAM で cyclic 境界を使う際の注意点：
- inlet/outlet パッチの面の法線ベクトルが反対向きで、面の形状が完全に一致している必要がある
createPatchユーティリティで周期パッチを作成した後、checkMeshでエラーがないことを確認meanVelocityForceの方向ベクトルが流れ方向と一致しているか確認
基本的な問題こそ丁寧に検証する姿勢が大事なんですね。
その通り。Hagen-Poiseuille 流が正しく解けないなら、それ以上複雑な問題を解く資格がないということだ。CFDの「健康診断」として定期的にベンチマークを回すことを勧める。
「Re数が正しいのに解が合わない」のよくある罠
層流管内流れの検証で「Re数は2000以下にしたのに速度プロファイルが放物線にならない」という相談がたまにあります。大抵の原因は「壁の粗さモデルが有効になっている」か「出口境界条件が逆流を許していない」のどちらかです。特に後者は盲点で、出口でゼロ勾配条件を使っても、数値的な揺らぎで流れが逆流すると収束が壊れます。もう一つの罠は「メッシュが粗すぎて放物線を十分解像できていない」ケース。断面方向に最低10要素は欲しいところを、5要素で「なんとなく形になってる」と思い込む。層流だから簡単、という先入観が落とし穴を招きます。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——層流管内流れ（Hagen-Poiseuille）の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告