Scalable Preconditioners for Structure Preserving Discretizations of Maxwell Equations in First Order Form

Phillips, Edward Geoffrey; Shadid, John N.; Cyr, Eric C.

doi:10.1137/17M1135827

Title: Scalable Preconditioners for Structure Preserving Discretizations of Maxwell Equations in First Order Form

Full Record
Other Related Research

Abstract

Here, we report multiple physical time-scales can arise in electromagnetic simulations when dissipative effects are introduced through boundary conditions, when currents follow external time-scales, and when material parameters vary spatially. In such scenarios, the time-scales of interest may be much slower than the fastest time-scales supported by the Maxwell equations, therefore making implicit time integration an efficient approach. The use of implicit temporal discretizations results in linear systems in which fast time-scales, which severely constrain the stability of an explicit method, can manifest as so-called stiff modes. This study proposes a new block preconditioner for structure preserving (also termed physics compatible) discretizations of the Maxwell equations in first order form. The intent of the preconditioner is to enable the efficient solution of multiple-time-scale Maxwell type systems. An additional benefit of the developed preconditioner is that it requires only a traditional multigrid method for its subsolves and compares well against alternative approaches that rely on specialized edge-based multigrid routines that may not be readily available. Lastly, results demonstrate parallel scalability at large electromagnetic wave CFL numbers on a variety of test problems.

Authors:

Phillips, Edward Geoffrey ^[1]; Shadid, John N. ^[1]; Cyr, Eric C. ^[1]

Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)

Publication Date:: Tue May 01 00:00:00 EDT 2018

Research Org.:: Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)

Sponsoring Org.:: USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR); USDOE National Nuclear Security Administration (NNSA)

OSTI Identifier:: 1441484

Report Number(s):: SAND-2018-4332J
Journal ID: ISSN 1064-8275; 662668

Grant/Contract Number:: AC04-94AL85000; NA0003525

Resource Type:: Accepted Manuscript

Journal Name:: SIAM Journal on Scientific Computing

Additional Journal Information:: Journal Volume: 40; Journal Issue: 3; Journal ID: ISSN 1064-8275

Publisher:: SIAM

Country of Publication:: United States

Language:: English

Subject:: 97 MATHEMATICS AND COMPUTING; Maxwell equations; preconditioning; compatible discretization

Citation Formats


                    Phillips, Edward Geoffrey, Shadid, John N., and Cyr, Eric C. Scalable Preconditioners for Structure Preserving Discretizations of Maxwell Equations in First Order Form.  United States: N. p., 2018. 
Web.  doi:10.1137/17M1135827.

Copy to clipboard


                    Phillips, Edward Geoffrey, Shadid, John N., & Cyr, Eric C. Scalable Preconditioners for Structure Preserving Discretizations of Maxwell Equations in First Order Form.  United States.  https://doi.org/10.1137/17M1135827

Copy to clipboard


                    Phillips, Edward Geoffrey, Shadid, John N., and Cyr, Eric C. Tue .  
"Scalable Preconditioners for Structure Preserving Discretizations of Maxwell Equations in First Order Form".  United States.  https://doi.org/10.1137/17M1135827.  https://www.osti.gov/servlets/purl/1441484.

Copy to clipboard


                    
@article{osti_1441484,

  title        = {Scalable Preconditioners for Structure Preserving Discretizations of Maxwell Equations in First Order Form},

  author       = {Phillips, Edward Geoffrey and Shadid, John N. and Cyr, Eric C.},

  abstractNote = {Here, we report multiple physical time-scales can arise in electromagnetic simulations when dissipative effects are introduced through boundary conditions, when currents follow external time-scales, and when material parameters vary spatially. In such scenarios, the time-scales of interest may be much slower than the fastest time-scales supported by the Maxwell equations, therefore making implicit time integration an efficient approach. The use of implicit temporal discretizations results in linear systems in which fast time-scales, which severely constrain the stability of an explicit method, can manifest as so-called stiff modes. This study proposes a new block preconditioner for structure preserving (also termed physics compatible) discretizations of the Maxwell equations in first order form. The intent of the preconditioner is to enable the efficient solution of multiple-time-scale Maxwell type systems. An additional benefit of the developed preconditioner is that it requires only a traditional multigrid method for its subsolves and compares well against alternative approaches that rely on specialized edge-based multigrid routines that may not be readily available. Lastly, results demonstrate parallel scalability at large electromagnetic wave CFL numbers on a variety of test problems.},

  doi          = {10.1137/17M1135827},

  journal      = {SIAM Journal on Scientific Computing},

  number       = 3,

  volume       = 40,

  place        = {United States},

  year         = {Tue May 01 00:00:00 EDT 2018},

  month        = {Tue May 01 00:00:00 EDT 2018}

}

Copy to clipboard

Journal Article:

Free Publicly Available Full Text

Accepted Manuscript (DOE)

Publisher's Version of Record

https://doi.org/10.1137/17M1135827

Other availability

Search WorldCat to find libraries that may hold this journal

Citation Metrics:

Cited by: 4 works

Citation information provided by
Web of Science

Save / Share:

Export Metadata

Save to My Library

Similar Records in DOE PAGES and OSTI.GOV collections:

On the Preconditioning of a Newton-Krylov Solver for a High-Order reconstructed Discontinuous Galerkin Discretization of All-Speed Compressible Flow with Phase Change for Application in Laser-Based Additive Manufacturing

Thesis/Dissertation Weston, Brian T.

This dissertation focuses on the development of a fully-implicit, high-order compressible ow solver with phase change. The work is motivated by laser-induced phase change applications, particularly by the need to develop large-scale multi-physics simulations of the selective laser melting (SLM) process in metal additive manufacturing (3D printing). Simulations of the SLM process require precise tracking of multi-material solid-liquid-gas interfaces, due to laser-induced melting/ solidi cation and evaporation/condensation of metal powder in an ambient gas. These rapid density variations and phase change processes tightly couple the governing equations, requiring a fully compressible framework to robustly capture the rapid density variations ofmore »« less
https://doi.org/10.2172/1361587

Full Text Available
Preconditioning a Newton-Krylov solver for all-speed melt pool flow physics

Journal Article Weston, Brian ; Nourgaliev, Robert ; Delplanque, Jean -Pierre ; ... - Journal of Computational Physics

n this paper, we introduce a multigrid block-based preconditioner for solving linear systems arising from a Discontinuous Galerkin discretization of the all-speed Navier-Stokes equations with phase change. The equations are discretized in conservative form with a reconstructed Discontinuous Galerkin (rDG) method and integrated with fully-implicit time discretization schemes. To robustly converge the numerically stiff systems, we use the Newton-Krylov framework with a primitive-variable formulation (pressure, velocity, and temperature), which is better conditioned than the conservative-variable form at low-Mach number. In the limit of large acoustic CFL number and viscous Fourier number, there is a strong coupling between the velocity-pressure systemmore »« less
Cited by 3
https://doi.org/10.1016/j.jcp.2019.07.045

Full Text Available
Scalability of preconditioners as a strategy for parallel computation of compressible fluid flow

Technical Report Hansen, G A

Parallel implementations of a Newton-Krylov-Schwarz algorithm are used to solve a model problem representing low Mach number compressible fluid flow over a backward-facing step. The Mach number is specifically selected to result in a numerically {open_quote}stiff{close_quotes} matrix problem, based on an implicit finite volume discretization of the compressible 2D Navier-Stokes/energy equations using primitive variables. Newton`s method is used to linearize the discrete system, and a preconditioned Krylov projection technique is used to solve the resulting linear system. Domain decomposition enables the development of a global preconditioner via the parallel construction of contributions derived from subdomains. Formation of the global preconditionermore »« less
https://doi.org/10.2172/237435

Full Text Available
Efficient discontinuous Galerkin implementations and preconditioners for implicit unsteady compressible flow simulations

Journal Article Franciolini, Matteo ; Fidkowski, Krzysztof J. ; Crivellini, Andrea - Computers and Fluids

This work presents and compares efficient implementations of high-order discontinuous Galerkin methods: a modal matrix-free discontinuous Galerkin (DG) method, a hybridizable discontinuous Galerkin (HDG) method, and a primal formulation of HDG, applied to the implicit solution of unsteady compressible flows. The matrix-free implementation allows for a reduction of the memory footprint of the solver when dealing with implicit time-accurate discretizations. HDG reduces the number of globally-coupled degrees of freedom relative to DG, at high order, by statically condensing element-interior degrees of freedom from the system in favor of face unknowns. The primal formulation further reduces the element-interior degrees of freedommore »« less
Cited by 8
https://doi.org/10.1016/j.compfluid.2020.104542

Full Text Available
A semi-implicit, energy- and charge-conserving particle-in-cell algorithm for the relativistic Vlasov-Maxwell equations

Journal Article Chen, Guangye ; Chacon, Luis ; Yin, Lin ; ... - Journal of Computational Physics

Conventional explicit electromagnetic particle-in-cell (PIC) algorithms do not conserve discrete energy exactly. Time-centered fully implicit PIC algorithms can conserve discrete energy exactly, but may introduce large dispersion errors in the light-wave modes. This can lead to intolerable simulation errors where accurate light propagation is needed (e.g. in laser-plasma interactions). In this study, we selectively combine the leap-frog and Crank-Nicolson methods to produce an exactly energy- and charge-conserving relativistic electromagnetic PIC algorithm. Specifically, we employ the leap-frog method for Maxwell's equations, and the Crank-Nicolson method for the particle equations. The semi-implicit formulation still features a timestep CFL, but facilitates exact globalmore »« less
https://doi.org/10.1016/j.jcp.2020.109228

Full Text Available

Similar Records