冲击射流换热数值模拟技术研究概述

介绍了冲击射流的特点,综合评述了对国内外冲击射流换热的试验和数值模拟研究进展.     

动静圆柱坐标系下N-S方程的分区算法及子域干涉处理研究

耦合了相对和绝对圆柱坐标系下三维Navier-Stokes方程,统一了动子域和静子域的控制方程;在子域计算中,采用高分辨率的三阶ENN格式及LU-SGS隐式解法;在子域干涉处理中,研究并对比了动静子域干涉面的重合处理与重叠处理;最终发展了一种新的既高效又准确的分区算法。通过对某带进口导叶的三级跨声速轴流压气机叶栅内三维紊流流场的数值模拟,验证了该分区算法的有效性,计算结果与实验数据吻合较好。     

修正曲率及雷诺正应力对分离流的影响

应用Ｓｔｒａｗｎ＆．Ｋｌｉｎｅ的准同时粘性／无粘流匹配模型，于动量积分布方程中加入了曲率和雷诺正应力两个二阶修正项，计算了二维不可压湍流附面层分离流地。文中首先分析推导了新的曲率修正公式，给出了两个算例，比较了修正项对计算结果的影响，并与实验做了对比。     

马赫数对侧压式高超音速进气道及等直隔离段三维内流场的影响的数值分析

用 N-S方程和 RNG k-ε紊流模型计算了 RBCC用侧压式高超音速进气道三维内流场 ,重点分析了马赫数对流场的影响。从出口截面上气流参数的均匀程度来看 ,Ma3 时比 Ma6时的更均匀 ;从压缩气流的主要部位来看 ,Ma6情况下在支板前肩点之后到喉部截面之前的这个等宽度收敛通道内 ,而 Ma3 时主要是在前肩点之前的收敛通道内 ;从对气流的压缩程度来看 ,Ma6比 Ma3 更大 ;从唇口激波的强度来看 ,Ma6比 Ma3 更强 ;从进气道的流量捕获率来看 ,Ma6比 Ma3 更高 ;高超音速进气道内靠近侧壁以及侧壁和顶板的角区附近是最有可能出现亚音速的区域。     

超燃冲压发动机燃烧室的准一维计算与分析

在考虑有限速率化学反应的准一维Euler方程基础上，通过增加截面面积变化、壁面摩擦和添质的源项，发展了适用于超燃燃烧室性能分析的准一维计算方法。依次以中国空气动力研究与发展中心（CARDC）和日本国家航空与航天实验室（NAL）的氢燃料燃烧室模型作为验证算例，分别采用传统的一步反应模型和发展的有限速率反应模型，模拟了燃烧室流场，并基于NAL燃烧室，计算分析了不同当量比和进口压强对燃烧室流动特性的影响。结果表明：两种方法都能得到与实验数据吻合良好的结果；和一步反应模型相比，有限速率反应模型不仅可以更细致地捕捉流场细节，而且能够初步分析化学非平衡效应的影响；对于NAL燃烧室，当量比≥0.6时，压强随当量比的升高而增大，当达到1.0时，反压已推进隔离段，且推进速度随当量比增大而增加；进口压强不大于110.444kPa时，反压随进口压强增大而升高，且当反压不小于82.833kPa时，反压被隔离在等直段燃烧室入口处；过小的当量比和过大的进口压强均会导致燃烧室出口马赫数严重下降，甚至出现亚声速出流状态。     

两相湍流的非线性k-ε-kp模型

由两相湍流的代数应力模型，提出一种两相湍流的非线性k-ε-kp模型，可以合理地模拟各向异性较强的旋流，浮力流等两相流动，具备二阶矩模型的优点，又比二阶矩模型简单。文中推导出气相、颗粒相的雷诺应力和两相脉动速度关联的非线性应力应变关系式。这些代数式和两相各自的湍动能k，kp以及两相脉动关联湍动能kpg的方程联立，就构成非线性k-ε-kp模型。本文将该模型用于模拟突扩湍流两相流动，给出两相时均速度场及雷诺应力各分量，并且将模拟结果和PDPA实验数据以及二阶矩模型的模拟结果对照。结果表明，该模型预报各向异性两相湍流的能力和二阶矩模型的能力接近，但是计算量比二阶矩模型的小得多。     

跨声速开槽壁超临界翼型洞壁干扰的N—S方程模拟

本文采用二维可压缩非定常Ｎ－Ｓ方程模拟了带有上下开槽壁的跨声速二维管风洞中超临界翼型绕流的流场。网格生成中内层采用解双曲型偏微分方程，外层采用代数生成法。开槽壁采用只有压力梯度没有顶板运动的Ｃｏｕｅｔｔｅ流动的解析解近似。推导出槽间流动与开闭比的三次方关系。上下壁边界层可由无滑移条件直接模拟，也可用平板边界层近似。本文对驻室压力作了几种近似处理，结果表明，所给Ｎ－Ｓ方程的边界条件与风洞实验多接近一     

解决欧拉方程奇异性的方法探讨

提出了一种解决欧拉方程奇异性的工程算法,从飞行器运动方程入手,推导出具体计算公式,并进行了数值计算验证。该方法无论是理论分析还是工程计算都证明是可行的,没有发现原理缺陷和方法误差,不会产生奇异性,在飞行器仿真、模拟和姿态控制中具有普遍的应用和参考价值。     

虚构解方法程序验证

对求解偏微分方程问题的程序的严格验证,一直以来由于其精确解有限,较难进行,针对这种情况,本文给出了一种新的程序验证方法,虚构解方法,该方法旨在解决大型科学计算程序如何较严格地进行程序验证这一问题;该方法通过构造虚构解,修改原控制方程,然后通过对计算结果的分析处理达到对程序的精确验证;文中给出了一组通用的虚构解的选取办法,解决了虚构解方法验证过程中确定虚构解的问题;同时还给出了运用该方法对非结构二维Euler计算程序进行验证的具体例子,网格收敛分析结果表明该方法是一种有效的程序验证方法,并且具有较好的通用性。     

液体火箭发动机推力室冷却通道流动与传热数值研究

采用气固耦合算法对液体火箭发动机推力室再生冷却通道的流动与传热过程进行了三维湍流流动与传热数值模拟，冷却工质为氢气，其密度、导热系数、动力粘度随着温度和压力而变化。应用大涡模拟及标准k-ε双方程模型两种湍流模型分别进行数值模拟，详细揭示了再生冷却通道固体区和流体区内的速度场和温度场，并在不同的计算网格数目下对两种湍流模型的计算结果进行了对比。结果表明，在相同的网格条件下，标准k-ε双方程模型与实验数据的吻合精度比大涡模拟模型更好，且满足工程计算精度。随着网格数的增加，大涡模拟的计算精度逐渐得到改善。