不同矢量偏角下二元喷管内流特性的数值研究

应用多重网格MacCormack显式格式对不同矢量角下二元喷管的内部流动进行了模拟分析。计算时, 在细网格上采用了雷诺平均N-S方程与k-ε湍流模型方程;在粗网格上则采用Euler方程。数值计算给出了较清晰的流谱, 壁面压力分布的计算也与实验结果相吻合。通过计算结果分析了矢量喷管几何参数对流场的影响, 并对矢量喷管的推力特性进行了数值研究。      

壁面抽吸条件下的爆震波传播特性

为研究静止气中壁面抽吸对爆震波传播特性的影响,采用数值模拟方法研究了多孔抽吸壁面条件下爆震波的流场结构、传播速度等特性的变化规律。结果表明多孔抽吸壁面对爆震波传播特性有两方面的影响。首先,流场与多孔壁面发生碰撞产生弧形激波,弧形激波对爆震波的横波结构造成直接破坏,导致靠近多孔壁面区域的爆震波产生明显的速度亏损,爆震波强度削弱甚至熄爆;其次,多孔壁面作用导致了流场不稳定性增强,流场与孔板碰撞产生高温高压点,对爆震波的传播有促进作用,尤其在临界条件下,可能导致爆震波熄爆后重新起爆。在保持抽吸压力和边界条件不变情况下,对不同活性及不同抽吸距离时的爆震波传播特性进行研究,发现预混气活性降低、抽吸距离增长时,爆震波结构和传播速度受壁面抽吸影响增强。在加长抽吸距离条件下,随着预混气活性的降低,存在3种爆震波传播现象,即自持传播、熄爆后重新起爆以及完全熄爆。将相应工况下的实验研究与数值模拟结果进行对比,验证了数值模拟结果的正确性。     

膨胀波对斜劈稳定斜爆震波特性影响的数值研究

为探究燃烧室内可能存在的膨胀波对斜爆震稳定燃烧的影响，基于块结构自适应网格加密的AMROC程序，求解多组分可压缩化学反应流Euler方程，研究了膨胀波对斜劈诱导斜爆震流场的影响。发现无限长斜劈模型中，斜爆震波面角度是缓慢增加的，膨胀波影响下，斜爆震波面角度明显下降。从爆震波面角度变化中可以明显看出膨胀波的影响范围和爆震波的衰减程度。来流的低静压不会改变膨胀波影响范围，但容易导致爆震波衰减直至发生解耦，因此高马赫数低静压来流条件下斜爆震燃烧的稳定性和燃烧效率需要着重考虑。此外，使用普朗特-梅耶膨胀波基本原理对膨胀区进行分析，发现爆燃区流场参数与理论值吻合较好，且壁面附近前马赫线角度与近似膨胀波前沿较为接近；基于此，发展了一种定性评估斜劈末端膨胀波影响范围的手段，在膨胀波前马赫线角度基础上适当增加4°～10°，可以近似得到膨胀波的前沿位置。     

CFD技术在带动力飞机气动设计中的应用

应用数值模拟的方法,研究了飞机带动力发动机的流场模拟,基于此计算分析了其性能特性.采用基于雷诺时均的N-S方程,动力方程中的雷诺应力通过SST二方程湍流模型和自动壁面处理计算;同时,采用守恒隐式有限体积法.离散方程使用有限高精度对流程序类似于Barth和Jesperson进行计算.质量流动通过Rhie和Chow算法压力-速度耦合进行计算.     

环形爆震波聚焦起爆

本文对三种带不同结构形式的环形预爆管的爆震发生器进行了数值模拟,来研究环形爆震波聚焦起爆现象及其气动特性。数值计算采用多组分理想气体详细的化学反应机理、二维轴对称非定常流动的Navier-Stokes方程来模拟化学动力学和流体动力学过程。研究发现用低的点火能量对环形预爆管中的燃料和氧化剂点火产生层流火焰,层流火焰在狭窄管壁的作用下完成爆燃向爆震转捩,形成环形爆震波,环形爆震波在聚焦腔经过反射、汇聚,最终引爆主爆震室中的可燃混合物;并对稳定爆震的距离和强度进行了对比研究。通过分析对环形爆震波衍射、反射、聚焦的气体动力学特性及复杂波系发展规律有了更深的认识,为进一步试验提供参考。     

稀释氩气对分叉管内爆震波绕射的影响

为了研究稀释气体Ar的浓度对爆震波从直管往分叉管绕射传播和对直管中爆震波传播的影响,采用三阶TVD迎风格式和Strang-splitting算子分裂法,对H2/O2/Ar混合物爆震燃烧进行了二维数值计算。计算结果表明:在298K,6670Pa初始条件下,四种组分比条件下,主爆震管下游由于分叉管绕射稀疏波的影响,都会产生不同程度的解耦,但是通过壁面马赫反射后,都能恢复平面爆震;在H2/O2/Ar体积比为2:1:7时,爆震波绕射进入分叉管后,爆震波解耦成爆燃火焰波;体积比为2:1:0,2:1:1和2:1:4条件下,水平爆震管中的爆震波,绕射进入垂直分叉管后,通过不同次数壁面马赫反射,能够二次起爆,解耦的爆震波恢复成为平面自持爆震波。     

冲压发动机进气道压力振荡过程的数值研究

针对来流马赫数3.5,0°迎角飞行条件下冲压发动机的简化模型,数值模拟了进气道整流罩开启过程的非定常流动。计算采用AUSMPW有限体积格式,基于结构与非结构的混合网格技术,求解了非定常N-S方程。计算结果表明,在进气道整流罩打开以后,超声速气流会在燃烧室通道内形成振幅衰减较慢的低频剧烈振荡。在同样计算条件下采用Euler方程进行对比性研究,数值实验发现对这类大雷诺数、强非定常性的问题,Euler方程和N-S方程得到的结果相差不大。     

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

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

受侧向膨胀影响的爆震波传播过程研究

基于Fay流体扩张理论,运用实验观测、数值模拟和理论分析等方法研究了侧向膨胀对爆震波传播过程的影响.详细分析了爆震波在侧向膨胀段的流场结构,由于气体可压缩壁面限制性作用较弱,爆震波强度受到不同程度衰减,波面发生一定变形,横波数量明显减小且横波间距增大.对比研究了预混气活性对爆震波流场结构和传播过程影响,实验发现预混气活性较高时,即ER=1.0,爆震波能够抵御侧向膨胀衰减影响,能够自持传播;预混气活性较低时,即ER =0.7,其抵御侧向膨胀能力较弱,速度亏损变大,波面变形严重,爆震波解耦而熄爆.若要维持爆震波在气体可压缩壁面稳定传播,需要一定活性来抵御侧向膨胀影响.     

轴对称喷管内外流场与结构温度场耦合计算

利用有限体积的对流项二次迎风插值格式和重整化群(RNG)k-ε湍流模型,二层增强型壁面函数,同时利用球形谐波法考虑热辐射的影响,以灰气体加权模型(WSGGM)确定气体介质的辐射性质,求解N-S方程、热传导方程、考虑吸收-发射性气体介质的辐射传输方程。采用流固耦合的流动与换热模型,流场与结构温度场互为边界条件交换数据,实现了流场解算与温度场解算的耦合数值分析。采用此计算模型对静止的轴对称收-扩喷管进行了数值模拟,计算得到的喷管壁温与试验数据吻合良好。在此基础上,进一步开展了巡航状态下轴对称收-扩喷管内外流场与结构温度场的耦合数值研究。     

连续旋转爆轰波在无内柱圆筒内的数值模拟

在对同轴圆管内的连续旋转爆轰的已有研究成果进行总结分析后,提出无内柱的连续旋转爆轰的燃烧室模型.利用结合了一步化学反应模型的圆柱坐标系下的Euler方程三维数值模拟了在这种燃烧室内爆轰波的传播行为.数值结果表明:爆轰波由初始起爆的单向传播最终自动收敛稳定到若干个典型爆轰波头的流场结构.在当前进气条件下,流量可以达到450kg/m2s,其基于反应物的比冲约为1900s.与同以H₂/Air作为燃料时同轴连续旋转爆轰发动机(RDE)燃烧室模型的2000s相比,基于反应物的比冲十分接近,但由于内部有可燃气以爆燃形式燃烧等原因,有5%的损失.      

基于OpenFOAM的三维H2/Air连续旋转爆轰流场数值模拟

为进一步研究旋转爆轰流场特征,基于开源计算流体动力学软件OpenFOAM,采用9组分19步的基元化学反应模型,对H2/Air连续旋转爆轰流场进行了三维数值模拟,得到了旋转爆轰波稳定传播时燃烧室内部流场的详细结构,研究了燃烧室头部激波的传播特性,分析了旋转爆轰燃烧室的压力增益性能。结果表明:旋转爆轰波后的第一道反射激波在由燃烧室外壁面向内壁面传播过程中反射激波的高度增加并在靠近内壁面附近与滑移线交汇形成局部高温高压区域;旋转爆轰波在外壁面位置处相位约落后于内壁面0.003rad～0.15rad,其相位差随燃烧室曲率差的增加而增大;燃烧室头部反射激波数目受到曲率差和进气总压的影响,燃烧室曲率差增大,反射激波数目减少,进气总压增大,反射激波数目增多;燃烧室压力增益保持在0.3以上,在进气总压一定的条件下,压力增益随着燃烧室曲率差的增大有增加的趋势。研究结果揭示了三维旋转爆轰流场的精细结构和燃烧室头部激波的传播规律。     

燃烧室宽度对煤油旋转爆震波传播模态的影响—“爆震专刊”

为研究基于煤油的旋转爆震波的传播特性,以煤油和含氧量40%的富氧空气作为燃料和氧化剂,基于燃烧室外径均为100 mm的无内柱燃烧室和燃烧室宽度分别为32 mm、26 mm和20 mm的环形燃烧室开展了对比实验。不同氧化剂流量下,共观察到四种燃烧波模态,分别为爆燃模态、准稳定爆震模态、双波对撞模态和稳定旋转爆震模态。无内柱燃烧室中,氧化剂流量较低时无法维持旋转爆震波的稳定传播,出现爆燃模态和准稳定爆震模态;当氧化剂流量超过120 g/s时,可以得到旋转爆震模态,旋转爆震波峰值压力超过0.7 MPa,平均传播速度为1750 m/s。对于环形燃烧室,旋转爆震波的传播速度仅为1245~1465 m/s,明显低于无内柱燃烧室中的传播速度。随环形燃烧室宽度减小,对应旋转爆震波模态的工况范围更窄,传播速度更慢。在本研究对应的工况范围内,增大燃烧室宽度,更有利于基于煤油的旋转爆震波的稳定传播。     

旋转爆轰流场的数值模拟

基于带化学反应的二维Euler方程,采用氢气-空气的9组分19步基元反应简化模型,对充有当量比的氢气-空气预混气和空气的环形旋转爆轰流场,从点火燃烧到发展成旋转爆轰的过程进行了数值模拟。根据数值结果分析了波后流场中爆轰产物受激波、高温和离心力等作用而挤向外壁,形成有利于充入燃料,实现持续稳定旋转爆轰的流场特征。还讨论流场中爆轰波、激波与间断面和内外壁面反射或折射,从而形成多个激波相交的波系特征。为认识和理解旋转爆轰流场,开展旋转爆轰的实验研究等具有现实的指导意义。     

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

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

旋转爆轰的三维数值模拟

旋转爆轰利用爆轰波在燃烧室的供气端传播,而爆轰产物从另一端排出。为了揭示其宏观特性,采用带化学反应的三维Euler方程,在贴体坐标系中,利用带限制函数的波传播算法,对圆环形燃烧室内的旋转爆轰进行了数值研究,讨论了爆轰传播过程中的波系结构。数值结果表明,即使仅在燃烧室内壁面的部分区域充入可燃气,爆轰波也能够以旋转方式在该区域稳定传播。     

燃烧室宽度对液态燃料旋转爆轰发动机影响实验研究

为了研究燃烧室宽度对液态燃料旋转爆轰发动机工作特性的影响,搭建了气液两相旋转爆轰实验系统,以汽油/富氧空气为工质,氢气/氧气预爆轰管作为点火装置,在不同燃烧室宽度下开展了一系列实验研究,分析了爆轰波的起爆过程,以及燃烧室宽度对爆轰波传播特性与发动机推力性能的影响。实验结果表明：点火后,燃烧室内需要经过一个爆燃转爆轰过程才能形成自持传播的爆轰波;爆轰波在不同燃烧室宽度下均以双波对撞模态传播,对应的波速分布在850~1025m/s内,随着当量比增加,波速整体呈增加趋势;当燃烧室宽度减小,波速整体有所降低;不同燃烧室宽度下推力性能存在显著差异,其中燃烧室宽度在16.5mm下,发动机的推力和燃料比冲要明显低于11.5mm和9mm的;随着燃烧室宽度减小,内外壁面边界层在流场中的作用更为突出,降低了发动机推力的稳定性。     

连续旋转爆轰发动机参数特性的三维数值模拟

采用一步化学反应模型,基于任意坐标系下Euler控制方程,对连续旋转爆轰发动机(CRDE)进行了三维数值模拟研究, 详细分析了来流总压对CRDE参数特性的影响.研究发现:随着来流总压的增大,爆轰波峰值随之增大,但是燃烧室头部Laval型喷注段的壅塞比基本不变.不同来流总压下,可燃气体在燃烧室头部均以亚声速入射.随着来流总压的增大,燃烧室内流场的平均压强增大,但是其平均轴向流速基本不变.CRDE的流量、推力和比冲均随着来流总压的增加而变大;但是,不同来流总压下,CRDE的流量与Laval型喷注段最大流量之比不变,并且此比值约等Laval型喷注段的壅塞比.Laval型喷注段最大流量与来流总压成正比,因此以上分析从理论上进一步解释了CRDE的流量与来流总压成正比的原因.      

氙离子推力器阳极工质分配方式优化

使用两种数值方法研究微尺度下惰性气体的流动现象,验证表明壁面速度滑移条件的N-S（Navier-Stokes）方法能够较好地模拟离子推力器气体分配环内工质流动特性.为了降低氙工质注入放电室过程中气流速度,并改善其周向均匀性,针对现有气体分配环构型提出符合4级环切场离子推力器放电特性的改进方案:增加2级供气通道,进行双侧45°角开孔并增加气孔数目.数值计算表明:结构优化后氙气进入放电室周向均匀性提高37%,气流速度降低32%；推力器装配优化方案气体分配环进行实验,大束流工况下束离子电离能耗自183W/A下降至167W/A,离子推力器放电室性能得到了提高.