2026《温差发电常物性模型研究》3700字.docxVIP

温差发电常物性模型研究

目录

TOC\o1-3\h\u23621温差发电常物性模型研究 1

108091.1常物性模型温差发电基本理论 1

154591.2常物性模型输出功率的计算 3

60211.3常物性模型GUI图形界面 5

491.4ANSYS下输出性能分析 11

实际应用中的温差发电装置是由许多温差发电单元通过特定的连接方式所构成的，因此如果要研究温差发电装置，首先需要对温差发电单元有深入的了解，那本章就主要借助仿真软件MATLAB和热电耦合仿真软件ANSYS对温差发电的常物性模型展开研究。

1.1常物性模型温差发电基本理论

实际的温差发电过程较为复杂，需要考虑的因素较多，为了简化分析得出一些一般性的结论，那需要在搭建常物性模型之前作出下列假设[28]：

(1)近似认为电偶臂上下端面的温度和所给定的冷、热源的温度相同，一旦温度给定后，端面的温度保持不变而且分布均匀。

(2)P、N型材料组成的电偶臂尺寸相同，并且它们的热电参数：赛贝克系数、导热系数、电阻率是恒定值。

(3)不考虑各接触面间的接触效应，同时因为所研究的温度范围不是很大，汤姆逊效应较为微弱，所以也可以忽略不计。

在温差发电的过程中，温度场和电场是相伴而生的，可以用如下方程组描述二者关系[29]：

E=??μ

j=1ρE?

q=πj

w=q+

式中各参数大都与温度和位移有联系，具体含义如下：

E：电场强度，单位：V/m；

μ：费米能级，单位：eV；

e：电子电荷量，单位：C；

φ：电场电位，单位：V；

J:电流密度，单位：A/m2;

ρ:电阻率，单位：Ω·m；

α：赛贝克系数，单位：V/K；

?T：T的散度记法；

q：热流密度，单位：J/(m2·s)

π：珀尔帖系数，单位：W/A；

λ：热导率，单位：W/(m·K)；

w：总能流密度，单位：J/(m2·s)；

通过求解式子（3-1）~（3-4）可以较为精确的得到稳定工作条件下温度和能流的分布变化情况，前提是需要满足以下条件[30]：

{div?(j

式子(1.5)中，w=q+μ?e

上式中相关参数的含义如下：

k(x,T)：温差发电单元热电偶沿长度方向的热导率，随着温度的变化而发生改变；

ρ(x,T)：温差发电单元热电偶沿长度方向的电阻率，随着温度的变化而发生改变；

α(x,T)：温差发电单元热电偶的赛贝克系数，其大小与温度和位移有关；

IS

所以总的来看，式子(1.6)表示了热传导效应、焦耳热效应和赛贝克效应组合后的变化情况，现在有许多研究人员都对此展开了研究，由于赛贝克系数、热导率和电阻率都是随温度和位置发生变化的量，因而求解并不容易。

1.2常物性模型输出功率的计算

温差发电单元热源温度设定以后，热量会在电偶臂间流通，如果不加以控制，热端温度会上升而冷端温度会下降，最终会影响到温差发电装置的输出性能。为了让温差发电装置稳定正常工作，除了要在陶瓷基板两侧设置相应的集热、散热装置外，还要连续不断地在热端提供热量，使热电偶两端温差保持恒定值，稳定的对外输出恒定功率。

将热电材料的赛贝克系数、电阻率、热导率视为恒定值，并用平均温度处的取值来代替，可以得到温差发电单元温差为ΔT时，常物性模型下的输出电动势：

E=αPN—

式中参数的含义如下：

αP—、αN—

Th、T

考虑温差发电单元的内阻，则流过的电流为：

I=ER+R

其中R、RL分别为温差发电单元的内阻和负载电阻，单位为Ω，假设令m=RL/R，则由式子(3-8）也可以写成如下式子：

I=αPN—

可以推导出温差发电单元输出功率为：

P=I2?

进一步推导出若m=1时，则输出功率最大，为：

Pmax=E

由此可以算出单位时间内温差发电单元冷、热端的珀尔帖热分别为：

Qpc=π

Qph=π

由于温度梯度而产生的传导热为：

QT=KT

常物性模型下假设温差发电单元四周绝热性能良好，所以可以假定焦耳热传送到电偶臂冷、热端，由此可以得到冷热端总的热流表达式分别为：

Qc=α

Qh=α

由式子（1.15）、（1.16）可以知道，不论是电偶臂冷端还是热端的总热流，它们都是由电