Simulation of Maxwell's Equations on GPU Using a High-Order Error-Minimized Scheme. (8th March 2017)

• Record Type: Journal Article
• Authors: Sheu, Tony W. H.; Wang, S. Z.; Li, J. H.; Smith, Matthew R.
• Is Part Of: Communications in computational physics. Volume 21:Number 4(2017:Apr.)
• Page Start: 1039
• Page End: 1064
• Publication Date: 2017-03-08
• Subjects: 35Q61, -- 65M06, -- 65Y05; Hamiltonian systems, -- symplecticity, -- GPU, -- 3-D Maxwell's equations, -- dispersion relation equation; Mathematical physics -- Data processing -- Periodicals; Physics -- Data processing -- Periodicals; 530.150285
• DOI: 10.4208/cicp.OA-2016-0079
• Languages: English
• ISSNs: 1815-2406

Abstract

Abstract: In this study an explicit Finite Difference Method (FDM) based scheme is developed to solve the Maxwell's equations in time domain for a lossless medium. This manuscript focuses on two unique aspects – the three dimensional time-accurate discretization of the hyperbolic system of Maxwell equations in three-point non-staggered grid stencil and it's application to parallel computing through the use of Graphics Processing Units (GPU). The proposed temporal scheme is symplectic, thus permitting conservation of all Hamiltonians in the Maxwell equation. Moreover, to enable accurate predictions over large time frames, a phase velocity preserving scheme is developed for treatment of the spatial derivative terms. As a result, the chosen time increment and grid spacing can be optimally coupled. An additional theoretical investigation into this pairing is also shown. Finally, the application of the proposed scheme to parallel computing using one Nvidia K20 Tesla GPU card is demonstrated. For the benchmarks performed, the parallel speedup when compared to a single core of an Intel i7-4820K CPU is approximately 190x.