固体火箭发动机内流场分区耦合数值计算

采用分区耦合方法计算固体火箭轴对称燃烧室与喷管流场。对于低速的燃烧室流场，选用不可压流的Ｎ－Ｓ方程描述并用ＳＩＭＰＬＥＣ方法数值求解；对于高Ｒｅ数的喷管流场，则采用Ｅｕ－ｌｅｒ方程描述并用ＳＣＭ方法求解。计算时用燃烧室出流为喷管流场提供入口参数，同时用喷管流场压强分布反馈影响燃烧室流动状况。对耦合边界条件处理方法进行了探讨。对典型的侧壁加质燃烧室与喷管流场进行了计算，计算结果揭示了单独喷管流场计算难以反映的喷管收敛段近壁区的低速区域，与已有的燃烧室流场实验结果一致并反映了燃烧室与喷管流场之间的联系，较好地模拟了流动中的物理现象。     

波瓣喷管引射-混合器的数值研究与验证

对基于CFD数值预测波瓣喷管引射 混合器引射流量比的计算方法进行了探讨,在计算过程中,主流进口采用速度边界条件,二次流进口采用总压压力边界条件,混合流出口采用静压压力边界条件,两者均设置为环境大气压力,与相关实验数据的对比验证表明计算结果与实验结果仅相差10%左右;同时通过改变混合管结构参数,得到了混合管结构参数对引射 混合特性的影响规律,进一步揭示了波瓣喷管有利于强化引射 混合的内在机理。在波瓣喷管出口对应于波谷区域存在一个相对低的静压区,对于引射 混合器系统存在一个最优的性能设计点,计算结果符合物理过程本质。     

气动塞式喷管底部二次流特性的数值模拟

发展了一套可以对底部有二次流的气动塞式喷管的性能进行预示的数值方法，采用二阶精度的NND格式求解二维层流N-S方程，对底部二次流的入口边界条件进行了特殊处理。以一种实验用的直排式气动塞式喷管为对象，针对两种不同的背压条件，对二次流流量分别为主流的0%，1%，3%和5%的几种工况研究了塞式喷管底部二次流流量变化对喷管性能影响。计算结果表明，底部二次流的加入使得气动塞式喷管性能有比较明显的提高。数值方法可以用于气动塞式喷管的设计和性能预报。     

结构参数对两相冲压发动机喷管性能影响分析

为了分析气液两相冲压发动机的喷管参数对推力及效率的影响,分别针对不同参数条件下喷管中泡状气液两相流场,基于双流体模型,采用变步长的Runge-Kutta法进行数值模拟。并对喷管中临界泡状流进行分析,重点研究了喷管喉部与入口面积比、半扩张角和半收缩角在各种气液质量流率比条件下对喷管推力及推进效率的影响,计算结果表明,推力随着喷管喉部与入口面积比增大而减小,气液质量流率比较大时,效率随其增大而增大;半扩张角增大时,推力和效率同时减小;半收缩角对喷管推力及效率影响均不大。在给定的喷管参数下,气液质量流率比在一定程度上增大时,喷管推力及效率同时增大。     

解耦N-S/DSMC方法计算推力器真空羽流的边界条件研究

准确模拟羽流流场是合理评估空间推力器真空羽流效应的基础。因羽流流场同时包含连续和稀薄两种流动机理,通常采用解耦方式的Navier-Stokes(N-S)方程和直接模拟蒙特卡罗(DSMC)混合方法对其进行模拟。为保证解耦N-S/DSMC方法应用于羽流计算时的准确性并尽量提高其计算效率,对计算中的DSMC入口和喷管壁面等边界条件设置问题开展了研究。通过与文献中低总压工况的羽流试验数据比较,确定了合理设置DSMC入口边界位置、壁面反射类型、壁温等边界条件的方法。针对将此方法应用于实际推力器工况计算时效率过低的问题,通过多种数值试验,表明从喷管出口处开始沿KnGL为0.05的等值线作为DSMC入口界面可同时保证仿真精度和较高的计算效率;并证明实际工况下壁温设置对羽流场仿真结果影响不大。     

波瓣喷管引射-混合器的数值研究与验证

对基于Navier-Stokes方程理论预测波瓣喷管引射一混合器引射流量比的计算方法进行了探讨，在计算过程中．主流进口采用速度边界条件，二次流进口采用总压压力边界条件，混合流出口采用静压压力边界条件，两者均设置为环境大气压力，与相关实验数据的对比验证表明计算结果与实验结果仅相差10％左右；同时通过改变混合管结构参数．得到了该参数对引射一混合特性的影响规律，进一步揭示了波瓣喷管有利于强化引射一混合的内在机理，计算结果符合物理过程本质。研究表明，本文的计算方法可以有效地预测引射流量比和揭示混合流场特性。     

波瓣引射-混合器特性参数影响的数值研究和验证

运用不可压缩流动Navier—Stokes方程,对波瓣喷管引射-混合器的流场和引射特性进行了三维数值研究。在计算过程中,主流进口采用质量流量边界条件,二次流进口和混合流出口采用压力边界条件,均设置为环境大气压力。与相关的实验结果对比表明本文的计算方法可以有效地预测引射流量比和混合流场。针对混合管的结构参数开展了系列研究,获得了混合管截面比和长径比对于引射系数和热混合效率的影响趋势。     

一维喷管化学反应流的吉尔方法研究

建立了一维喷管化学反应流计算模型,并进行了数值求解,给出了喷管内一维化学反应流计算结果。采用吉尔方法解决控制方程组的刚性问题获得满意效果。计算中解决了化学反应流计算中初值难于给定和喉部条件难于确定等问题。并通过对氢氧发动机喷管的分析计算,证明了这种方法的正确性。     

喷管碳基喉衬的烧蚀

一、计算方法和结果 喷管材料的烧蚀主要是热和化学作用的结果。热量来自碳与燃气中的氧化元素之间的异相化学反应,以及对流和幅射传热。本文中的解析模型用来预测喷管温度分布及烧蚀率。在处理联结核心流及喷管壁面的边界层流动时,把传热传质数值作为壁内导热计算的输入数据。边界层的边界条件是由一维无粘定常流的解来确定的。在材料响应的分析中介入了幅射边界条件,在UTC数据上引入了一个有效辐射系数。用有限元计算方法求解瞬态热传导方程来确定喷管壁内温度分布随时间的变化关系。     

斜切反喷管性能计算

本文提出了斜切反喷管流场计算方法,内部点用MacCormack二步显格式数值解控制方程,入口边界和固壁边界的参数用物理边界条件和可应用的沿双特征线的相容性方程计算.固壁边界计算中采用了激波装配技术和开关格式,有效地捕获了流场内的激波系列,计算中得到的激波情况和物理分析一致.本文在流场计算的基础上,进行了反喷管的性能计算.本文提出的方法可应用于反喷管的性能预估.     

超声速飞行器侧向喷流干扰流场传统数值模拟方法的误差分析

使用CFD方法,分别就真实喷管边界和简化喷口边界,计算超声速飞行器侧向喷流干扰流场,研究边界条件对干扰流场及气动力的影响.使用k-ε湍流模型封闭雷诺平均N-S方程,利用非结构网格对流场进行空间离散.通过对比,计算结果与实验值吻合良好,证明该方法具有一定可靠性.进一步研究表明喷流边界条件对喷流干扰流场具有一定影响:相对于简化喷口边界,真实喷管边界喷流出口的非均匀性导致喷口上游分离涡和激波位置较为靠前,从而引起附加气动力和力矩的变化;由于摩擦阻力的作用,真实喷管静推力存在损失;喷流压比为500时,总法向力和总俯仰力矩在两种边界条件之间的误差分别为8.21%和22.4%,误差较大.在进行侧向喷流干扰流场的精确计算时,需要考虑边界条件的影响.      

抛物线喷管型面参数对流动分离影响的数值模拟

为了研究抛物线喷管型面的3个参数——初始膨胀圆弧半径、入口角和出口角的变化对喷管流动分离过程的影响规律，模拟了不同参数组合下的抛物线喷管流场参数，获得了喷管流动位置、分离模态随集气室压强增加的变化过程.结果表明:这3个参数能够影响发动机起动过程中喷管流场的发展过程.在一定范围内分离模态转换喷管压比（NPR）与初始膨胀圆弧半径呈正相关；通过合理设定入口角和出口角，可以推迟和缩短受限激波分离过程，为解决彭冠侧向载荷问题提供了可能的方向.      

一种强剪切空气雾化喷嘴的流场和喷雾

采用数值和实验方法对一种强剪切空气雾化喷嘴的流场和喷雾特性进行研究.采用粒子动态分析仪(PDA)实验研究了不同气液比下强剪切空气雾化喷嘴的速度场和喷雾场,并对该喷嘴的流场进行了数值计算,所得的计算结果与实验数据比较吻合.研究结果表明:强剪切空气雾化喷嘴的设计合理,两股旋流空气在喷嘴内部完成强烈的动量交换,并对燃油进行了强剪切.强剪切空气雾化喷嘴在宽调节比范围内能够完成良好的雾化,当气液比变化10倍,液雾的索太尔平均直径分布均处于20~60μm之间.随着气液比的增大,线平均索太尔平均直径逐渐减小,推导了适用于该喷嘴的索太尔平均直径计算模型,可为强剪切空气雾化喷嘴设计提供参考.      

(高)超声速流动试验技术及研究进展

近年来,与高速飞行器相关的(高)超声速流动受到了极大的关注。这类流动所具有的非定常性、强梯度和可压缩性对试验方法和风洞设计技术提出了挑战。超声速纳米示踪平面激光散射(NPLS)技术是由作者所在团队研发的非接触光学测试技术。它能够以较高的空间分辨率来揭示超声速三维流场的一个瞬态剖面的时间解析的流动结构。介绍了NPLS技术以及基于NPLS开发的密度场测量、雷诺应力测量和气动光学波前测量等方法,并回顾了这些技术在超声速边界层、超声速混合层、超声速压缩拐角、激波/边界层相互作用和光学头罩绕流等流动中的应用,清晰地再现了边界层、混合层、激波等典型流场结构及其时空演化特性。另外,为了模拟和研究高空大气条件下边界层自然转捩和超声速混合层的转捩特性,介绍了高超声速静风洞、超-超混合层风洞的设计技术以及层流化喷管的设计方法。     

串装双燃速药柱发动机的内流场计算

对一维的加质燃烧室流场、轴对称的二维喷管内流场和高燃速药柱燃烧完成后的燃烧室二维内流场以及三维药柱几何做了一体化计算。应用轴对称的N－S方程、TVD有限体积法的显式MacCormack格式，对燃烧室后段和长尾喷管跨声速流场进行求解，得到了双燃速内弹道性能和轴对称二维流场中燃气参数的分布。结果具有较高的精度，并为进一步的燃烧室和喷管热结构分析提供数值依据。     

考虑附面层影响的二元混压式进气道设计

采用等激波强度的方法,考虑附面层修正,设计了一种飞行马赫数Ma=3.0的二元混压式进气道.通过数值仿真,模拟了激波-边界层的相互影响,研究了附面层抽吸对内流场的影响,获得了进气道内部复杂的流场分布,以及不同背压下进气道的起动特性.计算表明所设计的进气道性能较好,附面层抽吸对稳定正激波有明显的作用,提高了进气道抗反压能力.给出的方法可用于二元混压式进气道的初步设计和验证.      

离心式喷嘴内流场特性的数值模拟

基于单流体模型及气液两相混合分数的概念，建立了离心式喷嘴内流场的数学模型，并进行了数值计算，计算结果给出了离心式喷嘴内部液膜与气涡共存的流场结构及气液交界面的几何形状，计算得到的壁面压力分布、出口处液膜内轴向速度分布以及雾角大小都与实验结果吻合的很好，与常用的跟踪气液界面的ALE方法相比，本文建立的数值模拟方法具有简单，计算量小和易于实现的特点。     

叶片式预旋喷嘴出口流场实验研究

叶片式预旋喷嘴具有尺寸小,落后角大的特点。为了详细研究小尺寸预旋喷嘴的预旋性能,采用五孔探针对叶片式预旋喷嘴的出口流场进行了实验研究。测量了Ma=0.2,0.3时喷嘴出口的压力分布、速度分布和出口气流角度分布,实验获得了喷嘴的落后角和预旋效率,并进行了与实验工况相同的数值计算。通过实验获得的总压云图以及速度云图,可以发现叶片式预旋喷嘴的端壁二次流损失、尾迹损失严重,有明显的边界层分离现象。Ma=0.2时,喷嘴Re数为5.76×104,落后角2.84°,实验测得的预旋效率为0.73;Ma=0.3时,喷嘴Re数为1.06×105,预旋效率提高至0.77。实验模型端壁的影响使预旋效率实验结果偏低6.5%左右。数值结果与实验测得各参数符合较好:数值结果与测得的喷嘴出口截面平均总压、静压偏差在1%以内;出气速度、周向速度以及出气角度与实验结果偏差在4%以内。数值计算表明,叶片式预旋喷嘴的预旋效率基本不受压比影响,随Re数增大先增大后基本不变,最后基本稳定在0.85。