基于WENO重构的高阶移动网格动理学格式.pdfVIP

第30卷第4期 计 算 物 理 V01．30．No．4 CHINESE OFCOMPUTATIONALPHYSICS 2013年7月 JOURNAL July，2013 文章编号：1001．246X(2013)04-0509-06 基于WENO重构的高阶移动网格动理学格式 徐喜华1， 倪国喜2 (1．中国工程物理研究院研究生部，北京100088；2．计算物理实验室，北京应用物理与计算数学研究所，北京100088) 摘要：提出一种基于WENO重构的高阶(至少三阶)移动网格动理学格式．利用流体力学方程的积分形式得到 移动网格上离散格式，再利用自适应移动网格方法移动网格，进而得到网格速度，利用WENO重构得到高阶插值 多项式，最后使用时间方向上精确的动理学数值方法构造数值通量，得到移动网格单元上新的物理量．数值实验 表明这种格式同时具有高精度、高分辨率的特点． 关键词：可压流体力学方程；WENO重构；动理学格式；自适应移动网格方法 中图分类号：0351 文献标识码：A 0 引言 在计算流体力学中，高精度、高分辨率格式是一个热点课题，在很多流体力学问题的数值模拟中，高精 度数值方法能得到问题的真实解．近几年来流行的高阶格式主要有Essential Discontinuous Volume Essential methods(SV)‘3 Galerkin(DG)‘2|、Spectral 1及Weighted 的模板信息，在WENO格式通过引入变化的加权因子，使格式在光滑区解的截断误差阶数有进一步的提 高，而在间断区附近仍然保持ENO格式良好的分辨能力． 数值方法的精度除与格式本身有关以外，还与网格的分布有关，在现有格式的基础上提高数值解的精 度，自适应移动网格技术是一种有效的方法，可消除由于物理量变化剧烈而引起的解的震荡，降低数值耗 散，在提高精度的同时大幅度减少节点数，提高计算效率，这方面已取得了许多成果¨’71．例如Tang¨1等 人提出一种用基于变分原理的自适应移动网格方法，使网格分布自适应地跟踪间断，并给出了物理量的守 3等人发展了一套移动网格上的动理学格式．该方 恒的插值方法，提高数值模拟的精度，在此基础上，NiL8 法结合了ALE的思想，摒弃物理量在新旧网格之间的重映，从而减少了数值误差． 本文给出基于WENO重构的高阶(至少三阶)移动网格动理学格式．首先在移动框架下得到流体力学 方程的离散格式，利用自适应移动网格方法得到节点速度，对单元物理量进行WENO重构，再对高阶插值 点构造动理学数值通量．数值实验表明此格式不仅在光滑区域有很高分辨率，达到高阶精度，保证网格的 质量，并且在激波附近有很高的分辨率． 1 移动网格上的流体力学方程的离散格式 在控制体门(￡)上边界速度为U。的无粘可压缩流体力学方程的积分形式为 收稿日期：2012—10—18；修回日期：2013—0l一25 基金资助项目 作者简介：徐喜华(1986一)，男，博士生，从事计算流体力学方法研究，E-mail：xihuaxu@gmail．com 万方数据 510 计 算 物 理 第30卷 ——————————————————————————————————————————————————一 瓤㈤pd力=一kp(U一¨’础， 杀k删n=一如u(u一㈦‘础一b础， 。’ 未L删n=一kpE(U—Vg)刊s—kp¨蒯s， 度．当U。：0，方程组(1)是欧拉坐标的积分形式．当U。=u，则为拉