[发明专利]一种用于抑制预冷器结霜的模拟计算方法

说明书

本发明公开了一种用于抑制预冷器结霜的模拟计算方法，可以实现对预冷器结霜全过程的运算，揭示湿空气近冷壁面结霜现象的本质及影响因素，使得湿空气结霜过程的研究更加完善；通过计算可得，在一定时间范围内，增大管间距，减小管径，降低湿空气进口流速，提高管壁壁温或者降低湿空气进口相对湿度可以有效抑制预冷器极低温管束表面霜层的形成，起到抑霜作用。

技术领域

本发明提供了预冷器结霜模拟运算技术领域，尤其是一种用于抑制预冷器结霜的模拟计算方法。

背景技术

结霜现象普遍存在于制冷、低温存储、空气调节、航空航天等领域。因此，掌握具体结霜过程以及合理解释结霜机理对于工程应用来说具有十分重要的研究意义。近年来，随着航空航天事业的不断进步与发展，对于结霜现象的研究又有了新的需求。航空航天设备表面与周围环境之间温差较大(大于100K)，多数情况设备会在大温差工况下运转并发生快速结霜现象。以预冷器为例，主流气体所含水蒸气的结霜现象是设备运行过程中一个较为严重的问题。结霜现象的存在会导致预冷器主流压力损失增大以及热交换量降低等危害，甚至引起传热失效，更有甚者还会在低空加速阶段造成灾难性事故。扩大预冷器管路间距虽然可以控制主流压力损失增大，但仍存在热交换性能下降以及预冷器尺寸增加等问题。过去有学者进行了大量关于热交换器结霜现象的研究，但大多以冷冻和空调用热交换器为研究对象，主要用来解决小温差长时间结霜问题，这与预冷器大温差短时间结霜现象的本质有所不同。在目前的研究成果中，关于预冷器结霜过程的研究十分有限。对于预冷器结霜现象的深入研究有利于其抑霜方法的探索，进而削弱结霜带给预冷器的危害。总结前人研究成果可知，虽然目前已有许多公式与试验数据较为符合，但是其结果精确度不高，应用范围具有一定局限性，无法表征预冷器的局部结霜特性；对于预冷器结霜试验而言，由于试验条件和技术手段的局限性，加之预冷器结霜过程发生的尺度较小，难以进行试验观测，更难以对其中参数进行测量，目前通过试验方法很难观测该结霜过程，因此需要借助数值模拟手段来辅助研究。鉴于结霜机理的复杂性以及结霜方式的多样性，对于预冷器结霜过程的研究而言，实际结霜过程的预测一直是研究中的难题。湿空气在大温差环境下的结霜方式不同于普通的小温差结霜过程，主要由湿空气壁面结霜和湿空气近冷壁面结霜两部分组成。现有的数学模型一般未考虑湿空气近冷壁面结霜，所有关于此类问题的数值模拟研究都是近似将该过程简化为湿空气壁面结霜来进行，这一简化使得预冷器结霜过程的模拟结果有可能产生较大误差。

发明内容

为了更好地实现预冷器结霜全过程，揭示湿空气近冷壁面结霜现象的本质及影响因素，完善湿空气结霜过程的研究，本发明提供了一种用于抑制预冷器结霜的模拟计算方法。

本发明给出了一种用于抑制预冷器结霜的模拟计算方法，将改进的多组分焓法格子Boltzmann相变模型与LBM-LGA冰粒沉积模型耦合，使用密度分布函数和温度分布函数分别模拟流体的速度场和温度场，通过求解固相体积分数追踪相变过程中各相的体积分布情况，使用固相体积分数修正焓值以及流场演化方程迁移步，采用改进的多组分伪势多相流格子Boltzmann模型解决双组分两相流大密度比问题，并考虑相变潜热和辐射换热对相变过程的影响，利用LBM-LGA冰粒沉积模型模拟形成的冰粒沉积过程。

为解决上述技术问题，本发明创造采用的技术方案是：

一种用于抑制预冷器结霜的模拟计算方法，主要包括如下步骤：

(步骤一)根据物理模型计算所需的计算域和网格尺寸，并进行网格划分，如图1所示；

(步骤二)对模型计算域的密度、温度、速度以及固相体积分数进行初始化；

(步骤三)输出计算域的流场、温度场、辐射强度以及辐射热流量；

流场演化方程为：

其中，为组分σ的粒子分布函数，为组分σ的平衡态分布函数，τ_σ为组分σ的松弛时间；