On the Active Flux scheme for hyperbolic PDEs with source terms

The Active Flux scheme is a Finite Volume scheme with additional point values distributed along the cell boundary. It is third order accurate and does not require a Riemann solver: the continuous reconstruction serves as initial data for the evolution of the points values. The intercell flux is then obtained from the evolved values along the cell boundary by quadrature. This paper focuses on the conceptual extension of Active Flux to include source terms, and thus for simplicity assumes the homogeneous part of the equations to be linear. To a large part, the treatment of the source terms is independent of the choice of the homogeneous part of the system. Additionally, only systems are considered which admit characteristics (instead of characteristic cones). This is the case for scalar equations in any number of spatial dimensions and systems in one spatial dimension. Here, we succeed to extend the Active Flux method to include (possibly nonlinear) source terms while maintaining third order accuracy of the method. This requires a novel (approximate) operator for the evolution of point values and a modified update procedure of the cell average. For linear acoustics with gravity, it is shown how to achieve a well-balanced / stationarity preserving numerical method.

[1]  J. Greenberg,et al.  A well-balanced scheme for the numerical processing of source terms in hyperbolic equations , 1996 .

[2]  Christian Klingenberg,et al.  A General Well-Balanced Finite Volume Scheme for Euler Equations with Gravity , 2016 .

[3]  Timothy A. Eymann,et al.  Multidimensional Active Flux Schemes , 2013 .

[4]  Wasilij Barsukow,et al.  Stationarity preserving schemes for multi-dimensional linear systems , 2018, Math. Comput..

[5]  J. Christensen.  Doctoral thesis , 1970 .

[6]  D. Fan,et al.  On the Acoustic Component of Active Flux Schemes for Nonlinear Hyperbolic Conservation Laws , 2017 .

[7]  P. Roe,et al.  Investigations of a new scheme for wave propagation , 2015 .

[8]  Christian Klingenberg,et al.  A well‐balanced scheme to capture non‐explicit steady states in the Euler equations with gravity , 2016 .

[9]  Eitan Tadmor,et al.  Well-balanced schemes for the Euler equations with gravitation: Conservative formulation using global fluxes , 2017, J. Comput. Phys..

[10]  J. P. Berberich,et al.  High order discretely well-balanced finite volume methods for Euler equations with gravity - without any à priori information about the hydrostatic solution , 2020, ArXiv.

[11]  C. Caldwell.  Mathematics of Computation , 1999 .

[12]  Emmanuel Audusse,et al.  A Fast and Stable Well-Balanced Scheme with Hydrostatic Reconstruction for Shallow Water Flows , 2004, SIAM J. Sci. Comput..

[13]  Philip L. Roe,et al.  Third-order active-flux scheme for advection diffusion: Hyperbolic diffusion, boundary condition, and Newton solver , 2016 .

[14]  M. Mulholland,et al.  The University of Michigan. , 2003, Archives of surgery.

[15]  Christian Klingenberg,et al.  High order well-balanced finite volume methods for multi-dimensional systems of hyperbolic balance laws , 2019, 1903.05154.

[16]  Eric F Darve,et al.  Author ' s personal copy A hybrid method for the parallel computation of Green ’ s functions , 2009 .

[17]  Christiane Helzel,et al.  A New ADER Method Inspired by the Active Flux Method , 2019, Journal of Scientific Computing.

[18]  Christian Klingenberg,et al.  Second Order Finite Volume Scheme for Euler Equations with Gravity which is Well-Balanced for General Equations of State and Grid Systems , 2018, Communications in Computational Physics.

[19]  Alfredo Bermúdez,et al.  Upwind methods for hyperbolic conservation laws with source terms , 1994 .

[20]  Otto R. Spies,et al.  Tables of integral transforms, volume 2: edited by A. Erdelyi. 451 pages, 16 × 24 cm. New York, McGraw-Hill Book Co., Inc., 1954. Price, $8.00. , 1955 .

[21]  Wasilij Barsukow,et al.  The Active Flux Scheme for Nonlinear Problems , 2020, Journal of Scientific Computing.

[22]  Randall J. LeVeque,et al.  Balancing Source Terms and Flux Gradients in High-Resolution Godunov Methods , 1998 .

[23]  Christian Klingenberg,et al.  A Well-Balanced Scheme for the Euler Equation with a Gravitational Potential , 2014 .

[24]  Randall J. LeVeque,et al.  Tsunami modelling with adaptively refined finite volume methods* , 2011, Acta Numerica.

[25]  B. Vanleer,et al.  Towards the ultimate conservative difference scheme. IV. A new approach to numerical convection , 1977 .

[26]  C. Klingenberg,et al.  Well-balanced treatment of gravity in astrophysical fluid dynamics simulations at low Mach numbers , 2021, Astronomy & Astrophysics.

[27]  Philip L. Roe,et al.  Active Flux Schemes for Systems , 2011 .

[28]  Philip L. Roe,et al.  The Active Flux Scheme on Cartesian Grids and Its Low Mach Number Limit , 2018, Journal of Scientific Computing.

[29]  Christian Klingenberg,et al.  A Second Order Well-Balanced Finite Volume Scheme for Euler Equations with Gravity , 2015, SIAM J. Sci. Comput..

[30]  Siddhartha Mishra,et al.  A well-balanced finite volume scheme for the Euler equations with gravitation - The exact preservation of hydrostatic equilibrium with arbitrary entropy stratification , 2016 .