Finite element and NURBS approximations of eigenvalue, boundary-value, and initial-value problems

Abstract We study the spectral approximation properties of finite element and NURBS spaces from a global perspective. We focus on eigenfunction approximations and discover that the L 2 -norm errors for finite element eigenfunctions exhibit pronounced “spikes” about the transition points between branches of the eigenvalue spectrum. This pathology is absent in NURBS approximations. By way of the Pythagorean eigenvalue error theorem, we determine that the squares of the energy-norm errors of the eigenfunctions are the sums of the eigenvalue errors and the squares of the L 2 -norm eigenfunction errors. The spurious behavior of the higher eigenvalues for standard finite elements is well-known and therefore inherited by the energy-norm errors along with the spikes in the L 2 -norm of the eigenfunction errors. The eigenvalue pathology is absent for NURBS. The implications of these results to the corresponding elliptic boundary-value problem and parabolic and hyperbolic initial-value problems are discussed.

[1]  R. D. Richtmyer,et al.  Difference methods for initial-value problems , 1959 .

[2]  M. Okrouhlík,et al.  Complex wavenumber Fourier analysis of the B-spline based finite element method , 2014 .

[3]  G. M.,et al.  Partial Differential Equations I , 2023, Applied Mathematical Sciences.

[4]  Jesús Ildefonso Díaz Díaz,et al.  ON THE COMPLEX GINZBURG–LANDAU EQUATION WITH A DELAYED FEEDBACK , 2006 .

[5]  T. Hughes,et al.  ISOGEOMETRIC ANALYSIS: APPROXIMATION, STABILITY AND ERROR ESTIMATES FOR h-REFINED MESHES , 2006 .

[6]  Thomas J. R. Hughes,et al.  n-Widths, sup–infs, and optimality ratios for the k-version of the isogeometric finite element method , 2009 .

[7]  Victor M. Calo,et al.  The role of continuity in residual-based variational multiscale modeling of turbulence , 2007 .

[8]  R. D. Richtmyer,et al.  Difference methods for initial-value problems , 1959 .

[9]  Alessandro Reali,et al.  Isogeometric Analysis of Structural Vibrations , 2006 .

[10]  Jintai Chung,et al.  A Time Integration Algorithm for Structural Dynamics With Improved Numerical Dissipation: The Generalized-α Method , 1993 .

[11]  P. Pinsky,et al.  Complex wavenumber Fourier analysis of the p-version finite element method , 1994 .

[12]  Mark Ainsworth,et al.  Discrete Dispersion Relation for hp-Version Finite Element Approximation at High Wave Number , 2004, SIAM J. Numer. Anal..

[13]  Thomas J. R. Hughes,et al.  Improved numerical dissipation for time integration algorithms in structural dynamics , 1977 .

[14]  G. Strang,et al.  An Analysis of the Finite Element Method , 1974 .

[15]  Thomas J. R. Hughes,et al.  Isogeometric Analysis: Toward Integration of CAD and FEA , 2009 .

[16]  Alessandro Reali,et al.  AN ISO GEOMETRIC ANALYSIS APPROACH FOR THE STUDY OF STRUCTURAL VIBRATIONS , 2006 .

[17]  Alessandro Reali,et al.  Duality and unified analysis of discrete approximations in structural dynamics and wave propagation : Comparison of p-method finite elements with k-method NURBS , 2008 .

[18]  John A. Evans,et al.  Discrete spectrum analyses for various mixed discretizations of the Stokes eigenproblem , 2012 .