hp-Adaptive composite discontinuous Galerkin methods for elliptic eigenvalue problems on complicated domains

In this paper we develop the a posteriori error estimation of hp-adaptive discontinuous Galerkin composite finite element methods (DGFEMs) for the discretization of second-order elliptic eigenvalue problems. DGFEMs allow for the approximation of problems posed on computational domains which may contain local geometric features. The dimension of the composite finite element space is independent of the number of geometric features. This is in contrast with standard finite element methods, as the minimal number of elements needed to represent the underlying domain can be very large and so the dimension of the finite element space. Computable upper bounds on the error for both eigenvalues and eigenfunctions are derived. Numerical experiments highlighting the practical application of the proposed estimators within an automatic hp-adaptive refinement procedure will be presented

Spline-based Rayleigh-Ritz methods for the approximation of the natural modes of vibration for flexible beams with tip bodies

Rayleigh-Ritz methods for the approximation of the natural modes for a class of vibration problems involving flexible beams with tip bodies using subspaces of piecewise polynomial spline functions are developed. An abstract operator theoretic formulation of the eigenvalue problem is derived and spectral properties investigated. The existing theory for spline-based Rayleigh-Ritz methods applied to elliptic differential operators and the approximation properties of interpolatory splines are useed to argue convergence and establish rates of convergence. An example and numerical results are discussed

Detection of holes in an elastic body based on eigenvalues and traces of eigenmodes

We consider the numerical solution of an inverse problem of finding the shape and location of holes in an elastic body. The problem is solved by minimizing a functional depending on the eigenvalues and traces of corresponding eigenmodes. We use the adjoint method to calculate the shape derivative of this functional. The optimization is performed by BFGS, using a genetic algorithm as a preprocessor and the Method of Fundamental Solutions as a solver for the direct problem. We address several numerical simulations that illustrate the good performance of the method.info:eu-repo/semantics/publishedVersio

On Finite Element Methods for 2nd order (semi–) periodic Eigenvalue Problems

We deal with a class of elliptic eigenvalue problems (EVPs) on a rectangle Ω ⊂ R^2 , with periodic or semi–periodic boundary conditions (BCs) on ∂Ω. First, for both types of EVPs, we pass to a proper variational formulation which is shown to fit into the general framework of abstract EVPs for symmetric, bounded, strongly coercive bilinear forms in Hilbert spaces, see, e.g., [13, §6.2]. Next, we consider finite element methods (FEMs) without and with numerical quadrature. The aim of the paper is to show that well–known error estimates, established for the finite element approximation of elliptic EVPs with classical BCs, hold for the present types of EVPs too. Some attention is also paid to the computational aspects of the resulting algebraic EVP. Finally, the analysis is illustrated by two non-trivial numerical examples, the exact eigenpairs of which can be determined

Schnelle Löser für partielle Differentialgleichungen

[no abstract available

Extraction of Displacement Fields in Heterogeneous Media Using Optimal Local Basis Functions

The Multiscale Spectral Generalized Finite Element Method (MS-GFEM) was developed in recent work by Babuska and Lipton. The method uses optimal local shape functions, optimal in the sense of the Kolmogorov n-width, to approximate solutions to a second order linear elliptic partial differential equation with L-infinity coefficients. In this dissertation an implementation of MS-GFEM over a two subdomain partition of unity is outlined and several numerical experiments are presented. The method is applied to compute local fields inside high contrast particle suspensions. The method\u27s performance is evaluated for various examples with different contrasts between reinforcement particles and matrix material. The numerical experiments are shown to agree with a new theoretical estimate that shows the convergence rate is independent of the elastic properties of particles and matrix materials. A new domain decomposition method based on MS-GFEM is presented. Numerical computations using this iterative method are discussed and the theoretical convergence rate is provided. It is shown that the convergence rate is given by the same near-exponential bound given for MS-GFEM. A systematic method for identifying the worst case load amongst all boundary loads of a fixed energy is introduced. Here the worst case load delivers the largest fraction of input energy into a prescribed subdomain of interest. This leads to an eigenvalue problem, for which the largest eigenvalue is the maximum fraction of energy which concentrates in the subdomain. The associated eigenfunctions are the worst case solutions. These eigenfunctions are related back to the MS-GFEM shape functions and numerical results are presented for several different geometries