Refined isogeometric analysis for generalized Hermitian eigenproblems

09/17/2020
by   Ali Hashemian, et al.
0

We use the refined isogeometric analysis (rIGA) to solve generalized Hermitian eigenproblems (Ku=λ Mu). The rIGA framework conserves the desirable properties of maximum-continuity isogeometric analysis (IGA) discretizations while reducing the computation cost of the solution through partitioning the computational domain by adding zero-continuity basis functions. As a result, rIGA enriches the approximation space and decreases the interconnection between degrees of freedom. We compare computational costs of rIGA versus those of IGA when employing a Lanczos eigensolver with a shift-and-invert spectral transformation. When all eigenpairs within a given interval [λ_s,λ_e] are of interest, we select several shifts σ_k∈[λ_s,λ_e] using a spectrum slicing technique. For each shift σ_k, the cost of factorization of the spectral transformation matrix K-σ_k M drives the total computational cost of the eigensolution. Several multiplications of the operator matrices (K-σ_k M)^-1 M by vectors follow this factorization. Let p be the polynomial degree of basis functions and assume that IGA has maximum continuity of p-1, while rIGA introduces C^0 separators to minimize the factorization cost. For this setup, our theoretical estimates predict computational savings to compute a fixed number of eigenpairs of up to O(p^2) in the asymptotic regime, that is, large problem sizes. Yet, our numerical tests show that for moderately-sized eigenproblems, the total computational cost reduction is O(p). Nevertheless, rIGA improves the accuracy of every eigenpair of the first N_0 eigenvalues and eigenfunctions. Here, we allow N_0 to be as large as the total number of eigenmodes of the original maximum-continuity IGA discretization.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
12/28/2021

Performance of Refined Isogeometric Analysis in Solving Quadratic Eigenvalue Problems

Certain applications that analyze damping effects require the solution o...
research
04/19/2017

A Model Order Reduction Algorithm for Estimating the Absorption Spectrum

The ab initio description of the spectral interior of the absorption spe...
research
06/05/2020

A combined XFEM phase-field computational model for crack growth without remeshing

This paper presents an adaptive strategy for phase-field simulations wit...
research
10/29/2020

A non-asymptotic version of Cressie's refined continuity correction for the binomial distribution

In this paper, we prove a non-asymptotic version of the refined continui...
research
08/01/2019

The Tropical Division Problem and the Minkowski Factorization of Generalized Permutahedra

Given two tropical polynomials f, g on R^n, we provide a characterizatio...
research
05/23/2023

Basis Function Encoding of Numerical Features in Factorization Machines for Improved Accuracy

Factorization machine (FM) variants are widely used for large scale real...

Please sign up or login with your details

Forgot password? Click here to reset