非结构网格上欧拉方程的自适应逆风隐式有限元解法

本文发展了一种解二维欧拉方程的隐式逆风有限元格式。空间近似应用的是二阶精度逆风方法；时间近似上根据数值通量函数的线性化处理导出了有限元的隐格式。为减少数值振荡，引入了基于特征变量的限制因子。在经自适应处理的非结构网络上，本文对跨、超、高超声速绕流进行了数值模拟。计算实践表明发展的隐式有限元格式与显格式相比具有较高的计算效率和较强的适应性。     

高超声速热化学非平衡钝体绕流数值模拟

对高超声速热化学非平衡钝体绕流进行研究,发展了高温热化学非平衡条件下的高超声速气动热数值模拟方法。采用二阶精度迎风TVD格式的N-S方程有限体积法多块结构网格求解器,考虑了空气5组分及化学非平衡效应的影响,对组分生成源项及振动能源项采用了点隐式的处理,提高了计算效率。对球锥模型和RAM-CII飞行实验等算例进行了数值气动热模拟,分析了不同非平衡模型在热化学非平衡条件下流场的影响。计算结果表明建立的数值模拟方法具有较高的精度。     

高超声速压缩拐角湍流计算的收敛性分析

应用基于shear stress transport(SST)湍流模型的隐式紧耦合方法对高超声速压缩拐角流动进行了数值计算.考察了差分格式、限制器以及湍流时间步长等7个因素对湍流计算过程稳定性和收敛性的影响.结果表明:差分格式和限制器对收敛性有很大影响,总体而言耗散性越大的差分格式和限制器收敛性越好;增加湍流时间步长对收敛性影响很小;适当降低湍流变量最低限制值有利于湍流方程的收敛;增大Courant-Friedrichs-Lewy(CFL)数和隐式内迭代次数可显著加速收敛,但均存在最优值,超过最优值则加速效果不明显.      

近空间高超声速流场通量分裂型DSMC-IP方法

针对近空间高超声速流场的特点,采用Van Leer格式对直接模拟蒙特卡洛-信息保存(DSMC-IP)方法质量守恒方程中的计算格式进行改进.将局部马赫数作为分裂通量的标准,并重构单元分界面两侧的左右输运通量,使得计算格式具有通量分裂的特点,解决了DSMC-IP方法在高超声速流场计算中的应用问题.将改进后的通量分裂型DSMC-IP方法引入非结构网格中,对二维近空间高超声速流场进行数值模拟.计算结果表明:通量分裂型DSMC-IP方法所得出的数值结果与实验值及参考值符合较好,明显降低了直接模拟蒙特卡洛(DSMC)方法所带来的统计耗散.当来流气体的稀薄程度增加时,其非平衡效应也更加明显,而通量分裂型DSMC-IP方法的计算结果与参考值相差均在10%以内,良好地反映了非平衡条件下的流场特征,验证了通量分裂型DSMC-IP方法的可行性和有效性.      

高超声速化学非平衡粘性激波层方程的数值解

本文讨论了高超声速粘性激波层方程数值计算时差分格式引起的物理失真问题。具体分析了全隐格式格式粘性的影响,并作了数值试验。为了验证隐式结果的可靠性。在超声速激波风洞中测量了钝锥的表面压力分布,并与计算结果作了比较,两者基本一致。 本文采用隐式有限差分法数值计算了高超声速化学非平衡粘性激波层绕细长球锥的流动。计算时采用连续方程和法向动量方程耦合求解的方法以解决细长体远后身区计算中的问题。应用网格技术和加强系数矩阵主对角元素优势的方法提高了化学非平衡流计算的雷诺数范围。文中给出了高超声速化学非平衡流的计算结果,并与其它文献的结果作了比较。     

更准确、更精确、更高效——高超声速流动数值模拟研究进展

从准度、精度和效率3方面回顾了近几十年来高超声速流动数值模拟研究的进展。在物理模型方面,介绍了高超声速数值模拟中高温气体效应、稀薄气体效应以及湍流效应的建模与模拟,基于雷诺平均Navier-Stokes(RANS)方程重点对现阶段较为关注的高超声速边界层转捩的模式理论研究进行了介绍。在空间离散算法方面,主要介绍了高超声速数值模拟中常用的二阶精度迎风格式以及高阶精度格式的发展及其应用。在时间推进方面,主要回顾了隐式时间推进方法的发展及其应用。在误差和不确定度估计方面,主要介绍了其概念、来源以及常用的分析方法,同时给出了迭代误差估计、Richardson外插法以及敏感性导数方法等初步研究结果。最后,讨论了高超声速流动数值模拟中下一步需关注的问题。     

高超声速飞行器多层次结构强度分析方法

针对高超声速飞行器结构应力损伤问题,采用有限元子模型法对高超声速飞行器结构强度进行分析研究。有限元分析依据模型复杂情况分为3个层次进行:第1层次采用较为稀疏的网格,获得二级子模型边界上各点的位移和力分布;第2层次将二级子模型网格划分细密,并用Fastener单元模拟群钉连接结构,获得较为准确的应力计算结果及钉载分配;第3层次选取钉载最大的部位建立三级子模型,引入渐进损伤子程序对危险部位进行损伤分析。文中采用ABAQUS子模型法,结合Fastener单元及UMAT(User-defined Material Mechanical Behavior)对高超声速飞行器结构强度由整体到局部进行了有限元分析,解决了常规有限元分析法网格多、计算困难等难题。     

非结构网格高超声速化学非平衡MHD流场数值模拟

 在非结构网格上对考虑化学非平衡效应的二维高超声速磁流体绕钝头体流动进行了数值模拟。控制方程由二维理想磁流体动力学（MHD）方程和组元连续方程两部分组成,化学动力学模型为5组元17反应模型。MHD方程空间离散采用AUSM格式,时间推进采用显式5步龙格-库塔格式,并通过弱耦合的方式与化学反应控制方程结合在一起。计算模型为二维钝头体,外加磁场为偶极子场,磁场源位于钝头体内部。在高超声速来流条件下,对有、无磁场干扰,是否考虑化学反应下的4种工况进行了数值计算,得到了满意的结果,并与有限的参考文献进行了对比。结果表明本文发展的方法能准确地模拟考虑化学非平衡效应的高超声速MHD流场。     

基于人工粘性的间断Galerkin有限元方法

发展了一种基于人工粘性的间断Galerkin有限元方法,作出改进以增强适应性,并向非均匀网格推广.选取了典型算例对方法进行验证.一维激波管算例表明,改进的方法在保持计算结果高分辨率的同时,能够更好地抑制非物理振荡.分析得知,当人工粘性方法用于二阶DG格式时,所得计算结果的数值耗散较大,而当格式精度大于二阶时,采用人工粘性方法所得的结果的分辨率较高.通过计算圆柱高超声速粘性绕流,将三阶DG格式与三阶MUSCL格式和五阶WENO格式的结果进行对比,结果表明,该人工粘性方法对于高超声速流动计算也具有一定的优势.     

二维带动力吸气式高超声速飞行器绕流的PNS-NS混合求解

为了提高吸气式高超声速飞行器绕流的求解效率,采用空间推进求解抛物化PNS(parabolizedNavier-Stokes)方程和时间迭代求解Navier-Stokes(N-S)方程的混合计算流体动力学(CFD)方法来求解高超声速飞行器整机绕流.在超声速占主导的流动区域采用空间推进求解抛物化N-S方程的方法,在亚声速和分离区采用时间迭代求解N-S方程的方法.对于求解二维带化学反应的吸气式高超声速飞行器绕流,混合CFD方法和完全时间迭代方法相比,可得到同等准确的数值模拟结果,并且求解效率提高了数倍.      

一种基于三维粘性流场的乘波体生成方法

为了便于乘波体的优化设计工作,本文实现了一种基于任意三维粘性流场,在激波面上截取乘波体前缘线来生成乘波体外型的设计方法,并与NASA Langley研究中心的CFL3D软件的算例结果进行了比较检验,误差满足工程计算要求.由于在流场计算中直接计入了粘性,而不像大多数的研究方法那样通过无粘计算再进行粘性修正,而且对于高超声速乘波体设计,可以在程序中考虑真实气体效应、化学反应等因素,因此本文的设计方法更接近工程实际,更具有推广价值.     

计算机翼跨音速绕流的加罚Galerkin有限元法

1.引言 全速位Galerkin有限元法是跨音速有限元法中重要方法之一。处理超临界流时,为了自动捕捉到激波,通常采用上风技术,但因此导致了质量守恒条件得不到完全满足。文献[1]从Bateman变分原理出发,应用加罚函数方法,找出了强制满足质量守恒条件在Galerkin有限元方程中所应加入的附加项。本文将文献[1]的方法应用到超临界跨音速绕机翼流动,克服了一般Galerkin有限元法中由于采用上风技术而出现的激波不够陡和位     

粘性不可压流动压力修正的有限元计算

将ＳＩＭＰＬＥ方法压力修正思想应用于有限元方法，对粘性不可压缩流动的ＳＵＰＧ有限元方程，建立将压力与速度分离，进行压力修正迭代的计算格式；对典型算例进行多个Ｒｅ数的层流，湍流计算，显示此方法具有良好的收敛性和精度。     

高超声速湍流直接数值模拟技术

概述了近年来国内外高超声速湍流直接数值模拟(DNS)技术方面最新的研究进展,主要集中在高精度、高鲁棒性数值方法方面,同时也介绍了近年来典型的高超声速湍流DNS算例。在数值方法方面,主要介绍了高精度激波捕捉格式以及保持计算稳定的数值技术,重点是WENO格式及高阶保单调格式的最新进展。在高超声速湍流DNS算例方面,介绍了压缩性影响、壁温影响、真实气体效应以及高超声速转捩等方面的DNS研究。此外,还简要介绍了作者开发的可压缩高精度计算流体力学软件OpenCFD。     

高超声速稀薄流的气粒多相流动DSMC算法建模研究

基于直接模拟Monte Carlo(DSMC)方法,构造适用于DSMC算法的固态和液态颗粒碰撞、聚合和分离模型,发展稀薄条件下双向耦合作用的气粒多相流的DSMC算法,在此基础上初步实现高超声速稀薄流环境中的气粒多相喷流流场数值模拟.算例结果表明该方法能为稀薄过渡区气粒多相流动提供一种新的应用研究手段.     

高分辨率有限差分—有限元混合方法及其在气动热计算中的应用

在文献［１］结合ＮＮＤ格式思想＾［２］提出的有限元格式基础上，给出一种有限差分－有限元混合方法。通过求解完全Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，得到了高超声速情况下钝头体轴对称绕流的较满意结果，同时讨论了网格Ｒｅ数对驻点热流的影响。     

复合材料叠层板壳材料非线性有限元分析

 本文采用通用于复合材料叠层厚、薄板壳的Heferosis三维退化曲壳元,利用Hill屈服准则和层状模型法,研制了针对金属基复合材料叠层板壳进行弹塑性分析的超参元程序。算例结果表明,叠层板壳分析方法和所编程序是可靠的,精度良好。     

高超声速复杂外形热流的NND有限元数值计算方法研究

有限元法是一种加权积分的方法,这个特点和权函数的特性使得有限元法在计算诸如热流(与温度的一阶导数相关)等物理量时有其内在的优越性.本文研究了高超声速复杂外形热流的NND有限元数值计算方法,并对高超声速球钝锥外形和带翼飞行器外形的绕流进行了数值模拟,算例结果令人满意.