ENO守恒插值(重映)方法及其在流体计算中的应用 Conservative Interpolation(Remapping)Algorithm Based on ENO Interpolation and Application in Computational Fluid Dynamics.pdfVIP

计算物咖理PHYOFCOMPUONALPHYSICS i1 CHINESEJOURNAL 第巧善鬻期年 月TATIN。，j2008 。 [文章编号】1001．246x(2008)06．0641．08 ENO守恒插值(重映)方法及其在 流体计算中的应用 王永健1， 赵 宁2， 王春武1， 王东红1 【摘要】 Non．oscillatory)守恒插值方法的基础上，分析和研究现今流体力学计算中涉及的 在ENO(Essentially 几类网格技术：重叠网格技术、自适应加密技术和运动网格技术．基于ENO插值多项式构造的重映方法具有良 好的守恒性，可以有效保证数据传递中物理量的总体守恒．提出该类守恒插值方法在以上几种网格技术中的一 些应用前景，并给出一些数值算例． [关键词】 ENO插值多项式；守恒插值；重叠网格；自适应加密；运动网格 [中图分类号】0241．5 [文献标识码】 A O 引言 1974年Hirt等人提出了ALE(Arbitrary 网格重分步；③重映步乜。J，即守恒插值．由于二维流体运动中可能出现严重的扭曲现象，因而会造成 换成尽可能规整的新网格．把物理量从啪ge网格上重映到新的调整后的网格上，就是重映过程．它实际 上是利用旧网格上的物理量插值得到新网格上的物理量，通常要求插值过程应具备守恒性和单调性．即插值 过程不引入新的极值，也不改变重映量的积分值． 重映方法主要可分为两类，基于单元相交(CIB／DC，cell．intersection．based face．based (SFB／DC，simplified 较大，常用于Lagrange计算至网格变形较严重时进行重分、重映的情况№“1；后者效率较高，但需要新、旧网 格具有相同的拓扑结构，且对新、旧网格间的差异程度有一定要求，适用于两次重分、重映间隔的Lagrange计 算步数较少，甚至每个Lagrange计算步后都进行重分、重映，对网格变形进行微调的情况№’9’10|． 无论哪类重映方法，为了实现高精度计算，通常都用物理量的分片线性分布代替分片常数分布，这就需 要保证插值过程不产生新的极值点，往往会导致在大梯度区域插值的精度降低，影响计算结果．ENO插值的 能保证插值的单调性，但却允许出现非常微小的振荡，且具有一致高阶精度的优点．本文将ENO插值思想引 人重映过程，在梯度变化较大的区域，选取物理量分布最光滑的点插值，分别给出了基于ENO插值方法的两 类守恒重映方法¨引：基于单元相交的ENO守恒重映算法和基于面通量的ENO守恒重映算法． 守恒插值随ALE方法一起提出，但它的应用并不局限于此．计算流体力学中的重叠网格技术、运动网格 技术中网格数据的传递就需要利用守恒插值．典型的应用问题包括：多体分离、抖振颤振、飞机结冰问题计 算[I耵等．运动网格方法n弘191在一个时间推进步后会根据解的情况，对计算网格进行疏密度的调整；网格自适 应加密方法汹地1在一个时间推进步后，会根据解的情况，对计算网格进行加密或合并处理；重叠网格方法在 [收稿日期】2007—07～04；【修回日期]2007一lI一14 院青年创新基金资助项目 [作者筒介】王永健(1979一)。男．江苏太仓，博士生，从事计算流体力学方面的研究． 万方数据万方数据 计 算 物 理 第25卷 一个时间推进步后各层网格会产生相对运动．这些方法都涉及到新旧网格间的物理量的传递．一般的插值并 不能保证计算区域内物理量的守恒，而且会降低计算精度，以及无法捕捉激波的正确位置㈨，所以对于守恒 性要求很高的问题，守恒插值是唯一的选择． 1 ENO插值多项式构造 ENO插值通过自适应方式选取插值点，虽不能保证插值的单调性，但却允许出现非常微小的振荡，且具 有一致高阶精度的优点．构造的过程如下：设构造控制体0，上的插值多项式，记n；(i=1，2，3)上的网格平 均值为Ui．设在单元ni上，函数，(菇，Y)