两种非结构动网格方法特性研究.pdfVIP

第十三届全国空气弹性学术交流会会议论文集 两种非结构动网格方法特性研究 高传强, 张伟伟, 蒋跃文, 叶正寅 (西北工业大学 翼型叶栅空气动力学国家重点实验室，陕西 西安 710072) 摘 要：基于虚拟结构的弹簧法和基于代数方法的径向基函数 (RBFs)插值方法是两种主要的非结构动网格方法。 本文对两种方法的变形特性进行较详细的定量分析研究。以NACA0012翼型的俯仰运动为例，比较RBFs方法和线 性弹簧法的CPU计算时间和最大变形能力，分析两种方法的网格质量随俯仰角度和运动步数的关系，并给出典型 状态下网格单元质量云图。结果表明，基于RBFs插值的网格变形方法计算效率较高，变形能力较强且能保持较 高的网格质量，是一种高效的动网格方法。 关键字：径向基函数；动网格；弹簧法；非结构网格 引 言 [1] 非定常流动数值模拟中的运动边界的问题是计算流体力学的热点和难点之一 ，运动边界的处 理方法主要是网格变形技术。目前，非结构网格的网格变形方法主要有基于虚拟结构的弹簧法和基 [2] 于代数方法的径向基函数(RBFs)插值方法 。 弹簧法由Batina[3]在研究绕振荡翼型的Euler 流动时提出。他将计算域内的每一条网格线都看作 一根弹簧，整个网格则看成一个弹簧系统。弹簧的刚度系数仅由网格线的长度决定，所以也叫线性 弹簧法或标准弹簧法。Farhat[4]将扭转因子引入弹簧刚度系数，建立了扭转弹簧法。Blom[5]在研究二 维可动边界问题时，通过引入角度约束，提出了半扭转弹簧法。以上改进方法虽然一定程度上提高 了网格的变形能力，但是并没有从根本上解决弹簧法的计算效率较低和数据存储复杂等问题。 径向基函数（RBFs ）插值，作为一种代数插值方法，在各种散乱数据插值方面能取得满意的结 [6] [7] 果 ，近年来开始用于动态网格变形。Boer 首次将RBFs应用于网格变形技术，其基本原理是：运 用RBFs对物面边界节点的位移进行插值，然后利用构造出来的RBFs序列将物面位移效应光滑地分散 到整个网格区域。插值过程仅需要网格节点的坐标，而无需节点间连接信息，所以数据结构简单。 Rendall和 Allen [8-9]提出的基于贪心算法的物面节点精简方法，较大的提高了计算效率。王刚[10]在此 [11] 基础上进一步提出一种更高效的数据精简方法。林言中 等人分析不同基函数和紧支半径对网格质 量和计算效率的影响。 以上研究均只针对某一种方法的介绍或改进，对其网格变形特性做了较充分的研究，但是缺乏 两类方法的直接详尽对比。本文基于以上两类典型的非结构网格变形方法，以NACA0012翼型的二 维非结构网格为例，针对网格变形过程中主要的变形效率，最大变形能力和网格质量等特性做了较 详细的定量分析，总结两种方法在非结构网格变形方面的优缺点。 1.两种网格变形技术的基本方法及网格质量定义 1．1 线性弹簧法 i 对于顶点弹簧法，网格变形过程中网格节点的受力始终认为等于初始状态所受合力。对于节点 ， 假设计算域中共有N i 个节点与之相连，则合力表示为： 518 第十三届全国空气弹性学术交流会会议论文集 N i