各向异性网格自适应计算在超燃模拟中的应用

采用各向异性网格自适应求解技术,将其应用于DLR超燃冲压发动机燃烧室中的超燃模拟.开展了3个算例,包括采用滑移壁面条件的燃烧室冷流场模拟、采用无滑移壁面条件的冷流场模拟及采用无滑移壁面条件的反应流模拟.模拟中,各向异性网格自适应计算捕捉到了如激波、射流、边界层、火焰面等具有各向异性特征的大梯度区域,并利用各向异性网格进行了很好的加密.对比利用各向同性网格的初始流场计算,各向异性网格自适应计算使基于滑移条件、无滑移条件的冷流计算及反应流计算的网格单元数量分别下降了36.2%,36.4%和36.8%,有效降低了计算规模,而且流场大梯度区域的计算结果更准确,辨析度更好.结果表明:对于像超燃这类具有各向异性特征的问题,各向异性网格系统比各向同性网格系统有更好的计算效率及准确性,同时也表明基于Mach数场构造的各向异性网格系统可以有效应用于超燃计算.      

一种三维并行自适应计算方法

基于串行网格划分软件METIS与并行化消息传递编程接口（ MPICH2）对现有串行自适应程序进行简单的并行化改造,给出了一种三维可压缩无粘流数值模拟的并行自适应方法。首先利用单个进程调用METIS,串行划分网格;然后对所有进程并行计算以获得初始网格下的流场解;再次利用单个进程对整个流场运用自适应方法进行局部网格加密并调用METIS串行划分网格;最后全部进程在流场初始解的基础上继续并行计算,以获得自适应网格下的流场解。数值模拟算例验证了此方法的可靠性与高效性。     

自适应湍流耦合建表燃烧模型的振荡燃烧数值模拟

航空发动机及燃气轮机等动力装备燃烧室广泛采用的贫燃燃烧方式经常遇到破坏性的非定常热声耦合振荡燃烧问题。非定常振荡燃烧数值预测是一个长期的研究热点和难题。发展了针对振荡燃烧的耦合直接求解数值模拟方法,包括优化的动态模型参数的高精度自适应湍流模型（SATES）,耦合可压缩的详细化学反应建表FGM燃烧模型。选取的3种湍流燃烧模型包括有限速率模型（W1）及火焰面密度封闭方法中的Zimont(W2)和Fureby(W3)2种褶皱因子模型。针对经典的LIMOUSINE燃烧室多个部分预混振荡燃烧工况开展了数值研究,发现自适应湍流模拟框架下的3种燃烧模型均准确预测到了振荡燃烧的振荡频率,与试验相比,误差<6.4%;对于振荡燃烧压力脉动振幅的预测结果,有限速率模型（W1）和Zimont(W2)模型结果显著大于试验值,误差>380%;Fureby(W3)模型结果与试验值吻合较好,误差<17.9%。表明振荡燃烧的数值预测对不同的湍流及湍流燃烧模型具有较高的敏感性。不同的工况结果表明,振荡燃烧存在完全振荡模态和过渡模态,完全振荡模态中数值预测的特征主频在燃烧室上下游多个位置趋于一致;过渡振荡...     

不同PDF模型在超燃冲压发动机数值模拟中的应用

湍流与燃烧相互作用对复杂超声速燃烧流场产生影响,基于设定型PDF(ProbabilityDensityFunction)模型对其不同分布模型进行研究。考虑温度和组分联合概率密度函数,采用β-PDF和δ-PDF方法对HF2直连式发动机模型进行模拟,并与该发动机风洞试验数据进行对比。结果表明:无论是定量还是定性分析,两种PDF方法的预测结果与实验结果吻合好,说明湍流与化学反应相互作用不能忽略;采用β分布的算例壁面压强分布结果与实验结果吻合更好,燃烧效率更高,更能够捕捉湍流对于化学反应的影响。     

采用不同建表方法的火焰面模型在燃烧室中的应用研究

为加深对航空发动机燃烧室中湍流燃烧过程的理解,采用不同建表方法的火焰面模型对航空发动机模型燃烧室内的湍流燃烧过程进行数值模拟,包括层流火焰面数据库的构建和反应进度变量的PDF类型两个方面。其中,层流火焰面数据库的构造方法包括基于扩散火焰的FPV和基于预混火焰的FGM模型,反应进度变量的PDF类型包括δ和β分布。LISA和KHRT模型分别用于模拟液膜和液滴的破碎过程,非平衡Langmuir-Knudsen模型用于模拟液滴的蒸发过程。LISA模型得到的液膜破碎距离约为4.6mm,液滴直径在文氏管出口下游迅速减小到10μm左右,并在头部出口下游附近完全蒸发。通过与相干反斯托克斯喇曼散射(CARS)和可调谐二极管激光吸收光谱(TDLAS)测量温度的对比,验证了FPV和FGM模型的精度,并表明在流动变化较大的位置FPV模型具有更高的精度,而其他位置FGM模型具有更高的精度,采用β分布作为反应进度变量PDF的模型,可以有效提高温度的预测进度,而且主燃区内的误差基本都在5%以内。此外采用β分布作为反应进度变量PDF的FGM模型,可以更好地描述未燃混合物被回流燃气点火的过程,而且反应进度变量的PDF类型比层流火焰面数据库构建方法的影响更为显著。     

自适应多重网格法在跨声速势流计算中的应用

用全近似格式求解了轴向大扰动，径向小扰动的跨声速势流方程。并采用自适应多重网格技术控制粗细网格层间的迭代转换，采用嵌套迭代技术较好地选取了初值。数值计算表明：自适应多重网格法可使计算效率提高８倍以上。     

网格自适应策略在高升力构型计算中的应用

提出了一种基于多尺度网格误差分析的新型混合网格自适应局部加密策略。在初始网格基础上,通过多重网格计算,比较次细网格和细网格的计算结果,揭示出网格敏感区域,然后对这些网格敏感区域加密,计算自适应加密后的网格得到最终计算结果。相比于常规网格自适应策略,提出的新策略特别适合多种流动现象并存的高升力构型绕流流场计算中的敏感区域判断。NASA高升力全展构型的数值试验充分表明这种新策略的有效性。     

自适应网格局部加密方法在可压缩流场模拟中的应用

本文基于M.J.Berger和J.Oliger发展的自适应网格局部加密算法,它采用多个组合网格的思想,用于求解双曲型方程组。它结合显式有限差分格式,采用Richardson外推技术自动进行局部截断误差估计,对精度低的区域产生新细网格局部加密,或去除无需再加密的旧细网格,以最小的运算量达到给定的精度要求。网格可一层层加密下去,按层次覆盖,每个网格是有任意方向的矩形均匀网格。这套算法独立于求解所用的差分格式,很容易和各种格式结合。 我们从M.J.Berger和J.Oliger程序出发,参考M.J.Berger和P.Colella,对守恒型方程组,实现了粗细网格交接面上数值通量守恒,从而能计算间断解。对网格边界处理、误差估计和网格重新生成等做了一些改进,并且可以灵活地在任何感兴趣区域指定加密。 我们用不定常Euler方程组计算了带前台阶的二维管道中Mach-3流动,采用MacCormack显式附加人工粘性差分格式。共用了四层网格结构:Δ=1/10、1/20、1/80和1/320。其中基本网格是Δ=1/10;到台阶上第一次激波反射区域从Δ=1/20加密至Δ=1/80,计算到定常所需运算量仅是同等均匀网格△=1/80的28％;我们还用Δ=1/320网格对台阶拐角和切向间断处指定加密,结果表明可明显改善台阶上的激波Mach反射,切向间断清晰可辨,所需的计算量是上述均匀网格Δ=1/     

新的参考模型自适应控制方法在航空发动机中的应用

本文针对航空发动机模型的结构和参数不能明确已知,提出一种阶次低、鲁棒性较好的参考模型自适应控制法,即只需知道系统的相对阶次,且控制器结构简单,易于实现,适应性强。      

平面激光诱导荧光技术在超声速燃烧中的应用

使用283.553nm的紫外激光激发对温度不敏感的氢氧基Q1(8)线,拍摄代表氢氧基浓度分布的受激荧光发射图像,以此来分析超声速燃烧火焰的结构。用这种方法研究了不同凹腔长深比、凹腔后缘倾角和不同燃料喷注方案对超声速燃烧火焰结构的影响。观察到了凹腔与激波对点火和火焰稳定的作用,凹腔长深比和后缘倾角对燃料穿透、混合和燃烧的不同扰动作用,以及火焰中存在明显的湍流结构和分区燃烧现象,并进行了简要分析。结果显示PLIF是超声速燃烧研究的有力工具。     

双燃式冲压发动机中富油燃气射流的超燃研究

在空气流量1.2kg/s左右的地面连管试验台上，研究了双册槽结构的超燃室中多股高温富油燃气射流的超燃特性。试验结果表明，超燃点火可靠，火焰稳定，超燃效率可达0.8以上。因此超燃新方案是可行的。     

基于气动斜坡/燃气发生器的超燃燃烧室实验

为探索超燃冲压发动机燃烧室中的新的火焰稳定技术,提出了一种新型被动式燃料掺混增强技术—气动斜坡与燃气发生器组合燃料喷注技术,并在北航直联式超燃试验台对这种新型组合喷注器开展了超声速燃烧的试验研究。在模拟飞行马赫数5(燃烧室入口Ma=2),进行了冷流试验,获得了喷注器附近流场的纹影图像。本文设计了4种气动斜坡喷注单元,以乙烯为燃料,在约1kg/s试验气流中开展了多级喷注单元组合的超声速燃烧试验,在当量比0.78～1.22范围内实现了稳定的燃烧,经冲量分析法计算得到不同组合结构的燃烧效率为0.54～0.72。试验结果验证了这种新方案作为超然冲压发动机火焰稳定装置的可行性。     

超声速燃烧凹槽火焰稳定的研究动态

高超声速气流在燃烧室的停留时间非常短，使得超燃冲压发动机燃烧室内燃料与气流的混合及其燃烧变得非常困难。目前集燃料喷射、混合及火焰稳定为一体的凹槽有望改进这一情况而受到了普遍关注。介绍了超声速气流流过凹槽的自激振荡及其控制、停留时间和阻力等特性，总结、分析了近几年凹槽在超声速气流中增强混合及燃烧的实验与数值研究，指出了在把凹槽用于提高超燃冲压发动机性能方面存在的问题和进一步研究的方向。     

超声速进气道喉部附面层抽吸

为研究超声速进气道喉部之后流场激波附面层干扰,采用FLUENT软件模拟了单楔角进气道在设计工况下流动情况。通过分析,提出进气道喉部抽吸。计算了三种抽吸缝大小下进气道喉部之后流场,计算结果表明,喉部抽吸能使激波稳定于喉部,通过抽吸能改善喉部之后流场状况,提高进气道性能,少量抽气不改变流场结构,加大抽气量,使喉部之后激波串转变成正激波,正激波之后流场不分离,进气道出口性能参数提高显著。     

氢/碳氢燃料超声速燃烧的数值模拟

采用隐式上、下三角分解（Lower-Upper Decomposition）的算法，结合对组分连续方程中化学反应源项的点隐式处理，求解了多组分的N－S方程组，得到了二维和三维条件下的超声速掺混及燃烧的数值模拟结果，通过与水平喷射氢气及预燃煤油等算例的验证，计算结果与实验数据基本符合。在此基础上进行了垂直喷射氢气的三维超声速燃烧的数值模拟，计算结果显示了氢气在超声速空气中掺混、燃烧的过程，该程序可进一步用来模拟超燃冲压发动机燃烧室内的复杂流场。     

凹槽火焰稳定器结构对煤油超声速燃烧的影响

采用混合分数平衡化学反应模型和离散液滴模型，研究了不同长深比和后缘倾角的凹槽在液体煤油双模态燃烧中的火焰稳定特性。结果表明，集燃料喷射、混合及火焰稳定为一体的凹槽，在所模拟的工作过程中增强了燃料与空气的混合和燃烧；并得到最优的凹槽结构；混合分数平衡化学反应模型和离散液滴模型能准确地预测煤油在双模态燃烧室内的喷雾燃烧。     

凹腔结构对超声速燃烧室中横向燃料喷流流动与燃烧的影响

在超声速燃烧室的壁面上安装具有不同长深比、深度以及后壁倾角的凹腔,在凹腔上游壁面横向喷注燃料射流。分别利用掺有丙酮蒸气的氦气和氢气模拟不燃烧和燃烧两种情况下的横向燃料喷流,利用平面激光诱导荧光技术分别对流场中不同截面上的丙酮和氢氧基进行成像,同时对不燃烧情况下的喷流流场进行了数值仿真。研究发现:凹腔结构主要通过影响由凹腔后壁高压区向凹腔前壁传播高压扰动的行为来影响凹腔内部的静压分布,从而影响燃料的流动过程;凹腔长深比减小、深度或后壁倾角增大都有助于高压扰动的前传,导致燃料不易向凹腔内部偏转以及进入凹腔的燃料质量减少;凹腔长深比和深度的增加可增大凹腔内低速回流区的范围,有助于增强凹腔内部的燃烧以及火焰稳定能力;凹腔后壁倾角对燃料燃烧的影响不显著。