Abstract:We present a new higher-order accurate finite difference explicit jump Immersed Interface Method (HEJIIM) for solving two-dimensional elliptic problems with singular source and discontinuous coefficients in the irregular region on a compact Cartesian mesh. We propose a new strategy for discretizing the solution at irregular points on a nine point compact stencil such that the higher-order compactness is maintained throughout the whole computational domain. The scheme is employed to solve four problems embedded with circular and star shaped interfaces in a rectangular region having analytical solutions and varied discontinuities across the interface in source and the coefficient terms. We also simulate a plethora of fluid flow problems past bluff bodies in complex flow situations, which are governed by the Navier-Stokes equations; they include problems involving multiple bodies immersed in the flow as well. In the process, we show the superiority of the proposed strategy over the EJIIM and other existing IIM methods by establishing the rate of convergence and grid independence of the computed solutions. In all the cases our computed results extremely close to the available numerical and experimental results.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: solving two - dimensional elliptic problems with singular source terms and discontinuous coefficients in irregular regions. Specifically, the author proposes a new high - order compact explicit jump interface method (HEJIIM) for solving elliptic partial differential equations with complex boundary conditions on Cartesian grids.
### Problem Background and Challenges
1. **Complexity of Interface Problems**
- Interface problems are of great significance in many applications, such as two - phase flow, multiphase flow, bubble formation, free - surface flow, etc.
- The main difficulty lies in dealing with the discontinuity at the interface. Traditional methods such as body - fitted grids are difficult to generate high - quality grids and are computationally expensive in the case of moving boundaries.
- Cartesian grids are easy to generate, but special treatment is required at the interface to maintain high precision.
2. **Limitations of Existing Methods**
- Peskin's immersed boundary method (IBM) can only achieve first - order accuracy and is limited to the case of continuous solutions.
- Subsequent improved methods such as IIM (Immersed Interface Method), although improving the accuracy, still face challenges when dealing with high - order jump conditions.
- Other methods such as Ghost Fluid Method (GFM) and Matched Interface Boundary (MIB) methods can achieve high - order accuracy in some cases, but usually need to introduce virtual points or complex correction terms.
### Main Contributions of the Paper
1. **Proposing a New Method**
- The author proposes a new high - order compact explicit jump finite - difference immersed interface method (HEJIIM), which can maintain high - order accuracy on a nine - point compact stencil.
- By modifying the explicit jump immersed interface strategy of Wiegmann and Bube, fourth - order accuracy is ensured throughout the computational domain.
2. **Verification and Application**
- This method is used to solve four test problems in rectangular regions with circular and star - shaped interfaces, and its high accuracy and convergence are verified.
- It is further applied to simulate the steady - state flow problems绕过 obstacles, which are controlled by Navier - Stokes equations, including the case of multiple objects immersed in the flow field.
- The results show that this method has a good consistency between the analytical solutions and the experimental results.
### Mathematical Model
Consider the conservative form of a two - dimensional elliptic equation:
\[
(\beta u_x)_x+(\beta u_y)_y+\kappa(x, y)u = f(x, y)+\sigma\delta((x - x^*)(y - y^*)),\quad(x, y)\in\Omega,\quad(x^*, y^*)\in\Gamma
\]
where \(\Gamma\) is the interface embedded in the rectangular region \(\Omega = [x_0, x_f]\times[y_0, y_f]\), and the coefficients \(\beta(x, y)\), \(\kappa(x, y)\) and \(f(x, y)\) may be non - smooth or discontinuous at the interface, causing the solution and its derivatives to also be discontinuous on the interface.
To deal with these discontinuities, two physical jump conditions need to be defined on the interface:
\[
[u]=\lim_{(x,y)\to\Gamma^+}u(x,y)-\lim_{(x,y)\to\Gamma^-}u(x,y)=u^+-u^-=\hat{S}
\]
\[
[\beta u_n]=\lim_{(x,y)\to\Gamma^+}\beta u_n(x,y)-\lim_{(x,y)\to\Gamma^-}\beta u_n(x,y)=\beta^+\frac{\partial u^+}{\partial n}-\beta^-\frac{\partial u^-}{\partial n}=\sigma