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

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

高雷诺数槽道湍流的壁面模化大涡模拟研究

选取基准壁湍流的槽道流动,研究了多种模型的壁面模化大涡模拟.模型包括经典的大涡模拟、Spalart-Allmaras、分离涡模拟和一种动态混合模型.基于摩擦速度的雷诺数范围为395~12000,采用3组粗糙网格,流向和展向维数分别同取37,49和65,法向维数保证y+(1)~1.主要研究平均速度、雷诺切应力分布、详细分析了各模型的特性差异并展示了相应的湍流结构.研究表明:在高雷诺数粗糙网格下,大涡模拟失去求解精度,分离涡模拟出现对数律不匹配,动态混合模型的计算接近直接数值模拟,其对数率区可解应力约占雷诺切应力的93%,边界层外层可解应力约占99%.这说明合适的混合模型可以在经济成本下保证计算精度,具有解决实际问题的潜力.      

不同钝体宽度下的驻涡燃烧室排放试验

为了研究空气流量分配对驻涡燃烧室对排放特性的影响,了解对驻涡燃烧室内污染物生成的过程及其影响因素,设计了一个能够改变中心钝体宽度、仅凹腔供油的驻涡燃烧室.在常压下对该驻涡燃烧室进行了排放特性试验,进口温度保持200℃.试验中,燃烧室进口马赫数为0.15~0.3.影响排放的因素主要包括雾化质量、凹腔当量比以及与进口马赫数相关的驻留时间等.总体来说雾化质量、凹腔当量比的提高对降低CO和HC的排放是有利的,但是这会使NOx排放增加.在低凹腔当量比时,CO排放曲线变化下降比较平缓,甚至出现上升趋势,而HC排放曲线比较陡峭.这是由于HC的消耗速度比CO消耗速度快,随着凹腔当量比的增加,供油压力提高,燃油雾化粒径变小,燃油蒸发时间缩短,使HC排放快速减少,中间产物CO大量产生而来不及消耗.凹腔当量比进一步上升时,由于燃烧温度的提高,使得CO排放快速减少.在燃烧室内燃烧过程中,NOx的形成和消耗是非常复杂的过程,目前只能作一些定性的分析,而CO和HC的反应过程相对简单.通过对不同钝体槽宽下,具有相似凹腔前壁流量的工况的比较,发现CO和HC的形成主要受凹腔内工作状况影响,而NOx的形成过程更复杂,主流也对其产生着重要的影响.      

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

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

含激波的湍流流动高精度大涡数值模拟方法

针对采用亚格子模型进行含激波的湍流流动模拟时会面临激波附近的精度损失问题,考虑从通过亚格子模型以及数值模拟方法两方面的改进来实现湍流流动大涡模拟的精度提高.大涡模拟采用了Yee及Sj(o)green (2009)提出的高阶低耗散方法.该方法采用自适应的流场探测器以控制计算中所需区域的数值耗散,并考虑对动力学模型采用在激波位置使用Sagaut和Germano(2005)提出的单边亚格子过滤器和(或)直接禁用亚格子项等方法加以改进.对于标准的马赫数1.5和3条件下的激波-湍流干扰问题,上述新方法相较于全区域采用亚格子模型的方法均表现出了相似的精度提升.同时实现的数值精度改进方案采用了Harten的亚单元分辨过程来定位和锐化激波,并在精确激波位置附近的网格点处采用了单边测试滤波.     

火花放电合成射流与超声速来流相互干扰特性数值模拟研究

为了指导火花放电式合成射流激励器在超声速流动控制中的应用,数值模拟研究了火花放电合成射流与超声速来流的相互干扰特性。研究表明火花放电式合成射流在超声速流场中产生强烈扰动,产生较强的激波结构;随着射流的喷出,激励器上游分离区和流场中激波呈先增强后减弱的趋势,激波由弓形激波逐渐弱化为斜激波,并且随着放电能量的增加射流与主流的动量通量比不断增大,射流的干扰和控制能力显著增强。由于超声速流的较大惯性及其对腔内气体的引射作用,激励器的腔体回填速率大幅下降、回填时间明显增长,使得激励器的工作频率受到很大限制。     

基于特征线追踪的气动反设计

随着当今飞行器飞行速度的提高,三维超声速内流道的设计需求日益增加。提出一种基于特征线追踪的气动反设计方法,以解决三维超声速粘性流场设计问题。为了验证特征线追踪方法的设计能力,将其应用于二维超声速喷管和三维方转圆流道设计中。结果表明,设计所得流道内部不存在集中的膨胀波和压缩波,获得了流动均匀的菱形试验区,出口处马赫数误差低于0.5%。特征线追踪方法能够直接设计超声速粘性流场,避免了传统的边界层修正技术引入的设计误差,实现了既能满足预设流场参数分布,又能保证入口和出口形状的三维异形截面超声速流道设计。     

四代机发展之辨

最近一段时间以来,四代机的话题被媒体炒得很热。去年,美国抛弃了原来的四代论,取而代之提出了一个五代论的观点,在2007年第6期的《航空知识》上,我和傅前哨、徐德康两位同志也一起讨论了这个话题。此外,俄罗斯、日本、韩国     

人类何时才能突破高超声速屏障

高超声速飞行在广义上可泛指在大气层中马赫数大于5的飞行。在这个意义上,飞船、航天飞机(包括美国的X-37B)、导弹的再入段,都有一段是高超声速飞行,包括它们的上升段,也有一段是高超声速飞行,但它们都没有进行高超声速巡航,遇到的问题相对简单。在本文中我们着重讨论用吸气式发动机在大气层中进行高超声速巡航的飞行器,以及用火箭助推到高空然后滑翔     

管流转捩对比态定律的讨论：层,紊流并存流态转捩现象

转捩对管内流动和换热有很大影响,Ｎ．Ｔ．Ｏｂｏｔ曾在总结归纳圆管、矩形管、三角形管和环形管实验结果的基础上提出对比态定律,并断言：对任意截面管道,只要知道临界摩擦因子和临界雷诺数,根据对比态定律,就可以计算其流动和平均换热特性。本文通过对复杂截面通道流动的研究,指出管流转捩的对比态定律只适用于单一流态管流,而不适用于层、紊流并存管流,文中还分析了层、紊流并存现象及其产生的原因,它的起始点和终止点的判定,以及流动和换热特性计算的处理方法。