2026《温差发电变物性模型分析综述》4400字.docxVIP

PAGE40

温差发电变物性模型分析综述

目录

TOC\o1-3\h\u28735温差发电变物性模型分析综述 1

71291.1求解变物性模型电偶臂冷热端温度 1

319421.2求解变物性模型发电性能参数 3

319351.3常物性模型和变物性模型的对比分析 5

1.1求解变物性模型电偶臂冷热端温度

在分析变物性模型的过程中，采用的是有限元分析方法即将电偶臂平均分成许多小单元，在Matlab中使用fsolve函数编程求解得到电偶臂冷热端温度的具体数值，进而研究发电性能。考虑陶瓷基板、导流片以及相互间的接触效应，可以得到如下热导分布示意图，如图1.1所示。

图1.1热导分布示意图

图1.1中Ta、Tb分别表示的是电偶臂热端与冷端温度，K1为陶瓷基板的热导，K2为导流铜片的热导，K3、K4分别为热端和冷端的接触热导，为了简化计算，对变物性模型先作出下列假设[35]：

1.陶瓷基板外表面几乎无散热即温度恒定并且温度分布均匀；

2.P、N型材料分布均匀，热量沿长度方向传递，且二者电偶臂结构相同；

3.为了简化计算，假设陶瓷基板和导流片的热导率不发生变化，并且将热阻的数值设置相同，考虑二者接触效应的影响，忽略导流铜片和电偶臂间的接触效应；

1.稳态运行情况下，不考虑汤姆逊效应；

假设陶瓷基板冷热端的温度分别为Th和Tc，单位时间内，电偶臂热端得到的热量为Qh,传递到冷端时的热量为Qc,则：

Qh=Kh

Qc=Kc

由热导串并联公式可以得到温差发电单元冷、热端总热导表达式分别为：

Kc=K

Kh=K

假设焦耳效应产生的热量传递到电偶臂热端的系数是f,传递到电偶臂冷端的系数是1?f

f=KhK

可得如下方程：

Qh=α

Qc=α

进一步可得如下方程组：

K1Th

由式子（1.8）可以知道，电偶臂的数学模型是由非线性方程构造起来的，可以使用最小二乘法编程求解其冷热端温度，为后文进一步研究变物性模型的输出性能打下基础。

1.2求解变物性模型发电性能参数

1.根据前文赛贝克效应的定义可得温差为?T时输出的电动势为：

E=TbT

上式中αPN

假设将电偶臂的长度L平均划分为了n-1层，每层电阻的大小为ri，则具体表达式为：

ri=ΔLρ

倘若不考虑导流片的内阻，那么电偶臂的总电阻可以近似为所有被划分的电阻之和，即：

Rleg=i=1

那么每次迭代循环的热电偶总电阻和回路电流分别为：

R=i=1n?RN?leg+R

2.温差发电变物性模型下的输出功率为：

Pout=E

式中E是温度和热电材料赛贝克系数的函数，并且根据表达式（1.10）可以知道热电偶内阻是长度、电阻率以及电偶臂截面面积的函数，所以式子（1.14）可以进一步写成如下函数形式：

Pout=fT

3.温差发电变物性模型下的热电转换效率为：

η=Pout?

式中：

η：温差发电单元变物性模型的热电转换效率，没有量纲；

Pout

?：温差发电单元变物性模型热端吸收的热流量，单位：W；

从第三章常物性模型的分析中，可以知道电偶臂热端的热流量由珀尔帖热、传导热以及焦耳热构成。在对变物性模型下热流量传递分析的过程中，假定产生的焦耳热只有一半传递到热端，由于沿热电偶臂方向上的温度是呈非线性分布的并且导热系数的数值也会随温度的变化而改变，所以变物性模型下热端吸收的热量为：

Qh=IT

1.3常物性模型和变物性模型的对比分析

利用前文所建立的GUI界面计算常物性模型下温差发电单元的输出电压、电流、功率以及转换效率，再根据1.1和1.2小节计算出变物性模型下这些参数的具体数值，最后将这些输出值与ANSYS中的值进行对比，说明这两种模型的准确性和适用性。

研究变物性模型的输出性能之前，首先需要选定好所用的热电材料和温