吸气式高超声速飞行器多参数灵敏度分析

为实现吸气式高超声速飞行器多参数情况下的灵敏度分析及参数分类,降低设计的复杂度,在吸气式高超声速飞行器参数化建模的基础上,首先采用正交试验设计生成样本,再通过计算流体力学(CFD)进行高精度气动力性能计算,最后运用方差分析法进行气动性能的参数灵敏度分析,在运用较少的试验样本点的情况下,即可完成多参数、复杂构型条件下气动性能的参数灵敏度分析。结果表明,该方法可以正确地分析出参数对飞行器气动性能的敏感程度,得到参数对气动性能的影响规律。同时,通过灵敏度分析的计算样本,还可以初步选出气动性能较优的飞行器构型,为地面试验和优化设计提供参考。     

射流对高超声速进气道起动性能的影响

为探索利用射流技术降低进气道起动马赫数的可行性,对二元高超声速进气道二维流场进行了数值模拟,通过对比不同工况的流场结构、流量系数及总压恢复系数,分析了射流对高超声速进气道的作用效果,并研究了射流速度、压强及倾角对进气道起动性能的影响。分析结果表明:施加射流,激波与进气道边界层原有干扰形式发生改变,是降低进气道起动马赫数的主要原因。研究还表明,增大射流速度利于提高控制效果,但持续增大射流速度,会造成隔离段反压升高,并且这一现象与射流压强相关,降低射流压强能使进气道起动的射流速度区间扩大,同时在不同射流倾角下,上述规律表现一致。该研究揭示了进气道起动能力随射流参数变化的系统性规律,可用于指导工程设计及优化。     

基于两级直角网格结构的三维DSMC算法研究

结合直角网格计算的高效率以及表面非结构网格对飞行器几何外形的精细描述,本文建立了基于两级直角坐标网格结构的DSMC数值模拟算法。一级网格为计算区域离散后的背景网格,二级网格为独立的碰撞网格和取样网格,是在背景网格的基础上根据流场性质的变化经过自适应后获得的亚网格,碰撞分子对在亚网格内选取,以保证选取的分子是"最相邻的"。计算的Orion飞船外形从自由分子流到稀薄过渡流的气动力系数,与DS3V和DAC软件计算结果的一致性非常好,验证了本文算法的可靠性和有效性。同时对过渡区的流动特征和表面气动参数的分布规律进行了分析。     

适用于混合网格的并行GMRES+LU-SGS方法

给出了一种适用于混合网格的并行无矩阵GMRES+LU-SGS隐式时间格式。首先采用LU-SGS方法迭代若干步以获得一个合适的初场,然后切换到GMRES方法在每一时间步内近似求解,并将LU-SGS方法作为其预处理器。为加速收敛,将CFL数随着残差的降低逐步放大;为减少存储量和计算量,通量Jacobian采用无矩阵处理。在保证与串行执行一致的前提下,采用基于共享内存的OpenMP方法实现了并行计算,并通过对网格的分组避免了内存争夺。算例验证表明,方法极大地提高了计算收敛效率,并行结果与串行结果完全一致,计算结果与实验结果吻合较好。     

Roe格式中不同类型熵修正性能分析

针对Roe格式中较为流行的3类熵修正:Muller型、Harten-Yee型和Harten-Hyman型熵修正,从理论分析入手, 辅以Euler方程为控制方程的激波管问题,前台阶流动和运动激波的双马赫反射3个数值实验,对不同熵修正的性能做了深入研究,得出如下结论: Muller型和Harten-Yee型熵修正方法作用于激波和非物理的膨胀激波2种情况:激波情况下引入的数值耗散,能根本改善"Carbuncle"现象,膨胀激波时其数值粘性也足以使膨胀激波得以耗散;而类Harten-Hyman型熵修正只对膨胀过程起修正作用,对激波情况无效,不能改善"Carbuncle"现象,该类熵修正的数值粘性较小,不足以使膨胀过程正确求解;直接采用 δ 值代替特征值的特征值修正方法效果好于传统的特征值修正方法.     

层流平板摩擦阻力的数值计算

为提高层流摩擦阻力数值计算的精度,采用数值求解N-S方程的方法,对层流平板摩擦阻力进行了计算.计算中采用4种先进上风格式,同时配合以14种不同的壁面法向网格分布和5种常用限制器,并将计算结果与布拉修斯解进行对比.结果表明,van Leer的FVS(Flux Vector Splitting)格式数值耗散大,对网格分布和限制器敏感,不适合用于摩擦阻力计算.限制器和计算格式的粘性分辨率越高,摩擦阻力的计算精度越高.在计算量允许的情况下,可以减小壁面法向第1层网格高度,以提高计算精度,但应同时兼顾网格长宽比.对于摩擦阻力计算,边界层内应至少布置4个网格点.      

双三角翼前缘剖面形状对涡运动的影响

采用数值模拟和理论分析相结合的方法,研究了前缘剖面形状对双三角翼涡运动的影响,分析了前缘剖面形状对三角翼、双三角翼涡运动影响的不同机理 :对三角翼,尖前缘可以形成组织最好的涡结构,但对于双三角翼,圆前缘生成的旋涡结构较靠近翼面,涡结构紧密,诱导能力较强,可以形成有利的涡涡干扰,使内翼涡通过剪切层向外翼涡输入涡量更加容易,合并涡变得更加稳定,推迟了涡破裂,而且由于涡较靠近翼面,因而可以产生较高的非线性涡升力,这同传统的认识是不一致的。     

气动热CFD计算的格式效应研究

通过对N-S方程的求解计算了二维圆柱绕流及三维绕钝双锥高超声速流动过程中的气动热问题.对于求解方程,在空间离散的格式上,选用了3种上风格式——Roe的FDS,Van Leer的FVS和AUSM+格式以及一种中心差分格式,将各种格式所得的热流结果与实验结果进行了比较.对于AUSM+格式,还应用了minmod limiter 和 双 minmod limiter 两种限制器,并对结果进行了对比.      

翼型的气动最优化设计方法和反设计方法

本文将数值优化方法同计算流体动力学(CFD)相结合,形成两种跨声速翼型的气动设计方法,即:最优化方法和反设计方法.采用求解Euler方程的有限体积法计算流场,通过结合优化算法,从正反两个方向对跨声速翼型进行气动优化设计.实例证明它们是翼型气动设计的有效方法,有较高的工程应用价值.     

两类激波捕捉格式的性能分析

考虑了两类典型的激波捕捉格式:特征形式的MUSCL(Monotone Upstream-centred Schemes for Conservation Laws)格式和WENO(Weighted Essentially Non-Oscillatory)格式.MUSCL格式在作特征变换时使用了局部线性化的思想,并且针对波系的性质施加相应的限制器;通过逐维重构实现有限体积法的WENO格式.针对一维、二维和三维Euler系统进行数值实验.在一维和二维的情况下,特征形式的MUSCL格式在接触间断的捕捉上具有较明显的优势,而对于激波的捕捉则差别不大.对于三维问题则是WENO格式对流场的分辨更精细.最后对上述结果给出解释,并且提出了可能的改进方法.