4.7 Article

POD-based model order reduction with an adaptive snapshot selection for a discontinuous Galerkin approximation of the time-domain Maxwell's equations

Journal

JOURNAL OF COMPUTATIONAL PHYSICS
Volume 396, Issue -, Pages 106-128

Publisher

ACADEMIC PRESS INC ELSEVIER SCIENCE
DOI: 10.1016/j.jcp.2019.05.051

Keywords

Time-domain Maxwell's equations; Discontinuous Galerkin method; Model order reduction; Proper orthogonal decomposition; Adaptive snapshot selection algorithm

Funding

  1. NSFC [61772003]
  2. Science Strength Promotion Programme of UESTC

Ask authors/readers for more resources

In this work we report on a reduced-order model (ROM) for the system of time-domain Maxwell's equations discretized by a discontinuous Galerkin (DG) method. We leverage previous results on proper orthogonal decomposition (POD) [1,2], in particular for the wave equation [3], to propose a POD-based ROM with an adaptive snapshot selection strategy where the snapshots are produced by a high order discontinuous Galerkin time-domain (DGTD) solver. The latter is formulated on an unstructured simplicial mesh, and combines a centered scheme for the definition of the numerical fluxes of the electric and magnetic fields at element interfaces with a second order leap-frog (LF2) time scheme for the time integration of the associated semi-discrete equations. The POD-based ROM is established by projecting (Galerkin projection) the global semi-discrete DG scheme onto a low-dimensional space generated by the POD basis vectors. Inspired from the approach followed in [2,3], we derive error bounds for the POD-based ROM that is adapted to our particular modeling and discretization settings. The adaptive snapshot selection algorithm exploits the results of this analysis to measure the control error. A snapshot choosing rule aiming at keeping the error estimate close to a target selection error tolerance is proposed, which is similar to the standard rules found in adaptive time-stepping ordinary differential equations (ODEs) solvers. An incremental singular value decomposition (ISVD) algorithm is used to update the SVD on-the-fly when a new snapshot is available. The purpose of this adaptive selection strategy is to save memory without storing snapshots, while producing a smaller error. Numerical experiments for the 2-D time-domain Maxwell's equations nicely illustrate the performance of the resulting POD-based ROM with adaptive snapshot selection. (C) 2019 Elsevier Inc. All rights reserved.

Authors

I am an author on this paper
Click your name to claim this paper and add it to your profile.

Reviews

Primary Rating

4.7
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

No Data Available
No Data Available