Nonconforming Immersed Finite Element Methods for Interface Problems

Files

TR Number

Date

2013-05-04

Journal Title

Journal ISSN

Volume Title

Publisher

Virginia Tech

Abstract

In science and engineering, many simulations are carried out over domains consisting of multiple materials separated by curves/surfaces. If partial differential equations (PDEs) are used to model these simulations, it usually leads to the so-called interface problems of PDEs whose coefficients are discontinuous. In this dissertation, we consider nonconforming immersed "nite element (IFE) methods and error analysis for interface problems.

We "first consider the second order elliptic interface problem with a discontinuous diffusion coefficient. We propose new IFE spaces based on the nonconforming rotated Q1 "finite elements on Cartesian meshes. The degrees of freedom of these IFE spaces are determined by midpoint values or average integral values on edges. We investigate fundamental properties of these IFE spaces, such as unisolvency and partition of unity, and extend well-known trace inequalities and inverse inequalities to these IFE functions. Through interpolation error analysis, we prove that these IFE spaces have optimal approximation capabilities.

We use these IFE spaces to develop partially penalized Galerkin (PPG) IFE schemes whose bilinear forms contain penalty terms over interface edges. Error estimation is carried out for these IFE schemes. We prove that the PPG schemes with IFE spaces based on integral-value degrees of freedom have the optimal convergence in an energy norm. Following a similar approach, we prove that the interior penalty discontinuous Galerkin schemes based on these IFE functions also have the optimal convergence. However, for the PPG schemes based on midpoint-value degrees of freedom, we prove that they have at least a sub-optimal convergence. Numerical experiments are provided to demonstrate features of these IFE methods and compare them with other related numerical schemes.

We extend nonconforming IFE schemes to the planar elasticity interface problem with discontinuous Lam"e parameters. Vector-valued nonconforming rotated Q1 IFE functions with integral-value degrees of freedom are unisolvent with appropriate interface jump conditions. More importantly, the Galerkin IFE scheme using these vector-valued nonconforming rotated Q1 IFE functions are "locking-free" for nearly incompressible elastic materials.

In the last part of this dissertation, we consider potential applications of IFE methods to time dependent PDEs with moving interfaces. Using IFE functions in the discretization in space enables the applicability of the method of lines. Crank-Nicolson type fully discrete schemes are also developed as alternative approaches for solving moving interface problems.

Description

Keywords

Immersed Finite Element, Elliptic Interface Problems, Cartesian Mesh, Nonconforming Rotated Q1 Finite Element, Error Analysis, E

Citation