双声源驱动热声制冷机特征时间的熵产分析

陈　浩，吴　锋，何秋石，蒋智杰，田一泽

武汉工程大学理学院，湖北 武汉 430074

双声源驱动热声制冷机特征时间的熵产分析

陈浩，吴锋*，何秋石，蒋智杰，田一泽

武汉工程大学理学院，湖北 武汉 430074

在有限时间热力学理论的基础上，对双声源驱动热声制冷机的回热器声场进行分析，以特征时间为特征参数来表征热声系统的换能效应，寻求特征时间与回热器熵产之间的最优关系以及其他因素对熵产率的影响.结果表明：熵产率随着特征时间的减少，到达一个极小点后再平缓增大，通过调节阻抗比，寻求到一个最佳阻抗比值所对应的最小熵产率，存在一个最优特征时间使得系统的不可逆熵产最小，此时不可逆能量损失最小，有效的提高了双声源驱动热声制冷机效率，有利于双声源驱动系统运行.

双声源驱动；热声制冷机；回热器；特征时间；熵产

1　引言

热声制冷机［1-3］是基于热声效应原理的一种新型制冷机.它具备使用寿命长、可靠性高等特点，主要使用对环境无害的气体（例如氦气、氮气等）作为工质，且无运动部件，能有效地消除机械磨损等常规制冷机所存在的问题，因具备上述优点而使其在最近几年备受关注.

热声系统性能取决于整个系统的声场.而在双声源热声系统中，可以在回热器边界处调制出任意声场，并在文献［4］中给出了声场调制方法，该方法操作简单，只需要调节两个声源就能够调制出所需任意稳定频率、稳定压力振幅和稳定相角的回热器边界声场.同时该系统也可以作为双声源驱动热声制冷机的理论模型.

熵产［5］是衡量不可逆热力学过程的指标.熵产由不可逆传热和流体的黏性流动所引起，这两部分熵产生对某些参数来说并非总是同向变化的，从而存在某个数值使系统的总熵产为最小，并意味着系统的可用能量损失最小.熵产分析可以揭示系统过程的不可逆程度，从而反映系统能量利用的效率.由于熵产分析和有限时间热力学分析的方法有着本质的联系，本文对双声源驱动热声制冷机的回热器内流动与换热过程的熵产表达式进行推导，分析双声源驱动热声制冷机特征时间及阻抗比对熵产率的影响.

2　双声源驱动热声制冷机回热器声场

在双声源驱动热声制冷机中，流相工质的压力和振荡速度的一阶波动量［5］可表示为

fk和 fν是回热器的截面平均导热分布函数和截面平均黏性分布函数，对于平行板型回热器，fk和 fν关于特征时间［7］ωτ的表达式为：

其中ωτk、ωτν分别表示热弛豫特征时间和黏性弛豫特征时间，它们是评价整个热声系统性能参数的重要指标，同时可表征整个热声系统的换能效应.热声系统的特征时间ωτ还决定了回热器的特征尺寸（侵入层厚度δ）和谐振管长度.

根据热声理论，回热器中流相工质的一阶振荡温度为：

式（1）～（3）描述了双声源热声制冷机中回热器的声场.当双声源驱动热声制冷机稳定工作时，流相工质中存在稳定的沿回热器纵向（x方向）传播的声波，若令回热器两边界在谐振管中的坐标分别为zA和zB，两边界面处声压的一阶波动量为pA和 pB，则有边界条件：p（zA）=pA，p（）zB=pB，代入式（1），得到C1和C2的表达式，分别为

3　双声源驱动热声制冷机特征时间的熵产分析

如在考虑黏性时，热声系统的能量守恒方程可写为

式（5）中，E、q、V和φ分别为流相工质的比内能、热流量、比体积和黏性耗散函数.

令流相工质与回热器壁面之间的换热符合Fourier导热定律，由可得

将式（6）代入式（5）中，得到

联立式（7）和Gibbs方程Tds=dE+pdV，并略去高阶小量，得到式（8）中，S为单位体积的熵，ρ0为流体密度为横向（y方向）梯度.

在式（8）中，左边第一项为单位体积流体的熵增率，第二项为熵流；右边第一项为黏性不可逆性所引起的熵产率［8-11］，第二项为热传导不可逆性引起的熵产率.因此单位体积流体的不可逆性熵产率可写为

在小扰动条件下，黏性耗散函数φ可写成

将式（10）代入式（9）中，并取时间平均得

将式（2）和（3）代入式（11），并取实部，得到回热器任意一位置处的时均熵产率为

而Z、y0分别表示阻抗比和板间距，“*”表示复共轭.

很明显，式（12）表明回热器的时均熵产率与y有关，即时均熵产率沿横向截面有一分布，对式（12）取截面平均，得到回热器截面分布熵产率为

4　数值算例分析

4.1特征时间对熵产率的影响

图1　ṠiR与ωτk的关系（y0=0.001 5 m）Fig.1　Effects ofωτkon the values ofṠiR（y0=0.001 5 m）

图2　ṠiR与ωτν的关系（y0=0.001 5 m）Fig.2　Effects ofωτνon the values ofṠiR（y0=0.001 5 m）

图3　ṠiR与ωτk的关系（dT0/dx=1 000 Κ/m）Fig.3　Effects ofωτkon the values ofṠiR（dT0/dx=1 000 Κ/m）

图4　ṠiR与ωτν的关系（dT0/dx=1 000 Κ/m）Fig.4　Effects ofωτνon the values ofṠiR（dT0/dx=1 000 Κ/m）

回热器截面分布熵产率ṠiR随特征时间ωτk、ωτν的变化都有这样一个过程：ṠiR随着特征时间ωτk、ωτν急剧减小，到达一个极小点以后再平缓增大.在特征时间ωτk、ωτν值较小时，熵产率ṠiR较大，即不可逆熵产较大，此时热声系统不可逆程度较大，能量损失较大，能量转换效率较低，不利于双声源热声系统的运行.在急剧减小变化以后，且变化相对平缓，即有一最优特征时间ωτk、ωτν使ṠiR为最小值，换句话说，就是存在一个最优特征时间ωτk、ωτν使得系统的不可逆熵产最小，此时不可逆能量损失最小.

4.2阻抗比对熵产率的影响

以阻抗比Z为横坐标，回热器截面分布熵产率ṠiR为纵坐标，分别取ωτk和ωτν为0.5、1.5、2.5，其他参数与4.1节相同，可以得到图5和图6.

图5　ṠiR与Z的关系（ωτk=0.5，1.5，2.5）Fig.5　Effects of Z on the values ofṠiR（ωτk=0.5，1.5，2.5）

图6　ṠiR与Z的关系（ωτν=0.5， 1.5， 2.5）Fig.6　Effects of Z on the values ofṠiR（ωτν=0.5，1.5， 2.5）

由图5和6可知，回热器截面分布熵产率ṠiR随阻抗比Z先减小后增大，同时存在一个Z值使ṠiR最小，这是导热不可逆性与粘滞不可逆性耦合的结果.通过调节阻抗比Z，可以寻求一个最佳的阻抗比值所对应的最小熵产率，使熵产率最小可提高热声系统回热器的性能，有利于双声源热声k系统的工作运行.

5　结　语

在双声源驱动热声制冷机回热器流相工质的声压和速度分布表达式的基础上，结合特征时间为特征参数，通过理论计算得到回热器中关于特征时间的声场表达式.

对双声源驱动热声制冷机回热器熵产进行分析，将特征时间与回热器熵产率的优化关系进行理论分析及数值计算.数值计算结果表明，熵产率随着特征时间急剧减小，达到一极小值点后，平缓增大，为了有利于双声源驱动系统运行，应选取熵产率较小且增大趋势平缓的时刻所对应的特征时间值，同时还考查了阻抗比对熵产率的影响.

［1］吴锋，李青，郭方中，等.热声理论的研究进展［J］.武汉工程大学学报，2012，34（1）：1-6. WU F，LI Q，GUO F Z，et al.Advance in thermoacoustic theory［J］.Journal ofWuhaninstitutetechnology，2012，34（1）：1-6.

［2］鄂青，沈秋婉，陈秋霞.热声制冷机［J］.真空与低温，2009，15（33）：117-119. E Q，SHEN Q W，CHEN Q X.Thermoacoustic refrigerator［J］.Vacuum and cryogenics，2009，15（33）：117-119.

［3］吴锋，汪拓，费锦华，等.热声网络的辛矩阵分析［J］.热科学与技术，2013，12（4）：283-289. WU F，WANG T，FEI J H，et al.Symplectic matrix analysis of thermoacoustic network［J］.Journal of thermal science and technology，2013，12（4）：283-289.

［4］李雷，谢秀娟，周立华，等.双声源法调制回热器边界声场的理论分析［J］.低温工程，2009，172（6）：30-35. LI L，XIE Y J，ZHOU L H，et al.Theoretical analysis for modulating acoustic boundary conditions of a regenerator by double acoustic drivers［J］.Cryogenics，2009，172（6）：30-35.

［5］孟繁孔，李志信.驻型波热声制冷机熵产分析［J］.工程热物理学报，2006，27（3）：508-510. MENG F K，LI Z X.Entropy generation analysis of standing wave thermoacoustic refrigerator［J］.Journal of engineering thermophysics，2006，27（3）：508-510.

［6］ 周立华，谢秀娟，李雷，等.双声源驱动热声系统声场调制研究［J］.工程热物理学报，2012，33（1）：27-31. ZHOU L H，XIE Y J，LI L，et al.Study of modulating acoustic field in the thermoacoustic system by double acoustic drivers［J］.Journal of engineering thermophysics，2012，33（1）：27-31.

［7］ 张春萍.热声核特征参数实验研究及高频微型热声实验装置的研制［D］.武汉：华中科技大学，2011.

［8］ 姚寿广.对流热质传递过程的熵产分析及优化［J］.华东船舶工业学院学报，1993，7（3）：101-107. YAO S G.Analysis and optimization of entropy production in the process of convective heat transfer［J］.Journal of hechinormal college，1993，7（3）：101-107.

［9］BEJAN A.The thermodynamic design of heat and mass transfer process and devices［J］.Heat and fluid flow，1987，8（4）：398-406.

［10］BEJAN A.Entropy generation through heat and fluid flow［M］.New york：Wiley，1982.

［11］XIE Y J，GAO G，LI Q.Analysis of entropy generation rate inside the stack of standing-wave thermoacoustic refrigerator［J］.Journal of engineering thermophysics，2012，43（9）：609-612.

本文编辑：陈小平

Entropy Production of Characteristic Time for Thermoacoustic Refrigerator by Double Acoustic Drivers

CHEN Hao，WU Feng，HE Qiushi，JANG Zhijie，TIAN Yize
School of Science，Wuhan Institute of Technology，Wuhan 430074，China

The regenerator sound field of thermoacoustic refrigerator driven by dual sound source was analyzed on the basis of the theory of finite time thermodynamics.The characteristic time was used as parameters to characterize the transducer effect of thermoacoustic system，to obtain the optimal relationship between the characteristic time and regenerator entropy production，and explore the influence of other factors on the entropy production rate.The result shows that with the decrease of the characteristic time，the entropy production rate reaches a very small point and then increases again.By adjusting the impedance ratio，the minimum entropy production rate corresponding to a best impedance ratio is sought.The existence of an optimal characteristic time makes the irreversible entropy production of the system minimum.At this point，the irreversible energy loss is minimal，which effectively improves the efficiency of the thermoacoustic refrigerator driven by dual sound source，and is conducive to the system driven by double sound source.

dual sound source driver；thermoacoustic refrigerator；regenerator；characteristic time；entropy production

O302；TK11

A

10.3969/j.issn.1674-2869.2016.03.016

1674-2869（2016）03-0283-05

2016-03-02

国家自然科学基金（51176143）；武汉工程大学研究生教育创新基金（CX2014029）

陈浩，男，硕士研究生.E-mail：chenhao644681298@163.com

吴锋，男，博士，教授.E-mail：wufeng@mail.wit.edu.cn