Burman, E; Puppi, R; (2022) Two mixed finite element formulations for the weak imposition of the Neumann boundary conditions for the Darcy flow. Journal of Numerical Mathematics , 30 (2) pp. 141-162. 10.1515/jnma-2021-0042.

Preview

Text Burman_10.1515_jnma-2021-0042.pdf Download (650kB) | Preview

Abstract

We propose two different discrete formulations for the weak imposition of the Neumann boundary conditions of the Darcy flow. The Raviart-Thomas mixed finite element on both triangular and quadrilateral meshes is considered for both methods. One is a consistent discretization depending on a weighting parameter scaling as O(h-1), while the other is a penalty-type formulation obtained as the discretization of a perturbation of the original problem and relies on a parameter scaling as O(h-k-1), k being the order of the Raviart-Thomas space. We rigorously prove that both methods are stable and result in optimal convergent numerical schemes with respect to appropriate mesh-dependent norms, although the chosen norms do not scale as the usual L2-norm. However, we are still able to recover the optimal a priori L2-error estimates for the velocity field, respectively, for high-order and the lowest-order Raviart-Thomas discretizations, for the first and second numerical schemes. Finally, some numerical examples validating the theory are exhibited.

Download activity - last month

Export as JSONCSVXML

Download activity - last 12 months

Export as JSONXMLCSV

Downloads by country - last 12 months

Export as JSONCSVXML

Archive Staff Only

View Item