A stochastic collocation method for the second order wave equation with a discontinuous random speed

In this paper we propose and analyze a stochastic collocation method for solving the second order wave equation with a random wave speed and subjected to deterministic boundary and initial conditions. The speed is piecewise smooth in the physical space and depends on a finite number of random variables. The numerical scheme consists of a finite difference or finite element method in the physical space and a collocation in the zeros of suitable tensor product orthogonal polynomials (Gauss points) in the probability space. This approach leads to the solution of uncoupled deterministic problems as in the Monte Carlo method. We consider both full and sparse tensor product spaces of orthogonal polynomials. We provide a rigorous convergence analysis and demonstrate different types of convergence of the probability error with respect to the number of collocation points for full and sparse tensor product spaces and under some regularity assumptions on the data. In particular, we show that, unlike in elliptic and parabolic problems, the solution to hyperbolic problems is not in general analytic with respect to the random variables. Therefore, the rate of convergence may only be algebraic. An exponential/fast rate of convergence is still possible for some quantities of interest and for the wave solution with particular types of data. We present numerical examples, which confirm the analysis and show that the collocation method is a valid alternative to the more traditional Monte Carlo method for this class of problems.

[1]  Daniel M. Tartakovsky,et al.  Numerical Methods for Differential Equations in Random Domains , 2006, SIAM J. Sci. Comput..

[2]  N. Cutland,et al.  On homogeneous chaos , 1991, Mathematical Proceedings of the Cambridge Philosophical Society.

[3]  Bruno Després,et al.  Uncertainty quantification for systems of conservation laws , 2009, J. Comput. Phys..

[4]  Gianluigi Rozza,et al.  A reduced basis hybrid method for the coupling of parametrized domains represented by fluidic networks , 2012 .

[5]  Christoph Schwab,et al.  Convergence rates for sparse chaos approximations of elliptic problems with stochastic coefficients , 2007 .

[6]  Dongbin Xiu,et al.  The Wiener-Askey Polynomial Chaos for Stochastic Differential Equations , 2002, SIAM J. Sci. Comput..

[7]  R. Tempone,et al.  ON THE OPTIMAL POLYNOMIAL APPROXIMATION OF STOCHASTIC PDES BY GALERKIN AND COLLOCATION METHODS , 2012 .

[8]  R. Ghanem,et al.  Uncertainty propagation using Wiener-Haar expansions , 2004 .

[9]  H. Kreiss,et al.  Some Mathematical and Numerical Questions Connected with First and Second Order Time-Dependent Systems of Partial Differential Equations , 2001, gr-qc/0106085.

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

[11]  Nigel P. Weatherill,et al.  EQSM: An efficient high quality surface grid generation method based on remeshing , 2006 .

[12]  Raul Tempone,et al.  Analysis and implementation issues for the numerical approximation of parabolic equations with random coefficients , 2009 .

[13]  R. Ghanem,et al.  Stochastic Finite Elements: A Spectral Approach , 1990 .

[14]  Houman Owhadi,et al.  Numerical homogenization of the acoustic wave equations with a continuum of scales , 2006 .

[15]  Lars Hr̲mander,et al.  The Analysis of Linear Partial Differential Operators III: Pseudo-Differential Operators , 1985 .

[16]  Fabio Nobile,et al.  A Stochastic Collocation Method for Elliptic Partial Differential Equations with Random Input Data , 2007, SIAM Rev..

[17]  Peter Orbanz.  Probability Theory II , 2011 .

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

[19]  Dongbin Xiu,et al.  Galerkin method for wave equations with uncertain coefficients , 2008 .

[20]  L. Rodino.  Linear Partial Differential Operators in Gevrey Spaces , 1993 .

[21]  Guang Lin,et al.  Predicting shock dynamics in the presence of uncertainties , 2006, J. Comput. Phys..

[22]  J. Hosken,et al.  Ricker wavelets in their various guises , 1988 .

[23]  T. A. Zang,et al.  Spectral Methods: Fundamentals in Single Domains , 2010 .

[24]  G. Karniadakis,et al.  Long-Term Behavior of Polynomial Chaos in Stochastic Flow Simulations , 2006 .

[25]  A. Arias.  A measure of earthquake intensity , 1970 .

[26]  R. Ghanem,et al.  Multi-resolution analysis of wiener-type uncertainty propagation schemes , 2004 .

[27]  Alexandre Ern,et al.  Intrusive projection methods with upwinding for uncertain nonlinear hyperbolic systems , 2009 .

[28]  L. Hörmander.  The analysis of linear partial differential operators , 1990 .

[29]  Dongbin Xiu,et al.  High-Order Collocation Methods for Differential Equations with Random Inputs , 2005, SIAM J. Sci. Comput..

[30]  Julia Charrier,et al.  Strong and Weak Error Estimates for Elliptic Partial Differential Equations with Random Coefficients , 2012, SIAM J. Numer. Anal..

[31]  Marcel Oliver,et al.  Analyticity of Solutions for a Generalized Euler Equation , 1997 .

[32]  I. Babuska,et al.  Solving elliptic boundary value problems with uncertain coefficients by the finite element method: the stochastic formulation , 2005 .

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

[34]  P. Erdös,et al.  Interpolation , 1953, An Introduction to Scientific, Symbolic, and Graphical Computation.

[35]  Heinz-Otto Kreiss,et al.  Difference Approximations for the Second Order Wave Equation , 2002, SIAM J. Numer. Anal..

[36]  Christiaan C. Stolk,et al.  On the modeling and inversion of seismic data , 2000 .

[37]  Ronald L. Wasserstein,et al.  Monte Carlo: Concepts, Algorithms, and Applications , 1997 .

[38]  H. Kreiss,et al.  Initial-Boundary Value Problems and the Navier-Stokes Equations , 2004 .

[39]  Hermann G. Matthies,et al.  Galerkin methods for linear and nonlinear elliptic stochastic partial differential equations , 2005 .

[40]  Erich Novak,et al.  High dimensional polynomial interpolation on sparse grids , 2000, Adv. Comput. Math..

[41]  Tao Tang,et al.  Convergence Analysis for Stochastic Collocation Methods to Scalar Hyperbolic Equations with a Random Wave Speed , 2010 .

[42]  George Em Karniadakis,et al.  Stochastic modeling of random roughness in shock scattering problems: Theory and simulations , 2008 .

[43]  D. Xiu,et al.  Modeling Uncertainty in Steady State Diffusion Problems via Generalized Polynomial Chaos , 2002 .

[44]  K. Brown,et al.  Graduate Texts in Mathematics , 1982 .

[45]  Fabio Nobile,et al.  A Sparse Grid Stochastic Collocation Method for Partial Differential Equations with Random Input Data , 2008, SIAM J. Numer. Anal..

[46]  BabuskaIvo,et al.  A Stochastic Collocation Method for Elliptic Partial Differential Equations with Random Input Data , 2007 .

[47]  Alexandre Ern,et al.  Roe solver with entropy corrector for uncertain hyperbolic systems , 2010, J. Comput. Appl. Math..

[48]  G. M.,et al.  Partial Differential Equations I , 2023, Applied Mathematical Sciences.

[49]  Alexandre Ern,et al.  Intrusive Galerkin methods with upwinding for uncertain nonlinear hyperbolic systems , 2010, J. Comput. Phys..

[50]  Fabio Nobile,et al.  An Anisotropic Sparse Grid Stochastic Collocation Method for Partial Differential Equations with Random Input Data , 2008, SIAM J. Numer. Anal..