TY - JOUR
T1 - Domain uncertainty quantification in computational electromagnetics
AU - Aylwin, Ruben
AU - Jerez-Hanckes, Carlos
AU - Schwab, Christoph
AU - Zech, Jakob
N1 - Funding Information:
This work was supported in part by Fondecyt Regular 1171491, Conicyt-PFCHA/Doctorado Nacional/2017-21171791, SAM ETH Zurich-Puc Open Seed Fund 201603, and the Swiss National Science Foundation under grant SNF149819. The fourth author was supported by the Swiss National Science Foundation under Early Postdoc.Mobility Fellowship 184530.
Funding Information:
∗Received by the editors January 18, 2019; accepted for publication (in revised form) December 6, 2019; published electronically February 19, 2020. https://doi.org/10.1137/19M1239374 Funding: This work was supported in part by Fondecyt Regular 1171491, Conicyt-PFCHA/Doctorado Nacional/2017-21171791, SAM ETH Zurich–Puc Open Seed Fund 201603, and the Swiss National Science Foundation under grant SNF149819. The fourth author was supported by the Swiss National Science Foundation under Early Postdoc.Mobility Fellowship 184530. †School of Engineering, Pontificia Universidad Católica de Chile, 7820436 Santiago, Chile (rdaylwin@uc.cl). ‡Faculty of Engineering and Sciences, Universidad Adolfo Ibañez, 7941169 Santiago, Chile (carlos.jerez@uai.cl). §Seminar for Applied Mathematics, ETH Zürich, 8092 Zürich, Switzerland (schwab@math.ethz.ch). ¶Department of Aeronautics and Astronautics, Massachusetts Institute of Technology, Cambridge, MA 02139 (jzech@mit.edu).
Publisher Copyright:
© 2020 Society for Industrial and Applied Mathematics and American Statistical Association.
PY - 2020
Y1 - 2020
N2 - We study the numerical approximation of time-harmonic, electromagnetic fields inside a lossy cavity of uncertain geometry. Key assumptions are a possibly high-dimensional parametrization of the uncertain geometry along with a suitable transformation to a fixed, nominal domain. This uncertainty parametrization results in families of countably parametric, Maxwell-like cavity problems that are posed in a single domain, with inhomogeneous coefficients that possess finite, possibly low spatial regularity, but exhibit holomorphic parametric dependence in the differential operator. Our computational scheme is composed of a sparse grid interpolation in the high-dimensional parameter domain and an Hcurl -conforming edge element discretization of the parametric problem in the nominal domain. As a stepping-stone in the analysis, we derive a novel Strang-type lemma for Maxwell-like problems in the nominal domain, which is of independent interest. Moreover, we accommodate arbitrary small Sobolev regularity of the electric field and also cover uncertain isotropic constitutive or material laws. The shape holomorphy and edge-element consistency error analysis for the nominal problem are shown to imply convergence rates for multilevel Monte Carlo and for quasi-Monte Carlo integration, as well as sparse grid approximations, in uncertainty quantification for computational electromagnetics. They also imply expression rate estimates for deep ReLU networks of shape-to-solution maps in this setting. Finally, our computational experiments confirm the presented theoretical results.
AB - We study the numerical approximation of time-harmonic, electromagnetic fields inside a lossy cavity of uncertain geometry. Key assumptions are a possibly high-dimensional parametrization of the uncertain geometry along with a suitable transformation to a fixed, nominal domain. This uncertainty parametrization results in families of countably parametric, Maxwell-like cavity problems that are posed in a single domain, with inhomogeneous coefficients that possess finite, possibly low spatial regularity, but exhibit holomorphic parametric dependence in the differential operator. Our computational scheme is composed of a sparse grid interpolation in the high-dimensional parameter domain and an Hcurl -conforming edge element discretization of the parametric problem in the nominal domain. As a stepping-stone in the analysis, we derive a novel Strang-type lemma for Maxwell-like problems in the nominal domain, which is of independent interest. Moreover, we accommodate arbitrary small Sobolev regularity of the electric field and also cover uncertain isotropic constitutive or material laws. The shape holomorphy and edge-element consistency error analysis for the nominal problem are shown to imply convergence rates for multilevel Monte Carlo and for quasi-Monte Carlo integration, as well as sparse grid approximations, in uncertainty quantification for computational electromagnetics. They also imply expression rate estimates for deep ReLU networks of shape-to-solution maps in this setting. Finally, our computational experiments confirm the presented theoretical results.
KW - Bayesian inverse problems
KW - Computational electromagnetics
KW - Finite elements
KW - Shape holomorphy
KW - Sparse grid quadrature
KW - Uncertainty quantification
UR - http://www.scopus.com/inward/record.url?scp=85085339980&partnerID=8YFLogxK
U2 - 10.1137/19M1239374
DO - 10.1137/19M1239374
M3 - Article
AN - SCOPUS:85085339980
SN - 2166-2525
VL - 8
SP - 301
EP - 341
JO - SIAM-ASA Journal on Uncertainty Quantification
JF - SIAM-ASA Journal on Uncertainty Quantification
IS - 1
ER -