在非结构动网格上追踪多介质界面的ALE方法

在处理带有介质分界面的多介质流动问题时,把介质边界当作一种特殊的边界,并且在该边界上的点都有着双重的定义,分别对应边界两侧的不同介质状态,通过利用HLLC方法求解任意拉格朗日-欧拉方程(ALE),以及运用界面跟踪和动网格方法分别对一维和二维多介质流进行了数值模拟,并得到了满意的结果.     

再入体实验模型热化学非平衡绕流流场的数值分析

利用数值求解一个和多个振动温度热化学非平衡Navier-Stokes方程的CFD计算程序,对爆轰风洞球锥试验模型的热化学非平衡绕流流场进行了数值模拟,分析了分子组分振动温度在全流场(头身部流场和底部流场)中的分布规律.研究结果表明:(1)计算的压力、电子数密度以及流场的光辐射数据与试验数据符合较好;(2)在再入体身部,大多数分子的振动温度峰值大大高于平动温度的峰值;在再入体近尾出现振动温度冻结并大大高于平动温度的现象;CO2的两个振动态的振动温度分布非常接近平动温度分布.     

基于雷达散射截面约束的气动优化设计

采用拉格朗日乘子法优化设计了雷达散射截面约束条件下的锥形融合气动外形.拉格朗日乘子法中的极小化问题采用动态演化的优化设计方法求解.该方法是一种基于非定常演化的优化设计方法,即在求解非定常流动支配方程的时候同时履行优化过程,较其它基于定常解的优化方法具有高得多的计算效率.其中的雷达散射截面通过求解非结构的笛卡儿网格上的时域麦克斯韦方程来得到,而升阻比则通过求解锥形流方程来计算.通过优化设计,得到了M∞=8.0时,升阻比为4.98,雷达散射截面只有1.66m2的锥形融合气动外形.     

基于非结构重叠网格的二维N-S方程求解与应用研究

本文发展了非结构重叠网格技术.通过构造任意单元类型网格重叠的网格间边界确定算法、宿主单元搜索算法以及网格间插值算法,把非结构重叠网格推广应用于求解NS方程,并结合湍流模型模拟二维湍流流场.N-S方程的求解采用有限体积法,通量计算采用Osher格式,湍流模型采用SA模型.为验证本文方法的正确性和有效性,模拟了RAE2822翼型和某三段翼型绕流湍流流场,结果表明采用非结构重叠网格模拟复杂外形湍流流场不但网格生成效率高,而且容易推广应用于具有相对运动的非定常问题,具有较高的工程应用价值.     

飞行器外挂投放数值模拟

采用非结构重叠网格Euler方程的CFD方法数值模拟了外挂投放.该方法在重叠网格     

非失速二维振荡叶栅非定常流动数值模拟研究

进一步完善了非定常雷诺平均N-S方程求解程序,借助快速网格生成方法,对周期性俯仰运动下NACA0012叶栅的非失速绕流问题进行了数值模拟研究,并分析了振荡叶栅的迟滞效应以及俯仰运动各个参数对迟滞效应的影响.数值模拟结果与实验结果对比表明,本文发展的数值模拟方法是可靠的;结果还表明,随着振幅,折合频率的增加迟滞效应增强,平均攻角对迟滞效应无影响;保持无量纲参数--振幅和折合频率的乘积--恒定时,迟滞效应相同.     

基于非结构网格的不可压N-S方程多矩有限体积法

提出了一种基于三角形及四面体非结构网格的有限体积法(FVM),用以鲁棒且精确地求解不可压粘性流动问题。与传统的FVM方法仅将体积分平均值(VIA)作为计算变量的做法不同,本文提出的方法将VIA及点值(PV)同时作为计算变量并在每个迭代步进行计算更新。VIA以通量形式进行计算以确保数值守恒,PV可以通过控制方程的不同形式进行求解更新,无需守恒,因此可以采用非常高效的方法进行求解。将PV作为增加的变量使得紧致网格模板得以实现更高阶精度的重构,而且由此获得的数值模型对于非结构网格变得更鲁棒。本文针对二维/三维的三角形/四面体非结构网格提出了数值格式,给出了几个基准测试算例,验证了本文提出的数值方法在采用非结构网格求解不可压粘性流动问题时的精确性和鲁棒性。     

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

结合直角网格计算的高效率以及表面非结构网格对飞行器几何外形的精细描述,本文建立了基于两级直角坐标网格结构的DSMC数值模拟算法.一级网格为计算区域离散后的背景网格,二级网格为独立的碰撞网格和取样网格,是在背景网格的基础上根据流场性质的变化经过自适应后获得的亚网格,碰撞分子对在亚网格内选取,以保证选取的分子是     

高超声速进气道边界层强制转捩试验

最优增长扰动理论是近年来高超声速边界层非模态失稳现象研究中的主要理论之一,可被看作一种对非模态扰动增长的工程化估计。本文以高超声速钝楔边界层为研究对象,采用最优增长理论和二次失稳分析方法,对钝楔边界层的非模态失稳特性进行了研究。结果表明,在平楔面各个工况的有限计算域中无法找到不稳定模态扰动,非模态的最优增长扰动可获得较大的能量增益。研究还发现,最优扰动的能量增益随马赫数的增加逐渐减小,而壁面曲率对最优扰动的增长起促进作用。但是,由于不同曲率条件下对应的最大能量增长扰动的展向波数存在一定差异,因此变曲率的等熵压缩面不存在统一的变化规律。进一步,研究还采用非线性抛物化稳定性方程研究了零频率最优扰动的非线性演化,并以最优扰动形成的条纹边界层的二次稳定性为基本流,采用线性理论进行了二次失稳分析,发现最优增长扰动形成的条纹结构具有较强的二次失稳增长率,这有利于下游边界层转捩成湍流。该研究结果对设计高速飞行器前体进气道的强制转捩装置有一定参考意义。

     

火星进入器高超声速气动力/热研究综述

以火星大气条件下飞行器非平衡流场研究为背景,采用态-态模型对定容静止CO₂系统自30 km火星大气条件（181 K）瞬时升温后的热非平衡过程开展数值模拟。考虑CO₂分子的对称拉伸、弯曲和非对称拉伸三个振动模态共201个振动能级,微观过程包括：振动-平动（VT）能量交换过程引起的同一模态内部的和不同模态间的振动能级跃迁,振动-振动（VV）能量交换过程引起的同一模态内部振动能级跃迁。分析不同振动模态的能级分布演化特点,并考察对应微观过程细节探究内在原因,结果表明：VT过程起支配作用,对由初始低温突然升温后保持定温的气体非平衡过程,VV过程的贡献可以忽略;CO₂分子的弯曲振动模态具有最快的激发速率,平衡时也具有最高的粒子数占比,升温至相似文献