《一种基于径向基函数的非结构混合网格变形技术》.pdfVIP

2011年 10月 西 北 工 业 大 学 学 报 0ct． 2011 第29卷第5期 JournalofNorthwesternPolytechnicalUniversity Vo1．29 No．5 一 种基于径 向基函数的非结构混合网格变形技术 王 刚，雷博琪，叶正寅 (西北工业大学 翼型、叶栅空气动力学国防科技重点实验室，陕西 西安 710072) 摘 要：网格变形技术被广泛运用于气动外形优化设计或气动弹性力学仿真之 中。非结构混合网格 由于拓扑结构比较复杂，其网格变形的难度较大。本章发展了一套基于径向基函数的非结构混合网 格变形技术。其基本原理是运用径向基函数对物面边界网格节点的位移进行插值，然后利用构造出 来的径向基函数插值序列将物面的位移效应光滑地分散到整个网格区域的节点上。为了提高网格变 形效率，需要在物面位移的径向基函数插值过程中进行数据精简，为此，文中提出了一种基于贪心法 逐级选择径向基函数空间子集来实现插值逼近的数据精简算法。选择NACA0012翼型和三维LANN 机翼的非结构混合网格的典型变形问题作为数值算例，对该方法的实用性进行了验证。计算结果表 明该方法具备较高的计算效率，对大尺度变形问题的适应性很好，变形后计算网格的质量仍然可以得 到有效地保证。 关 键 词：径向基函数，非结构混合网格，网格变形，贪心法 中图分类号：V211．3 文献标识码 ：A 文章编号：1000-2758(2011)05-0783-06 在气动弹性数值仿真、非定常流动数值模拟以 渐进处理。 及气动外形优化设计等方面的研究工作和工程实践 近年来，采用径向基函数插值的网格变形技术 中J’2J，经常需要使用网格变形技术对当前计算网 逐渐兴起，目前已经在流固耦合数值模拟中获得了 格进行适当变换，从而在保留当前网格拓扑结构的 成功的应用 ．8J。其基本原理是运用径向基函数对 前提下，使计算网格能够准确反映物面边界的形状 物面边界网格节点的位移进行插值 (或拟合)，然后 变化。对结构网格而言，其拓扑结构相对比较简单， 利用构造出来的径向基函数序列将物面的位移效应 一 般采用简单的代数插值方法 就可以满足常规 光滑地分散到整个网格区域的节点上。对于物面网 的网格变形需要。而非结构网格的网格变形要比结 格节点数 目较多的三维问题而言，直接采用所有物 构网格复杂许多。特别是对于粘性流动数值模拟中 面节点进行径向基函数插值的计算量很大。Renda． 使用的非结构混合网格，网格变形方法必须处理的 1l等 采用贪心法在物面节点位移的拟合过程中进 拓扑结构更加复杂。目前针对非结构网格提出的网 行数据精简，以便尽量减小径向基函数函数空间的 格变形方法主要分三类 ：弹簧法 j、弹性体方法_5 维数，从而获得较高的计算效率。 和Delaunay图法 。弹簧法需要构建和求解大型 本文以非结构混合网格为背景，对基于径向基 矩阵，计算量比较大，而且主要是针对四面体非结构 函数插值的网格变形技术进行了研究，发展出了一 网格。弹性体方法对任意网格类型都适用，并能很 套实用的网格变形方法。以网格变形方法的实用性 好地保证变形后的网格质量，但是其计算量非常大， 为基本出发点，对基于径向基函数的物面边界位移 难 以在非定常流动计算中多次反复使用。Delaunay 插值条件方程的系数矩阵特性进行了分析，在此基 图法的效率和鲁棒性都优于弹簧法 ，其缺点是：在处 础上提出了一种新的基于贪心法的数据精简方法， 理大变形问题时，Delaunay图可能由于边界节点位 显著提高了网格变形的效率。为了验证本文方法的 置的改变而破坏，因此需要用连续小变形过程进行 有效性和适应性，选取 RAE282