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 sub-domain, a Dirichlet-Neumann 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 multi-point 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 initial-boundary value problem for the time-dependent Maxwell--Schr\"{o}dinger equations in the Coulomb gauge. We first prove the global existence of weak solutions to the equations. Next we propose an energy-conserving 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 time-step 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 time-dependent Maxwell-Schr\"{o}dinger equations under the Lorentz gauge. An alternating Crank-Nicolson 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 H1-norm provided that the exact solution has sufficient regularity, but also numerical results indicate that their convergence rates in both the H1-norm and the L2-norm 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.