基于背景网格簇的动网格生成方法

提出了基于背景网格簇实现网格变形的新方法。背景网格簇由任一内边界节点和远场边界角点生成的背景网格构成,并随边界运动而变化。通过保持计算网格节点在每一个背景网格中所在单元的面积坐标不变(三维时保持体积坐标不变)而求出一组期望位置,加权平均该组期望位置,确定出该计算网格节点的新位置。其中权值与内边界节点和计算网格节点距离的倒数相关。新方法引入的背景网格簇可以有效地改进单个背景网格的不足之处。算例结果表明:基于背景网格簇的动网格生成方法实现简单,与弹簧近似法相比,新方法因不需要迭代求解方程组而非常高效,且拥有更强的网格变形能力;与Delaunay图映射法相比,该方法背景网格的单元个数极少,因此易于定位,且不会出现背景网格单元交叉的现象,网格变形能力更强,变形后的网格质量更好。     

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

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

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

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

国家数值风洞工程结构网格变形程序SGDP V1.0开发与应用

对结构网格变形程序SGDP V1.0的理论方法、组成架构、主要功能、运行流程、运行性能、变形算例和工程应用进行了简要的介绍。SGDP V1.0是国家数值风洞NNW项目共性基础库网格变形模块下实现的第1款模块化程序,在吸收前期研究代码的基础上,通过代码重构、接口规范、功能扩展等技术实现开发。目前,SGDP V1.0程序依托RBF-TFI方法实现结构网格变形,采用模块化的代码结构和扁平化的参数设计,可以实现4种不同类型运动下的变形网格生成。SGDP V1.0程序每百万节点网格的读取时间约为0.01s、写出时间约为0.008s、网格更新时间约为0.1s,可在数秒内实现千万量级网格的变形更新,具有非常高的运行效率。网格变形算例和工程应用结果表明,SGDP V1.0程序能够处理复杂拓扑网格的变形问题,且有很好的网格质量维持能力,从而能够获得更加精确的数值模拟数据。     

基于高阶谐波平衡的跨声速颤振高效预测方法

跨声速流场激波及其诱导的附面层分离等非线性因素导致跨声速颤振边界很难被准确预测,尤其是目前工程常用的偶极子格网法,在跨声速时该方法的预测精度大幅下降。在雷诺平均Navier-Stokes方程流场求解器的框架内,利用结构模态建立广义结构运动方程,利用径向基函数建立模态振型的插值方法,结合径向基函数和无限插值两种网格变形方法的优点实现高效高鲁棒性网格变形方法,从而实现颤振时间推进分析流程,利用国际颤振标模AGARD445.6机翼验证程序在跨声速颤振边界预测中的可靠性。然而,时域方法在气动/结构反复迭代,需要耗费大量的计算资源和时间。为了提高颤振预测效率,基于高阶谐波平衡（HOHB）方法快速获得广义力影响系数矩阵,利用该矩阵建立频域模态位移和气动力之间的关系,实现高效颤振频域分析方法。通过二维翼型和三维机翼算例进行验证,结果表明：在不对计算精度产生明显影响的前提下,HOHB方法能够提高颤振预测效率约6倍。     

基于ProCAST仿真文件的精铸模具反变形方法

基于ProCAST数值模拟网格文件,提出了一种精铸模具反变形方法.首先提取所需网格文件中节点信息,将无序节点进行匹配变为有序节点.网格文件简化法与快速排序法的结合提高了排序匹配的效率.建立了模具反变形系统,经验证,所采用的方法取得了良好的效果.     

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

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

基于网格变形技术的涡轮叶片变形传递

为实现学科耦合,需要在涡轮叶片多学科设计优化的气动模型和结构模型之间传递温度、气压和变形等载荷.研究结构变形向气动网格传递的实现方法.参数空间插值方法解决耦合面网格不一致时的变形插值问题;网格变形技术用来调整气动网格内部节点的位置,保证变形传递后的网格质量.某涡轮叶片变形传递实例的表明:该方法能够用于多学科设计优化问题中的变形传递.      

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

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

二维任意域内基于节点的局部网格生成算法

凸域内基于节点的局部网格生成算法,克服了基于节点的有限元方法的网格生成可能产生的不一致性。将该基于节点的局部网格生成算法的适用范围拓展到二维任意域。另外,提出了通过使用约束Delaunay路径来划分任意域的区域划分算法,该算法使得在并行实现网格生成的过程中各处理器之间无需通信,从而大大提高了节点给定情形下有限元方法网格生成的并行效率。     

RBFs-MSA Hybrid Method for Mesh Deformation

Simulating unsteady flow phenomena involving moving boundaries is a challenging task,one key requirement of which is a reliable and fast algorithm to deform the computational mesh.Radial basis functions(RBFs) interpolation is a very simple and robust method to deform the mesh.However,the number of operations and the requirement of memory storage will be increased rapidly as the number of grid nodes increases,which limits the application of RBFs to three-dimensional(3D) moving mesh.Moving submesh approach(MSA) is an efficient method,but its robustness depends on the method used to deform the background mesh.A hybrid method which combines the benefits of MSA and RBFs interpolation,which is called RBFs-MSA,has been presented.This hybrid method is proved to be robust and efficient via several numerical examples.From the aspect of the quality of deforming meshes,this hybrid method is comparable with the RBFs interpolation;from the aspect of computing efficiency,one test case shows that RBFs-MSA is about two orders of magnitude faster than RBFs interpolation.For these benefits of RBFs-MSA,the new method is suitable for unsteady flow simulation which refers to boundaries movement.     

非结构动网格生成的弹簧-插值联合方法

 提出了基于弹簧法的两点改进,以解决边界发生大变形时的非结构网格变形问题。为了使边界运动引起的网格变形能更好地由边界传递到内部网格中,提出了一种新的基于Delaunay网格插值的弹簧倔强系数逐层改进方法。为了提高弹簧法的计算效率,引入背景网格和直接插值方法,提出了弹簧-插值法。弹簧-插值法首先生成计算域的背景网格(粗网格),然后由弹簧法求解边界运动引起的背景网格变形,最后利用变形后的背景网格直接插值得到计算网格的变形。算例结果表明：改进后的方法一方面有效地提高了动网格的变形能力和变形后的网格质量;另一方面通过降低弹簧法的求解规模,显著地提高了动网格的变形效率。     

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

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

An efficient large-scale mesh deformation method based on MPI/OpenMP hybrid parallel radial basis function interpolation

An efficient MPI/OpenMP hybrid parallel Radial Basis Function (RBF) strategy for both continuous and discontinuous large-scale mesh deformation is proposed to reduce the computational cost and memory consumption. Unlike the conventional parallel methods in which all processors use the same surface displacement and implement the same operation, the present method employs different surface points sets and influence radius for each volume point movement, accompanied with efficient geometry searching strategy. The deformed surface points, also called Control Points (CPs), are stored in each processor. The displacement of spatial points is interpolated by using only 20–50 nearest control points, and the local influence radius is set to 5–20 times the maximum displacement of control points. To shorten the searching time for the nearest control point clouds, an Alternating Digital Tree (ADT) algorithm for 3D complex geometry is designed based on an iterative bisection technique. Besides, an MPI/OpenMP hybrid parallel approach is developed to reduce the memory cost in each High-Performance Computing (HPC) node for large-scale applications. Three 3D cases, including the ONERA-M6 wing and a commercial transport airplane standard model with up to 2.5 billion hybrid elements, are used to test the present mesh deformation method. The robustness and high parallel efficiency are demonstrated by a wing deflection case with a maximum bending angle of 45° and more than 80% parallel efficiency with 1024 MPI processors. In addition, the availability for both continuous and discontinuous surface deformation is verified by interpolating the projecting displacement with opposite directions surface points to the spatial points.     

基于径向基函数的网格变形及非线性气动弹性时域仿真研究

为开展非线性气动弹性研究,基于非线性结构有限元软件NASTRAN和自主研制的多块结构化计算流体力学(CFD)求解器,开发了一套基于计算流体力学/计算结构动力学(CFD/CSD)耦合求解方法的气动弹性时域仿真程序.该程序采用径向基函数(RBF)交换两套求解器之间的数据并进行网格变形.为提高RBF方法的效率,构造了基于多次插值的空间待插值点精简算法.在多次插值过程中,每次插值的对象为上次插值的误差,并同时限制插值区域,以此实现了空间待插值网格数的精简.数个网格变形的算例表明该方法可支持大变形运动,并且具有较高的计算效率.采用此程序开展了AGARD 445.6机翼颤振计算、大展弦比机翼的静气动弹性计算与切尖三角翼极限环振荡(LCO)现象的动气动弹性仿真,结果揭示了当机翼展弦比较大或者响应幅值较大时,结构非线性对于气动弹性有显著影响.     

Delaunay graph-based moving mesh method with damping functions

The fluid–structure interaction and aerodynamic shape optimization usually involve the moving or deforming boundaries, thus the dynamic mesh techniques are the key techniques to cope with such deformation. A novel dynamic mesh method was developed based on the Delaunay graph in this paper. According to the Delaunay graph, the mesh points were divided into groups. In each group, a factor ranging from 0 to 1 was calculated based on the area/volume ratio. By introducing a proper function for this factor, this method can control the mesh quality with high efficiency. Several test cases were compared with other dynamic mesh methods regarding mesh quality and CPU time, such as radial basis function method and Delaunay graph mapping method.     

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

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

基于非结构背景网格的动网格方法研究

借鉴弹簧动网格方法和Delaunay背景图方法,提出一种基于非结构背景网格的动网格生成方法。先用少量的物面控制点生成相对较粗的非结构背景网格,再用弹簧法将动边界的变化量作用到该背景网格,然后利用计算网格与背景网格的映射关系插值生成对应的动网格,最后对物面网格点及其附近受影响的密网格点进行局部网格修正,获得最终计算所需的动网格。所提方法与直接弹簧法比较,由于弹簧法作用的背景网格远粗于实际的计算网格,提高了计算效率,此外,由于网格的运动受控于受弹簧法作用的背景网格,边界变形或运动能较为均匀地传递到内部网格上,使网格单元过渡光滑,提高了动网格的质量。给出若干算例,展示了方法的效率和可行性。     

多块协调变形的网格变形技术及其应用

基于多块网格协调变形的理论,提出了适用于叶轮机内流多物边界的更具一般性的网格变形技术,先获得块边界处的变形量,然后插值得到块内网格各点的变形量。通过二维振荡叶栅和带叶尖间隙的三维跨音风扇叶片的非定常绕流算例,证明网格变形技术在叶轮机气弹稳定性数值模拟中的应用是有效的,能够在保证网格高精度的条件下,提高网格变形速度。