The Poisson equation is central in numerous physics and engineering applications, such as computational fluid dynamics and acoustic wave propagation, where efficient and accurate solutions are essential. This study focuses on the numerical solution of the 2D Poisson equation with Dirichlet boundary conditions using a fourth-order compact Implicit Finite Difference scheme. Finite difference methods, particularly high-order schemes, are advantageous for solving the Poisson equation due to their efficiency and suitability for structured grids. To address the computational demands of large-scale problems, we incorporate domain decomposition and the Multicolor Successive Over Relaxation method, facilitating parallel computation. Through numerical experiments, we demonstrate that our approach significantly enhances both accuracy and computational efficiency when compared to traditional second-order methods.
Abstract The Poisson equation is central in numerous physics and engineering applications, such as computational fluid dynamics and acoustic wave propagation, where efficient and accurate [...]
This work forms the foundation for addressing high-order immersed interface methods to solve interface problems and enables us to conduct in-depth examination of this theory. Here, we focus on the introduction a fourth-order finite-difference formulation to approximate the second-order derivative of discontinuous functions. The approach is based on the combination of a high-order implicit formulation and the immersed interface method. The idea is to modify the standard schemes by introducing additional contribution terms based on jump conditions. These contributions are calculated only at grid points where the stencil intersects with the interface. Here, we discuss the issues of implementing the one-dimensional Poisson equation and the heat conduction equation with discontinuous solutions as a three-point stencil for each grid point on the computational domain. In both cases, the resulting discretization approach yields a tridiagonal linear system with matrix coefficients identical to those employed for smooth solutions. We present several numerical experiments to verify the feasibility and accuracy of the method. Thus, this high-order method provides an attractive numerical framework that can efficiently lead to the solution to more complex problems.
Abstract This work forms the foundation for addressing high-order immersed interface methods to solve interface problems and enables us to conduct in-depth examination of this theory. [...]
This paper introduces a sixth-order Immersed Interface Method (IIM) for addressing 2D Poisson problems characterized by a discontinuous forcing function with straight interfaces. In the presence of this discontinuity, the problem exhibits a non-smooth solution at the interface that divides the domain into two regions. Here, the IIM is employed to compute the solution on a fixed Cartesian grid. This method integrates necessary jump conditions resulting from the interface into the numerical schemes. In order to achieve a sixth-order method, the proposed approach combines implicit finite differences with the IIM. The proposed scheme is efficient because the matrix arising from discretization remains the same as in the smooth problem, and changes are made to the resulting linear system by introducing new terms on the right side. These supplementary terms account for the discontinuities in the solution and its derivatives, with calculations restricted near the interface. The paper demonstrates the accuracy of the proposed method through various numerical examples.
Abstract This paper introduces a sixth-order Immersed Interface Method (IIM) for addressing 2D Poisson problems characterized by a discontinuous forcing function with straight interfaces. [...]
Multicolor Parallel Fourth-Order Implicit Finite Difference for Solving the 2D Poisson Equation 
