A novel surface-derivative-free of jumps AIIM with triangulated surfaces for 3D Helmholtz interface problems

Abstract

Triangular surface-based 3D IIM (Immersed Interface Method) algorithms face major challenges due to the need to calculate surface derivative of jumps. This paper proposes a fast, easy-to-implement, surface-derivative-free of jumps, augmented IIM (AIIM) with triangulated surfaces for 3D Helmholtz interface problems for the first time, which combines the simplified AIIM with domain decomposed and embedding techniques. The computational domain is divided into sub-domains along the interface and the solutions of sub-domains are continuously extended into larger regular domains by embedding. The jumps in normal derivative of solution along the interfaces in the extended domains are introduced as unknowns to impose the original jump relations. The original problem is simplified into Helmholtz interface problems with constant coefficients by coupling them with the augmented equation, which is then solved using fast simplified AIIM. This approach eliminates the need to compute surface derivatives of jumps, making implementation of 3D IIM based on triangulated surfaces fairly simple. Numerical results demonstrate that the algorithm is efficient and can achieve the overall second-order accuracy.