共役热传递分析
共役热传递的理论基础
概要
老师！今天是讲共役热传递分析的话题吧？那是什么东西呢？
固体和流体的热传递同时求解。电子设备冷却和涡轮叶片的热设计。界面处的温度热流通量连续条件。
等等，固体和流体的热传递同时解，也就是说可以在这样的情况下也使用吗？
支配方程
哇～，共役热传递分析的记述话题，超级有意思！请告诉我更多。
离散化方法
这个方程，在计算机上实际是怎样求解的呢？
使用有限元法（FEM）进行空间离散化。组装单元刚度矩阵，构建整体刚度方程。
矩阵求解算法
矩阵求解算法，具体来说是什么意思呢？
用直接法（LU分解、Cholesky分解）或迭代法（CG法、GMRES法）求解联立方程组。大规模问题中带预处理的迭代法很有效。
| 求解法 | 分类 | 内存使用量 | 适用规模 |
|---|---|---|---|
| LU分解 | 直接法 | O(n²) | 小～中规模 |
| Cholesky分解 | 直接法（对称正定值） | O(n²) | 小～中规模 |
| PCG法 | 迭代法 | O(n) | 大规模 |
| GMRES法 | 迭代法 | O(n·m) | 大规模、非对称 |
| AMG预处理 | 预处理 | O(n) | 超大规模 |
也就是说，在有限元法的部分省略的话，之后会很痛苦吧。我铭记于心！
商用工具中的实现
那么，进行共役热传递分析的话用什么软件呢？
| 工具名 | 开发方/现在 | 主要文件格式 |
|---|---|---|
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
| COMSOL Multiphysics | COMSOL AB | .mph |
| ANSYS Mechanical（原ANSYS Structural） | Ansys Inc. | .cdb, .rst, .db, .ans, .mac |
供应商系列与产品整合历史
各个软件的成立历史，相当具有戏剧性吗？
Ansys Fluent
接下来是Ansys Fluent的话题。内容是什么呢？
由Fluent Inc.开发。2006年被Ansys收购。基于非结构化网格的通用CFD求解器。
现在的所属：Ansys Inc.
Simcenter STAR-CCM+
接下来是Simcenter STAR的话题。内容是什么呢？
由CD-adapco开发。2016年被Siemens收购并整合到Simcenter品牌下。多面体网格是其特点。
现在的所属：Siemens Digital Industries Software
听到这里，才终于理解了为什么开发历史很重要！
COMSOL Multiphysics
请给我讲讲「COMSOL Multiphysics」！
1986年在瑞典成立。最初作为MATLAB联动的FEMLAB开始，后来改名为COMSOL。在多物理场方面有强项。
现在的所属：COMSOL AB
哇～，关于开发的话题，超级有意思！请告诉我更多。
文件格式与互操作性
在不同软件之间交换数据时有什么要注意的呢？
| 格式 | 扩展名 | 种别 | 概要 |
|---|---|---|---|
| STEP | .stp/.step | 中立CAD | ISO 10303符合的3D CAD数据交换格式。形状+PMI对应。 |
| IGES | .igs/.iges | 中立CAD | 初期的CAD数据交换规范。曲面数据的互换性有课题。正在向STEP过渡。 |
| CGNS | .cgns | CFD数据 | CFD通用符号系统。CFD结果的标准交换格式。 |
| VTK | .vtk/.vtu | 可视化 | 可视化工具包格式。在ParaView等中使用。 |
在不同求解器之间转换模型时，需要注意单元类型的对应关系、材料模型的兼容性、荷载边界条件的表示差异。特别是高阶单元或特殊单元（粘聚单元、用户定义单元等）在求解器之间往往无法直接转换。
明白了……格式看似简单，但实际上非常深奥啊。
实务注意事项
教科书里没有的「现场智慧」这样的东西有吗？
网格收敛性的确认、边界条件的妥当性检证、材料参数的敏感性分析非常重要。
哎呀，共役热传递分析真的很深奥啊……但是多亏了老师的说明，我理解得差不多了！
嗯，加油啊！实际动手做是最好的学习方法。有不明白的地方随时问我。
Biot数诉说的「固体和流体哪边占主导」
共役热传递的本质浓缩在Biot数（Bi = hL/k）中。若Bi≪1，则固体内部的温度分布均匀，流体侧的热传递为速率限制——电子器件的散热器就是这种情况。若Bi≫1，则固体内部出现大的温度梯度，固体的热传导为速率限制——厚陶瓷壁的炉子就是这种情况。有趣的是「CHT真正必要的是Bi≒1的场合」，两者相抗衡时，同时求解固体和流体才有意义。在Bi≪1或Bi≫1的情况下，仅详细分析各侧并给定边界条件的「解耦方法」往往也能获得足够精度。
共役热传递的数值计算方法
数值方法详解
具体怎样用算法来求解共役热传递分析呢？
离散化的定式化
用形状函数 $N_i$ 来近似未知量：
用公式表示就是这样。
基本方程的离散形式
用公式表示就是这样。
嗯……光看式子我感觉不到…… 它表示什么呢？
连续体的支配方程离散化后，得到以下代数方程组：
这里 $[K]$ 是整体刚度矩阵（或等价的系统矩阵），$\{u\}$ 是未知节点变量向量，$\{F\}$ 是外力向量。
啊，明白了！连续体的支配方程就是这样离散化的啊。
单元技术
我听过「单元技术」这个词，但可能没真正理解……
| 单元类型 | 次数 | 节点数(3D) | 精度 | 计算成本 |
|---|---|---|---|---|
| 四面体1次 | 线性 | 4 | 低（剪切锁定） | 低 |
| 四面体2次 | 二次 | 10 | 高 | 中 |
| 六面体1次 | 线性 | 8 | 中 | 中 |
| 六面体2次 | 二次 | 20 | 非常高 | 高 |
| 棱柱 | 线性/二次 | 6/15 | 中～高 | 中 |
积分方案
积分方案，具体来说是什么意思呢？
听到这里，才终于理解了为什么单元类型很重要！
收敛性和稳定性
收敛失败的话，首先应该检查什么呢？
收敛速度：二次单元中误差以 $O(h^2)$ 的阶减小（光滑解的情况）
明白了……网格细分看似简单，但实际上非常深奥啊。
求解器设置建议
具体怎样用算法来求解共役热传递分析呢？
| 参数 | 推荐值 | 备注 |
|---|---|---|
| 迭代法收敛判定 | $10^{-6}$ | 残差范数基准 |
| 预处理方法 | ILU(0) 或 AMG | 根据问题规模 |
| 最大迭代次数 | 1000 | 未收敛时重新审视设置 |
| 内存模式 | In-core | 尽可能 |
单一方法
将全部物理场作为1个联立方程组同时求解。对强耦合稳定，但实现复杂且内存消耗大。
分割法（分离迭代法）
独立求解各物理场，在界面处进行数据交换。实现容易，能利用现有求解器。适用于弱耦合。
界面数据转写
最近邻法（最简单但精度低）、投影法（保守）、RBF插值（对网格非一致性强）。保守性与精度的平衡很重要。
子迭代
在各耦合步内进行充分迭代，确保界面条件的整合性。残差基准应根据各物理场的典型值进行缩放。
Aitken缓和
自动调整耦合迭代的缓和系数。防止过缓和导致的发散，同时加速收敛的自适应方法。
稳定性条件
注意added mass效应（流体-结构耦合中结构密度≈流体密度的场合）。不稳定时应用Robin型界面条件或IQN-ILS法。
共役热传递的实务应用
实践指南
老师，请给我讲讲「实践指南」！
解说共役热传递分析的实务分析流程和注意事项。
老师的说明好懂！共役热传递分析实务的疑虑一扫而光。
分析流程
从最初一步开始教我！应该从什么开始？
1. 前处理（Pre-processing）
- CAD数据的导入及形状简化
- 材料特性的定义
- 网格生成（单元类型尺寸的确定）
- 边界条件及荷载条件的设置
2. 求解（Solving）
- 求解器设置（求解法、收敛基准、输出控制）
- 任务投入及计算执行
- 收敛监控
3. 后处理（Post-processing）
- 结果的可视化（位移、应力、其他物理量）
- 结果的检证及妥当性确认
- 报告制作
网格生成的最佳实践
怎样判断网格的好坏呢？
单元品质指标
请给我讲讲「单元品质指标」！
| 指标 | 理想值 | 允许范围 | 影响 |
|---|---|---|---|
| 纵横比 | 1.0 | < 5.0 | 精度下降 |
| 雅可比比率 | 1.0 | > 0.3 | 单元退化 |
| 翘曲 | 0° | < 15° | 精度下降 |
| 歪斜度 | 0° | < 45° | 收敛性恶化 |
| 锥形比 | 0 | < 0.5 | 精度下降 |
网格密度的确定
网格密度的确定，具体来说是什么意思呢？
边界条件设置指南
边界条件，听说如果这里出错的话全都完蛋……
啊，明白了！过度拘束注意就是这样的仕组啊。
按商用工具的实现步骤
有各种软件吧？请给我讲讲各个的特点！
| 工具名 | 开发方/现在 | 主要文件格式 |
|---|---|---|
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
| COMSOL Multiphysics | COMSOL AB | .mph |
| ANSYS Mechanical（原ANSYS Structural） | Ansys Inc. | .cdb, .rst, .db, .ans, .mac |
Ansys Fluent
接下来是Ansys Fluent的话题。内容是什么呢？
由Fluent Inc.开发。2006年被Ansys收购。基于非结构化网格的通用CFD求解器。
现在的所属：Ansys Inc.
Simcenter STAR-CCM+
接下来是Simcenter STAR的话题。内容是什么呢？
由CD-adapco开发。2016年被Siemens收购并整合到Simcenter品牌下。多面体网格是其特点。
现在的所属：Siemens Digital Industries Software
先生的说明好懂！工具名的疑虑一扫而光。
常见失败和对策
初学者容易犯的失败模式有吗？事先想了解！
| 症状 | 原因 | 对策 |
|---|---|---|
| 计算不收敛 | 网格品质不良、不适当的边界条件 | 网格改善、拘束条件重新检视 |
| 应力异常大 | 应力奇异点、网格依赖 | 奇异点回避、局部网格细分化 |
| 位移不现实 | 材料常数错误、单位系混乱 | 输入数据确认 |
| 计算时间过长 | 不必要的细分化、无效率的求解法 | 网格最优化、并行计算 |
质量保证检查清单
教科书里没有的「现场智慧」这样的东西有吗？
哎呀，共役热传递分析真的很深奥啊……但是多亏了老师的说明，我理解得差不多了！
嗯，加油啊！实际动手做是最好的学习方法。有不明白的地方随时问我。
燃气轮机叶片冷却孔——CHT改变了设计过程
现代喷气发动机的高压涡轮叶片在燃烧气体温度（1400～1600℃）超过叶片材料融点（镍超合金：约1300℃）的状态下使用。这由叶片内部复杂的冷却空气流路实现，从叶片表面打开的数百个小孔喷出冷气形成膜冷却。这种设计中固体（叶片）和流体（冷却空气燃烧气）的耦合=共役热传递分析是必需的，GE和Rolls-Royce从1970年代开始用内部代码实施。现在用ANSYS Fluent的尺度分析模拟（SRS）可以分析冷却孔周围的喷流干涉了。
共役热传递的软件比较
商用工具比较
有各种软件吧？请给我讲讲各个的特点！
解说支持共役热传递分析的主要商用CAE工具的功能比较，及各产品的历史背景。
等等，支持共役热传递分析，也就是说可以在这样的情况下也使用吗？
支持工具列表
那么，进行共役热传递分析的话用什么软件呢？
| 工具名 | 开发方/现在 | 主要文件格式 |
|---|---|---|
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
| COMSOL Multiphysics | COMSOL AB | .mph |
| ANSYS Mechanical（原ANSYS Structural） | Ansys Inc. | .cdb, .rst, .db, .ans, .mac |
Ansys Fluent
接下来是Ansys Fluent的话题。内容是什么呢？
由Fluent Inc.开发。2006年被Ansys收购。基于非结构化网格的通用CFD求解器。
现在的所属：Ansys Inc.
Simcenter STAR-CCM+
接下来是Simcenter STAR的话题。内容是什么呢？
由CD-adapco开发。2016年被Siemens收购并整合到Simcenter品牌下。多面体网格是其特点。
现在的所属：Siemens Digital Industries Software
听到这里，才终于理解了为什么开发历史很重要！
COMSOL Multiphysics
请给我讲讲「COMSOL Multiphysics」！
1986年在瑞典成立。最初作为MATLAB联动的FEMLAB开始，后来改名为COMSOL。在多物理场方面有强项。
现在的所属：COMSOL AB
ANSYS Mechanical（原ANSYS Structural）
请给我讲讲「ANSYS Mechanical」！
1970年由Swanson Analysis Systems Inc. (SASI) 开发。基于APDL（Ansys参数化设计语言）。
现在的所属：Ansys Inc.
啊，明白了！开发就是这样的仕组啊。
功能比较矩阵
预算和时间都有限，最划算的是哪个？
| 功能 | Fluent | Star-CCM+ | COMSOL | Ansys Mechanical |
|---|---|---|---|---|
| 基本功能 | ○ | ○ | ○ | ○ |
| 高级功能 | ○ | ○ | ○ | △ |
| 自动化/脚本 | ○ | ○ | ○ | ○ |
| 并行计算 | ○ | ○ | ○ | ○ |
| GPU对应 | △ | △ | △ | ○ |
转换时的风险
转换时的风险，具体来说是什么意思呢？
啊，明白了！不同工具间的模型转换就是这样的仕组啊。
许可证形式
我听过「许可证形式」这个词，但可能没真正理解……
| 工具 | 许可证 | 特点 |
|---|---|---|
| 商用FEA | 节点锁定/浮动 | 高额但有官方支持 |
| OpenFOAM | GPL | 无偿但支持有偿 |
| COMSOL | 节点锁定/浮动 | 按模块购买 |
| Code_Aster | GPL | EDF开发的OSS求解器 |
选择指南
最后应该选哪个，判断标准告诉我？
共役热传递分析工具选择时需考虑以下几点：
哎呀，共役热传递分析真的很深奥啊……但是多亏了老师的说明，我理解得差不多了！
嗯，加油啊！实际动手做是最好的学习方法。有不明白的地方随时问我。
Siemens Simcenter和ANSYS——CHT有什么不同
共役热传递分析的代表工具Simcenter STAR-CCM+和ANSYS Fluent都具有高功能，但擅长领域略有不同。STAR-CCM+通过「多面体网格」对复杂形状的适应力和从CAD到分析的工作流程自动化能力很强——正因为此在汽车冷却系统设计中广泛应用。Fluent的物性数据库和化学反应模型丰富，通过UDF的扩展性强——在燃气轮机和化学工艺领域的采用率高。明确应用目的选择很重要，「先用知名的Fluent」这样的选择方式可能之后会后悔。
共役热传递的前沿研究
前沿话题和研究动向
共役热传递分析领域将来会怎样发展呢？
看一下共役热传递分析领域的最新研究动向和先进方法。
最新的数值方法
接下来是最新数值方法的话题。内容是什么呢？
嗯……光看式子我感觉不到…… 它表示什么呢？
高性能计算（HPC）的对应
| 并行化方法 | 概要 | 适用求解器 |
|---|---|---|
| MPI（领域分割） | 分布式内存型。大规模问题的标准 | 全主要求解器 |
| OpenMP | 共有内存型。节点内并行 | 多数求解器 |
| GPU（CUDA/OpenCL） | GPGPU活用。特别是显式法有效 | LS-DYNA、Fluent等 |
| 混合 MPI+OpenMP | 节点间+节点内并行 | 大规模HPC环境 |
共役热传递的故障排除
故障排除
明白了。那么共役热传递分析相关的处理做好的话，首先就没问题了吧？
常见错误和对策
老师也在共役热传递分析中通过调试熬过夜吗？（笑）
1. 收敛失败
收敛失败，具体来说是什么意思呢？