.png){kind=icon}
.png){kind=icon}
Return
A FFT accelerated fourth order finite difference method for solving three-dimensional elliptic interface problems
.png){kind=icon}
OA
AI
DOI：10.1016/j.jcp.2023.111924.png)
.png){kind=icon}
.png){kind=icon}
Share.png){kind=icon}
Save
Abstract
En 中文
In this paper, a fourth order augmented matched interface and boundary (AMIB) method is proposed for solving a three-dimensional elliptic interface problem which involves a smooth material interface inside a cuboid domain. On the boundary of the cuboid domain, the fourth order AMIB method can handle different types of boundary conditions, including Dirichlet, Neumann, Robin and their mix combinations in fictitious value generation. Moreover, zero-padding solutions are introduced so that the fast Fourier transform (FFT) algorithm is still valid near the boundary. In dealing with the interior interface, a fourth order ray-casting matched interface and boundary (MIB) scheme is proposed, which enforces the jump conditions along the normal direction for calculating fictitious values. Comparing with the existing MIB scheme, the ray-casting scheme naturally bypasses the corner issue and becomes more robust in handling complex geometry. Based on fictitious values generated near interface and boundary, the fourth order central difference can be corrected at various irregular points including corner points, by introducing Cartesian derivative jumps as auxiliary variables. This gives rise to an enlarged linear system, which can be efficiently solved by the Schur complement procedure together with the FFT inversion of the discrete Laplacian. Extensive numerical experiments have been carried to test the proposed ray-casting AMIB method for numerical accuracy, efficiency, and robustness in corner treatment. The numerical results demonstrate that the ray-casting AMIB scheme not only maintains a fourth order of accuracy in treating various interfaces and boundaries for both solutions and solution gradients, but also attains an overall efficiency on the order of O (n3 log n) for an x n x n uniform grid.(c) 2023 Elsevier Inc. All rights reserved.
Keywords:
Elliptic interface problem
High order central difference schemes
Fast Fourier transform (FFT)
Matched interface and boundary method
(MIB)
Mixed boundary conditions
Gradient recovery
AI Summary
Key information extracted from the uploaded paper, including a brief overview, abstract, background, key highlights, visual analysis, and future outlook.
Journal
IF:
3.8
Papers:
1.5W
Citations:
7.4W
