モノリシック（一体型）CHT解法

根本的な違いは、支配方程式の離散化と解法の段階にあります。セグリゲート法は流体と固体の方程式を別々のソルバーで交互に解き、境界でデータを交換します。一方、モノリシック法は最初から一つの連立方程式系として定式化します。例えば、非定常熱伝導と流体のエネルギー方程式を、界面での温度と熱流束の連続条件を内部境界条件として直接組み込んだ形です。これにより、界面での反復収束遅延や、小さなタイムステップでの不安定性が本質的に軽減されます。

良い質問です。変数は統一します。多くの実装では、全体の計算領域（流体+固体）で「温度」を一次変数として採用します。流体領域ではナビエ-ストークス方程式と連続の式に加え、エネルギー方程式を温度形式で解く。固体領域では熱伝導方程式を解く。これらを一つのマトリックスに組み上げると、提示されたようなブロック構造になります：

はい、界面要素（例えば流体と固体に跨がる要素）の離散化において、弱形式の境界積分項を適切に扱うことで実現されます。具体的には、界面ガンマ上で、

まさにそこが実装の肝です。全体マトリックスは非対称かつ大規模になります。直接法は現実的ではないため、Krylov部分空間法（例えばGMRES）と前処理の組み合わせが必須です。具体的な前処理として、流体ブロック

それがモノリシック法の大きな利点の一つです。一つの連立方程式系として扱うため、陰的積分法（例えば後退差分公式BDF2）を全体に適用できます。これにより、流体のCFL条件や固体のフーリエ数に縛られず、物理現象に基づいたタイムステップを選択可能になります。例えば、自動車エンジン部品の熱応答解析では、燃焼サイクル（数ms）とブロックの熱浸透（数秒）を同時に扱う必要がありますが、モノリシック法ならタイムステップを1msから10ms程度で設定し、安定して解けます。Abaqus/Standardの共役熱伝導解析も、同様の陰的時間積分を採用しています。

必ずしも一致メッシュは必要ありません。非一致メッシュを扱うためには、界面での積分において、一方のメッシュ（例えば流体側）の節点での温度や熱流束を、他方のメッシュ（固体側）の形状関数を使って補間する「マッピング」技術を用います。これにより、

まずは「物性値の桁数」を確認してください。特に熱伝導率と体積熱容量です。空気の熱伝導率は約0.026 W/mK、鋼は約50 W/mKで、3桁以上違います。モノリシック法ではこれがマトリックス内に直接入るため、条件数が悪化し、ソルバーが発散する原因になります。対策として、無次元化を行うか、ソフトウェアのスケーリングオプションを有効にします。ANSYS Mechanical APDLでCHTを行う場合、`MP`コマンドで物性を入力する際、単位系の一貫性を厳密にチェックする必要があります。

その通りです。最低限、以下の2点を定量チェックします。

界面を含めて一括で行うべきです。なぜなら、モノリシック法の精度は界面近傍のメッシュ解像度に強く依存するからです。具体的な手順としては、まず境界層メッシュ（流体側）と、固体の熱浸透深さを考慮したメッシュを適切に生成します。その後、全体のメッシュサイズを系統的に細かくし（例えば基底サイズを1.0, 0.7, 0.5倍）、監視対象（例えば界面の平均熱伝達率や固体最高温度）の変化が例えば2%以内に収まることを確認します。Altair AcuSolveのチュートリアルでは、CHTのメッシュ依存性検証として、界面近傍の要素アスペクト比を5以下に保つことを推奨しています。

真のモノリシック解法を実装しているのは、現状では限られています。COMSOL Multiphysicsは「係数型偏微分方程式」インターフェースで流体と固体の方程式を一つの系として記述するため、本質的にモノリシックです。一方、ANSYS Fluentの「Coupled Thermal Boundary Condition」は、エネルギー方程式のみを界面で結合し、運動量方程式とは分離して解く「部分的モノリシック」と言えます。見分け方は、ソルバーの出力ログを見て、流体と固体の温度場が「1つの線形ソルバー呼び出し」で更新されているか、それとも「外部/内部反復」として別々に呼び出されているかを確認することです。

OpenFOAMの`chtMultiRegionFoam`は、デフォルトでは**セグリゲート法**です。流体ソルバー（`buoyantSimpleFoam`等）と固体ソルバー（`solidTemperature`）を別々に呼び出し、界面で温度と熱流束を交換する「外部反復」ループを回しています。ただし、開発バージョンや派生ソルバーには、エネルギー方程式部分をモノリシックに解く実験的な実装も存在します。完全なモノリシック法を求めるなら、`fe40`や`CalculiX`のような有限要素法ベースのオープンソースコードで、自分でカスタム要素を実装する必要があるでしょう。

それら単体では「外部のCFDソフトとの連携」という形になります。Abaqus/Standardには「共役熱伝導解析」機能がありますが、これは複数の固体領域間の熱伝導を扱うもので、流体は含みません。流体を含むCHTを行うには、Abaqus CFD（旧Dassault Systèmesの`PowerFLOW`連携）や、ANSYSの場合は`System Coupling`を使ってFluent（流体）とMechanical（固体）を双方向に結合させます。これはソフトウェア間のセグリゲート法であり、計算コストとデータ転送のオーバーヘッドが大きいという欠点があります。

最も多い原因は「初期場の不整合」です。特に、固体の初期温度を一様な室温（20°C）に設定した場合、流入する高温流体（例えば500°C）と界面で極端な温度勾配が生じ、最初のタイムステップで計算された熱流束が物理的にあり得ない巨大な値になり、ソルバーが破綻します。対策は、流体の流入温度に近い値で固体の初期温度を設定するか、最初の数十ステップは非常に小さなタイムステップ（例えば1e-6秒）でゆっくりと物理場を立ち上げる「ランプ関数」を適用することです。Siemens Star-CCM+の`Field Function`を使えば、このような初期化プロファイルを簡単に設定できます。

これは、流体と固体の熱的インピーダンス（熱拡散率の比）が極端に異なる場合に起こる、数値的な「スティフネス」の問題です。空気（熱拡散率 α ≈ 2e-5 m²/s）と銅（α ≈ 1e-4 m²/s）ならまだ良いですが、空気と樹脂（α ≈ 1e-7 m²/s）だと問題が顕著になります。対策は二つあります。1) 時間積分スキームをより安定なものに変える（例えば、BDF1からBDF2へ）。2) ソルバーの収束基準を厳しくする。特に、COMSOLでは「相対許容誤差」をデフォルトの0.01から0.001以下に引き下げることで、この振動が抑制されることが多いです。根本的には、界面近傍のメッシュを、より熱拡散率の低い（=熱が伝わりにくい）材料側で細かくする必要があります。

その現象は、ソルバーが「局所的な」質量・運動量・エネルギーの保存則の残差は小さくしているが、「グローバルな」熱バランスが達成されていない状態を示しています。原因は主に二つ考えられます。1) **緩和係数が強すぎる**: 特に固体の熱伝導方程式に強いアンダーリラクゼーションをかけていると、更新量が小さすぎて全体の熱平衡に至るまでの反復回数が膨大になります。2) **境界条件の漏れ**: 固体の外部断熱境界が完全でない、または流体の流出境界で熱流束がゼロでない（対流流出など）。対策として、グローバルな熱収支（流入熱量 = 固体の内部発熱 + 系からの損失）をモニターし、これが収束するまで計算を続けるようにします。ANSYS Fluentでは、`Surface Monitor`で境界の積分熱流束を監視できます。

それはマッピング誤差による「数値的な熱流束の漏れ」の典型的な症状です。許容できません。非一致メッシュを用いる場合、以下の点を確認してください。