基于虚网格的格心ADT搜索法

针对重叠网格寻点问题,提出了一种直接在格心网格下操作,基于虚网格的ADT(Alternating Digital Tree)搜索方法.通过对计算网格进行扩展,建立虚边界网格,解决了格心网格因覆盖区域不完整而无法直接建立ADT数据结构的问题.搜索结果即为包含解的格心单元集合,可直接对结果列表遍历以得到合理贡献单元,故完全摒弃了可靠性差的Stencil Walk方法.由虚网格的定义,使寻点在边界附近的处理更为灵活,可以准确给出边界附近贡献单元的有效信息,同时简化了虚网格体系的构建.扩展了重叠边界类型,构建的搜索空间能完整覆盖网格范围,解决了对称面重叠问题.算例研究表明:该方法可靠性好,边界处理能力强,有效提高了重叠网格方法对大型复杂网格的解算能力.     

新型重叠网格洞面优化方法及其应用

高度自动化的洞面优化方法是重叠网格关键技术之一。通过对割补法的分析,针对其洞面切割后洞边界位置不确定性及不可预知性问题,提出将其与物面距离优化准则相结合的新型混合洞面优化方法。该方法针对各物体间体网格,采用物面距离优化准则,将洞边界位置控制在物体中间的位置,使重叠区域更合理、可预测;针对背景网格,将其物面距离设置为较大值,采用物面距离优化准则将其落入物体网格内部的单元一并挖去;针对物面处重叠的网格,采用割补法对其进行切割和填补,保证物面及近壁面网格得以有效的重叠。该方法在保证高适用性的前提下提高了网格重叠质量,自动化程度高,无需人工干预。3个典型复杂流动算例计算结果与实验结果吻合良好,网格重叠区域流场变量传递正确,等值线过渡光滑,流场刻画准确,证明该方法准确可靠。     

高精度多维限制器的性能分析

目前常用的限制器大都是基于一维构造,无法在多维情况下保证物理量的单调特性进而导致非物理振荡.为弥补传统方法的这一构造缺陷,多维限制器(MLP)通过多维修正使单元通量值介于周围相邻单元通量的最大值和最小值之间,在保证求解精度的情况下有效避免了多维振荡.基于一维激波管、无黏涡及激波边界层干扰等算例,对高精度MLP的特性进行了研究分析.结果显示:3阶MLP在连续和间断区域均可有效地避免多维振荡;与高阶WENO(Weighted Essentially Non-Oscillatory)方法相比,3阶MLP不仅算法简单、易于实现,还可显著提高求解的精度、保单调性及收敛性.因此可用于工程及科学研究的复杂流动,具有较好的应用前景.      

界面化、自动化的战术导弹CFD系统平台

介绍了研究、开发成功的完全界面化、网格生成和CFD计算自动化的战术导弹CFD系统 平台.该软件平台主要基于Visual Basic和Fortran的混合编程,同时大量运用了数据库集 成技术,大大方便了工程技术部门人员的使用.该平台已经在多个航天部门中得到成功的运 用,并取得了可靠的结果.      

DSMC方法中的自适应当地时间步长法

理论分析表明,Laux提出的DSMC(Direct Simulation of Monte Carlo)方法中的当地时间步长法尽管能够显著缩短流场达到稳定所需的CPU计算时间,提高DSMC程序的运行效率,却存在仿真分子运动和碰撞计算复杂,并需要耗费额外计算机内存的缺憾.对Laux的方法中仿真分子的运动处理时机提出了改进,并改变了其碰撞抽样方法,从而简化了仿真分子的运动和碰撞计算处理,避免了额外的计算机内存消耗.应用改进后的自适应时间步长法,对圆柱的稀薄气体绕流进行了采用和未采用改进的自适应当地时间步长法的对比计算.结果表明,改进后的自适应当地时间步长法能明显缩短流场达到稳定所需的计算时间,对流场模拟结果产生的影响却非常小.      

高超声速再入钝头体表面热流计算

基于跟踪流线的轴对称比拟法,采用纯工程算法、Euler数值计算与边界层内工程算法相结合的方法,对高超声速再入钝头体的表面热流进行了计算,并将计算结果与风洞实验数据进行了对比,两者吻合较好,验证了两种工程算法在计算高超声速飞行器热环境方面的正确性.将两种工程算法与数值求解N-S方程进行对比,表明工程算法在迎风面的热流计算方面有较高的精度,节约了计算时间,很好地满足高超声速飞行器概念研究和设计的需要.      

DES方法对返回舱气动热环境数值模拟

脱体涡模拟方法(DES,Detached Eddy Simulation)由于结合了湍流模型(RNAS,Reynolds Averaged Navier-Stokes)和大涡模拟(LES,Large Eddy Simulation)两者各自优势,是模拟非定常大分离流动的有效方法.采用基于SA(Spalart-Allmaras)湍流模型的DES方法,对高超声速返回舱外形进行数值模拟,计算所得的返回舱表面压力及热流分布与实验结果吻合,证实了DES方法的优越性.最后,在DES方法中加入可压缩修正,结果证实在高超声速流动中,可压缩修正方法在一定程度上能够提高原始模型的预测性能.     

喷流干扰气动热数值模拟的若干影响因素

工程上对喷流气动热进行近似计算时,常常对真实的全化学反应流近似处理,这时需要了解喷流气体和喷管外形的近似方法对计算结果的影响规律。采用数值试验对比分析的方法,探讨了近似计算的2个影响因素:其一是使用量热完全的空气喷流近似燃气喷流进行数值模拟时,喷管几何形状对热流分布的影响;其二是使用空气喷流近似燃气喷流进行数值模拟的方法,分析不同的喷流热力学参数匹配方法对喷流形态及热流分布的影响。通过对不同方法进行对比计算,揭示了各个方法对喷流干扰气动热计算的影响规律,可以用于指导工程分析。      

超音速底部流动的大涡模拟研究

超音速底部流动是CFD的难点问题之一。采用基于可压缩流动的大涡模拟方法求解N-S方程,对超音速底部流动进行了模拟研究,分析了其流动结构。计算结果表明：得到的速度型、底面压力分布与实验值吻合,优于雷诺平均方法的结果。     

圆弧翼型跨声速流动的动态模态分析

跨声速翼型的激波周期性自激振荡会给机翼结构带来附加的脉动载荷,从而加剧飞行器表面结构的疲劳损伤。使用动态模态分解（DMD）方法研究了跨声速下绕厚度18%的对称双圆弧翼型的压力脉动场,分析了DMD提取的各阶主模态的频率特征、压力脉动的空间分布以及压力脉动随激波振荡的时间演化过程,并使用DMD模态进行流场重构。结果表明,DMD方法能准确捕捉流场各特征频率的模态,第1阶模态是激波抖振的主频,在激波的自激振荡过程中占主导作用,前7阶模态重构的流场损失函数降低至4%以内,误差主要分布于激波间断处。