A polytopal method for the Brinkman problem robust in all regimes
Abstract
In this work we develop a discretisation method for the Brinkman problem that is uniformly well-behaved in all regimes (as identified by a local dimensionless number with the meaning of a friction coefficient) and supports general meshes as well as arbitrary approximation orders. The method is obtained combining ideas from the Hybrid High-Order and Discrete de Rham methods, and its robustness rests on a potential reconstruction and stabilisation terms that change in nature according to the value of the local friction coefficient. We derive error estimates that, thanks to the presence of cut-off factors, are valid across the all regimes and provide extensive numerical validation.
- Publication:
-
Computer Methods in Applied Mechanics and Engineering
- Pub Date:
- May 2023
- DOI:
- 10.1016/j.cma.2023.115981
- arXiv:
- arXiv:2301.03272
- Bibcode:
- 2023CMAME.409k5981D
- Keywords:
-
- Mathematics - Numerical Analysis;
- 65N30;
- 65N08;
- 76S05;
- 76D07
- E-Print:
- doi:10.1016/j.cma.2023.115981