Directly Coupled (Monolithic) Conjugate Heat Transfer

Simply put, solving the governing equations for fluid and solid as a single simultaneous equation system is direct coupling (monolithic method). Partitioned coupling (partitioned method) is an approach where the solid solver and fluid solver run separately, exchanging temperature and heat flux at the interface like a game of catch.

Yes, in partitioned coupling, an outer iteration is required: "calculate temperature on the solid side → pass interface heat flux to fluid → fluid calculates temperature → return interface temperature to solid." This back-and-forth repeats until convergence, which increases computation time and can diverge under poor conditions.

Direct coupling eliminates the need for this interface exchange; the solid and fluid temperature fields are determined simultaneously in a single solve of the simultaneous equations. Since temperature continuity and heat flux conservation at the interface are algebraically guaranteed, there are no convergence worries.

There are pros and cons. Direct coupling requires a dedicated solver and the solid and fluid meshes must be integrated. Partitioned coupling is more flexible when you want to combine existing structural and CFD solvers. However, for problems involving thin solid layers like chip cooling on electronic boards, the Biot number tends to become large, and direct coupling often offers better accuracy and efficiency.

First, let's write the heat equations for the solid domain $\Omega_s$ and fluid domain $\Omega_f$.

Solid side (heat conduction only):

Fluid side (convection + diffusion):

In direct coupling, after discretizing these two using the Finite Element Method (FEM) or Finite Volume Method (FVM), they are assembled into a single block matrix system with shared interface nodes. It looks something like this:

$\Gamma$ is a subscript representing the interface between solid and fluid. Nodes on the interface belong to both solid and fluid, and the temperature $\mathbf{T}_\Gamma$ is shared. This is the mechanism that algebraically guarantees "temperature continuity."

In practice, it's often written more compactly as a 2x2 matrix:

Here, $\mathbf{K}_{sf}$ and $\mathbf{K}_{fs}$ are the interface coupling matrices (off-diagonal coupling blocks), representing the exchange of heat through the interface. The difference with partitioned coupling is that these coupling terms don't exist there; instead, data is exchanged via outer iterations.

There are two conditions that must be satisfied at the interface $\Gamma$. First, temperature continuity (Dirichlet condition):

And heat flux conservation (Neumann condition):

$n$ is the normal direction vector at the interface. In direct coupling, since interface nodes are shared, temperature continuity is automatically satisfied by node coincidence. Heat flux conservation is algebraically guaranteed within the weak form (variational formulation) because the interface integrals cancel between adjacent elements.

To use an analogy, partitioned coupling is like "two rooms exchanging temperature over the phone," while direct coupling is like "removing the wall to make one room." In a game of telephone, information discrepancies (interface errors) occur, but if there's no wall, there's no possibility for discrepancy in the first place.

Good question. In partitioned coupling, the "solid → fluid → solid → ..." iteration becomes a Gauss-Seidel-like sequential update. If the ratio of thermal properties between solid and fluid (especially $k_s/k_f$ or $\rho c_p$ ratio) is extreme, the spectral radius of this iteration can exceed 1, leading to divergence in some cases.

In direct coupling, since all degrees of freedom are solved simultaneously, using implicit time integration makes it unconditionally stable. Particularly:

• Cases with large thermal conductivity differences like $k_s \gg k_f$ (metal solid + air fluid)
• Cases with thin solid layers where the Biot number $\text{Bi} = hL/k_s$ is $O(1)$ or greater
• Transient analyses with large time steps

Under these conditions, partitioned coupling iterations struggle to converge, but direct coupling can solve them without issue.

There's a theoretical result shown by Giles (1997): the interface condition number for partitioned coupling is:

Larger $\kappa_\Gamma$ leads to slower convergence. For a combination like steel ($k_s = 50$ W/(m$\cdot$K)) and air ($k_f = 0.026$), the ratio is about 2000 times. The biggest advantage of direct coupling is that it is not affected by this $\kappa_\Gamma$.

The basics are the same as regular FEM/FVM. Discretize the solid and fluid domains separately to create element matrices, then assemble them into the global matrix. The key point is the numbering of interface nodes.

For FEM, the weak form for a solid element is:

For fluid elements, the convection term $\int \rho_f c_{p,f} N_i (\mathbf{u} \cdot \nabla N_j) \, d\Omega$ is added to this.

The trick in assembly is to assign elements on the interface to both the solid and fluid domains. By numbering the interface nodes so they appear only once in the global matrix, temperature continuity is automatically guaranteed. In actual solvers, care must be taken to avoid double-counting the DOFs (degrees of freedom) of shared nodes.

Sharp observation. CFD solvers like STAR-CCM+ or Fluent are FVM-based, but their CHT functionality also handles the solid domain within the same FVM framework. The mechanism ensures heat flux conservation at interface cell faces. In FVM, the "face flux" between adjacent cells corresponds to the interface coupling.

Specifically, the heat flux at an interface face is calculated as:

using a harmonic average. Here, $\delta_s$, $\delta_f$ are the distances from the interface to the solid and fluid cell centers, respectively. Distributing this to the corresponding rows in the solid and fluid matrices completes the monolithic matrix.

Exactly. The monolithic matrix includes all DOFs of solid+fluid, so it's very large. The solution method is chosen based on problem scale:

It's extremely effective. AMG (Algebraic Multigrid) automatically creates a hierarchy of "coarse approximations" and "fine approximations" of the matrix, simultaneously damping low-frequency and high-frequency error components. The monolithic matrix for CHT tends to have a poor condition number because material properties differ by orders of magnitude between solid and fluid, but AMG can significantly mitigate this effect.

In practice, setting the AMG convergence criterion to a residual below $10^{-6}$ is standard.

The most important requirement is mesh conformity (conformal mesh) at the interface. The nodes (FEM) or faces (FVM) on the solid side and fluid side of the interface must match. Some solvers support non-conformal interfaces, but accuracy tends to suffer.

Practical points to note:

• Boundary layer mesh on solid walls: Sufficient prism layers (inflation layers) are needed near the interface on the fluid side. Mesh design based on $y^+$ is essential.
• Solid interior: Place at least 3 or more element layers in thin wall sections where temperature gradients are steep.
• Mesh size ratio: The element size ratio between solid and fluid sides at the interface should ideally be within a factor of 5.

That's true. However, with tools like STAR-CCM+ that can "generate solid+fluid meshes in one go within a single software," creating a conformal mesh is relatively easy. In fact, partitioned coupling can sometimes be more troublesome due to the effort required to "align meshes between two solvers."

有限要素法（FEM）による空間離散化を使うんだ。要素剛性マトリクスを組み立て、全体剛性方程式を構築する。

弱形式（変分形式）への変換を行い、試験関数と形状関数を用いてGalerkin法による定式化を使うんだ。要素タイプの選択（低次要素 vs. 高次要素、完全積分 vs. 低減積分）は解の精度と計算コストのトレードオフに直結するんだよ。

直接法（LU分解、Cholesky分解）または反復法（CG法、GMRES法）により連立方程式を解く。大規模問題では前処理付き反復法が効果的なんだ。

Fluent Inc.が開発。2006年にAnsysが買収。非構造格子ベースの汎用CFDソルバー。

現在の所属: Ansys Inc.

CD-adapcoが開発。2016年にSiemensが買収しSimcenterブランドに統合。ポリヘドラルメッシュが特徴。

現在の所属: Siemens Digital Industries Software

1986年スウェーデンで設立。MATLAB連携のFEMLABとして開始、後にCOMSOLに改名。マルチフィジックスに強み。

現在の所属: COMSOL AB

異なるソルバー間でモデルを変換する際は、要素タイプの対応関係、材料モデルの互換性、荷重・境界条件の表現差異に注意が必要になるんだ。特に高次要素や特殊要素（コヒーシブ要素、ユーザー定義要素等）はソルバー間で直接変換できない場合が多い。

メッシュ収束性の確認、境界条件の妥当性検証、材料パラメータの感度分析がすごく大事なんだ。

• メッシュ依存性の検証: 少なくとも3水準のメッシュ密度で収束性を確認
• 境界条件の妥当性: 物理的に意味のある拘束条件の設定
• 結果の検証: 理論解、実験データ、既知ベンチマーク問題との比較

うん、いい調子だよ！ 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。

形状関数 $N_i$ を用いて未知量を近似:

これを数式で表すとこうなるよ。

これを数式で表すとこうなるよ。

連続体の支配方程式を離散化すると、以下の代数方程式系が得られる:

ここで $[K]$ は全体剛性マトリクス（または同等のシステムマトリクス）、$\{u\}$ は未知節点変数ベクトル、$\{F\}$ は外力ベクトルなんだ。

• 完全積分: 全ての項を正確に積分。剛性過大評価の傾向（ロッキング）
• 低減積分: 積分点数を削減。計算効率向上だが、アワーグラスモード発生のリスク
• 選択的低減積分 (B-bar法): 体積項と偏差項を分離して積分。ロッキング回避

• h-refinement: メッシュを細分化（要素サイズ h を小さく）して精度向上
• p-refinement: 要素の多項式次数を上げて精度向上
• hp-refinement: h と p を同時に最適化

収束速度: 二次要素で $O(h^2)$ のオーダーで誤差が減少（滑らかな解の場合）

直接カップリング法の実務的な解析フローと注意点を解説する。

1. 前処理 (Pre-processing)

• CADデータのインポートと形状簡略化
• 材料特性の定義
• メッシュ生成（要素タイプ・サイズの決定）
• 境界条件と荷重条件の設定

2. 求解 (Solving)

• ソルバー設定（解法、収束基準、出力制御）
• ジョブ投入と計算実行
• 収束モニタリング

3. 後処理 (Post-processing)

• 結果の可視化（変位、応力、その他の物理量）
• 結果の検証と妥当性確認
• レポート作成

• 応力集中部: 最低3層以上の要素を配置
• 応力勾配の大きい領域: 要素サイズを周囲の1/3〜1/5に
• 荷重印加点近傍: 局所細分化
• 遠方領域: 粗いメッシュで計算効率を確保

• 過拘束に注意: 剛体移動の拘束は6自由度のみ
• 対称条件の活用: 計算規模の削減
• 荷重の等価分配: 集中荷重 vs. 分布荷重の選択

Fluent Inc.が開発。2006年にAnsysが買収。非構造格子ベースの汎用CFDソルバー。

現在の所属: Ansys Inc.

CD-adapcoが開発。2016年にSiemensが買収しSimcenterブランドに統合。ポリヘドラルメッシュが特徴。

現在の所属: Siemens Digital Industries Software

• メッシュ収束性を3水準以上で確認したか
• 力の釣り合い（反力合計）を検証したか
• 結果が物理的に妥当な範囲か確認したか
• 既知の理論解またはベンチマーク問題と比較したか

うん、いい調子だよ！ 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。

直接カップリング法に対応する主要な商用CAEツールの機能比較と、各製品の歴史的背景を詳述する。

Fluent Inc.が開発。2006年にAnsysが買収。非構造格子ベースの汎用CFDソルバー。

現在の所属: Ansys Inc.

CD-adapcoが開発。2016年にSiemensが買収しSimcenterブランドに統合。ポリヘドラルメッシュが特徴。

現在の所属: Siemens Digital Industries Software

1986年スウェーデンで設立。MATLAB連携のFEMLABとして開始、後にCOMSOLに改名。マルチフィジックスに強み。

現在の所属: COMSOL AB

1970年にSwanson Analysis Systems Inc. (SASI) が開発。APDL（Ansys Parametric Design Language）ベース。

現在の所属: Ansys Inc.

• 要素タイプの非互換: ソルバー固有要素は中立フォーマットで表現不可
• 材料モデルの差異: 同名でも内部実装が異なる場合がある
• 境界条件の再定義: 多くの場合、手動での再設定が必要
• 結果データの比較: 出力変数の定義（節点値 vs. 要素値、積分点値）に差異

直接カップリング法のツール選定においては以下を考慮:

• 解析規模: 数万〜数億DOFへのスケーラビリティ
• 物理モデル: 必要な構成則・要素タイプの対応状況
• ワークフロー: CADとの連携、自動化の容易さ
• コスト: 初期投資 + 年間保守 + 教育コスト
• サポート: 技術サポートの質とレスポンス

うん、いい調子だよ！ 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。

直接カップリング法における最新の研究動向と先進的手法を見ていこう。

• 等幾何解析 (IGA): NURBS基底関数を直接使用し、CAD-CAE間のシームレスな連携を実現
• 粒子法 (SPH, MPM): メッシュフリー手法による大変形・破壊の追跡
• 位相場法 (Phase-Field): 界面の暗示的表現による複雑な界面追跡
• 機械学習支援: サロゲートモデル、物理インフォームドニューラルネットワーク (PINN)

症状: ソルバーが指定反復回数内に収束せず異常終了

考えられる原因:

• メッシュ品質の不足（過度に歪んだ要素）
• 材料パラメータの不適切な設定
• 不適切な初期条件
• 非線形性が強すぎる（荷重ステップの不足）

対策:

• メッシュ品質チェックを実施（アスペクト比、ヤコビアン）
• 材料パラメータの単位系を確認
• 荷重を複数ステップに分割（サブステップ数の増加）
• 収束判定基準の緩和（ただし精度に注意）

症状: 応力/変位/温度等が物理的に非現実的な値

考えられる原因:

• 境界条件の誤設定
• 単位系の混在（SI単位と工学単位の混同）
• 不適切な要素タイプの選択
• 応力特異点の存在

対策:

• 反力の合計を確認（力の釣り合い）
• 単位系の一貫性を確認
• 要素タイプの適切性を再検討
• 特異点除去またはサブモデリング

症状: 計算が想定時間の何倍もかかる

対策:

• メッシュの粗密分布の最適化
• 対称性の活用（1/2, 1/4モデル）
• ソルバー設定の最適化（反復法、前処理の選択）
• 並列計算の活用

症状: Out of Memory エラー

対策:

• アウトオブコア解法の使用
• メッシュ規模の削減
• 64bit版ソルバーの使用確認
• メモリ割り当ての増加

• FATAL 2012: 特異剛性マトリクス → 拘束条件の見直し
• USER WARNING 5291: 要素品質不良 → メッシュ修正
• SYSTEM FATAL 3008: メモリ不足 → MEM設定の調整

• Excessive distortion: 要素の過大変形 → NLGEOM確認、メッシュ改善
• Zero pivot: 拘束不足 → 境界条件追加
• Time increment too small: 収束失敗 → ステップ設定見直し