Stationarity preserving nodal Finite Element methods for multi-dimensional linear hyperbolic balance laws via a Global Flux quadrature formulation
Fuente:
arXiv
Enregistré dans:
| Auteurs principaux: | , , |
|---|---|
| Format: | Preprint |
| Publié: |
2025
|
| Sujets: | |
| Accès en ligne: | |
| Tags: |
Ajouter un tag
Pas de tags, Soyez le premier à ajouter un tag!
|
| _version_ | 1866914072616763392 |
|---|---|
| author | Barsukow, Wasilij Ricchiuto, Mario Torlo, Davide |
| author_facet | Barsukow, Wasilij Ricchiuto, Mario Torlo, Davide |
| contents | We consider linear, hyperbolic systems of balance laws in several space dimensions. They possess non-trivial steady states, which result from the equilibrium between derivatives of the unknowns in different directions, and the sources. Standard numerical methods fail to account for this equilibrium, and include stabilization that destroys it. This manifests itself in a diffusion of states that are supposed to remain stationary. We derive new stabilized high-order Finite Element methods based on a Global Flux quadrature: we reformulate the entire spatial operator as a mixed derivative of a single quantity, referred to as global flux. All spatial derivatives and the sources are thus treated simultaneously, and our methods are stationarity preserving. Additionally, when this formulation is combined with interpolation on Gauss-Lobatto nodes, the new methods are super-convergent at steady state. Formal consistency estimates, and strategies to construct well-prepared initial data are provided. The numerical results confirm the theoretical predictions, and show the tremendous benefits of the new formulation. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2510_02928 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Stationarity preserving nodal Finite Element methods for multi-dimensional linear hyperbolic balance laws via a Global Flux quadrature formulation Barsukow, Wasilij Ricchiuto, Mario Torlo, Davide Numerical Analysis We consider linear, hyperbolic systems of balance laws in several space dimensions. They possess non-trivial steady states, which result from the equilibrium between derivatives of the unknowns in different directions, and the sources. Standard numerical methods fail to account for this equilibrium, and include stabilization that destroys it. This manifests itself in a diffusion of states that are supposed to remain stationary. We derive new stabilized high-order Finite Element methods based on a Global Flux quadrature: we reformulate the entire spatial operator as a mixed derivative of a single quantity, referred to as global flux. All spatial derivatives and the sources are thus treated simultaneously, and our methods are stationarity preserving. Additionally, when this formulation is combined with interpolation on Gauss-Lobatto nodes, the new methods are super-convergent at steady state. Formal consistency estimates, and strategies to construct well-prepared initial data are provided. The numerical results confirm the theoretical predictions, and show the tremendous benefits of the new formulation. |
| title | Stationarity preserving nodal Finite Element methods for multi-dimensional linear hyperbolic balance laws via a Global Flux quadrature formulation |
| topic | Numerical Analysis |
| url | https://arxiv.org/abs/2510.02928 |