A Fourier penalty method for solving the time-dependent Maxwell's equations in domains with curved boundaries

DOI10.1016/J.JCP.2015.11.031zbMATH Open1351.78058OpenAlexW1520113090MaRDI QIDQ729324FDOQ729324

Authors Ryan Galagusz, D. Shirokoff, Jean-Christophe Nave 

Publication date 20 December 2016

Published in Journal of Computational Physics (Search for Journal in Brave)

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

zbMATH Keywords

Maxwell equationsFourier continuationFourier methodsactive penalty methodsharp mask function

Mathematics Subject Classification ID

Electromagnetic theory (general) (78A25) Spectral, collocation and related methods for initial value and initial-boundary value problems involving PDEs (65M70) Spectral, collocation and related methods applied to problems in optics and electromagnetic theory (78M22)

Abstract: We present a high order, Fourier penalty method for the Maxwell's equations in the vicinity of perfect electric conductor boundary conditions. The approach relies on extending the smooth non-periodic domain of the equations to a periodic domain by removing the exact boundary conditions and introducing an analytic forcing term in the extended domain. The forcing, or penalty term is chosen to systematically enforce the boundary conditions to high order in the penalty parameter, which then allows for higher order numerical methods. We present an efficient numerical method for constructing the penalty term, and discretize the resulting equations using a Fourier spectral method. We demonstrate convergence orders of up to 3.5 for the one-dimensional Maxwell's equations, and show that the numerical method does not suffer from dispersion (or pollution) errors. We also illustrate the approach in two dimensions and demonstrate convergence orders of 2.5 for transverse magnetic modes and 1.5 for the transverse electric modes. We conclude the paper with numerous test cases in dimensions two and three including waves traveling in an irregular waveguide, and scattering off of a windmill-like geometry.

Recommendations

• Numerical investigation of a boundary penalization method for Maxwell equations
• Interior penalty method for the indefinite time-harmonic Maxwell equations
• A Legendre pseudospectral penalty scheme for solving time-domain Maxwell's equations
• Characteristic-based, time-dependent Maxwell equation solvers on a general curvilinear frame
• On the Finite Element Approximation of a 2D Maxwell Eigenvalue Problem in a Domain with Curved Boundaries
• On the immersed interface method for solving time-domain Maxwell's equations in materials with curved dielectric interfaces
• A Finite-Volume Method for the Maxwell Equations in the Time Domain
• Numerical solution of Maxwell's equations in axisymmetric domains with the Fourier singular complement method
• A Time-Domain Finite Element Method for Maxwell's Equations
• Finite-Element Domain Approximation for Maxwell Variational Problems on Curved Domains

Cites Work

• scientific article; zbMATH DE number 2212009 (Why is no real title available?)
• A Fourier spectral method for the Navier-Stokes equations with volume penalization for moving solid obstacles
• A PseudoSpectral Procedure for the Solution of Nonlinear Wave Equations with Examples from Free-Surface Flows
• A characteristic based volume penalization method for general evolution problems applied to compressible viscous flows
• A fictitious domain model for the Stokes/Brinkman problem with jump embedded boundary conditions
• A mathematical analysis of the PML method
• A perfectly matched layer for the absorption of electromagnetic waves
• A sharp-interface active penalty method for the incompressible Navier-Stokes equations
• A stable and conservative interface treatment of arbitrary spatial accuracy
• A volume penalization method for incompressible flows and scalar advection-diffusion with moving obstacles
• An optimal penalty method for a hyperbolic system modeling the edge plasma transport in a tokamak
• Computation of turbulent flow past an array of cylinders using a spectral method with Brinkman penalization
• Efficient and stable perfectly matched layer for CEM
• Eulerian–Lagrangian grid coupling and penalty methods for the simulation of multiphase flows interacting with complex objects
• High-order unconditionally stable FC-AD solvers for general smooth domains. I: Basic elements
• High-order unconditionally stable FC-AD solvers for general smooth domains. II: Elliptic, parabolic and hyperbolic PDEs; theoretical considerations
• Long time behavior of the perfectly matched layer equations in computational electromagnetics
• On the Long-Time Behavior of Unsplit Perfectly Matched Layers
• Penalization for non-linear hyperbolic system
• Revisiting and extending interface penalties for multi-domain summation-by-parts operators
• Simulation of confined magnetohydrodynamic flows with Dirichlet boundary conditions using a pseudo-spectral method with volume penalization
• Spectral Methods for Time-Dependent Problems

Cited In (5)

• Title not available (Why is no real title available?)
• FDTD schemes for Maxwell's equations with embedded perfect electric conductors based on the correction function method
• A Legendre pseudospectral penalty scheme for solving time-domain Maxwell's equations
• Numerical investigation of a boundary penalization method for Maxwell equations
• Bounded error schemes for the wave equation on complex domains

Uses Software

• Eigtool
• DistMesh

This page was built for publication: A Fourier penalty method for solving the time-dependent Maxwell's equations in domains with curved boundaries

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