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

Abstract

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 This work was supported by the King Abdullah University of Science and Technology (AEA project "Bayesian earthquake source validation for ground motion simulation"), the VR project "Effektiva numeriska metoder för stokastiska differentialekvationer med tillämpningar", and the PECOS center at ICES, University of Texas at Austin (Project Number 024550, Center for Predictive Computational Science). The second author was partially supported by the Italian grant FIRB-IDEAS (Project no. RBID08223Z) "Advanced numerical techniques for uncertainty quantification in engineering and life science problems".

A stochastic multiscale method for the elastodynamic wave equation arising from fiber composites

“…In this section, we briefly review the stochastic collocation method for computing the statistical moments of the solution to the micro problem (4), where Y ∈ Γ ⊂ R N is a vector of N independent random variables [5,29,22]. The stochastic collocation method consists of three main steps.…”

“…Two typical examples of sparse grids include total degree and hyperbolic cross sparse grids. We refer to [5,29,22] for more details.…”

“…Recently, it is shown in [22,21] that the solution of stochastic hyperbolic problems, unlike the solution of elliptic and parabolic problems, is not in general analytic with respect to the random variables. The wave solutions may posses high regularity for particular types of smooth data.…”

“…We note that both spatial mesh sizes ∆x 1 and ∆x 2 and time step ∆t are of order h and related by a proper CFL condition. Also Note that the constant in the term O(h r ) depends on Q(u h ) and is uniform with respect to Y , see [22].…”

A wavelet Galerkin approximation of Fredholm integral eigenvalue problems with bidimensional Haar functions

“…This representation is obtained from the eigenvalues and eigenfunctions of a homogeneous Fredholm integral equation of the second kind whose kernel is given by the covariance function of the random process. This approach has been widely used in the parametrization of parameters for elasticity problems, heat and mass transfer, fluid mechanic and acoustic [1,6,9].…”