Abstract
This article presents a finite element method on a fixed mesh for solving a group of inverse geometric problems for recovering the material interface of a linear elasticity system. A partially penalized immersed finite element method is used to discretize both the elasticity interface problems and the objective shape functionals accurately regardless of the shape and location of the interface. Explicit formulas for both the velocity fields and the shape derivatives of IFE shape functions are derived on a fixed mesh and they are employed in the shape sensitivity framework through the discretized adjoint method for accurately and efficiently computing the gradients of objective shape functions with respect to the parameters of the interface curve. The shape optimization for solving an inverse geometric problem is therefore accurately reduced to a constrained optimization that can be implemented efficiently within the IFE framework together with a standard optimization algorithm. We demonstrate features and advantages of the proposed IFE-based shape optimization method by several typical inverse geometric problems for linear elasticity systems.
Original language | English |
---|---|
Article number | 109123 |
Pages (from-to) | 1-24 |
Number of pages | 24 |
Journal | Journal of Computational Physics |
Volume | 404 |
DOIs | |
Publication status | Published - 1 Mar 2020 |
Keywords
- Discontinuous Lamé parameters
- Elasticity systems
- Immersed finite element methods
- Inclusions reconstruction
- Inverse problems
- Shape optimization
ASJC Scopus subject areas
- Numerical Analysis
- Modelling and Simulation
- Physics and Astronomy (miscellaneous)
- General Physics and Astronomy
- Computer Science Applications
- Computational Mathematics
- Applied Mathematics