An Acceleration Technique for the Augmented IIM for 3D Elliptic Interface Problems

An Acceleration Technique for the Augmented IIM for 3D Elliptic Interface Problems

Year:    2021

Author:    Changjuan Zhang, Zhilin Li, Xingye Yue

Numerical Mathematics: Theory, Methods and Applications, Vol. 14 (2021), Iss. 3 : pp. 773–796

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.

You do not have full access to this article.

Already a Subscriber? Sign in as an individual or via your institution

Journal Article Details

Publisher Name:    Global Science Press

Language:    English

DOI:    https://doi.org/10.4208/nmtma.OA-2020-0112

Numerical Mathematics: Theory, Methods and Applications, Vol. 14 (2021), Iss. 3 : pp. 773–796

Published online:    2021-01

AMS Subject Headings:   

Copyright:    COPYRIGHT: © Global Science Press

Pages:    24

Keywords:    3D elliptic interface problem augmented IIM fast Poisson solver directional least squares interpolation.

Author Details

Changjuan Zhang

Zhilin Li

Xingye Yue

  1. A FFT accelerated fourth order finite difference method for solving three-dimensional elliptic interface problems

    Ren, Yiming

    Zhao, Shan

    Journal of Computational Physics, Vol. 477 (2023), Iss. P.111924

    https://doi.org/10.1016/j.jcp.2023.111924 [Citations: 5]