。 图1为某时刻的温度分布，图2为某时...)
comsol仿真，方向为热电发电器/温差发电器(TEG)。 图1为某时刻的温度分布，图2为某时刻的电势分布，图3为输出功率随时间的变化。
智能手表背面的硬币大小装置持续供电的秘密，藏在塞贝克效应与有限元分析的碰撞中。打开COMSOL的瞬态热电耦合模块，我们先在几何建模里画了个三明治结构——两片陶瓷基板夹着碲化铋半导体，金属电极像夹心饼干的奶油层般精准定位。
材料库调取参数时特别要注意各向异性设置：
material = "Bi2Te3"; sigma_n = 1.2e5; % n型电导率[S/m] sigma_p = 0.8e5; % p型电导率 Seebeck_n = -200e-6; % 塞贝克系数[V/K] Seebeck_p = 220e-6;这对搞反了正负号会导致电势场倒挂，就像把电池正负极接反一样尴尬。边界条件设置阶段，热端用了个随时间波动的温度函数：
temp_hot = 310 + 5*sin(t/10); // 310K基准叠加周期波动这模拟了实际工况中热源不稳定的状态，比如汽车尾气管的温度脉动。求解器跑起来后，图1的温度云图显示出明显的阶梯状分布，半导体臂两侧的温差达到67K，热流在p-n结连接处形成漩涡状分布。
电势分布（图2）呈现对称的太极图样，p型区的红色电势云与n型区的蓝色区域在电极处交汇。通过后处理提取的功率曲线（图3）暴露了有趣现象——虽然输出电压跟随温度变化，但功率峰值却有0.5秒的滞后，这其实是材料热容导致的缓冲效应：
// 瞬态计算核心方程 dT/dt = (∇·(k∇T) - J·E)/ρCp; ∇·(σ∇V) = -∇·(σS∇T);第一个方程里的ρCp项就像热水袋储热，延缓了温度传导。当把半导体臂截面从2x2mm改为3x3mm时，电流提升了2.25倍但内阻也增大，这中间的平衡点需要参数化扫描来找，就像在迷宫里找最短逃生路线。
有个坑新人容易踩：网格划分时在电极-半导体界面用了过粗的网格，导致接触电阻计算偏差20%。后来改用边界层网格，像给手机贴膜般在关键界面布置了五层渐进式网格，结果立刻顺滑多了。这种细节处理能力，往往是区分老司机和小白的关键指标。
