熱解析の収束失敗

いい質問だ。収束しないとは、離散化した支配方程式の残差が、設定した許容値（例えば1.0e-6）以下に減らない状態だ。物理的には、モデル内のどこかで熱平衡が成り立つ解が存在しない、または計算過程で発散していることを意味する。例えば、発熱量が放熱能力を大幅に上回り、温度が無限に上昇するような非現実的な状況を計算しようとしている可能性がある。

非定常熱伝導の支配方程式は、エネルギー保存則から導かれる。体積Vを持つコントロールボリュームで考えると、蓄積項 = 流入項 + 発生項だ。これを微分形で書くと、

となる。ソルバーは、この式を離散化した代数方程式系の残差ノルムを各反復で監視している。Ansys Mechanical APDLでは`CNVTOL, HEAT`コマンドでこの許容値を設定する。

非線形性が強くなるからだ。例えば、ある絶縁材料の`k`が温度の2乗に比例して増加する場合、方程式は`k(T)`の項を含む非線形方程式になる。ソルバーは反復計算で解を探すが、初期推定値と実際の解が離れすぎていると、反復ごとの解の更新量が発散したり、振動したりする。Abaqusでは、このような非線形問題に対してデフォルトで最大10回の反復（`*CONTROLS, ANALYSIS=DISCONTINUOUS`で変更可能）を行うが、そこで収束しなければ失敗と判定される。

これは重要だ。非定常解析では、ある時間ステップ`Δt`の中での平衡を求める「反復」を行う。この時間ステップ`Δt`が大きすぎると、一つのステップ内での温度変化が急激になり、反復計算が収束しなくなる。まずいじるべきは「時間増分」、つまり`Δt`を小さくすることだ。Ansys Fluentのデフォルトの時間ステップは流れの時間スケールから決まるが、熱伝導のみの場合は、材料の熱拡散率αから推定される拡散時間スケール

根本的に異なる。直接法（如きMUMPS, PARDISO）は離散化された連立一次方程式を有限回の演算で「正確に」解く。ここで失敗するのはメモリ不足や特異行列（条件が定義されていない節点がある等）が主因だ。一方、反復法（如き共役勾配法：CG）は残差を反復ごとに減らしていく。これが失敗するのは、問題が悪条件（材料特性の比が極端など）であるか、許容反復回数が少なすぎる場合だ。COMSOLでは、定常熱解析で「分離」ソルバー（反復法）を使う場合、プレコンディショナーの設定（幾何学的マルチグリッドなど）が収束性を大きく左右する。

反復計算における解の更新量を制御する係数だ。反復iからi+1への温度更新を

以下の順で確認せよ。1. **境界条件の漏れ**: 全ての外面に熱的境界条件（温度、熱流束、対流、放射）が定義されているか。未定義面は「断熱」とみなされるが、意図しない断熱面がないか。2. **材料特性の入力ミス**: 熱伝導率`k`に1.0e-6のような極端に小さい値や、密度`ρ`に0が入っていないか。単位系（SI vs mm単位系）の不一致はないか。3. **接触熱抵抗の設定**: 部品間の接触面で「熱接触」が定義されていないと、熱が伝わらず温度が局所的に発散する。

次に、物理的に解が存在するか疑え。例えば、密閉筐体内のCPUの発熱が10W、筐体外面の自然対流（熱伝達率5 W/m²K）のみで放熱している場合、平衡温度は数百℃に達し、材料特性がその温度で大きく変化するかもしれない。まず、発熱量と放熱量のバックオブザエンベロープ計算をせよ。発熱量`Q_gen`が放熱量`Q_loss = h A (T_s - T_amb)`を常に上回る設定なら、物理的に定常解は存在しない。発熱量を半分にするか、強制空冷（hを20 W/m²K以上に）を検討する必要がある。

その時間で何らかの「事象」が発生しているからだ。具体的には、1. **相変化（融解・凝固）**: 潜熱の吸収・放出により熱容量が実質的に無限大になり、ソルバーが平衡温度を見つけられない。2. **温度依存材料特性の不連続点**: データテーブルの入力間隔が粗く、ある温度を境に`k`の値がジャンプしている。3. **接触状態の変化**: 熱接触のギャップが閉じる/開くなど、システムの剛性が劇的に変化する。対策としては、事象が起こる時間付近で自動時間ステップ制御を有効にし、最大反復回数を増やす（例えば25回に）ことだ。

「Analysis Settings」の以下の項目だ。1. **Solver Type**: プログラム制御（通常は直接法）から「Iterative (PCG)」に変えてみる。大規模モデルで有効。2. **Nonlinear Controls**: 「Force Convergence」と「Momentum Convergence」の許容値（Tolerance）をデフォルトの0.5%から1%や2%に緩和する。物理的に妥当か結果を確認する必要がある。3. **Line Search**: 「On」にすると更新方向を調整し、強非線形問題で収束を助ける。4. **Stabilization**: 「Energy」スタビライゼーションを「On」にし、係数を1.0e-3など小さな値から試す。これは数値的なダンピングとして働く。

根本的に異なる。**Standard**はニュートン・ラフソン法による反復で平衡を求める。収束しない場合は、`*CONTROLS`パラメータで、許容残差(`c1`)、最小反復回数(`i1`)、最大反復回数(`i2`)を調整する。例えば`*CONTROLS, PARAMETERS=TIME INCREMENTATION`で`, i1=8, i2=15`のように増やす。一方、**Explicit**は動的陽解法で、時間積分を行うため「収束」の概念自体がない。熱伝導も、熱伝導行列を対角化した（lumped）熱容量行列を用いて陽的に解く。安定性は安定時間ステップ

COMSOLの分離ソルバーは、物理場ごとに方程式を順次解く反復法（Segregated法）だ。収束失敗時は、1. **ソルバーシーケンスの調整**: 「完全従属ソルバー」に切り替える。これは全物理場を連立して解く直接法で、確実だがメモリ消費大。2. **緩和係数の調整**: 分離ソルバーの「高非線形性」設定を使用するか、手動で「ダンピング係数」を0.5以下に下げる。3. **初期値の改善**: 「補助掃引」を使用し、あるパラメータ（例えば発熱量）を0から最終値まで段階的に増加させながら解を求める。これにより、各ステップで良い初期値が得られる。特に放射熱伝達を含む問題で有効だ。

それは最後の手段だ。反復回数が30回でも収束しないということは、解に向かって収束する「傾向」すら見えていない可能性が高い。まず、ソルバー出力（.logファイルや.outputファイル）を確認し、残差ノルムの履歴を見よ。単調減少しているなら回数を増やす価値がある（Ansysでは`NEQIT`コマンド）。しかし、振動したり、一定値で停滞したりしているなら、根本原因（境界条件、材料特性、接触）を疑うか、緩和係数を下げるべきだ。回数だけ増やすと、無駄な計算時間を浪費するだけだ。

これは重大なエラーで、通常は材料特性の入力ミスか、極端な非線形性による数値的不安定性が原因だ。熱伝導マトリックスの対角項は、その節点からの熱の流れ出しやすさを示す正の値になるはずだ。これが負になるのは、1. **熱伝導率`k`に負の値が入力されている**（あり得ない）。2. **温度依存性のデータテーブルで、外挿が起こり`k`が負になっている**。3. **放射の形状係数（View Factor）の計算が不正確**で、実質的に負の熱抵抗が生じている。まず、全ての材料の`k`が全ての温度範囲で正であることを確認せよ。Abaqusでは、`*CONDUCTIVITY`の温度依存データの範囲を広げ、`*EXPANSION`と混同していないか確認する。

出る。これは陽解法（Explicit）だけでなく、陰解法の自動時間ステップ制御でも発生する。ソルバーが収束するために必要な時間ステップ`Δt_req`が、ユーザーが設定した最小時間ステップ`Δt_min`よりも小さくなったことを意味する。原因は、1. **急激な負荷変化**: ステップ関数で発熱がONになるなど。対策は、負荷を線形に立ち上げる（Amplitude定義）。2. **接触状態の急変**: ギャップ閉鎖の瞬間。接触剛性を適切に設定する。3. **極端に小さい要素が1つある**: 熱拡散の安定条件

放射の熱流束が温度の4乗