An FDTD scheme on a face-centered-cubic (FCC) grid for the solution of the wave equation

A method is proposed to improve the numerical dispersion characteristics for simulations of the scalar wave equation in 3D using the FDTD method. The improvements are realized by choosing a face-centered-cubic (FCC) grid instead of the typical Cartesian (Yee) grid, which exhibits non-physical distortions of the wavefront due to the FD stencil. FCC grids are the logical extension of hexagonal grids in 2D, and have been shown previously to provide optimal sampling of space based on close packing of spheres (highest density). The difference equations are developed for the wave equation on this alternative grid, and the dispersion relationship and stability for grids of equal and non-equal aspect ratios are derived. A comparison is made between FCC and Cartesian formulations, based upon having an equal volume density of gridpoints in each method (i.e. the computational storage requirements of each method would be the same for the same simulated space). The comparison shows that the FCC grid exhibits a much more isotropic dispersion relation than the Cartesian grid of equivalent density. Furthermore, for an equivalent density, the FCC method has a more relaxed stability criterion by a factor of approximately 1.35, resulting in a further reduction in computational resources.

[1]  John B. Schneider,et al.  An FDTD method for analysis of scattering from rough fluid-fluid interfaces , 2001 .

[2]  Masahiro Sato Application of the finite-difference time-domain method and total-field/scattered field formulation to scattering phenomena in solids , 2007 .

[3]  Jun Zhang,et al.  A fourth-order compact difference scheme on face centered cubic grids with multigrid method for solving 2D convection diffusion equation , 2003, Math. Comput. Simul..

[4]  Gary Cohen Higher-Order Numerical Methods for Transient Wave Equations , 2001 .

[5]  Waymond R. Scott,et al.  Elastic waves interacting with buried land mines: a study using the FDTD method , 2002, IEEE Trans. Geosci. Remote. Sens..

[6]  David Middleton,et al.  Sampling and Reconstruction of Wave-Number-Limited Functions in N-Dimensional Euclidean Spaces , 1962, Inf. Control..

[7]  Jinquan Xu,et al.  Multidimensional Finite Differencing (MDFD) with Hypersphere-Close-Pack Grids , 2005 .

[8]  B. Finkelstein,et al.  The spectral order of accuracy: A new unified tool in the design methodology of excitation-adaptive wave equation FDTD schemes , 2009, J. Comput. Phys..

[9]  Alfredo Bautista,et al.  An explicit finite differences scheme over hexagonal tessellation , 2001, Appl. Math. Lett..

[10]  Tang Xiaohong,et al.  The construction of low-dispersive FDTD on hexagon , 2005 .

[11]  Yen Liu,et al.  Fourier Analysis of Numerical Algorithms for the Maxwell Equations , 1993 .

[12]  Mikko Karttunen,et al.  Stencils with isotropic discretization error for differential operators , 2006 .

[13]  Maarten van Walstijn,et al.  Wideband and Isotropic Room Acoustics Simulation Using 2-D Interpolated FDTD Schemes , 2010, IEEE Transactions on Audio, Speech, and Language Processing.

[14]  Allen Taflove,et al.  Computational Electrodynamics the Finite-Difference Time-Domain Method , 1995 .

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