高温圆管内多孔介质耦合换热的辐射作用

通过数值模拟圆管内填充多孔介质的辐射对流耦合换热,研究高温固体骨架辐射效应对温度分布及换热的影响.基于局部非热平衡模型分别建立流体和固相能量方程,采用蒙特卡罗法求解固体骨架的辐射换热;对不考虑热辐射引起的温度场偏差和管壁发射率以及多孔结构参数的影响进行讨论.结果表明:多孔固体骨架的辐射效应对入口段温度场的影响明显,不考虑辐射将导致较大偏差,壁温为1500K时最大偏差为16%.管壁发射率对温度场影响较小,壁温为1500K时影响小于3%.孔隙率或孔径增大,壁面辐射热流密度比例增加,辐射效应体现明显.      

三维等离子体MHD气动热环境数值模拟

电磁流动控制技术是一个多学科交叉融合的重要研究方向,在高超声速飞行器气动特性优化、气动热环境减缓、边界层转捩和等离子体分布等流动控制方面显示出广阔的应用前景。考虑高超声速飞行器绕流流场中发生的离解、复合、电离和置换等化学反应,气体分子振动能激发以及化学非平衡效应,耦合电磁场作用并基于低磁雷诺数假设,通过数值模拟求解三维非平衡Navier-Stokes流场控制方程和Maxwell电磁场控制方程,建立磁场与三维等离子体流场耦合数值模拟方法及程序,采用典型算例进行考核。在此基础上,开展不同条件下磁场对再入三维等离子体流场以及气动热环境影响分析。研究表明:建立的高超声速飞行器的等离子体流场与磁场耦合计算方法及程序,其数值模拟结果与文献符合,外加磁场使飞行器头部弓形激波外推,磁场强度越强,激波面外推距离越大;不同磁场强度环境下,流场中温度峰值大小略有变化,变化幅度较小;磁场对绝大部分区域的热流有减缓作用,作用的大小与飞行高度、马赫数以及磁场的配置紧密相关;当前的计算条件下,飞行的高度越高,磁场的作用越明显。     

电离对高超声速热化学非平衡气动热环境的影响

高超声速飞行,激波后高温气体会发生电离,飞行器气动热环境复杂。5组元（N₂,O₂,NO,O,N）、7组元（N₂,O₂,NO,O,N,NO+,e-）和11组元（N₂,O₂,NO,O,N,N₂+,O₂+,NO+,O+,N+,e-）热化学反应采用Gupta化学反应模型,分别数值研究电离作用对高超声速热化学非平衡气动热环境影响。本文分析了不同催化壁面条件下,高超声速热化学非平衡电离流场气动热环境特性。电离作用对激波离体距离和气动力载荷的影响很小。5组元热化学非平衡不考虑电离作用,流场温度和壁面热流密度偏大。11组元热化学平衡强电离流场温度最低;7组元热化学非平衡弱电离流场NO+和e-生成量过低;11组元热化学反应能对热化学非平衡电离流场气动力和热流密度载荷可靠预测。壁面催化作用会增大壁面热流密度,但它对高超声速热化学非平衡电离流场温度和气动力载荷的影响很小。      

真实气体效应下高马赫数内转进气道特性研究

为初步研究高马赫数内转进气道在真实气体效应下的工作特性,首先设计额定工作状态Ma=12的高超声速内转进气道,再结合不同气体模型对其进行数值模拟。研究结果表明:化学非平衡气体在流场结构、工作性能和气动加热方面与热完全气体较为相近,与热化学非平衡气体存在一定差别。离解反应发生在边界层内和低速涡流区内,热化学非平衡气体的离解反应程度比化学非平衡气体大。在隔离段内激波反射处,相比完全气体,化学反应气体的静温降低了2000～2500K。高热流区在上壁面喉道位置与下壁面激波反射点位置附近,温度较高的等温壁面、热化学非平衡气体均可降低壁面热流密度,不同壁面条件对隔离段出口性能参数影响较为明显。真实气体效应、壁面温度对隔离段涡流区的影响较为复杂,有待进一步研究。     

基于支板燃烧室的喷管化学非平衡效应

采用有限体积法全隐式格式和代用燃料C₁₂H₂₃的10组元13步化学反应Arrhenius有限速率模型研究煤油燃料超燃冲压发动机单边膨胀喷管(SERN, Single Expansion Ramp Nozzle)内的化学非平衡流动,通过建立支板燃烧室——喷管模型有效解决了单边膨胀喷管模型的"入口薄层"问题.计算结果表明,整个单边膨胀喷管内,流动呈现化学非平衡效应,喷管入口附近区域尤其显著;非平衡流动喷管性能明显高于冻结流动,随发动机当量混合比 ε增加,非平衡流动的喷管推力系数和升力系数相对冻结流动的百分比增量 δ不断升高,当 ε=0.8时,推力系数百分比增量 δ_CF达到9.41%,升力系数百分比增量 δ_CY达到16.39%,化学非平衡效应对煤油燃料超燃冲压发动机尾喷管性能的影响不可忽略.     

非平衡电弧加热射流光谱学诊断方法：第二部分：氩,氮混合射流

本文是非平衡电弧加热射流光谱学诊断方法实验研究工作的第二部分，给出了氩、氮混合射流的结果，证明亚声速射流中确有非平衡囊兜现象存在。     

提高Spalart-Allmaras湍流模型对分离模拟能力的研究

为了提高Spalart-Allmaras湍流模型对分离的模拟能力,提出了两种改进方法.第一种是根据模型特点而提出的提高对雷诺正应力的模拟精度,第二种是基于湍流特性的理论分析而提出的对模型生成与耗散的关系进行修正.研究结果表明:第一种方法,由于没有明显提高对雷诺正应力的模拟精度,从而并没有提高模型对分离流动的模拟能力;第二种方法,由于改进了模型对强非平衡湍流的模拟能力,从而明显地提高了模型对分离流动的模拟精度.      

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

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

面向稀薄流非线性本构预测的机器学习方法

稀薄非平衡流域内连续介质假设已经失效,主要围绕Boltzmann方程及模型方程对稀薄非平衡流开展理论与计算研究,统一气体动理论格式（UGKS）是其中一种代表性方法。在稀薄非平衡流数值模拟中,Navier-Stokes （N-S）方程连续介质假设已经失效,不能有效描述流场非平衡特征。UGKS方法虽然计算精度高,但速度空间离散导致计算效率低下,多维高速条件下数值计算难以开展。基于数据驱动的思想,在N-S方程与UGKS方法的研究基础上发展出了一种稀薄非平衡流非线性本构关系求解方法（DNCR）。该方法以N-S与UGKS求解器获得的流场数值模拟计算结果作为训练数据集,基于流场特征参数采用极端随机树算法生成机器学习模型,对预测流场中线性黏性应力项与热流项进行非线性修正,并耦合非线性本构关系求解宏观守恒方程得到目标状态稀薄非平衡流动数值解。针对DNCR方法中所采用的机器学习方法-极端随机树模型,通过二维顶盖驱动方腔流算例对高维非线性建模涉及的特征参数选取、参数调优开展了相关验证工作,选取若干典型状态对极端随机树模型的泛化性能开展研究,并评估了相关模型与方法的计算精度与计算效率。     

燃烧多相流的介尺度动理学建模研究进展

基于求解Navier-Stokes方程组的传统计算流体力学已经在诸多领域取得了巨大的成功,但在航空、航天、微流控等领域也遇到了新的瓶颈与挑战。其原因分为2个方面：①物理建模层面的问题;②离散格式带来的数值精度和稳定性问题。微尺度燃烧等一系列燃烧新概念的研究表明,特征更加丰富但以前知之甚少的热力学非平衡行为蕴含着大量待开发的物理功能。物理模型合理和具备相应功能是数值仿真研究的前提;物理建模层面的问题无法通过数值精度的提高来解决。本文从物理建模与复杂物理场分析角度,介绍了非平衡燃烧系统离散玻尔兹曼建模方法（DBM）的研究进展。DBM是非平衡统计物理学粗粒化建模理论在流体力学领域的具体应用之一,是相空间描述方法在离散玻尔兹曼方程形式下的进一步发展。它选取一个视角,研究系统的一组动理学性质,因而要求描述这组性质的动理学矩在模型简化中保值;以该组动理学矩的独立分量为基,构建相空间,使用该相空间和其子空间来描述系统的非平衡行为特征;研究视角和建模精度随着研究推进而调整。借助DBM可以研究反应过程中不同自由度内能之间的不平衡和相互转换等Navier-Stokes模型无法模拟的动理学过程。在内爆和外爆过程中,几何汇聚与发散效应等效于一个"外场力",在球心处系统始终处于热力学平衡态;在冯·纽曼压强峰处,系统不是偏离平衡最远,而是在平衡态附近;在冯·纽曼峰后反应区以外,Chapmann-Jouguet理论值、Zeldovich-Neumann-Doering （ZND）理论值和DBM结果相互验证;在反应区内DBM结果与ZND结果一致;在冯·纽曼压强峰前的压缩阶段,DBM模拟结果在物理上更合理。在冲击压缩过程中,相对于其他自由度,压缩波所在自由度上的内能先增加,因而这一自由度上的内能总是朝着正向偏离其平衡态值,而横向自由度上的内能总是朝着逆向偏离其平衡态值。在二流体模型视角下,反应物和产物朝着相反的方向偏离热力学平衡态。