热声发动机系统稳定性分析 汤珂.docVIP

中国工程热物理学会 工程热力学与能源利用 学术会议论文 编号：091131 热声发动机系统稳定性分析 汤珂，林小钢，张玙，金滔，陈国邦 浙江大学制冷与低温研究所，杭州市浙大路38号，邮编：310027 （Tel：0571E-mail：ktang@zju.edu.cn） 摘要： 热声发动机系统的自激振荡起振过程是在远离平衡态条件下发生的自组织过程，是决定热声发动机能够正常运行的关键过程。本文根据非平衡态热力学，针对起振过程开展热声系统热力学稳定性研究。将Glansdorff-Prigogine稳定性判据引入到热声系统稳定性分析中，借助CFD软件，分析了两种定态的稳定性，两种定态分别对应了温度梯度小于和大于临界温度梯度的情况。验证了非平衡态热力学稳定性分析在热声发动机自激振荡起振现象研究中的可行性。 关键词： 热声发动机 非平衡态热力学 热力学稳定性 自激振荡 0 前言 近年来，热声发动机由于无运动部件、稳定性好、可利用热能驱动等优点，备受学术界及工业界的关注，在许多方面都取得了突破性进展。与此同时，仍存在着许多基本问题值得深入研究，自激振荡起振问题就是其中之一。 由于自激振荡起振过程是决定热声发动机能否正常运行的关键过程，国内外许多科研院校针对热声系统的起振和消振过程开展了相关研究工作，包括数值模拟和实验观测。然而目前的研究工作主要针对起振的动力学过程进行[1-9]，即模拟或者观测微小压力振动如何发展并最终达到稳定压力振荡的过程，包括频率模态的转换，以及各种结构形式、结构参数、运行参数、甚至操作方法等等对起振动力学过程的影响，对于起振过程的热力学前提很少涉及。 热声系统的起振过程是在远离平衡的情况下发生的非线性动态过程，是依靠能量耗散得以生成和维持的耗散结构的建立过程，是多种不可逆过程耦合作用的结果。经典热力学针对孤立系统的熵增原理，以及近平衡区线性非平衡态热力学的最小熵产生原理均无法应用于起振过程的研究。根据非线性非平衡态热力学，诸如起振过程这样在远离平衡态的条件下发生的自组织过程的前提条件是系统的热力学不稳定性，即当系统处于热力学不稳定状态，同时结合由边界条件等共同确定的系统动力学特性，才可能形成具有时空有序的耗散结构[10-11]。可见，热力学不稳定性是起振的基本前提，而动力学特性是 起振过程，即具体形成怎样的耗散结构（即系统声场）的决定因素。 在热声系统的热力学稳定性分析方面，哥伦比亚Rivera-Alvarez[12]等人利用非平衡热力学研究中的分歧曲线法、极限环等手段，研究了热声系统板叠温度梯度小于、大于临界温度梯度两种情况下，温度、压力、速度等扰动随时间的发展趋势，据此来判断热声系统的稳定性，并从中推导了临界温度梯度的计算公式，其结果同Swift[13]的推导结果相一致。随后墨西哥Ovando、Huelsz[14]等人利用Rivera-Alvarez所推导的方法进行了在热声板叠出施加磁场的磁流体热声系统的稳定性分析。虽然他们应用了非平衡热力学研究的诸多方法，但其本质上还是从研究压力、速度等动力学稳定性因素方面进行分析的。 本文从热力学角度出发，将基于非平衡态热力学的Glansdorff-Prigogine稳定性判据[10](以下简称G-P判据)引入到热声系统起振现象研究，利用CFD软件，考察给定定态（参考态）在系统中被施加了特定的扰动之后，相对于参考态的熵的二阶偏离(超熵)随时间的变化趋势，然后结合G-P判据判断参考态的稳定性。考察了两个定态（参考态）的稳定性，通过这两种不同定态的稳定性的分析，有望为热声系统发生起振过程的临界热力学条件问题的探究提供一些有益的参考。 计算模型 1.1 G-P判据 根据非平衡态热力学，对于远离平衡态的非线性非平衡区定态的稳定性存在G-P判据[10]： 稳定 介稳定 不稳定 其中为超熵产，定义如下： （1） 也就是说，当针对某定态（参考态）引入微小扰动后，计算得到，则随着时间的推移，微小扰动将逐渐消失，系统回复到之前的定态，即系统为稳定；如果，则意味着微小扰动将随时间发展增大，系统不能回复原定态，而是向新的定态发展；而表示该定态为介于稳定和不稳定之间的临界状态。以上三种状态分别对应图1中的三条超熵曲线，根据式（1）可知，图1中超熵曲线的斜率即超熵产。 图1 非平衡定态稳