各向异性网格自适应计算在超燃模拟中的应用

采用各向异性网格自适应求解技术,将其应用于DLR超燃冲压发动机燃烧室中的超燃模拟.开展了3个算例,包括采用滑移壁面条件的燃烧室冷流场模拟、采用无滑移壁面条件的冷流场模拟及采用无滑移壁面条件的反应流模拟.模拟中,各向异性网格自适应计算捕捉到了如激波、射流、边界层、火焰面等具有各向异性特征的大梯度区域,并利用各向异性网格进行了很好的加密.对比利用各向同性网格的初始流场计算,各向异性网格自适应计算使基于滑移条件、无滑移条件的冷流计算及反应流计算的网格单元数量分别下降了36.2%,36.4%和36.8%,有效降低了计算规模,而且流场大梯度区域的计算结果更准确,辨析度更好.结果表明:对于像超燃这类具有各向异性特征的问题,各向异性网格系统比各向同性网格系统有更好的计算效率及准确性,同时也表明基于Mach数场构造的各向异性网格系统可以有效应用于超燃计算.      

基于径向基函数与混合背景网格的动态网格变形方法

发展了一种基于径向基函数(RBF)插值和混合背景网格映射的动态网格变形方法。混合背景网格由运动物面附近的各向异性网格和偏离运动物面的各向同性网格组成。物面变形或运动后,包含在各向异性背景网格中的所有物面网格点的位移由指定的物面运动规律得到,而各向同性背景网格点的位移通过径向基函数插值获得。然后变形后的计算网格利用计算网格与混合背景网格之间固定的映射关系代数插值得到。各向同性网格距离物面一定的距离,消除了物面保形对径向基函数插值的要求,从而可以减少径向基函数基点的数量,以提高插值的效率和可靠性。最后,通过二维和三维算例对该方法的变形网格质量、参数影响、网格拓扑依赖性、复杂外形变形能力及网格变形效率进行了测试。测试结果表明,该方法能够高效生成大变形下高质量的变形网格。     

一种适用于三维混合网格的GMRES加速收敛新方法

为提高流场计算收敛效率,发展了一套适用于三维混合网格Naiver-Stokes方程求解的并行广义最小残差（GMRES）隐式时间推进方法。该方法由科学计算可移植扩展工具包（PETSc）中的Krylov子空间求解器实现,线性方程系统中的系数矩阵直接以显式给出以提高算法的稳定性。为进一步提高GMRES方法的收敛速度,对非结构网格的序号进行了重排序,使得系数矩阵的非零元素尽量向主对角线靠近。利用所发展的GMRES方法,完成了对ONERA-M6机翼、AIAA阻力预测会议通用研究模型（CRM）等算例的计算,计算结果与试验结果吻合良好。通过与其他隐式推进方法进行比较,对算法的收敛特性进行了研究。结果表明,所发展的GMRES方法计算更加稳定,残差下降速度相对LU-SGS（Lower-Upper Symmetric Gauss-Seidel）方法更快,尤其是气动力系数向着收敛解逼近的速度更加明显,提高了计算效率。     

基于驻留时间比的近地空间碎片环境建模

为分析近地空间碎片的分布规律,提出了一种以碎片在空间网格内驻留时间为基础的碎片环境统计建模方法.该方法利用多项式拟合和求根方法统计碎片在空间网格内的停留时间,获取模型基础数据,并据此采用多项式预测、插值和时间序列分析等技术,综合分析空间碎片的分布与演化规律.给出了一个基于双行根数(TLE,Two Line Elements)数据的建模实例,该实例通过了ORDEM2000模型的对比验证,并获得了一些更精细的近地空间碎片环境特征.所得建模方法和分析结论可为长期运行的近地航天器轨道设计、碰撞风险评估及防护等提供技术支撑.     

三维曲线坐标系N－S方程通用形式和数字解

运用张量分析理论，比较简捷地导出了三维任意曲线坐标系上椭圆型Ｎ－Ｓ方程的通用形式。在对这些方程进行数值离散时，采用了有限控制体积的非交错网格设置，该算法以ＳＩＭＰＬＥ为基础，为了抑制压力振荡场的出现，在压力修正方程中对压力梯度的离散附加一项１—δｘ的差分。作为算例，本文给出了三维９０°强曲率弯曲管道中气流的不可压缩紊流状态的数值计算结果，并与相关的实验数据进行了对比，结果是满意的。     

格栅翼减阻特性研究

为了探索减少格栅翼阻力的方法和途径 ,进行了超声速M∞ =2 .5 2 1下格栅翼的边框几何形状和尺寸以及格栅翼茎厚度、格栅几何形状对格栅翼阻力特性影响的风洞实验。结果显示 ,格栅翼的边框对格栅翼的阻力影响最大 ,选择合适的边框厚度和剖面形状可以有效地减少格栅翼的阻力     

短舱及离散精度对螺旋桨桨叶载荷分布的影响

采用多块结构网格生成技术,分别生成螺旋桨加中心体和螺旋桨加短舱两种构型的动态面搭接网格。在相同计算状态下,分别对螺旋桨加中心体和螺旋桨加短舱流场进行不同时空离散精度的数值模拟。与螺旋桨加中心体构型的计算结果进行对比,发现短舱的出现较大程度改变了螺旋桨桨叶载荷分布,并且采用不同时空离散精度,载荷分布的改变程度不同。     

武器分离及舱门开启过程数值模拟研究

基于非结构重叠网格技术,发展了一套能较好模拟空腔流动和武器分离问题的数值计算方法.通过模拟典型空腔流动和WPFS模型分离轨迹,并与试验结果比较,说明本文方法可较好模拟空腔复杂流动和武器分离轨迹.采用发展的数值方法对某内埋弹舱舱门开启过程和内埋武器分离轨迹进行了数值模拟,研究了舱门运动的非定常效应对舱门气动特性的影响.研究表明:一定条件下,舱门运动产生的非定常效应对舱门气动特性影响较大,舱门和舱门控制系统设计时必须考虑;在本文给定条件下,弹舱内的内、外侧导弹均能安全分离.通过研究,为内埋弹舱的舱门设计和内埋武器分离安全评估提供了技术手段和依据.     

基于非结构/混合网格的高阶精度DG/FV混合方法研究进展

DG/FV混合方法因其具有紧致、易于推广获得高阶格式及相比同阶精度DG方法计算量、存储量小等优点,自提出以来已成功应用于一维、二维标量方程和Euler/N-S方程的求解。综述了DG/FV混合方法的研究进展,重点介绍了DG/FV混合方法的空间重构算法、针对RANS方程的求解方法、隐式时间离散格式、数值色散耗散及稳定性分析、计算量理论分析,并给出了系列粘性流算例的计算结果,包括用于验证混合方法数值精度的库埃特流,以及方腔流、亚声速剪切层、低速平板湍流、NACA0012翼型湍流绕流等。数值计算结果表明DG/FV混合方法达到了设计的精度阶,且相比同阶DG方法计算量减少约40%,而隐式方法能大幅提高定常流的收敛历程,较显式Runge-Kutta的收敛速度提高1~2个量级。     

多面体网格在静气动弹性计算中的应用

介绍了多面体网格技术在静气动弹性问题中的应用.基于M6算例和径向基函数插值法,探究了多面体网格的CFD计算效率和变形能力,使用AGARD445.6机翼进行了静气动弹性数值计算.结果表明,多面体网格具有计算效率高、变形能力强、计算结果准确等优点,可适用于静气动弹性问题的数值计算.