A boundary penalization technique to remove outliers from isogeometric analysis on tensor-product meshes

Abstract

We introduce a boundary penalization technique to improve the spectral approximation of isogeometric analysis (IGA). The technique removes the outliers appearing in the high-frequency region of the approximate spectrum when using the $C^{p-1},p$-th ($p\ge 3$) order isogeometric elements. We focus on the classical Laplacian (Dirichlet) eigenvalue problem in 1D to illustrate the idea and then use the tensor-product structure to generate the stiffness and mass matrices for multiple dimensional problems. To remove the outliers, we penalize the product of the higher-order derivatives from both the solution and test spaces at the domain boundary. Intuitively, we construct a better approximation by weakly imposing features of the exact solution. Effectively, we add terms to the variational formulation at the boundaries with minimal extra computational cost. We then generalize the idea to remove the outliers for the isogeometric analysis to the Neumann eigenvalue problem (for $p\ge 2$). The boundary penalization does not change the test and solution spaces. In the limiting case when the penalty goes to infinity, we perform the dispersion analysis of $C^2$ cubic elements for Dirichlet eigenvalue problem and $C^1$ quadratic elements for Neumann eigenvalue problem. We obtain the analytical eigenpairs for the resulting matrix eigenvalue problems. Numerical experiments show optimal convergence rates for the eigenvalues and eigenfunctions of the discrete operator.

30 pages

Cited by in corpus (13)

• A mathematical theory for mass lumping and its generalization with applications to isogeometric analysis
• Towards higher-order accurate mass lumping in explicit isogeometric analysis for structural dynamics
• Critical time-step size analysis and mass scaling by ghost-penalty for immersogeometric explicit dynamics
• Variationally consistent mass scaling for explicit time-integration schemes of lower- and higher-order finite element methods
• A variational approach based on perturbed eigenvalue analysis for improving spectral properties of isogeometric multipatch discretizations
• SoftIGA: soft isogeometric analysis
• A theoretical study on the effect of mass lumping on the discrete frequencies in immersogeometric analysis
• Application of optimal spline subspaces for the removal of spurious outliers in isogeometric discretizations
• Nonlinear dynamic analysis of shear- and torsion-free rods using isogeometric discretization and outlier removal
• Outlier removal for isogeometric spectral approximation with the optimally-blended quadratures
• Analytical solutions to some generalized and polynomial eigenvalue problems
• Outlier-free isogeometric discretizations for Laplace eigenvalue problems: closed-form eigenvalue and eigenvector expressions
• SoftFEM: revisiting the spectral finite element approximation of second-order elliptic operators