直接耦合（整体式）共轭传热

简单来说，将流体和固体的控制方程作为一个联立方程组同时求解的就是直接耦合（整体法）。分离型（分区法）则是将固体求解器和流体求解器分开运行，在界面上交换温度和热通量的方式。

是的，分离型需要“固体侧计算温度 → 将界面热通量传给流体 → 流体计算温度 → 将界面温度返回固体”这样的外部迭代。这个过程需要来回多次直到收敛，所以计算时间长，条件不好时还会发散。

直接耦合则不需要界面上的数据交换，一次联立方程求解就能同时确定固体和流体的温度场。界面的温度连续性和热通量守恒在代数上得到保证，所以不用担心收敛问题。

各有利弊。直接耦合需要专用求解器，并且必须将固体和流体网格一体化。如果想组合使用现有的结构求解器和CFD求解器，分离型则更灵活。不过，对于包含薄固体层的问题（如电子基板上的芯片冷却），毕渥数容易变大，直接耦合通常在精度和效率上都更好。

首先，写出固体区域 $\Omega_s$ 和流体区域 $\Omega_f$ 各自的热方程。

固体侧（仅热传导）:

流体侧（对流+扩散）:

直接耦合是将这两个方程用有限元法（FEM）或有限体积法（FVM）离散化后，组装成共享界面节点的单一分块矩阵系统。会变成这样：

$\Gamma$ 是表示固体与流体界面的下标。界面上的节点同时属于固体和流体，共享温度 $\mathbf{T}_\Gamma$。这就是代数上保证“温度连续性”的机制。

实际上，也常用更紧凑的2行2列形式表示：

这里的 $\mathbf{K}_{sf}$ 和 $\mathbf{K}_{fs}$ 是界面耦合矩阵（非对角耦合块），表示通过界面的热量传递。分离型则没有这个耦合项，取而代之的是通过外部迭代进行数据交换，这就是区别。

在界面 $\Gamma$ 上需要满足两个条件。首先是温度连续性（Dirichlet条件）：

然后是热通量守恒（Neumann条件）：

$n$ 是界面的法线方向向量。直接耦合由于共享界面节点，温度连续性通过节点一致自动满足。热通量守恒则通过弱形式（变分公式）中相邻单元间的界面积分相互抵消，在代数上得到保证。

打个比方，分离型像是“两个房间之间通过电话互相告知温度”，而直接耦合则是“拆掉墙壁变成一个房间”。电话传话游戏会产生信息偏差（界面误差），但如果没有墙壁，就根本不会产生偏差。

问得好。分离型中“固体→流体→固体→…”的迭代是类似Gauss-Seidel的逐次更新。当固体和流体的热物性比（特别是 $k_s/k_f$ 或 $\rho c_p$ 的比值）极端时，这种迭代的谱半径可能超过1，导致发散。

直接耦合同时求解所有自由度，因此如果使用隐式时间积分，就能实现无条件稳定。特别是：

• $k_s \gg k_f$（金属固体 + 空气流体）这种导热率差异大的情况
• 薄固体层中毕渥数 $\text{Bi} = hL/k_s$ 达到 $O(1)$ 以上的情况
• 非稳态分析中时间步长较大的情况

在这些条件下，分离型迭代难以收敛，但直接耦合可以顺利求解。

Giles (1997) 提出了一个理论结果，分离型的界面条件数（interface condition number）为：

$\kappa_\Gamma$ 越大，收敛越慢。钢（$k_s = 50$ W/(m$\cdot$K)）和空气（$k_f = 0.026$）的组合，比率可达约2000倍。直接耦合则不受这个 $\kappa_\Gamma$ 的影响，这是其最大的优点。

基本和通常的FEM/FVM一样。分别离散化固体区域和流体区域生成单元矩阵，再组装成整体矩阵。关键在于界面节点的编号。

FEM的情况下，固体单元的弱形式为：

流体单元则在此基础上增加对流项 $\int \rho_f c_{p,f} N_i (\mathbf{u} \cdot \nabla N_j) \, d\Omega$。

组装的技巧在于，将界面上的单元同时归属于固体和流体两个区域。只要确保界面节点在整体矩阵中只出现一次进行编号，温度连续性就能自动保证。在实际求解器中，需要注意共享节点的自由度（DOF）不要被重复计算。

很敏锐的观察。STAR-CCM+、Fluent等CFD求解器虽然是基于FVM的，但其CHT功能也在同一个FVM框架下处理固体区域。其机制保证了界面单元面上的热通量守恒。在FVM中，相邻单元间的“面通量”相当于界面耦合。

具体来说，界面面的热通量计算如下：

这里使用了调和平均。$\delta_s$, $\delta_f$ 是界面到固体、流体单元中心的距离。将其分配到固体矩阵和流体矩阵的对应行中，整体矩阵就完成了。

没错。整体矩阵包含固体+流体的全部自由度，所以非常大。根据问题规模选择解法：

非常有效。AMG（代数多重网格）自动将矩阵的“粗略近似”和“精细近似”分层，同时衰减低频和高频误差分量。CHT的整体矩阵由于固体和流体的物性值差异巨大，条件数容易变差，但AMG可以大幅缓解这种影响。

实际应用中，通常将AMG的收敛判据设为残差 $10^{-6}$ 以下。

最重要的是界面处的网格匹配（共形网格）。界面处固体侧和流体侧的节点（FEM）或面（FVM）必须一致。有些求解器也支持非共形界面，但精度容易下降。

实际注意事项：

• 固体壁面的边界层网格：流体侧界面附近需要有足够的棱柱层（膨胀层）。必须根据 $y^+$ 进行网格设计
• 固体内部：在温度梯度大的薄壁部分，至少配置3层以上的单元
• 网格尺寸比：界面处固体侧与流体侧的单元尺寸比建议在5倍以内

是的。不过，像STAR-CCM+这样“在一个软件内统一生成固体+流体网格”的工具，创建共形网格相对容易。相反，分离型有时在“协调两个求解器间的网格”上更麻烦。

有限要素法（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: 収束失敗 → ステップ設定見直し