Search

This paper proposes a low-cost, penalty parameter-free, and pressure-robust Stokes solver based on the enriched Galerkin (EG) method with a discontinuous velocity enrichment function. The EG method employs the interior penalty discontinuous Galerkin (IPDG) formulation to weakly impose the continuity of the velocity function. However, despite its advantage of symmetry, the symmetric IPDG formulation requires a lot of computational effort to choose an optimal penalty parameter and compute different trace terms. To reduce such effort, we replace the derivatives of the velocity function with its weak derivatives computed by the geometric data of elements. Therefore, our modified EG (mEG) method is a penalty parameter-free numerical scheme that has reduced computational complexity and conserves the optimal convergence orders. Moreover, we achieve pressure robustness for the mEG method by employing a velocity reconstruction operator on the load vector on the right-hand side of the discrete system. The theoretical results are confirmed through numerical experiments with two- and three-dimensional examples. [ABSTRACT FROM AUTHOR]