自然对流-传导耦合
自然对流-传导耦合的理论基础
概述
老师！ 今天讲自然对流-传导耦合是吧？那是什么东西？
浮力驱动流与固体传导的耦合。密闭空间的自然对流。建筑保温设计。Rayleigh数和Nusselt数。
哦~浮力驱动流与固体传导的话，真的超有意思！ 请继续给我讲讲。
控制方程
离散化方法
这些方程在计算机上如何实际求解呢？
采用有限元法（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 | 可视化 | Visualization Toolkit格式。用于ParaView等。 |
在不同求解器间转换模型时，需要注意单元类型的对应关系、材料模型的兼容性、荷载和边界条件的表达差异。特别是高次单元或特殊单元（内聚单元、用户定义单元等）往往无法在求解器间直接转换。
明白了…格式看起来很简单，但其实深不见底呢。
实务注意事项
教科书里没讲的「现场智慧」有什么吗？
网格收敛性验证、边界条件合理性检验、材料参数的敏感性分析真的很重要。
自然对流-传导耦合的全貌我掌握了！ 明天开始就在实务中用这些意识试试。
好！你的状态很不错！ 实际动手操作才是最好的学习，有问题随时来问我。
Rayleigh-Bénard对流——「混沌」的入口
水平面从下方加热时产生的Rayleigh-Bénard对流，是热流动向混沌转变的著名入口。当Rayleigh数（Ra=GrPr）小于1708时只有传导，超过1708时出现规则的柱状对流，进一步增大时出现周期振荡，最终进入乱流混沌。这一「从有序到混乱的转变」是Lorenz吸引子研究的起源，也是气候科学中研究大气和海洋混沌行为的原始模型。现代CFD通过DNS可以解析Ra数超过10⁸的Rayleigh-Bénard乱流，但大气的Ra数（约10²⁰）即使用超级计算机也难以完全解析。
自然对流-传导耦合的数值计算方法
数值方法详解
具体用什么算法来求解自然对流-传导耦合呢？
离散化公式
使用形状函数 $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 | 尽可能使用 |
单一法
将所有物理场作为一个联立方程组同时求解。对强耦合问题稳定，但实现复杂、内存占用大。
分割法（分离迭代法）
各物理场独立求解，在界面处进行数据交换。实现容易，可复用现有求解器。适用于弱耦合。
界面数据转移
最邻近法（最简单但精度低）、射影法（保守性好）、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
听你的解释分明了！ 工具的迷茫都消散了。
常见失败与对策
初学者常犯的错误有什么规律吗？ 想事先知道！
| 症状 | 原因 | 对策 |
|---|---|---|
| 计算不收敛 | 网格品质不良、边界条件不当 | 改进网格、重审约束条件 |
| 应力异常大 | 应力奇点、网格相关性 | 避免奇点、局部网格细化 |
| 位移非物理 | 材料常数错误、单位系混乱 | 确认输入数据 |
| 计算时间过长 | 不必要的细化、求解低效 | 网格最优化、并行计算 |
质量保证检查清单
教科书里没讲的「现场智慧」有什么吗？
自然对流-传导耦合的全貌我掌握了！ 明天开始就在实务中用这些意识试试。
好！你的状态很不错！ 实际动手操作才是最好的学习，有问题随时来问我。
电力变压器的自然冷却——油浸变压器的热设计
支撑电力网的油浸变压器（最大容量1000 MVA以上）由线圈铜损和铁心铁损产生热，通过绝缘油的自然对流进行冷却。油的粘度温度依赖性很强（40℃→80℃约降为1/4），线圈周围的高温油上升，外部冷却翅片冷却的低温油下降，形成恒定的热虹吸式自然对流。变压器设计中「最高点油温」（TOT：Top Oil Temperature）和「最高点卷线温度」（HST：Hot Spot Temperature）由IEC 60076规范规定，用自然对流-传导耦合CFD预测这些值已成为设计品质保证的标准工程。
自然对流-传导耦合的软件比较
商用工具比较
有各种各样的软件吧？ 各自的特色给我讲讲！
自然对流-传导耦合支持的主要商用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 Parametric Design Language）。
现属：Ansys Inc.
啊！原来如此！开发就是这样的机制啊！
功能比较矩阵
预算和时间都有限，性价比最高的是哪个呢？
| 功能 | Fluent | Star-CCM+ | COMSOL | Ansys Mechanical |
|---|---|---|---|---|
| 基础功能 | ○ | ○ | ○ | ○ |
| 高级功能 | ○ | ○ | ○ | △ |
| 自动化/脚本 | ○ | ○ | ○ | ○ |
| 并行计算 | ○ | ○ | ○ | ○ |
| GPU支持 | △ | △ | △ | ○ |
转换时的风险
转换时的风险具体是什么意思？
啊！原来如此！不同工具间的转换就是这样的机制啊！
许可证形式
听过「许可证形式」这个词，但我可能理解得不够透彻…
| 工具 | 许可证 | 特色 |
|---|---|---|
| 商用FEA | 节点锁定/浮动 | 高价但有官方支持 |
| OpenFOAM | GPL | 免费但支持付费 |
| COMSOL | 节点锁定/浮动 | 按模块购买 |
| Code_Aster | GPL | EDF开发的开源求解器 |
选择指南
预算和时间都有限，到底选哪个呢？ 判断标准是什么？
自然对流-传导耦合工具选择考虑以下因素：
自然对流-传导耦合的全貌我掌握了！ 明天开始就在实务中用这些意识试试。
好！你的状态很不错！ 实际动手操作才是最好的学习，有问题随时来问我。
FloTHERM和Icepak——电子设备自然对流冷却用什么工具比较合适
专门针对电子设备自然对流冷却解析的有名工具是Mentor公司的FloTHERM（现Siemens旗下）和ANSYS Icepak。两者都支持PCB·电子部件的发热设置的便利性，以及网络模型和紧凑型热模型（CTM）。不同之处在于，FloTHERM的「SmartParts」部件库很充实，电子设备设计者容易上手，而Icepak则与ANSYS生态系统的联动很强（Mechanical进行热应力解析等）。两者的GUI都很丰富，与OpenFOAM的自然对流求解器（buoyantSimpleFoam）相比，非CFD专家也能在设计阶段使用，获得好评。
自然对流-传导耦合的先端研究
先端话题与研究动向
自然对流-传导耦合领域今后会怎么发展？
自然对流-传导耦合的最新研究动向和先进手法探索。
最新的数值方法
下面是最新数值方法的介绍是吧？内容如何？
嗯…仅看公式我还是不太明白… 这些表示什么呢？
高性能计算 (HPC) 的支持
| 并行化方法 | 概述 | 应用求解器 |
|---|---|---|
| MPI (领域分割) | 分布式内存型。大规模问题的标准 | 全主要求解器 |
| OpenMP | 共享内存型。节点内并行 | 多数求解器 |
| GPU (CUDA/OpenCL) | GPGPU活用。特别是显式法有效 | LS-DYNA, Fluent等 |
| 混合 MPI+OpenMP | 节点间+节点内并行 | 大规模HPC环境 |
自然对流-传导耦合的故障处理
故障处理
常见错误与对策
老师也有过自然对流-传导耦合的通宵调试吗？（笑）
1. 收敛失败
收敛失败具体是什么意思？
症状：求解器未在指定迭代次数内收敛，异常终止
可能原因：
- 网格品质不足（过度歪曲的单元）
- 材料参数设置不当
- 初始条件不当
- 非线性过强（荷载步数不足）
对策：
- 执行网格品质检查（宽高比、雅可比行列式）
- 确认材料参数单位系
- 分步加荷（增加子步数）