可压涡卷对并的燃烧反应放热效应分析

通过以高精度迎风/对称紧致差分杂交算法对组元输运守恒方程和含源项流动Navier-Stokes方程的全耦合求解,在无放热、低放热和高放热三种情形下,对可压平面混合层中大尺度涡卷对并所受反应放热效应现象进行了直接数值模拟研究.对比显示,高放热下无涡并出现,亚谐扰动将基频涡卷雏形撕裂,直接形成大尺度亚谐涡卷.文中还分析了不同反应放热量对混合层厚度、产物生成率以及流场脉动行为特征等方面的影响.     

超声速剪切层的混合问题

用数值模拟和物理分析相结合的方法,研究了入口速度剖面具有阶跃型的超声速混合层的混合机理和特征.指出在上游起始阶段,混合层产生的展向涡是平行发展的,近似二维情况,涡的运动先是稳定的,后经非线性不稳定性和分叉,并进一步发展为两涡合并或两合并涡的再合并.随着向下游的增加,由于展向出现物理量的三维效应,展向涡要弯曲,并且沿其轴向要分叉演化,产生一个或多个极限环,这就开始改变二维涡合并的发展规律,产生流向和法向的旋涡.如果进一步走向下游,展向涡有的要产生螺旋型和泡型破裂,并且破裂涡与上游来的涡要混掺、缠绕,形成中间夹有小涡的新的大涡相干结构,相应流向、法向涡进一步非线性增强,流动完全改变了二维混合规律,变成真正的三维混合.由于涡的破裂带有随机和间歇性,这种混合过程也具带有随机和间歇性.通过本文对概率密度分布、分数维和间歇因子的计算,证明在流动出口处,基本上已达到了湍流.文中还给出了转捩发生的特征和机理.     

二维超声速剪切层的非线性失稳过程分析

利用直接数值模拟方法，模拟了二维超声速剪切流动中的小扰动波随时间发展在空间上的传播过程，由此分析其稳定性，通过相图分析，李雅普诺夫指数分析和分维分析，本文得到了二维超声速自由剪切层流动由线性到非线性失稳过程的数学特征，结果分析表明，在本文的计算条件下，二维超声速混合层流动在失稳过程中出现了间歇现象，并进入了混沌运动的初期阶段。     

大涡模拟超声速两相流混合层的非线性分析

为研究超声速两相流混合机理,采用大涡模拟(LES)方法数值仿真了超声速混合层内液滴两相流流场结构.用二维盒式滤波法处理可压缩气相N-S方程,得到亚格子形式的LES控制方程;用轨道模型仿真液滴运动,并用单液滴蒸发模型模拟液滴蒸发过程.利用伪邻近法和互信息量法确定了相空间重构的嵌入维数和时间延迟,对混合层压力时间序列进行了相空间重构.分析了仿真结果中混合层的涡结构特性和液滴对流场的反作用效果,并讨论了混合层失稳过程的非线性特征.发现混合层发展的不同阶段与其系统的非线性特征值密切相关.      

三维效应对超声速湍流燃烧流动大涡模拟的影响研究

本文对空间发展Mac=0.61超声速湍流燃烧混合层流动开展了二维/三维大涡模拟（LES）研究,探讨了三维效应在超声速湍流燃烧流动的大涡模拟中的重要性。研究发现,相较于三维LES计算结果,二维LES大幅低估了超声速湍流燃烧混合层的增长率,湍动能的二维LES结果显著偏大而雷诺剪切应力显著偏小。因而,二维LES难以准确模拟超声速湍流燃烧混合层的增长率和湍流脉动特性。二维流动中漩涡拉伸、压缩机制失效,致使二维LES的漩涡结构与三维LES结果显著不同：二维LES计算所得流场呈现清晰的旋涡生成、发展、并对和破碎等演化现象,而三维LES的湍流结构呈不规则状。湍流与燃烧间存在强烈的耦合相互作用,由于二维LES难以准确模拟湍流特性,燃烧场进而在二维LES下呈现一系列特殊的、不同于真实燃烧流动情况的数值现象。漩涡运动带动燃料和氧化剂掺混,掺混界面为火焰前锋。由于漩涡结构不同,致使二维LES的火焰前锋面面积不同于三维LES结果,导致其难以获得准确的燃料效率,最终导致二维LES计算获得的统计平均火焰温度显著偏低。     

超声速混合层液滴两相流动的大涡模拟

为研究超声速气流中液滴与气流的混合及液滴蒸发对混合的影响,采用大涡模拟(LES)方法数值仿真了超声速混合层内液滴两相流流场结构,气相流场采用亚格子(SGS)模型和切应力输运(k-ωSST)湍流模型,液相模拟采用轨道模型和单液滴蒸发模型。在混合层前缘入口处均匀持续地投放液滴,并在液滴入口处下方添加非周期小扰动,并观察液滴蒸发过程对该小扰动造成的影响。分析了入口小扰动在流场中不同的发展情况,发现液滴的蒸发过程使混合层厚度增加并加速混合层的发展,对气相流场扰动较强,可能导致流动失稳,对混合过程有很大的促进作用。     

超声速平面剪切层多模态混合扰动数值分析

对Ma_c=1.2二维空间发展超声速平面自由剪切层进行扰动模态及流动结构数值分析.使用线性稳定理论和扰动直接数值模拟方法对超声速剪切层中不同流向站位的流动剖面进行扰动参数分析和局部扰动直接数值模拟,获得了剪切层不同站位处的扰动特性.通过特定单频扰动下的直接数值模拟,在剪切层的不同区域内,无相差地促发了慢模态、混合模态和快模态等多种声涡模态,得到了特定单频扰动促发的多模态无相差不稳定波.      

二维超声速混合层流动稳定性的数值分析与并行计算

本文用DNS方法对两股空气的二维超声速混合流动(Mc=0.5)的稳定性问题进行了计算模拟.通过对采集的流动参数在时空两个方面进行的分析,发现流动失稳后出现了间歇结构.有限差分计算采用了并行计算技术,在微机并行环境下可高效率地完成大规模计算任务.     

三维超声速不等密度平面混合层的小激波形成与演变

基于大涡模拟方法与高精度混合格式,对两种密度比(1.0与0.138)条件下的高对流马赫数(Mac=1.1)三维空间发展超声速平面混合层进行了数值模拟。数值结果清晰描述了扰动激励下超声速平面混合层的演变过程,讨论了三维混合层中小激波的结构及其形成过程。分析了混合层两侧密度变化对混合层失稳及小激波形成的影响,发现当上、下层密度不等时,混合层失稳加速,小激波强度减弱,而密度较低侧的小激波会更弱。     

高对流Mach数三维混合层转捩特性分析及小激波结构模拟

采用空间大涡模拟方法对超/超混合、超/亚混合两类三维可压缩平面混合层转捩及其全湍流流场进行了研究,认为混合层动量厚度饱和点可作为流场转捩完成的标志.计算所得到的线性扰动波激励下流场转捩拟序结构与随机扰动下自然失稳结构以及文献结果进行了比对,其结果是一致的,表明了引入线性扰动激励来研究流场转捩结构是合理的.同时,本文还在较高对流Mach数流动下得到了三维流场动态小激波结构,其分布具有非对称特性,且形状与实验及直接数值模拟结果相似.不同条件混合层转捩计算表明:高对流Mach数下混合层转捩以Λ涡结构的形成和发展为主导机制,受扰动及对流条件的影响Λ涡结构不尽相同,某些情况下流场出现二维与三维涡结构共存现象.充分发展湍流区域,流场脉动速度分量量级相同,湍流压缩效应随着对流Mach数提高而明显增强.     

涡流器燃烧室头部两相反应流数值模拟

在任意曲线坐标系中数值模拟带涡流器燃烧室头部的三维两相反应流场,其中,采用偏微分方程法生成块结构网格;气相在Euler坐标系下采用修正的κ-ε双方程模型处理紊流特性,用修正的EDM-Arrhenius紊流燃烧模型模拟化学反应速率,用六通量模型考虑热辐射的影响;液相处理采用在Lagrange坐标系下颗粒群轨道模型.在非交错网格系下采用SIMPLE算法及混分差分格式求解离散方程,采用单元内颗粒源法(PSIC)算法进行两相之间的耦合,计算所得的旋流数与试验数据相当接近.结果表明涡流器几何形状对燃烧室流场影响较大,涡流器的优化设计可以提高燃烧性能.     

不可压湍流反应流直接模拟和RANS反应模型检验

用谱方法对三维不可压槽道湍流反应流动进行了直接数值模拟,得到了温度和质量分数耦合的瞬态数据库.结果显示温度与质量分数脉动在近壁区都有条带结构.推导了雷诺平均方法中湍流二阶矩反应模型中所求关联量的精确输运方程,发现耗散项是化学反应率系数-质量分数关联量封闭的关键,耗散项需要考虑化学反应影响.基于数据库的统计结果,对关联量模型方程中各项进行了先验研究,发现产生项和耗散项的贡献最大,扩散项和反应项的贡献较小,化学反应对各项大小和分布形状有明显影响.在算例中,直接模拟统计得到的化学反应率系数-质量分数关联和用代数二阶矩模型的模拟值很接近,说明ASOM模型具有一定的合理性.      

超燃燃烧室流场计算方法比较分析

采用准一维流动模型、轴对称方程、三维数值模拟方法结合13组分、32反应氢燃料反应机理计算了等直段+扩张段燃烧室内超声速燃烧问题.结果表明准一维流动模型与三维数值模拟方法的计算结果反应了主要气动参数变化规律;除预燃激波串区域外,计算数据与试验数据吻合良好;准一维流动模型计算快速、适应性强;三维数值模拟方法给出了详细流场细节,所需计算时间长.针对两种模型特点,在超燃燃烧室设计过程中应有机结合协同工作.      

二维凹槽超声速湍流流动数值模拟

本文采用时间、空间均为二阶精度的NND差分格式,以及Menterk ωSST二方程湍流模型,数值模拟了二维开式凹槽超声速粘性流动。用组分扩散方程计算了凹槽内粒子驻留时间。本文给出了M∞=3,长深比L/D=3,后壁倾斜角为90°和30°两种几何形状凹槽的计算结果,与相应的实验符合甚好。计算亦表明,该模型较高地估计了涡粘性,从而使流动更稳定。     

回流环形燃烧室反应流场数值模拟

采用自编计算程序,在任意曲线坐标系下对带扩压器、涡流杯、火焰筒、内外环冷却通道与弯管的回流环形燃烧室三维两相化学反应整体流场进行计算.数值预测火焰筒头部结构参数对燃烧流场的影响.所得计算结果与实验数据符合较好,结果表明:采用的数学模型和计算方法合理,该软件可供回流燃烧室优化设计与研制使用.      

BL模式在超声速湍流边界层中应用的检验

湍流模型是目前解决工程中湍流问题不可缺少的工具.然而它终究是建立在简化与假设的基础上的,其准确性必须经过实际的检验.对超声速边界层,由于缺乏系统的实验数据,检验工作很难做.本文以BL模型为例,通过可压缩湍流边界层的三个典型算例,与直接数值模拟(DNS)的结果作对比.分析了其成功及不足之处,以及影响其精度的若干因素,并检验了用BL模型预测转捩位置的有效性.     

回流燃烧室燃烧过程的三维数值模拟

在三维任意曲线坐标系下数值模拟回流燃烧室火焰筒内两相燃烧过程,采用RNG k-ε模型模拟紊流粘性,EBU-Arrhenius模型模拟燃烧反应速率、离散坐标模型以及六通量模型考虑辐射传热,液相采用颗粒轨道模型,气相采用SIMPLE算法求解,并用PSIC算法考虑气液两相之间的相互作用的影响,计算得到燃烧室内速度、温度等各气流参数分布.通过将计算与实验结果对比表明,计算方法可靠,离散坐标模型优于六通量辐射模型,更适用于模拟火焰筒内两相燃烧流场.      

超声速高度欠膨胀冲击射流的大涡模拟

冲击射流广泛应用于短距起飞垂直降落飞行器(SVTOL)等航空航天领域.本文采用大涡模拟方法对高度欠膨胀的超声速冲击射流的流场进行了数值模拟.本文数值模拟得到了高度欠膨胀冲击射流流场中的激波结构和内外剪切层中不同尺度的涡结构.数值结果观察到了马赫盘的振荡,以及在斜激波、马赫盘及大尺度涡结构的共同作用下,射流内外剪切层之间的环形激波的生成与消失的周期过程.并对流场内剪切层的涡结构的演化进行了研究,数值结果显示内剪切层的大尺度涡结构的形成与马赫盘的振荡相关,在内外剪切层的作用下形成了壁射流区内外交错的涡结构.     

固液火箭冲压发动机内两相反应流场数值计算

为了计算固液混合式火箭冲压发动机补燃室内的三维反应流场，用块隐式法求解气相Navier-Stokes方程组，用连续介质模型和k-ε-Ap模型计算颗粒相的湍流流动与蒸发过程，用修正的k-ε-g模型描述燃料的燃烧。为了分析发动机设计参数对反应流场的影响，用不同的条件进行计算，并由此分析了补燃室几何结构和液体燃料切始颗粒直径对燃烧效率的影响。算例表明，计算方法有效可行。数值结果能够反映流场结构、液体燃料的蒸发和两种燃料的燃烧过程。