ANSYS CFX — CAE用語解説
ANSYS CFX
ANSYS CFXとANSYS Fluentって同じANSYS社のCFDソフトですよね。どっちを使えばいいんですか?
理論と物理
基本概念と支配方程式
CFXで「定常」と「非定常」を選ぶ基準は何ですか?時間変化があるかどうかだけでは判断が難しい気がします。
良い質問だ。定性的な判断だけでは不十分で、無次元数が基準になる。例えば、ストローハル数(Strouhal number)
RANS方程式で出てくる「渦粘性係数」って、具体的にどう決まるんですか?物性値のように固定なのでしょうか。
いや、流れ場の状態によって空間的に大きく変動する計算量だ。最も基本的な零方程式モデルである「混合距離モデル」では、壁からの距離など幾何学的な量で決める。一方、CFXでよく使われる
「圧力基盤」と「密度基盤」という言葉を聞きます。CFXはどちらですか?なぜそれが重要なのですか。
CFXは「圧力基盤」ソルバーだ。これは、連続の式と運動量方程式をカップリングさせ、圧力補正方程式(ポアソン型の方程式)を解く手法を中核としている。一方、Fluentのデフォルトは「密度基盤」で、質量、運動量、エネルギーの保存式を直接的に解く。圧力基盤は、非圧縮性から圧縮性流れまで一貫したアルゴリズムで扱え、特に低速流れや複雑な物理モデル(多相流、化学反応)のカップリングに強い。マッハ数0.3以下の自動車の外気解析やポンプ内部流れでは、CFXのこの特性が有利に働くことが多い。
「収束判定」の残差を見ていますが、10^-4とか10^-5まで下がればOKという感覚でいます。これに物理的な意味はあるんですか?
残差の値そのものに絶対的な物理的意味はない。離散化された方程式の「不平衡さ」の尺度だ。10^-4を一つの目安にするのは、多くの工業的問題でこれ以上下げても、抗力係数や効率などの注目する物理量の変化が1%以内に収まる経験則があるからだ。しかし、重要なのは「モニタポイント」の値が定常状態に達しているかどうかだ。例えば翼型解析なら、揚力係数Clの時系列をプロットし、その変動が±0.001以内で安定することを確認する。残差が低くても、モニタポイントが振動している場合は解が発振している可能性がある。
数値解法と実装
離散化とソルバー設定
CFXの「高分解能」と「二風上差分」のスキーム選択で迷います。マニュアルでは高分解能を推奨とありますが、二風上差分は使わない方がいいのですか?
「高分解能」は、本質的には二風上差分をベースにした「ブレンドファクター」を各要素・各変数ごとに自動計算し、精度を最大化する手法だ。一方、純粋な「二風上差分」は1次の精度で数値拡散が大きく、特にせん断層や渦の解像度が悪い。ではなぜオプションにあるか? 計算が非常に安定で、荒いメッシュでもとりあえず流れ場の傾向を得たい初期段階や、収束が難しい複雑な問題の「とっかかり」として使うことがある。ただし最終結果には使うべきではない。高分解能は実質的に2次精度に近づけ、工業的精度を満たす標準的な選択肢だ。
「タイムステップ」の大きさはどう決めればいいですか?小さすぎると計算時間が、大きすぎると発散しそうで心配です。
非定常解析のキーパラメータだ。物理的な現象を捉えるために必要な「時間分解能」と、数値的に安定な「CFL条件」の両方を考慮する。例えば、羽根車が1回転する時間が0.01秒で、1回転を20ステップで捉えたいなら、タイムステップは0.0005秒だ。一方、CFL数は
「多重グリッド」オプションがありますが、これは何をしているんですか?ON/OFFで計算時間が大きく変わります。
ソルバーの収束を劇的に加速する手法だ。通常の計算(細かいグリッド)では、長波長の誤差成分がなかなか減衰しない。多重グリッドは、元の細かいメッシュ(Fine Grid)と、仮想的に作った粗いメッシュ(Coarse Grid)を何段階も用意する。粗いメッシュでは長波長の誤差を素早く減衰させ、その情報を細かいメッシュにフィードバックする。CFXでは圧力補正方程式のソルバーに対して適用される。デフォルトでONだが、極めて複雑な幾何形状や、不安定な流れの初期計算では、一度OFFにして安定化させ、その後ONにするというテクニックも使われる。計算時間は30〜70%短縮できることが多い。
実践ガイド
ワークフローとチェックリスト
解析を始める前の「前処理チェックリスト」のようなものはありますか?いきなり計算を回して失敗するのが多いです。
必ず以下の5項目を確認せよ。1. **単位系の一貫性**: CFXはkg, m, s, K, JのSI単位系が内部標準。入力値がこれと一致しているか。圧力がPaか、粘度がPa・sか。2. **ドメインのタイプ**: 「流体」ドメインのみか、「固体」ドメインが混在する共役伝熱(Conjugate Heat Transfer)か。3. **境界条件の物理的妥当性**: 流入速度と流入温度は同時に指定できるが、圧力出口と温度は同時には指定できない(どちらかは計算される)。4. **インターフェース**: 回転/静止域の「フローズンロータ」インターフェース設定は正しいか。5. **初期条件**: 特に圧力場は「ゼロ」ではなく、出口静圧など現実的な値で初期化する。
メッシュの品質指標で「Orthogonal Angle」と「Aspect Ratio」を見ますが、許容される基準値はありますか?
CFXのメッシュ適合格判定では明確な基準がある。**Orthogonal Angle(直交角)**: 理想的には90度。10度を下回る要素が多数あると、拡散項の計算精度が著しく低下する。最低でも20度以上をキープしたい。**Aspect Ratio(アスペクト比)**: 要素の長辺と短辺の比。一般的な流れ場では100以下が望ましい。しかし、境界層メッシュのように流れ方向に極端に細長い要素では、アスペクト比が1000を超えることもある。その場合は、壁に垂直な方向の要素サイズが適切か(y+の値で判断)がより重要だ。品質指標は「絶対的な合格点」ではなく、「問題を特定するための目安」と捉えよ。
計算が途中で「浮動小数点例外」で止まります。最初に疑うべきポイントはどこですか?
まず、エラーメッセージの直前に出力されている「警告」を確認せよ。具体的には、1. **物性値の極端な変化**: 特に「理想気体」モデルを使う高速流れで、局所的に温度や圧力が計算上ゼロや負になり、密度が定義できなくなる。2. **悪いメッシュ**: 体積がゼロまたは負の要素(特に非構造メッシュの生成失敗)。3. **不適切な境界条件**: 流入速度が高すぎて、出口条件と矛盾し、計算領域内の圧力が物理的にあり得ない値になる。最初の対策は、**「二風上差分」スキームと「自動タイムステップ」で、より小さなCFL数(初期値5など)で再計算**することだ。これで収束すれば、その後、高分解能スキームに戻せる。
ソフトウェア比較
Ansys内部での位置づけ
AnsysにはCFXとFluentの2つのCFDソルバーがあります。なぜ2つもあるんですか?どちらを使うべきか迷います。
歴史的経緯(別会社の買収)もあるが、現在は用途である程度棲み分けがされている。**CFX**は圧力基盤、要素中心の有限体積法、強カップリングソルバーが特徴。ターボ機械(ポンプ、ファン、タービン)の解析には専用の「ターボモード」と「フローズンロータ」インターフェースが強く、業界でデファクトスタンダードだ。また、粒子追跡や多相流(特に欧拉法的アプローチ)も堅牢。**Fluent**は密度基盤(圧力基盤も選択可)、セル中心の有限体積法、分離型ソルバーが主流。自動車外気解析、燃焼、化学反応、複雑な動き(6DOF)のモデリングに強みがある。迷ったら、ターボ機械と複雑な多相流はCFX、その他はFluent、と覚えておけば大体間違いない。
「Workbench」と「CFX-Pre/Post」の関係がわかりません。どちらを使えばいいですか?
**Ansys Workbench**は統合プラットフォームだ。ここで「CFX」コンポーネントを起動すると、その中でCFX-Pre(前処理)、CFX-Solver Manager(求解制御)、CFX-Post(後処理)が連携して動く。Workbenchの利点は、Geometry (SpaceClaim/DM) や Meshing (Fluent Meshing/Tetra) とデータがシームレスに連携し、パラメータスタディや構造解析(Mechanical)との連成が容易なこと。一方、**スタンドアロンのCFX-Pre/Post**は、Workbenchを介さず直接起動する。特にCFX-Preは、ターボ機械用のテンプレートや高度な式(CEL)の編集など、細かい設定に熟練者が好む傾向がある。初心者はWorkbenchのフローに沿うのが無難で、特にメッシュをAnsys Meshingで作るならWorkbench必須だ。
他のCAEソフト、例えばSiemensのSTAR-CCM+やAltairのAcuSolveと比べて、CFXの強みは何ですか?
最大の強みは「**ターボ機械解析のエコシステム**」の完成度の高さだ。CFXには「ターボモード」があり、ブレードの幾何学データから自動的に流路メッシュを生成し、フローズンロータ、ミキシングプレーン、ステージインターフェースを簡単に設定できる。性能曲線(圧力-流量、効率-流量)の自動計算や、ISO 1940などのバランス規格に基づくアンローディング解析も標準装備。STAR-CCM+も強力だが、ターボ専門家の間ではCFXのノウハウと信頼性は別格だ。また、Ansysエコシステム内での**多物理場連成**(Maxwellとの電磁流体連成、Mechanicalとの流体構造連成FSI)も強み。AcuSolveはGPGPUによる高速計算が売りだが、物理モデルのライブラリの豊富さではAnsysに分がある。
トラブルシューティング
よくあるエラーと対策
「The density ratio has exceeded 1.0e+05」というエラーが出ます。これは何が原因で、どう直せばいいですか?
これは「理想気体」モデルを使用した圧縮性流れ解析で頻発するエラーだ。局所的な流速が音速を超えたり(マッハ数>1)、逆に極端に低速になったりして、連続の式を解く過程で密度の比が計算上の限界を超えたことを意味する。対策は以下の順で試せ。1. **初期条件の見直し**: 全圧、静圧、温度を現実的な値に設定。特に全体の初期圧力を、最も低い出口静圧より高く設定する。2. **境界条件の緩和**: 流入境界を「亜音速流入」から「入口質量流量」や「入口総圧」に変更し、出口を「亜音速流出」(静圧指定)にする。3. **ソルバー設定**: 「二風上差分」スキームと、非常に小さなCFL数(1以下)で計算を再開し、流れ場が安定してから高分解能に戻す。4. **メッシュ改善**: 流速や圧力が急変する領域(ノズル喉部など)のメッシュを細かくする。
非定常計算で「Global Time Step is too large.」と出て止まります。タイムステップを小さくするしかないですか?
直接的な対策はタイムステップを小さくすることだが、根本原因を探る必要がある。このエラーは、あるタイムステップ内での物理量(圧力、速度)の変化がソルバーの内部変数の許容範囲を超えた時に出る。原因は、1. **急激な物理現象**: バルブの急閉など。対策は現象に合った十分に小さいタイムステップを設定する。2. **悪いメッシュ**: 極端に小さい要素や歪んだ要素があると、そこでの計算が不安定になる。メッシュ品質を改善する。3. **不適切な境界条件/初期条件**: 現実とかけ離れた条件から計算を始めている。まず定常計算(RANS)で妥当な流れ場を得て、それを非定常の初期条件として使う「定常→非定常」のワークフローが有効だ。また、CFXの「自動タイムステップ」機能は、このエラーを回避するためにステップを自動で調整してくれるので活用せよ。
計算は最後まで終わるのですが、後処理で壁面近くの流速が異常に高かったり、逆にゼロになっていたりします。何が起きているのでしょうか?
これは典型的な「**壁面関数とメッシュ解像度の不一致**」による問題だ。CFXのデフォルトの壁面処理は「スケーラブルな壁面関数」で、これは最初の要素の中心が壁面から
並列計算で「A partition has encountered an error.」と出て、計算が一部のコアだけで進みます。どう対処すればいいですか?
関連トピック
なった
詳しく
報告