二维欧拉方程在非结构网格上的自适应有限元算法

发展了一种易于和自适应网格加密技术相结合的非结构三角形网格生成技术。这种方法的特点是域内节点可任意规定;根据欧拉方程解的误差估计,在流场某些区域增加节点后,能重新快速生成新的具有高质量三角形单元的非结构网格。欧拉方程采用Ｇａｌｅｒｋｉｎ空间离散和Ｒｉｃｈｔｍｙｅｒ两步显式时间推进相结合的有限元法求解,并以沿流线熵的变化量作为物理判据进行网格自适应加密,提高解的精度。     

脉冲发动机中金属膜片式隔舱动态破坏过程研究

为了实现固体火箭发动机的多脉冲启动功能,设计了一种金属膜片式隔舱(PSD)结构,利用断裂力学理论和圆板大挠度理论建立了金属膜片的设计公式,根据燃烧室的使用要求确定了打开压强为2.2 MPa的金属膜片尺寸,利用脆性断裂模型模拟了膜片的破坏过程,得到其打开压强为2.15 MPa,有限元数值计算结果与公式计算结果符合较好。设计了金属膜片打开破坏的单项试验,5次试验平均打开压强为2.10 MPa,数值模拟的破坏形式与试验结果一致。通过区间估计,计算了膜片真实平均打开压强置信度为90%的置信区间为[1.90 MPa,2.30 MPa],该范围可以较好满足脉冲发动机的使用要求。     

机翼积冰数值模拟

采用计算流体力学（CFD）方法对三段翼和MS-317后掠翼进行了积冰数值模拟研究,基于欧拉两相流理论,应用分区算法,对空气流场和水滴流场进行了数值求解,得到水滴收集率.求解3-D积冰模型得到积冰量.将MS-317的冰形计算结果与实验数据以及软件LEWICE3D的计算结果进行了对比.结果显示:三维积冰模型具有较好的精度,虽然预测的冰形与实验数据有一些差异,最大厚度方面,误差在13.1%以内,但是积冰的类型与体积和实验数据基本一致.      

大气温度对燃气涡轮发动机性能的影响及其修正方法

本文通过对燃气涡轮发动机在不同进气温度下的变比热特性计算,分析了按相似理论导出的常规性能换算公式的误差,并提出了较精确而实用的换算方法.本文还就发动机循环参数及工作状态等对换算公式的影响作了初步分析.     

三维可压缩流欧拉方程组的隐式时间推进法

本文提出了求解三维定常无粘可压缩流场问题的一种隐式时间推进法。在计算空间以逆变速度分量为未知变量的欧拉方程组为控制方程,从而简化壁面边界的处理。在Beam-Warming的AF-方法基础上,用对角化及矢通量分裂等方法提高了计算速度及解的精度和稳定性。作为算例,对伴有激波的方形收放喷管跨声速流动进行了计算,以显示方法的有效性。     

燃油直接喷射入横向气流中雾化机理数值研究

采用商业软件Fluent对燃油直接喷射入横向气流中的破碎雾化过程进行数值模拟,选用欧拉-欧拉方法中的VOF（volume of fluid）耦合level-set方法模拟气液两相流动.数值模拟了横向气流中燃油直接喷射射流破碎雾化过程及其相应流场信息,分析了流场结构对液柱破碎雾化过程的影响,获得了横向气流中射流破碎点位置以及射流穿透深度,并与试验数据进行了对比.数值计算结果表明:①采用VOF耦合level-set算法能够较好地模拟射流喷入横向气流中燃油破碎雾化过程;②由于反向旋转的涡对的存在,使得液膜和液滴沿着展向（y向）输运,对于液滴的二次雾化起到促进作用;③射流液柱总是在喷嘴下游约8倍直径处开始破碎,且破碎点位置取决于喷嘴几何特性.      

四参数全范围Forman裂纹扩展速率曲线p-da/dN-△K曲线拟合技术

裂纹扩展分析软件NASA／FLAGRO VERSION 2．0，采用的Forman裂纹扩展速率公式，是包括门槛区和快速扩展区的全范围裂纹扩展公式。采用四参数Forman公式可对裂纹扩展的整个过程进行描述，Fonnan公式考虑了应力幅值△K和应力比R对裂纹扩展速率的影响。通过对试验数据的拟和，可以得到Forman的四个参数C、n、p、q，进而得到各种不同材料的da／dN-△K-R曲面，再代入不同的应力比R，即可求得不同R下的da／dN-△K曲线。此外，为计算具有可靠度的疲劳裂纹扩展寿命，需采用具有可靠度的裂纹扩展速率表达式，即p-da／dN-△K表达式。     

基于非结构网格的多段翼型流场数值模拟

将非结构网格生成技术应用于多段翼型流场数值模拟，针对多连通域可以快捷、迅速地生成非结构网格。通过控制背景网格节点的网格生成参数，进行网格的局部加密，提高多段翼型流场分析的精度，应用有限体积法求解欧拉方程，应用各种加速收敛措施提高流场求解速度，并将数值计算结果与风洞试验结果进行了对比分析，证明该方法具有型号设计应用价值。     

旋翼桨叶翼型流场和跨音速翼涡干扰噪声研究

基于桨涡干扰流动的噪声场分析了NACA0012翼型的远场气动噪声.采用Euler方程计算了NACA0012、NACA23012、OA207三种翼型的绕流场,选取Kirchhoff面选为包围整个翼型的一层网格,通过得到NACA0012翼型的Kirchhoff面上在桨尖涡干扰流动下的近场解,根据Kirchhoff方法求解远场声压.样例计算分析体现出:在计算桨涡干扰噪声时,旋翼CFD技术与Kirchhoff方法相结合的良好效果.     

应用于跨声速旋翼气动声学计算的两种时域方法的比较

高速脉冲噪声(HSI)与流动非线性效应密切相关,是跨声速旋翼的主要声源之一。采用理论分析和数值实验相结合的方法,对预测HSI噪声的FfowcsWilliams-Hawkings方法(简称FW-H方法)和Kirchhoff方法进行了比较和分析。理论分析从两种方法的控制方程出发,采用广义函数方法对方程源项进行整理和比较,证实在线性区域Kirchhoff方程可视为FW-H方程的有效近似。数值实验以UH-1H悬停旋翼跨声速气动噪声计算为例,由三维Euler方程数值模拟提供近场气动数据,并在同一控制面上分别采用FW-H方法和Kirchhoff方法对远场噪声进行定量预测。计算结果表明,如果适当构造控制面,2种方法均能有效地预测跨声速旋翼非线性HSI噪声。进一步研究了控制面位置及控制面上物理参数时间导数计算精度对计算结果的影响,发现控制面位置对Kirchhoff方法计算结果具有决定性的影响,而控制面上物理参数时间导数的计算精度对FW-H方法计算结果具有重要影响。