Detailed analysis of the effects of stencil spatial variations with arbitrary high-order finite-difference Maxwell solver

DOI10.1016/J.CPC.2015.11.009zbMATH Open1351.78057arXiv1507.05572OpenAlexW1652551217MaRDI QIDQ2374022FDOQ2374022

Authors: H. Vincenti, Jean-Luc Vay 

Publication date: 14 December 2016

Published in: Computer Physics Communications (Search for Journal in Brave)

Abstract: Due to discretization effects and truncation to finite domains, many electromagnetic simulations present non-physical modifications of Maxwell's equations in space that may generate spurious signals affecting the overall accuracy of the result. Such modifications for instance occur when Perfectly Matched Layers (PMLs) are used at simulation domain boundaries to simulate open media. Another example is the use of arbitrary order Maxwell solver with domain decomposition technique that may under some condition involve stencil truncations at subdomain boundaries, resulting in small spurious errors that do eventually build up. In each case, a careful evaluation of the characteristics and magnitude of the errors resulting from these approximations, and their impact at any frequency and angle, requires detailed analytical and numerical studies. To this end, we present a general analytical approach that enables the evaluation of numerical discretization errors of fully three-dimensional arbitrary order finite-difference Maxwell solver, with arbitrary modification of the local stencil in the simulation domain. The analytical model is validated against simulations of domain decomposition technique and PMLs, when these are used with very high-order Maxwell solver, as well as in the infinite order limit of pseudo-spectral solvers. Results confirm that the new analytical approach enables exact predictions in each case. It also confirms that the domain decomposition technique can be used with very high-order Maxwell solver and a reasonably low number of guard cells with negligible effects on the whole accuracy of the simulation.

Full work available at URL: https://arxiv.org/abs/1507.05572

Recommendations

• Analysis of spatial high-order finite difference methods for Maxwell's equations in dispersive media
• Efficiency of the perfectly matched layer with high-order finite difference and pseudo-spectral Maxwell solvers
• Choice of the perfectly matched layer boundary condition for frequency-domain Maxwell's equations solvers
• Optimized three-dimensional FDTD discretizations of Maxwell's equations on Cartesian grids
• κ-Space Transfer Function Design of Discrete Operators: Application To Maxwell's Time-Domain Equations

zbMATH Keywords

perfectly matched layersdomain decomposition technique3D electromagnetic simulationseffects of stencil truncation errorspseudo-spectral Maxwell solververy high-order Maxwell solver

Mathematics Subject Classification ID

Finite difference methods applied to problems in optics and electromagnetic theory (78M20) Maxwell equations (35Q61) Multigrid methods; domain decomposition for initial value and initial-boundary value problems involving PDEs (65M55)

Cites Work

• A perfectly matched layer for the absorption of electromagnetic waves
• Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media
• Rational Chebyshev spectral methods for unbounded solutions on an infinite interval using polynomial-growth special basis functions
• Asymmetric perfectly matched layer for the absorption of waves.
• Efficiency of the perfectly matched layer with high-order finite difference and pseudo-spectral Maxwell solvers
• A domain decomposition method for pseudo-spectral electromagnetic simulations of plasmas
• High-Order Finite Differences and the Pseudospectral Method on Staggered Grids
• Numerical stability analysis of the pseudo-spectral analytical time-domain PIC algorithm
• On the elimination of numerical Čerenkov radiation in PIC simulations
• Closed-form expressions for the finite difference approximations of first and higher derivatives based on Taylor series
• Mathematical proof of closed form expressions for finite difference approximations based on Taylor series

Cited In (6)

• Two-step perfectly matched layer for arbitrary-order pseudo-spectral analytical time-domain methods
• X-dispersionless Maxwell solver for plasma-based particle acceleration
• A generalized massively parallel ultra-high order FFT-based Maxwell solver
• A dual grid geometric electromagnetic particle in cell method
• \textsc{Smilei} : a collaborative, open-source, multi-purpose particle-in-cell code for plasma simulation
• A hybrid nodal-staggered pseudo-spectral electromagnetic particle-in-cell method with finite-order centering

Uses Software

• Mathematica

This page was built for publication: Detailed analysis of the effects of stencil spatial variations with arbitrary high-order finite-difference Maxwell solver

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2374022)