High order two-grid finite difference methods for interface and internal layer problems

by   Zhilin Li, et al.

Second order accurate Cartesian grid methods have been well developed for interface problems in the literature. However, it is challenging to develop third or higher order accurate methods for problems with curved interfaces and internal boundaries. In this paper, alternative approaches based on two different grids are developed for some interface and internal layer problems, which are different from adaptive mesh refinement (AMR) techniques. For one dimensional, or two-dimensional problems with straight interfaces or boundary layers that are parallel to one of the axes, the discussion is relatively easy. One of challenges is how to construct a fourth order compact finite difference scheme at boarder grid points that connect two meshes. A two-grid method that employs a second order discretization near the interface in the fine mesh and a fourth order discretization away from the interface in the coarse and boarder grid points is proposed. For two dimensional problems with a curved interface or an internal layer, a level set representation is utilized for which we can build a fine mesh within a tube |φ( x) | ≤δ h of the interface. A new super-third seven-point discretization that can guarantee the discrete maximum principle has been developed at hanging nodes. The coefficient matrices of the finite difference equations developed in this paper are M-matrices, which leads to the convergence of the finite difference schemes. Non-trivial numerical examples presented in this paper have confirmed the desired accuracy and convergence.


page 1

page 2

page 3

page 4


Compact 9-Point Finite Difference Methods with High Accuracy Order and/or M-Matrix Property for Elliptic Cross-Interface Problems

In this paper we develop finite difference schemes for elliptic problems...

Elastic wave propagation in curvilinear coordinates with mesh refinement interfaces by a fourth order finite difference method

We develop a fourth order accurate finite difference method for the thre...

Geometric multigrid method for solving Poisson's equation on octree grids with irregular boundaries

A method is presented to include irregular domain boundaries in a geomet...

A ghost-point based second order accurate finite difference method on uniform orthogonal grids for electromagnetic scattering around PEC

We propose a finite difference method to solve Maxwell's equations in ti...

New Numerical Interface Scheme for the Kurganov-Tadmor second-order Method

In this paper, we develop a numerical scheme to handle interfaces across...

A SBP-SAT FDTD Subgridding Method Using Staggered Yee's Grids Without Modifying Field Components

A summation-by-parts simultaneous approximation term (SBP-SAT) finite-di...

Alternating direction implicit time integrations for finite difference acoustic wave propagation: Parallelization and convergence

This work studies the parallelization and empirical convergence of two f...

Please sign up or login with your details

Forgot password? Click here to reset