基于非结构重叠网格的二维外挂物投放模拟

本文发展了基于非结构重叠网格的动态网格技术,并提出了新的相邻单元搜索算法。非结构重叠网格集非结构网格和重叠网格的优点于一身,适合处理复杂外形和多体相对运动问题。新的相邻单元搜索算法在碰到物面时能继续搜索下去,搜索效率大大提高。非定常Euler方程的求解,采用有限体积法,时间上采用双时间步长推进。为验证算法的有效性,先对绕振荡NACA0012翼形的流动进行了模拟;然后通过耦合求解刚体三自由度运动方程,模拟了二维外挂物投放过程,结果表明非结构重叠网格在处理多体相对问题时具有优越性。     

二维非定常欧拉方程的有限元解法

本文通过求解时间相关的欧拉方程，对振荡翼型的非定常二维无粘可压缩绕流作了数值模拟。求解方程时采用了有限体积积分的方法，并用ＶａｎＬｅｅｒ逆风通量分裂格式来计算控制体通量。隐式格式的采用，大大提高了计算效率；同时，通过插值使得空间离散达到了二阶精度。另外，本文采用基于非结构三角形网格的动态网格来处理刚性物体的移动，从而顺利地解决了复杂几何外形物体的非定常绕流的问题。文章最后给出了ＮＡＣＡ００１２翼型绕四分之一弦线作简谐俯仰振动的计算结果。一个周期内的瞬时压强分布以及升力和力矩系数都与实验结果吻合良好，     

非结构网格上N—S方程解法研究及在多段翼型绕流中的应用

在非结构网格上求解N－S方程时，靠近物面的网格伸展比相当大，用有限体积法进行显式求解时不易收敛并且雷诺数有一定的限制。本文采用的人工耗散项中含有网络的方向尺度，因而在不同区域选取不同：在粘性区内该项比物理粘性项小得多，在粘性区外与Jameson的人工耗散项相同。生成非结构网格采用的是层推进技术与传统的阵面推进法相结合的方法，这样能保证在粘中性区内网格有一定的方向规则，而在粘性区外网格完全非结构化。     

用DES方法数值模拟喷流开启/关闭时的流场

针对喷流在开启/关闭时的气动特性进行了数值模拟。控制方程为Navier-Stokes方程,其求解采用了基于面的有限体积方法,该离散方法适合任意网格单元。空间离散采用了Osher格式。同时也采用了基于Spalart-Allmaras一方程模型的DES模型。对非定常过程的处理采用了双时间步长。在加速收敛的过程中采用了区域分裂的并行技术,利用METIS进行网格分区,MPI库函数作为并行环境,并采用了非阻断通信来进行区域间的信息交换。对导弹在来流马赫数为6.15时开启/关闭喷流时的流场变化进行了数值模拟,并与实验结果进行了比较。结果显示,该方法对处理此类问题非常有效。     

混合网格在求解三维N—S方程中的应用

用三棱柱／四面体组成的混合网格求解三维Navier-Stokes方程，在物面附近采用三陵柱网格，而其他地方则采用四面体网格，与单一类型网格相比，混合网格易于调节粘性层内网点数目，并能减少存储量，提高解的效率，Navier-Stokes方程的求解采用有限体积法，对人工耗散项进行了改进，使得求解的结果更接近真实的物理解，本文对上述方法开展了数值实验研究。结果表明，用三棱柱／四面体组成的混合网格求解N－S方程是非常有效的。     

二维非结构网格DSMC方法的实现及其应用

研究了非结构网格DSMC方法实现的过程，将Bird位置元方案中的子网格思想引入到非结构网格上来，只存储子网格的总体标识号，做到了利用较少的计算网格提高分子的分辨率和计算精度；提出了将面积元坐标搜索算法与交替数字二叉树搜索算法(ADT)相结合来跟踪模拟分子在网格之间的迁移；使用ADT方法判别分子与物面是否相交，节约了计算时间；最后对过渡流域高超声速绕流进行了数值试验．结果表明这种方法的可行性、高效性及算法的通用性。     

基于基因算法与博弈论的气动高升力优化(英文)

讨论了基于基因算法与博弈论的组合优化算法在高升力气动优化问题中的应用。引入基因算法与博弈论相结合的分布式组合算法 ,可以将复杂的优化问题分解为几个简单的局部优化问题。文中论述了组合优化算法的构造方法 ,并应用于高升力多段翼型气动优化。与传统基因算法的数值计算结果进行了比较 ,表明本文构造的方法具有高效收敛性及强的鲁棒性 ,可广泛应用于先进气动设计问题。     

三维非结构混合网格区域分解及数据处理

讨论了三维非结构混合网格区域分解技术及子分区显示技术。在三维复杂外形的流场计算中,为了得到较高的计算效率及保证数值模拟的准确性,经常要采用混合网格及并行计算技术。本文针对三维非结构混合网格的自动分区及子网格的数据处理进行了细致的讨论,编写了相应的软件。对三维实际外形非结构混合网格的处理结果表明,本文发展的方法及软件系统具有通用性强、人为干涉参数少、自动化程度高等特点,可作为通用软件推广使用。     

非结构网格高超声速绕流数值模拟及涵道构型减阻特性分析

在非结构网格上对二维高超声速化学非平衡粘性绕流进行了数值模拟，并应用此方法对涵道构型高超声速减阻特性进行了数值分析。空间离散采用VanLeer逆风通量分裂格式，时间推进采用显式的Runge—Kutta格式。化学非平衡动力学模型为五组元五反应模型，对化学反应源项进行了点隐式格式处理，温度场的计算采用牛顿迭代法。分别对二维类弹头体涵道构型、二维半圆柱模型的高超声速绕流流场进行了数值模拟，得到了数值结果，并与完全气体绕流计算结果进行了对比分析。结果表明本文方法可以应用于高超声速绕流计算，并且涵道构型具有高超声速减阻特性。     

不同飞行工况下的某高超飞行器流道特征及气动力特性研究

对一种类似于X-43A飞行器的一体化构形进行了全三维数值仿真模拟,将前体下表面侧沿程静压分布的计算结果与实验结果进行了对比,二者吻合较好,表明数值模拟方法正确,结果可信.在此基础上研究了不同飞行工况下飞行器流道特征及气动力特性的变化.研究结果表明:①研究范围内来流马赫数的变化对全流道的流动结构和全机气动力特性有着一定影响,使得飞行器全机的升阻比略有下降;②随着飞行攻角的改变,全流道的流动结构和升力系数变化显著,而阻力系数的变化并不明显;③侧滑角导致了不对称的流动结构,但在本文的研究范围内其进气道以及全机的气动特性的不利影响并不明显.