【技术】微型低温节流制冷器结构优化设计
J-T节流制冷器被广泛用于红外探测、热成像技术和冷冻手术刀等设备的制冷设备。在有限结构内达到冷量最大化，是节流制冷器在小型化发展中面临的问题之一。作为红外寻的系统的辅助设备，微型节流制冷器的工作性能和结构参数都受到严格限制。在有限结构内为红外芯片的工作环境实现冷量最大化，是节流制冷器在小型化发展中面临的问题之一。
Hampson型换热器作为节流制冷器中的主要部件，其结构参数对制冷器的性能有至关重要的影响。本文中研究的换热器主要由翅片毛细盘管、芯轴和外壳部分组成，基本结构如图1所示。高压高温流体通过入口进入换热器毛细管内，在流动过程中被外侧回流的低温低压流体预冷，流出换热器的高压高温流体经过节流元件的节流作用，变为低温低压流体；节流后的流体吸收蒸发腔内的热量，通过芯轴和外壳间的环形通道排出。
在进行数值计算时，采用有限体积法划分网格，螺旋翅片管的网格与高压流体的网格划分方法一致，芯轴和外壳的网格与低压流体的网格划分方法一致。具体的网格设定、模型简化以及边界条件参照文献。待求解微元的基本控制方程为：
式中：m为质量流量，g·s－1；x为空间步长，m；r为流体密度，kg·m－3；u为流体流速，m·s－1；p为流体压力，Pa；τw为壁面剪切应力，kg·m－1·s－2；ps为湿周，m；A为管道通流面积，m2；h为流体焓值，J·kg－1；a为对流换热系数，W·m－2·K－1；Tf为流体温度，K；Ts为固体温度，K；k为固体导热系数，W·m－1·K－1；As为固体导热截面积，m2；Qrad为单位长度辐射传热量，J·m－1。
公式(1)～(3)分别为流体的连续性方程、动量方程和能量方程，公式(4)为固体的能量方程。公式(1)～(4)中的摩擦因子和换热系数相关关联式参照由Timmerhaus等人总结的经验公式。
节流制冷器的结构优化问题，即为数学函数求解最优值问题。对于本文中的结构优化设计，是致力于分析换热器中所有几何参数变化时对换热器性能甚至制冷器性能的影响，且这些几何参数间对最终的结果影响作用复杂，因此本文中将各几何因素对性能的影响假设为目标函数，通过对此目标函数求解最优值来确定最优的几何结构。考虑到遗传算法全局搜索能量强，能够很好地处理约束、跳出局部最优，得到全局的最优解，本文选择遗传算法作为结构设计优化方法。
3.1单因素性能影响分析
在对换热器结构进行单因素影响分析时，假定换热器芯轴的尺寸为固定值，其余参数受到所选单因素参数的影响。对换热器冷端的条件设定为冷流体入口为饱和气态，换热器冷端的冷量为冷端冷流体入口能量与热流体的出口能量之差。图2为换热器中的单因素几何结构对换热器冷端冷量、熵产、冷端入口压力和入口温度的影响。
从图中可以看出，各单因素参数的变化对换热器的熵产、冷端入口温度影响并不显著。螺旋圈数（即毛细管长度）、螺旋直径、翅厚和翅高增加时，换热器冷端的冷量都随之增加，但翅间距增加导致换热器冷端冷量的减少。各参数变化时，对换热器冷端入口压力的影响比较显著。螺旋圈数（即毛细管长度）、翅厚增加，冷端入口压力呈现明显增加趋势。在相同质量流量下，毛细管长度的增加增大了冷流体的回流流程，而翅片厚度的增加则加大了冷流体流动的阻力，这两者作用的结果就是冷流体侧入口压力增加，这也导致了冷端温度的升高。
翅高和螺旋直径增加，换热器冷端入口压力降低，同时，冷端入口温度减小，这有利于实现换热器在冷端温度较低时满足较大的冷量需求。翅片节距的增加，导致换热器冷端冷量、进气压力急剧下降。以上所有参数中，螺旋圈数与翅厚对换热器的影响变化规律和螺旋直径与翅高对换热器的影响变化规律比较类似。因此本文中以翅片参数（翅片高度、翅片节距、翅片厚度）为例进行优化，并对最优结构进行系统的校核和验证。
3.2换热器最优结构
根据前文所述，翅片参数的待求解区域分别是：翅片高度为0.1～0.5mm、翅片厚度为0.1～0.6mm、翅片节距为0.1～0.6mm。制冷工质为氩气，换热器的进气端工况为：进气压力18MPa、进气温度300K。流量设定为0.3～0.6g/s。
以换热器的熵产最小和冷端冷量最大为优化目标，得到各个工况下的最优结构如表3中所示。从表中结果可知，当进气工况和流量确定时，换热器冷端冷量最大值和系统的熵产最小值出现在翅片厚度最小值和翅片高度最大值。对比不同质量流量的工况，可以发现，翅片节距基本为翅片厚度的2倍左右。因此，在进行换热器的结构设计时，对于肋片参数可以将翅片节距设定为翅片厚度的2倍，以达到较好的换热效果。
3.3最优结构校核
从换热器中流出的高压流体，经过节流装置的等焓节流过程成为低温低压流体。节流制冷器中采用的节流孔结构与制冷器系统的质量流量可以通过公式(5)确定。
根据节流前的流体状态，可以得到确定流量下的节流孔直径。
当进口工况确定，满足上述表3中工况的节流孔直径分别为0.079mm、0.086mm、0.093mm、0.106mm、0.113mm、0.121mm。以表3中的翅片参数和小孔参数组合为新的制冷器，假定气库容积为300cm3，当气库初始压力为18MPa时，得到的流量变化与降温曲线如图3、图4所示，结构1～6分别对应表3中工况下得到的最优结构。
从图3中可以看出，在制冷器启动过程中，随着气库放气，制冷器的流量会出现短暂的上升趋势，随后很快转变为下降。节流孔直径逐渐增大，系统流量的上升幅度也较大。在制冷器运行一段时间后，不同结构对应的流量关系发生逆转，节流孔较大的结构，质量流量相对较小。这是由于初始阶段，节流前流体的压力还比较高，质量流量受小孔面积的影响比较大，等制冷器运行一段时间后，流量大的结构中的气库气体减少迅速，导致整个节流前流体的压力也相对较低，因此呈现出节流孔大的结构质量流量较小的结果。
如图4中降温曲线所示，节流孔较大的结构对应的冷端温度值较高、节流孔小的结构对应的冷端温度之较低，但不同结构的降温时间并没有确定的变化规律。节流孔最小的结构1进入恒定冷端温度区间时间最长，对应的冷端温度最低；结构3进入恒定冷端温度区间的时间最短，对应的冷端温度值比结构1和结构2略高。当冷端温度进入稳定阶段后，温度值基本保持稳定。随着制冷器工作时长增加，气库内体积变小、质量流量减少，节流孔最大的结构7冷端温度最先出现上升，脱离恒定温度区间的时间也最短。
图5为气库容积300cm3时，不同结构下对应的实时冷量变化。从冷端的冷量数值来看，冷端温度进入恒定后，冷量开始逐渐上升，在25s左右，所有结构对应的冷量值达到最大。此时，结构3的冷量值最大，结构2、4、5的冷量最大值比较接近，结构1和结构7的冷量最大值比较接近。随后，冷量进入迅速衰减期，在25～50s时间内，相同时刻的冷量大小分别对应的结构为结构3、2、1、4、5、6；在50s以后，相同时刻的冷量大小对应的结构分别为结构1、2、3、4、5、6。图5的冷量值呈现出这种结果，主要是受到系统流量和冷端温度的双重影响。
本文中建立了J-T节流制冷器一维稳态数值模型，对制冷器中换热器的几何结构进行了单因素性能影响分析，结果表明换热器的翅片参数（肋高、肋厚和肋间距）和螺旋参数（螺旋直径、螺旋圈数）对换热器的性能具有明显的影响。利用遗传算法主要对节流制冷器中的换热器翅片参数（肋高、肋厚和肋间距）进行了优化，结果表明：肋厚小、肋高大时，换热器能达到较大的冷端冷量和较小的熵产，且肋间距与肋厚之间存在2倍的最优结构。在气库容积为300cm3，初始压力为18MPa的工况下，肋厚0.1mm、肋节距0.222mm、肋高0.492mm、小孔直径0.093mm的结构组合，能够实现制冷器快速降温，能较长时间将冷端温度维持在90K以下，且冷端能提供的冷量较大。