A locally-conservative proximal Galerkin method for pointwise bound constraints
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_ | 1866908540286795776 |
|---|---|
| author | Fu, Guosheng Keith, Brendan Masri, Rami |
| author_facet | Fu, Guosheng Keith, Brendan Masri, Rami |
| contents | We introduce the first-order system proximal Galerkin (FOSPG) method, a locally mass-conserving, hybridizable finite element method for solving heterogeneous anisotropic diffusion and obstacle problems. Like other proximal Galerkin methods, FOSPG finds solutions by solving a recursive sequence of smooth, discretized, nonlinear subproblems. We establish the well-posedness and convergence of these nonlinear subproblems along with stability and error estimates under low regularity assumptions for the linearized equations obtained by solving each subproblem using Newton's method. The FOSPG method exhibits several advantages, including high-order accuracy, discrete maximum principle or bound-preserving discrete solutions, and local mass conservation. It also achieves prescribed solution accuracy within asymptotically mesh-independent numbers of subproblems and linear solves per subproblem iteration. Numerical experiments on benchmarks for anisotropic diffusion and obstacle problems confirm these attributes. Furthermore, an open-source implementation of the method is provided to facilitate broader adoption and reproducibility. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2412_21039 |
| institution | arXiv |
| publishDate | 2024 |
| record_format | arxiv |
| spellingShingle | A locally-conservative proximal Galerkin method for pointwise bound constraints Fu, Guosheng Keith, Brendan Masri, Rami Numerical Analysis 35J86, 49J40, 65N30 We introduce the first-order system proximal Galerkin (FOSPG) method, a locally mass-conserving, hybridizable finite element method for solving heterogeneous anisotropic diffusion and obstacle problems. Like other proximal Galerkin methods, FOSPG finds solutions by solving a recursive sequence of smooth, discretized, nonlinear subproblems. We establish the well-posedness and convergence of these nonlinear subproblems along with stability and error estimates under low regularity assumptions for the linearized equations obtained by solving each subproblem using Newton's method. The FOSPG method exhibits several advantages, including high-order accuracy, discrete maximum principle or bound-preserving discrete solutions, and local mass conservation. It also achieves prescribed solution accuracy within asymptotically mesh-independent numbers of subproblems and linear solves per subproblem iteration. Numerical experiments on benchmarks for anisotropic diffusion and obstacle problems confirm these attributes. Furthermore, an open-source implementation of the method is provided to facilitate broader adoption and reproducibility. |
| title | A locally-conservative proximal Galerkin method for pointwise bound constraints |
| topic | Numerical Analysis 35J86, 49J40, 65N30 |
| url | https://arxiv.org/abs/2412.21039 |