毕业惯性管的功传输特性研究.docx

惯性管的功传输特性研究 摘要 阶梯活寒脉管制冷机是一种新型的可回收膨胀功的脉管制冷机。其特点是在不增加移动部件的 前提下，利用压缩腔和膨胀腔的中压力波的相位差回收膨胀功。本文利用ANSYS Fluent对阶梯活塞脉 管制冷机中的惯性管进行模拟。得到了惯性管中许多基本的流动规律和现象，并研究不同频率下惯性 管的输入功变化以及惯性管传输功的效率。结果显示，压力波在惯性管中传播时相位有明显的偏移。 在计算范围内，惯性管的功传输效率可达0. 93o当工作在共振频率时，其功传输效率也可达到0.89。 关键词脉管制冷机；CED；惯性管；功传输效率 0前言 脉管制冷机以其高可靠性、结构简单、长寿命等特点成为空间应用低温制冷机主要 形式。但是现在广泛应用的双向进气型和惯性管型脉管制冷机存在着脉管热端功无法回 收的问题。其热端的功以热量的形式散失在调相机构中［1］。为了解决这个问题，许多研 究者提出不同的解决方法，浙江大学甘智华等人提出了功回收级联型脉管制冷机［2］；屮 科院理化所王晓涛等人利用推移活塞回收膨胀功的混合脉管制冷机实现了 24.2%的相对 卡诺效率，是目前该温区下单级脉管制冷机的最高效率［3］。 |与此同吋，Zhu提出了阶梯活塞结构的脉管制冷机［4-5］,如图1所示。阶梯活塞与 阶梯气缸形成两个工作腔1和20工作腔2后面依次连接热端换热器3、冋热器4、冷端 换热器5、脉管6、锥形管7、惯性管8、锥形管9和气库10。气库10与工作腔1连接。 冷端换热器5的功通过惯性管7和气库10传输至工作腔1回收。其回收功的原理是压 力波经过惯性管的传输后相位发生180度的变化，文献中已经通过数值计算对该种制冷 机进行模拟，结果显示，与同尺寸的惯性管脉管制冷机相比，该制冷机的比卡诺效率可 相差5%,证明该制冷机是nJ以提高效率的。在该制冷机中，惯性管同时起着调节相位 和传输功两个作用，因此，惯性管传输功的效率直接关系到有多少功能够被回收利用。 如果惯性管传输功的效率很低，这种结构虽然理论上很好，但无实用价值，如杲惯性管 传输功的效率很高，这与斯特林制冷机或推移活塞脉管制冷机相比，将是一个结构大为 简化的新型热机。因此惯性管是这种制冷机成功与否的核心部件，对英进行研究有重要 的意义。 2 A R A 丄u 图1阶梯活塞脉管制冷机示意图 基金项目：国家自然科学基金面上项0 屮科院理化所罗二仓研究员课题组采用多个完全相同的热声机和脉管制冷机提出 了一种完全无运动的新型热声制冷技术[6],其中热声机Z间连接的谐振管同样有功传输 的功能；浙江大学甘智华等人提出的级联式脉管制冷机中，实验证明英两个脉管制冷机 之间的连接线同样具冇功传输能力，但其传输管的优化目标是调节压力波和质量流的相 位，并未对功传输效率做深入研究。本文利用计算流体力学(CFD)的方法研究惯性管 功传输特性，以及压力波在惯性管中传输时的相位偏移现象。 1 CFD计算设置 近些年利用商业软件CFD进行脉管制冷机整机的二维数值模拟已经获得了许多成 功的经验⑺先而且，二维的数值计算可以得到许多一维数值计算无法获得的信息。本 文中，利用商业软件ANSYS Fluent 16」进行数值计算卩叭 计算采用非稳态二阶方法， k-3湍流模型，压力基求解器，PRESTO!压力差值方法，PISO速度-压力耦合方式，密 度、动量、湍流模型参数以及能量方程均采用二阶迎风格式。利用G“mbit[ll]划分网格。 模型如图2所示。该模型将图1中的工作腔2-脉管-惯性管-气库■工作腔1部分提取出来。 采用二维轴对称模型。尺寸信息见表1。其屮工作腔2和工作腔1分别作为压缩腔和膨 胀腔。即压缩腔向管内输入功，膨胀腔冋收功。通过用户自定义函数(UDF)來控制移 动边界的运动，模拟工作腔活塞的运动。计算模型网格数为30900,工质为氨气，充气 压力为3.0MPa,所有壁而均为等温300K。 表1模型尺寸 脉管 ①50mmX100mm 惯性管 ①30mmX1400mm 气库 2L 压缩腔活塞的运动方程为 v = v0sin(eot) 膨胀腔活塞的运动方程为 V = Vosin(6)t + 7T) 其中，t;是活塞移动的瞬时速度，％为活塞运动速度的幅值，3为角频率，t为时间。 这样可以得到运动方向完全相反的两个活塞面，但是工作腔的容积变化规律是一 致的。因此，当压缩腔处于高压时，膨胀腔压力处于低压，即压缩腔输入功时，膨胀腔 输岀功。 压缩腔扫气容积 T yc = I Ac ? %dt 膨胀腔扫气容积 其中，Sc为压缩腔的截面积, 扫气容积比 其中，Sc为压缩腔的截面积, 扫气容积比 人为膨胀腔的截面积，丁为周期。 ye ￡ =— yc 压缩腔的PV功定义为 Wc = PcdVc 膨胀腔的PV功定义为 率变化时其功