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

DOI10.1016/J.CPC.2015.11.009MaRDI QIDQ2374022zbMATH OpenOpenAlexFDO

Authors H. Vincenti, Jean-Luc Vay

Publication date 14 December 2016

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

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

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)

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.

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

Cites work

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

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

Describes a project that uses

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)