High order ADER-DG method with local DG predictor for solutions of differential-algebraic systems of equations

Abstract

A numerical method ADER-DG with a local DG predictor for solving a DAE system has been developed, which was based on the formulation of ADER-DG methods using a local DG predictor for solving ODE and PDE systems. The basis functions were chosen in the form of Lagrange interpolation polynomials with nodal points at the roots of the Radau polynomials, which differs from the classical formulations of the ADER-DG method, where it is customary to use the roots of Legendre polynomials. It was shown that the use of this basis leads to A-stability and L1-stability in the case of using the DAE solver as ODE solver. The numerical method ADER-DG allows one to obtain a highly accurate numerical solution even on very coarse grids, with a step greater than the main characteristic scale of solution variation. The local discrete time solution can be used as a numerical solution of the DAE system between grid nodes, thereby providing subgrid resolution even in the case of very coarse grids. The classical test examples were solved by developed numerical method ADER-DG. With increasing index of the DAE system, a decrease in the empirical convergence orders p is observed. An unexpected result was obtained in the numerical solution of the stiff DAE system -- the empirical convergence orders of the numerical solution obtained using the developed method turned out to be significantly higher than the values expected for this method in the case of stiff problems. It turns out that the use of Lagrange interpolation polynomials with nodal points at the roots of the Radau polynomials is much better suited for solving stiff problems. Estimates showed that the computational costs of the ADER-DG method are approximately comparable to the computational costs of implicit Runge-Kutta methods used to solve DAE systems. Methods were proposed to reduce the computational costs of the ADER-DG method.

98 pages, 44 figures, 21 tables. arXiv admin note: text overlap with arXiv:2409.09933

References in corpus (19)

• A Posteriori Subcell Limiting of the Discontinuous Galerkin Finite Element Method for Hyperbolic Conservation Laws
• Very High Order $\PNM$ Schemes on Unstructured Meshes for the Resistive Relativistic MHD Equations
• ADER-WENO Finite Volume Schemes with Space-Time Adaptive Mesh Refinement
• Space-time adaptive ADER discontinuous Galerkin finite element schemes with a posteriori sub-cell finite volume limiting
• High order direct Arbitrary-Lagrangian-Eulerian schemes on moving Voronoi meshes with topology changes
• High Order Space-Time Adaptive WENO Finite Volume Schemes for Non-Conservative Hyperbolic Systems
• Arbitrary-Lagrangian-Eulerian discontinuous Galerkin schemes with a posteriori subcell finite volume limiting on moving unstructured meshes
• ExaHyPE: An Engine for Parallel Dynamically Adaptive Simulations of Wave Problems
• Space-time adaptive ADER-DG schemes for dissipative flows: compressible Navier-Stokes and resistive MHD equations
• Conformal and covariant Z4 formulation of the Einstein equations: strongly hyperbolic first-order reduction and solution with discontinuous Galerkin schemes
• Space-time adaptive ADER discontinuous Galerkin schemes for nonlinear hyperelasticity with material failure
• A comparison of high order explicit Runge-Kutta, extrapolation, and deferred correction methods in serial and parallel
• On the Eigenvalues of the ADER-WENO Galerkin Predictor
• Approximate Taylor methods for ODEs
• An Efficient ADER-DG Local Time Stepping Scheme for 3D HPC Simulation of Seismic Waves in Poroelastic Media
• Efficient iterative arbitrary high order methods: an adaptive bridge between low and high order
• Space-time adaptive ADER-DG finite element method with LST-DG predictor and a posteriori sub-cell WENO finite-volume limiting for simulation of non-stationary compressible multicomponent reactive flows
• Space-time adaptive ADER-DG finite element method with LST-DG predictor and a posteriori sub-cell ADER-WENO finite-volume limiting for multidimensional detonation waves simulation
• Arbitrary high order ADER-DG method with local DG predictor for solutions of initial value problems for systems of first-order ordinary differential equations

Cited by in corpus (1)

• Improving the local solution of the DG predictor of the ADER-DG method for solving systems of ordinary differential equations and its applicability to systems of differential-algebraic equations