非结构混合网格超声速化学反应流场的数值模拟研究.pdfVIP

342 非结构混合网格超声速化学反应流场的数值模拟 王兰吴颖川乐嘉陵 621000) (中国空气动力研究与发展中心吸气式高超声速技术实验室(AHL)，四川绵阳 摘要：为了解决超燃发动机和高超声速飞行器的设计中遇到的结构网格生成费时费力造成的瓶颈，吸气式 高超声速技术实验室AHL发展了基于非结构、混合网格的并行隐式流场解算器AHL_UNS3D。该解算器用 MPI实现并行，采用无矩阵LU·SGS方法耦合求解三维或二维轴对称的Euler或N．s方程及组分连续方程， 燃料带凹槽超燃冲压发动机燃烧室和支板喷氢超燃冲压发动机整机的流动，取得了很好的结果。 引言 以吸气式超燃冲压发动机为动力装置的高超声速飞行潜在的军事和经济效益吸引着世界各国对超 用阶段。然而在超燃冲压发动机的研究中，通过试验获得的流场信息非常有限，而且试验非常昂贵。 CFD能在较短时间内给出超燃冲压发动机燃烧室详细的流场结构。因此，用功能强大、高效率的CFD 软件平台进行大规模计算，结合有限的试验数据来进行研究是一种很好的途径，尤其在设计阶段，可以 大大节约人力物力财力。CARDC的吸气式高超声速技术实验室AHL发展了基于非结构、混合网格用 MPI并行的隐式流场解算器AHLUNS3D。我们已经成功实现了对无粘和层流反应流的模拟【【l】】，并应 用到各种高速湍流流动和多组分扩散流的流场模拟中【【2】】，各种算例的成功证明AHLUNS3D可以求解 从亚音速到超音速的各种可压缩流动，包括完全气体的无粘流、层流和湍流及多组分混合气体的化学反 应流场。 AHL 计算法，计算区域简化为从喷孔中心线到距喷孔中心一半孔间距的一片区域，隐含了侧向对称假定，忽 略了侧壁影响。对带凹槽氢燃料超燃冲压发动机燃烧室流动和支板喷氢超燃冲压发动机整机进行了数值 模拟，得到了冷流和化学反应流场壁面压力分布和对称面上的物理量云图，压力分布与实验值基本吻合， 表明AHLUNS3D已经具备模拟超声速氢气反应流场的能力，可以为超燃冲压发动机设计研究提供可 靠的数据。 1数值求解方法 1．1控制方程及其离散 在以加为边界的控制面力上的积分形式控制方程如下： 导OffV+OdV+fF。(Q，元)豳=f(面(Q)+j(Q)矽V+f丘(Q，5)dS÷ I ，元)豳=I(面(Q)+s(Q)矽+1只(Q， t…t) ’’n m n m 矢量，S(Q)和面(Q)分别为湍流模型和化学反应源项。采用格点格式有限体积法，对单元f积分： L纠 巧ia幼-+∑(厂iy～)△墨．，=(St+q)巧+∑(工．驴·～)△墨，， “‘ J=ribt” J=曲t 上式中，K为单元f的体积，nb(i)表示单元f的相邻单元，下标“耖”表示单元f和单元j『的交 界面，％，△S，分别界面∥的法线矢量和面积。 1．2无矩阵LU．SGS方法和边界条件 对无粘和粘性通量项及湍流模型源项和化学反应源项进行隐式处理，令豆为余项，方程可表示成： 鲁=一亩(．，邑(∥吲哦)+∥+∥+管(．，邑(厂1删·～飚√=∥ ◇ 式处理采用赵慧勇论文中的方法【【1211。 由(3)得LU．SGS左右扫描删： △研=D卅限一o·5薹[心(踢嘞)一成△踢】晦 (5) AQ,=△研一0．5D_∑(△F(g，坳)-pAAQj]ASU jEU 谱半径几=旷‘nl+a+ns茜h。三和【，分别对应 式中，对角矩阵。=晦+。．5∑J(O△岛几十巧或)， -ff?-Yn i的编号小于f和大于f的邻点的集合。△F(Q，，，zff)的计算采用无矩阵方法【【4】Ⅲ5】】． 舻(咖盯)：丝丝丝趣型s：1．0e．6 ‘