hexagon logo

圧縮性流体解析の際、FE106で計算が止まってしまいました。対処方法を教えてください。

回答​
圧縮性流体解析の際、現実ではありえない負の密度が出力されるとFE106で解析は中断いたします。この場合、以下の対処方法をお試しください。
 
----------------------------------------
メッシュの品質を確認する
----------------------------------------
密度が負になった節点を確認し、その周辺の要素の品質を「h-ratio」および「Shape Factor」で確認してください。もし、品質の悪い要素が見つかった場合は該当箇所のオクタントを細かくして再度メッシュを作成してください。
 
----------------------------------------
緩和係数を小さくする (定常計算の場合)
----------------------------------------
流速・圧力・温度に関して、緩和係数をデフォルト値から0.1刻み程度で徐々に下げていき計算が安定するかお試しください。
 
----------------------------------------
慣性不足緩和を小さくする (定常計算の場合)
----------------------------------------
流速・圧力・温度に関して、慣性不足緩和を“1”程度に下げて解析すると安定するかお試しください。もし安定した場合、収束が早くなるよう徐々に慣性不足緩和を上げていってください。
 
----------------------------------------
流入流出境界条件(流速、圧力)をサイクル毎に緩やかに変化させる(定常計算の場合)
----------------------------------------
解析開始時の急激な流れの発生は計算の不安定性に結びつくため、 変数テーブル機能を使い、数百サイクルをかけて徐々に目的とする状態に遷移させると安定して解析が行える場合があります。
 
----------------------------------------
初期値を設定する (定常計算の場合)
----------------------------------------
圧縮性流体解析ではデフォルトで調和関数から初期値が与えられますが、調和関数を利用せず初期値を直接設定すると安定に解析出来る場合があります。
 
----------------------------------------
倍精度版ソルバーを使用する
----------------------------------------
倍精度版ソルバーを使用すると、数値誤差が蓄積しにくくなり計算の安定性に結びつきます。