[发明专利]一种火箭发动机真空干扰羽流场仿真方法

摘要

一种火箭发动机真空干扰羽流场仿真方法，步骤为：(1)采用CFD计算流体力学方法计算火箭发动机喷流在10Pa环境压力下的喷流干扰流场；(2)在流场未受干扰的喷流核心区域选取密度等值面作为三维DSMC计算的喷流初始边界；(3)进行喷流干扰流场的三维DSMC计算，实现所述火箭发动机真空干扰羽流场仿真。该方法克服了发动机真空干扰羽流场中同时存在连续流、过渡流和自由分子流，不能采用单一数值方法进行仿真预示的难题，结合了直角网格的高效率和物面三角形网格对复杂边界的精确描述，提高预示精度和计算效率，成功实现发动机真空干扰羽流场CFD/DSMC相结合的仿真预示。

主权项

一种火箭发动机真空干扰羽流场仿真方法，其特征在于步骤如下：(1)根据火箭发动机喷管燃烧室出口参数和喷管型面，采用CFD计算流体力学方法计算火箭发动机喷流在10Pa环境压力下的喷流干扰流场；(2)在流场未受干扰的喷流核心区域选取喷流密度最小的密度等值面，将密度等值面上的流场参数作为三维DSMC计算的喷流初始边界；所述流场参数包括燃气密度、速度和温度；(3)进行喷流干扰流场的三维DSMC计算，实现所述火箭发动机真空干扰羽流场仿真；所述进行喷流干扰流场的三维DSMC计算具体为：(3.1)划分三维DSMC计算区域的直角网格，确定流场区域边界；对物面边界划分三角形非结构网格；所述流场区域边界包括出口边界、入口边界、物面边界和对称面边界；(3.2)进行流场初始化，流场初始化包括：(a)计算所述物面边界三角形非结构网格的外法向矢量、形心位置和三角形面积；(b)将所述物面边界三角形非结构网格放入喷流干扰流场区域的直角网格，确定并标识飞行器放入计算区域后占据的流场网格、物面三角形面元与流场网格相交的物面边界网格，建立物面边界网格与其相关物面三角形面元的链表关系；(c)设置所有网格内的模拟分子数量，给所有模拟分子编号；(d)在喷流干扰流场区域内，设置模拟分子的运动参数，运动参数包括模拟分子的初始位置和运动速度；(e)在入口边界处，按照入口边界条件确定在时间步长Δt内进入喷流干扰流场区域的新模拟分子数量以及运动参数；所述入口边界即为三维DSMC计算的喷流初始边界；所述入口边界条件包括喷流初始边界的燃气密度、速度和温度；(3.3)对于喷流干扰流场区域内的模拟分子，按照其自身的运动速度和初始位置运动一个时间步长Δt后，确定模拟分子到达的新位置；(3.4)判断模拟分子运动后到达的新位置与喷流干扰流场区域边界的关系；如果该新位置到达对称面边界，则模拟分子在对称面边界上作镜面反射；如果该新位置到达物面边界，则按完全热调节的漫反射模型，确定反射之后的模拟分子的运动参数；如果该新位置到达出口边界外，则将模拟分子作逸出处理，消除该分子的编号；(3.5)重新为所有处于计算区域内的模拟分子以及新进入计算区域的模拟分子进行编号，并记录所有模拟分子的运动参数；(3.6)在每个喷流干扰流场区域的直角网格内选取可能的碰撞分子对，并由随机数判断真正发生碰撞的分子，计算碰撞后的分子速度及内能；(3.7)重复执行步骤(3.3)到(3.6)，运行N个时间步长之后，N为正整数，判断喷流干扰流场是否达到稳定状态，如果达到稳定状态，便对喷流干扰流场区域的所有直角网格内的模拟分子的数量以及运动参数进行统计计算，得到喷流干扰流场的燃气密度、速度和温度；否则返回步骤(3.3)；所述稳定状态是指喷流干扰流场中的模拟分子总数的变化率不大于5％。