双脉冲发动机内两相流场数值模拟研究.pdfVIP

中国航天第三专业信息网第27届年会航天动力技术发展与应用学术会议论文集 双脉冲发动机内两相流场数值模拟 孙长宏牛禄张伟青 (上海航天技术研究院810所，上海200125) 摘要：采用欧拉一拉格朗日法对双脉冲发动机隔板打开后发动机内燃气二维轴对称两相流 动进行数值模拟，研究隔舱开口尺寸对脉冲发动机内流场的影响以及凝聚相射流对I脉冲绝热层 的烧蚀影响。理论分析表明：隔舱出口处存在一个回流区，开口尺寸较小时，燃烧室内颗粒相可 能被回流带到上游沉积起来，强化局部抉热和烧蚀，导致I脉冲发动机过热，该定性分析结果与 脉冲试验发动机热试后I脉冲壳体严重过熟的现象吻合． 主题词：双脉冲发动机；两相流场：数值模拟 1引言 固体双脉冲发动机(以下简称双脉冲发动机)是指在燃烧室内有两套独立的装药和点火系统，两套 装药间隔点火工作，间隔时间不受限制，工作过程中共用一个喷管。双脉冲发动机II脉冲点火工作时 I脉冲绝热层完全暴露在高温燃气下，与常规单室双推力发动机相比I脉冲绝热工况恶劣，发动机热 试车后I脉冲壳体过热现象严重[见图1]；此外，I、II脉冲之间采用隔舱结构进行功能隔离，当II 脉冲发动机按导弹指令点火工作时，隔舱结构能可靠打开，而II脉冲工作时隔舱开口大小对I脉冲发 动机内流场以及颗粒轨迹有一定的影响。本文采用欧拉一拉格朗日法对双脉冲发动机不同隔舱开口直 径时的两相流场进行了数值模拟，并进行结果分析，为发动机的设计提供了参考意见。 图1 双脉冲发动机热试验后壳体过热 孙长宏等：双脉冲发动机内两相流场数值模拟 2数学和物理模型 2．1物理模型： 一 发动机内流场为一个动态流动过程，随着燃面的退移，燃气流动通道的形状发生很大变化。由于 目前工程软件难以实现三维燃面的退移，计算不同隔舱开13直径的流场时假设I、II脉冲药柱已经烧 完，在保证计算精度的前提下，对发动机内腔形状作了适当简化。 为节约计算资源，考虑到发动机结构的对称性，本文采用二维轴对称模型，将燃气简化为单相理 想气体，不考虐化学反应和热辐射，不考虑侵蚀燃烧效应的影响和内壁面烧蚀；假设Al疵为球形，不 考虑Al疵颗粒的燃烧、蒸发和破碎，粒子的体积比小于0．1。 2．2数学模型 型，近壁处理采用非平衡壁面函数法。 2．2．1气相控制方程 固体火箭发动机的可压缩两相内流场，气相基本方程的通用形式为：’ 妄(删+毒(p叫)。毒卜毒j+s+％ ㈤ 其中：妒为任意独立变量， r。为输运系数，邑为气体源项，％为气体与颗粒相互作用的源项· 2．2．2颗粒相控制方程 连续方程： M=I。体％刚=corot ‘(2) 其中：一为通道截面积，Lk为颗粒流速垂直于流管截面的分量。 在i方向上的动量方程： 警=(ul一％舨托 (3) 其中：7t为颗粒动力驰豫时问。 颗粒相的能量方程： 鲁=(Q^一Q一级)／(豫气) (4) 式中：魏为异相反应热，幺为颗粒辐射热，—‰为女种颗粒每个颗粒的质量。 当颗粒和固体壁面碰撞时，采用Tabakoff等人提出的颗粒碰撞前后的速度公式”： (5) 等=1．0-0．4159崩～o．4994fl-2+o．292群 等乩o-2．12”3．0775，?“l卢? 4lO 中图航天第三专业信息网第27届年会航天动力技术发展与应用学术会议论文集 式中：％l、屹2和VrI、V，2分别为颗粒碰撞前后的法向和切向速度，届是颗粒碰撞前的速度和 壁面切向之问的夹角。根据Vn、K2就可求出颗粒碰撞后的角度屈： 展：arctg(vz) (6)