A discrete action principle for electrodynamics and the construction of explicit symplectic integrators for linear, non-dispersive media

Jeffrey M. McMahon, Stephen K. Gray, George C Schatz

Research output: Contribution to journalArticle

4 Citations (Scopus)

Abstract

In this work, we derive a discrete action principle for electrodynamics that can be used to construct explicit symplectic integrators for Maxwell's equations. Different integrators are constructed depending on the choice of discrete Lagrangian used to approximate the action. By combining discrete Lagrangians in an explicit symplectic partitioned Runge-Kutta method, an integrator capable of achieving any order of accuracy is obtained. Using the von Neumann stability analysis, we show that the integrators greatly increase the numerical stability and reduce the numerical dispersion compared to other methods. For practical purposes, we demonstrate how to implement the integrators using many features of the finite-difference time-domain method. However, our approach is also applicable to other spatial discretizations, such as those used in finite element methods. Using this implementation, numerical examples are presented that demonstrate the ability of the integrators to efficiently reduce and maintain a minimal amount of numerical dispersion, particularly when the time-step is less than the stability limit. The integrators are therefore advantageous for modeling large, inhomogeneous computational domains.

Original languageEnglish
Pages (from-to)3421-3432
Number of pages12
JournalJournal of Computational Physics
Volume228
Issue number9
DOIs
Publication statusPublished - May 20 2009

Fingerprint

integrators
Electrodynamics
electrodynamics
Runge Kutta methods
Finite difference time domain method
Convergence of numerical methods
Maxwell equations
Finite element method
numerical stability
Runge-Kutta method
finite difference time domain method
Maxwell equation
finite element method

Keywords

  • Dispersion
  • Electrodynamics
  • FDTD
  • Lagrangian
  • Runge-Kutta
  • Stability
  • Symplectic integrator

ASJC Scopus subject areas

  • Computer Science Applications
  • Physics and Astronomy (miscellaneous)

Cite this

A discrete action principle for electrodynamics and the construction of explicit symplectic integrators for linear, non-dispersive media. / McMahon, Jeffrey M.; Gray, Stephen K.; Schatz, George C.

In: Journal of Computational Physics, Vol. 228, No. 9, 20.05.2009, p. 3421-3432.

Research output: Contribution to journalArticle

@article{d520d98e76e8408ea770915369eb5bf1,
title = "A discrete action principle for electrodynamics and the construction of explicit symplectic integrators for linear, non-dispersive media",
abstract = "In this work, we derive a discrete action principle for electrodynamics that can be used to construct explicit symplectic integrators for Maxwell's equations. Different integrators are constructed depending on the choice of discrete Lagrangian used to approximate the action. By combining discrete Lagrangians in an explicit symplectic partitioned Runge-Kutta method, an integrator capable of achieving any order of accuracy is obtained. Using the von Neumann stability analysis, we show that the integrators greatly increase the numerical stability and reduce the numerical dispersion compared to other methods. For practical purposes, we demonstrate how to implement the integrators using many features of the finite-difference time-domain method. However, our approach is also applicable to other spatial discretizations, such as those used in finite element methods. Using this implementation, numerical examples are presented that demonstrate the ability of the integrators to efficiently reduce and maintain a minimal amount of numerical dispersion, particularly when the time-step is less than the stability limit. The integrators are therefore advantageous for modeling large, inhomogeneous computational domains.",
keywords = "Dispersion, Electrodynamics, FDTD, Lagrangian, Runge-Kutta, Stability, Symplectic integrator",
author = "McMahon, {Jeffrey M.} and Gray, {Stephen K.} and Schatz, {George C}",
year = "2009",
month = "5",
day = "20",
doi = "10.1016/j.jcp.2009.01.019",
language = "English",
volume = "228",
pages = "3421--3432",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",
number = "9",

}

TY - JOUR

T1 - A discrete action principle for electrodynamics and the construction of explicit symplectic integrators for linear, non-dispersive media

AU - McMahon, Jeffrey M.

AU - Gray, Stephen K.

AU - Schatz, George C

PY - 2009/5/20

Y1 - 2009/5/20

N2 - In this work, we derive a discrete action principle for electrodynamics that can be used to construct explicit symplectic integrators for Maxwell's equations. Different integrators are constructed depending on the choice of discrete Lagrangian used to approximate the action. By combining discrete Lagrangians in an explicit symplectic partitioned Runge-Kutta method, an integrator capable of achieving any order of accuracy is obtained. Using the von Neumann stability analysis, we show that the integrators greatly increase the numerical stability and reduce the numerical dispersion compared to other methods. For practical purposes, we demonstrate how to implement the integrators using many features of the finite-difference time-domain method. However, our approach is also applicable to other spatial discretizations, such as those used in finite element methods. Using this implementation, numerical examples are presented that demonstrate the ability of the integrators to efficiently reduce and maintain a minimal amount of numerical dispersion, particularly when the time-step is less than the stability limit. The integrators are therefore advantageous for modeling large, inhomogeneous computational domains.

AB - In this work, we derive a discrete action principle for electrodynamics that can be used to construct explicit symplectic integrators for Maxwell's equations. Different integrators are constructed depending on the choice of discrete Lagrangian used to approximate the action. By combining discrete Lagrangians in an explicit symplectic partitioned Runge-Kutta method, an integrator capable of achieving any order of accuracy is obtained. Using the von Neumann stability analysis, we show that the integrators greatly increase the numerical stability and reduce the numerical dispersion compared to other methods. For practical purposes, we demonstrate how to implement the integrators using many features of the finite-difference time-domain method. However, our approach is also applicable to other spatial discretizations, such as those used in finite element methods. Using this implementation, numerical examples are presented that demonstrate the ability of the integrators to efficiently reduce and maintain a minimal amount of numerical dispersion, particularly when the time-step is less than the stability limit. The integrators are therefore advantageous for modeling large, inhomogeneous computational domains.

KW - Dispersion

KW - Electrodynamics

KW - FDTD

KW - Lagrangian

KW - Runge-Kutta

KW - Stability

KW - Symplectic integrator

UR - http://www.scopus.com/inward/record.url?scp=62049084024&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=62049084024&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2009.01.019

DO - 10.1016/j.jcp.2009.01.019

M3 - Article

AN - SCOPUS:62049084024

VL - 228

SP - 3421

EP - 3432

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

IS - 9

ER -