热电冷却(珀尔帖元件)
热电冷却(珀尔帖元件)的理论基础
概述
先生! 今天是热电冷却(珀尔帖元件)的讲座吧? 这是什么东西?
珀尔帖效应的主动冷却。温度控制。性能系数(COP)的优化。
先生的讲解很清楚! 关于珀尔帖效应的疑惑已解开。
支配方程
等等,热电冷却是这样的情况吗?也能用?
离散化方法
这个方程在计算机上如何实际求解?
使用有限元法(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 |
| 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 Mechanical"!
1970年由Swanson Analysis Systems Inc. (SASI) 开发。基于APDL(Ansys参数化设计语言)。
现在所属:Ansys Inc.
Ansys Fluent
接下来是Ansys Fluent的话题吧。讲讲内容?
由Fluent Inc.开发。2006年被Ansys收购。基于非结构网格的通用CFD求解器。
现在所属:Ansys Inc.
听到这里,终于明白为什么开发历程这么重要了!
Simcenter STAR-CCM+
接下来是Simcenter STAR的话题吧。讲讲内容?
由CD-adapco开发。2016年被西门子收购,整合到Simcenter品牌。多面体网格是特色。
现在所属:Siemens Digital Industries Software
哇,开发历程的故事太有趣了! 继续听下去吧。
文件格式和互操作性
在不同软件之间互相转移数据时有什么需要注意的?
| 格式 | 扩展名 | 种类 | 概述 |
|---|---|---|---|
| STEP | .stp/.step | 中立CAD | 符合ISO 10303的3D CAD数据交换格式。支持形状+PMI。 |
| IGES | .igs/.iges | 中立CAD | 早期CAD数据交换规范。曲面数据兼容性存在问题。正在逐步迁移到STEP。 |
在不同求解器间转换模型时,要注意单元类型的对应关系、材料模型的兼容性、荷载和边界条件的表达差异。特别是高阶单元或特殊单元(内聚单元、用户定义单元等),往往不能在求解器间直接转换。
原来如此……格式看似简单,其实有很多深层的东西。
实务上的注意事项
教科书里没写的"现场智慧"有吗?
网格收敛性确认、边界条件妥当性验证、材料参数的敏感性分析特别重要。
是啊,进展顺利! 只有实际动手才是最好的学习。有不明白的地方随时问我吧。
塞贝克系数超过300μV/K
热电转换是塞贝克效应(1821年发现)、珀尔帖效应(1834年发现)、汤姆逊效应的三种现象总称。决定转换效率的性能指数ZT = α²σT/κ(α:塞贝克系数、σ:电导率、κ:热导率)在室温Bi₂Te₃上长期以来上限约为1.0,但2014年MIT研究(发表于Science)中SnSe单晶达到ZT=2.6,震惊业界。SnSe等新材料的塞贝克系数最大值报告已超过300μV/K。
热电冷却(珀尔帖元件)的数值计算方法
数值方法的详细内容
具体用什么算法来求解热电冷却(珀尔帖元件)?
离散化的公式化
使用形状函数 $N_i$ 来近似未知量:
用公式表示就是这样。
基本方程式的离散形式
用公式表示就是这样。
嗯,光看公式有点想象不出来… 这到底表示什么呢?
连续体的支配方程离散化后,得到下面的代数方程组:
这里$[K]$是全局刚度矩阵(或等价的系统矩阵),$\{u\}$是未知节点变量向量,$\{F\}$是外力向量。
哎呀,原来是这样! 连续体的支配方程离散化就是这个道理啊。
单元技术
"单元技术"这个词听过,但可能理解得不够深…
| 单元类型 | 阶数 | 节点数(3D) | 精度 | 计算成本 |
|---|---|---|---|---|
| 四面体1阶 | 线性 | 4 | 低(剪切锁定) | 低 |
| 四面体2阶 | 二阶 | 10 | 高 | 中 |
| 六面体1阶 | 线性 | 8 | 中 | 中 |
| 六面体2阶 | 二阶 | 20 | 非常高 | 高 |
| 棱柱体 | 线性/二阶 | 6/15 | 中~高 | 中 |
积分方案
积分方案具体是什么意思?