English

Stabilization of isogeometric finite element method with optimal test functions computed from $L_2$ norm residual minimization

Numerical Analysis 2024-11-26 v1 Numerical Analysis

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 L2L_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 L2L_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.

Keywords

Cite

@article{arxiv.2411.15565,
  title  = {Stabilization of isogeometric finite element method with optimal test functions computed from $L_2$ norm residual minimization},
  author = {Marcin Łoś and Tomasz Służalec and Maciej Paszyński and Eirik Valseth},
  journal= {arXiv preprint arXiv:2411.15565},
  year   = {2024}
}

Comments

22 pages, 14 figures, 12 tables