Numerical Analysis and Scientific Computing Seminar
High Order Accurate Compact Finite Difference Scheme for Wave Equation in General Shaped Domains using Simple Meshes
Speaker: Eli Turkel, Tel Aviv University
Location: Warren Weaver Hall 1302
Date: Nov. 2, 2018, 10 a.m.
We consider the wave equation in a general shaped domain perhaps including general interfaces. We consider two fourth order accurate algorithms that use only simple (Cartesian or Polar) grids. Both approaches rely on the difference potential approach.
In one approach we discretize implicitly in time with fourth order accuracy which leads to a modified Helmholtz equation. (2D or 3D with a variable speed of sound) for the upper time level. This modified Helmholtz equation is then approximated by a fourth order, in space, scheme. Difference potentials are then used to solve for the general shaped domain using Cartesian grids.
The other approach employs the difference potentials directly in 3D+time and again the shape of the boundary does not have to conform to the grid, on which the difference potentials are computed. A fourth order discretization is used in both time and space. Furthermore, it takes advantage of Huygens' principle. F\So for long simulation times, the method offers sub-linear complexity, i.e., appears asymptotically cheaper than a typical explicit scheme.
Results are presented for both variations.