一种不利高过载诱发的喷雾冷却表面液膜畸变的建模方法
【技术领域】
[0001] 本发明涉及一种不利高过载诱发的喷雾冷却表面液膜畸变的建模方法。
【背景技术】
[0002] 喷雾冷却作为一种高效的换热方式，在航空航天领域的大功率电子系统中有十分 广阔的应用前景。电子设备在机载环境下将面临不同方向上的过载影响。以飞机为例，机 动飞行过程中所受载荷比水平直线稳定飞行时大很多，而且飞机机动性要求越高，所承受 的过载也越大。飞机过载系数是指飞机所受除重力之外的外力总和与飞机重量之比，战斗 机在俯冲时的过载系数可以达到9~10。机动飞行动作不同则所受的过载力大小和方向也 不相同。喷雾液滴撞击热源面后将在壁面上形成一层薄液膜，该液膜对喷雾冷却的换热性 能有重要作用。过载力的存在将会引发表面液膜畸变，改变液膜在热源面上的分布和流动 形态，进而影响喷雾冷却换热效果。对于过载系数较大的高过载环境，这一作用效果将更为 严重，且有可能对换热造成不利影响。因此，针对喷雾冷却应用于航空航天领域的可行性， 有必要研究不利高过载对喷雾冷却表面液膜畸变的影响。
[0003] 目前研究重力因素对喷雾冷却性能的影响多以实验方法为主，常用的方法有两 种。一是利用飞机的抛物飞行，结合调整热源面的朝向，为喷雾冷却创造超重和失重环境。 二是利用落塔装置实现微重力环境的模拟。两种方法可以获得的实验时间很短，如飞机的 抛物飞行实验约为20s，落塔实验只有几秒。相关实验可以直观地反映重力因素对喷雾冷却 换热的影响，但由于实验时间很短且实验条件有差异，实验现象观察不充分，结论存在不一 致性，不利于非常重力环境下喷雾冷却换热机理的研究。
[0004] 另一方面，针对喷雾冷却换热机理已开展大量的数值模拟研究，多以单个或有限 个液滴撞击壁面或表面液膜的流动过程为研究对象，观察不同初始条件对结果的影响。喷 雾冷却是大量液滴源源不断撞击壁面后形成液膜的流动过程，有限个液滴撞击并不能很好 反映出热源表面液膜的形成和流动形态。目前较为常用的计算模型为VOF多相流动模型。 VOF模型多被用于计算常重力环境下的液滴撞击热源面问题，几乎没有被用于高过载环境。
[0005] VOF模型是一种多相流模型，通过求解单一的动量方程和追踪控制体内流体的体 积分数来模拟两种或多种互不相融的流体的流动。当在模型中引入新的一相，则会加入该 相的体积分数作为新的变量，所有相的体积分数之和为1。若各相的体积分数已知，则所有 变量和属性由各相共享且以均值来体现。
[0006] 质量守丨旦方程：▽
[0007] 体积分数方程：
[0009] 物质属性方程：P = Σ a i p i CN 105117558 A 说明书 2/5 页
[0012] 其中，$为速度，α为体积分数，t为时间，S为质量源项，P为密度，m为质量传 递，Φ为混合属性，P为静态压力，μ为动力粘度，,为动力源项。
【发明内容】
[0013] 根据本发明的一个方面，提供了一种不利高过载诱发的喷雾冷却表面液膜畸变的 建模方法，其特征在于包括：提出过载加速度的模拟方法；建立喷雾状态模拟方法；建立不 利高过载导致液膜堆积现象发生的判断方法；建立不利高过载导致液膜破碎现象发生的判 断方法。
[0014] 根据本发明的另一个方面，提供了一种高过载环境下液膜畸变的确定方法，其特 征在于：
[0015] 提出过载加速度的模拟方法，在VOF多相流动模型中引入过载加速度，建立过载 环境下的VOF计算模型。
[0016] 建立喷雾状态模拟方法，将喷雾模型简化为多列液滴撞击热源面模型，以三个液 滴参数控制工质流量和喷雾疏密程度；
[0017] 建立不利过载导致液膜堆积现象发生的判断方法，通过对热源壁面节点对应的液 膜厚度进行检验，判断液膜堆积现象发生的位置；
[0018] 建立不利过载导致液膜破碎现象发生的判断方法，通过对热源壁面节点对应的液 膜厚度进行检验，判断液膜断裂现象发生的位置。
[0019] 说明书附图
[0020] 图1显示了根据本发明的一个实施例的方法的组成步骤。
[0021] 图2为根据本发明的一个实施例的模型实例。
[0022] 图3为根据本发明的一个实施例的喷雾状态简化图。
[0023] 图4为根据本发明的一个实施例的不利过载导致液膜堆积现象发生的判断流程。
【具体实施方式】
[0024] 本发明提供了一种不利高过载诱发的喷雾冷却表面液膜畸变的建模方法，解决实 验条件限制导致的实验时间过短及现象观察不充分等问题，对常重力条件下的VOF模型加 以修正，建立过载条件下的计算模型，并以液膜堆积和破碎的判断方法为依据，确定不利过 载引发液膜畸变现象发生的位置，研究过载影响。
[0025] 根据本发明的一个方面，建立了过载条件下的数值计算模型；对喷雾状态进行模 拟，其中对喷雾模型进行了简化，以三个液滴参数控制工质流量和喷雾疏密程度；对液膜畸 变现象发生进行判断，确定在一定过载条件和喷雾状态下，液膜堆积和破裂现象发生的位 置。
[0026] 图1显示了根据本发明的一个实施例的一种不利高过载诱发的喷雾冷却表面液 膜畸变的建模方法，包括：
[0027] (1)建立过载环境下的算法模型
[0028] 以常重力环境下的VOF模型为基础，将参考系建立在水平加速运动的热源面上， 引入平移惯性力，建立过载环境下的算法模型。确定计算区域，进行时间及空间上的网格划 分，设置边界条件。（图1的步骤101)
[0029] (2)对喷雾状态进行模拟
[0030] 将喷雾模型简化为多列液滴撞击壁面模型。以液滴直径、液滴速度和液滴间距为 描述喷雾状态的三个重要参数，控制工质流量和喷雾疏密程度。（图1的步骤102)
[0031] (3)判断液膜畸变
[0032] 根据计算的液膜厚度，判断液膜畸变，确定由不利高过载引发的液膜堆积和破裂 现象发生的位置。（图1的步骤103)
[0033] 下面结合二维实例介绍本发明研究液膜畸变方法的具体实施步骤。如图2所示， 计算区域宽度为1〇_，高度为5_，空间步长Δ X和Δ y均为〇. 1_，网格数量为5 X IO3个。 左右两端为压力出口，顶端为压力入口，底部为热源面壁面。边界大气压力为〇· IMPa，温度 为27°C。空气为主相，液态水为第二相，水的表面张力系数为0. 073N/m，水与壁面的静态接 触角为90°。
[0034] (1)过载环境下算法模型的建立
[0035] 在地面惯性系中热源面以加速度a。运动，将坐标系建立在热源壁面上，则确定了 一个非惯性系。在此非惯性系中引入惯性加速度a，其大小与a。相等，方向与a。相反。据 此，液滴撞击以加速度a。运动的热源面问题，转换为了处于惯性力场a中的液滴撞击固定 热源面问题。
[0037] 过载环境下VOF模型中的动量方程变为如下形式
[0039] 其中P为密度；
[0040] μ为动力粘度；
[0041] 交为速度；
[0042] t为时间；
[0043] P为静态压力；
[0044] f为重力加速度；
[0045] 5为引入的惯性加速度；
[0046] J5为表面张力产生的动量源相。
[0047] (2)模拟喷雾状态
[0048] 喷雾状态直接影响液膜在表面上的流动形态，与多个参数有关。根据本发明的一 个实施例，将喷雾简化为多列液滴撞击热源面问题，对于某一喷雾状态，假设所有液滴直 径、液滴间距和液滴速度固定，且液滴撞击热源面时的速度垂直于热源面，如图3所示。其 中，喷雾过程中液滴源源不断以一定初始条件撞击热源面。把液