应用于轴对称流的流图变换中的两个问题

本文导得了轴对称流的位函数φ与流函数ψ的偏导数之间的一个关系式,此式与平面流的关系式是类似的,它是Molenbroek-Chaplygin变换式的推广。本文还导得了轴对称流中的导数的表达式。该式虽然与平面流的形式不一样,但也可以统一起来,它也是平面特征线关系式的推广。     

轴对称羽流流场数值模拟和实验验证

针对在轨航天器羽流效应对航天器整体寿命的影响,在"DFH"系列卫星姿控发动机喷管羽流实验研究的基础上,建立喷管的轴对称模型,选取矩形计算区域,设定一定的边界条件和计算条件,采用轴对称单组分直接模拟蒙特卡罗(DSMC,Direct Simulation Monte Carlo)方法对相应实验状态的羽流流场进行数值模拟研究.通过对比数值模拟结果和实验测量结果,表明DSMC数值模拟羽流流场的计算结果符合实际流动规律,实验验证了所编制的轴对称DSMC数值模拟程序的正确性和有效性.      

绕大钝头短体外形跨音速大迎角湍流流动计算

采用数值模拟求解N-S方程方法计算与分析了绕大钝头短体外形的跨音速大迎角时的湍流流动.数值方法中选用Van Leer 分裂格式离散无粘通量项,并构造了一种Limiter函数以保证TVD(Total Variation Diminishing)性质,采用中心差分格式来离散粘性通量项.计算中采用了B-L湍流模型.算例定性地表明绕流的三维分离具有准周期性的非定常特性.     

进气旋流对推力喷管性影响的研究

介绍了模拟近代大涵道比发动机排气状况的旋流模型,和在喷管落压比1.2到6.5范围内,进气旋流对轴对称喷管,矩形收-扩喷管,楔形二元喷管,和单边膨胀斜板形喷管等性能和射流流态影响等的一些实验结果。实验数据指出:进气旋流可使喷管出口射流柱中的旋涡加强,射流与外界气流间的混合区扩大,中心核心流区缩小,但它对轴对称喷管出口流态的影响很小。进气旋流同时对喷管推力系数和流量系数均起不利影响,增大进气旋流角度会使喷管推力系数和流量都下降,特别是当旋流角超过某一“临界”值后二者降低很急剧,并且二元喷管降低的程度比轴对称喷管更严重,因之在具体应用情况中,保留多大发动机的旋流,要针对具体情况综合考虑各种因素来选定。     

轴对称跨音速流中无激波超音速区的设计

本文直接从音速线出发计算了一个旋成体跨音速绕流超音速泡中的特征线网。提出了从音速线开拓到一条超音速线的方法,计算了无旋轴对称流的特征线网,发现其特点是不能跨越对称轴线,并对其原因作了初步的讨论。本文的计算结果可应用于物面修形。     

气动塞式喷管实验与数值模拟研究

提出了塞锥型面设计的特征线法及示例结果,得出了侧喷管粘性跨音速和超音速内流场数值计算结果,及塞锥外流和尾流的Euler方程解,总结了塞式喷管工作特性的固体推进剂燃气模拟实验系统、实验结果及主要结论.塞式喷管火箭发动机是一个复杂而有挑战性的研究方向,对于发展新一代先进天地往返运输系统有重要作用.      

使用Gao-Yong湍流方程计算翼体角偶流动

Gao-Yong湍流方程基于侧偏统计平均方法保留了湍流脉动量的一阶统计信息,并引入加权漂移速度对称性及正交各向异性,导出漂移流连续方程、动量方程和机械能方程,最后依据湍流物理的唯象论,使用动量传输链概念模化封闭了整个方程组.大量算例验证了此方程对广泛范围的复杂湍流问题的适定性.为了进一步的工程应用,将其加入开源计算流体力学软件OpenFOAM(Open Field Operation and Manipulation)中,并就翼体角偶流动进行了数值模拟.翼体角偶流动具有三维分离、马蹄涡等十分复杂的流动特征,通过与实验测量结果的对比研究表明,Gao-Yong湍流方程不论在定性还是定量上,都成功捕捉到这类复杂流动现象的主要特征,得到了令人满意的结果.     

引射混合器内流场的数值分析

计算了以波瓣为主喷管的三维引射掺混流场；分析了包括复杂型面和引射二次流在内的边界条件处理；再现了流向涡的产生，计算结果与实验结果作了对比分析，吻合程度较好。     

超音速喷流DPIV瞬时速度场实验测量

介绍了采用数字式互相关粒子图像测速系统(DPIV)在测量超音速喷流实验中的应用;实验给出了设计马赫数Ma=1.5的小型拉瓦尔喷管在不同总压和反压比条件下,喷流速度场、流线、涡量分布等定量信息.实验结果显示出流场中激波前后流体速度,涡量分布的明显变化,波系的结构和不同于一维管流的流动特性.     

火焰筒头部流场研究之二——二元轴对称等温廻流区计算

本文将火焰筒头部的等温流场,用二元轴对称不可压定常流来处理并用数值方法解出初步结果。基本方法是用涡量流函数法,双方程粘度模型,采用交替方向隐含格式差分求解,在边界及迭代初值提出自己的看法,初步计算结果与[10]之实验定性一致。     

多段翼型贴体网格的生成方法

本文给出了一种多段翼型绕流的分区贴体网格生成方法。采用泊松方程的解来给出贴体网格的坐标,泊松方程中的源项通过靠近边界的网格间距和正交性来决定,本文对分区网格的区域边界的选择作了讨论,采用生成的分区贴体网格计算了绕多段翼型位流的解,其结果与实验所得是很一致的。     

三维抛物化N-S方程全隐格式矢通量分裂求解

采用全隐格式的矢通量分裂算法,在通用曲线坐标系内求解三维无因次抛物化强守恒型N-S(Navier-Stokes)方程.所建立的计算程序,可计算较为复杂边界的三维流场,如求解亚、跨、超音速共存的流场,分区模拟内外流共存流动及其混合剪切尾流.给出了一轴对称亚跨超音速喷管内流场的计算结果,包括计算网格,流矢分布,压力、温度和马赫数等值线分布.也给出了一喷管内外流及其混合尾流的计算结果,包括计算网格,流矢分布,无因次压力、密度和温度等值线分布.     

矢量喷流对战斗机绕流和气动特性影响的研究

在双三角翼、椭圆机身的冀身融合体气动布局中配置二元矢量喷管,本文对该布局形式的矢量喷流和主流绕流之间的干扰流场及其对机翼气动特性的影响进行了系统的实验研究。结果表明,矢量喷流对主流及其气动特性的干扰效应可以分为三类:矢量喷流对附流的干扰效应仅限于喷口附近,主要表现为引射或阻塞作用;矢量喷流对前缘涡的干扰效应起稳定旋涡的作用,这样对稳定涡系的干扰影响很弱,但对破裂涡绕流的干扰效应较强,且随喷流矢量角的增大而增强。     

固体运载火箭变轨发动机喷管气流分离研究

固体运载火箭变轨发动机喷管在工作过程中可能产生气流分离问题,为研究气流分离对喷管性能的影响,开展了理论计算与数值模拟分析。通过分析获得了气流分离点位置、推力系数、喷管壁面的压强、对流换热系数、温度分布。结果表明:地面推力系数是真空推力系数的73.3%,喷管气流分离影响了发动机能量转换;气流分离后喷管壁面压强、对流换热系数、温度存在跃变现象,从而会对喷管扩张段产生不利影响。该分析为进一步研究固体火箭发动机高空喷管通过地面试验性能预示高空性能及喷管扩张段热防护设计提供参考。     

战术导弹横向喷流数值模拟

为了研究喷流对导弹不同升力面的干扰影响,通过数值求解N-S方程来模拟超声速外流场中横向喷流的干扰流场,采用分块对接网格和"O"型网格技术,精确模拟喷口截面及弹翼形状,生成高质量的贴体计算网格.通过几种升力面和弹身组合的超声速横向喷流干扰流场的数值模拟,计算的马赫数为2和4.5,计算结果显示具有升力面的喷流干扰比光弹身的喷流干扰大.研究和分析了喷口附近流场的涡系结构和波系结构,并将导弹喷流干扰力放大因子及喷流引起的压心变化的计算结果与目前国外公布的试验值进行比较分析,二者显示了良好的一致性.     

喷管两相流耦合数值模拟与粒子热增量预测

含有氧化铝粒子的两相流是固体火箭发动机喷管流场的重要特征.在有限体积方法框架下,采用基于热增量试验数据的粒子壁面反弹模型以及基于粒子轨道的单元内颗粒源(PSIC,Particle Source in Cell)两相流耦合算法,对喷管内两相流流场及粒子撞击产生的壁面热增量进行了计算和分析,研究了氧化铝粒子尺寸对粒子轨道分布和喷管壁面热增量分布的影响规律.研究结果表明:喷管扩张段内粒子稀疏区域范围随粒子直径增加而增大;粒子热增量只分布于喷管收缩段内,粒子直径越大,产生的壁面热增量越强.     

带副翼的翼身组合体绕流的 Euler方程解

提出了一种求解带副翼偏转的翼身组合体绕流的Euler 方程计算方法.将对接分区网格与分区求解算法相结合,有效地求解了绕此外形的复杂流动.提出了一种满足通量守恒的内边界耦合条件.数值方法中选用Van Leer 分裂格式离散无粘通量项,并构造了一种限制器(Limiter)函数以保证TVD性质.数值算例表明本文方法是求解带操纵面偏转的翼身组合体绕流的有效方法.     

直升机发动机入口绕流的数值模拟

对直升机发动机入口处流场采用有限体积法,通过生成贴体结构化网格求解Navier-Stokes(N-S)方程, 进而对MI-171直升机低空悬停和前飞状态进行数值模拟,计算了该型直升机发动机进气道处的流场特性.然后以广义尾流入手,并引入尖涡涡核作用的半经验修正的旋翼自由尾流模型对该机发动机进气道处流场进行计算分析.对沿发动机轴线剖面以及垂至于该轴线剖面的分析结果表明,CFD(Computational Fluid Dynamics)数值模拟方法与旋翼自由尾流分析方法均可以对涡诱导流场进行预测.     

一种跨声速定常流场求解加速方法

跨声速定常流场的隐式求解相当于使用牛顿迭代法求解一个非线性方程组。为满足牛顿迭代收敛性的要求,通常需要对所求解问题进行全局化处理。在同伦延拓的框架内,提出了一种基于拉普拉斯算子的方程延拓方法,提高了定常流场隐式求解收敛速度。针对定常流场通常初始化为均匀来流的特点,一方面利用拉普拉斯算子的椭圆性加快边界条件信息向流场内部的传播,另一方面利用拉普拉斯算子的线性和正定性改善延拓问题的正则性,综合两者增加拟牛顿算法的稳定性,提高可用CFL数,最终达到提高流场求解效率的目的。由于流场问题的复杂性和非线性,难以通过理论分析得出先验的最优非线性求解策略。因此,通过无黏NACA0012翼型、湍流RAE2822翼型和三维ONERA M6机翼等算例的数值实验,研究了拉普拉斯项参数对收敛效率的影响,给出了效率较优的参数组合,验证了本文方法在跨声速情况下相对于经典伪时间推进法可以节约20%以上的CPU计算时间。