高精度ENO格式在射流数值模拟中的应用

本文分析了ENO格式的特点，并应用于全NS方程的迁移项和压力顶，模拟了轴对称和三维射流。首先对轴对称自由射流为算例比较了高阶TVD格式和ENO格式，证实了ENO格式捕获复杂激波波系的能力；其次计算了四喷管湍射流干扰流场，给出了相应的计算结果；最后计算了四股射流流入圆管内的流场，给出了物理量沿壁面和截面的分布。     

增强型双喉道射流推力矢量喷管的流动特性试验

对一种增强型双喉道射流推力矢量喷管开展了内部流动特性的试验研究,获得了其在不同次流压比状态下的内流结构和沿程静压分布.试验结果显示:在基准双喉道矢量喷管尾部附加扩张段后,能够以2.8%的次流消耗率获得超过20°的平均气流偏角,这表明通过附加扩张段来增加喷管矢量角的设计概念是可行的.在凹腔内,增强型双喉道射流推力矢量喷管的静压分布规律与基准双喉道矢量喷管一致,但在附加的扩张段内,下壁面的压强要明显高于上壁面,这正是其推力矢量角得到显著增大的原因.随着次流压比的增加,喷管获得的推力矢量角单调增加,但是喷管附加扩张段的矢量增强效果基本维持不变.      

逆流矢量喷管流动自身导致的非定常现象

基于数值方法,研究了逆流推力矢量喷管中由于流动自身导致的非定常现象。结果表明,当抽吸背压下降至某个范围时(50.66~60.8 kPa),主流在接近外套壁面的某个位置处呈现出周期性振荡,流动参数(抽吸二次流量、推力矢量角、合成推力系数)也以小振幅振荡,但是同向二次流量却没有振荡。深入的流动机理分析表明,此时在逆流剪切层中,分离涡周期性地形成和发展,没有呈现出流动失稳,同时剪切层中的速度比又远大于定常流动下的值,因此流动表现为自激振荡特征。     

次流喷射控制推力矢量喷管的流场数值模拟

为了探讨次流喷射控制二维推力矢量喷管的性能,用时间推进求解Ｎ－Ｓ方程的方法数值模拟了喷管的流场,数值格式用中心有限体积格式,借助于当地时间步长和隐式残差光顺技术加速收敛,得到了在不同二次喷流压力下推力矢量喷管内流场的变化特征。计算结果与实验结果比较有满意的一致性。研究表明,应用次流喷射控制主流流动可以实现较大的推力矢量转折,但是,二次喷流必须具有足够的压力值。     

火箭喷管推力偏心高精度计算软件

利用火箭喷管从亚声速、跨声速到超声速的非对称流动数学模型建立了计算喷管推力偏心的理论模型，并采用高精度有限体积TVD格式求解数值方程，可以对火箭发动机喷管的推力偏心特性进行评估与设计，并集成为Nvc软件。该软件完成了从喷管网格生成、数值计算到推力偏心分析全过程，界面友好，使用方便，大大提高了设计效率，通过大量的计算与试验结果的比较，计算结果可以满足工需要，减少试验次数，节约成本，具有广泛的用途。     

一种矢量增强型双喉道射流推力矢量喷管的数值模拟

为克服双喉道射流矢量喷管矢量角偏小的缺点,提出了一种矢量增强型双喉道矢量喷管的设计概念:在喷管尾部增加一扩张段,利用流体的附壁效应使主流在扩张段中进一步偏转,从而获取更大的矢量角.首先对设计概念的可行性进行了仿真分析,而后对扩张段的设计规律进行了研究.结果表明,在喷管尾部附加扩张段可显著强化其推力矢量性能,使矢量角达到20°以上,但也导致了一定的推力损失.在研究范围内,扩张段扩张角、扩张段长度、扩张段型线等设计参数对喷管的矢量效率、推力系数以及内部流态均有着显著影响,而在扩张段开缝则可以作为一种抑制尾喷流过膨胀的有效措施.若将内凹型扩张段与开缝方案相结合,仅需消耗2.8%的次流便可获得24.12°的推力矢量角和0.929的推力系数.      

斜切反喷管粘性流场数值模拟

为了在级间分离期间提供反推力，许多固体火箭发动机前端都装有一组斜切反管，由于反管气动型面具有尖点，并在超音速区有台阶，喷管内存在一系列激波，并伴有流动分离现象。从雷诺平均的非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程出发，利用时间相关法，采用ＭａｃＣｏｒｍａｃｋ两步显格式，结合Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型，数值模拟了斜切反喷管流场。     

WENO格式在隐态问题中的应用

高阶、高分辨率的WENO格式大大加强了对于复杂流场结构的分辨能力和计算精度，但计算时间较TVD格式有了明显的增加。近似隐式分裂的LU方法能够加快求解稳态问题的收敛速度。本文结合WENO格式和LU方法实现了求解稳态问题的有效方法，它使得求解稳态问题在收敛速度和计算精度上都有较大的提高。通过算例的数值分析可以看出，上述方法不但收敛速度有了明显的提高，而且对于复杂流场现象仍然具有良好的分辨能力。     

二元双喉道射流推力矢量喷管流动参数影响的数值研究

采用数值模拟方法研究了不同流动参数对二元双喉道射流推力矢量喷管(Dual-throat fluidic Thrust-vectoring Nozzle,DTN)内流特性和推力矢量控制效果的影响。结果表明,DTN在非推力矢量情况下,NPR在3~4范围时,推力系数较大,达到0.968,而流量系数较小,仅为0.93;NPR再增大,推力系数迅速下降。在推力矢量情况下,落压比一定时,随着次流流量比的增加,推力矢量角增加,而流量系数、推力系数、推力矢量效率减小;次流流量比一定时,随着落压比的增加,推力矢量角减小,系统推力系数先增加后减小,流量系数略微增加。     

轴对称矢量喷管结构的计算机模拟及运行分析

采用工作站与微机混合建模的方式,利用基于微机的动画及图像处理软件,对轴对称矢量喷管热态试验件的真实运动过程的进行计算机三维运动模拟,验证运动机构方案的合理性,可行性,结构尺寸,连接方式,型面设计的正确性,分析运动干涉情况,确定主要构件的运行轨迹,为进一步的工程设计提供了依据。     

基于非结构嵌套网格方法的旋翼地面效应数值模拟

 建立了一个基于非结构嵌套网格的流场求解器,用来精确模拟复杂的旋翼近地流场,为更好地分析地面效应(IGE)对旋翼气动特性的影响提供一套计算方法。在该求解器中,控制方程采用惯性坐标系下的非定常N-S方程,空间方向上采用二阶迎风格式,并用背景网格的一个面模拟地面作用,以方便地面边界条件的处理。应用所建立的模型,首先针对有实验结果可供对比的旋翼无地面效应（OGE）流场进行了数值模拟,以验证计算方法。然后着重计算了IGE下流场中的桨尖涡空间位置和旋翼拉力增益,并对旋翼在小速度前飞状态下的近地流场及地面涡形成过程进行了模拟和分析。在此基础上,得出了一些有意义的结论。     

11组元轴对称热化学非平衡流场的数值研究

本文从轴对称热化学非平衡Ｎ－Ｓ方程出发,采用双温度、十一组元反应气体模型,利用隐式ＮＮＤ有限差分格式的时间预处理技术数值求解了ＦＩＲＥＩＩ飞船热化学非平衡流场。文中同时给出了七组元模型的计算结果,并进行了结果比较。     

单边膨胀矢量喷管气动和红外特性研究

在实验数据验证的基础上,通过计算流体力学/红外辐射(CFD/IR)综合的方法,对不同落压比、不同几何矢量角下的单边膨胀喷管(SERN)进行分析。研究结果表明:喷管无几何矢量动作下,低落压比下的单膨胀边上过度膨胀是造成喷管推力性能急剧下降的原因;喷管在负矢量角下,过度膨胀加剧,推力性能降低;随着喷管几何矢量角绝对值的增加,矢量推力增加,但推力系数减小,喷管几何矢量角在±25°、喷管落压比在3~6的研究范围内,喷管推力系数最低为0.88左右,最高达0.98;喷管几何矢量角为5°时,喷流红外辐射强度最大,喷管矢量角偏离5°的程度越大,尾焰红外辐射强度越低,但是空间分布规律不变。随着喷管几何矢量角的改变,喷管整体红外辐射强度的空间分布规律发生改变,几何矢量角为负时,辐射强度值大的探测角度向下方移动,几何矢量角为正时,喷管整体红外辐射较强的位置分布在上方,由单膨胀边高温壁面以及喷管内腔的可视面积决定。     

固体火箭发动机流体喉部推力矢量特性

针对采用水作为二次流工质的流体喉部进行了冷流实验及数值模拟研究.研究了该种固体火箭发动机流体喉部的一般规律,包括不同二次流射流方式,不同二次流流量下流体喉部的扼流性能,推力偏角及推力效率,数值模拟及实验结果吻合较好.结果表明:扼流性能与二次流的注射位置、注射角度及流量比有关,且随二次流/主流流量比的增大而增大.喉部二次流喷射能有效的调节有效喉部面积进而调节推力大小,当流量比为0.4时,最大有效喉部面积比为0.8;扩张段二次流喷射能有效调节推力方向,当流量比为0.4时,最大推力偏角为20°;喉部二次流与扩张段二次流入射位置存在相位差可有效降低喉部与扩张段二次流干扰.      

球面收敛二元扩张矢量喷管热射流特征的实验研究

为了获得球面收敛二元扩张矢量喷管热射流特征，在小型热射流实验台架上运用红外热像显示和流场测量，在总温673～823K、喷管落压比(NPR)1．46～2．25状态下，对不同喉道宽高比、俯仰角和偏航角的热射流特征进行了实验研究。在本研究范围内，NPR对矢量喷管上壁面附面层流动影响不大；在扩张上壁附面层分离的情况下．壁面附近的流线被上抬，形成“马鞍”形异型总压分布；喷121出121下游的激波呈现一连串的“葫芦”状，随喷管矢量角的增大，流场的不对称性越明显；在收扩喷管喉道宽高比（ARi）为2时，热射流被分成两股。     

叶轮机械三维不可压Euler方程的直接解法

本文用近似因式分解交替方向隐式格式直接求解叶轮机械中原参数三维不可压Euler方程组。在连续方程中引入“拟压缩性”之后,将原来类型不确定的控制方程变成双曲型方程,给定初边值之后用时间推进法求得稳定解,对一个单级压气机转子内部流场进行了计算,结果与实验值比较符合。     

ENO-AUSMPW格式对激波涡干扰的数值模拟

高阶ENO格式在复杂流场计算中具有重要的地位,它理论上可以达到任何阶的精度。本文采用二维高阶ENO插值方法,构造了高阶AUSMPW格式,并对包含复杂激波结构的激波涡干扰流场进行了数值模拟。计算结果表明了高阶ENO插值下的AUSMPW格式具有较高的激波分辨率和较低的数值耗散。     

悬停状态倾转旋翼/机翼干扰流场及气动力的CFD计算

基于一套高效通用的多层运动嵌套网格技术,建立了适合悬停状态倾转旋翼机旋翼/机翼气动干扰特性分析的高效混合计算流体力学（CFD）方法。在倾转旋翼/机翼贴体网格区采用可压雷诺平均Navier-Stokes（RANS）方程作为主控方程,过渡/背景网格区选用Euler方程,湍流模型选用Spalart-Allmaras模型。时间推进上采用双时间推进格式进行非定常求解,并在方法中运用了SPMD（Single Program Multiple Data）模式的并行加速技术。在此基础上,首先分别采用UH-60A直升机旋翼及XV-15倾转旋翼机旋翼作为数值算例,验证了CFD方法的有效性。然后着重对倾转旋翼/机翼的非定常干扰流场及气动力分布特征进行了数值研究,模拟得到与实际情况相符的“喷泉效应”干扰现象。计算结果表明,干扰作用使得倾转旋翼相对于孤立旋翼拉力数值减小了3%,但总的拉力系数损失达到了17%,证明悬停状态下气动干扰对飞行器气动性能有重要影响。     

矢量推力喷流对飞行器气动载荷影响的研究

矢量推力喷流对飞行器绕流具有干扰作用,引起飞行器气动载荷变化,对飞行器操纵性也有潜在影响。这种气动干扰和载荷变化是矢量推力技术研究的一个重要方面。与普通喷流相比,矢量推力喷流具有喷流偏转的特点,喷流偏转对气动干扰有重要贡献。因此本文采用多块搭接网格模拟了矢量推力飞行器无粘绕流流场,通过力系数和力矩系数对比,探讨了飞行器气动载荷随喷流偏转角的变化规律,这些规律具有理论和实际意义。