共役熱伝達(CHT) — トラブルシューティング集
より充実した内容を conjugate-ht.html でご覧いただけます。
界面温度が不連続になる
先生、CHT解析で界面の固体側と流体側で温度が跳んでいるんですけど…
それは界面マッピングの設定ミスが最も疑わしい。まず界面がちゃんとcoupled wallとして認識されているか確認しよう。Fluentなら境界条件でwallタイプがcoupledになっているか、STAR-CCM+ならContact Interfaceが正しく生成されているかをチェック。
確認したけどちゃんとcoupledでした。他に原因は?
non-conformal interfaceで、マッピング先が見つからないfaceが存在する場合がある。Fluentのmesh interfaceでorphan faceが報告されていないか確認。STAR-CCM+ではContact設定のTolerance値が小さすぎて界面が部分的にしか接続されていないケースがある。
エネルギー残差が下がらない
運動量や連続の残差は下がるのに、energyの残差だけ高止まりします。
CHTでよくある症状だ。原因と対策を整理しよう。
| 原因 | 対策 |
|---|---|
| 固体-流体の熱容量比が極端に大きい | エネルギー方程式の緩和係数を0.85程度に下げる |
| 初期温度が実態と大きく乖離 | 初期温度を予想される分布に近づける |
| メッシュが界面近傍で粗すぎる | prism layerを追加、または壁面第一層を薄くする |
| 輻射モデル併用時の非線形性 | 輻射のunder-relaxationを0.5程度に下げる |
擬似非定常法(pseudo transient)を使えばいいって聞いたんですが?
FluentのPseudo Transient設定は定常CHTの収束改善にかなり効果的だ。Time Scale FactorをAutoにしておけば、ソルバーが自動的に適切な擬似時間刻みを選んでくれる。特に発熱量が大きい問題で有効だよ。
壁面熱流束が実験値と合わない
CFDで得た壁面熱流束が実験値の半分しかないんです。
まず乱流モデルの妥当性を疑おう。標準k-εは衝突噴流領域や再付着点付近で熱伝達を過小評価する傾向がある。SST k-ωモデルやRealizable k-εモデルに切り替えてみよう。$y^+$が十分小さい(1以下)ことも再確認すること。
乱流モデルを変えても改善しない場合は?
固体側の熱伝導率が正しいか、界面の接触熱抵抗が考慮されているか、流体の物性値(特にPrandtl数に効く粘度と熱伝導率)が正確かを順にチェックする。また、2D軸対称で解いている場合に実は3D効果が重要なケースもある。最終手段としてLES(Large Eddy Simulation)を検討する価値があるよ。
CHT解析で「温度が変になった」ときの一番多い原因
CHT解析のトラブル相談で圧倒的に多いのが「固体の温度が異様に高い(または低い)」という症状だ。原因の大半は熱伝導率の単位ミスだ。材料の熱伝導率をW/(m·K)で入力するべきところをW/(mm·K)で入力するだけで、値が1000倍違ってくる。CGS単位とSI単位が混在するような古いデータベースを流用したときに起きやすい。CHTで答えがおかしいと思ったらまず単位をチェック——経験豊富なエンジニアほど「あるある」と頷く基本の落とし穴です。
関連トピック
なった
詳しく
報告