Stabilization of isogeometric finite element method with optimal test functions computed from $L_2$ norm residual minimization
Abstract
We compare several stabilization methods in the context of isogeometric analysis and B-spline basis functions, using an advection-dominated advection\revision{-}diffusion as a model problem. We derive (1) the least-squares finite element method formulation using the framework of Petrov-Galerkin method with optimal test functions in the $L_2$ norm, which guarantee automatic preservation of the \emph{inf-sup} condition of the continuous formulation. We also combine it with the standard Galerkin method to recover (2) the Galerkin/least-squares formulation, and derive coercivity constant bounds valid for B-spline basis functions. The resulting stabilization method are compared with the least-squares and (3) the Streamline-Upwind Petrov-Galerkin (SUPG)method using again the Eriksson-Johnson model problem. The results indicate that least-squares (equivalent to Petrov-Galerkin with $L_2$-optimal test functions) outperforms the other stabilization methods for small P\'eclet numbers, while strongly advection-dominated problems are better handled with SUPG or Galerkin/least-squares.
- Publication:
-
arXiv e-prints
- Pub Date:
- November 2024
- arXiv:
- arXiv:2411.15565
- Bibcode:
- 2024arXiv241115565L
- Keywords:
-
- Mathematics - Numerical Analysis
- E-Print:
- 22 pages, 14 figures, 12 tables