[发明专利]一种飞行器飞行中遭遇降雨时的数值模拟方法

摘要

本发明涉及一种飞行器飞行中遭遇降雨时的数值模拟方法，包括步骤：生成飞行器流场空间计算网格；确定飞行条件和降雨条件；求解空气相流体气动性能数据和空气相流体壁面剪切应力；求解飞行器表面局部雨滴收集率和飞行器表面局部雨滴撞击质量流量；求解雨滴撞击力系数增量；取飞行器流场空间计算网格中的翼型表面贴体网格作为表面水膜控制体，设定翼型表面同时存在气流附着点、气流转捩点和气流再附点，依次求解所有表面水膜控制体的水膜厚度；求解雨滴溅射模型的四个参数；求解等效的水膜层表面粗糙度；根据水膜层厚度更新飞行器几何形状，重构飞行器流场空间计算网格；获得更新空气相流体气动性能数据；获得降雨引起的飞行器气动性能损失数据。

主权项

一种飞行器飞行中遭遇降雨时的数值模拟方法，包括以下步骤：1)根据目标飞行器的几何外形，生成飞行器流场空间计算网格；2)根据目标飞行器的实际使用条件，确定飞行条件和降雨条件；3)根据飞行器流场空间计算网格、飞行条件和降雨条件，进行空气相‑雨滴相两相流体流动模拟计算，获得流场数据；4)根据流场数据获得空气相流体气动性能数据和空气相流体壁面剪切应力；5)利用动量定理建立模拟雨滴对飞行器表面撞击效应的计算式，求解飞行器表面局部雨滴收集率和飞行器表面局部雨滴撞击质量流量；6)根据飞行器表面局部雨滴撞击质量流量，求解雨滴撞击力系数增量；7)取飞行器流场空间计算网格中的翼型表面贴体网格作为表面水膜控制体，采用QUICK格式离散步骤4)得到的空气相流体壁面剪切应力，得到各表面水膜控制体边界处的空气相流体壁面剪切应力；8)设定翼型表面同时存在气流附着点、气流转捩点和气流再附点，分别记气流附着点、气流转捩点和气流再附点所在的表面水膜控制体为i0控制体、i1控制体和i2控制体，并初始化i0控制体、i1控制体和i2控制体的位置；9)求解i0控制体水膜厚度；10)求解i0控制体左侧界面方向所有表面水膜控制体的水膜厚度；11)检测i1控制体的存在性，若存在，更新其位置；否则，不予更新；12)求解i0控制体与i1控制体之间的表面水膜控制体的水膜厚度；13)检测i2控制体的存在性，若存在，更新其位置；否则，不予更新；14)求解i2控制体水膜厚度；15)求解i1控制体与i2控制体之间表面水膜控制体的水膜厚度；16)求解i2控制体右侧边界方向所有表面水膜控制体的水膜厚度；17)整理计算得到的所有表面水膜控制体的水膜厚度；18)建立雨滴溅射模型，求解雨滴溅射模型的四个参数：无量纲溅射凹坑高度实际修正值、无量纲溅射凹坑半径、无量纲溅射凹坑形成时间和无量纲凹坑间隔距离；19)根据雨滴溅射模型的四个参数，求解等效的水膜层表面粗糙度；20)根据步骤17)得到的所有表面水膜控制体的水膜层厚度更新飞行器几何形状，根据带有附着水膜界面的、新的飞行器几何形状重构飞行器流场空间计算网格；21)引入步骤19)得到的等效的水膜层表面粗糙度，根据重构的飞行器流场空间计算网格，以及飞行条件和降雨条件，再次进行空气相‑雨滴相两相流体流动模拟计算，获得更新流场数据，根据更新流场数据获得更新空气相流体气动性能数据；22)将步骤21)得到的更新空气相流体气动性能数据与步骤4)得到的空气相流体气动性能数据相减，再和步骤6)得到的雨滴撞击力系数增量数据叠加，获得降雨引起的飞行器气动性能损失数据。