求解叶轮机械三维粘性可压流的初步研究

给出了用“参数多项式方法”求解叶轮机械三维粘性可压缩流的基本思路、方程和方法，采用圆柱坐标系下的三次参数多项式。中还应用此方法以一个透平和一个压气机叶片为例计算了三维粘性流场。     

粘性跨音速喷管流场计算

本文继用隐式近似因子分解法成功地计算无粘跨音速喷管流场之后,用同样的方法,结合任意曲线坐标系,通过非定常方程在相当长时间后的定常解,针对矩形截面喷管和轴对称喷管粘性跨音速流场,求解了可压缩层流薄层N-S方程,获得了核心流,特别是边界层中的流动参数.对于矩形截面和轴对称两种喷管进行了计算,其结果和实验数据相当一致.将本文的方法应用到粘性两相流动计算,可望较多的节省机时.     

微可压缩模型(SCM)与可压缩N-S方程数值计算对比研究

本文开展了微可压缩模型与可压缩N-S方程之间的对比研究.通过对层流平板边界层和椭球高雷诺数湍流绕流的数值模拟,比较了微可压缩模型和可压缩N-S方程之间的计算精度和效率.计算结果表明,微可压缩模型和可压缩N-S方程具有相同的计算精度,但是微可压缩模型具有更高的计算效率,内存和CPU时间都节省20%左右.     

再生冷却推力室耦合传热数值模拟

对大推力液体火箭发动机再生冷却推力室内部燃气、室壁和再生冷却剂进行了耦合传热数值计算.采用二维轴对称N-S方程描述推力室内部燃气的湍流流动与传热,对冷却剂流动采用简化的一维模型,通过室壁的偶合传热采用一维传热模型.N-S方程的求解采用贴体坐标系下的有限容积法,速度和压强的耦合采用可压缩的SIMPLE算法.湍流的模拟采用可压缩的标准κ-ε模型,辐射传热采用热流模型计算.研究表明,本文方法可较好地模拟燃气的二维流动,同时能快速计算壁面热流密度、壁温和冷却套温升,计算结果对大推力发动机推力室的设计具有一定的指导意义.     

可压流的龙格-库塔时间步进求解法及其应用

研究了用经典四阶龙格－库塔法计算可压缩流场的可行性,并将其用到了电弧喷射推力器内部等离子体流场的数值求解中。应用情况表明,在结合了局域时间步长、隐式残值光滑加速收敛措施后,本格式能够成功地计算比较复杂的可压缩气体流场以及等离子体流场。有关计算结果揭示了气体流经电弧喷射推力器通道但无电流时形成的纯气动流场以及有电流通过时形成的等离子体流场的丰富的结构和一些重要的影响因素,为研究其过程机制提供了依据。     

粘性不可压缩流动三维复杂流场分块耦合求解

提出了一个粘性不可压缩流动三维复杂流场的分块耦合求解方法。可解决科学和工程计算中的各种稳态和瞬态粘性不可压缩流动计算问题。     

径向开槽药柱发动机内流场数值模拟

采用时间相关的显示ＭａｃＣｏｒｍａｃｋ预估－校正两步格式和可压缩全Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，计算了轴对称径向开槽药柱发动机内流场。计算中采用了贴体坐标系统和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型。边界参数计算采用参考平面上的特征线法和物理边界条件。通过大量算例的分析表明，方法有效，计算结果合理。当有集中质量加入时，燃烧室头部压强升高，喷射强度和开槽口的倾斜方向对头部压强也有较大的影响。     

用于粘性可压缩流动数值计算的SIMPLE方法

对SIMPLE方法加以推广,使之适用于可压缩粘性流动的数值计算。采用非交错网格技术,推出三维任意曲线坐标系下可压缩形式SIMPLE方法的计算方程。对亚音速、跨音速及超音速等5个流场进行了数值计算并与有关文献及实验数据进行比较。      

弱可压缩流体与可压缩流体模型的管路水击研究

为提升液体推进系统模型对可压缩流体的仿真能力,针对供应管路的瞬变流动,采用一维有限体积法,使用Modelica语言基于MWorks平台分别开发了适用于补燃循环发动机仿真模型库的弱可压缩流体与可压缩流体的管路模型,不同于控制方程中显含声速的弱可压缩流体模型,可压缩流体模型通过将守恒方程与真实流体的SRK型状态方程耦合求解来计算压力波的传播。在管路供应系统中,分别以液氢、液氧及液甲烷为工质进行关阀水击的仿真研究,结果表明:使用SRK状态方程避免了可压缩流体模型对内能拟合公式的依赖,减轻了建模难度;可压缩流体模型对液氧及液甲烷的仿真结果较为准确,最大误差不超过2%,对液氢的最大计算误差约9%;相对于可压缩流体模型,弱可压缩流体模型计算得到的水击压力在振幅及频率上均偏大;可压缩流体模型计算得到的流体密度波动不超过3%,因此工程实践中处理水击问题时将液体假设为弱可压具有一定的合理性。     

可压缩多介质流动问题的高精度数值模拟方法

在可压缩流动问题的数值模拟领域,激波的高分辨率计算已取得重要进展.但是包含物质界面的可压缩多介质流动的数值模拟还存在诸多数值挑战,主要表现为界面处数值耗散过大和非物理振荡等问题.界面处流体性质的不连续性是造成可压缩多介质流动问题物理建模与数值方法困难的主要原因.为了建立一套高效的可压缩多介质流动问题的高精度数值模拟方案...     

发动机进气流动三维瞬态数值模拟研究

采用任意拉格朗日-欧拉法(ALE法)和所提出的动态网格生成新方法,实现了车用发动机进气流动过程的三维瞬态可压缩湍流数值模拟,湍流模型采用经压缩修正的κ-ε双方程模型.计算结果揭示了发动机进气过程的三维瞬态流动特性,为发动机进气系统的优化设计提供了理论依据.     

有限分析方法求解粘性可压缩紊流

本文将有限分析方法进一步推广到粘性可压缩紊流流场的数值研究, 在贴体座标网格上数值分析了单列和串列叶栅内粘性可压缩紊流流场。计算表明, 有限分析方法具有数值精度高, 稳定性好和收敛速度快等优点, 适于粘性可压缩紊流流场数值计算。      

不可压缩粘弹性有限元法及其对推进剂药柱应力分析的应用

由于固体推进剂具有近似不可压缩的特点 ,利用位移法有限元求解这类材料的力学问题不可能得到满意的结果。本文根据 Herrmann变分原理导出一种适用于不可压缩和近似不可压缩粘弹性物质的新的本构关系式 ,介绍了相应有限元公式 ,并开发了用于平面问题分析的计算机程序 VFAPINP。计算结果表明 ,本方法和程序可以用于所有泊松比的粘弹性问题的计算 ,尤其适用于固体推进剂药柱的应力分析。     

固体火箭发动机喷管气固两相流动的数值模拟

对颗粒相采用颗粒轨道模型,气相求解可压缩N-S方程组,计算方法采用显式Runge-Kutta时间推进法与有总变差衰减(TVD)性质的高精度MUSCL-Roe格式;自主开发了曲线坐标系下二维轴对称可压缩N-S方程组的解算器Solve2D,研究了固体火箭发动机喷管中颗粒相对流场的影响以及不同尺寸颗粒运动规律.结果表明:颗粒相对流场的影响主要表现在喷管喉部以及扩张段,和单相流场相比,沿轴线马赫数减小,且颗粒尺寸越小减少得越多;沿轴线气相温度升高,且颗粒尺寸越小温度升高越多;颗粒尺寸越小,无粒子区越小;颗粒越大与收缩段壁面碰撞越剧烈,无粒子区越大.      

固体推进剂药柱松弛模量随机粘弹性有限元分析

为了研究性能参数的分布对药柱结构完整性的影响，以不可压缩粘弹增量有限元和摄动法为基础，采用局部平均方法对随机场进行离散，发展了一种松弛模量随机粘弹有限元方法，讨论了模量的随机性对结构分析的影响。对药柱进行随机模拟，程序实施简单，计算效率和精度较高。     

计算流体动力学及其发展

本文介绍计算流体动力学的四个发展阶段和数值求解 Navier-Stokes 方程组的基本差分方法及困难;较详细地叙述了不可压缩流动的逐次逼近法、稳定化法和 MAC 法,可压缩流动的 Lax-wendroff法和 ICE 法。     

三维机翼可压缩层湍流边界层计算

 从三维时间相依可压缩边界层动量积分方程和平均流动能积分方程出发,将二维可压缩层、湍流边界层积分方程算法推广到求解有限翼展后掠机翼的边界层流动。利用四步Runge-Kutta时间步进方案数值求解积分方程组,并利用当地时间步长加速迭代收敛。文中分析了数值方法的稳定性与收敛性,并考查了横向耗散项对计算结果的影响。算例表明,能获得令人满意的三维机翼定常可压缩层、湍流边界层的计算结果。     

用BGK格式对平面可压缩自由混合剪切流动的数值模拟

本文介绍了ＢＧＫ－Ｂｏｌｔｚｍａｎｎ格式的构造与用有限体积法计算网络内宏观量的方法。文中还介绍了可压缩平面混合剪切层计算域在不同出口马赫数下的边界条件。然后将上述方法应用于四种不同情况的可压缩平面混合剪切层的数值模拟。还详细介绍了计算结果,并与其他方法所得的的结果以及实验数据进行比较。表明本文的结果较好,方法是有效的。     

计算流体动力学及其发展

本文介绍了计算流体动力学的四个发展阶段和数值求介Navier-Stokes方程组的基本差分方法及困难;较详细地叙述了不可压缩流动的逐次逼近法、稳定化法和MAC法,可压缩流动的Lax-Wendroff法和ICE法。     

涡扇发动机低转速部件特性扩展和风车状态性能模拟

借助涡扇发动机风车特性研究的最新成果,开展了某型涡扇发动机的风车性能数值模拟研究.首先,基于不可压流理论,附以可压缩效应修正的指数法发展了补充慢车转速以下部件特性的计算模型和计算程序;其次,根据流量连续、压力平衡、功率平衡、燃烧室进出口总温相等及转速相等等约束条件建立了风车状态时涡扇发动机高低压转子及整机的共同工作方程;最后发展了具有一定计算精度的涡扇发动机风车状态特性计算模型和程序,并以某型涡扇发动机为例,计算分析了不同飞行条件下发动机的风车特性,从计算结果看,其趋势是合理的.