三维自适应非结构网格的Euler方程解

 将 Ausm+ 迎风格式应用于三维非结构网格中求解 Euler方程。对单元变量进行重构以获得空间高阶精度,对时间域采用多步龙格库塔法推进,并采用了当地时间步长和隐式残差光顺技术来加速收敛。采用多点择优推进阵面法生成复杂曲面的三角形网格,利用推进阵面法生成四面体网格。采用网格自适应技术对网格进行局部加密,以减少总体网格数目,从而提高计算效率。最后给出了绕 ONERA M6机翼的跨音速流动及绕麻雀 导弹的超音速流动算例,结果表明了本方法的有效性。     

一种通量加权型紧致格式

将通量加权的思想引入到紧致格式中,构造了一个传统方法与紧致格式混合组成的加权型差分格式.利用该格式与二阶TVD格式分别计算了一维方波、组合波问题和一些黎曼问题,如Sod问题和Shu问题,以及一维定常激波问题.计算结果的比较表明加权格式无论在捕捉各种间断,还是在分辨各种复杂波系上,都具有较大的优势,并与精确解非常吻合.      

稀薄流航天器鼻锥迎风凹腔气动力和气动热性能研究

本文采用直接模拟蒙特卡罗(DSMC)方法,对高超声速稀薄流中航天器鼻锥迎风凹腔气动力与气动热性能进行了数值研究。得到了鼻锥外壁面、凹腔侧壁面以及凹腔底面的热流密度分布,分析了不同凹腔深宽比对鼻锥冷却效率以及凹腔腔体内气体参数的影响;以深宽比为1的凹腔为基准,研究了凹腔唇口钝化半径对航天器气动热与气动力的影响。数值结果表明,高超声速稀薄流中迎风凹腔能够降低鼻锥外壁面的热流密度;当凹腔深宽比达到1之后,凹腔侧壁面热流变化趋于一致,热流密度最低点的轴向位置不随深宽比改变,且凹腔底部热流很小;凹腔近底部气体均由稀薄流转化为连续流,腔内气体压力不断振荡;唇口钝化没有明显优势,虽然可以降低鼻锥峰值热流,但是会带来严重的气动力性能下降。     

改进的二维三阶半离散中心迎风格式

对二维三阶半离散中心迎风格式中的权函数给出了简化改进. 在保持格式精度的基础上, 改进后的权函数在二维情况下具有更加简单直接的结构而且严格非负. 该改进方法得到的格式仍然具有半离散中心迎风格式的优点, 同时保持了重构函数的非振性. 时间离散采用保持强稳定性的三阶Runge-Kutta方法, 并利用四阶Lax-Wendroff (L-W) 格式计算磁流体算例中的磁场散度. 用该修正格式计算了二维磁流体数值算例, 得到高精度无振荡的结果, 验证了此方法的有效性.      

飞行器鼻锥凹腔-发散组合冷却数值模拟

尖锐鼻锥冷却方案是可复用式航天飞行器研究领域一个十分重要的课题。传统发散冷却虽然可以有效降低鼻锥结构温度,但是由于驻点外极高的热流、压力,会出现驻点冷却效果差的问题。迎风凹腔结构是一种针对鼻锥驻点区域的减阻防热方案,尖锐唇口的分流作用可以使附近压力、热流降低。因此,提出一种新型冷却结构——凹腔-发散组合冷却,利用迎风凹腔结构对驻点的强化冷却解决发散冷却中驻点难以冷却的问题。以楔形鼻锥为物理模型,对发散冷却、迎风凹腔结构和凹腔-发散冷却3种冷却结构进行数值模拟,并和无冷却的纯鼻锥结构进行对比。结果表明,与传统发散冷却相比,使用凹腔-发散组合冷却可以使结构温度峰值下降16.8%;与没有冷却的纯鼻锥模型相比,鼻锥头部圆弧段表面平均温度降幅可达64%,证实了这种新型冷却结构的可行性和高效性。     

几种典型迎风格式的分析与比较

重点分析了AUSM 格式的构造方法,同时将其与Van Leer, Roe格式等构造方法进行比较和分析,通过求解一维激波管问题,以及二维和三维跨音速流场的计算比较了AUSM 格式,VanLeer格式,Roe格式和中心格式的数值模拟结果和网格适应性,以及在经典四步Runge-Kutta显式时间推进和LU-SGS隐式时间推进下上述四种格式的计算效率.获取了对于基础研究和工程应用很有意义的结论.     

非结构网格上的TV-van Leer混合迎风格式

为了改进迎风型对流通量格式的激波捕捉计算稳定性,对迎风格式产生数值激波不稳定现象的原因进行分析,并且在分析结果基础上构造新型迎风型格式。为此,首先设计了数值算例,分析非结构网格上迎风格式出现数值激波不稳定现象的主要影响因素。结果表明,在一定条件下,网格和流动状态共同作用导致产生不合理的数值结果,并影响计算的稳定。其中,网格间的动量通量对数值激波不稳定现象影响显著。因此,基于TV格式构造了一种自适应调节动量通量耗散项的混合格式,从而在激波附近引入耗散,抑制扰动的产生和传播。数值计算结果显示新型混合格式在非结构网格离散条件下能够稳定地捕捉激波,同时混合格式保持了TV格式原有的黏性流动计算分辨率。     

贴体曲线坐标系中恢复函数的一种构造方法及应用

本文在结构网格下，从三点迎风紧致逼近出发提出了一种适合于有限体积离散的恢复函数生成办法，在光滑区它具有三阶精度，并且在捕捉激波时有较高的激波分辨率。典型的几个内流算例表明：用该方法得到的恢复函数去计算数值能量时能得到与实验数据较接近的数值解。     

基于迎风格式数值模拟三角翼涡流场

本文应用基于有限体积方法的高精度迎风格式模拟了三角翼在大迎角条件下的绕流流场。选用不同几何外形的三角翼，由计算结果可以清楚的看到前缘自由剪切层的发展，涡的形成及其复杂的相互作用。     

凹腔尺寸对迎风凹腔与逆向喷流组合热防护系统性能的影响

针对高超声速飞行器热防护系统(TPS)的设计,对迎风凹腔与逆向喷流组合热防护系统展开研究.在数值方法实验验证的基础上,通过求解Navier-Stokes方程得到了带组合热防护系统的鼻锥的流场结构以及壁面热流分布.验证了组合热防护系统的有效性.在逆向喷流条件不变的情况下,进一步研究了凹腔的尺寸变化对其防热能力的影响.研究发现:凹腔的直径越小,深度越深,气动加热值越低.自由来流与逆向喷流形成的回流区在减少鼻锥的气动加热上起到关键的作用.相对于凹腔深度的变化,鼻锥壁面的气动加热更敏感于凹腔直径的变化.