叶轮机械三维不可压Euler方程的直接解法

本文用近似因式分解交替方向隐式格式直接求解叶轮机械中原参数三维不可压Euler方程组。在连续方程中引入“拟压缩性”之后,将原来类型不确定的控制方程变成双曲型方程,给定初边值之后用时间推进法求得稳定解,对一个单级压气机转子内部流场进行了计算,结果与实验值比较符合。     

求解Euler方程的无网格与有限体积混合算法研究

研究了无网格方法与有限体积法相结合的、求解Euler方程的混合算法.与单一的无网格方法相比较,由于在大部分计算区域采用了有限体积法,使得发展的算法在运算效率上能与有限体积法相当;同时在物体附近嵌入了无网格区,使得在处理几何外形上更具灵活性.有关算法涉及的无网格与网格区域搭接处理,通过在交界面处引入辅助卫星点和网格单元构成边界信息,实现了区域间的流动信息传递.Euler方程的空间导数分别在两区域用有限体积法和无网格法离散近似,时间方向都采用四步显式Runge-Kutta格式推进求解,数值模拟了喷管内流和绕翼型外流,并分别与整体有限体积法和整体无网格方法进行了比较.算例展示出,提出的混合算法能有效捕捉激波间断,且两区域等值线过渡光滑,算法效率如预期与有限体积法相当,表明该方法是可行的.     

冷却环带喷注结构对煤油超临界液膜的影响研究

针对液氧/煤油补燃发动机液膜冷却过程,建立了超临界条件下的液膜冷却模型,分析了冷却环带喷注结构对局部流动和冷却效果的影响.在超临界条件下,煤油和周围燃气为同种流体,用同一组方程来描述其流动与传热过程,对3种典型冷却环喷注结构的流动进行了数值模拟.结果表明:冷却环带的出口角度、台阶结构和喷射角对局部流动和传热有显著影响,出口角度小于90°的台阶式结构可减小局部回流从而有利于推力室的热防护.      

含内热源球壳的复合传热特性及其表面温度均匀性研究

应用计算流体动力学方法,对带内热源的空腔球壳所涉及的含导热、自然对流和辐射换热的耦合传热问题进行了数值建模,得到了球壳空腔内、外侧流体的温度场、速度场。在数值模拟基础上,提出并试验验证了能有效提高球壳外表面温度均匀性的电加热绳方案。     

感性耦合夹层等离子体隐身天线罩电磁散射分析

设计了一种石英夹层感性耦合等离子体（ICP）隐身天线罩模型,采用有限元与Z变换时域有限差分（ZT-FDTD）联合仿真的方法,建立了ICP放电的流体模型,得到不同气压及功率下与电磁散射相关的电子密度空间分布,在此基础上建立Z变换时域有限差分法模型,对石英夹层等离子体隐身天线罩的宽频段后向散射进行计算,同时利用微波干涉法及XFDTD软件对算法及程序进行验证。结果表明：感性耦合等离子体可产生较高密度等离子体,能有效实现雷达散射截面（RCS）的减缩,在气压为2 Pa时,碰撞衰减较弱,等离子体密度分布较均匀,衰减带宽集中在等离子体振荡频率附近,功率增加会使衰减峰值向高频方向移动,气压为20 Pa时,碰撞衰减增强,且等离子体密度分布有较大梯度,衰减带宽有效增加,同时RCS曲线的波动特性加强。     

超临界压力下碳氢燃料裂解与流动传热模拟的快速算法

为了快速预测发动机冷却通道内碳氢燃料在考虑热裂解时的流动传热特性,基于正癸烷热裂解反应机理,建立了一套模拟正癸烷裂解吸热和超临界压力传热现象的快速算法。采用三维物性库查表算法计算裂解反应混合物的热物性,同时对组分输运方程进行简化,简化后只需求解1个组分输运方程。通过与现有的实验和数值结果进行比较,检验了快速算法的计算效率和可靠性。结果表明,快速算法与求解全组分输运方程的算法精度相当,但计算效率提升了约20倍。最后采用该算法对三维矩形截面冷却通道内的超临界压力正癸烷裂解与传热过程进行数值模拟,进一步考察了本文快速算法的计算精度及其工程应用价值。     

剪切稀化非牛顿射流撞击液膜破碎直接数值模拟

液态射流撞击是液体火箭推进系统中广泛采用的一种燃料雾化方法,其破碎特征直接影响燃料最终的掺混及燃烧效率。采用直接数值模拟（DNS）工具,研究了低雷诺数（Re_l=41）和中等韦伯数（We_l=163）条件下剪切稀化非牛顿射流撞击液膜破碎的问题,着重分析了对角液膜的三维结构、破碎特征和非牛顿特性等。研究结果表明：在所研究的射流参数下,该非牛顿撞击液膜破碎属于Open Rim类别,破碎过程具有三维特性并伴随液丝与边缘的融合、液丝向液滴的转变等时域流动特征。液体的总表面积随时间不断增长,但单位表面积随液膜破碎的发生而下降,液膜扩张半角随时间逐渐增加并趋于恒定值43°,而后部液膜的长度几乎不随时间发生变化。此外,撞击液膜表现出明显的剪切稀化特性,液体内部最低黏性系数仅为零剪切黏性系数的1/5。     

高超声速飞行器宽速域翼型/机翼设计与分析

对于现代高超声速飞行器的设计而言,除了需要保证高超声速的性能外,还必须兼顾满足工程需求的亚跨超声速特性。首先,采用雷诺平均Navier-Stokes（RANS）方程流动求解器,结合基于Kriging模型的代理优化算法,开展了高超声速飞行器宽速域翼型的优化设计研究,设计出了一种下表面具有双"S"形特征的新翼型。综合性能评估结果表明,该翼型相比于常规的高超声速翼型,在跨声速和高超声速下具有更加优良的气动特性;其跨声速状态下的升阻比达到78.9,高超声速状态下的升阻比达到5.94,能够实现宽速域内良好的综合气动性能。其次,开展了仿德国"桑格尔号"（SANGER）空天飞机运载机机翼的气动特性研究,对配置宽速域翼型与常规高超声速翼型的机翼进行了气动力特性综合对比分析。结果表明,配置新翼型的机翼在宽速域范围内整体气动性能更优,说明所设计的宽速域翼型在三维机翼上也具有一定的实用价值。     

基于非结构网格流场超大规模并行计算

大规模并行的计算流体力学已成为现代航空工业研发的核心手段之一。基于非结构混合网格和有限体积法,发展了适用于工业级复杂外形气动计算的并行流动数值模拟方法。文中首先介绍了紧致数值离散格式、基于Metis的分布式多核系统网格分区技术、并行边界虚拟单元技术和MPI并行实现等相关算法。采用网格量相对较小的旋成体构型绕流模型对比分析多核并行计算结果与单核计算结果以验证并行计算的正确性,比较了不同并行规模下并行效率和残差收敛情况。然后通过对上亿网格单元的运输机复杂构型绕流模拟,开展并行效率的测试,结果表明,本文方法并行加速性能高,直到多达18816核并行效率都保持在80%以上。     

基于深度神经网络的粒子图像测速算法

粒子图像测速(PIV)作为一种流体力学实验技术,能够从流体图像中获取全局、定量的速度场信息。随着人工智能技术的发展,设计用于粒子图像测速的深度学习技术具有广泛的应用前景和研究价值。借鉴在计算机视觉领域用于运动估计的光流神经网络,采用人工合成的粒子图像数据集进行监督学习训练,从而获得适用于流体运动估计的深度神经网络模型,并且能够高效地提供单像素级别分辨率的速度场。文中采用人工合成的湍流流场粒子图像进行初步实验评估,并讨论PIV神经网络的隐藏层输出和内在原理,同时将训练而成的深度神经网络模型与传统的相关分析法、光流法对比;随后进行射流流场测速实验,验证深度神经网络PIV的实用性。实验结果表明,文中提出的基于深度神经网络的粒子图像测速在精度、分辨率、计算效率上具有优势。