求解可压缩流的高精度非结构网格WENO有限体积法.pdfVIP

计 算 物 理 Vo1．28．NO．5 第28卷 第 5期 CHINESE JOURNAL OFCOMPUTATIONAL PHYSICS Sep．，2011 2011年 9月 文章编号 ：1001-246X(2011)05-06330·8 求解可压缩流的高精度非结构网格 WENO有限体积法 雷国东 ， 李万爱 ， 任玉新 (1．中国航空研究院，北京 100012； 2．清华大学航天航空学 院，北京 100084) 摘 要 ：提 出一种基于最小二乘重构和 WENO限制器 的非结构网格高精度有限体积方法 ．用 中心网格的某些邻 居网格建立重构多项式，给出一定的原则搜索和存储足够多的邻居网格 以建立重构多项式，采用最小二乘法求 解重构多项式的系数．用一种通用的方法控制重构邻居个数 ，以减少存储和计算，采用 WENO限制器和旋转 Riemann求解器 以达到统一的高精度并且抑制守恒律方程求解中的非物理振荡．为检验上述算法，以基于节点的 梯度重构，BathandJesperson限制器 的二阶算法为基准 ，给出三阶和 四阶格式与二阶格式 以及高阶格式若干经典 算例计算结果的对 比和分析． 关键词：最 ／bZ-乘重构 ；高精度；非结构网格；WENO有限体积法；旋转 Riemann求解器 中图分类号：035 文献标识码 ：A 0 引言 WENOFV和 DG等数值格式适合于非结构 网格高精度可压缩流计算 ¨，为 了采用较少的模板抑制 DG 方法的非物理振荡 ，Shu，Luo等提 出HermiteWENO重构解 决这一问题 ；Capdeville 提 出将 Hermite WENO应用于 FVM，以减少重构模板中的邻居个数 ，这种方法类似于 DG，但是 比WENOFVM方法 ，增加了 变量导数方程的求解．而且增加了变量导数 的第一类边界条件．对于特定 问题可行 ，但是对于欧拉方程是难 以给定的．WENOFVM对强间断具有的良好鲁棒性 ，可以采用较大 的CFL数推进 ，较小 的内存消耗 ；而 DG 方法计算量较大 (采用 WENO，HermiteWENO，计算量更大)，稳定性条件更为苛刻 ，随着精度的增加 ，必须采 用越来越小的CFL数，由于必须存储流场各阶导数信息，存储量相对较大 ．一般结构网格上的WENO格 式 采用选定的若干模板，构造 出若干满足局部场信息的多项式，根据多项式的振荡程度选择相应的权系 数，加权构造出能够抑制振荡 的最终重构多项式．其 中模板的选择 ，多项式的重构过程在结构化 网格上进行 比较简单 ，邻居网格 的搜索、存储 比较容易 ；而对于非结构网格，由于其邻居网格 的不规则性，必须采用一 定的原则建立邻居关系．对于每一个非结构网格，其邻居的个数并不确定，为保证不出现欠定问题，邻居模板 个数应该足够多，这时可以采用最小二乘法构造重构多项式 -s]．Mandal，Subramanian。。提 出采用采用加权 最 ／b乘重构 、采用各种限制器、二阶精度的FVM；Dumbser 提出采用 ADER方法 ，将物理坐标转化到计算 坐标 ，最 ／bZ-乘重构的任意阶高精度 FVM，应用于多维问题． 本文采用类似的加权最小二乘重构应用于多维高精度(大于等于三阶)问题，但是不采用文 [5—7]中物 理坐标向计算坐标 的转化 ，所有积分和重构均在物理坐标 (直角坐标)下进行 ，采用 WENO限制器和旋转 Riemann求解器来抑制守恒律方程计算 中出现 的非物理振荡 ，提 出一种通用的邻居 网格选择方法和限制中 心选择方法 ，以保证邻居 网格选择 的有效性和紧凑性．本方法充分利用 FVM对强间断的鲁棒性，最 ／b乘高 阶重构的高分辨率 ，WENO限制器 的稳定性和不降低分辨率特性．三角形网格上 的体积分采用 Hammer积 分，FV方法的体积分对于没有源项的控制方程仅仅出现在最小二乘重构时的参数计算 中，参数结果仅仅与 网格相关 ，这种计算是一次性 的；而反复出现的网格面积分 ，即数值通量的计算 ，可以通过事先给出的积分点 收稿 日期 ：2010一l0—19；修 回日期 ：2011—02—16 作者简介 ：雷国东(1980一)，男 ，湖北嘉鱼 ，主要