000013660 001__ 13660
000013660 005__ 20161114170126.0
000013660 04107 $$aeng
000013660 046__ $$k2011-05-25
000013660 100__ $$aMotamed, M.
000013660 24500 $$aAnalysis and Computation of Second Order Stochastic Hyperbolic Equations

000013660 24630 $$n3.$$pComputational Methods in Structural Dynamics and Earhquake Engineering
000013660 260__ $$bNational Technical University of Athens, 2011
000013660 506__ $$arestricted
000013660 520__ $$2eng$$aFollowing the works on stochastic elliptic and parabolic problems [1,2], we consider the problem of numerically approximating statistical moments of the solution of linear second order hyperbolic partial differential equations with random coefficients and deterministic initial and boundary data. The coefficients are assumed to depend on a finite number of random variables. In particular, we study the wave field in multi-layered media, where the number of layers determines the number of random variables. We show that the problem is well-posed under suitable assumptions on the data. We propose and analyze a stochastic collocation method to numerically solve the problem. In this method, the problem is first discretized in space and time using a deterministic numerical method, such as the finite element or finite difference method, to obtain a semi-discrete problem. Next, the semi-discrete problem is collocated in the zeros of suitable tensor product orthogonal polynomials (Gauss points) in the probability space. This leads to the solution of uncoupled deterministic problems as in the Monte Carlo approach. 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 different 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, in general, the rate of convergence is only algebraic. An exponential rate of convergence is still possible for 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.

000013660 540__ $$aText je chráněný podle autorského zákona č. 121/2000 Sb.
000013660 653__ $$a

000013660 7112_ $$aCOMPDYN 2011 - 3rd International Thematic Conference$$cIsland of Corfu (GR)$$d2011-05-25 / 2011-05-28$$gCOMPDYN2011
000013660 720__ $$aMotamed, M.$$iNobile, F.$$iTempone, R.
000013660 8560_ $$ffischerc@itam.cas.cz
000013660 8564_ $$s10052$$uhttps://invenio.itam.cas.cz/record/13660/files/454.pdf$$yOriginal version of the author's contribution as presented on CD, section: MS 31 Uncertainty and Reliability in Computational Structural Dynamics.
000013660 962__ $$r13401
000013660 980__ $$aPAPER