Consistent modeling of boundaries in acoustic finite-difference time-domain simulations.

The finite-difference time-domain method is one of the most popular for wave propagation in the time domain. One of its advantages is the use of a structured staggered grid, which makes it simple and efficient on modern computer architectures. A drawback, however, is the difficulty in approximating oblique boundaries, having to resort to staircase approximations. In many scattering problems this means that the grid resolution required to obtain an accurate solution is much higher than what is dictated by propagation in a homogeneous material. In this paper zero boundary data are considered, first for the velocity and then the pressure. These two forms of boundary conditions model perfectly rigid and pressure-release boundaries, respectively. A simple and efficient method to consistently model curved rigid boundaries in two dimensions was developed in Tornberg and Engquist [J. Comput. Phys. 227, 6922-6943 (2008)]. Here this treatment is generalized to three dimensions. Based on the approach of this method, a technique to model pressure-release surfaces with second order accuracy and without additional restriction on the timestep is also introduced. The structure of the standard method is preserved, making it easy to use in existing solvers. The effectiveness is demonstrated in several numerical tests.

[1]  Björn Engquist,et al.  A new type of boundary treatment for wave propagation , 2006 .

[2]  R. J. Kruhlak,et al.  Simple conformal methods for finite-difference time-domain modeling of pressure-release surfaces , 1998 .

[3]  Dick Botteldooren,et al.  ACOUSTICAL FINITE-DIFFERENCE TIME-DOMAIN SIMULATION IN A QUASI-CARTESIAN GRID , 1994 .

[4]  B. Engquist,et al.  Numerical approximations of singular source terms in differential equations , 2004 .

[5]  A. Cangellaris,et al.  Analysis of the numerical error caused by the stair-stepped approximation of a conducting boundary in FDTD simulations of electromagnetic phenomena , 1991 .

[6]  Allen Taflove,et al.  Three-dimensional contour FDTD modeling of scattering from single and multiple bodies , 1993 .

[7]  Allen Taflove,et al.  Finite-difference time-domain modeling of curved surfaces (EM scattering) , 1992 .

[8]  R. Sorrentino,et al.  A simple way to model curved metal boundaries in FDTD algorithm avoiding staircase approximation , 1995 .

[9]  Rolf Schuhmann,et al.  Conformal FDTD-methods to avoid time step reduction with and without cell enlargement , 2007, J. Comput. Phys..

[10]  Raj Mittra,et al.  A MODIFIED LOCALLY CONFORMAL FINITE-DIFFERENCE TIME-DOMAIN ALGORITHM FOR MODELING THREE-DIMENSIONAL PERFECTLY CONDUCTING OBJECTS , 1998 .

[11]  Thomas Rylander,et al.  Stable FEM-FDTD hybrid method for Maxwell's equations , 2000 .

[12]  K. Yee.  Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

[13]  R. Mittra,et al.  Time-domain (FE/FDTD) technique for solving complex electromagnetic problems , 1998 .

[14]  J. Hesthaven,et al.  Convergent Cartesian Grid Methods for Maxwell's Equations in Complex Geometries , 2001 .

[15]  Björn Engquist,et al.  Consistent boundary conditions for the Yee scheme , 2008, J. Comput. Phys..

[16]  Ruey-Beei Wu,et al.  Hybrid finite-difference time-domain modeling of curved surfaces using tetrahedral edge elements , 1997 .

[17]  B. Gustafsson.  The convergence rate for difference approximations to mixed initial boundary value problems , 1975 .

[18]  Julius G. Tolan,et al.  Locally conformal method for acoustic finite-difference time-domain modeling of rigid surfaces. , 2003, The Journal of the Acoustical Society of America.

[19]  R. Mittra,et al.  A locally conformal finite-difference time-domain (FDTD) algorithm for modeling three-dimensional perfectly conducting objects , 1997 .

[20]  Heinz-Otto Kreiss,et al.  Difference Approximations for the Second Order Wave Equation , 2002, SIAM J. Numer. Anal..

[21]  D. Botteldooren.  Finite‐difference time‐domain simulation of low‐frequency room acoustic problems , 1995 .

[22]  Modeling sea surface scattering by the time‐domain finite‐difference method , 1996 .

[23]  A finite-difference time-domain solution to scattering from a rough pressure-release surface , 1997 .