高超声速流场与结构温度场一体化计算方法

在对国内外流场与结构温度场一体化计算方法的不足进行细致讨论的基础上,提出了一种高超声速流场与结构温度场一体化计算方法.采用统一的积分方程组作为气动加热和结构传热物理过程的控制方程,对整个物理场进行统一的迎风格式有限体积方法离散,给出了流场与结构交界面上温度、温度梯度及导热系数等参数的计算方法.在时间推进方面,定常状态采用多步龙格库塔迭代格式,非定常状态则采用双时间步长方法.采用发展的一体化算法对二维圆管模型的气动加热和结构传热问题进行了数值模拟.结果表明:2s时驻点物面温度为390K,与其他文献的误差在3.1K范围内;稳定时驻点物面温度为647K.      

分级型进化算法及其在翼型跨声速气动优化中的应用

提出了变精度模型的分级型进化算法并对初始翼型为NACA0012的二维翼型进行了跨声速流动下的形状增升优化。借鉴自然进化中不同进化阶段个体适应度评估环境不同的机理，构造了分级精度模型以加速优化过程。对翼型进行了给定气动条件下的形状增升优化，给出了优化结果，并与传统基因算法及单精度模型算法结果进行了对比。     

低马赫数流动中的预处理Euler方程的收敛特性

对预处理后Euler方程在低马赫数流动数值模拟中的收敛特性进行了细致研究。首先采用量级分析法研究了Weiss-Smith预处理方法对低速情形下Euler方程收敛性的影响;然后针对低速情形下预处理后的Euler方程组中各方程的收敛性,提出了一个改进的预处理矩阵。数值模拟结果表明：Weiss-Smith预处理方法可以很好地改善Euler方程在低速流动时的收敛特性,并且改进后的预处理矩阵在连续方程的收敛特性方面效果明显。     

三维非结构混合网格高超声速流场并行计算

在非结构混合网格上对三维高超声速化学非平衡粘性绕流进行了基于PC-Cluster的分布式并行数值模拟.本采用区域分裂思想,研究了三维非结构混合网格区域自动分解技术,并以此为基础对高超声速化学非平衡绕流进行了并行数值计算.控制方程为多组分N-S方程,空间离散采用有限体积格心格式,时间推进为显式Runge-Kutta格式.化学非平衡动力学模型为七组元带电离反应模型,对化学反应源项进行了点隐式处理,温度场的计算采用牛顿迭代法.在PC-Cluster上对三维双椭球模型的高超声速绕流流场进行了基于区域分解技术的并行数值模拟,所得数值结果与参考献中的结果作了对比验证。     

基于基因算法与博奕论的翼型跨音速Euler方程气动优化

将基因算法与博奕论中的 Nash平衡相结合 ,构建了多目标基因优化算法—— Nash基因算法 (NashGAs) ,并对初始翼型为 NACA0 0 1 2二维翼型进行给定跨音速流动下的形状增升优化。计算中应用 Bézier曲线对翼型形状进行参数化 ,避免了非流线型的不合理形状产生 ;采用有限元非结构网格 AUSM+通量分裂格式Euler方程数值解进行个体适应度值评估 ;应用动态网格技术调整计算网络 ,节约了 CPU机时 ;最后给出了优化结果     

黏性边界层修正对前体/进气道一体化乘波布局气动性能影响分析

针对基于流线追踪生成的乘波进气道在内流部分存在较大黏性误差的问题,以构建具有消波隔离段的内外流一体化乘波气动布局为目的,发展了一种多级波系锥导或吻切锥乘波体的边界层黏性修正方法.在基于特征线法的多级前体与截短Busemann进气道一体化乘波布局设计方法基础上,耦合了高超声速可压缩轴对称流动的冯?卡门动量方程积分方法,通...     

三维横向流喷与超声速主流干扰流场的数值模拟

本文对三维隅角机翼声速喷流与超声速主流的干扰流扬进行了数值模拟。三维欧拉方程的求解采用非结构网络有限体积伽辽金法（Finite Volume Galerkin Method）。引入了总体结点积分域的概念,简化了从单元矩阵到总体矩阵的汇总过程。通量的分裂采用Osher格式,通过外差使其由一阶精度上升为二阶精度。发展了一种基于线化流量的逆风非结构网格隐式有限元格式以提高求解精度及效率。最后给出了三维隅角机翼流场的算例。     

高超声速化学非平衡绕流分布式并行算法

讨论了非结构网格高超声速化学非平衡绕流流场的分布式并行数值模拟方法。控制方程采用带化学反应源项的N-S方程,数值离散格式采用Jameson有限体积法。化学反应模型为五组元五反应模型,没有考虑电离效应,对化学反应源项进行了点隐式处理,温度场的计算采用牛顿迭代法。通过调用METIS库函数将整体计算区域进行分裂,使用PVM并行机制在PCs-Cluster机群上对二维钝头体高超声速化学非平衡绕流进行了分布式并行计算,得到了结果,并对算法的加速特性进行了对比分析。     

一种新的用于MHD方程逆风格式的Jacobian矩阵分裂方法

针对理想MHD方程,提出了一种新的基于MacCormack算法的雅可比矩阵分裂方法,克服了原有方法稳定性差的问题,并成功地应用于理想MHD方程的求解.控制方程在非结构混合网格上进行空间离散,其中对流项采用本文发展的逆风向量分裂格式,并引入了双曲型磁场散度清除技术,时间推进为显式5步龙格-库塔方法.对MHD激波管流动和带均匀磁场干扰的二维高超声速钝头体绕流流场进行了数值模拟,得到了与参考文献相吻合的数值结果,表明本文发展的数值分裂方法可以有效地捕捉MHD流场的流动特征,并且具有比MacCormack方法更高的稳定性和计算精度.     

电离对高超声速稀薄流飞行器气动热影响

将电离反应模型扩展到（direct simulation Monte Carlo,DSMC）方法中,研究了电离反应效应对高超声速稀薄流飞行器气动热的影响特性.针对稀薄流场中电子出现带来的实际困难,引入“捆绑法”思想处理电子在流场中的运动,并给出了电离反应模型及电离反应处理方法.在以RAM-C Ⅱ飞行器外形为例对增加了电离反应的DSMC代码进行验证的基础上,以“星尘号”探测器外形为研究对象,针对不同飞行高度下5组元混合气体模型（无电离）和11组元混合气体模型（含电离）的化学非平衡流动开展了数值模拟,细致分析和对比了电离反应效应对探测器气动热的影响规律.研究结果表明:采用的电离反应处理方案能够模拟带电离反应的高超声速化学非平衡稀薄流动.在飞行高度为60km时电离反应对探测器气动热的影响最为强烈,使探测器的驻点热流密度降低了5.12%,电离反应对探测器气动热的影响随气体稀薄程度增加而减弱.