A High Performance Computing and Sensitivity Analysis Algorithm for Stochastic Many-Particle Wave Scattering

We present a high performance computing framework for quantifying uncertainties in the propagation of acoustic waves through a stochastic media comprising a large number of three-dimensional particles. We subsequently describe an efficient postprocessing approach using our framework to statistically quantify the sensitivity of the uncertainties with respect to input parameters that govern the stochasticity in the model. The stochasticity arises through the random positions and orientations of the component particles in the media. Simulation even for a single deterministic three-dimensional configuration is inherently difficult because of the large number of particles; the stochasticity leads to a larger dimensional model involving three spatial variables and additional stochastic variables, and accounting for uncertainty in key parameters of the input probability distributions leads to prohibitive computational complexity. In the first part of our paper we describe a high performance computing framework f...

[1]  R. García-Pelayo,et al.  Multiple Scattering , 2001 .

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

[3]  Ivan G. Graham,et al.  A high-order algorithm for obstacle scattering in three dimensions , 2004 .

[4]  Thomas Wriedt,et al.  Light scattering by single erythrocyte: Comparison of different methods , 2006 .

[5]  Josef Dick,et al.  Walsh Spaces Containing Smooth Functions and Quasi-Monte Carlo Rules of Arbitrary High Order , 2008, SIAM J. Numer. Anal..

[6]  Omar M. Knio,et al.  Spectral Methods for Uncertainty Quantification , 2010 .

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

[8]  J. Hesthaven,et al.  Adaptive sparse grid algorithms with applications to electromagnetic scattering under uncertainty , 2011 .

[9]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[10]  Mustafa Kuzuoglu,et al.  A coordinate transformation approach for efficient repeated solution of Helmholtz equation pertaining to obstacle scattering by shape deformations , 2014, Comput. Phys. Commun..

[11]  Elisabeth Ullmann,et al.  Further analysis of multilevel Monte Carlo methods for elliptic PDEs with random coefficients , 2012, Numerische Mathematik.

[12]  Tony Warnock,et al.  An Algorithm for Uniform Random Sampling of Points In and On a Hypersphere , 1985, Inf. Process. Lett..

[13]  Dongbin Xiu,et al.  A Stochastic Collocation Algorithm with Multifidelity Models , 2014, SIAM J. Sci. Comput..

[14]  Benjamin Stamm,et al.  A reduced basis method for electromagnetic scattering by multiple particles in three dimensions , 2012, J. Comput. Phys..

[15]  Stuart C. Hawkins,et al.  A high-order algorithm for multiple electromagnetic scattering in three dimensions , 2009, Numerical Algorithms.

[16]  Jukka Sarvas,et al.  UNIFIED DERIVATION OF THE TRANSLATIONAL ADDITION THEOREMS FOR THE SPHERICAL SCALAR AND VECTOR WAVE FUNCTIONS , 2008 .

[17]  Ralf Hiptmair,et al.  Convergence analysis with parameter estimates for a reduced basis acoustic scattering T-matrix method , 2012 .

[18]  A. Lacis,et al.  Multiple Scattering of Light by Particles: Radiative Transfer and Coherent Backscattering , 2006 .

[19]  Josef Dick,et al.  Explicit Constructions of Quasi-Monte Carlo Rules for the Numerical Integration of High-Dimensional Periodic Functions , 2007, SIAM J. Numer. Anal..

[20]  A. Ditkowski,et al.  Wave scattering by randomly shaped objects , 2012 .

[21]  J. Conoir,et al.  Multiple Scattering: Interaction of Time-Harmonic Waves with N Obstacles , 2006 .

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

[23]  Mahadevan Ganesh,et al.  An efficient $$\mathcal {O}(N)$$O(N) algorithm for computing $$\mathcal {O}(N^2)$$O(N2) acoustic wave interactions in large $$N$$N-obstacle three dimensional configurations , 2014 .

[24]  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.

[25]  Max Gunzburger,et al.  A Multilevel Stochastic Collocation Method for Partial Differential Equations with Random Input Data , 2014, SIAM/ASA J. Uncertain. Quantification.

[26]  M. Mishchenko Electromagnetic Scattering by Particles and Particle Groups: An Introduction , 2014 .

[27]  Thomas Wriedt,et al.  Simulation of light scattering by biconcave Cassini ovals using the nullfield method with discrete sources , 2005 .

[28]  Claude Jeffrey Gittelson,et al.  Sparse tensor discretizations of high-dimensional parametric and stochastic PDEs* , 2011, Acta Numerica.