×

A new parallelization strategy for solving time-dependent 3D Maxwell equations using a high-order accurate compact implicit scheme. (English) Zbl 1106.78012

Summary: With progress in computer technology there has been renewed interest in a time-dependent approach to solving Maxwell equations. The commonly used Yee algorithm (an explicit central difference scheme for approximation of spatial derivatives coupled with the Leapfrog scheme for approximation of temporal derivatives) yields only a second-order of accuracy. On the other hand, an increasing number of industrial applications, especially in optic and microwave technology, demands high-order accurate numerical modelling. The standard way to increase accuracy of the finite difference scheme without increasing the differential stencil is to replace a 2nd-order accurate explicit scheme for approximation of spatial derivatives with the 4th-order accurate compact implicit scheme. In general, such a replacement requires additional memory resources and slows the computations. However, the curl-based form of Maxwell equations allows us to construct an effective parallel algorithm with the alternating domain decomposition (ADD) minimizing the communication time. We present a new parallel approach to the solution of three-dimensional time-dependent Maxwell equations and provide a theoretical and experimental analysis of its performance.

MSC:

78M20 Finite difference methods applied to problems in optics and electromagnetic theory
65Y05 Parallel numerical computation
65M06 Finite difference methods for initial value and initial-boundary value problems involving PDEs

Software:

ScaLAPACK; MPI
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Yee, IEEE Transactions on Antennas and Propagation 14 pp 302– (1966)
[2] Kashdan, Journal of Scientific Computing 27 pp 75– (2006)
[3] Kashdan, Journal of Computational Physics (2006)
[4] Hesthaven, Computer Modelling in Engineering and Sciences 5 pp 395– (2003)
[5] Spotz, Numerical Methods for Partial Differential Equations 17 pp 657– (2001)
[6] High-order methods. In Advances in Computational Electrodynamics: The Finite-Difference Time-Domain Method, Chapter 2, (ed.). Artech House: Boston, MA, 1998; 63–110.
[7] High-order accurate methods for Maxwell equations. Ph.D. Thesis, Tel Aviv University, 2004.
[8] . Comparisons of staggered and non-staggered schemes for Maxwell’s equations. 12th Annual Review of Progress in Applied Computational Electromagnetics, Monterey, CA, March 1996; 1122–1131.
[9] . Fourth order method for Maxwell’s equations on a staggered mesh. In IEEE Antennas and Propagation Society International Symposium, vol. 4, Montréal, Canada, July 1997; 2156–2159.
[10] Carpenter, Journal of Computational Physics 108 pp 541– (1993)
[11] Fourth order accurate compact implicit method for the Maxwell equations. Ph.D. Thesis, Tel Aviv University, 1998.
[12] . Computational Electrodynamics: The Finite-Difference Time-Domain Method (2nd edn). Artech House: Boston, MA, 2001.
[13] , , , . Progress in cray-based algorithms for computational electromagnetics. In Progress in Electromagnetics Research, vol. 7, Chapter 2, (eds). EMW Publishing: Cambridge, MA, 1993.
[14] Gedney, International Journal on Numerical Modelling: Electronic Networks, Devices and Fields 8 pp 265– (1995)
[15] Varadarajan, IEEE Microwave Guided Wave Letters 4 pp 144– (1994)
[16] . A parallel compact multidimensional numerical algorithm with aeroacoustics applications. ICASE Report No. 99-34, NASA Langley Research Center, 1999.
[17] . A multi-domain tri-diagonal solver. Technical Report TAE No. 753, Technion–Israel Institute of Technology, 1995.
[18] . Boundary and interface conditions for high order finite difference methods applied to the Euler and Navier–Stokes equations. ICASE Report No. 98-19, NASA Langley Research Center, 1998.
[19] , , . Numerical Linear Algebra for High-Performance Computers. SIAM Publishers: Philadelphia, PA, U.S.A., 1998. · Zbl 0914.65014
[20] . Parallel Computers 2: Architecture, Programming and Algorithms. Institute of Physics Publishing: Kent, 1988. · Zbl 0641.68001
[21] , , , . MPI: The Complete Reference. The MIT Press: Cambridge, MA, 1996.
[22] Efficient implementation of a 3-dimensional ADI method on the iPCS/860. In Proceedings of Supercomputing Conference, Portland, OR, U.S.A., 1993; 102–111.
[23] Foster, Journal of Atmospheric and Oceanic Technology 13 pp 1031– (2002)
[24] ScaLAPACK Users’ Guide. SIAM Publishers: Philadelphia, PA, U.S.A., 1997.
[25] , , . Solution of alternating-line processes on modern parallel computers. 28th International Conference on Parallel Processing, Fukushima, Japan, 1999; 208–215.
[26] , , , , , . Tera-scalable algorithms for variable-density hydrodynamics with spectral accuracy. In Proceedings of Supercomputing Conference, Seattle, WA, U.S.A., 2005; 601–617.
[27] . Performance Optimization of Numerically Intensive Codes. SIAM Publishers: Philadelphia, PA, U.S.A., 2001. · Zbl 0985.68006
[28] , . Validation of Dimemas communication model for MPI collective operations. 7th European PVM/MPI Users’ Group Meeting on Recent Advances in Parallel Virtual Machine and Message Passing Interface, 2000; 39–46.
[29] . Overview of Recent Supercomputers. http://www.top500.org/ORSC
[30] , , , , , . LogP: towards a realistic model of parallel computation. 4th ACM SIGPLAN Symposium on Principles and Practice of Parallel Programming, 1993.
[31] Alexandrov, Journal of Parallel and Distributed Computing 44 pp 71– (1997)
[32] Dimemas: performance prediction for message passing applications. http://www.cepba.upc.es/dimemas
[33] Wang, Microwave and Optical Technology Letters 37 pp 8– (2003)
[34] Berenger, Journal of Computational Physics 114 pp 185– (1994) · Zbl 0814.65129
[35] Turkel, Applied Numerical Mathematics 27 pp 533– (1998)
[36] Gedney, IEEE Transactions on Antennas and Propagation 44 pp 1630– (1996)
[37] Time-domain methods for the Maxwell’s equations. Ph.D. Thesis, Royal Institute of Technology, Sweden, 2001.
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.