Abstract: | A new fast algorithm based on the augmented immersed interface method
and a fast Poisson solver is proposed to solve three dimensional elliptic interface
problems with a piecewise constant but discontinuous coefficient. In the new approach, an augmented variable along the interface, often the jump in the normal
derivative along the interface is introduced so that a fast Poisson solver can be utilized. Thus, the solution of the Poisson equation depends on the augmented variable
which should be chosen such that the original flux jump condition is satisfied. The
discretization of the flux jump condition is done by a weighted least squares interpolation using the solution at the grid points, the jump conditions, and the governing
PDEs in a neighborhood of control points on the interface. The interpolation scheme
is the key to the success of the augmented IIM particularly. In this paper, the key
new idea is to select interpolation points along the normal direction in line with the
flux jump condition. Numerical experiments show that the method maintains second order accuracy of the solution and can reduce the CPU time by 20-50%. The
number of the GMRES iterations is independent of the mesh size. |