三角形网格CE／SE方法在带锥形尾喷管两相脉冲爆轰发动机流场计算中的应用（可编辑）.docVIP

三角形网格CE／SE方法在带锥形尾喷管两相脉冲爆轰发动机流场计算中的应用 ? 格式简单、精度高、捕捉间断能力强等优点,所以近年 来/ 方法成为国内外学者研究的焦点之一。 ? 如等人 在最初提出的 / 方法的基础上进 ? 如一 行了改善,提高了它的计算精度和稳定性。如 一等人 在 的基础上将 / 方法进
一 一 步发展,针对粘性流采用长扁型网格发展了/一/ 的几种计算格式,并进行了算例验证和效率分析。 一 / 等人 运用矩形网格 / 方法对两相脉冲 午 爆轰发动机内流场进行了数值研究。李听等人 运 用矩形网格/ 方法求解流体动力学 ? 方程 与麦克斯韦方程,实现了对电磁轨道炮中等离子体电;枢的数值研究。刘建文等人。 考虑到流场的不规则性运用非结构网格 / 方法对单相多管脉冲爆轰 发动机模型进行了数值研究。 本文在 的基础上推导了带源项的结构化 三角形网格 / 方法计算格式,这方面的内容在 式中, 一 为二维流动; 一 为轴对称流动。其中下 国内还未见报道。而此结构化三角形网格 / 方 标 、 分别表示气相和液相; 、 为体积分数比,满 法与矩形网格 / 方法相比无需坐标变换,网格 足% 一 ; 、 、 、 分别是密度、轴向速度、径向 分布均匀;与非结构网格 / 方法相比推导过程速度、压力。 是单位总能, ? ÷ “ 。 、 简化,计算格式简单。本文应用此 / 方法数值 厶 研究带锥形扩张尾喷管的两相脉冲爆轰发动机的内 、 、 的含义和表达式见参考文献 。 外流场。分析尾喷管内不同填充状态时的内外 三角形网格 / 方法 流场及推进性能。 . 三角形 / 方法中守恒元与求解元的确定 理论模型/ 方法 。’ 将整个空间一时间计算区域划 计算采用轴对称两相燃烧转爆轰控制方程 : 分为若干个求解元。在每个求解元内,假设流场的变 量是连续的,并可以用 级数展开,穿过相邻求 一 ? . 解元的边界,流场的变量可以不连续。而在每个网格点对应的守恒元上,空间一时间的积分通量是守恒的。 考虑如图 所表示的交替网格,每个三角形的中『 ; 心点用实心圆和空心圆分别表示其交替网格点。令 ? 一 空间 表示 。上的所有网格点, ,, 的集合。 令空心圆点表示 时刻,实心点表示 “ 时刻。如图 所示, 表示 时刻,““ 【表示 / 时刻,表示 时 刻。那么点 对应的求解元 为, 、 和 ;对应的守恒元 为四棱柱 , 和 ? 一 一; ; . 计算格式
户 在方程中,根据对时间和空间统一处 理的思想 ,令 一 , : , 一 作为三维 空 气 动 力 学 学 报 第 卷 仕 台 乜上坷 瓦进 仃 胖 ,侍 剑 二用 彤/ 方法求解格式如下: 当 一,,或当 一, ? , , ?;走一 , ,时: 萎翼篓萎萋 , 二维三角形交替网格定义 ~ 号 , 图 二维三角形网格 / 方法中时空几何图形 一?. /△ , 卜 欧氏空间 。的 个坐标,且令 :: 一 ; 一 , , ,利用高斯定理对其积分有: 卜 卜 ,一『 卜 ,卜矿 : 十 式中,是 在 。上任意时空区域的边界, 一 , , 是时间一空间通量流密度矢量。 ÷一 十一 ,卜 一 一 引入雅克比系数矩阵 、 : 一 , 。 因 卜.为 和 均是 的齐次函数,所以 ? , ? 。 当 一 ,一 ,或当 一 , 一 一 , , 这样 进一步表示为:一 , , 。 ?; 一 , ,时: 根据求解元和守恒元的定义,对于任意的 , ,
∈, , , , ,
、 , ,
、 , , 可用其相 一 一 卜 ,一 应的离散量, ,; , , 、, , ; ,, 、~ ,一 , , ; , ,扎 近似代替。利用泰勒级数展开有: , , ; , , 一:,:,一 ~Ⅲ 十, ,一 ,
一 一 ,, 芝 ,一 一 所以方程用离散通量表示如下:, 中 ? △ 刁 : 为了方便计算引入一个新的坐标系,如图 。新 号, △,. 坐标系下取单位长度:△刁一、// 和△ 一 ?一/叫。 ,其中叫一』 , 。 △ △ : 和 。。求解如下 : 【擘 :, 一 , 若 “’一 , ,图 , 和 ,町 坐标系 ∑ 九。 。。 /∑ 九。, 其它. 一, 若 “ 一一 , , 新旧坐标系变换关系为: 一,其中一 ‘ “一奎九。 /奎九 其它 。雅克此矩阵 在新坐标系中对应为 、 ; 、 在新坐标系中对应为 、 。 式中: 是一个可调的参数,通常取 或 。第 期 王研艳等:三角形网格 / 方法在带锥形尾喷管两相脉冲爆轰发动