
We present a new fixed mesh algorithm for solving a class of interface
inverse problems for the typical elliptic interface problems. These interface
inverse problems are formulated as shape optimization prob lems whose
objective functionals depend on the shape of the interface. Regardless of the
location of the interface, both the governing partial differential equations
and the objective functional are discretized optimally, with respect to the
involved polynomial space, by an immersed finite element (IFE) method on a
fixed mesh. Furthermore, the formula for the gradient of the descritized
objective function is de rived within the IFE framework that can be computed
accurately and efficiently through the discretized adjoint procedure. Features
of this proposed IFE method based on a fixed mesh are demonstrated by its
applications to three representative interface inverse problems: the interface
inverse problem with an internal measurement on a subdomain, a
DirichletNeumann type inverse problem whose data is given on the boundary, and
a heat dissipation design problem.

We construct and analyze a group of immersed finite element (IFE) spaces
formed by linear, bilinear and rotated Q1 polynomials for solving planar
elasticity equation involving interface. The shape functions in these IFE
spaces are constructed through a group of approximate jump conditions such that
the unisolvence of the bilinear and rotated Q1 IFE shape functions are always
guaranteed regardless of the Lam\`e parameters and the interface location. The
boundedness property and a group of identities of the proposed IFE shape
functions are established. A multipoint Taylor expansion is utilized to show
the optimal approximation capabilities for the proposed IFE spaces through the
Lagrange type interpolation operators.

In this paper, we consider the initialboundary value problem for the
timedependent MaxwellSchr\"{o}dinger equations in the Coulomb gauge. We
first prove the global existence of weak solutions to the equations. Next we
propose an energyconserving fully discrete finite element scheme for the
system and prove the existence and uniqueness of solutions to the discrete
system. The optimal error estimates for the numerical scheme without any
timestep restrictions are then derived. Numerical results are provided to
support our theoretical analysis.

In this paper we study the numerical method and the convergence for solving
the timedependent MaxwellSchr\"{o}dinger equations under the Lorentz gauge.
An alternating CrankNicolson finite element method for solving the problem is
presented and the optimal error estimate for the numerical algorithm is
obtained by a mathematical inductive method. Numerical examples are then
carried out to confirm the theoretical results.

This article presents new immersed finite element (IFE) methods for solving
the popular second order elliptic interface problems on structured Cartesian
meshes even if the involved interfaces have nontrivial geometries. These IFE
methods contain extra stabilization terms introduced only at interface edges
for penalizing the discontinuity in IFE functions. With the enhanced stability
due to the added penalty, not only these IFE methods can be proven to have the
optimal convergence rate in the H1norm provided that the exact solution has
sufficient regularity, but also numerical results indicate that their
convergence rates in both the H1norm and the L2norm do not deteriorate when
the mesh becomes finer which is a shortcoming of the classic IFE methods in
some situations. Trace inequalities are established for both linear and
bilinear IFE functions that are not only critical for the error analysis of
these new IFE methods, but also are of a great potential to be useful in error
analysis for other IFE methods.