Well-balanced Arbitrary-Lagrangian-Eulerian finite volume schemes on moving nonconforming meshes for the Euler equations of gas dynamics with gravity

被引:46
|
作者
Gaburro, Elena [1 ]
Castro, Manuel J. [2 ]
Dumbser, Michael [1 ]
机构
[1] Univ Trento, Dept Civil Environm & Mech Engn, Via Mesiano 77, I-38123 Trento, Italy
[2] Univ Malaga, Dept Math Stat & Appl Math, Campus Teatinos, E-29071 Malaga, Spain
基金
欧盟地平线“2020”; 欧洲研究理事会;
关键词
accretion; accretion discs; convection; hydrodynamics; instabilities; methods: numerical; NONCONSERVATIVE HYPERBOLIC SYSTEMS; SHALLOW-WATER EQUATIONS; ORDER ADER SCHEMES; CONSERVATION-LAWS; HYDROSTATIC RECONSTRUCTION; NUMERICAL SCHEMES; KINETIC SCHEME; COMBINED FORCE; SOURCE TERMS; SLIDE LINES;
D O I
10.1093/mnras/sty542
中图分类号
P1 [天文学];
学科分类号
0704 ;
摘要
In this work, we present a novel second-order accurate well-balanced arbitrary Lagrangian-Eulerian (ALE) finite volume scheme on moving nonconforming meshes for the Euler equations of compressible gas dynamics with gravity in cylindrical coordinates. The main feature of the proposed algorithm is the capability of preserving many of the physical properties of the system exactly also on the discrete level: besides being conservative for mass, momentum and total energy, also any known steady equilibrium between pressure gradient, centrifugal force, and gravity force can be exactly maintained up to machine precision. Perturbations around such equilibrium solutions are resolved with high accuracy and with minimal dissipation on moving contact discontinuities even for very long computational times. This is achieved by the novel combination of well-balanced path-conservative finite volume schemes, which are expressly designed to deal with source terms written via non-conservative products, with ALE schemes on moving grids, which exhibit only very little numerical dissipation on moving contact waves. In particular, we have formulated a new HLL-type and a novel Osher-type flux that are both able to guarantee the well balancing in a gas cloud rotating around a central object. Moreover, to maintain a high level of quality of the moving mesh, we have adopted a nonconforming treatment of the sliding interfaces that appear due to the differential rotation. A large set of numerical tests has been carried out in order to check the accuracy of the method close and far away from the equilibrium, both, in one-and two-space dimensions.
引用
收藏
页码:2238 / 2262
页数:25
相关论文
共 50 条