基于图剖分的多块结构网格负载平衡方法

负载平衡是影响并行计算性能的重要因素。针对多块结构网格,给出了一种改进的多层次图剖分负载平衡方法。该方法设计了新的网格剖分算法,采用改进的子块分裂方法与图剖分算法的循环调用实现结构对接网格剖分,并通过建立不同物体重叠网格间的连接关系,实现了结构重叠网格的负载平衡。采用2个典型算例对方法进行了对比验证,数值结果表明,子块分裂方法对剖分结果具有重要影响,采用循环调用算法及改进的子块分裂方法能有效地实现计算负载均衡及通信量优化,同时显著减少了网格块数及因虚网格导致的内存需求,有利于提高并行效率。该负载平衡方法与网格拓扑无关,适用于多块结构对接网格及重叠网格,且整体型剖分方式对于多块结构重叠网格具有更好的剖分效果。     

Development process of muzzle flows including a gun-launched missile

Numerical investigations on the launch process of a gun-launched missile from the muzzle of a cannon to the free-flight stage have been performed in this paper. The dynamic overlapped grids approach are applied to dealing with the problems of a moving gun-launched missile. The high-resolution upwind scheme(AUSMPW+) and the detailed reaction kinetics model are adopted to solve the chemical non-equilibrium Euler equations for dynamic grids. The development process and flow field structure of muzzle flows including a gun-launched missile are discussed in detail.This present numerical study confirms that complicated transient phenomena exist in the shortly launching stages when the gun-launched missile moves from the muzzle of a cannon to the freeflight stage. The propellant gas flows, the initial environmental ambient air flows and the moving missile mutually couple and interact. A complete structure of flow field is formed at the launching stages, including the blast wave, base shock, reflected shock, incident shock, shear layer, primary vortex ring and triple point.     

基于超网格的重叠网格守恒插值方法

保证重叠网格边界数据插值的守恒性是计算流体力学面临的一大难题。基于超网格方法与格心格式有限体积法,发展了一种新型混合重叠网格守恒插值方法。详细研究了用网格切割技术在重叠网格边界构造局部超网格,用网格求交算法合理地扩大贡献单元模板,建立了一种适用于任意网格类型的隐式并行重叠网格守恒插值方法,以超网格为媒介可实现重叠网格二阶精度的守恒插值。数值结果表明,本文方法对二阶分布的流场变量具有严格的守恒性,相比三线性插值方法和逆向距离权插值方法,本文方法减小了数值误差,提高了重叠网格边界的插值精度,加快了计算收敛速度,改善了重叠区域网格尺度相差较大时流场的光滑性和连续性。     

高效非结构网格变形与流场插值方法

非结构动网格在含动边界的非定常流动模拟中应用十分广泛。为了提高网格变形的计算效率和对复杂外形的适用性,发展了高效的K近邻-径向基函数（KNN-RBF）动网格方法。并采用高效的交界面节点一致性方法实现了网格变形的并行化。由于边界的运动和变形形式多种多样,在自主开发的计算流体力学（CFD）程序中实现了动网格的用户自定义函数（UDF）功能模块,可以灵活地设置此类问题中的边界运动。在模拟大变形问题时,采用了全局网格重构的策略对网格质量较差的区域进行处理。提出了快速查找新旧网格最近点的两级KD（K-Dimensional）树方法,并根据最近点在新旧网格间进行流场插值。从而建立了高效健壮的并行流场插值方法。通过对若干算例进行模拟,验证了所发展方法的正确性和健壮性。     

动基座近舰面流场数值模拟

航母尾迹流场对舰载机的着舰有较大影响,所以需要对其流场特点进行研究,分析不同状态下舰载机气动特性的变化。采用嵌套网格技术对航母处于垂荡状态下无人机的着舰进行了模拟。首先,利用SFS2舰船进行数值计算,验证了舰船流场的数值模拟方法。然后,对比了单独无人机定常与非定常计算结果,表明所建立的嵌套网格适用于无人机流场的模拟。接着,对航母单相流和两相流的流场结果进行了分析,结果显示甲板下方的流动对甲板上方流场没有大的影响。因此,忽略了水的影响只对航母在空气流场中的特性进行研究,结果表明航母尾迹非定常特性明显,静止航母下滑轨迹上的速度均处于周期性波动状态,且波动幅值随着远离航母而逐渐衰减;而在垂荡情况下航母尾迹变得更加紊乱,水平方向速度波动的周期性减弱,但垂向速度的波动幅值进一步增大。对于静止航母,无人机在不同时刻着舰气动特性的变化也存在差异;当航母处于垂荡状态时,无人机的升力和俯仰力矩在短时间内会有更大的波动。     

基于非结构网格的燃面推进算法

根据惠更斯波传播原理,构建了一种在通用CFD软件的非结构网格系统上直接计算燃面推进的数值方法。该方法将燃面对应的计算网格面作为一系列平面波源,利用网格几何关系对边界节点奇异性进行判定,通过所有平面波源影响区域的叠加求解新的燃面位置。数值实现过程中,根据节点奇异性和波源影响区域以三维矢量运算得到节点的推进矢量。结合网格的重分,在通用CFD软件环境中实现了燃面的推进,所用数值算法具有良好的精度。另外,对三维动态内弹道模拟中燃面计算的稳定性和精度等问题进行了研究。     

洞映射方法的研究和改进

针对传统洞映射方法存储大的缺点,对"广义封闭"的概念进行扩展,提出了最小洞映射方法,该方法允许挖洞曲面结束于网格截断面,有效缩小了洞映射区域,节省了存储开销.对适用于广义封闭的洞映射单元识别方法进行了分析,指出了用物面信息判断映射单元属性可靠性较差.发展了一种新的特别适合广义封闭问题的识别方法:"Inverse mark",使用计算网格结点信息自动识别洞外单元,再作为初始点在网格内部推进.研究表明:"Inverse mark"方法计算效率高,可靠性好,自动化程度高,有效提高了重叠网格方法对缝隙等局部复杂结构的适应性     

混合重叠网格插值方法的改进及应用

混合重叠网格间合理插值是保证流场计算正确的基础之一。本文针对重叠区网格尺寸匹配度较差时插值误差较大问题,发展了一种新型混合重叠网格插值方法,通过使用二阶精度插值、按单元尺度区分和扩充模板,改善了插值区网格尺寸匹配度较差时的插值精度。本文方法适用于任意单元类型的混合网格重叠,对各类单元处理透明,实现简单。计算结果表明,采用本文方法在网格重叠区流场变量传递正确有效,插值区网格尺寸匹配度较差时,相比原始方法,等值线过渡更为光滑,变量经过插值区耗散更小,计算与试验值符合更好。     

网格之间流场信息传递的高精度方法和验证

理论分析了插值方法引入误差的原因,提出一种利用时间推进方法和移动网格技术实现网格之间流场信息传递的思想,并构造计算方法。在介绍新的信息传递算法原理的基础上,根据三维情况下传递过程需要满足的时间推进步长稳定性条件,推导出最小推进步数的要求。最后,通过对激波诱导物体运动的非定常流动模拟,比较了信息传递算法和线性插值、加权最小二乘插值方法,计算结果表明,与插值方法相比,信息传递过程没有引入误差,能有效抑制非物理振荡,保持流场计算精度。     

多体分离数值模拟及方法验证

利用自主研发的混合网格Euler／RANs方程解算软件WoF90以及六自由度动力学与运动学方程解算程序,采用脚本语言技术实现网格的自动生成,建立了一套可用于多体分离问题的准定常Euler流动非结构网格方法。通过AEDC算例对非结构准定常Euler方程方法解决多体分离问题进行了计算研究,验证了方法的可行性。