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

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

基于动态非结构重叠网格法的直升机前飞非定常流场数值模拟研究

 发展了非结构重叠网格技术，提出适合于直升机旋翼运动的动态非结构重叠网格算法，该算法使得在采用非结构重叠网格模拟直升机流场时实现程序自动处理网格，不需任何人工干预。流动控制方程采用非定常可压缩Euler方程，空间离散采用Jameson有限体积法，时间推进采用双时间步长法，为加速收敛，在双时间步长内循环迭代中采用预处理技术。为验证本发展的预处理解算器以及动态非结构重叠网格的算法的可信性和有效性，开展一系列的数值实验，模拟了不同速度范围的NACA0012翼型绕流,Caradonna旋翼悬停和前飞流场以及Robin直升机机身-旋翼组合的前飞非定常流场，结果表明采用本文算法得到的结果是可信的，网格处理是自动高效的，可以应用于直升机全机流场模拟，为直升机设计提供可信指导数据。     

基于非结构网格的非定常流数值模拟方法研究

本文给出了两种基于非结构网格的非定常流动问题的数值模拟方法.一种是基于非结构网格的重叠网格方法,该方法结合了非结构网格方法和重叠网格法的优点,适合于处理复杂外形飞行器的非定常运动问题.另一种是基于Denauley图映射的动弹网格方法,该方法具有较强处理粘性网格的变形运动的能力.本文通过多个数值算例(三维外挂物投放、直升机机身+旋翼前飞、后缘襟翼大偏角运动以及机翼弹性变形等)对这些方法的正确性和有效性进行了验证,结果表明这两类非结构网格的非定常流数值模拟方法是精确高效的.     

基于非结构重叠网格的二维N—S方程求解与应用研究

本文发展了非结构重叠网格技术.通过构造任意单元类型网格重叠的网格间边界确定算法、宿主单元搜索算法以及网格间插值算法,把非结构重叠网格推广应用于求解NS方程,并结合湍流模型模拟二维湍流流场.N-S方程的求解采用有限体积法,通量计算采用Osher格式,湍流模型采用SA模型.为验证本文方法的正确性和有效性,模拟了RAE2822翼型和某三段翼型绕流湍流流场,结果表明采用非结构重叠网格模拟复杂外形湍流流场不但网格生成效率高,而且容易推广应用于具有相对运动的非定常问题,具有较高的工程应用价值.     

用动态非结构重叠网格法模拟三维多体相对运动绕流

 发展了动态非结构重叠网格技术。本文提出了基于阵面推进的相邻单元搜索算法和适合于高阶空间离散格式的网格间插值关系定义方法，前者提高了非结构重叠网格在模拟多体间具有大幅相对运动的非定常问题时网格处理的效率，而后者保证了网格间边界附近流场求解精度；本文非定常Euler方程的求解采用有限体积法，通量计算采用Osher格式，时间积分采用双时间步长推进。为验证本文提出的方法，分别对定常M6机翼绕流和机翼下外挂物投放过程的非定常绕流进行了模拟,结果表明本文提出的方法提高了非结构重叠网格在模拟多体间具有大幅相对运动问题时的效率，保证了网格边界处的精度，具有较高的工程应用价值。     

基于非结构动态网格的翼型积冰过程数值模拟

讨论了积冰气象条件下翼型积冰过程的非结构动态网格数值模拟方法。外部无黏绕流流场采用基于非结构网格的分布式并行计算技术,空间离散采用VanLeer迎风有限体积格式,时间推进为5步显式龙格 库塔方法。采用四阶龙格 库塔方法求解过冷水滴运动轨迹方程,使用了一种合理、简单的物面水积冰量计算方法,提高了计算效率。使用非结构动态网格技术对积冰后的翼型网格进行快速调整。在对NACA0012翼型积冰形状进行校核验证后,对某型飞机在积冰气象条件下飞行时翼型积冰过程进行了数值模拟,得到了与参考文献结果相符的冰形特性。     

基于非结构动态网格的翼型积冰过程数值模拟

讨论了积冰气象条件下翼型积冰过程的非结构动态网格数值模拟方法.外部无黏绕流流场采用基于非结构网格的分布式并行计算技术,空间离散采用VanLeer迎风有限体积格式,时间推进为5步显式龙格-库塔方法.采用四阶龙格-库塔方法求解过冷水滴运动轨迹方程,使用了一种合理、简单的物面水积冰量计算方法,提高了计算效率.使用非结构动态网格技术对积冰后的翼型网格进行快速调整.在对NACA0012翼型积冰形状进行校核验证后,对某型飞机在积冰气象条件下飞行时翼型积冰过程进行了数值模拟,得到了与参考文献结果相符的冰形特性.     

带控制律导弹投放数值模拟

基于非结构重叠网格技术,耦合舵面控制律模块,采用非定常雷诺平均NS方程,建立了模拟控制律作用下导弹投放分离动态特性的CFD数值方法.非结构重叠网格生成技术以“物面相交”方法确定洞边界,在多套网格重叠部分采用算术平均进行物理量的计算传递.采用机翼/挂架/带舵外挂物模型,以及不同的舵面控制规律展开了对比计算,分析了结果.并比较了某导弹有无控制律作用时从某战斗机分离的动态特性,研究了增加控制律的必要性.     

并行化非结构重叠网格隐式装配技术

重叠网格技术广泛应用于复杂外形及多体分离问题的数值模拟,但是超大规模重叠网格装配的鲁棒性和效率依然是实际应用中的瓶颈。基于格心型的有限体积格式,发展了一种并行化的非结构重叠网格隐式装配技术。该技术采用节点的壁面距离作为插值边界的判断准则,并通过物理边界推进确定出活跃区域。此外,基于网格分区策略实现了隐式装配技术的并行化,为了能够适应超大规模的计算网格,在每一个网格子区内单独建立用于查询的交替数字树（ADT）数据结构,整体网格点在各个子区内进行查询,并通过信息通讯接口（MPI）传递查询结果并进行逻辑判断。该并行化的非结构重叠网格隐式装配技术具有较高的自动化程度,通过多体装配算例以及机翼外挂物分离算例,证明了该方法的实用性。     

旋翼/机身干扰非定常流场数值模拟(英文)

基于非结构动态嵌套网格方法对直升机旋翼/机身干扰非定常流场进行了三维Euler方程数值模拟。嵌套网格系统由包含旋翼的运动区网格和包含机身的静止区网格组成,两区网格均来自于同一套初始的整体非结构网格,人工边界以及初始的插值贡献单元对应关系在网格分区时同时确定。在非定常计算中两区网格始终存在一定的重叠区域,无需重新建立人工边界。高效的贡献单元搜索算法可以根据初始的插值贡献单元对应关系很快更新每一个真实时间步上的贡献单元对应关系。为了准确计入旋翼桨叶的挥舞和变距变形,在运动区网格中引入线性弹簧方法。采用非定常Euler方程和双时间推进方法模拟了美国佐治亚理工学院旋翼/机身干扰实验外形非定常流场,计算结果与实验值以及文献的计算值相吻合。     

柔性扑翼非定常流场的数值计算方法

提出一种将Delaunay图映射网格变形技术和非结构嵌套网格方法结合使用的策略,解决网格变形和嵌套网格单独用于柔性扑翼流场计算时需要网格再生的问题。该方法为嵌套网格中的每个嵌于背景网格的贴体非结构网格生成Delaunay背景图;每个时间步,根据扑翼的运动和变形规律移动背景图,再根据网格点和背景图的映射关系移动网格点,之后自动完成嵌套边界的定义和插值关系的建立。为方便嵌套关系的建立,嵌套网格进行分层管理。也研究了一种内存消耗少、效率较高的搜索算法,以及格心格式和格点格式统一的边界拓宽算法。非定常可压缩Navier-Stokes方程在非结构的动态网格上用有限体积法离散,并用预处理的双时间步推进、隐式LU-SGS迭代求解。几个扑翼算例的结果表明,该方法充分利用了Delaunay图映射网格变形方法的高效率,同时也发挥了嵌套网格处理大幅运动的优势;用于既有整体大幅扑动又有局部小变形的柔性扑翼流场计算,可取得令人满意的精度和效率。     

大幅运动复杂构形扑翼动态网格生成的一种新方法

 基于Delaunay图映射的动态网格生成方法无需迭代计算,效率高,稳定性好。但对大幅运动复杂构形的动边界问题,背景图容易交叉,重新生成背景图和重新定位网格点信息不仅费时而且会导致网格质量的严重下降。提出一种双重Delaunay图映射的动网格生成方法,通过在初始背景图中添加辅助点,生成虚实两种背景图和虚实两种映射关系。分别根据虚映射关系和虚背景图、实映射关系和实背景图,移动辅助点和网格点。几个复杂构形的扑翼算例表明,双重图映射方法多付出极少的内存代价即可避免背景图交叉引起的问题,提高了动网格生成的效率和质量,增强了处理大变形复杂动边界问题的强健性。     

基于径向基函数改进的Delaunay图映射动网格方法

为了提高大变形下动网格生成的效率和质量,基于径向基函数插值技术对原始Delaunay图映射动网格方法进行了改进。首先通过带多项式的径向基函数插值方法计算出背景网格远场边界点的位移;然后将背景网格位移插值到计算网格;最后利用衰减函数将计算网格远场位移衰减为零。基于矩形旋转网格变形实例比较了改进方法与原始Delaunay图映射动网格方法之间的差别,并研究了径向基基点数目和衰减函数类型对变形网格质量的影响。矩形旋转网格变形实例说明该方法能够精确恢复出Delaunay背景网格的旋转特性。NACA 0012翼型、NLR 7301两段翼和M6机翼网格变形算例进一步证明,通过添加控制点,该方法能够不重构背景网格实现大变形下高质量动网格的生成。     

基于非结构网格的气动弹性数值方法研究

基于非结构网格和网格变形技术,研究了耦合求解流体和弹性结构体互相作用问题的数值方法.气动弹性模型采用了耦合求解可压缩雷诺平均的Navier-Stokes方程和线性结构动力学方程的方法.为了能处理在复杂几何域内的问题,采用了基于混合单元的非结构网格有限体积方法.开发的方法在二维翼型简谐振动进行了验证并应用于AGARD 445.6翼型颤振边界的数值模拟.      

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

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

气动弹性计算中网格变形方法研究进展

网格变形是气动弹性计算中实现计算网格随运动边界变形的主要方法。在总结目前网格变形方法发展现状的基础上,对近几年常用的网格变形方法,即弹簧法、弹性体法、超限插值法、Delaunay背景网格法、径向基函数插值法和温度体法等做了简要的总结。根据各方法构建模型的不同,将它们分成物理模型法、数学插值法和混合方法3类,简要介绍了各方法的基本思想和研究进展,重点比较了各方法的网格变形特性（变形能力、变形质量和变形效率）和优缺点。总结了当前流场边界与结构边界的数据传递方法,对当前气动弹性计算中遇到的网格变形的难点问题作了简要评述并对未来网格变形方法的发展方向进行了探讨。     

基于Delaunay图映射的动弹网格技术及其应用

在多块混合网格的基础上发展了一套基于Delaunay图映射方法的动弹网格技术,实现了针对气动弹性计算的网格变形运动方法,并使用基于多块混合网格的RANS方程的有限体积法耦合结构静平衡方程,针对大展弦比无人机在实际使用中可能出现的显著的静气动弹性问题进行了计算、分析和研究。研究结果表明,所发展的动弹网格方法为粘性混合网格的变形提供了较好的解决方案;RANS方程解算器具有较高的精度和较强的适用性,具备求解粘性流体与结构耦合等复杂航空气动弹性问题的能力。     

两种典型非结构网格变形方法特性对比研究

非定常流动仿真中常常遇到运动边界的情形,基于物理模型的线性弹簧法和基于数学方法的径向基函数(RBFs)插值法是实现计算网格随运动边界变化的两种主要的非结构网格变形方法。以NACA0012翼型的俯仰运动为例,对上述两种方法的变形特性进行较详细地定量分析研究,比较线性弹簧法和RBFs方法的CPU计算时间和最大变形能力,分析两种方法的网格质量随俯仰角度和运动步数的关系,并给出典型状态下网格单元质量云图和流场计算结果。结果表明:基于RBFs插值的网格变形方法计算效率较高,变形能力较强且能保持较高的网格质量,是一种高效的网格变形方法。