Back and forth error compensation and correction method for linear hyperbolic systems with application to the Maxwell's equations

Abstract We study the Back and Forth Error Compensation and Correction (BFECC) method for linear hyperbolic systems and in particular for the Maxwell's equations. BFECC has been applied to schemes for scalar advection equations to improve their stability and order of accuracy. Similar results have been established in this paper for linear hyperbolic systems with constant coefficients. We apply BFECC to the central difference scheme, Lax-Friedrichs scheme and a combination of them for the Maxwell's equations and obtain second order accurate schemes with large CFL numbers (greater than 1 in one or two dimensions). The method is further applied to schemes on non-orthogonal unstructured grids. The new BFECC schemes for the Maxwell's equations operate on a single non-staggered grid and are simple to implement on unstructured grids. Numerical examples are given to demonstrate the effectiveness of the new schemes.

[1]  S. Gedney,et al.  Full wave analysis of microwave monolithic circuit devices using a generalized Yee-algorithm based on an unstructured grid , 1996 .

[2]  Ignacio Llamas,et al.  Simulation of bubbles in foam with the volume control method , 2007, SIGGRAPH 2007.

[3]  Ronald Fedkiw,et al.  An Unconditionally Stable MacCormack Method , 2008, J. Sci. Comput..

[4]  Chaowei Hu,et al.  No . 98-32 Weighted Essentially Non-Oscillatory Schemes on Triangular Meshes , 1998 .

[5]  E. Süli,et al.  A convergence analysis of Yee's scheme on nonuniform grids , 1994 .

[6]  Raj Mittra,et al.  Modeling three-dimensional discontinuities in waveguides using nonorthogonal FDTD algorithm , 1992 .

[7]  Ignacio Llamas,et al.  Advections with Significantly Reduced Dissipation and Diffusion , 2007, IEEE Transactions on Visualization and Computer Graphics.

[8]  Chi-Wang Shu,et al.  Locally divergence-free discontinuous Galerkin methods for the Maxwell equations , 2004, Journal of Computational Physics.

[9]  Pingwen Zhang,et al.  Discontinuous Galerkin methods for dispersive and lossy Maxwell's equations and PML boundary conditions , 2004 .

[10]  P. Frederickson,et al.  Higher order solution of the Euler equations on unstructured grids using quadratic reconstruction , 1990 .

[11]  Peter D. Lax,et al.  On the stability of difference approximations to solutions of hyperbolic equations with variable coefficients , 1961 .

[12]  Yingjie Liu,et al.  Back and forth error compensation and correction methods for semi-lagrangian schemes with application to level set interface computations , 2006, Math. Comput..

[13]  C. Makridakis,et al.  Time-discrete finite element schemes for Maxwell's equations , 1995 .

[14]  Moysey Brio,et al.  Overlapping Yee FDTD Method on Nonorthogonal Grids , 2009, J. Sci. Comput..

[15]  D. Komatitsch,et al.  An unsplit convolutional perfectly matched layer improved at grazing incidence for the seismic wave equation , 2007 .

[16]  T. Dupont,et al.  Back and forth error compensation and correction methods for removing errors induced by uneven gradients of the level set function , 2003 .

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

[18]  W. Steen.  Absorption and Scattering of Light by Small Particles , 1999 .

[19]  Pierre Degond,et al.  On a finite-element method for solving the three-dimensional Maxwell equations , 1993 .

[20]  Zhiming Chen,et al.  Finite Element Methods with Matching and Nonmatching Meshes for Maxwell Equations with Discontinuous Coefficients , 2000, SIAM J. Numer. Anal..

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

[22]  J. Hesthaven,et al.  Nodal high-order methods on unstructured grids , 2002 .

[23]  Mengping Zhang,et al.  An analysis of and a comparison between the discontinuous Galerkin and the spectral finite volume methods , 2005 .

[24]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous galerkin finite element method for conservation laws. II: General framework , 1989 .

[25]  Endre Süli,et al.  Error estimates for Yee's method on non-uniform grids , 1994 .

[26]  J. Douglas Faires,et al.  Numerical Analysis , 1981 .

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

[28]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[29]  John B. Schneider,et al.  Understanding the Finite-Difference Time-Domain Method , 2011 .

[30]  Gang Bao,et al.  Finite element approximation of time harmonic waves in periodic structures , 1995 .

[31]  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 .

[32]  Ignacio Llamas,et al.  FlowFixer: Using BFECC for Fluid Simulation , 2005, NPH.