Quasi-Monte Carlo finite element analysis for wave propagation in heterogeneous random media

We propose and analyze a quasi-Monte Carlo (QMC) algorithm for efficient simulation of wave propagation modeled by the Helmholtz equation in a bounded region in which the refractive index is random and spatially heterogenous. Our focus is on the case in which the region can contain multiple wavelengths. We bypass the usual sign-indefiniteness of the Helmholtz problem by switching to an alternative sign-definite formulation recently developed by Ganesh and Morgenstern (Numerical Algorithms, 83, 1441-1487, 2020). The price to pay is that the regularity analysis required for QMC methods becomes much more technical. Nevertheless we obtain a complete analysis with error comprising stochastic dimension truncation error, finite element error and cubature error, with results comparable to those obtained for the diffusion problem.

[1]  Jared Wunsch,et al.  Sharp High-Frequency Estimates for the Helmholtz Equation and Applications to Boundary Integral Equations , 2015, SIAM J. Math. Anal..

[2]  I. Sloan.  Lattice Methods for Multiple Integration , 1994 .

[3]  Jeffrey Galkowski,et al.  Optimal constants in nontrapping resolvent estimates and applications in numerical analysis , 2018, Pure and Applied Analysis.

[4]  R. Gantner.  Dimension Truncation in QMC for Affine-Parametric Operator Equations , 2016 .

[5]  Mahadevan Ganesh,et al.  High-order FEM-BEM computer models for wave propagation in unbounded and heterogeneous media: Application to time-harmonic acoustic horn problem , 2016, J. Comput. Appl. Math..

[6]  Frances Y. Kuo,et al.  Application of Quasi-Monte Carlo Methods to Elliptic PDEs with Random Diffusion Coefficients: A Survey of Analysis and Implementation , 2016, Foundations of Computational Mathematics.

[7]  F. Ihlenburg.  Finite Element Analysis of Acoustic Scattering , 1998 .

[8]  Junshan Lin,et al.  An Efficient Numerical Method for Acoustic Wave Scattering in Random Media , 2014, SIAM/ASA J. Uncertain. Quantification.

[9]  Frances Y. Kuo,et al.  High-dimensional integration: The quasi-Monte Carlo way*† , 2013, Acta Numerica.

[10]  HELMUT HARBRECHT,et al.  On the quasi-Monte Carlo method with Halton points for elliptic PDEs with log-normal diffusion , 2016, Math. Comput..

[11]  Frances Y. Kuo,et al.  Multi-level quasi-Monte Carlo finite element methods for a class of elliptic partial differential equations with random coefficients , 2012, 1208.6349.

[12]  M. Ganesh,et al.  A coercive heterogeneous media Helmholtz model: formulation, wavenumber-explicit analysis, and preconditioned high-order FEM , 2019, Numerical Algorithms.

[13]  Ivan G. Graham,et al.  Stability and finite element error analysis for the Helmholtz equation with variable coefficients , 2019, Math. Comput..

[14]  Francisco-Javier Sayas,et al.  An overlapping decomposition framework for wave propagation in heterogeneous and unbounded media: Formulation, analysis, algorithm, and simulation , 2019, J. Comput. Phys..

[15]  Josef Dick,et al.  Multi-level higher order QMC Galerkin discretization for affine parametric operator equations , 2014, 1406.4432.

[16]  O. R. Pembery,et al.  The Helmholtz equation in heterogeneous media: A priori bounds, well-posedness, and resonances , 2018, Journal of Differential Equations.

[17]  Frances Y. Kuo,et al.  Circulant embedding with QMC: analysis for elliptic PDE with lognormal coefficients , 2018, Numerische Mathematik.

[18]  R. Kress,et al.  Inverse Acoustic and Electromagnetic Scattering Theory , 1992 .

[19]  Frances Y. Kuo,et al.  Analysis of quasi-Monte Carlo methods for elliptic eigenvalue problems with stochastic coefficients , 2018, Numerische Mathematik.

[20]  Théophile Chaumont-Frelet,et al.  On high order methods for the heterogeneous Helmholtz equation , 2016, Comput. Math. Appl..

[21]  W. McLean.  Strongly Elliptic Systems and Boundary Integral Equations , 2000 .

[22]  Harald Niederreiter,et al.  Random number generation and Quasi-Monte Carlo methods , 1992, CBMS-NSF regional conference series in applied mathematics.

[23]  O. R. Pembery,et al.  The Helmholtz Equation in Random Media: Well-Posedness and A Priori Bounds , 2018, SIAM/ASA J. Uncertain. Quantification.

[24]  Frances Y. Kuo,et al.  Higher Order QMC Petrov-Galerkin Discretization for Affine Parametric Operator Equations with Random Field Inputs , 2014, SIAM J. Numer. Anal..

[25]  Mahadevan Ganesh,et al.  A Sign-Definite Preconditioned High-Order FEM, Part I: Formulation and Simulation for Bounded Homogeneous Media Wave Propagation , 2017, SIAM J. Sci. Comput..

[26]  Frances Y. Kuo,et al.  Higher order QMC Galerkin discretization for parametric operator equations , 2013, 1309.4624.

[27]  C. Schwab,et al.  Quasi-Monte Carlo integration for affine-parametric, elliptic PDEs: local supports imply product weights , 2016 .

[28]  Hélène Barucq,et al.  Stability analysis of heterogeneous Helmholtz problems and finite element solution based on propagation media approximation , 2016, Math. Comput..

[29]  S. Nicaise,et al.  Wavenumber explicit convergence analysis for finite element discretizations of general wave propagation problems , 2019, IMA Journal of Numerical Analysis.

[30]  Josef Dick,et al.  Higher Order Quasi-Monte Carlo Integration for Holomorphic, Parametric Operator Equations , 2014, SIAM/ASA J. Uncertain. Quantification.

[31]  J. Nédélec.  Acoustic and electromagnetic equations , 2001 .

[32]  R. Caflisch.  Monte Carlo and quasi-Monte Carlo methods , 1998, Acta Numerica.