振荡管内非定常流动的数值模拟.pdf

第 卷 第 期 福州大学学报（自然科学版） 29 3 Vol . 29 No. 3 年 月 （ ） 2001 6 Journal of f uzhou University Natural science Jun . 2001 文章编号： （ ） 1000 - 2243 2001 03 - 0109 - 04 振荡管内非定常流动的数值模拟 1 2 李学来 ，朱 彻 （ 福州大学化学工程系，福建 福州 ； 大连理工大学化工机械系，辽宁大连 ） 1. 350002 2 . 116012 摘要：建立了振荡管内非定常流动与传热的数学模型，分别采用Godunov 格式和控制容积法对气动方程及传 热方程进行离散，得到了管内压力波形、气流速度及温度随时间的变化历程. 所得管内压力波在波形和相位上 与实验结果吻合较好；排气初期的气流速度最大且温度最低；管内出现强反射激波. 提出的求解Riemann 问题 迭代初值可使迭代快速、稳定收敛. 关键词：压力波制冷机；非定常流动；激波； 格式；数值模拟 Godunov 中图分类号： ； 文献标识码： TG051. 5 TK124 A ［ ］ 压力波制冷是近 年发展起来的一种气体膨胀制冷技术 1 - 6 ，这种制冷机具有结构简单、造价低、 30 操作维护方便、除气源压力能外不需消耗外加能量、制冷部分无转动部件，抵抗两相流侵蚀的能力强等优 点，因而在石油、化工、能源、航空航天及科学研究等领域中有广阔的应用前景，并已在气体分离场合得 ［］ 到应用7 . 压力波制冷机的能量转换和分离过程全部在振荡管内完成，管内非定常流动、摩擦及传热的耦 合以及波系间相互作用的复杂性和管内气流速度及温度（脉动周期达 量级）的高速变化性，使得其内部 ms 流场的研究无论从理论上还是实验上都较困难. 要揭示内部流场细节，数值模拟有望成为一条可行的途径. 振荡管内的流动为具有周期性边界条件、有摩擦和传热的非定常流动，不但会出现如膨胀波波头及 波尾这类弱间断，而且还存在激波和接触面这两种强间断. 为了处理流场中的强间断，计算流体力学领域 中先后构造出许多优秀格式，如 格式、 格式、 格式、 格式及 Godunov MacCormark Beam - Warming TVD ［ ］ NND 格式8 - 11 . Godunov 格式是由Riemann 问题的近似分析解而建立起来的，适合原有问题的物理特 性，构造简单，稳定性好，故本文采用这一格式. ! 控制方程 振荡管的基本结构见图 管内的流动不但伴随着非定常摩擦和传热，还存在波系之间、波与接触面 1 . 之间