クエット流れ — トラブルシューティングガイド
トラブルシューティング
クエット流れの計算でハマりやすいポイントは何ですか？
シンプルな問題だからこそ、間違いが明確に見える。
線形分布にならない
原因1: 圧力勾配が残っている。入出口BCの設定ミスで意図せぬ圧力差が生じている
- 対策: 周期境界条件を使うか、出口で $p = 0$ を明示的に固定
原因2: 移動壁の速度方向が間違い。壁面法線方向に速度成分があるとno-penetration条件と矛盾
- 対策: 移動壁速度は壁面接線方向のみに設定
Taylor-Couette渦が出ない
Ta > 1708 にしたのにTaylor渦が現れません。
原因1: 初期条件が完全な軸対称Couette解のまま→微小擾乱を加えていない
- 対策: 初期速度場にランダム擾乱（振幅0.1%程度）を重畳
原因2: 軸方向ドメインが渦波長と整合していない
- 理論波長: $\lambda \approx 2d$（$d = R_o - R_i$: ギャップ幅）
- 対策: $L_z = n \times 2d$（$n$: 整数）に設定
原因3: 2D軸対称計算にしている。Taylor渦は3D構造なので、3D計算が必要
移動壁面での圧力振動
非構造格子でcollocated配置（圧力・速度が同一セル）の場合、移動壁面近傍で圧力のチェッカーボード振動が起こることがある。
- 対策: Rhie-Chow補間が正しく機能しているか確認。壁面近傍のメッシュを均一にする
粘度計シミュレーションの注意点
回転粘度計のシミュレーションで気をつけることは？
- 端面効果: 実際の粘度計は有限長なので、上下端面の影響がある。Guard ring がない場合は端面を含めた3D計算が必要
- 非ニュートン流体: せん断速度がギャップ内で一様でないため、粘度のせん断速度依存性が結果に影響する。ギャップ比 $R_o/R_i < 1.1$ が広ギャップ補正不要の目安
粘度計の結果とシミュレーションを比較するには端面効果を考慮しないといけないんですね。
F1と空力の戦い
F1マシンは時速300kmで走ると、車重と同じくらいのダウンフォース（下向きの空力的な力）を発生します。つまり理論上、天井に貼り付けて走れる！ チームは数千CPU時間のCFDシミュレーションを毎週実行し、フロントウィングの角度を0.1°単位で最適化しています。F1はCAEの技術力がそのまま順位に直結する世界です。
トラブル解決の考え方
デバッグのイメージ
CFDのデバッグは「水道管の詰まり修理」に似ている。まず「どこで詰まっているか」（どの残差が下がらないか）を特定し、次に「何が詰まっているか」（メッシュ品質？境界条件？乱流モデル？）を調べ、最後に「どう直すか」（メッシュ修正？緩和係数？）を判断する。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——クエット流れの問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
CFDメッシュの品質管理や乱流モデルの選定に悩む時間を、もっと創造的な設計作業に使えたら。 — Project NovaSolverはそんな実務者の声から生まれました。
CAEの未来を、実務者と共に考える
Project NovaSolverは、クエット流れにおける実務課題の本質に向き合い、エンジニアリングの現場を支える道具づくりを目指す研究開発プロジェクトです。
プロジェクトの最新情報を見る →