脉冲爆震发动机快起爆的二维数值模拟

为研究脉冲爆震发动机的点火过程，用二维欧拉方程和Korobeinikov两步爆震模型对爆震波快起爆过程进行了数值模拟研究。用高温高压未燃气体点火方式来模拟试验中的电火花塞点火。研究表明：电火花塞的能量不足以直接点燃爆震波，爆震波是在经历了一系列反射激波的相互作用后最终建立的。激波．爆震波转捩是快点火过程，是爆震波建立的一个非常重要的机理。在设计点火装置时，火花塞要尽可能靠近壁面和拐角，以便产生高强度的反射激波。同时还要控制反射激波的方向，使之向未燃气体中顺利传播，才能起到很好的点火效果。     

热点触发爆震波的数值研究

运用时空守恒元和解元(CE/SE)法模拟爆震波由热点触发的过程,分析通过热点触发爆震波的机理以及热点触发爆震的影响因素。数值研究表明,热点内部均匀温度分布和线性温度分布有着不同的加速机理。触发爆震波过程受到热点区域大小、形状及温度分布等因素的影响,半球形的热点形状最易起爆;热点区域越大,温度分布线性化越高,越容易形成爆震。     

二维流场、热、结构一体化数值模拟

流场、热、结构问题,在传统的计算中,通常是先算流场、再算热响应和热结构问题,这显然是人为地分割了本身统一的、连续的过程,使得计算结果存在较多不足。本文给出一种把流场、热、结构耦合起来进行一体化数值模拟的方法。算例证明,本方法较好地模拟了动态的物理过程。     

旋转爆震波发动机二维数值模拟

采用一种改进的化学非平衡流解耦方法对Euler反应流方程解耦处理,对流项采用五阶WENO格式离散,化学反应源项的刚性采用简化的隐式方法处理,时间步进采用二阶精度的Runge-Kutta方法,对H2/Air预混气旋转爆震发动机内流场进行了二维数值模拟。模拟结果给出了不同发动机尺寸条件下的详细起爆过程,结果表明当发动机尺寸小于临界直径时无法成功起爆;详细分析了流场结构和爆震波形状,旋转爆震波的传播速度与理论预测值吻合;性能分析结果表明在喷注总压低于燃烧室平均压力时仍可实现推进剂喷注,没装尾喷管的情况下发动机比冲达到176.5s。     

旋转冲压压缩转子二维进气流道数值研究

借鉴二维超声速进气道的设计方法,设计了一种旋转冲压压缩转子的二维进气流道,并采用二维雷诺平均N-S方程和Spalart-Allmaras湍流模型对其流场进行了数值仿真,研究了转速、背压对二维进气流道中波系结构、内部流动特性和性能的影响.计算结果表明:所提出的旋转冲压压缩转子二维进气流道设计方法是可行的;二维进气流道中产生的激波系与二维超声速进气道中产生的激波系相类似,所不同的是二维进气流道中产生的激波为弯曲激波,而二维超声速进气道中产生的激波为平直激波;转速和背压对二维进气流道的性能有较大影响.      

旋转爆震发动机燃烧室气膜出流影响数值研究

针对旋转爆震发动机（RDE）壁面的高热负荷问题,开展旋转爆震发动机燃烧室壁面气膜冷却的数值仿真,探究气膜出流与爆震波、斜激波和燃烧室流场之间的相互作用以及气膜对壁面的冷却特性。研究结果表明:爆震波对气膜的压缩和冷却孔的堵塞作用明显,气膜对爆震波整体的传播特性影响较小。受爆震波和燃烧室流场的影响,气膜出流存在周期性的摆动情况,这在一定程度上影响了壁面的冷却效果。在爆震波覆盖的壁面区域,峰值壁温下降程度有限,但时均壁温的降幅超过26.9%;在斜激波覆盖区域,随着冷气量的增加,峰值壁温和时均壁温的降幅超过32.5%和51.3%,气膜对该区域壁面的冷却效果更加明显。      

环形爆震室中火焰加速数值模拟及验证

对带不同数量孔板的环形爆震室进行数值模拟,并通过试验对数值计算进行验证,来研究火焰加速现象、爆燃向爆震转变过程和不同当量比下起爆距离.数值计算采用二维轴对称非定常Navier-Stokes方程来模拟流体动力学过程.研究发现用较低的点火能量对火焰混合区的可燃气点火产生低速火焰,低速火焰向环形爆震室射流并改变方向向出口传播,火焰在孔板的阻碍作用以及火焰诱导激波和反射波的加速作用下,由层流变为湍流,湍流火焰与其诱导激波相互加强,最终引爆未燃混气;还对爆震波在孔板区的传播过程进行了分析,对不同当量比下的火焰速度和起爆距离进行了模拟研究.      

二维旋转湍流边界层内气固两相流动的数值模拟

本文用修正的ＤＯＲＯＤ－ＦＥＭ方法对不可压流体二维旋转湍流边界层进行了数值模拟，得到了与实验数据较符合的边界层形状因子和壁面摩擦系数的计算结果。采用单向耦合（Ｏｎｅ－ｗａｙ Ｃｏｕｐｌｉｎｇ）方法，计算了固粒在湍流边界层和无粘流中的运动轨迹，给出了壁面被固粒撞击后的冲蚀磨损分布，分析了近壁流场性态与壁面冲蚀磨损的关系。     

旋转爆震发动机数值研究

采用9组分19步的基元反应模型,忽略黏性、热传导和扩散等输运效应,对以2H2+O2+4N2为反应混合物的旋转爆震发动机进行了二维模拟。采用分段填充方法起爆旋转爆震波,研究结果表明,分段填充方法可以有效的形成沿同一方向传播的爆震波;详细分析了内部流场结构,旋转爆震波会使流场呈现周期性变化,斜激波和接触间断面是造成流场在出口平面波动的主要原因。旋转爆震发动机具有较高的比冲,基于燃料的平均比冲为49131.5 N.s/kg。     

二维爆震波数值模拟

针对爆震发动机的多维效应,利用CFD计算软件对单管脉冲爆震发动机进行了二维的数值模拟。以氢气和氧气混合物为例,仿真了爆震波在管内的传播过程以及传出管外后的流场情况。模拟所得的爆震波参数与文献中的计算结果和实验值基本一致,验证了模拟的正确性。爆震波传出爆震管后退化成一道激波,对该激波的流动情况进行了分析。     

圆盘旋转爆震燃烧室中爆震波传播模态

为了研究圆盘结构下旋转爆震波传播特性,通过求解14组分19步CH4/O2反应的欧拉方程,数值模拟了不同燃料喷射温度、压力、燃烧室内外直径比值和燃烧室径向长度并分析了其对旋转爆震波传播模态的影响.研究表明,只有预混气喷射温度为500～900K,喷射压力为0.5～3.5MPa的条件下能够在固定燃烧室内形成连续的旋转爆震波....     

点火位置对圆盘结构下旋转爆震波起爆过程的影响

为分析点火位置改变对旋转爆震波(RDW)起爆过程的影响,在圆盘形旋转爆震发动机上进行点火实验,研究了不同点火位置、质量流率条件下旋转爆震波的建立过程及工作特性.结果表明,不同点火位置下,RDW的起爆过程皆经历点火器放电、爆燃转爆震以及稳定旋转爆震3个阶段.点火位置靠近喷注面附近时,RDW起爆过程中的无序缓燃模式缩短,起...     

旋转爆轰波与涡轮平面叶栅相互作用数值模拟

由于旋转爆轰燃烧室具有自增压特性,可提高热力循环效率,因此将旋转爆轰燃烧室应用于燃气轮机可进一步提高系统的性能。基于非稳态雷诺时均Navier-Stokes方法,采用剪切应力输运k-ω湍流模型,建立旋转爆轰燃烧室与涡轮平面叶栅耦合计算模型,研究旋转爆轰燃烧室内的复杂波系与涡轮叶片的相互作用,分析涡轮叶栅对高频爆轰压力振荡的抑制作用。结果表明：旋转爆轰燃烧室内的燃气在涡轮叶栅内加速,并且在斜激波后的局部区域马赫数的增加更为明显。斜激波与涡轮静转子叶片的前缘、压力面、吸力面以及尾缘相互作用,由于旋转爆轰波不同的传播方向,使得斜激波与静子叶片呈相互垂直或平行,进而形成两种不同的波系结构。涡轮叶栅对高频压力振荡存在明显的抑制作用,涡轮叶栅上下游高频压力振荡幅值的衰减率达到80%以上。研究结果展示了旋转爆轰波作用下涡轮叶栅内复杂波系结构特征,并对基于爆轰燃烧推进技术的应用提供了一定的理论基础。     

旋转爆轰波中多波流动模式的数值研究

为了研究旋转爆轰燃烧室内复杂的波系结构,特别是爆轰波多波传播模式的影响因素,采用两步诱导-放热总包反应模型,对简化的二维旋转爆轰波进行数值模拟。在不考虑环形爆轰燃烧室曲率的情况下,研究了进口总温和周向尺寸对单波及多波流场结构的影响。结果表明,一定范围内进口总温的增加会使得爆轰波的波头数目增加,且双波流场结构中存在双波对撞和双波同向两种传播形式,双波同向传播和三波同向传播之间则稳定存在着两组双波对撞流场,且有对应的温度范围;对于双波对撞模式,持续减小周向尺寸能够使得流场结构转变成单波模式;周向尺寸的增加则会使得双波对撞向双波同向传播模式转变。     

考虑燃料雾化的气液两相连续旋转爆轰数值模拟

为研究气液两相连续旋转爆轰发动机燃烧室内爆轰波的传播特性,以汽油为燃料,富氧空气为氧化剂,建立了欧拉-拉格朗日模型进行二维数值仿真,其中气相方程采用时空守恒元与求解元方法求解,液相方程采用标准四阶龙格库塔法求解。在两相旋转爆轰模型还考虑了液滴雾化破碎过程。计算结果表明：起爆后形成的初始爆轰波经过初始燃料填充区域后逐渐衰减,随后入口附近新生成的压力波经过一系列发展形成了自持稳定传播的旋转爆轰波;旋转爆轰波的传播模态受燃料与氧化剂的喷注压力和氧化剂填充比影响,在不同工况下旋转爆轰波呈现出4种传播模态,即稳定单波模态、稳定双波模态、不稳定双波模态和不稳定单波模态;在双波模态工况下,燃烧室内初始只形成1个爆轰波,后由入口附近局部爆炸产生的压力波发展为新的爆轰波,转化为双波模态后爆轰波的强度略有下降,但燃烧室整体推力更加平稳。     

连续旋转爆震波与涡轮静子叶栅相互作用数值研究

为了研究旋转爆震燃烧与涡轮部件组合的工作特性,对旋转爆震波与涡轮静子叶栅的相互作用过程进行了数值模拟,考虑了不同传播方向的影响,详细分析了爆震波与涡轮叶栅相互作用机理。结果表明,爆震波顺着叶盆方向传播时,在叶栅的叶盆处出现高温区,逆向传播时,同时在叶盆和叶背处都出现高温区;并且顺向传播时产生的反射波强度更大。分析了涡轮进出口压力和温度的变化过程,发现涡轮对压力的波动有一定抑制作用,顺向和逆向传播的爆震波经过涡轮叶栅后压力变化幅值分别下降了68%和57%。得到了对于当前叶栅构型,顺向传播的爆震波总压损失为11.03%,而逆向传播的爆震波总压损失为6.7%。     

扇形脉冲爆震燃烧室起爆过程数值模拟

为了探究不同障碍物结构对扇形脉冲爆震燃烧室起爆过程的影响,对采用3种典型孔板障碍物的扇形脉冲爆震燃烧 室,截取其中间截面进行了起爆过程的数值模拟,研究了障碍物型式、障碍物堵塞比及障碍物间距对扇形脉冲爆震燃烧室爆燃转 爆震（DDT）距离的影响规律。结果表明：障碍物堵塞比对扇形脉冲爆震燃烧室DDT距离的影响最大,且在堵塞比为0.35~0.70时, 堵塞比越高DDT距离越短;在3种型式的障碍物中,在低堵塞比下,截面为前掠三角形的障碍物缩短DDT距离的效果最好,在高堵 塞比下,截面为后掠三角形的障碍物缩短DDT距离的效果最好;在3种障碍物间距中,当间距分别等于1倍和1.2倍的爆震室当量 直径时,扇形脉冲爆震燃烧室的DDT距离相差不大,且均短于障碍物间距等于0.8倍爆震室当量直径时的DDT距离。     

激光支持爆轰波的仿真研究

采用了爆轰波的C-J起爆模型,计算了激光支持爆轰波初始参数,得到激光功率密度与爆轰波波面初始参数的关系.利用流体力学的方法建立爆轰渡在靶与约束层双重约束作用下的二维辐射膨胀模型,用Fluent软件进行仿真研究,得到了激光冲击强化后不同时间对应的爆轰波压力和轴向速度分布,仿真结果与激光冲击强化GH742试件后表面残余应力分布的实验结果一致.