Computation of moments for Maxwell's equations with random interfaces via pivoted low-rank approximation

Abstract The aim of this paper is to compute the mean field and variance of solutions to three-dimensional Maxwell's equations with random interfaces via shape calculus and pivoted low-rank approximation. Based on the perturbation theory and shape calculus, we characterize the statistical moments of solutions to Maxwell's equations with random interfaces in terms of the perturbation magnitude via the first order shape-Taylor expansion. In order to capture oscillations with high resolution close to the interface, an adaptive finite element method using Nedelec's third order edge elements of the first kind is employed to solve the deterministic Maxwell's equations with the mean interface to approximate the expectation of solutions. For the second moment computation, an efficient low-rank approximation of the pivoted Cholesky decomposition is proposed to compute the two-point correlation function to approximate the variance of solutions. Numerical experiments are presented to demonstrate our theoretical results.

[1]  Shaozhong Deng,et al.  An upwinding boundary condition capturing method for Maxwell's equations in media with material interfaces , 2007, J. Comput. Phys..

[2]  Kai Zhang,et al.  Multi-level Monte Carlo weak Galerkin method for elliptic equations with stochastic jump coefficients , 2016, Appl. Math. Comput..

[3]  Christoph Schwab,et al.  Karhunen-Loève approximation of random fields by generalized fast multipole methods , 2006, J. Comput. Phys..

[4]  Stefan Heldmann,et al.  An octree multigrid method for quasi-static Maxwell's equations with highly discontinuous coefficients , 2007, J. Comput. Phys..

[5]  J. Zolésio,et al.  Introduction to shape optimization : shape sensitivity analysis , 1992 .

[6]  Roger Ghanem,et al.  Adaptive polynomial chaos expansions applied to statistics of extremes in nonlinear random vibration , 1998 .

[7]  J. Nédélec.  Mixed finite elements in ℝ3 , 1980 .

[8]  Reinhold Schneider,et al.  Sparse second moment analysis for elliptic problems in stochastic domains , 2008, Numerische Mathematik.

[9]  Gang Bao,et al.  A Robust Numerical Method for the Random Interface Grating Problem via Shape Calculus, Weak Galerkin Method, and Low-Rank Approximation , 2018, J. Sci. Comput..

[10]  Lung-an Ying,et al.  Partial differential equations and the finite element method , 2007, Math. Comput..

[11]  Ralf Hiptmair,et al.  Shape derivatives in differential forms I: an intrinsic perspective , 2013 .

[12]  J. Zou,et al.  An adaptive edge element method and its convergence for a Saddle‐Point problem from magnetostatics , 2012 .

[13]  Kai Zhang,et al.  A weak Galerkin method for diffraction gratings , 2017 .

[14]  Eli Turkel,et al.  A High-Order Accurate Method for Frequency Domain Maxwell Equations with Discontinuous Coefficients , 2006, J. Sci. Comput..

[15]  Ralf Hiptmair,et al.  Convergence analysis of finite element methods for H(curl; Ω)-elliptic interface problems , 2012, Numerische Mathematik.

[16]  Shan Zhao,et al.  High-order FDTD methods via derivative matching for Maxwell's equations with material interfaces , 2004 .

[17]  D. Xiu,et al.  An efficient spectral method for acoustic scattering from rough surfaces , 2007 .

[18]  Liying Zhang,et al.  Preservation of physical properties of stochastic Maxwell equations with additive noise via stochastic multi-symplectic methods , 2014, J. Comput. Phys..

[19]  Raino A. E. Mäkinen,et al.  Introduction to shape optimization - theory, approximation, and computation , 2003, Advances in design and control.

[20]  J. Schöberl,et al.  High order Nédélec elements with local complete sequence properties , 2005 .

[21]  Haijun Wu,et al.  Edge Element Methods for Maxwell's Equations with Strong Convergence for Gauss' Laws , 2014, SIAM J. Numer. Anal..

[22]  Daniel M. Tartakovsky,et al.  Stochastic analysis of transport in tubes with rough walls , 2006, J. Comput. Phys..

[23]  Raúl Tempone,et al.  Galerkin Finite Element Approximations of Stochastic Elliptic Partial Differential Equations , 2004, SIAM J. Numer. Anal..

[24]  Chang-Yeol Jung,et al.  Maxwell solutions in media with multiple random interfaces , 2014 .

[25]  J. Nédélec.  A new family of mixed finite elements in ℝ3 , 1986 .

[26]  Weiying Zheng,et al.  An Adaptive FEM for a Maxwell Interface Problem , 2016, J. Sci. Comput..

[27]  Liying Zhang,et al.  A stochastic multi-symplectic scheme for stochastic Maxwell equations with additive noise , 2014, J. Comput. Phys..

[28]  George E. Karniadakis,et al.  Adaptive Generalized Polynomial Chaos for Nonlinear Random Oscillators , 2005, SIAM J. Sci. Comput..

[29]  Peter Benner,et al.  Uncertainty Quantification for Maxwell's Equations Using Stochastic Collocation and Model Order Reduction , 2013 .

[30]  Jingzhi Li,et al.  Optimal shape for a nozzle design problem using an arbitraryLagrangian–Eulerian finite element method , 2014 .

[31]  Peter Monk,et al.  Finite Element Methods for Maxwell's Equations , 2003 .

[32]  Zhiming Chen,et al.  Finite Element Methods with Matching and Nonmatching Meshes for Maxwell Equations with Discontinuous Coefficients , 2000, SIAM J. Numer. Anal..

[33]  Helmut Harbrecht,et al.  First order second moment analysis for stochastic interface problems based on low-rank approximation , 2013 .

[34]  Kai Zhang,et al.  Multi-level Monte Carlo weak Galerkin method with nested meshes for stochastic Brinkman problem , 2018, J. Comput. Appl. Math..

[35]  D. Xiu,et al.  Modeling uncertainty in flow simulations via generalized polynomial chaos , 2003 .

[36]  H. Harbrecht,et al.  On the low-rank approximation by the pivoted Cholesky decomposition , 2012 .

[37]  Claudio Canuto,et al.  A fictitious domain approach to the numerical solution of PDEs in stochastic domains , 2007, Numerische Mathematik.

[38]  Shaozhong Deng,et al.  On the immersed interface method for solving time-domain Maxwell's equations in materials with curved dielectric interfaces , 2008, Comput. Phys. Commun..

[39]  Peter Monk,et al.  A finite element method for approximating the time-harmonic Maxwell equations , 1992 .

[40]  M. Delfour,et al.  Shapes and Geometries: Analysis, Differential Calculus, and Optimization , 1987 .