一维双曲守恒律的龙格-库塔控制体积间断有限元方法

计算物理 ›› 2009, Vol. 26 ›› Issue (4): 501-509.

一维双曲守恒律的龙格-库塔控制体积间断有限元方法

陈大伟¹, 蔚喜军²

1. 中国工程物理研究院研究生部, 北京 100088;
2. 北京应用物理与计算数学研究所计算物理实验室, 北京 100088

收稿日期:2008-03-12 修回日期:2008-07-29 出版日期:2009-07-25 发布日期:2009-07-25
作者简介:陈大伟(1982-),男,吉林公主岭,硕士生,从事计算流体力学方面的研究,北京2101信箱100088.
基金资助:
国家自然科学基金(10471011,10771019)资助项目

RKCVDFEM for One-dimensional Hyperbolic Conservation Laws

CHEN Dawei¹, YU Xijun²

1. Graduate School of China Academy of Engineering Physics, Beijing 100088, China;
2. Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Beijing 100088, China

Received:2008-03-12 Revised:2008-07-29 Online:2009-07-25 Published:2009-07-25

摘要/Abstract

摘要： 给出数值求解一维双曲守恒律方程的新方法——龙格-库塔控制体积间断有限元方法(RKCVDFEM),其中空间离散基于控制体积有限元方法,时间离散基于二阶TVB Runge-Kutta技术,有限元空间选取为分段线性函数空间.理论分析表明,格式具有总变差有界(TVB)的性质,而且空间和时间离散形式上具有二阶精度.数值算例表明,数值解收敛到熵解并且对光滑解的收敛阶是最优的,优于龙格-库塔间断Galerkin方法(RKDGM)的计算结果.

关键词: 双曲守恒律, 龙格-库塔技术, 控制体积有限元方法

Abstract: Runge-Kutta control volume discontinuous finite element method (RKCVDFEM) is proposed to solve numerically hyperbolic conservation laws,in which space discretization is based on control volume finite element method (CVFEM) while time discretization is based on a second order accurate TVB Runge-Kutta technique.Piecewise linear function space is chosen as finite element space.The scheme is total variation bounded (TVB) and is formally second order accurate in space and time.Numerical examples show that numerical solution converges to the entropy solution,and order of convergence is optimal for smooth solution.Compared with numerical results of Runge-Kutta discontinuous Galerkin method (RKDGM) those of RKCVDFEM are better.

Key words: hyperbolic conservation laws, Runge-Kutta technique, control volume finite element method

中图分类号:

陈大伟, 蔚喜军. 一维双曲守恒律的龙格-库塔控制体积间断有限元方法[J]. 计算物理, 2009, 26(4): 501-509.

CHEN Dawei, YU Xijun. RKCVDFEM for One-dimensional Hyperbolic Conservation Laws[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2009, 26(4): 501-509.

[1]	郑素佩, 建芒芒, 封建湖, 翟梦情. 保号WENO-AO型中心迎风格式[J]. 计算物理, 2022, 39(6): 677-686.
[2]	郭子滔, 冯仁忠. 一种高精度的修正Hermite-ENO格式[J]. 计算物理, 2019, 36(2): 141-152.
[3]	赵丰祥, 潘亮, 王双虎. 基于非结构四边形网格的WENO格式[J]. 计算物理, 2018, 35(5): 525-534.
[4]	程晓晗, 聂玉峰, 蔡力, 封建湖. 基于移动网格的熵稳定格式[J]. 计算物理, 2017, 34(2): 175-182.
[5]	任炯, 封建湖, 刘友琼, 梁楠. 求解双曲守恒律方程的高分辨率熵相容格式[J]. 计算物理, 2014, 31(5): 539-551.
[6]	冯涛, 蔚喜军, 安恒斌, 崔霞, 吴迪, 李珍珍. Lax-Wendroff时间离散的自适应间断有限元方法求解三维可压缩欧拉方程[J]. 计算物理, 2013, 30(6): 791-798.
[7]	吴迪, 蔚喜军, 徐云. 局部时间步长间断有限元方法求解三维欧拉方程[J]. 计算物理, 2011, 28(1): 1-9.
[8]	罗力, 封建湖, 唐小娟, 向量. 求解双曲守恒律方程的高分辨率熵稳定格式[J]. 计算物理, 2010, 27(5): 671-678.
[9]	吴迪, 蔚喜军. 自适应间断有限元方法求解三维欧拉方程[J]. 计算物理, 2010, 27(4): 492-500.
[10]	徐云, 蔚喜军. 自适应间断有限元方法求解双曲守恒律方程[J]. 计算物理, 2009, 26(2): 159-168.
[11]	蔚喜军, 符鸿源, 常谦顺. 用有限元方法求解双曲守恒律[J]. 计算物理, 1999, 16(5): 457-466.
[12]	朱庆勇, 马延文. 求解双曲型守恒律方程的高精度迎风紧致群速度控制法[J]. 计算物理, 1998, 15(5): 531-536.