距离倒数加权插值法最优参数探究.pdfVIP

第十五届全国计算流体力学会议论文集 距离倒数加权插值法最优参数研究 周璇，李水乡 (北京大学，工学院力学与空天技术系，北京，100871) 摘要：距离倒数加权法(InverseDistance weiglltingIDW)是目前应用较广的插值方法之一，是高效、简单、 直接的插值方法。本文基于IDW插值法模型和动网格应用实例，研究了基于IDW的网格交形方法中的最优 参数问题a研究了不同网格类型、变形结构边界、运动形式下网格变形能力和网格变形质量与参数选取的关 系。本文以二维基本的平动和转动变形为例，根据实际网格变形效果，得到了参数选取与网格类型和结构边 界无关的一般性结论，以及不同变形形式下参数的取值范围。不论是转动变形还是平动变形，平动项参数取 值范围为[1，2]时，网格变形能力和网格质量均较优。而转动项参数的取值与变形形式相关，转动变形中转动 项参数取值为2；而平动变形时可忽略转动项，即转动项参数趋于无穷大。 关键词：距离倒数插值法、网格变形、动网格、最优参数、插值法 变形能力。类似于弹簧法的弹性体法【8】把网格变形 0引言 当成连续介质力学问题来求解，通过求网格变形能 的最小值来确定网格节点位移。弹性体法的网格变 在工程实际问题的数值模拟中，常常需要求解 形能力和网格质量都较高，但实现复杂、计算量大。 边界运动等非定常流动问题。如自由表面问题、多 Stem等【9】同时采用了多种改进方法：网格变形由 体分离问题、强迫振动问题及流固耦合问题等。在 单元大小和形状的改变决定；使网格形状对弹性刚 求解这类非定常问题时需要改变流场计算网格以 度的影响大于网格单元大小对弹性刚度的影响；引 适应边界的运动。动网格方法是解决这类问题的难 点和关键，尤其对于多体作相对运动或边界做不规 结合弹性体法和局部网格质量改进，不仅增加小单 则变形运动的非定常问题[1】。非结构网格能够方便 元的弹性模量，而且单独处理网格质量较差单元。 地生成复杂外形的流场网格，可通过流场中大梯度 但这些改进在提高网格变形能力时都不可避免地 区域自适应来提高对间断(如激波等)的分辨率，而 增加了计算量和算法复杂性。 且网格分区以及并行计算比采用结构网格更加直 近年来发展较快的网格变形方法是代数法，代 接。主要的非结构网格变形方法可分为虚拟结构法、 数法仅根据网格节点的空闯位置来移动节点，不考 偏微分方程法和代数法[2]。 虑节点邻接关系和网格类型，是一种高效和直接的 最直接的网格变形方法是弹簧近似法(弹簧 网格变形方法。主要的代数法可分为基于Delaunay 法)[3】，把整个计算网格认为是一张达到平衡的弹 背景网格的直接插值法[111、径向基函数插值法B2】 簧网。通过求解椭圆型平衡方程，将边界变形所形 和距离函数插值法[13】等。Liu等提出的基于 成的扰动传递到整个网格中。弹簧法简单易于实现， 但需要迭代计算求解方程组，而且变形量有限。不 时间步生成Delaunay网格作为背景网格，边界运动 少学者通过增加扭转弹簧【4】、虚拟弹簧f5】或结合边 对网格的扰动由Delaunay网格的变形实现。通过 界改进[6】等方法，约束节点位移，提高弹簧法的网 Delaunay网格和计算网格间的映射关系，得到最后 格变形能力。周璇等[7】结合改进后的弹簧法和直接 的动网格。直接插值法不需要迭代计算，因此效率 插值法，不仅提高了网格变形效率，也提高了网格 很高。但Delaunay网格具有局部性，不能将边界扰 539 周璇，李水乡：距离