热电冷却（珀尔帖效应）
热电冷却（珀尔帖效应）的理论基础
概述
老师！今天是热电冷却（珀尔帖元件）的话题对吧？这是什么样的东西？
利用塞贝克效应·珀尔帖效应的电子冷却。传感器冷却、小型恒温槽适用。
支配方程
等等等等，热电冷却的话，也可以用在这样的情况吗？
离散化手法
这个方程怎样在计算机上实际求解？
采用有限元法（FEM）进行空间离散化。组装要素刚性矩阵，构建全局刚性方程。
矩阵求解算法
矩阵求解算法具体是什么意思？
通过直接法（LU分解、Cholesky分解）或迭代法（CG法、GMRES法）求解连立方程。大规模问题中预处理迭代法很有效。
| 求解法 | 分类 | 内存使用量 | 适用规模 |
|---|---|---|---|
| LU分解 | 直接法 | O(n²) | 小～中规模 |
| Cholesky分解 | 直接法（对称正定） | O(n²) | 小～中规模 |
| PCG法 | 迭代法 | O(n) | 大规模 |
| GMRES法 | 迭代法 | O(n·m) | 大规模·非对称 |
| AMG预处理 | 预处理 | O(n) | 超大规模 |
也就是说，在有限元法那一步马虎的话，之后会很吃亏对吧。我记住了！
商用工具中的实现
那，做热电冷却（珀尔帖元件）需要什么样的软件？
| 工具名 | 开发商/现状 | 主要文件格式 |
|---|---|---|
| Ansys Mechanical (旧ANSYS Structural) | Ansys Inc. | .cdb, .rst, .db, .ans, .mac |
| COMSOL Multiphysics | COMSOL AB | .mph |
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
供应商系统和产品整合经历
各个软件的成长过程，是不是有些戏剧性的故事？
Ansys Mechanical (旧ANSYS Structural)
请给我讲讲「Ansys Mechanical」！
1970年由Swanson Analysis Systems Inc. (SASI)开发。基于APDL（Ansys参数化设计语言）。
现属：Ansys Inc.
COMSOL Multiphysics
请给我讲讲「COMSOL Multiphysics」！
1986年在瑞典成立。始于MATLAB联动的FEMLAB，后改名为COMSOL。多物理场方面强。
现属：COMSOL AB
听到这里，终于理解为什么要关注开发历程了！
Ansys Fluent
接下来讲ANSYS Fluent的话题对吧。内容是什么？
由Fluent Inc.开发。2006年被Ansys收购。基于非结构网格的通用CFD求解器。
现属：Ansys Inc.
哦～，开发的故事，真的超级有意思！让我多听一些。
文件格式和相互可操作性
不同软件间交换数据时有什么要注意的？
| 格式 | 扩展名 | 种类 | 概述 |
|---|---|---|---|
| STEP | .stp/.step | 中立CAD | 符合ISO 10303的3D CAD数据交换格式。支持形状+PMI。 |
| IGES | .igs/.iges | 中立CAD | 早期CAD数据交换标准。曲面数据的兼容性有课题。向STEP过渡。 |
| VTK | .vtk/.vtu | 可视化 | 可视化工具箱格式。用于ParaView等。 |
在不同求解器间转换模型时，需要注意要素类型的对应关系、材料模型的兼容性、荷载·边界条件的表现差异。特别是高次要素或特殊要素（内聚要素、用户定义要素等）在求解器间往往无法直接转换。
明白了……格式看起来很简单，但实际上有很深的内涵啊。
实务注意事项
教科书里没有的「现场智慧」有吗？
网格收敛性确认、边界条件的妥当性验证、材料参数的灵敏度分析特别重要。
好的，很有劲！实际动手练习是最好的学习。有不懂的地方随时问我。
珀尔帖冷却的最大温度差是70℃
珀尔帖元件（TEC模块）通过通电在冷却侧（冷端）和发热侧（热端）产生温度差。Bi₂Te₃系TEC的最大温度差ΔTmax ≈ ZT²/2 = 约65～75℃（热端25℃固定时）是物理上限，超越这一冷却通过多段堆积（级联）结构实现。二段级联时ΔTmax ≈ 100～110℃，三段时130℃以上的温度差成为可能，已应用于光通信激光器恒温槽（-40～+85℃工作保证）。
热电冷却（珀尔帖效应）的数值计算手法
数值手法详细说明
具体用什么算法来求解热电冷却（珀尔帖元件）？
离散化的定式化
用形状函数$N_i$近似未知量：
用公式表示如下。
基础方程的离散形式
用公式表示如下。
嗯，只有公式的话摸不着头脑……这表示什么？
连续体的支配方程离散化后得到下列代数方程组：
这里[K]是全局刚性矩阵（或等效的系统矩阵），{u}是未知节点变量矢量，{F}是外力矢量。
哦，是这样！连续体的支配方程离散化的仕组就是这样啊。
要素技术
听说过「要素技术」，但可能没正确理解……
| 要素类型 | 次数 | 节点数(3D) | 精度 | 计算成本 |
|---|---|---|---|---|
| 四面体1次 | 线性 | 4 | 低（剪切locking） | 低 |
| 四面体2次 | 二次 | 10 | 高 | 中 |
| 六面体1次 | 线性 | 8 | 中 | 中 |
| 六面体2次 | 二次 | 20 | 非常高 | 高 |
| 棱柱 | 线性/二次 | 6/15 | 中～高 | 中 |
积分方案
积分方案具体是什么意思？
听到这里，终于明白为什么要素类型这么重要了！
收敛性和稳定性
收敛不了时，首先检查什么？
收敛速度：二次要素时误差以$O(h^2)$阶减小（光滑解的情况）
明白了……网格细化看似简单，实际非常深刻啊。
求解器设置建议
具体用什么算法来求解热电冷却（珀尔帖元件）？
| 参数 | 推荐值 | 备注 |
|---|---|---|
| 迭代法收敛判定 | $10^{-6}$ | 残差范数准则 |
| 预处理手法 | ILU(0) or AMG | 取决于问题规模 |
| 最大迭代次数 | 1000 | 未收敛时重新检查设置 |
| 内存模式 | In-core | 尽可能使用 |
线性要素 vs 2次要素
热传导分析中线性要素通常能获得足够精度。温度梯度急变区域（热冲击等）推荐2次要素。
热流量评估
从要素内温度梯度计算。与节点应力相同，往往需要光滑处理。
对流-扩散问题
佩克莱数高（对流支配）时需要逆风稳定化（SUPG等）。纯热传导问题无需。
非定常分析的时间步长
热扩散的特征时间$\tau = L^2 / \alpha$（$\alpha$：热扩散率）中设定充分小的步长。急剧温度变化时自动时间步长控制有效。
非线性收敛
温度依赖物性导致的非线性通常较缓和，Picard迭代（直接置换法）往往充分。放射的强非线性推荐牛顿法。
定常分析的判定
所有节点温度变化不超过阈值（$|\Delta T| / T_{max} < 10^{-5}$等）判定收敛。
热电冷却（珀尔帖效应）的实务应用
实践指南
老师，请给我讲讲「实践指南」！
热电冷却（珀尔帖元件）的实务分析流程和注意点进行讲解。
分析流程
从最开始教我！从什么开始？
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 Mechanical (旧ANSYS Structural) | Ansys Inc. | .cdb, .rst, .db, .ans, .mac |
| COMSOL Multiphysics | COMSOL AB | .mph |
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
Ansys Mechanical (旧ANSYS Structural)
请给我讲讲「Ansys Mechanical」！
1970年由Swanson Analysis Systems Inc. (SASI)开发。基于APDL（Ansys参数化设计语言）。
现属：Ansys Inc.
COMSOL Multiphysics
请给我讲讲「COMSOL Multiphysics」！
1986年在瑞典成立。始于MATLAB联动的FEMLAB，后改名为COMSOL。多物理场方面强。
现属：COMSOL AB
你的说明很容易理解！工具名的困惑解开了。
常见失败和对策
初心者常犯什么失误？想事先了解！
| 症状 | 原因 | 对策 |
|---|---|---|
| 计算不收敛 | 网格品质不足、不恰当的边界条件 | 网格改善、拘束条件审视 |
| 应力异常大 | 应力奇异点、网格依赖 | 奇异点回避、局部网格细分 |
| 位移不现实 | 材料常数误差、单位系不一致 | 确认输入数据 |
| 计算时间过长 | 不必要的细分、无效的求解 | 网格最优化、并行计算 |
质量保证检查表
教科书里没有的「现场智慧」有吗？
好的，很有劲！实际动手练习是最好的学习。有不懂的地方随时问我。
PCR装置的温度循环由珀尔帖担当
COVID-19诊断中使用的PCR（聚合酶链式反应）装置的温度循环（95℃→55℃→72℃重复30～40次）是由珀尔帖元件和PID控制实现的。为了满足温度变化速度±5℃/秒、精度±0.1℃的要求，美国Bio-Rad公司「C1000 Touch」采用8×24排列的TEC阵列独立控制结构。温度精度的关键是TEC直下的Pt100热电对传感器响应时间（100msec以下）和控制增益优化。
热电冷却（珀尔帖效应）的软件比较
商用工具比较
有很多软件对吧？各自的特点教一下！
热电冷却（珀尔帖元件）的对应主要商用CAE工具的功能比较和各产品的历史背景进行详述。
对应工具一览
那，做热电冷却（珀尔帖元件）需要什么样的软件？
| 工具名 | 开发商/现状 | 主要文件格式 |
|---|---|---|
| Ansys Mechanical (旧ANSYS Structural) | Ansys Inc. | .cdb, .rst, .db, .ans, .mac |
| COMSOL Multiphysics | COMSOL AB | .mph |
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
Ansys Mechanical (旧ANSYS Structural)
请给我讲讲「Ansys Mechanical」！
1970年由Swanson Analysis Systems Inc. (SASI)开发。基于APDL（Ansys参数化设计语言）。
现属：Ansys Inc.
COMSOL Multiphysics
请给我讲讲「COMSOL Multiphysics」！
1986年在瑞典成立。始于MATLAB联动的FEMLAB，后改名为COMSOL。多物理场方面强。
现属：COMSOL AB
听到这里，终于理解为什么要关注开发历程了！
Simcenter STAR-CCM+
接下来讲Simcenter STAR的话题对吧。内容是什么？
由CD-adapco开发。2016年被Siemens收购，整合到Simcenter品牌中。多面体网格是特色。
现属：Siemens Digital Industries Software
哦，是这样！开发故事就是这种仕组啊。
功能比较矩阵
预算和时间都有限，成本效益最高的是哪个？
| 功能 | Ansys Mechanical | COMSOL | Fluent | Star-CCM+ |
|---|---|---|---|---|
| 基本功能 | ◎ | ◎ | ◎ | ◎ |
| 高级功能 | ◎ | ◎ | ◎ | △ |
| 自动化/脚本 | ◎ | ◎ | ◎ | ◎ |
| 并行计算 | ◎ | ◎ | ◎ | ◎ |
| GPU对应 | △ | △ | △ | ◎ |
转换时的风险
转换时的风险具体是什么意思？
哦，是这样！不同工具间转换就是这种仕组啊。
许可形式
听说过「许可形式」，但可能没正确理解……
| 工具 | 许可 | 特点 |
|---|---|---|
| 商用FEA | 节点锁定/浮动 | 昂贵但有官方支持 |
| OpenFOAM | GPL | 免费但支持收费 |
| COMSOL | 节点锁定/浮动 | 按模块购买 |
| Code_Aster | GPL | EDF开发的开源求解器 |
选择指南
最后到底选哪个，给我判断标准！
在热电冷却（珀尔帖元件）工具选择中应考虑：
好的，很有劲！实际动手练习是最好的学习。有不懂的地方随时问我。
Laird Thermal Systems在精密TEC市场排名第一
精密温调用珀尔帖模块的全球最大制造商是英国Laird Thermal Systems（原Marlow Industries + Nextreme合并），在医疗·通信·国防领域强大。日本京瓷自2021年开始量产光通信模块用超薄TEC「TEC封装」，通过小型化（2.5mm见方）实现差别化。中国国内最大的TEC厂商广东克鲁克（Guangdong Kelk）正在扩大低成本汎用TEC市场份额。
热电冷却（珀尔帖效应）的先端研究
先端课题和研究动向
热电冷却（珀尔帖元件）的领域今后如何发展？
热电冷却（珀尔帖元件）领域的最新研究动向和先进手法进行讲解。
最新的数值手法
接下来讲最新数值手法的话题对吧。内容是什么？
只有公式的话摸不着头脑……这表示什么？
高性能计算 (HPC) 应对
| 并行化手法 | 概述 | 适用求解器 |
|---|---|---|
| MPI (领域分割) | 分布式内存型。大规模问题的标准 | 全主要求解器 |
| OpenMP | 共有内存型。节点内并行 | 多数求解器 |
| GPU (CUDA/OpenCL) | GPGPU活用。特别对显式法有效 | LS-DYNA, Fluent等 |
| 混合 MPI+OpenMP | 节点间+节点内并行 | 大规模HPC环境 |
热电冷却（珀尔帖效应）的故障排查
故障排查
常见错误和对策
老师也曾为热电冷却（珀尔帖元件）通宵debug过吗？（笑）
1. 收敛失败
收敛失败具体是什么意思？
症状：求解器在指定迭代次数内未收敛而异常终止
可能的原因：
- 网格品质不足（过度扭曲的要素）
- 材料参数设置不恰当
- 初始条件不恰当
- 非线性性过强（荷载步数不足）
对策：
- 进行网格品质检查（宽纵比、雅可比比）
- 确认材料参数单位系
- 将荷载分为多个步骤（增加子步数）
- 放松收敛判定准则（但要注意精度）
也就是说，在收敛失败的地方马虎的话，之后会很吃亏对吧。我记住了！
2. 非物理的结果
接下来讲非物理的结果的话题对吧。内容是什么？
症状：应力/位移/温度等出现物理上非现实的值
可能的原因：
- 边界条件设置错误
- 单位系混在（SI单位与工学单位混同）
- 要素类型选择不恰当
- 应力奇异点的存在
对策：
- 确认反力合计（力的平衡）
- 确认单位系的一致性
- 重新检讨要素类型的恰当性
- 奇异点除去或子模型分析