Multidomain pseudospectral time‐domain algorithm in curvilinear coordinates system

Abstract

Multidomain pseudospectral time‐domain (MPSTD) algorithm, previously proposed for a Cartesian coordinates system, is extended to Maxwell's equations in a curvilinear coordinates system and applied to electromagnetic wave problems. A novelty patching technique is derived to effectively implement the boundary conditions at the interface between two adjacent subdomains by using characteristic variables. A perfectly matched layer (PML) in a curvilinear coordinates system is constructed according to a well‐posed PML in a Cartesian coordinates system. Numerical results are presented to illustrate the efficiency and accuracy of the method. © 2007 Wiley Periodicals, Inc. Microwave Opt Technol Lett 49: 2618–2624, 2007; Published online in Wiley InterScience (www.interscience.wiley.com) DOI 10.1002/mop.22772

A parallel 3‐D staggered grid pseudospectral time domain method for ground‐penetrating radar wave simulation

“…Thus, the multidomain PSTD method can be an alternative potential for parallel PSTD implementations, especially for models with complicated geometries. However, comparing with the multidomain PSTD method, which often adopts Chebyshev PSTD method [Fan et al, 2002;Zhao and Liu, 2003;Shi and Liang, 2007] and is not suitable to have a big size of subdomain, our parallel PSTD algorithm is more flexible in the selection of subdomain size. Thus, our parallel PSTD algorithm is suitable for simulating large geophysical problems.…”

“…This may be the reason why only 16 grids are adopted as the size of subdomain in most multidomain Chebyshev PSTD simulations [Fan et al, 2002;Zhao and Liu, 2003;Shi and Liang, 2007]. In order to simulate a large problem using multidomain Chebyshev PSTD method, the number of subdomains should be significantly increased.…”

Stability and accuracy of Engquist–Majda absorbing boundary condition for pseudo spectral time domain method

“…This approach is similar to forcing incoming flux to zero as in an uncoupled Riemann problem. 9 The CVs are also used for information exchange (continuity) between the subdomains 10,11 and implementation of boundary conditions such as the outer boundary termination of the PML 12 and perfect electric conductor (PEC). 11 However, any CVs based implementations must be analyzed in detail for stability since the L-CPSTD method, being global in nature, is more sensitive to local treatments of ABCs.…”