上端壁翘曲对涡轮叶栅流场的影响

利用N－S方程求解程序，对采用不同端壁结构的三套大转角涡轮叶栅进行了数值模拟，探讨了上端壁翘曲对涡轮叶栅内气流流动的影响。计算结果表明,上端壁翘曲后，上半翼殿内的二次流流场发生了明显的变化，非轴对称端壁是控制叶栅内二次流流动的一种有效的方法。     

机翼和叶栅非定常流的Kutta条件

本文在现有的实验研究和数值计算的基础上，对经典Ｋｕｔｔａ条件在非定常流动中的适用性作了进一步探讨，并对文献〔１，２〕中提出的一种新的Ｋｕｔｔａ条件－Ｇｉｅｓｉｎｇ－Ｍａｓｋｅｌｌ（Ｇ－Ｍ）模型进行了分析与比较，本文给出了二维振荡机翼和叶栅绕流的变域变分有限元解法中Ｋｕｔｔａ条件的处理方法。本文旨在为机翼和叶栅非定常绕流的数值计算提供更合理可靠的尾缘边界条件。     

H型叶栅通道流计算网格的改进

提出一种壁面正交、尾迹区网格节点分布逐渐均匀化叶栅通道网格生成方法。为了研究壁面处网格正交性对数值计算结果的影响，推导出壁面附近流动参数一阶导数离散误差与网格线交角的数学关系。采用N-S方程分别对平板紊流附面层流和叶栅通道流进行计算分析。结果表明：壁面处网格正交性对附面层内流动计算预测精度影响较大，而对势流区流动影响很小；尾迹区切向网格节点逐渐均匀化分布能更准确计算尾迹区流动。     

平面叶栅跨音绕流时间推进法中特征相容条件在边界点上的应用

在平面叶栅定常跨音绕流计算中,采用时间推进法,由于基本方程组是非定常欧拉方程组,属于双曲型,可以应用特征理论。在边界点上,特征相容方程精确地刻画了流动参量间的关系,使得数值处理更精确地反映无粘流的物理模型。本文从基本方程组出发,推导出适合于边界点应用的特征相容方程组的通用形式;并应用特征理论,讨论了平面叶栅绕流问题中边界条件的提法;对于边界条件的几种不同的近似处理方法,进行了数值试验,比较了它们对数值结果的影响。     

弦向缝隙叶栅对边界层分离的控制

本文提出了一种在Ｓ１流面上进行势流与粘性流迭代求解弦向缝隙叶栅流场的计算方法，将Ｓ１流面上势流方程组的解与经Ｉ１１ｉｎｇｗｏｒｔｈ－Ｓｔｅｗａｒｔｓｏｎ变换后的边界层方程结合起来联立求解弦向缝隙叶栅流场。计算值同实验结果相当吻合，在正大攻角时，向缝隙叶栅能效地控制边界层分离。     

基于叶栅型预旋喷嘴的计算方法

随着发动机涡轮前进口温度越来越高,设计好的预旋结构,并准确地进行预旋结构的流体动力计算尤为重要。针对发动机中普遍使用的叶栅型预旋结构,采用商业软件CFX对某种叶栅型预旋喷嘴流动进行了数值模拟计算;根据2种叶栅型预旋喷嘴的试验结果,分析了叶栅型预旋喷嘴的流动特点,同时总结了1套关于叶栅型预旋结构的流体动力计算和预旋相关参数的计算方法,并编制了计算程序。结果表明:用该方法计算得到的预旋相关参数变化规律与真实物理特性相一致。     

一种新的适用于叶栅绕流计算的网格生成法及其应用

本文对平面叶栅气动计算网格生成方法作了简要评述。在此基础上,根据叶栅槽道流动的特点,提出一种以等比级数为基础,用几何方法构作H网格的数值方法。用这种网格进行流场计算的结果表明,所生成的网格在叶片前后缘区有很高的分辨能力,能够算出速度分布的细微结构。本文的计算结果与亚音速及跨音速叶栅速度图解符合得很好。     

叶栅参数对反推力装置气动性能影响规律

针对某型叶栅式反推力装置,开展了气动性能数值预测的研究工作。采用本文提出的数值计算方法,在不同落压比下反推力预测值与实验值最大相差8．2％。流场速度等值线分布展示反推力装置叶栅内部流动特征．叶栅叶盆处的旋涡和叶背处的失速是导致气动性能下降的不利因素。大量数值计算结果揭示了叶片进入角、出流角、稠度等结构参数对轴向反推力、流量系数、总压恢复系数的影响规律。在本研究范围内,叶片进入角54°、出流角135°、稠度1.3的叶栅几何结构综合气动性能最优。     

一种新型多维通量函数在求解平面叶栅跨音速流动的应用

采用一种改进的Ｒｏｅ格式对平面叶栅的无粘流动进行了数值模拟,并与Ｒｏｅ格式的计算结果进行了比较。采用四阶Ｒｕｎｇｅ－Ｋｕｔａ时间差分格式,空间三阶有限体积法。网格生成采用贴体网格近壁面修正的方法。计算结果表明,该格式在减少格式对网格的依赖性,抑制数值粘性方面具有明显的作用。     

有间隙涡轮叶栅壁面流谱实验与数值模拟研究

对具有 3 .6 %相对叶顶间隙的涡轮叶栅进行了壁面流动显示和三维流场数值模拟 ,分析了大叶顶间隙涡轮叶栅的壁面流动特点 ,从数值模拟上特别对尾缘附近的壁面流动结构进行了详细讨论。实验与计算结果表明叶栅尾缘附近的流动是非常复杂的 ,同时由于叶顶间隙的存在 ,上下尾缘附近的壁面流谱明显不同     

用全位势方程计算侧滑机翼的跨声速绕流特性

本文对大后掠小展弦比细长机翼，由垂直于机翼纵轴的每一横流截面生成Ｏ型网格，从而形成对机翼流场的Ｈ－Ｏ型网格，用守恒型全位势方程，差分方法和隐式近似因式分解ＡＦ２迭代算法计算侧滑机翼的流场。     

跨音速喷管流动的数值模拟

数值方法采用的差分格式是具有高分辨率和快速稳定收敛性质的全变差衰减(TVD)格式。喷管壁面形状和反压可任意给定,对于典型轴对称与二维喷管(即壁面由收缩角为45°,扩张角为15°的两条直线用圆弧光滑连接而成。圆弧的无量纲半径等于0.625,其参考长度为喉道半宽),进行了有激波与无激波流动的数值模拟。计算结果与实验结果以及其他计算结果符合良好,此法可推广到非定常跨音速喷管流动的计算,并可用于工程中喷管设计。     

水滴撞击特性的粒子统计法研究

基于水滴撞击特性研究的拉格朗日方法,提出了适合于考察非定常情况的粒子统计法.通过统计非定常流场中不同时刻撞击在防冰表面上不同区域的粒子数,获得防冰表面上粒子浓度分布随时间的变化情况,与基于染色法水滴撞击试验获得的色度分布具有直接对应关系;采用基于染色法的动态图像测量方法对某型发动机帽罩缩比模型开展了水滴撞击特性试验,对所提粒子统计法进行了验证.基于拉格朗日方法的粒子统计法获得的计算结果与试验结果符合较好,与采用拉格朗日轨迹法获得的结果具有较好的一致性.粒子统计法所模拟的物理过程与试验过程一致,在处理湍流随机性和来流液态水含量分布不均等方面具有优势.     

旋转离心叶轮内湍流的非结构化网格数值计算

给出了同位非结构化网格数值求解旋转离心叶轮内湍流流场的SIMPLE算法。网格划分采用改进的阵面推进法。控制方程采用通用的对流扩散方程并利用有限控制容积法离散。湍流模型采用标准的k-ε模型。形成的代数方程组用带有预条件器的共轭梯度平方法（CGS）求解。最后利用所给出的算法对一实验用的旋转离心叶轮进行计算。计算结果和实验数据基本相符,说明该算法合理,编程正确,可作为流场预测的理论工具。     

气体运动论BGK格式的翼型绕流数值模拟

以气体运动论BGK(Bhatnagaar-Gross-Krook)格式为基础,通过改进非均匀网格下的通量计算方法,把气体运动论BGK方法应用到二维贴体网格上模拟了可压缩的无黏和黏性流动。在重构阶段,以交界面周围4个网格单元为一组,先把每个网格单元的宏观量(密度、动量和总能量)转换在相对于交界面的坐标系下,然后使用van Leer限制器插值求出交界面的宏观量,在求解阶段再使用BGK理论求出流过交界面的通量。通过上述方法使贴体网格和均匀网格,在计算格式上一致,与其他求解二维贴体网格的气体运动论BGK方法相比,该方法相对简单,计算量减小。结合贴体网格通量计算和气体运动论BGK格式方法理论,分别对亚、跨、超声速下的NACA0012翼型绕流进行了数值模拟和对比研究,所得可压缩无黏流动的数值结果与北大西洋公约组织航天研究与发展咨询组(AGARD)的结果吻合良好。本文还模拟了低雷诺数状态下的可压缩黏性流动,所得结果与文献中的数值计算结果吻合良好,证明了该方法在数值模拟及激波捕捉方面具有鲁棒性,为在复杂流动中的应用提供了一种简单的途径。     

提取流管实验装置中声衬声阻抗的模态方法

发展了一种模态匹配方法进行有均匀流条件下敷设声衬的二维矩形流管中声传播的计算, 计算结果与NASA流管实验结果和CAA计算结果符合得很好.在此基础上, 着重介绍采用求极值的变尺度法优化目标函数和给定求解范围寻找目标函数最小值两种方法进行管内未知声阻抗的提取, 给出了算例验证, 并分析了随机误差和已知声压点数对计算结果的影响.      

基于笛卡尔网格的非定常气动力计算

在非结构自适应笛卡尔网格基础上构建了一种基于笛卡尔动网格的非定常气动力计算方法。为了适应非定常流场对网格品质要求较高的特点,网格生成采用了基于几何外形、模型表面曲率为基础的自适应生成及加密方法。采用显格式龙格库塔推进方法和中心格式空间离散的有限体积法对NACA0012翼型,三维机翼的跨音速刚性振动进行了数值模拟。对伪时间收敛精度控制,真实时间步长对结果的影响进行了讨论。将计算结果与试验数据进行了比较,得到了良好的效果。     

一种计算非定常二维流动的无网格算法

主要目的是发展一套求解非定常流动的无网格算法。计算区域的离散方面,提出了一种按区域进行填充布点的点云自动生成方法;发展了一种点云的运动技术来实现离散点云对物面边界的随体运动;在点云离散的基础上,采用最小二乘法求解矛盾方程的方法来求取空间导数,进而获得数值通量;采用双时间方法进行时间离散推进,其中物理时间迭代采用二阶隐式格式,伪时间迭代采用四步龙格一库塔显式格式,为了加速收敛,在伪时间迭代中采用了当地时间步长和隐式残值光顺等加速收敛措施。最后,利用本文算法模拟了NACA0012翼型和NACA64A010翼型的跨音速非定常流动,并将计算结果与实验测量结果进行了对比分析,验证了上述方法的正确性和实用性。     

基于DES模型的超燃气动斜坡/燃气发生器方案仿真研究

采用基于剪切应力输运(SST)k-ω两方程湍流模型的分离涡(DES)方法和雷诺时均Navier-Stokes(RANS)方法对超燃气动斜坡结合燃气发生器增进掺混方案进行数值仿真研究。通过对比DES仿真结果、RANS仿真结果与试验结果,发现DES仿真结果对流场内涡结构的捕捉和分辨能力强于RANS方法获得的结果。选取燃气发生器喷流形成的壁面涡流特征角为比较标准, DES仿真结果显示该角度为48°,RANS仿真结果为37°,而试验测得该角度为47°。可见DES仿真结果与试验结果更为接近,说明对于气动斜坡结合燃气发生器的超燃掺混方案,DES方法仿真结果是合理的。     

预测轴流压气机叶栅通道端壁附面层及叶片力亏损

本文将二次流理论与端壁附面层理论相关联,提出一种用于预测轴流压气机叶栅通道内的端壁附面层及叶片力亏损的新方法。附面层主流流动采用包含叶片力亏损的动量积分方程,横流流动采用双层速度分布模型,使用二次流分析和有限差分计算方法预测横流外层速度分布,主流流动和横流流动交替计算。用于两种高负荷压气机叶栅的预测结果表明,不仅端壁附面层发展的预测结果与实验结果吻合较好,而且叶片力亏损发展的预测值也与实验值一致性较好;此外,该方法不仅自动,而且还能较准确地预测S形横向流动的发展。