脉管制冷机的全三维数值模拟计算(下).pdfVIP

中国T程热物理学会 T程热力学’j能源利用 58 第f．_一扁学术会议 编号：05ll 脉管制冷机的全三维数值模拟计算(下) 丁文静何雅玲母 陶文铨 西安交通大学动力工程多相流国家重点实验室，能源’j动力工程学院，西安，710049 Tel：029-8266385lE．mail：valinghe@ma．1．xitu．edu．cn 摘要： 本文上篇阐述了脉管制冷机全三维数值计算程序的开发，详细地给出了物理模型和控制方程。 下篇通过计算，详细地分析脉管内的流场和温度场以及其他各主要参数的变化，以便于对脉管制冷机 的_T作过程和制冷机理有进一步的了解和认识，为脉管制冷机的结构改进和实用化奠定一定的基础。 关键词：温度场，温度滞后，流场，二次流，环形效应，制冷量 1前言 本文上篇详细地给出了脉管制冷机的物理模型、控制方程和数值方法。下篇详细分 析脉管制冷机的温度场、流场及其他主要参数的变化，探讨脉管制冷机的制冷机理。 2一个周期内截面平均温度随时间的变化 图1为脉管制冷机轴向不同截面上平均温度随时间的变化曲线，图中前半周期为压 缩过程，后半周期为膨胀过程。 图1一个周期内脉管轴向不同截商、F均温度随时问的变化曲线 从图l可以看出，脉管截面平均温度沿轴向的分布具有以卜．特点：(1)在靠近冷端 位置，温度随时间变化平缓，在靠近热端区域温度变化幅度较火。(2)脉管冷端温度在 一个计算周期内，有时高于冷端换热器温度，有时低于冷端换热器温度；从数值上可以 看出，温度低于冷端换热器的幅值较大，高于冷端换热器温度的幅值较小，这也是产生 周期制冷量的必要条件。(3)就常规流动而言，在压缩过程中，．：亡质应该在热端换热器 资助项目(No．G2000026303)。术通讯作者 683 内放出热量：在膨胀过程中，脉管冷头，L：质温度应该低丁．冷端换热器设定温度，此时在 冷端换热器内有冷量产生。而实际上在脉管制冷机压缩过程的起始阶段，脉管热端温度 仍低于热端换热器设定温度，没有热量传出；在膨胀过程的起始阶段，脉管冷端温度反 而高于冷端换热器温度，也没有冷量在冷端换热器产生。由图1可以清楚地看山，在脉 管内压缩放热和膨胀吸热效具有滞后性，其压缩和膨胀效果大约滞后1／4个周期才能体 现出来。 3纵截面速度场分布示意图 图2给出了脉管纵截面周期平均质量流量的分布。由图2可以看出，脉管周期平均 质量流量并不等于零，在全场存在比较复杂的流动分布，具体存在以下特点：(1)壁面 处速度和常规流动不同，常规流动壁面处速度较小，而脉管内的交变可压缩流动中壁面 处存在相对较大的轴向流速；(2)全场明显存在多个较大的回流区，称为二次质量流动， 这是弓f起脉管内流体扰动的主要原因。(3)脉管的轴向中间位置是反向流速的交汇点， 即周期静流量滞．1|==点，也是能量损失的最大部位。(4)轴向速度沿径向的分布比较复杂， 不同轴向位置分布不同，在同一轴向位置，轴向速度沿径向也多次改变流向。(5)壁面 处周期平均质量流餐和中心区域反向，且数值上较大。 n们 ∞ 至≮隧季}}芗于于于f于芗f}芗fff芗 一 、 ……7 7 ’ ’ - - · 一 一 - - 、 ～ 0：!≤；i圣三三；!蔓：jj；一一!ii三i ∞E、叱 ：j；；；三三三三；i—i‘i；：o；；；{誓；三三 - · - - · - - · 一 ’ 一 ： ‘_…～ … … ～