Saved in:
| Main Authors: | , , |
|---|---|
| Format: | Preprint |
| Published: |
2025
|
| Subjects: | |
| Online Access: | https://arxiv.org/abs/2505.24504 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866909629194174464 |
|---|---|
| author | Birken, Philipp Dedner, Andreas Klöfkorn, Robert |
| author_facet | Birken, Philipp Dedner, Andreas Klöfkorn, Robert |
| contents | Discontinuous Galerkin (DG) methods are promising high order discretizations for unsteady compressible flows. Here, we focus on Numerical Weather Prediction (NWP). These flows are characterized by a fine resolution in $z$-direction and low Mach numbers, making the system stiff. Thus, implicit time integration is required and for this a fast, highly parallel, low-memory iterative solver for the resulting algebraic systems. As a basic framework, we use inexact Jacobian-Free Newton-GMRES with a preconditioner.
For low order finite volume discretizations, multigrid methods have been successfully applied to steady and unsteady fluid flows. However, for high order DG methods, such solvers are currently lacking. %The lack of efficient solvers suitable for contemporary computer architectures inhibits wider adoption of DG methods. This motivates our research to construct a Jacobian-free precondtioner for high order DG discretizations. The preconditioner is based on a multigrid method constructed for a low order finite volume discretization defined on a subgrid of the DG mesh. We design a computationally efficient and mass conservative mapping between the grids. As smoothers, explicit Runge-Kutta pseudo time iterations are used, which can be implemented in parallel in a Jacobian-free low-memory manner.
We consider DG Methods for the Euler equations and for viscous flow equations in 2D, both with gravity, in a well balanced formulation. Numerical experiments in the software framework DUNE-FEM on atmospheric flow problems show the benefit of this approach. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2505_24504 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Jacobian-free Multigrid Preconditioner for Discontinuous Galerkin Methods applied to Numerical Weather Prediction Birken, Philipp Dedner, Andreas Klöfkorn, Robert Numerical Analysis Computational Physics 65M55, 65M60 Discontinuous Galerkin (DG) methods are promising high order discretizations for unsteady compressible flows. Here, we focus on Numerical Weather Prediction (NWP). These flows are characterized by a fine resolution in $z$-direction and low Mach numbers, making the system stiff. Thus, implicit time integration is required and for this a fast, highly parallel, low-memory iterative solver for the resulting algebraic systems. As a basic framework, we use inexact Jacobian-Free Newton-GMRES with a preconditioner. For low order finite volume discretizations, multigrid methods have been successfully applied to steady and unsteady fluid flows. However, for high order DG methods, such solvers are currently lacking. %The lack of efficient solvers suitable for contemporary computer architectures inhibits wider adoption of DG methods. This motivates our research to construct a Jacobian-free precondtioner for high order DG discretizations. The preconditioner is based on a multigrid method constructed for a low order finite volume discretization defined on a subgrid of the DG mesh. We design a computationally efficient and mass conservative mapping between the grids. As smoothers, explicit Runge-Kutta pseudo time iterations are used, which can be implemented in parallel in a Jacobian-free low-memory manner. We consider DG Methods for the Euler equations and for viscous flow equations in 2D, both with gravity, in a well balanced formulation. Numerical experiments in the software framework DUNE-FEM on atmospheric flow problems show the benefit of this approach. |
| title | Jacobian-free Multigrid Preconditioner for Discontinuous Galerkin Methods applied to Numerical Weather Prediction |
| topic | Numerical Analysis Computational Physics 65M55, 65M60 |
| url | https://arxiv.org/abs/2505.24504 |