A non-intrusive model order reduction approach for parameterized time-domain Maxwell's equations

Abstract

<p style='text-indent:20px;'>We present a non-intrusive model order reduction (NIMOR) approach with an offline-online decoupling for the solution of parameterized time-domain Maxwell's equations. During the offline stage, the training parameters are chosen by using Smolyak sparse grid method with an approximation level <inline-formula><tex-math id="M1">\begin{document}$ L $\end{document}</tex-math></inline-formula> (<inline-formula><tex-math id="M2">\begin{document}$ L\geq1 $\end{document}</tex-math></inline-formula>) over a target parameterized space. For each selected parameter, the snapshot vectors are first produced by a high order discontinuous Galerkin time-domain (DGTD) solver formulated on an unstructured simplicial mesh. In order to minimize the overall computational cost in the offline stage and to improve the accuracy of the NIMOR method, a radial basis function (RBF) interpolation method is then used to construct more snapshot vectors at the sparse grid with approximation level <inline-formula><tex-math id="M3">\begin{document}$ L+1 $\end{document}</tex-math></inline-formula>, which includes the sparse grids from approximation level <inline-formula><tex-math id="M4">\begin{document}$ L $\end{document}</tex-math></inline-formula>. A nested proper orthogonal decomposition (POD) method is employed to extract time- and parameter-independent POD basis functions. By using the singular value decomposition (SVD) method, the principal components of the reduced coefficient matrices of the high-fidelity solutions onto the reduced-order subspace spaned by the POD basis functions are extracted. Moreover, a Gaussian process regression (GPR) method is proposed to approximate the dominating time- and parameter-modes of the reduced coefficient matrices. During the online stage, the reduced-order solutions for new time and parameter values can be rapidly recovered via outputs from the regression models without using the DGTD method. Numerical experiments for the scattering of plane wave by a 2-D dielectric cylinder and a multi-layer heterogeneous medium nicely illustrate the performance of the NIMOR method.</p>