高超声速全机外形气动加热与结构传热快速计算方法

发展了一种无黏流场解与工程计算方法相结合的高超声速全机外形气动加热与结构传热快速计算方法。该计算方法结合了三维块结构网格无黏流场数值计算技术可处理复杂外形流动的优点与工程计算方法效率高的特点,将气动热的计算简化为绕飞行器的无黏外流(边界层以外)数值解和边界层内热流求解两个部分,同时耦合了防热结构传热计算模型、高温化学非平衡热效应估算方法以及弹道状态动态插值方法,可用于快速计算与分析三维复杂外形高超声速飞行器在弹道飞行状态下全机热环境参数、防热结构内温度场等随飞行时间的变化特性。以RAM-CⅡ、类Ⅹ-37B等典型高超声速飞行器为研究对象,在设定的飞行条件及热防护方案下,进行了气动加热与结构传热问题的求解,给出了全机表面热流密度与防热结构材料温度的时变特性。结果对比表明,所发展的方法具有快速、高效的特点,且计算精度可满足工程设计初期选型需求,可为高超声速飞行器的热防护系统初期设计及热环境特性快速计算分析提供技术支持。     

钝体外形气动加热与结构传热一体化数值模拟

为了提高气动加热与热防护结构传热多物理场数值模拟的稳态计算效率与计算精度,发展了一种基于有限体积法的气动加热与结构传热一体化数值计算方法。该方法将高速流场与结构温度场统一到同一物理场,基于统一的控制方程组,采用基于LU-SGS隐式时间迭代和自适应时间步长的有限体积方法进行求解,避开了传统气动加热与结构传热耦合求解方法在时间域内的所需繁琐数据交替迭代策略。对二维/三维钝体进行一体化数值计算分析,计算结果表明:二维钝体非稳态下,得到2s时圆管驻点温度最高达到390.2K,驻点热流密度和结构温度与参考文献和实验值吻合较好,证明了方法的可靠性和可行性。同时分析了三维钝体应用算例的流-固-热稳态计算特征,计算得到稳态时钝头体结构外壁表面最高温度达到535.6K,表明一体化计算方法可用于长航时飞行条件下的气动加热-结构传热多物理场耦合计算分析,为高速飞行器热防护结构设计与选材提供一定的理论与技术支持。     

高超声速飞行器热环境与结构传热的多场耦合数值研究

为了准确预测高超声速飞行器面临的严峻气动热/力环境以及结构的热力响应,发展了高超声速流动与结构传热耦合框架。采用分区求解方法,通过耦合界面的实时数据传递,实现了基于Navier-Stokes方程的高超声速化学非平衡计算流体力学（CFD）求解器与结构的热力全耦合有限元法（FEM）求解器的多场耦合计算,建立了高超声速飞行器的多场耦合数值分析方法。首先对经典高超声速圆柱绕流实验进行了耦合计算,结果与实验值吻合良好。然后针对典型的超高温陶瓷（UHTC）材料的耦合传热问题进行了数值研究,考虑热传导效应对气动热环境和结构热响应预测的影响,结果表明对于复杂外形且热导率相对较高的UHTC材料,结构内部热传导对热环境和表面温度分布的影响不可忽略。最后针对UHTC材料热物性（比热和热导率）非线性对高超声速流动传热过程的影响进行了研究,结果表明当比热和热导率处于合理的误差范围内时,材料表面温度响应对其变化并不敏感。     

高超声速飞行器多层复杂热防护结构气-固耦合快速热分析方法

为了实现高超声速飞行器多种复杂结构热防护系统的气-固耦合快速热分析,采用热网络法建立非稳态等效一维传热模型;对于具有弧度的热防护结构,提出了驻点和翼前缘热阻等效计算方法,并给出了修正计算公式;结合气动热环境工程算法,实现了对任意多层复杂防热结构外部气动加热与内部结构传热的快速耦合分析。分别对钝锥气动加热和高超声速二维圆管气-固耦合传热问题进行了模拟,得到了与实验符合较好的结果,且计算效率很高;并对Micro-X验证机的全过程进行了耦合热分析,结果表明多层防热结构具有很好的防热效果,显著降低了结构内部温度。和传统耦合算法相比,此算法可快速有效地分析模拟气-固耦合问题,满足高超声速飞行器热防护系统初始设计阶段的使用要求。     

高速飞行器综合热管理方案快速仿真平台

高速飞行器面临严峻的外部气动热和内部设备热环境,以热防护系统和热管理系统组成的综合热管 理方案决定该类飞行器总体方案的闭合性,因此需要在飞行器初步方案设计阶段对综合热管理系统设计方案 进行快速仿真计算和可行性评估。以电子设备舱为例,采用模块化设计和一维简化方法,构建一个将飞行器热 源、传热和热控等热环境因素进行模块化封装的快速仿真平台;通过仿真算例对该平台进行验证。结果表明: 该平台可以快速实现飞行器电子设备舱综合热管理方案的模型搭建、计算与评估。     

高超声速多体分离过程气动加热计算技术

发展了无粘外流解与气动热工程方法相结合的用于全机外形复杂流动气动加热计算技术,并将其与动态多体分离过程流场解相结合,提出了一种高超声速飞行器多体分离过程气动加热特性计算技术,然后耦合结构传热,实现了飞行器热防护系统中结构材料温度分布特性的数值模拟。采用本文方法,针对假定高超声速多体分离布局及飞行条件开展了分离过程中的结构温度分布特性与热流密度分布特性计算,分析并给出了组合体分离过程中部件干扰对气动热的影响特性。结果表明本文发展的计算技术可为高超声速飞行器的防热方案设计、气动热特性定性分析及热防护系统设计等提供技术支持。     

“Ⅳ型”激波干扰中流-热-固耦合问题一体化计算分析

针对高超声速进气道前缘"Ⅳ型"激波干扰产生的气动加热与结构传热多物理场耦合计算问题，发展了一种基于有限体积法的流-热-固一体化计算方法。该方法采用一体化控制方程组统一离散求解外部高速流场与内部结构温度场，规避了传统分区耦合算法在时间域内交替迭代的繁琐数据交换策略。另外，提出一种新的双温阻模型计算流-固交界面的物性参数以保证计算准确性，采用LU-SGS隐式时间迭代和自适应时间步长以提高计算效率。采用经典高超声速二维圆管流-热-固耦合算例对该一体化方法进行验证，计算结果与试验值和参考文献数据吻合较好，证明了该方法的可靠性和正确性。利用一体化方法对高超声速前缘"Ⅳ型"激波干扰流-热-固耦合问题进行定常/非定常计算与分析，给出了温度与热流的时变特性，计算结果表明，激波干扰作用产生的超声速"喷流"不断冲击壁面，使得壁面最大压力系数增大约9倍，壁面最大热流增大约4.7倍，给高速飞行器的热防护设计与选材带来严峻挑战。同时，也表明了一体化计算方法可以较好地用于长航时飞行条件下与复杂飞行环境下的高超声速热防护系统的热环境特性分析与综合性能评估。     

一种耐高温多层热防护结构的优化设计与性能

为解决高速飞行器飞行过程中剧烈的气动加热问题，以“高温防热层+隔热缓冲层+核心隔热层”顺序设计的一体化多层热防护结构的传热过程为研究对象，建立了高温环境下热防护结构内部一维非稳态导热-辐射耦合传热模型，通过数值模拟计算得到了高温环境下热防护结构各层的温度分布。利用不同热防护材料的隔热性能差异，针对构建的热防护结构，提出了在满足一定约束条件下，以轻质多层热防护结构总质量和总厚度为目标函数的优化设计方案，得到了多层结构的最优几何参数，并通过实验考核了优化后热防护结构的防隔热性能。实验表明：该结构可耐受1 473 K的高温1 800 s而背温不超过370 K。     

基于TLBM-FVM耦合的飞行器舱内热环境跨尺度预测方法

舱内热环境的有效预测是优化飞行器热控与防热设计、减小系统冗余并保障飞行器热安全的重要基础。在国家数值风洞（NNW）工程支持下,针对目前舱内热环境多尺度、多机制复合传热特点及其数值预测面临的精度与效率提升难题,发展了多区域协同推进的时空耦合模型及流/固界面的自适应分辨率识别算法,建立了基于热格子Boltzmann方法（TLBM）与有限体积法（FVM）相互耦合的舱内复合传热跨尺度预测方法,开展了典型飞行器仪器舱的综合热分析,验证了耦合方法的计算精度及效率。研究表明,相关方法可实现舱内热环境的局部精细化与整体大规模的协同模拟,用于开展整体自然对流与设备局部热量传递的多尺度数值模拟,掌握不同环境参数对舱内热质传递过程的影响规律,从而为飞行器热防护/热管理一体化设计提供重要技术支撑。     

不同热边界对超音速机翼气动热数值分析的影响

飞行器在超音速飞行时受到的气动加热效应给结构强度及热防护设计带来极大影响，且真实状态下的气动热环境需要考虑外流场与结构的耦合及内壁面边界条件的影响。采用S-A湍流模型求解Navier-Stokes方程，通过流场与固体壁面交界处的信息传递，实现外流场与结构场的耦合数值分析。针对三种不同翼型的超音速绕流气动加热进行耦合数值研究，对比翼型内壁面在不同热边界条件下的气动热效应，结果表明：不同翼型具有与气动力相似的的气动热效应；内壁面考虑对流换热的边界条件最接近真实；考虑机翼燃油箱满油时，三维机翼前缘驻点处热流密度最高可达4200w/m2。     

突然启动流动问题:从不可压到高超声速流动

突然启动流动问题在气动弹性、扑翼飞行和机动飞行中有重要应用价值。鉴于此,对突然启动流动问题涉及的流动与升力演化机制和理论分析方法进行统一介绍并指出目前的理论空白。比较性地介绍了升力随时间变化的原因、预测升力演化的理论方法和主要变化规律。从中发现,不可压缩和可压缩突然启动问题存在各自独立的研究方法与结论,因此进行统一介绍与分析对建立二者之间的关联有指导意义。同时指出,大迎角可压缩突然启动问题尚无理论分析方法和研究结果。作为一项补充研究,采用数值计算发现一个之前尚未报道的现象,即升力系数首先从较低值快速增加至一个峰值,接着快速下降,趋于定常值。该文综述介绍的方法可用于分析突然启动问题升力演化规律。     

径流及斜流压气机任意曲面叶型长短叶片的造型设计方法

发展了一套径流及斜流式叶轮机的任意曲面叶型长短叶片的造型方法,并在离心式压气机和斜流式压气机的设计中做了具体的应用。新发展了令叶片的厚度按三次样条函数沿流向分布的一种设计思想,得到了比较理想的叶片。同时在径流及斜流式压气机中使用分流叶片,引入了长短叶片的造型设计方法。通过CFD流场检验验证了这种叶片造型方法的可操作性和适用性。这为径流式以及斜流式叶轮机转子、静子的进一步气动优化设计奠定了基础。      

开式空腔气动声学特性及其流动控制方法

在高速风洞中对空腔流场气动声学特性进行了试验研究,采用剪切层扰流法对流场进行流动控制,空腔长深比为4.1。通过对空腔流场的脉动压力试验结果分析,研究了亚、跨声速条件下开式空腔流场的气动声学特性及气动噪声抑制效果,对开式空腔流场气动噪声形成机制及流动控制机理进行了分析。试验结果表明:开式空腔流场气动声学环境恶劣,最大总声压级(OSPL)高达177dB;开式空腔流场存在强烈的自持振荡,声压频谱曲线上存在多个不同模态的单调声;亚声速条件下,采用剪切层扰动法进行流动控制可导致空腔流场气动声学环境更加恶劣;跨声速时,采用剪切层扰动法进行流动控制使空腔流场气动声学环境明显改善。     

高比转速斜流叶轮流场分析

以两个高负荷、高比转速离心叶轮与斜流叶轮为例,采用数值模拟的方法,在级环境下分析了两种叶轮的总体性能,以及内部流动机理.计算结果表明在高比转速情况下,斜流叶轮有利于减弱叶尖泄漏损失,改善叶轮出口均匀性,提高叶轮效率和级效率.      

直升机旋翼下洗气流对排气喷流的影响

采用数值模拟方法,对旋翼下洗气流作用下的排气喷流流动特征进行了研究,分析了旋翼下洗气流速度和排气喷口方向对排气喷流流动以及排气系统引射能力的影响.研究结果表明:排气喷流受到旋翼下洗气流的作用而发生明显的向后机身下方以及旋翼转动方向的偏转,其偏转程度随旋翼下洗气流速度的增大而加剧;当排气喷口向上排气时,排气喷流在旋翼下洗气流作用下的偏转能够形成对后机身表面的撞击,排气系统的引射能力有微弱的降低,引射系数减小约0.01;而当排气喷口斜向上或侧向时,排气喷流对后机身未形成撞击,引射能力得到了一定程度的提升,引射系数最大增大0.12.      

燃气非平衡流再生冷却流动传热数值模拟

为准确预测液体火箭发动机推力室身部再生冷却换热状况,采用数值模拟方法,对燃气、推力室壁和超临界气氢进行三维流动和换热耦合计算。采用6组分9步反应动力学模型计算燃气的非平衡化学反应,采用DO模型计算燃气辐射换热,考虑超临界气氢物性随温度和压力的变化。获得了室壁温度场、燃气及冷却剂流场。结果表明,Redlich-Kwong方程、Peng-Robinson方程、Lucas法、TRAPP法能分别准确计算超临界氢的密度、定压比热容、粘度、导热系数,采用燃气非平衡流计算所得流场值更符合实际情况。     

固冲补燃室内流场冷流模拟及显示

为了获得对固冲补燃室内混台的认识,以气体为介质进行了补燃室内流场的模拟,通过在气流中播入粒子,利用激光片光成像获取了补燃室内的混合图像对获取的不稳定的连续图像进行处理,获得了可用的稳态混合图像一另外,对其进行伪彩色处理获得了令人满意的结果,与测温法获得的数据进行了比较,结果表明,处理结果是可靠的,这种冷流模拟显示的方法可以用于同冲补燃室内的混合研究。     

Recent progress of machine learning in flow modeling and active flow control

In terms of multiple temporal and spatial scales, massive data from experiments, flow field measurements, and high-fidelity numerical simulations have greatly promoted the rapid development of fluid mechanics. Machine Learning(ML) provides a wealth of analysis methods to extract potential information from a large amount of data for in-depth understanding of the underlying flow mechanism or for further applications. Furthermore, machine learning algorithms can enhance flow information and automat...     

Multi-phase flow effect on SRM nozzle flow field and thermal protection materials

Multi-phase flow effect generated from the combustion of aluminum based composite propellant was performed on the thermal protection material of solid rocket motor(SRM) nozzle.Injection of alumina(Al₂O₃) particles from 5% to 10% was tried on SRM nozzle flow field to see the influence of multiphase flow on heat transfer computations.A coupled,time resolved CFD(computational fluid dynamics) approach was adopted to solve the conjugate problem of multi-phase fluid flow and heat transfer in the solid rocket motor nozzle.The governing equations are discretized by using the finite volume method.Spalart-Allmaras(S-A) turbulence model was employed.The computation was executed on the different models selected for the analysis to validate the temperature variation in the throat inserts and baking material of SRM nozzle.Comparison for temperatures variations were also carried out at different expansion ratios of nozzle.This paper also characterized the advanced SRM nozzle composites material for their high thermo stability and their high thermo mechanical capabilities to make it more reliable simpler and lighter.      

反推气流对发动机进口流场影响的数值研究

通过求解三维雷诺平均的Navier-Stokes(N-S)方程,获得了某大型运输机降落滑跑过程中,在不同滑跑速度下的反推气流扰流流场细节.单台发动机反推气流扰流流场的计算结果表明,在任何状态下,反推气流都不会被发动机重新吸入.飞机/发动机一体化计算结果表明:随着相对来流马赫数的减小,反推气流被发动机重新吸入的可能性不断增大.当相对来流马赫数减小到0.1时,反推气流会被外侧发动机重新吸入,此时,发动机进口截面出现了明显的流场畸变,周向稳态总压畸变指数增加明显.当相对来流马赫数减小到0.05时,两台发动机都会吸入反推气流.当相对来流马赫数减小到0时,反推气流没有被重新吸入发动机,但是反推气流会干扰吸入发动机的自由流,降低自由流的总压,从而也会造成发动机进口的流场畸变.