Uniform stable conformal convolutional perfectly matched layer for enlarged cell technique conformal finite-difference time-domain method*

Based on conformal construction of physical model in a three-dimensional Cartesian grid, an integral-based conformal convolutional perfectly matched layer (CPML) is given for solving the truncation problem of the open port when the enlarged cell technique conformal finite-difference time-domain (ECT-CFDTD) method is used to simulate the wave propagation inside a perfect electric conductor (PEC) waveguide. The algorithm has the same numerical stability as the ECT-CFDTD method. For the long-time propagation problems of an evanescent wave in a waveguide, several numerical simulations are performed to analyze the reflection error by sweeping the constitutive parameters of the integral-based conformal CPML. Our numerical results show that the integral-based conformal CPML can be used to efficiently truncate the open port of the waveguide.

[1]  J. Bérenger.  Perfectly matched layer for the FDTD solution of wave-structure interaction problems , 1996 .

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

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

[4]  Raj Mittra,et al.  Frequency dependence of the constitutive parameters of causal perfectly matched anisotropic absorbers , 1996 .

[5]  J. Bérenger,et al.  Application of the CFS PML to the absorption of evanescent waves in waveguides , 2002, IEEE Microwave and Wireless Components Letters.

[6]  Q.H. Liu,et al.  Enlarged cells for the conformal FDTD method to avoid the time step reduction , 2004, IEEE Microwave and Wireless Components Letters.

[7]  Yongdong Li,et al.  UNIPIC code for simulations of high power microwave devices , 2009 .

[8]  Jianguo Wang,et al.  Truncation of open boundaries of cylindrical waveguides in 2.5-dimensional problems by using the convolutional perfectly matched layer , 2006, IEEE Transactions on Plasma Science.

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

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

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

[12]  王玥 Wang Yue,et al.  Technique of complex geometry modeling and grid generation for fully electromagnetic particle simulation , 2011 .

[13]  Jean-Pierre Berenger,et al.  Improved PML for the FDTD solution of wave-structure interaction problems , 1997 .

[14]  Jean-Pierre Berenger,et al.  An effective PML for the absorption of evanescent waves in waveguides , 1998 .

[15]  Yongdong Li,et al.  Three-dimensional parallel UNIPIC-3D code for simulations of high-power microwave devices , 2010 .

[16]  Jean-Pierre Berenger,et al.  Evanescent waves in PML's: origin of the numerical reflection in wave-structure interaction problems , 1999 .

[17]  Juan Chen,et al.  A New Method to Avoid Reduction of Time Step in the CP-FDTD Method , 2007, IEEE Transactions on Antennas and Propagation.