Structure preserving nodal continuous Finite Elements via Global Flux quadrature
Fuente:
arXiv
Saved in:
| Main Authors: | , , |
|---|---|
| Format: | Preprint |
| Published: |
2024
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866908361027485696 |
|---|---|
| author | Barsukow, Wasilij Ricchiuto, Mario Torlo, Davide |
| author_facet | Barsukow, Wasilij Ricchiuto, Mario Torlo, Davide |
| contents | Numerical methods for hyperbolic PDEs require stabilization. For linear acoustics, divergence-free vector fields should remain stationary, but classical Finite Difference methods add incompatible diffusion that dramatically restricts the set of discrete stationary states of the numerical method. Compatible diffusion should vanish on stationary states, e.g. should be a gradient of the divergence. Some Finite Element methods allow to naturally embed this grad-div structure, e.g. the SUPG method or OSS. We prove here that the particular discretization associated to them still fails to be constraint preserving. We then introduce a new framework on Cartesian grids based on surface (volume in 3D) integrated operators inspired by Global Flux quadrature and related to mimetic approaches. We are able to construct constraint-compatible stabilization operators (e.g. of SUPG-type) and show that the resulting methods are vorticity-preserving. We show that the Global Flux approach is even super-convergent on stationary states, we characterize the kernels of the discrete operators and we provide projections onto them. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2407_10579 |
| institution | arXiv |
| publishDate | 2024 |
| record_format | arxiv |
| spellingShingle | Structure preserving nodal continuous Finite Elements via Global Flux quadrature Barsukow, Wasilij Ricchiuto, Mario Torlo, Davide Numerical Analysis Numerical methods for hyperbolic PDEs require stabilization. For linear acoustics, divergence-free vector fields should remain stationary, but classical Finite Difference methods add incompatible diffusion that dramatically restricts the set of discrete stationary states of the numerical method. Compatible diffusion should vanish on stationary states, e.g. should be a gradient of the divergence. Some Finite Element methods allow to naturally embed this grad-div structure, e.g. the SUPG method or OSS. We prove here that the particular discretization associated to them still fails to be constraint preserving. We then introduce a new framework on Cartesian grids based on surface (volume in 3D) integrated operators inspired by Global Flux quadrature and related to mimetic approaches. We are able to construct constraint-compatible stabilization operators (e.g. of SUPG-type) and show that the resulting methods are vorticity-preserving. We show that the Global Flux approach is even super-convergent on stationary states, we characterize the kernels of the discrete operators and we provide projections onto them. |
| title | Structure preserving nodal continuous Finite Elements via Global Flux quadrature |
| topic | Numerical Analysis |
| url | https://arxiv.org/abs/2407.10579 |