Numerical Solution for Random Forced SPDE via Galerkin Finite Element Method
|
|
- Gerard Ezra Henry
- 6 years ago
- Views:
Transcription
1 Journal of mathematics and computer science 9 (014), 71-8 umerical Solution for Random Forced SE via Galerkin Finite Element ethod Rasoul aseri, A. alek 1 epartment of Applied athematics, Faculty of athematical Sciences, Tarbiat odares University, Article history: Received October 013 Accepted ecember 013 Available online ecember 013. O. Box , Tehran, Iran (naserir@modares.ac.ir& mala@modares.ac.ir) Abstract In this paper, we present a deterministic finite element approach for solving a random forced iffusion equation. Separation of random and deterministic variables is done by Karhunen- Loeve expansion. Truncating the Karhunen-Loeve expansion of the permeability field leads to a finite dimensional approximation of the problem. The problem is discretized, in spatial part, using the finite-element method and the polynomial chaos expansion in stochastic part. Finally, using Kronecker product preconditioner and thus, preconditioned conjugate gradient method the governed system of equation is solved. umerical experiments are presented for illustrating the theoretical results. Keywords: Stochastic partial differential equation, Karhunen-Loeve expansion, Wiener Chaos expansion, finite element method. 1. Introduction hysical problems, often, can be formulated as mathematical models which in many cases contain partial differential equations (Es). Uncertainty might plagues everything from modeling assumptions to experimental data. So, the differential coefficient and source functions in the model might represent as functions of the spatial domain and some sample space. In this case we deal with a stochastic partial differential equation (SE). In this paper, we consider the stochastic steady-state diffusion equation along with homogeneous irichlet boundary-value conditions. In stochastic engineering, the perturbation method is a popular technique, cf., e.g., [7, 10, 11, 13, 14, 16, 1]. Considering certain smoothness conditions, the random functions and operators involved in the differential equation are expanded in a Taylor series about their respective mean values. In this way, usually good results are obtained only for small deviations, cf., e.g., [11,16]. The eumann expansion series method is another approach. In this 1 Corresponding author. address: mala@modares.ac.ir.
2 R. aseri, A. alek / J. ath. Computer Sci. ( ), - approach, the inverse of the stochastic operator is approximated by its eumann series [4,18]. Recently, a method based on a spectral representation of the uncertainty is introduced by Ghanem and Spanos [9], where utilizes the Karhunen-Loeve expansion (KLE) of correlated random functions, cf. [1]. In order to allow for other representations of the random processes, the spectral finite element scheme is generalized by Babuska and et al. [1,,3]. Here, we follow the Galerkin finite element method of stochastic diffusion problems. We use the kronecker product preconditioner presented by Ullmann [0] to clustering the eigenvalues of the global matrix gained at the end of variational computation. Thus, using preconditioned conjugate gradient method [17], the large scale system of equation is solved, which gives good computational results. The organization of this paper is as follow: In Section, Karhunen-Loeve expansion (KLE) and polynomial chaos expansion (CE) is presented. In Section 3, stochastic Galerkin method, kroneker product preconditioner and computation of statistics is proposed. In Section 4, some numerical result is illustrated. roblem formulation Our model problem is a steady state diffusion problem in a domain with inhomogeneous stochastic diffusion coefficient a x, ω :. a x, ω u x, ω = z x, ω in u x, ω = 0 1 where is the spatial domain, is probability space and a x, ω is the correlated random fields. Hence, the solutions u x, ω of the SE (1) are also random fields. As an important assumption for the stochastic diffusion equation (1), it is assumed that the random coefficient u x, ω satisfies the elliptic condition. That is, there exist a constant a min such that 0 < a min a x, ω ().1 Karhunen-Loeve Expansion (KLE) Consider a random field a(x, ω), x, with finite second order moment E[ a (x, ω)]dx < Assume that E a = a x. It is possible to expand a(x, ω), for a given orthonormal basis {ψ k } in L (), as a generalized Fourier series where a x, ω = a x + a k ω = a x, ω a k ω ψ k x dx, ψ k x (3) k = 1,, are random variables with zero means. It is important, now, to find a special basis {φ k } that makes corresponding a k uncorrelated: E[a i a j ] = 0 for al li j. enoting the covariance function of a x, ω by R x, y = E[a x, ω a y, ω ], the basis functions {φ k } should satisfy 7
3 R. aseri, A. alek / J. ath. Computer Sci. ( ), - E a i a j = φ i (x) dx R x, y φ j y dy = 0, i j Complete and orthonormality of {φ k } in L () follows that φ k (x) are eigenfunctions of R x, y : where λ j = E a j R x, y φ j y dy = λ j φ j x, j = 1,, (4) > 0. Indeed, by choosing basis functions φ k (x) as the solutions of the eigenproblem (4), the random variables a k ω will be uncorrelated. enoting θ k = a k / λ k, we have the following expansion: a x, ω = a x + λ k θ k ω φ k x (5) where θ k satisfy E[θ k ] = 0 and E[θ i θ j ] = δ ij. In the case that a x, ω is considered as a Gaussian process, a k ω, k = 1,, will be independent Gaussian random variables. The expansion (5) is known as the Karhunen-Loeve expansion (KLE) of the stochastic process a x, ω. Using the KLE (5), the stochastic process can be represented as a series of uncorrelated random variables. Since the basic functions φ k (x) are deterministic, the spatial dependence of the random process can be resolved by them. The KLE converges to the random process a x, ω in the mean square sense where lim E a x, ω a x, ω dx = 0, a = a x + λ k θ k φ k (6) is a finite term [1,15]. So, among all the Fourier expansions (3), the KLE with the basis function φ k (x) determined by the eigen problem (4) has the fastest convergence rate. The convergence rate of the KLE only depends on the smoothness of the covariance function, but not on the specific probability structure of the process being expanded. The KLE converges fast if the covariance function R x, y is very smooth. In this case, we only need to keep the leading order terms (quantified by the magnitude of λ k ) in the finite KLE and still capture most of the energy of the stochastic process a x, ω. For the finite KLE (6), we define its energy ratio as e = E a dx = E a dx λ k λ k If the eigenvalues λ k, k = 1,,, decay very fast, then the finite term KLE would be good approximation of the stochastic process [5,8].. olynomial Chaos Expansion (CE) There are problems, as the solution of a E with random inputs, that the covariance function of a random process u x, ω, x is not known. The solution of such problems can be represented using a polynomial chaos expansion (CE) given by u x, ω = u k x Ψ k ω (7) 73
4 R. aseri, A. alek / J. ath. Computer Sci. ( ), - where, the functions u k x are deterministic coefficients, ω is a vector of orthonormal random variables and Ψ k ω are multi-dimensional orthogonal polynomials with the following properties: Ψ 1 E Ψ 1 = 1, Ψ k = 0, k > 1, Ψ i Ψ j = h i δ ij (8) The convergence property of CE for a random quantity in L is ensured by the Cameron-artin theorem [5, ], i.e. u x, ω u k x Ψ k ω L 0. Hence, for the sake of computation, this convergence justifies a truncation of CE to a finite number of terms, u x, ω = u k x Ψ k ω (9) where, the value of is determined by the highest degree of polynomial d, used to represent u, and the number of random variables the length of ω with the formula + 1 = ( + d)!! d! []. Generally, the value of is the same as the number of uncorrelated random variables in the system or equivalently, the truncation length of the truncated KLE. Typically, the value of d is chosen by some heuristic method. Indeed, in the case of d = 1 and random variables, the KLE is a special case of the CE. 3 Stochastic Galerkin Suppose W i L ρi i with dimension ρ i for i = 1,, and V H 1 0 () with dimension. In addition; let {ψ i p n } i n=1 for i = 1,, be basis of Wi and {φ i } i=1 a basis of V. The finite dimensional tensor product space W 1 W V can be defined as the space spanned by the functions {ψ 1 n1,, ψ n φ i } for all n 1 {1,, p 1 },,n {1,, p } and i {1,, }. To simplify notation, let n denote a multiindex whose kth component n k {1,, p k } and I denote the set of such multi-indices, then the basis functions for the tensor product space have the form v ni x, ω = φ i (x)ψ n ω where Ψ n = k ψ nk ω k. Then, the Galerkin method is looking for a solution u W 1 W V such that ρ ω a x, ω u h x, ω. v ni x, ω dxdω = ρ ω z x, ω v ni x, ω dxdω for all n I and i = 1,,, which ρ(ω) is the integration weight. Equivalently,plugging in the explicit formula for the basic functions, equation (10) can be written as ρ ω Ψ n (ω) a x, ω u h x, ω. φ i x dxdω = ρ(ω)ψ n (ω) z x, ω φ i x dxdω (11) for all n I and i = 1,,. ow, considering the approximate solution as (10) 74
5 R. aseri, A. alek / J. ath. Computer Sci. ( ), - u h x, ω = and substituting (1) into the equation (11) yields u jm j =1 m I j =1 m I u jm φ j x Ψ m ω (1) ρ ω Ψ n ω Ψ m ω a x, ω φ j x. φ i x dxdω (13) = ρ ω Ψ n ω z x, ω φ i x dxdω for all m I and i = 1,,. Let = I = p k, then, a natural bijection between {1,, } and I can be defined. Thus, equation (13) can be written as u jm j =1 m =1 ρ(ω)ψ n ω Ψ m ω a x, ω φ j x. φ i x dxdω 14 = ρ(ω)ψ n (ω) z(x, ω)φ i x dxdω ow, if a and z have the following KLE: and Then, we can rewrite equation (14) as where u jm K 0 i,j j =1 m =1 a x, ω = a 0 x + y i a i (x) i=1 z x, ω = z 0 x + y i z i (x) ρ ω Ψ n ω Ψ m ω dω + i=1 K k i,j m =1 = (z 0 ) i ρ(ω)ψ n (ω)dω + (z k ) i ρ(ω)ψ n (ω)y k dω m=1 ρ ω Ψ n ω Ψ m ω y k dω 15 (K k ) i,j = a k x φ i x. φ j x dx w r a k x r φ i x r. φ j x r (16) and r z k i = z i x φ i x dx w r z i x r φ i x r for k = 1,, and i = 1,, [3]. Equations (16) and (17) are approximated by quadrature rule in spatial domain. ow, assume that there exist functions Ψ n, n = 1,, such that r (17) 75
6 R. aseri, A. alek / J. ath. Computer Sci. ( ), - and ρ(ω)ψ n ω Ψ m ω dω = δ mn ρ(ω)ψ n ω Ψ m ω y k dω = C kn δ mn. where, for m = n, the above integral is approximated using quadrature rule in random domain as follow: Then, equation (15) becomes or, equivalently, u mj j =1 u jm j =1 m =1 ρ(ω)ψ n ω Ψ n ω y k dω w s ρ ω s Ψ n ω s y k ω s (18) s (K 0 ) i,j + C kn (K k ) i,j δ mn = (z 0 ) i ρ(ω)ψ n (ω)dω + z k i ρ ω Ψ n ω y k dω (19) (K 0 ) i,j + C kn (K k ) i,j = (z 0 ) i ρ(ω)ψ n (ω)dω + (z k ) i ρ ω Ψ n ω y k dω (0) In equations (19) and (0) we have: and ρ(ω)ψ n (ω)dω = 0; for n 0, ρ(ω)ψ n (ω)y k dω = ρ(ω)ψ 0 (ω)ψ n (ω)y k dω = C kn δ 0n = C k0. Finally, equation (0) can be considered as the following block diagonal system: K 0 + C k1 K k 0 0 K 0 + C kp K k u 1 u = z 1 z p (1) or briefly, Ku = z. Once u has been computed, it can be post-processed to obtain meaningful information, such as the mean and variance of u h. efinitions of mean and variance of a random field (noting that a random field is a stochastic process) can be found in [14,19]. The expected value of u h is given by E(u h ) = ρ ω u jm j =1 m =1 φ j x Ψ m ω dω 76
7 R. aseri, A. alek / J. ath. Computer Sci. ( ), - = u jm φ j x ρ ω Ψ m ω dω j =1 m =1 = u j1 φ j x j =1 on account that Ψ 1 = 1 and the orthogonality conditions (8). Therefore, if x j is a node corresponding to the basis function φ j, E(u h (x j, ω)) = u j1, the variance of u h is given by V(u h ) = ρ(ω) u jm φ j x Ψ m ω dω E(u h ) Adopting the notation j =1 m =1 = u jm φ j x ρ ω Ψ m ω dω j =1 m =1 j =1 u k x = (u k ) j φ j x, one can compute the m th statistical moments of the solution as j =1 u j1 φ j x E u., x m E Ψ k. u k x m = E Ψ k1. Ψ km u k1 x u km x. 3.1 Kronecker product preconditioners. Rewriting equation (1) as: k=0 G k K k K 77 u = z the strategy is to choose = G K 0 as a preconditioner, such that G = argmin{h R : K H K 0 F } where, as mentioned in CE, is equals with the degree of polynomial chaos truncation, and. F denotes the Frobenius norm. The closed form of the solution can be written as follow [6,0], tr(k T k K 0 ) G = I + tr(k T G k. 0 K 0 ) q Which tr(k T k K 0 ) = [K k ] i,i [K 0 ] i=1 i,i and hence, the coefficients in above equality can be computed straightforward. In addition, since K and K 0 are symmetric and positive definite, so G and = G K 0 have also these properties. ow, using preconditioned ewton s conjugate gradient algorithm, we can find fast and accurate convergence of the solution [0]. In fact can be written as:
8 R. aseri, A. alek / J. ath. Computer Sci. ( ), - = tr(k T k K 0 ) 1 + C k1 tr(k T 0 K 0 ) So 1 K can be expressed as: K K = 0 tr(k T k K 0 ) C kp tr(k T 0 K 0 ) C k1 K k tr (K 1 + C T k K 0 ) k1 tr (K 0 T K 0 ) 0 K 0 + C kp K k tr(k 1 + C T k K 0 ) k1 0 0 K 0 + C kp K k tr (K 1 + C T k K 0 ) kp tr(k 0 T K 0 ) tr (K 0 T K 0 ) Finally, after preconditioning the problem, we are interested in solving the following equation: K 0 + C k1 K k z 1 0 tr (K 1 + C T k K 0 ) tr (K k1 tr (K T 0 K 0 ) u 1 + C T k K 0 ) k1 1 tr (K T 0 K 0 ) =. u z tr (K 1 + C T k K 0 ) k1 tr (K 0 T K 0 ) Using ewton s conjugate gradient algorithm we can find fast and accurate convergence of the solution. 4 Experimental results Example 1- oisson Equation: The first example is a oisson equation in two spatial dimensions with one constant random parameter. Let = { x, y : 1 x 1, 1 y 1} be the spatial domain of the problem and the sample space. Then we seek a solution u x, y, ω such that satisfies a ω u x, y, ω xx + u x, y, ω yy = ex p 64 x 1 + y 1 in u x, y, ω = 0 on where a ω is a random variable distributed uniformly over the interval [1,3]. First we express a in terms of a uniform random variable on [ 1,1], i.e. a = ξ + where ξ~u 1,1. Figure 1, illustrates approximated solution, standard deviation and expectation counter plot and global matrix sparsity, for 50*50 mesh points of spatial domain, d=3, (number of random variables) and =4 (degree of polynomials). 78
9 R. aseri, A. alek / J. ath. Computer Sci. ( ), - Fig 1.Illustration of the approximated solution and its mean value, standard deviation and global matrix sparsity for d=3, =4, 50*50 mesh points in Example 1. Example : consider the following SE. a x, y, ω u x, y, ω = z x, y x, y ε, ωε where = [ 1,1], ~U[ 1,1], the boundary condition is given as u = 1.1sign sin πx sin πy. and z x, y = sin πx sin πy + x. The random field a x, y, ω is characterized by its mean and covariance function E a = 10, R x 1, x = e (x 1,x ), x 1, x 1,1. The truncated KLE of a x, y, ω can be expressed as: a x, y, ω = 10 + λ j ω j φ j x. j =0 The eigenpairs λ j, ω j in truncated KLE solve the integral equation e (x 1,x ) φ j x d, x = λ j φ j, x 1 For this special case of the covariance function, we have explicit expression for λ j and φ j [9]. Let ω jeven and ω j odd solve the equations 1 ω jeven tan ω jeven = 0, ω j odd + tan ω j odd = 0 79
10 R. aseri, A. alek / J. ath. Computer Sci. ( ), - Then the even and odd indexed eigenfunctions are given by: φ jeven x = cos (ω j even x) φ j 1 + sin (ω jeven ) odd x = ω j even sin (ω j odd x) 1 sin (ω j odd ) ω jodd and corresponding eigenvalues are given by: λ jeven x = + 1 λ j odd x = ω jeven ω j odd + 1 We choose ω j = (ω 1,, ω ) T to be independent random variables uniformly distributed over the interval 1,1. First of all, we find numerical approximations of ω jeven and ω j odd with bisection method, and then, with the eigenpairs evaluated with these ω jeven and ω j odd, we construct the KLE of d =, = 6, which emphasis that = 8. Figure, illustrates mean value representation, standard deviation and expectation counter plot and global matrix sparsity, for 0*0 mesh points Fig. Illustration of the approximated solution and its mean value, standard deviation and global matrix sparsity for d=, =6, 0*0 mesh points in Example. Example 3: similar to Example, consider the following SE. a x, y, ω u x, y, ω = z x, y x, y ε, ωε where = [ 1,1], ~U[ 1,1], with the boundary condition u x, y, ω = 0 on and z x, y = π sin πx sin πy in 80
11 R. aseri, A. alek / J. ath. Computer Sci. ( ), - Figure 3, illustrates mean value representation, standard deviation and expectation counter plot and global matrix sparsity, for 50*50 mesh points of spatial domain, d=1, (number of random variables) and =3 (degree of polynomials) which emphasis that =4. Fig 3. umerical solution of example 3 with 50*50 mesh points, d=1, =3 5 Conclusion We described a deterministic finite element (FE) solution algorithm for the stochastic elliptic boundary value problem, whose coefficient are assumed to be random field. Separation of random and deterministic variables is achieved via Karhunen-Loeve Expansion (KLE). Finite element discretization utilized in both spatial space and probability space. Using the olynomial Chaos Expansion (CE) and constructing the solution based on orthogonal basis and truncating the KLE of the permeability field lead to find a finite dimensional approximation. Finally, since the resulted system of equation governed by discretizing the problem is large and sparse, we use Kronecker product preconditioner and ewton s conjugate gradient method to solve this system of equation more accurate and fast. umerical results presented for illustrating the theoretical results. References [1] I. Babuska, F. obile, and R. Tempone, A stochastic collocation method for elliptic partial differential equations with random input data, SIA J. umer. Anal., 45(3): (electronic), 007. [] I. Babuska, R. Tempone, and G. E. Zouraris, Solving elliptic boundary value problems with uncertain coefficients by the finite element method: the stochastic formulation, Comput. ethods Appl. ech. Engrg., 194 (1-16): , 005. [3] I. Babuska, R. Tempone, and G. E. Zouraris, Galerkin finite element approximations of stochastic elliptic partial differential equations, SIA J. umer. Anal., 4(): (electronic), 004. [4] A. T. Bharrucha-Reid, On random operator equations in Banach space, Bull, Acad. olon. 81
12 R. aseri, A. alek / J. ath. Computer Sci. ( ), - Sci., Ser. Sci. ath. Astr. hys. 7 (1959) [5] H. Elman,. Furnival, Solving the stochastic steady-state diffusion problem using multigrid, IA Journal of umerical Analysis 7, , (007). [6] O. G. Ernst and E. Ullmann, Stochastic Galerkin matrices, SIA J. atrix Anal. Appl., 31(4): , 010. [7] B. Fathi Vajargah, R. Asghari, Simulation ethod for Solving Stochastic differential Equations with Constant iffusion Coefficients, Journal of mathematics and computer Science 8 (014) 8-3. [8] R. G. Ghanem, R.. Kruger, umerical solution of spectral stochastic finite element systems, Comput. ethods Appl. ech. Engrg , [9] R. Ghanem,.. Spanos, Stochastic Finite Elements: A Spectral Approach, Springer, Berlin, [10]. heydari,. ahdavi azdeh and. Bayat, Stochastic Online Scheduling With reemption enalties, Journal of mathematics and computer Science 6 (013), [11] T. Hisada, S. akagiri, Role of the stochastic finite element method in structural safety and reliability, roc. 4th Int. Conf. on Struct. Safety and Reliability, Kobe, Japan, vol. I, , 1985,. [1] K. Karhunen, Uber lineare ethoden in der Wahrscheinlichkeitsrechnung, Ann. Acad. Sci. Fennicae. Ser. A. I. ath.-hys., 1947(37):79, [13]. Kleiber, T.. Hien, The Stochastic Finite Element ethod, Wiley, ew York, 199. [14] W. K. Liu, T. Belytschko, A. ani, Random field finite elements, Int. J. um. eth. Engrg , [15]. Lo eve, Fonctions aleatoires de second ordre, Revue Sci., 84:195-06, [16] S. akagiri, T. Hisada, Stochastic finite element method applied to structural analysis with uncertain parameters, roc. Intl. Conference on FE, August pp , 198. [17] J. ocedal and S. J. Wright, umerical Optimization, Springer Verlag, Berlin, Heidelberg, ew York, second edition, 006. [18]. Shinozuka, E. Lenoe, A probabilistic model for spatial distribution of material properties, Eng. Fract. ech , [19].. Spanos, R. Ghanem, Stochastic finite element expansion for random media, Eng. ech., ASCE 115, [0] E. Ullmann, A Kronecker product preconditioner for stochastic Galerkin finite element discretizations, 3():93-946, 010. [1] E.H. Vanmarcke,. Grigoriu, Stochastic finite element analysis of simple beams, J. Eng. ech., ASCE 109 (1983) []. Wiener, The homogeneous chaos, Amer. J. ath., 60: ,
Efficient 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 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 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 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 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 informationBlock-diagonal preconditioning for spectral stochastic finite-element systems
IMA Journal of Numerical Analysis (2009) 29, 350 375 doi:0.093/imanum/drn04 Advance Access publication on April 4, 2008 Block-diagonal preconditioning for spectral stochastic finite-element systems CATHERINE
More informationSTOCHASTIC FINITE ELEMENTS WITH MULTIPLE RANDOM NON-GAUSSIAN PROPERTIES
STOCHASTIC FIITE ELEMETS WITH MULTIPLE RADOM O-GAUSSIA PROPERTIES By Roger Ghanem, 1 Member, ASCE ABSTRACT: The spectral formulation of the stochastic finite-element method is applied to the problem of
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 informationResearch Article Numerical Optimal Control for Problems with Random Forced SPDE Constraints
ISR Applied athematics Volume 24, Article I 97435, pages http://d.doi.org/.55/24/97435 Research Article umerical Optimal Control for roblems with Random Forced SE Constraints R. aseri and A. alek epartment
More informationHierarchical Parallel Solution of Stochastic Systems
Hierarchical Parallel Solution of Stochastic Systems Second M.I.T. Conference on Computational Fluid and Solid Mechanics Contents: Simple Model of Stochastic Flow Stochastic Galerkin Scheme Resulting Equations
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 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 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 informationSchwarz Preconditioner for the Stochastic Finite Element Method
Schwarz Preconditioner for the Stochastic Finite Element Method Waad Subber 1 and Sébastien Loisel 2 Preprint submitted to DD22 conference 1 Introduction The intrusive polynomial chaos approach for uncertainty
More informationKarhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques
Institut für Numerische Mathematik und Optimierung Karhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques Oliver Ernst Computational Methods with Applications Harrachov, CR,
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 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 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 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 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 informationUNCERTAINTY ASSESSMENT USING STOCHASTIC REDUCED BASIS METHOD FOR FLOW IN POROUS MEDIA
UNCERTAINTY ASSESSMENT USING STOCHASTIC REDUCED BASIS METHOD FOR FLOW IN POROUS MEDIA A REPORT SUBMITTED TO THE DEPARTMENT OF ENERGY RESOURCES ENGINEERING OF STANFORD UNIVERSITY IN PARTIAL FULFILLMENT
More informationPARALLEL COMPUTATION OF 3D WAVE PROPAGATION BY SPECTRAL STOCHASTIC FINITE ELEMENT METHOD
13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 24 Paper No. 569 PARALLEL COMPUTATION OF 3D WAVE PROPAGATION BY SPECTRAL STOCHASTIC FINITE ELEMENT METHOD Riki Honda
More informationRandom Fields in Bayesian Inference: Effects of the Random Field Discretization
Random Fields in Bayesian Inference: Effects of the Random Field Discretization Felipe Uribe a, Iason Papaioannou a, Wolfgang Betz a, Elisabeth Ullmann b, Daniel Straub a a Engineering Risk Analysis Group,
More informationSpectral methods for fuzzy structural dynamics: modal vs direct approach
Spectral methods for fuzzy structural dynamics: modal vs direct approach S Adhikari Zienkiewicz Centre for Computational Engineering, College of Engineering, Swansea University, Wales, UK IUTAM Symposium
More informationCollocation based high dimensional model representation for stochastic partial differential equations
Collocation based high dimensional model representation for stochastic partial differential equations S Adhikari 1 1 Swansea University, UK ECCM 2010: IV European Conference on Computational Mechanics,
More informationKarhunen-Loève Expansions of Lévy Processes
Karhunen-Loève Expansions of Lévy Processes Daniel Hackmann June 2016, Barcelona supported by the Austrian Science Fund (FWF), Project F5509-N26 Daniel Hackmann (JKU Linz) KLE s of Lévy Processes 0 / 40
More informationMULTISCALE FINITE ELEMENT METHODS FOR STOCHASTIC POROUS MEDIA FLOW EQUATIONS AND APPLICATION TO UNCERTAINTY QUANTIFICATION
MULTISCALE FINITE ELEMENT METHODS FOR STOCHASTIC POROUS MEDIA FLOW EQUATIONS AND APPLICATION TO UNCERTAINTY QUANTIFICATION P. DOSTERT, Y. EFENDIEV, AND T.Y. HOU Abstract. In this paper, we study multiscale
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 informationarxiv: v1 [math.na] 18 Dec 2017
A Bramble-Pasciak conjugate gradient method for discrete Stokes problems with lognormal random viscosity arxiv:1712.06472v1 [math.na] 18 Dec 2017 Christopher Müller 1,2, Sebastian Ullmann 1,2, and Jens
More informationFast Algorithm for Computing Karhunen-Loève Expansion
Fast Algorithm for Computing Karhunen-Loève Expansion Qian-Yong Chen Abstract Karhunen Loève expansion has been proven to be an efficient way to approximate second-order random processes. But numerical
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 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 informationA Vector-Space Approach for Stochastic Finite Element Analysis
A Vector-Space Approach for Stochastic Finite Element Analysis S Adhikari 1 1 Swansea University, UK CST2010: Valencia, Spain Adhikari (Swansea) Vector-Space Approach for SFEM 14-17 September, 2010 1 /
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 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 informationNumerical methods for the discretization of random fields by means of the Karhunen Loève expansion
Numerical methods for the discretization of random fields by means of the Karhunen Loève expansion Wolfgang Betz, Iason Papaioannou, Daniel Straub Engineering Risk Analysis Group, Technische Universität
More informationA Preconditioned Recycling GMRES Solver for Stochastic Helmholtz Problems
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 6, No. 2, pp. 342-353 Commun. Comput. Phys. August 2009 A Preconditioned Recycling GMRES Solver for Stochastic Helmholtz Problems Chao Jin 1 and Xiao-Chuan
More informationarxiv: v2 [math.pr] 27 Oct 2015
A brief note on the Karhunen-Loève expansion Alen Alexanderian arxiv:1509.07526v2 [math.pr] 27 Oct 2015 October 28, 2015 Abstract We provide a detailed derivation of the Karhunen Loève expansion of a stochastic
More informationLecture 1: Center for Uncertainty Quantification. Alexander Litvinenko. Computation of Karhunen-Loeve Expansion:
tifica Lecture 1: Computation of Karhunen-Loeve Expansion: Alexander Litvinenko http://sri-uq.kaust.edu.sa/ Stochastic PDEs We consider div(κ(x, ω) u) = f (x, ω) in G, u = 0 on G, with stochastic coefficients
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 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 informationNotes Regarding the Karhunen Loève Expansion
Notes Regarding the Karhunen Loève Expansion 1 Properties of the Karhunen Loève Expansion As noted in Section 5.3 of [3], the Karhunen Loève expansion for a correlated second-order random field α(x, ω),
More informationNumerical Approximation of Stochastic Elliptic Partial Differential Equations
Numerical Approximation of Stochastic Elliptic Partial Differential Equations Hermann G. Matthies, Andreas Keese Institut für Wissenschaftliches Rechnen Technische Universität Braunschweig wire@tu-bs.de
More informationEfficient Observation of Random Phenomena
Lecture 9 Efficient Observation of Random Phenomena Tokyo Polytechnic University The 1st Century Center of Excellence Program Yukio Tamura POD Proper Orthogonal Decomposition Stochastic Representation
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 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 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 informationThe Conjugate Gradient Method
The Conjugate Gradient Method Classical Iterations We have a problem, We assume that the matrix comes from a discretization of a PDE. The best and most popular model problem is, The matrix will be as large
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 informationPolynomial Chaos and Karhunen-Loeve Expansion
Polynomial Chaos and Karhunen-Loeve Expansion 1) Random Variables Consider a system that is modeled by R = M(x, t, X) where X is a random variable. We are interested in determining the probability of the
More informationRational Chebyshev pseudospectral method for long-short wave equations
Journal of Physics: Conference Series PAPER OPE ACCESS Rational Chebyshev pseudospectral method for long-short wave equations To cite this article: Zeting Liu and Shujuan Lv 07 J. Phys.: Conf. Ser. 84
More informationSpectral Polynomial Chaos Solutions of the Stochastic Advection Equation
Spectral Polynomial Chaos Solutions of the Stochastic Advection Equation M. Jardak, C.-H. Su and G.E. Karniadakis Division of Applied Mathematics Brown University October 29, 21 Abstract We present a new
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 informationA Spectral Approach to Linear Bayesian Updating
A Spectral Approach to Linear Bayesian Updating Oliver Pajonk 1,2, Bojana V. Rosic 1, Alexander Litvinenko 1, and Hermann G. Matthies 1 1 Institute of Scientific Computing, TU Braunschweig, Germany 2 SPT
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 informationLecture 12: Detailed balance and Eigenfunction methods
Miranda Holmes-Cerfon Applied Stochastic Analysis, Spring 2015 Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility),
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 informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationQuantifying Uncertainty: Modern Computational Representation of Probability and Applications
Quantifying Uncertainty: Modern Computational Representation of Probability and Applications Hermann G. Matthies with Andreas Keese Technische Universität Braunschweig wire@tu-bs.de http://www.wire.tu-bs.de
More informationLecture 12: Detailed balance and Eigenfunction methods
Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility), 4.2-4.4 (explicit examples of eigenfunction methods) Gardiner
More informationA Stochastic Collocation based. for Data Assimilation
A Stochastic Collocation based Kalman Filter (SCKF) for Data Assimilation Lingzao Zeng and Dongxiao Zhang University of Southern California August 11, 2009 Los Angeles Outline Introduction SCKF Algorithm
More informationOn the Application of Intervening Variables for Stochastic Finite Element Analysis
On the Application of Intervening Variables for Stochastic Finite Element Analysis M.A. Valdebenito a,, A.A. Labarca a, H.A. Jensen a a Universidad Tecnica Federico Santa Maria, Dept. de Obras Civiles,
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 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 informationMicrostructurally-Informed Random Field Description: Case Study on Chaotic Masonry
Microstructurally-Informed Random Field Description: Case Study on Chaotic Masonry M. Lombardo 1 J. Zeman 2 M. Šejnoha 2,3 1 Civil and Building Engineering Loughborough University 2 Department of Mechanics
More informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
More informationH(div) Preconditioning for a Mixed Finite Element Formulation of the Stochastic Diffusion Problem 1
University of Maryland Department of Computer Science CS-TR-4918 University of Maryland Institute for Advanced Computer Studies UMIACS-TR-2008-15 H(div) Preconditioning for a Mixed Finite Element Formulation
More informationSENSITIVITY ANALYSIS IN NUMERICAL SIMULATION OF MULTIPHASE FLOW FOR CO 2 STORAGE IN SALINE AQUIFERS USING THE PROBABILISTIC COLLOCATION APPROACH
XIX International Conference on Water Resources CMWR 2012 University of Illinois at Urbana-Champaign June 17-22,2012 SENSITIVITY ANALYSIS IN NUMERICAL SIMULATION OF MULTIPHASE FLOW FOR CO 2 STORAGE IN
More informationRandom Vibrations & Failure Analysis Sayan Gupta Indian Institute of Technology Madras
Random Vibrations & Failure Analysis Sayan Gupta Indian Institute of Technology Madras Lecture 1: Introduction Course Objectives: The focus of this course is on gaining understanding on how to make an
More informationUncertainty quantification for flow in highly heterogeneous porous media
695 Uncertainty quantification for flow in highly heterogeneous porous media D. Xiu and D.M. Tartakovsky a a Theoretical Division, Los Alamos National Laboratory, Mathematical Modeling and Analysis Group
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 informationAnalysis and finite element approximations of stochastic optimal control problems constrained by stochastic elliptic partial differential equations
Retrospective Theses and issertations Iowa State University Capstones, Theses and issertations 2008 Analysis and finite element approximations of stochastic optimal control problems constrained by stochastic
More informationDynamic Data Driven Simulations in Stochastic Environments
Computing 77, 321 333 (26) Digital Object Identifier (DOI) 1.17/s67-6-165-3 Dynamic Data Driven Simulations in Stochastic Environments C. Douglas, Lexington, Y. Efendiev, R. Ewing, College Station, V.
More informationUncertainty Evolution In Stochastic Dynamic Models Using Polynomial Chaos
Noname manuscript No. (will be inserted by the editor) Uncertainty Evolution In Stochastic Dynamic Models Using Polynomial Chaos Umamaheswara Konda Puneet Singla Tarunraj Singh Peter Scott Received: date
More informationSampling and low-rank tensor approximation of the response surface
Sampling and low-rank tensor approximation of the response surface tifica Alexander Litvinenko 1,2 (joint work with Hermann G. Matthies 3 ) 1 Group of Raul Tempone, SRI UQ, and 2 Group of David Keyes,
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 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 informationAnalysis and Computation of Hyperbolic PDEs with Random Data
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,
More informationNovember 7, 2002 GALERKIN FINITE ELEMENT APPROXIMATIONS OF STOCHASTIC ELLIPTIC PARTIAL DIFFERENTIAL EQUATIONS
November 7, 00 GALERKIN FINITE ELEMENT APPROXIMATIONS OF STOCHASTIC ELLIPTIC PARTIAL IFFERENTIAL EQUATIONS IVO BABUŠKA, RAÚL TEMPONE AN GEORGIOS E. ZOURARIS Abstract. We describe and analyze two numerical
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 informationAND ELISABETH ULLMANN
Proceedings of ALGORITMY 25 pp. 1 1 COMPUTATIONAL ASPECTS OF THE STOCHASTIC FINITE ELEMENT METHO MICHAEL EIERMANN, OLIVER G. ERNST AN ELISABETH ULLMANN Abstract. We present an overview of the stochastic
More informationStochastic Design Criteria in Linear Models
AUSTRIAN JOURNAL OF STATISTICS Volume 34 (2005), Number 2, 211 223 Stochastic Design Criteria in Linear Models Alexander Zaigraev N. Copernicus University, Toruń, Poland Abstract: Within the framework
More informationThe Matrix Dyson Equation in random matrix theory
The Matrix Dyson Equation in random matrix theory László Erdős IST, Austria Mathematical Physics seminar University of Bristol, Feb 3, 207 Joint work with O. Ajanki, T. Krüger Partially supported by ERC
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 34: Improving the Condition Number of the Interpolation Matrix Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu
More informationStochastic structural dynamic analysis with random damping parameters
Stochastic structural dynamic analysis with random damping parameters K. Sepahvand 1, F. Saati Khosroshahi, C. A. Geweth and S. Marburg Chair of Vibroacoustics of Vehicles and Machines Department of Mechanical
More informationNumerical Analysis of Differential Equations Numerical Solution of Elliptic Boundary Value
Numerical Analysis of Differential Equations 188 5 Numerical Solution of Elliptic Boundary Value Problems 5 Numerical Solution of Elliptic Boundary Value Problems TU Bergakademie Freiberg, SS 2012 Numerical
More informationA Data-Driven Stochastic Method
A ata-riven Stochastic Method Mulin Cheng, Thomas Y. Hou, Pengchong Yan Applied and Computational Mathematics, California Institute of Technology, Pasadena, CA 96 Abstract Stochastic partial differential
More informationA Review of Preconditioning Techniques for Steady Incompressible Flow
Zeist 2009 p. 1/43 A Review of Preconditioning Techniques for Steady Incompressible Flow David Silvester School of Mathematics University of Manchester Zeist 2009 p. 2/43 PDEs Review : 1984 2005 Update
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 informationStochastic Finite Element Modeling and Its Applications in Heat Conduction and Elastostatics Problems
Stochastic Finite Element Modeling and Its Applications in Heat Conduction and Elastostatics Problems RICHARD J. AQUINO A master s dissertation submitted to the University of Wales in fulfillment of the
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 informationMultigrid and stochastic sparse-grids techniques for PDE control problems with random coefficients
Multigrid and stochastic sparse-grids techniques for PDE control problems with random coefficients Università degli Studi del Sannio Dipartimento e Facoltà di Ingegneria, Benevento, Italia Random fields
More informationStochastic Processes. M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno
Stochastic Processes M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno 1 Outline Stochastic (random) processes. Autocorrelation. Crosscorrelation. Spectral density function.
More informationDynamic response of structures with uncertain properties
Dynamic response of structures with uncertain properties S. Adhikari 1 1 Chair of Aerospace Engineering, College of Engineering, Swansea University, Bay Campus, Fabian Way, Swansea, SA1 8EN, UK International
More informationUncertainty Quantification of Two-Phase Flow in Heterogeneous Porous Media
Uncertainty Quantification of Two-Phase Flow in Heterogeneous Porous Media M.Köppel, C.Rohde Institute for Applied Analysis and Numerical Simulation Inria, Nov 15th, 2016 Porous Media Examples: sponge,
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 informationGeneralized Spectral Decomposition for Stochastic Non Linear Problems
Generalized Spectral Decomposition for Stochastic Non Linear Problems Anthony Nouy O.P. Le Maître Preprint submitted to Journal of Computational Physics Abstract We present an extension of the Generalized
More informationNumerical techniques for linear and nonlinear eigenvalue problems in the theory of elasticity
ANZIAM J. 46 (E) ppc46 C438, 005 C46 Numerical techniques for linear and nonlinear eigenvalue problems in the theory of elasticity Aliki D. Muradova (Received 9 November 004, revised April 005) Abstract
More informationPRECONDITIONING MARKOV CHAIN MONTE CARLO SIMULATIONS USING COARSE-SCALE MODELS
PRECONDITIONING MARKOV CHAIN MONTE CARLO SIMULATIONS USING COARSE-SCALE MODELS Y. EFENDIEV, T. HOU, AND W. LUO Abstract. We study the preconditioning of Markov Chain Monte Carlo (MCMC) methods using coarse-scale
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 information