Analysis and Computation of Hyperbolic PDEs with Random Data
|
|
- Lisa Richardson
- 5 years ago
- Views:
Transcription
1 Analysis and Computation of Hyperbolic PDEs with Random Data Mohammad Motamed 1, Fabio Nobile 2,3 and Raúl Tempone 1 1 King Abdullah University of Science and Technology, Thuwal, Saudi Arabia 2 EPFL Lausanne, Switzerland 3 Politecnico di Milano, Milan, Italy March 20, Introduction Hyperbolic partial differential equations (PDEs) are mathematical models of wave phenomena, with applications in a wide range of scientific and engineering fields such as electromagnetic radiation, geosciences, fluid and solid mechanics, aeroacoustics, and general relativity. The theory of hyperbolic problems, including Friedrichs and Kreiss theories, has been well developed based on energy estimates and the method of Fourier and Laplace transforms [8, 16]. Moreover, stable numerical methods, such as the finite difference method [14], the finite volume method [17], the finite element method [6], spectral methods [4], and the boundary element method [11], have been proposed to compute approximate solutions of hyperbolic problems. However, the development of the theory and numerics for hyperbolic PDEs has been based on the assumption that all input data, such as coefficients, initial data, boundary and force terms, and computational domain, are exactly known. There is an increasing interest in including uncertainty in these models and quantifying its effects on the predicted solution or other quantities of physical interest. The uncertainty may be due to either an intrinsic variability of the physical system (aleatoric uncertainty), or our ignorance or inability to accurately characterize all input data (epistemic uncertainty). For example, in earthquake modeling, seismic waves propagate in a geological region where, due to soil spatial variability and the uncertainty of measured soil parameters, both kinds of uncertainties are present. Consequently, the field of uncertainty quantification (UQ) has arisen as a new scientific discipline. UQ is the science of quantitative characterization, reduction, and propagation of uncertainties and the key for achieving validated predictive computations. The numerical solution of hyperbolic problems with random inputs is a new branch of UQ, relying on a broad range of mathematics and statistics groundwork with associated algorithmic and computational development, and aiming at accurate and fast propagation of uncertainties and the quantification of uncertain simulation-based predictions. The most popular method for solving PDEs in probabilistic setting is the Monte Carlo sampling, see for instance [7]. It consists in generating independent realizations drawn from the input mohammad.motamed@kaust.edu.sa fabio.nobile@epfl.ch raul.tempone@kaust.edu.sa 1
2 distribution and then computing sample statistics of the corresponding output values. This allows one to reuse available deterministic solvers. While being very flexible and easy to implement, this technique features a very slow convergence rate. In the last few years, other approaches have been proposed, which in certain situations feature a much faster convergence rate. They exploit the possible regularity that the solution might have with respect to the input parameters, which opens up the possibility to use deterministic approximations of the response function (i.e. the solution of the problem as a function of the input parameters) based on global polynomials. Such approximations are expected to yield a very fast convergence. Stochastic Galerkin [10, 22, 38, 3, 32] and Stochastic Collocation [2, 27, 26, 37] are among these techniques. Such new techniques have been successfully applied to stochastic elliptic and parabolic PDEs. In particular, it is shown that, under particular assumptions, the solution of these problems is analytic with respect to the input random variables [2, 25]. The convergence results are then derived from the regularity results. For stochastic hyperbolic problems, the analysis was not well developed until very recently, see [24, 23, 1]. In the case of linear problems, there are a few works on the one-dimensional scalar advection equation with a time- and space-independent random wave speed [36, 13, 31]. Such problems also possess high regularity properties provided the data live in suitable spaces. The main difficulty however arises when the coefficients vary in space or time and are possibly non-smooth. In this more general case, the solution of linear hyperbolic problems may have lower regularity than those of elliptic, parabolic and hyperbolic problems with constant random coefficients [24, 23, 1]. There are also recent works on stochastic nonlinear conservation laws, see for instance [18, 19, 28, 34, 35]. In the present notes, we assume that the uncertainty in the input data is parametrized either in terms of a finite number of random variables or more generally by random fields. Random fields can in turn be accurately approximated by a finite number of random variables when the input data vary slowly in space, with a correlation length comparable to the size of the physical domain. A possible way to describe such random fields is to use the truncated Karhunen-Loéve [20, 21] or polynomial chaos expansion [39]. We will address theoretical issues of well-posedness and stochastic regularity, as well as efficient treatments and computational error analysis for hyperbolic problems with random inputs. We refer to [24, 23, 1] for further details. 2 Problem Statement In this section, for simplicity, we consider the linear second order scalar acoustic wave equation with random wave speed and deterministic boundary and initial conditions. We motivate and describe the source of randomness and address the well-posedness of the problem. For extensions to the system of elastic wave equations see [23, 1]. Let D be a convex bounded polygonal domain in R d, d = 2, 3, and (Ω, F, P ) be a complete probability space. Here, Ω is the set of outcomes, F 2 Ω is the σ-algebra of events and P : F [0, 1] is a probability measure. Consider the stochastic initial boundary value problem (IBVP): find a random function u : [0, T ] D Ω R, such that P -almost everywhere in Ω, i.e. almost surely (a.s), the following holds u tt (t, x, ω) (a 2 (x, ω) u(t, x, ω) ) = f(t, x) in [0, T ] D Ω, u(0, x, ω) = g 1 (x), u t (0, x, ω) = g 2 (x) on {t = 0} D Ω, (1) u(t, x, ω) = 0 on [0, T ] D Ω. 2
3 Here, the solution u is the displacement, and t and x = (x 1,..., x d ) are the time and location, respectively, and the data f L 2 ((0, T ); L 2 (D)), g 1 H 1 0 (D), g 2 L 2 (D), (2) are compatible. The only source of randomness is the wave speed a which is assumed to be bounded and uniformly coercive, 0 < a min a(x, ω) a max <, almost everywhere in D, a.s. (3) Assumption (3) guarantees that the energy is conserved and therefore the stochastic IBVP (1) is well-posed [24]. In many wave propagation problems, the source of randomness can be described or approximated by only a small number of uncorrelated random variables. For example, in seismic applications, a typical situation is the case of layered materials where the wave speeds in the layers are not perfectly known and therefore are described by uncorrelated random variables. The number of random variables is therefore the number of layers. In this case the randomness is described by a finite number of random variables. Another situation is when the wave speeds in layers are given by random fields, which in turn are approximated by a truncated Karhunen-Loéve expansion. Hence, the number of random variables corresponds to the number of layers as well as the number of terms in the expansion. In this case the randomness is approximated by a finite number of random variables. This motivates us to make the following finite dimensional noise assumption on the form of the wave speed, a(x, ω) = a(x, Y 1 (ω),..., Y N (ω)), almost everywhere in D, a.s, (4) where N N + and Y = [Y 1,..., Y N ] R N is a random vector. We denote by Γ n Y n (Ω) the image of Y n and assume that Y n is bounded. We let Γ = N n=1 Γ n and assume further that the random vector Y has a bounded joint probability density function ρ : Γ R + with ρ L (Γ). We note that by using a similar approach to [2, 5] we can also treat unbounded random variables, such as Gaussian and exponential variables. The finite dimensional noise assumption (4) implies that the solution of the stochastic IBVP (1) can be described by only N random variables, i.e., u(t, x, ω) = u(t, x, Y 1 (ω),..., Y N (ω)). This turns the original stochastic problem into a deterministic IBVP for the wave equation with an N- dimensional parameter, which allows the use of standard finite difference and finite element methods to approximate the solution of the resulting deterministic problem u = u(t, x, Y ), where t [0, T ], x D, and Y Γ. Note that the knowledge of u = u(t, x, Y ) fully determines the law of the random field u = u(t, x, ω). Here, as an example of form (4), we consider a random speed a given by a(x, ω) = a 0 (x) + N a n (x, ω) I Dn (x), a n (x, ω) = Y n (ω) α n (x), α n C (D n ), (5) n=1 where, I is the indicator function. In this case, D is a heterogeneous medium consisting of N nonoverlapping sub-domains {D n } N n=1, {Y n } N n=1 are independent random variables, and {α n } N n=1 are smooth functions defined on sub-domains. The boundaries of sub-domains, which are interfaces of speed discontinuity, are assumed to be smooth. The ultimate goal is the prediction of statistical moments of the solution u or statistics of some given quantities of physical interest. As an example, we consider the following quantity, T Q(Y ) = u(t, x, Y ) φ(x) dx dt + u(t, x, Y ) ψ(x) dx, (6) 0 D 3 D
4 where u solves (1) and the mollifiers φ and ψ are given functions of x. 3 Non-intrusive Numerical Methods There are in general three types of methods for propagating uncertainty in PDE models with random inputs; intrusive, non-intrusive, and hybrid methods. Intrusive methods, such as perturbation expansion, intrusive polynomial chaos [39], and stochastic Galerkin [10], require extensive modifications in existing deterministic solvers. On the contrary, non-intrusive methods, such as Monte Carlo and Stochastic Collocation, are sample-based approaches. They rely on a set of deterministic models corresponding to a set of realizations, and hence require no modification to the existing deterministic solvers. Finally, hybrid methods are a mixture of both intrusive and non-intrusive approaches. Non-intrusive (or sample-based) methods are attractive in the sense that they require only a deterministic solver for computing deterministic models. In addition, since the deterministic models are independent, it is possible to distribute them onto multiple processors and perform parallel computation. The most popular non-intrusive technique is the Monte Carlo method. While being very flexible and easy to implement, this technique features a very slow convergence rate. The Multi-Level Monte Carlo method has been proposed to accelerate the slow convergence of Monte Carlo sampling [12]. Quasi Monte Carlo methods can be considered as well, that aim at achieving higher rates of convergence [15]. The stochastic collocation (SC) method is another non-intrusive technique, in which, regularity features of the quantity of interest with respect to the input parameters can be exploited to obtain a much faster or possibly spectral convergence. In SC, the problem(1) is first discretized in space and time, using a deterministic numerical method. The obtained semi-discrete problem is then collocated in a set of η collocation points {Y (k) } η k=1 Γ to compute the approximate solutions u h(t, x, Y (k) ), where h represent the discretization mesh/grid size. A global polynomial approximation is then built upon those evaluations, u h,η (t, x, Y ) = η u h (t, x, Y (k) ) L k (Y ), k=1 for suitable multivariate polynomial bases {L k } η k=1 such as Lagrange polynomials. Finally, using the Gauss quadrature formula, we can easily approximate the statistical moments of the solution. For instance, E[u(., Y )] E[u h,η (., Y )] = u h,η (., Y ) ρ(y ) dy = η u h (., Y (k) ) k=1 Γ Γ L k (Y ) ρ(y ) dy η u h (., Y (k) ) θ k. A key point in SC is the choice of the set of collocation points {Y (k) }, i.e. the type of computational grid in the N-dimensional stochastic space. A full tensor grid, based on cartesian product of mono-dimensional grids, can only be used when the number of stochastic dimensions N is small, since the computational cost grows exponentially fast with N (curse of dimensionality). Alternatively, sparse grids can reduce the curse of dimensionality. They were originally introduced by Smolyak for high dimensional quadrature and interpolation computations [30]. In the following we will briefly review the sparse grid construction. k=1 4
5 Let j Z N + be a multi-index containing non-negative integers. For a non-negative index j n in j, we introduce a sequence of one-dimensional polynomial interpolant operators U jn : C 0 (Γ n ) P p(jn)(γ n ) of degree p(j n ) on p(j n ) + 1 suitable knots. With U 1 = 0, we define the detail operator jn := U jn U jn 1. Finally, introducing a sequence of index sets I(l) Z N +, the sparse grid approximation of u : Γ V at level l reads u η(l,n) (., Y ) = N jn [u](., Y ). (7) j I(l) n=1 Furthermore, in order for the sum (7) to have some telescopic properties, which are desirable, we impose an additional admissibility condition on the set I [9]. An index set I is said to be admissible if j I, j e n I for 1 n N, j n 1, holds. Here, e n is the n-th canonical unit vector. To fully characterize the sparse approximation operator in (7), we need to provide the following: A level l N and a function p(j) representing the relation between an index j and the number of points in the corresponding one-dimensional polynomial interpolation formula U j. A sequence of sets I(l). Typical examples include total degree and hyperbolic cross grids. The family of points to be used, such as Gauss or Clenshaw-Curtis abscissae, [33]. The rate of convergence for SC depends on the stochastic regularity of the output quantity of interest, which may be the solution or some given functional of the solution. For example, a fast spectral convergence is possible for highly regular outputs. We will discuss this issue in more details in the next section. 4 Stochastic Regularity and Convergence of Stochastic Collocation As extensively discussed in [24, 23, 1], the error in the stochastic collocation method is related to the stochastic regularity of the output quantity of interest (the solution or a given functional of the solution). It has been shown that, under particular assumptions, the solution of stochastic elliptic and parabolic PDEs is analytic with respect to the input random variables [2, 25]. In contrast, the solution of stochastic hyperbolic PDEs may possess lower regularity. In this section, first, we will convey general stochastic regularity properties of hyperbolic PDEs through simple examples. Then, we state the main regularity results followed by the convergence results for SC. 4.1 Examples: General Stochastic Regularity Properties Example 1. Consider the 1D Cauchy problem for the scalar wave equation, u tt (t, x, y) y 2 u xx (t, x, y) = 0, in [0, T ] R, u(0, x, y) = g(x), u t (0, x, y) = 0, on {t = 0} R, 5
6 with g C0 (R). The wave speed is constant and given by a single random variable y. We want to investigate the regularity of the solution u with respect to y. For this purpose, we extend the parameter y into the complex plane and study the extended problem in the complex plane. It is well known that if the extended problem is well-posed and the first derivative of the resulting complexvalued solution with respect to the parameter satisfies the so called Cauchy-Riemann conditions, the solution can analytically be extended into the complex plane, and hence u will be analytic with respect to y. We therefore let y = y R + i y I, where y R, y I R, and apply the Fourier transform with respect to x to obtain û(t, k, y) = ĝ(k) ( e i y k t + e i y k t), 2 where û and ĝ are the Fourier transforms of u and g with respect to x, respectively. Hence, û(t, k, y) ĝ(k) e y I k t. Therefore, the Fourier transform of the solution û(t, k, y) grows exponentially fast in time, unless the Fourier transform of the initial solution ĝ(k) decays faster than e y I k t. In order for the Cauchy problem to be well-posed, g needs to belong to the Gevrey space G q (R) of order q < 1 [29]. For ɛ k 1/q such functions, we have ĝ(k) C e for positive constants C and ɛ, and hence the problem is well-posed in the complex strip Σ r = {(y R + i y I ) C : y I r}. Note that G 1 (R) is the space of analytic functions. If g G 1 (R), the problem is well-posed only for a finite time interval when t ɛ/r. This shows that even if the initial solution g is analytic, the solution u is not analytic with respect to y for all times in Σ r. Example 2. Consider the 1D Cauchy problem for the scalar wave equation in a domain consisting of two homogeneous half-spaces separated by an interface at x = 0, u tt (t, x, Y ) ( a 2 (x, Y ) u x (t, x, Y ) ) = 0, in [0, T ] R, x u(0, x, Y ) = g( x), u t (0, x, Y ) = a g ( x), on {t = 0} R, with g C0 (0, ). The wave speed is piecewise constant and a function of a random vector of two variables Y = [y, y + ], { y, x < 0, a(x, Y ) = y +, x > 0. By d Alembert s formula and the interface jump conditions at x = 0, the solution reads u(t, 0, Y ) = u(t, 0 +, Y ), y 2 u x (t, 0, Y ) = y 2 + u x (t, 0 +, Y ), (8) u(t, x, Y ) = { g(y t x) + y y+ y +y + g(y t + x), x < 0, 2 y y +y + g( y y + (y + t x)), x > 0. (9) Clearly, the solution (9) is infinitely differentiable with respect to both parameters y and y + in (0, + ). Note that the smooth initial solution u(0, x, Y ), which is contained in one layer with zero value at the interface, automatically satisfies the interface conditions (8) at time zero. Otherwise, if for instance the initial solution crosses the interface without satisfying (8), a singularity is introduced in the solution, and the high regularity result does not hold any longer. 6
7 In the more general case of multi-dimensional heterogeneous media consisting of sub-domains, the interface jump conditions on a smooth interface Υ between two sub-domains D I and D II are given by [u(t,., Y )] Υ = 0, [a 2 (., Y ) u n (t,., Y )] Υ = 0. (10) Here, the subscript n represents the normal derivative, and [v(.)] Υ is the jump in the function v across the interface Υ. In this general case, the high regularity with respect to parameters holds provided the smooth initial solution satisfies (10). The jump conditions are satisfied for instance when the initial data are contained within sub-domains. This result for Cauchy problems can easily be extended to IBVPs by splitting the problem to one pure Cauchy and two half-space problems. See [24] for more details. We note that the above high stochastic regularity result is valid only for particular types of smooth data. In real applications, the data are not smooth. Let us now consider a more practical datum. Example 3. Consider the 1D Cauchy problem for the scalar wave equation in a domain consisting of two homogeneous half-spaces separated by an interface at x = 0, u tt (t, x, y) ( a 2 (x, y) u x (t, x, y) ) = 0, in [0, T ] R, x u(0, x, y) = g(x), u t (0, x, y) = 0, on {t = 0} R, with g H0 1 (R) being a hat function with a narrow support supp g = [x 0 α, x 0 +α], with 0 < α 1, located at x 0 < α. The wave speed is piecewise constant and a function of a random variable y, { 1, x < 0, a(x, y) = y, x > 0. Similar to Example 2, we can obtain the closed form of the solution as 1 2 g(t + x), x < x 0 t 12 g( t + x) + 1 y 2 (1+y) u(t, x, y) = g( t x), x 0 t < x < y g( t + x y ), 0 < x < x 0 + y t 0, x > x 0 + y t Since g H 1 0, there is only one bounded derivative y u, and higher bounded y-derivatives do not exist. However, if we consider a quantity of interest Q(y) as in (6) with mollifiers φ, ψ C 0 (R) vanishing at x < 0, we can employ integration by parts and shift the derivatives on g to the mollifiers. Therefore, although the solution u has only one bounded y-derivative, the quantity of interest Q is smooth with respect to y. Remark 1. Immediate results of the above three examples are the following: 1. For the solution of the 1D Cauchy problem for the wave equation to be analytic with respect to the random wave speed at all times in a given complex strip Σ r with r > 0, the initial datum needs to live in a space strictly contained in the space of analytic functions, which is the Gevrey space G q (R) with 0 < q < 1. Moreover, if the problem is well-posed and the data are analytic, the solution may be analytic with respect to the parameter in Σ r only for a short time interval. 2. In a 1D heterogeneous medium with piecewise smooth wave speeds, if the data are smooth and the initial solution satisfies the interface jump conditions (10), the solution to the Cauchy problem is smooth with respect to the wave speeds. If the initial solution does not satisfy (10), the solution is not smooth with respect to the wave speeds. 3. In general, the solution u has only finite stochastic regularity. The stochastic regularity of functionals of the solution (such as mollified quantities of interest) can however be considerably higher than the regularity of the solution. 7
8 4.2 General Results: Stochastic Regularity We now state some regularity results in the more general case when the data satisfy the minimal assumptions (2). See [24, 23] for proofs. Theorem 1. For the solution of the stochastic IBVP (1) with data given by (2) and a random piecewise smooth wave speed satisfying (3) and (5), we have Y u C 0 (0, T ; L 2 (D)), Y Γ. Theorem 2. Consider the quantity of interest Q(Y ) in (6) and let the smooth mollifiers φ, ψ C 0 (D) vanish at the discontinuity interfaces. Then with the assumptions of Theorem 1, we have Q C (Γ). In practical applications, quantities of interest, such as Arias intensity, spectral acceleration, and von Mises stress, are usually nonlinear in u. In such cases, the high stochastic regularity property might not hold. However, one can perform a low-pass filtering (LPF) on the low regular solutions or quantities of interest by convolving them with smooth kernels such as Gaussian functions, K δ (x) = 1 ( 2 π δ) d exp ( x 2 2 δ 2 ), x R d. (11) Here, the standard deviation δ is inversely proportional to the maximum frequency that is allowed to pass. For instance, the filtered solution is given by the convolution Q δ (u)(t, x, Y ) = (u K δ )(t, x, Y ) = K δ (x x) u(t, x, Y ) d x. (12) The filtered solution (12) is of a type similar to the quantity of interest (6) with smooth mollifiers. However, the main difference here is the boundary effects introduced by the convolution. Therefore, in the presence of a compactly supported smooth kernel K C 0 (D) as mollifier, the quantity (12) will have high Stochastic regularity on a smaller domain x D δ D with dist( D δ, D) d > 0. We choose d so that K δ (x 0 ) with x 0 = d is essentially zero. This implies that for any x D δ, the support of K δ (x x) is essentially vanishing at D. We note that by performing filtering, we introduce an error u Q δ (u) which is proportional to O(δ 2 ). This error needs to be taken into account. We refer to [1] for a more rigorous error analysis of filtered quantities. 4.3 Convergence Results for Stochastic Collocation In order to obtain a priori estimates for the total error u u h,η W L 2 ρ (Γ), with W := L 2 (0, T ; L 2 (D)), we split it into two parts D ε := u u h,η u u h + u h u h,η =: ε I + ε II. (13) The first term ε I controls the convergence of the deterministic numerical scheme with respect to the mesh size h and is of order O(h r ), where r is the minimum between the order of accuracy of the finite element or finite difference method used and the regularity of the solution. Notice that the constant in the term O(h r ) is uniform with respect to Y. The second term ε II is derived as follows from the stochastic regularity results. See [24] for proofs. 8
9 Theorem 3. Consider the isotropic full tensor product interpolation. Then ε II C η s/n. (14) where the constant C = C s N n=1 max k=0,...,s D k Y n u h,l L (Γ;W ) does not depend on l. Theorem 4. Consider the Smolyak sparse tensor product interpolation based on Gauss-Legendre abscissas, and let u h,η be given by (7). Then for the solution u h with s 1 bounded mixed Y - derivatives, where the constant ε II C ( 1 + log 2 1 Cs N C = C s ρ 1/2 1 C s depends on u h, s, and N, but not on l. η ) 2 N log 2 η s ξ+log N, ξ = 1 + log 2 (1 + log2 1.5) 2.1, (15) N max d=1,...,n max 0 k 1,...,k d s Dk1 Y 1... D k d Y d u h L 2 (Γ;W ), Remark 2. It is possible to show that the semi-discrete solution u h can analytically be extended on the whole region Σ(Γ, τ) = {Z C N, dist(γ n, Z n ) τ, n = 1,..., N}, with the radius of analyticity τ = O(h) [24]. This can be used to show that for both full tensor and Smolyak interpolation, we will have a fast exponential decay in the error when the product h p(l) is large. As a result, with a fixed h, the error convergence is slow (algebraic) for a small l and fast (exponential) for a large l. Moreover, the rate of convergence deteriorates as h gets smaller. Remark 3. The main parameters in the computations include h, η(l), and possibly δ if filtered quantities are used. In order to find the optimal choice of the parameters, we need to minimize the computational complexity of the SC method, subject to the total error constraint ε = TOL. We refer to [1] for more details. References [1] I. Babuška, M. Motamed, and R. Tempone. A stochastic multiscale method for the elastodynamic wave equations arising from fiber composites Preprint. [2] I. Babuška, F. Nobile, and R. Tempone. A stochastic collocation method for elliptic partial differential equations with random input data. SIAM J. Numer. Anal., 45: , [3] I. Babuška, R. Tempone, and G. E. Zouraris. Solving elliptic boundary value problems with uncertain coefficients by the finite element method: the stochastic formulation. Computer Methods in Applied Mechanics and Engineering, 194: , [4] J. P. Boyd. Chebyshev and Fourier Spectral Methods. Springer, [5] J. Charrier. Strong and weak error estimates for elliptic partial differential equations with random coefficients. IMA J. Numer. Anal., 50: , [6] K. Eriksson, D. Estep, P. Hansbo, and C. Johnson. Computational differential equations. Studentlitteratur,
10 [7] G. S. Fishman. Monte Carlo: Concepts, Algorithms, and Applications. Springer- Verlag, New York, [8] K. O. Friedrichs. Symmetric hyperbolic linear differential equations. Comm. Pure Appl. Math., 7: , [9] T. Gerstner and M. Griebel. Dimension-adaptive tensor-product quadrature. Computing, 71:65 87, [10] R. G. Ghanem and P. D. Spanos. Stochastic finite elements: A spectral approach. Springer, New York, [11] W. Gibson. The Method of Moments in Electromagnetics. Chapman and Hall/CRC, [12] M. B. Giles. Multilevel Monte Carlo path simulation. Oper. Res., 56: , [13] D. Gottlieb and D. Xiu. Galerkin method for wave equations with uncertain coefficients. Commun. Comput. Phys., 3: , [14] B. Gustafsson, H. O. Kreiss, and J. Ologer. Time dependent problems and difference methods. Wiley-Interscience, [15] T. Y. Hou and X. Wu. Quasi-Monte Carlo methods for elliptic PDEs with random coefficients and applications. J. Comput. Phys., 230: , [16] H.-O. Kreiss. Initial boundary value problems for hyperbolic systems. Comm. Pure Appl. Math., 23: , [17] R. J. Leveque. Finite volume methods for hyperbolic problems. Cambridge University Press, [18] G. Lin, C.-H. Su, and G. E. Karniadakis. Predicting shock dynamics in the presence of uncertainties. J. Comput. Phys., 217: , [19] G. Lin, C.-H. Su, and G. E. Karniadakis. Stochastic modeling of random roughness in shock scattering problems: theory and simulations. Comput. Methods Appl. Mech. Engrg., 197: , [20] M. Loéve. Probability theory I. Grad. Texts in Math. 45. Springer-Verlag, New York, [21] M. Loéve. Probability theory II. Grad. Texts in Math. 46. Springer-Verlag, New York, [22] H. G. Matthies and A. Kees. Galerkin methods for linear and nonlinear elliptic stochastic partial differential equations. Computer Methods in Applied Mechanics and Engineering, 194: , [23] M. Motamed, F. Nobile, and R. Tempone. Analysis and computation of the elastic wave equation with random coefficients Submitted. [24] M. Motamed, F. Nobile, and R. Tempone. A stochastic collocation method for the second order wave equation with a discontinuous random speed. Numer. Math., 123: , [25] F. Nobile and R Tempone. Analysis and implementation issues for the numerical approximation of parabolic equations with random coefficients. Int. J. Numer. Meth. Engng., 80: ,
11 [26] F. Nobile, R. Tempone, and C. G. Webster. An anisotropic sparse grid stochastic collocation method for partial differential equations with random input data. SIAM J. Numer. Anal., 46: , [27] F. Nobile, R. Tempone, and C. G. Webster. A sparse grid stochastic collocation method for partial differential equations with random input data. SIAM J. Numer. Anal., 46: , [28] G. Poette, B. Després, and D. Lucor. Uncertainty quantification for systems of conservation laws. J. Comput. Phys., 228: , [29] L. Rodino. Linear partial differential operators in Gevrey spaces. World Scientific Publishing Co., Singapore, [30] S. A. Smolyak. Quadrature and interpolation formulas for tensor products of certain classes of functions. Doklady Akademii Nauk SSSR, 4: , [31] T. Tang and T. Zhou. Convergence analysis for stochastic collocation methods to scalar hyperbolic equations with a random wave speed. Commun. Comput. Phys., 8: , [32] R. A. Todor and C. Schwab. Convergence rates for sparse chaos approximations of elliptic problems with stochastic coefficients. IMA J. Numer. Anal., 27: , [33] L. N. Trefethen. Is Gauss quadrature better than Clenshaw-Curtis? SIAM Rev., 50:67 87, [34] J. Tryoen, O. Le Maitre, M. Ndjinga, and A. Ern. Intrusive projection methods with upwinding for uncertain nonlinear hyperbolic systems. J. Comput. Phys, 229: , [35] J. Tryoen, O. Le Maitre, M. Ndjinga, and A. Ern. Roe solver with entropy corrector for uncertain hyperbolic systems. Comput. Appl. Math., 235: , [36] X. Wang and G. E. Karniadakis. Long-term behavior of ploynomial chaos in stochastic flow simulations. Computer Methods in Applied Mechanics and Engineering, 195: , [37] D. Xiu and J. S. Hesthaven. High-order collocation methods for differential equations with random inputs. SIAM J. Sci. Comput., 27: , [38] D. Xiu and G. E. Karniadakis. Modeling uncertainty in steady state diffusion problems via generalized polynomial chaos. Computer Methods in Applied Mechanics and Engineering, 191: , [39] D. Xiu and G. E. Karniadakis. The Wiener-Askey polynomial chaos for stochastic differential equations. SIAM J. Sci. Comput., 24: ,
Numerische Mathematik
Numer. Math. (203) :493 536 OI 0.007/s002-02-0493-5 Numerische Mathematik A stochastic collocation method for the second order wave equation with a discontinuous random speed Mohammad Motamed Fabio Nobile
More informationPerformance Evaluation of Generalized Polynomial Chaos
Performance Evaluation of Generalized Polynomial Chaos Dongbin Xiu, Didier Lucor, C.-H. Su, and George Em Karniadakis 1 Division of Applied Mathematics, Brown University, Providence, RI 02912, USA, gk@dam.brown.edu
More informationStochastic Collocation Methods for Polynomial Chaos: Analysis and Applications
Stochastic Collocation Methods for Polynomial Chaos: Analysis and Applications Dongbin Xiu Department of Mathematics, Purdue University Support: AFOSR FA955-8-1-353 (Computational Math) SF CAREER DMS-64535
More informationBeyond Wiener Askey Expansions: Handling Arbitrary PDFs
Journal of Scientific Computing, Vol. 27, Nos. 1 3, June 2006 ( 2005) DOI: 10.1007/s10915-005-9038-8 Beyond Wiener Askey Expansions: Handling Arbitrary PDFs Xiaoliang Wan 1 and George Em Karniadakis 1
More informationSTOCHASTIC SAMPLING METHODS
STOCHASTIC SAMPLING METHODS APPROXIMATING QUANTITIES OF INTEREST USING SAMPLING METHODS Recall that quantities of interest often require the evaluation of stochastic integrals of functions of the solutions
More informationResearch Article A Pseudospectral Approach for Kirchhoff Plate Bending Problems with Uncertainties
Mathematical Problems in Engineering Volume 22, Article ID 7565, 4 pages doi:.55/22/7565 Research Article A Pseudospectral Approach for Kirchhoff Plate Bending Problems with Uncertainties Ling Guo, 2 and
More informationMulti-Element Probabilistic Collocation Method in High Dimensions
Multi-Element Probabilistic Collocation Method in High Dimensions Jasmine Foo and George Em Karniadakis Division of Applied Mathematics, Brown University, Providence, RI 02912 USA Abstract We combine multi-element
More informationUncertainty Quantification for multiscale kinetic equations with random inputs. Shi Jin. University of Wisconsin-Madison, USA
Uncertainty Quantification for multiscale kinetic equations with random inputs Shi Jin University of Wisconsin-Madison, USA Where do kinetic equations sit in physics Kinetic equations with applications
More informationUncertainty Quantification and hypocoercivity based sensitivity analysis for multiscale kinetic equations with random inputs.
Uncertainty Quantification and hypocoercivity based sensitivity analysis for multiscale kinetic equations with random inputs Shi Jin University of Wisconsin-Madison, USA Shanghai Jiao Tong University,
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
More informationQuasi-optimal and adaptive sparse grids with control variates for PDEs with random diffusion coefficient
Quasi-optimal and adaptive sparse grids with control variates for PDEs with random diffusion coefficient F. Nobile, L. Tamellini, R. Tempone, F. Tesei CSQI - MATHICSE, EPFL, Switzerland Dipartimento di
More informationNONLOCALITY AND STOCHASTICITY TWO EMERGENT DIRECTIONS FOR APPLIED MATHEMATICS. Max Gunzburger
NONLOCALITY AND STOCHASTICITY TWO EMERGENT DIRECTIONS FOR APPLIED MATHEMATICS Max Gunzburger Department of Scientific Computing Florida State University North Carolina State University, March 10, 2011
More informationIntroduction to Uncertainty Quantification in Computational Science Handout #3
Introduction to Uncertainty Quantification in Computational Science Handout #3 Gianluca Iaccarino Department of Mechanical Engineering Stanford University June 29 - July 1, 2009 Scuola di Dottorato di
More informationA Polynomial Chaos Approach to Robust Multiobjective Optimization
A Polynomial Chaos Approach to Robust Multiobjective Optimization Silvia Poles 1, Alberto Lovison 2 1 EnginSoft S.p.A., Optimization Consulting Via Giambellino, 7 35129 Padova, Italy s.poles@enginsoft.it
More informationUncertainty Quantification for multiscale kinetic equations with high dimensional random inputs with sparse grids
Uncertainty Quantification for multiscale kinetic equations with high dimensional random inputs with sparse grids Shi Jin University of Wisconsin-Madison, USA Kinetic equations Different Q Boltmann Landau
More informationImplementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs
Implementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs Roman Andreev ETH ZÜRICH / 29 JAN 29 TOC of the Talk Motivation & Set-Up Model Problem Stochastic Galerkin FEM Conclusions & Outlook Motivation
More informationStochastic Regularity of a Quadratic Observable of High Frequency Waves
Research in the Mathematical Sciences manuscript No. (will be inserted by the editor) Stochastic Regularity of a Quadratic Observable of High Frequency Waves G. Malenová M. Motamed O. Runborg Received:
More informationAdaptive Collocation with Kernel Density Estimation
Examples of with Kernel Density Estimation Howard C. Elman Department of Computer Science University of Maryland at College Park Christopher W. Miller Applied Mathematics and Scientific Computing Program
More informationSparse polynomial chaos expansions in engineering applications
DEPARTMENT OF CIVIL, ENVIRONMENTAL AND GEOMATIC ENGINEERING CHAIR OF RISK, SAFETY & UNCERTAINTY QUANTIFICATION Sparse polynomial chaos expansions in engineering applications B. Sudret G. Blatman (EDF R&D,
More informationAccuracy, Precision and Efficiency in Sparse Grids
John, Information Technology Department, Virginia Tech.... http://people.sc.fsu.edu/ jburkardt/presentations/ sandia 2009.pdf... Computer Science Research Institute, Sandia National Laboratory, 23 July
More informationUncertainty Quantification in Discrete Fracture Network Models
Uncertainty Quantification in Discrete Fracture Network Models Claudio Canuto Dipartimento di Scienze Matematiche, Politecnico di Torino claudio.canuto@polito.it Joint work with Stefano Berrone, Sandra
More informationResearch Article Multiresolution Analysis for Stochastic Finite Element Problems with Wavelet-Based Karhunen-Loève Expansion
Mathematical Problems in Engineering Volume 2012, Article ID 215109, 15 pages doi:10.1155/2012/215109 Research Article Multiresolution Analysis for Stochastic Finite Element Problems with Wavelet-Based
More informationFinal Report: DE-FG02-95ER25239 Spectral Representations of Uncertainty: Algorithms and Applications
Final Report: DE-FG02-95ER25239 Spectral Representations of Uncertainty: Algorithms and Applications PI: George Em Karniadakis Division of Applied Mathematics, Brown University April 25, 2005 1 Objectives
More informationUncertainty Quantification in Computational Models
Uncertainty Quantification in Computational Models Habib N. Najm Sandia National Laboratories, Livermore, CA, USA Workshop on Understanding Climate Change from Data (UCC11) University of Minnesota, Minneapolis,
More informationarxiv: v1 [math.na] 3 Apr 2019
arxiv:1904.02017v1 [math.na] 3 Apr 2019 Poly-Sinc Solution of Stochastic Elliptic Differential Equations Maha Youssef and Roland Pulch Institute of Mathematics and Computer Science, University of Greifswald,
More informationFast Numerical Methods for Stochastic Computations
Fast AreviewbyDongbinXiu May 16 th,2013 Outline Motivation 1 Motivation 2 3 4 5 Example: Burgers Equation Let us consider the Burger s equation: u t + uu x = νu xx, x [ 1, 1] u( 1) =1 u(1) = 1 Example:
More informationEFFICIENT STOCHASTIC GALERKIN METHODS FOR RANDOM DIFFUSION EQUATIONS
EFFICIENT STOCHASTIC GALERKIN METHODS FOR RANDOM DIFFUSION EQUATIONS DONGBIN XIU AND JIE SHEN Abstract. We discuss in this paper efficient solvers for stochastic diffusion equations in random media. We
More informationJ. Beck 1, F. Nobile 2, L. Tamellini 2 and R. Tempone 1
ESAIM: PROCEEDINGS, October 2011, Vol. 33, p. 10-21 C. Dobrzynski, T. Colin & R. Abgrall, Editors IMPLEMENTATION OF OPTIMAL GALERKIN AND COLLOCATION APPROXIMATIONS OF PDES WITH RANDOM COEFFICIENTS, J.
More informationarxiv: v2 [math.na] 8 Sep 2017
arxiv:1704.06339v [math.na] 8 Sep 017 A Monte Carlo approach to computing stiffness matrices arising in polynomial chaos approximations Juan Galvis O. Andrés Cuervo September 3, 018 Abstract We use a Monte
More informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationComparison of Clenshaw-Curtis and Leja quasi-optimal sparse grids for the approximation of random PDEs
MATHICSE Mathematics Institute of Computational Science and Engineering School of Basic Sciences - Section of Mathematics MATHICSE Technical Report Nr. 41.2014 September 2014 Comparison of Clenshaw-Curtis
More informationUncertainty Quantification in Computational Science
DTU 2010 - Lecture I Uncertainty Quantification in Computational Science Jan S Hesthaven Brown University Jan.Hesthaven@Brown.edu Objective of lectures The main objective of these lectures are To offer
More informationA sparse grid stochastic collocation method for elliptic partial differential equations with random input data
A sparse grid stochastic collocation method for elliptic partial differential equations with random input data F. Nobile R. Tempone C. G. Webster June 22, 2006 Abstract This work proposes and analyzes
More informationDownloaded 08/12/13 to Redistribution subject to SIAM license or copyright; see
SIAM J. MATRIX ANAL. APPL. Vol. 33, No. 1, pp. 195 209 c 2012 Society for Industrial and Applied Mathematics A POSTERIORI ERROR ANALYSIS OF PARAMETERIZED LINEAR SYSTEMS USING SPECTRAL METHODS T. BUTLER,
More informationSolving the stochastic steady-state diffusion problem using multigrid
IMA Journal of Numerical Analysis (2007) 27, 675 688 doi:10.1093/imanum/drm006 Advance Access publication on April 9, 2007 Solving the stochastic steady-state diffusion problem using multigrid HOWARD ELMAN
More informationSolving the steady state diffusion equation with uncertainty Final Presentation
Solving the steady state diffusion equation with uncertainty Final Presentation Virginia Forstall vhfors@gmail.com Advisor: Howard Elman elman@cs.umd.edu Department of Computer Science May 6, 2012 Problem
More informationASSESSMENT OF COLLOCATION AND GALERKIN APPROACHES TO LINEAR DIFFUSION EQUATIONS WITH RANDOM DATA
International Journal for Uncertainty Quantification, 11):19 33, 2011 ASSESSMENT OF COLLOCATION AND GALERKIN APPROACHES TO LINEAR DIFFUSION EQUATIONS WITH RANDOM DATA Howard C. Elman, 1, Christopher W.
More informationA study of hyperbolicity of kinetic stochastic Galerkin system for the isentropic Euler equations with uncertainty
A study of hyperbolicity of kinetic stochastic Galerkin system for the isentropic Euler equations with uncertainty Shi Jin and Ruiwen Shu September 19, 2018 Abstract We study the fluid dynamic behavior
More informationSparse Quadrature Algorithms for Bayesian Inverse Problems
Sparse Quadrature Algorithms for Bayesian Inverse Problems Claudia Schillings, Christoph Schwab Pro*Doc Retreat Disentis 2013 Numerical Analysis and Scientific Computing Disentis - 15 August, 2013 research
More informationIntroduction to Computational Stochastic Differential Equations
Introduction to Computational Stochastic Differential Equations Gabriel J. Lord Catherine E. Powell Tony Shardlow Preface Techniques for solving many of the differential equations traditionally used by
More informationQuadrature for Uncertainty Analysis Stochastic Collocation. What does quadrature have to do with uncertainty?
Quadrature for Uncertainty Analysis Stochastic Collocation What does quadrature have to do with uncertainty? Quadrature for Uncertainty Analysis Stochastic Collocation What does quadrature have to do with
More informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationICES REPORT A posteriori error control for partial differential equations with random data
ICES REPORT 13-8 April 213 A posteriori error control for partial differential equations with random data by Corey M. Bryant, Serge Prudhomme, and Timothy Wildey The Institute for Computational Engineering
More informationAN ADAPTIVE REDUCED BASIS COLLOCATION METHOD BASED ON PCM ANOVA DECOMPOSITION FOR ANISOTROPIC STOCHASTIC PDES
AN ADAPTIVE REDUCED BASIS COLLOCATION METHOD BASED ON PCM ANOVA DECOMPOSITION FOR ANISOTROPIC STOCHASTIC PDES Heyrim Cho & Howard Elman 2 Department of Mathematics, University of Maryland, College Park,
More informationA stochastic Galerkin method for general system of quasilinear hyperbolic conservation laws with uncertainty
A stochastic Galerkin method for general system of quasilinear hyperbolic conservation laws with uncertainty Kailiang Wu, Huazhong Tang arxiv:6.4v [math.na] 6 Jan 6 HEDPS, CAPT & LMAM, School of Mathematical
More informationCollocation methods for uncertainty quantification in PDE models with random data
Collocation methods for uncertainty quantification in PDE models with random data Fabio Nobile CSQI-MATHICSE, EPFL, Switzerland Acknowlegments: L. Tamellini, G. Migliorati (EPFL), R. Tempone, (KAUST),
More informationDownloaded 08/15/18 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 40, No. 4, pp. A2152 A2173 c 2018 Society for Industrial and Applied Mathematics STOCHASTIC DOMAIN DECOMPOSITION VIA MOMENT MINIMIZATION DONGKUN ZHANG, HESSAM BABAEE, AND GEORGE
More informationSparse Grids. Léopold Cambier. February 17, ICME, Stanford University
Sparse Grids & "A Dynamically Adaptive Sparse Grid Method for Quasi-Optimal Interpolation of Multidimensional Analytic Functions" from MK Stoyanov, CG Webster Léopold Cambier ICME, Stanford University
More informationSimulating with uncertainty : the rough surface scattering problem
Simulating with uncertainty : the rough surface scattering problem Uday Khankhoje Assistant Professor, Electrical Engineering Indian Institute of Technology Madras Uday Khankhoje (EE, IITM) Simulating
More informationPascal s Triangle on a Budget. Accuracy, Precision and Efficiency in Sparse Grids
Covering : Accuracy, Precision and Efficiency in Sparse Grids https://people.sc.fsu.edu/ jburkardt/presentations/auburn 2009.pdf... John Interdisciplinary Center for Applied Mathematics & Information Technology
More informationOptimal polynomial approximation of PDEs with stochastic coefficients by Galerkin and collocation methods
Optimal polynomial approximation of PDEs with stochastic coefficients by Galerkin and collocation methods Fabio Nobile MOX, Department of Mathematics, Politecnico di Milano Seminaire du Laboratoire Jacques-Louis
More informationLucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche
Scuola di Dottorato THE WAVE EQUATION Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche Lucio Demeio - DIISM wave equation 1 / 44 1 The Vibrating String Equation 2 Second
More informationarxiv: v1 [math.na] 14 Sep 2017
Stochastic collocation approach with adaptive mesh refinement for parametric uncertainty analysis arxiv:1709.04584v1 [math.na] 14 Sep 2017 Anindya Bhaduri a, Yanyan He 1b, Michael D. Shields a, Lori Graham-Brady
More informationGiovanni Migliorati. MATHICSE-CSQI, École Polytechnique Fédérale de Lausanne
Analysis of the stability and accuracy of multivariate polynomial approximation by discrete least squares with evaluations in random or low-discrepancy point sets Giovanni Migliorati MATHICSE-CSQI, École
More informationOverview. Bayesian assimilation of experimental data into simulation (for Goland wing flutter) Why not uncertainty quantification?
Delft University of Technology Overview Bayesian assimilation of experimental data into simulation (for Goland wing flutter), Simao Marques 1. Why not uncertainty quantification? 2. Why uncertainty quantification?
More informationMultilevel stochastic collocations with dimensionality reduction
Multilevel stochastic collocations with dimensionality reduction Ionut Farcas TUM, Chair of Scientific Computing in Computer Science (I5) 27.01.2017 Outline 1 Motivation 2 Theoretical background Uncertainty
More informationSolving the Stochastic Steady-State Diffusion Problem Using Multigrid
Solving the Stochastic Steady-State Diffusion Problem Using Multigrid Tengfei Su Applied Mathematics and Scientific Computing Advisor: Howard Elman Department of Computer Science Sept. 29, 2015 Tengfei
More informationStochastic collocation and mixed finite elements for flow in porous media
Stochastic collocation and mixed finite elements for flow in porous media Benjamin Ganis a, Hector Klie b, Mary F. Wheeler b, Tim Wildey b, Ivan Yotov a, and Dongxiao Zhang c a Department of Mathematics,
More informationELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS
ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationA MULTILEVEL STOCHASTIC COLLOCATION METHOD FOR PARTIAL DIFFERENTIAL EQUATIONS WITH RANDOM INPUT DATA
A MULTILEVEL STOCHASTIC COLLOCATION METHOD FOR PARTIAL DIFFERENTIAL EQUATIONS WITH RANDOM INPUT DATA A. L. TECKENTRUP, P. JANTSCH, C. G. WEBSTER, AND M. GUNZBURGER Abstract. Stochastic collocation methods
More informationSteps in Uncertainty Quantification
Steps in Uncertainty Quantification Challenge: How do we do uncertainty quantification for computationally expensive models? Example: We have a computational budget of 5 model evaluations. Bayesian inference
More informationAn Empirical Chaos Expansion Method for Uncertainty Quantification
An Empirical Chaos Expansion Method for Uncertainty Quantification Melvin Leok and Gautam Wilkins Abstract. Uncertainty quantification seeks to provide a quantitative means to understand complex systems
More informationA sparse grid stochastic collocation method for partial differential equations with random input data
A sparse grid stochastic collocation method for partial differential equations with random input data F. Nobile R. Tempone C. G. Webster May 24, 2007 Abstract This work proposes and analyzes a Smolyak-type
More informationAn operator splitting based stochastic Galerkin method for the one-dimensional compressible Euler equations with uncertainty
An operator splitting based stochastic Galerkin method for the one-dimensional compressible Euler equations with uncertainty Alina Chertock, Shi Jin, and Alexander Kurganov Abstract We introduce an operator
More informationPDEs, part 3: Hyperbolic PDEs
PDEs, part 3: Hyperbolic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Hyperbolic equations (Sections 6.4 and 6.5 of Strang). Consider the model problem (the
More informationSparse Grid Collocation for Uncertainty Quantification
Sparse Grid Collocation for Uncertainty Quantification John Burkardt Department of Scientific Computing Florida State University Ajou University 14 August 2013 http://people.sc.fsu.edu/ jburkardt/presentations/......
More informationUncertainty analysis of large-scale systems using domain decomposition
Center for Turbulence Research Annual Research Briefs 2007 143 Uncertainty analysis of large-scale systems using domain decomposition By D. Ghosh, C. Farhat AND P. Avery 1. Motivation and objectives A
More informationDinesh Kumar, Mehrdad Raisee and Chris Lacor
Dinesh Kumar, Mehrdad Raisee and Chris Lacor Fluid Mechanics and Thermodynamics Research Group Vrije Universiteit Brussel, BELGIUM dkumar@vub.ac.be; m_raisee@yahoo.com; chris.lacor@vub.ac.be October, 2014
More informationAlgorithms for Uncertainty Quantification
Algorithms for Uncertainty Quantification Lecture 9: Sensitivity Analysis ST 2018 Tobias Neckel Scientific Computing in Computer Science TUM Repetition of Previous Lecture Sparse grids in Uncertainty Quantification
More informationEvaluation of Non-Intrusive Approaches for Wiener-Askey Generalized Polynomial Chaos
49th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference 6t 7 - April 28, Schaumburg, IL AIAA 28-892 Evaluation of Non-Intrusive Approaches for Wiener-Askey Generalized
More informationTong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH
Consistency & Numerical Smoothing Error Estimation An Alternative of the Lax-Richtmyer Theorem Tong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH 43403
More informationIMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1
Computational Methods in Applied Mathematics Vol. 1, No. 1(2001) 1 8 c Institute of Mathematics IMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1 P.B. BOCHEV E-mail: bochev@uta.edu
More informationMultilevel accelerated quadrature for elliptic PDEs with random diffusion. Helmut Harbrecht Mathematisches Institut Universität Basel Switzerland
Multilevel accelerated quadrature for elliptic PDEs with random diffusion Mathematisches Institut Universität Basel Switzerland Overview Computation of the Karhunen-Loéve expansion Elliptic PDE with uniformly
More informationKernel-based Approximation. Methods using MATLAB. Gregory Fasshauer. Interdisciplinary Mathematical Sciences. Michael McCourt.
SINGAPORE SHANGHAI Vol TAIPEI - Interdisciplinary Mathematical Sciences 19 Kernel-based Approximation Methods using MATLAB Gregory Fasshauer Illinois Institute of Technology, USA Michael McCourt University
More informationAN EFFICIENT COMPUTATIONAL FRAMEWORK FOR UNCERTAINTY QUANTIFICATION IN MULTISCALE SYSTEMS
AN EFFICIENT COMPUTATIONAL FRAMEWORK FOR UNCERTAINTY QUANTIFICATION IN MULTISCALE SYSTEMS A Dissertation Presented to the Faculty of the Graduate School of Cornell University in Partial Fulfillment of
More informationPolynomial chaos expansions for sensitivity analysis
c DEPARTMENT OF CIVIL, ENVIRONMENTAL AND GEOMATIC ENGINEERING CHAIR OF RISK, SAFETY & UNCERTAINTY QUANTIFICATION Polynomial chaos expansions for sensitivity analysis B. Sudret Chair of Risk, Safety & Uncertainty
More informationMultiscale methods for time-harmonic acoustic and elastic wave propagation
Multiscale methods for time-harmonic acoustic and elastic wave propagation Dietmar Gallistl (joint work with D. Brown and D. Peterseim) Institut für Angewandte und Numerische Mathematik Karlsruher Institut
More informationAn inverse scattering problem in random media
An inverse scattering problem in random media Pedro Caro Joint work with: Tapio Helin & Matti Lassas Computational and Analytic Problems in Spectral Theory June 8, 2016 Outline Introduction and motivation
More informationRelaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni
Relaxation methods and finite element schemes for the equations of visco-elastodynamics Chiara Simeoni Department of Information Engineering, Computer Science and Mathematics University of L Aquila (Italy)
More informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
More informationJournal of Computational Physics
Journal of Computational Physics 229 (2010) 1536 1557 Contents lists available at ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp Multi-element probabilistic
More informationEfficient Solvers for Stochastic Finite Element Saddle Point Problems
Efficient Solvers for Stochastic Finite Element Saddle Point Problems Catherine E. Powell c.powell@manchester.ac.uk School of Mathematics University of Manchester, UK Efficient Solvers for Stochastic Finite
More informationTHE CHOICE OF AUXILIARY DENSITY FUNCTION IN STOCHASTIC COLLOCATION
IN STOCHASTIC COLLOCATION KELLIN RUMSEY & NEVIN MARTIN 12/06/2016 1 IN STOCHASTIC COLLOCATION 2 1. Introduction Stochastic collocation has been used in a wide range of application spaces as a way to efficiently
More informationCharacterization of heterogeneous hydraulic conductivity field via Karhunen-Loève expansions and a measure-theoretic computational method
Characterization of heterogeneous hydraulic conductivity field via Karhunen-Loève expansions and a measure-theoretic computational method Jiachuan He University of Texas at Austin April 15, 2016 Jiachuan
More informationStochastic Elastic-Plastic Finite Element Method for Performance Risk Simulations
Stochastic Elastic-Plastic Finite Element Method for Performance Risk Simulations Boris Jeremić 1 Kallol Sett 2 1 University of California, Davis 2 University of Akron, Ohio ICASP Zürich, Switzerland August
More informationAdaptive methods for control problems with finite-dimensional control space
Adaptive methods for control problems with finite-dimensional control space Saheed Akindeinde and Daniel Wachsmuth Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy
More informationA reduced-order stochastic finite element analysis for structures with uncertainties
A reduced-order stochastic finite element analysis for structures with uncertainties Ji Yang 1, Béatrice Faverjon 1,2, Herwig Peters 1, icole Kessissoglou 1 1 School of Mechanical and Manufacturing Engineering,
More informationNonclassical Particle Transport in Heterogeneous Materials
Nonclassical Particle Transport in Heterogeneous Materials Thomas Camminady, Martin Frank and Edward W. Larsen Center for Computational Engineering Science, RWTH Aachen University, Schinkelstrasse 2, 5262
More informationc 2004 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 6, No., pp. 578 59 c Society for Industrial and Applied Mathematics STOCHASTIC SOLUTIONS FOR THE TWO-DIMENSIONAL ADVECTION-DIFFUSION EQUATION XIAOLIANG WAN, DONGBIN XIU, AND GEORGE
More informationRobust optimization for windmill airfoil design under variable wind conditions
Center for Turbulence Research Proceedings of the Summer Program 2012 241 Robust optimization for windmill airfoil design under variable wind conditions By H. Tachikawa, D. Schiavazzi, T. Arima AND G.
More informationEFFICIENT SHAPE OPTIMIZATION USING POLYNOMIAL CHAOS EXPANSION AND LOCAL SENSITIVITIES
9 th ASCE Specialty Conference on Probabilistic Mechanics and Structural Reliability EFFICIENT SHAPE OPTIMIZATION USING POLYNOMIAL CHAOS EXPANSION AND LOCAL SENSITIVITIES Nam H. Kim and Haoyu Wang University
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationUpscaling Wave Computations
Upscaling Wave Computations Xin Wang 2010 TRIP Annual Meeting 2 Outline 1 Motivation of Numerical Upscaling 2 Overview of Upscaling Methods 3 Future Plan 3 Wave Equations scalar variable density acoustic
More informationarxiv: v1 [math.cv] 18 Aug 2015
arxiv:508.04376v [math.cv] 8 Aug 205 Saddle-point integration of C bump functions Steven G. Johnson, MIT Applied Mathematics Created November 23, 2006; updated August 9, 205 Abstract This technical note
More informationarxiv: v2 [math.na] 8 Apr 2017
A LOW-RANK MULTIGRID METHOD FOR THE STOCHASTIC STEADY-STATE DIFFUSION PROBLEM HOWARD C. ELMAN AND TENGFEI SU arxiv:1612.05496v2 [math.na] 8 Apr 2017 Abstract. We study a multigrid method for solving large
More informationNumerical Analysis and Methods for PDE I
Numerical Analysis and Methods for PDE I A. J. Meir Department of Mathematics and Statistics Auburn University US-Africa Advanced Study Institute on Analysis, Dynamical Systems, and Mathematical Modeling
More informationA new method for the solution of scattering problems
A new method for the solution of scattering problems Thorsten Hohage, Frank Schmidt and Lin Zschiedrich Konrad-Zuse-Zentrum Berlin, hohage@zibde * after February 22: University of Göttingen Abstract We
More informationLow-rank techniques applied to moment equations for the stochastic Darcy problem with lognormal permeability
Low-rank techniques applied to moment equations for the stochastic arcy problem with lognormal permeability Francesca Bonizzoni 1,2 and Fabio Nobile 1 1 CSQI-MATHICSE, EPFL, Switzerland 2 MOX, Politecnico
More informationSlow Growth for Gauss Legendre Sparse Grids
Slow Growth for Gauss Legendre Sparse Grids John Burkardt, Clayton Webster April 4, 2014 Abstract A sparse grid for multidimensional quadrature can be constructed from products of 1D rules. For multidimensional
More information