Numerical Approximation of Stochastic Elliptic Partial Differential Equations
|
|
- Abel Tyler
- 5 years ago
- Views:
Transcription
1 Numerical Approximation of Stochastic Elliptic Partial Differential Equations Hermann G. Matthies, Andreas Keese Institut für Wissenschaftliches Rechnen Technische Universität Braunschweig
2 Overview 2. Stochastic modell problem and stochastic PDE 2. Spatial discretisation of random fields (Karhunen-Loève) 3. Stochastic Galerkin Methods 4. Structure of Discrete Equations 5. Solver for linear SPDEs 6. Hierarchical Parallelisation 7. Numerical Example
3 Model Problem 3 Geometry Aquifier 2D Model simple stationary model of groundwater flow (κ(x) u(x)) = f(x) & b.c., x G R d κ(x) u(x) = g(x), x G, u hydraulic head, κ conductivity, f and g sinks and sources.
4 Stochastic Model 4 Uncertainty of system parameters e.g. κ = κ(x, ω) = κ(x) + κ(x, ω), f = f(x, ω), g = g(x, ω) are stochastic fields ω Ω = probability space of all realisations, with probability measure P. Assumption: < κ κ(x, ω) < κ. Possibilities: Transformation κ(x, ω) = φ(x, γ(x, ω)) := F κ(x) of Gaussian field γ with given 2nd order statistic. e.g. κ(x, ω) has marginal β(/2, /2)-distribution, if (with a(x) > ) κ(x, ω) = a(x) + b(x) arccos (π Φ ( γ(x, ω) )) or log-normal distribution (violates assumption) κ(x, ω) = a(x) + exp(γ(x, ω)) Φ(γ(x, ω)) x
5 Realisation of κ(x, ω) 5
6 Stochastic PDE and Variational Form 6 Insert stochastic parameters into PDE stochastic PDE (κ(x, ω) u(x, ω)) = f(x, ω), x G & b.c. on G Solution u(x, ω) is a stochastic field in tensor product form W S u(x, ω) = µ v µ (x)u (µ) (ω) W is a normal spatial Sobolev space, and S a space of random variables, e.g. S = L 2 (Ω, P ) if u has finite variance, otherwise a stochastic distribution. Well-posed variational formulation: Find u W S, such that v W S : a(κ; v, u) := v(x, ω) (κ(x, ω) u(x, ω)) dx dp (ω) = Ω Ω G [ G v(x, ω)f(x, ω) dx + G ] v(x, ω)g(x, ω) ds(x) dp (ω) =: f, v.
7 Theoretical matters 7 Ghanem & Spanos (99) Idea of stochastic Galerkin methods, first formulation. Holden et al. (996): Products between RVs as Wick product. Interpretation doubtful. Classical PDE in space. Existence via S-transform. Benth & Gjerde (998): Convergence rates for polynomial chaos approximation in variational setting. Use of Lax-Milgram. M. & Bucher (999): Variational formulation of Wick product SPDEs, equivalence with S-transform. Besold (2): Existence via Lax-Milgram for normal product SPDEs. Babuška et al. (22): p-version FE-convergence rates for solution of a restricted modell (finite dim. Prob. space of independent RVs). M. & Keese (23): Stability of Galerkin approximation for Karhunen-Loève and polynomial chaos expansion of κ. Nonlinear SPDEs via monotone operator theory.
8 Statistics 8 Desirable: Uncertainty Quantification: The goal is to compute statistics which are deterministic functionals of the solution: Ψ u (x) = E (Ψ(u(x, ))) := Ψ(u(x, ω)) dp (ω) Ω e.g.: ū(x) = E (u(x)), or var u (x) = E ( (ũ(x)) 2), or Pr{u(x) u } = E ( ) χ {u(x) u } Principal Approach:. Discretise spatial (and temporal) part (e.g. via finite elements, finite differences). 2. Approximate stochastic fields in finitely many independent random variables (RVs), stochastic discretisation on spatially continuous or discrete representation. 3. Compute solution/statistics: Via high-dimensional integration (e.g. Monte Carlo, Smolyak, FORM). Or approximate solution with response-surface, then integration stochastic Galerkin methods are a form of response surface computation
9 Example Solution 9 Geometry 2.5 flow = flow out 2 Sources 5 Dirichlet b.c Mean of solution 5 Realization of κ Variance of solution Realization of solution
10 Random Fields Consider region in space G, a random field is a RV κ x (ω) at each x G, alternatively for each ω Ω a random function a realisation κ ω (x) on G Covariance considered at different positions C κ (x, x 2 ) := E ( κ(x, ) κ(x 2, )) ( If κ(x) κ, and C κ (x, x 2 ) = c κ (x x 2 ), process is (weakly) homogeneous.) We have to deal with another RV at each point x. (uncountably many) Need to discretise spatial aspect (generalise Fourier representation). One possibility is the Karhunen-Loève Expansion. Can serve also for system (information) reduction. Need to discretise the random variables. One possibility is the Polynomial Chaos Expansion.
11 Karhunen-Loève Expansion KL mode KL mode Eigenvalue-Problem for Karhunen-Loève Expansion KLE. Other names: Proper Orth. Decomp.(POD), Sing. Value Decomp.(SVD), Prncpl. Comp. Anal.(PCA): C κ (x, y)g j (y) dy = κ 2 jg j (x) G gives spectrum {κ 2 j} and orthogonal KLE eigenfunctions g j (x) Representation of κ: κ(x, ω) = κ(x) + κ j g j (x)ξ j (ω) =: κ j g j (x)ξ j (ω) j= with centred, uncorrelated random variables ξ j (ω). Truncate after m largest eigenvalues optimal in variance expansion in m RVs. j=
12 Representation in Gaussian Random Variables 2 Represent {ξ j (ω)} in (uncorrelated ) independent Gaussian RVs θ m (ω) = (θ (ω),..., θ m (ω)) (one for each KLE function) on an m-dimensional space Θ (m) R m with Gaussian measure dγ m (θ m ) dγ m (θ m ) = (2π) m/2 exp( θ m 2 /2) dθ m m = (2π) /2 exp( θj 2 /2) dθ j j= Each function (RV) ξ j (θ m ) of these Gaussian RVs θ m (ω) (with finite variance) can be approximated arbitrarily well (in variance, i.e. the L 2 (Ω, P )-norm) with polynomials in (θ (ω),..., θ m (ω)) Wiener s Plynomial Chaos Realisation m = Realisation m = 2
13 Polynomial Chaos Expansion 3 Much more is actually true (polynomial chaos expansion PCE): Theorem[Norbert Wiener]: Any RV r(ω) L 2 (Ω, P ) (with finite variance) can be represented in orthogonal polynomials of Gaussian RVs {θ m (ω)} m= =: θ(ω): r(ω) = ϱ (α) H α (θ(ω)), with H α (θ(ω)) = h αj (θ j (ω)), α J j= where h l (ϑ) are the usual Hermite polynomials, and J := {α α = (α,..., α j,...), α j N, α := are multi-indices, where only finitely many of the α j are non-zero. In other words, the H α are an orthogonal basis in L 2 (Ω), with H α, H β = E (H α H β ) = α! δ αβ, where α! := j= (α j!). α j < } j=
14 Computational Approaches 4 Variational formulation discretised in space, e.g. via finite elements u(x, ω) = N(x) T u(ω): A(ω)[u(ω)] = f(ω), or A(θ)[u(θ)] = f(θ). The principal computational approaches are: Monte Carlo / Direct Integration Directly compute statistic by quadrature: Ψ u (x) = E (Ψ(, u(x, ))) = Θ Ψ(θ, N(x)T u(θ)) dγ (θ) by numerical integration. Perturbation Assume that stochastics is a small perturbation around mean value, do Taylor expansion and truncate usually after linear or quadratic term. Response Surface Try to find a functional fit u(θ) v(θ), then compute with v. Stochastic Galerkin Use an ansatz for the solution u(θ) = β u(β) H β (θ) in the stochastic dimension, then do Galerkin method to obtain u (β).
15 Monte Carlo / Direct Integration 5 Numerical integration of Ψ u proceeds as follows:. Select integration rule, i.e. sampling points {θ z } Z z= and corresponding weights w z. 2. At each θ z, first solve the equations for u(θ z ) for fixed realisation θ z, then evaluate the integrand Ψ(θ z, u(x, θ z )). Use of normal, existing software. 3. Approximate Ψ u (x) = Θ Ψ(θ, N(x)T u(θ)) dγ (θ) Z z= w zψ(θ z, N(x) T u(θ z )). For Monte Carlo (MC), the points θ z are random according to measure Γ, and weights are w z = /Z. For Quasi Monte Carlo (QMC), the points θ z are non-random according to number-theoretic low discrepancy series, and weights are still w z = /Z. For Product Gauss or sparse grid Smolyak rules, the points θ z and weights w z are known.
16 High-Dimensional Integration 6 We want E (ψ(θ)) = Θ m ψ(θ m ) dγ m (θ m ) = ψ(θ,..., θ m )dγ (θ ) dγ (θ m ) θ θ m Monte-Carlo (MC) with Z realisations and expected error ɛ: Pure Monte Carlo (MC) ɛ = O( ψ 2 Z /2 ) Quasi Monte Carlo (QMC) ɛ = O ( ψ BV Z (log Z) m) Quadrature-formulas (Integrand in C r (Θ m ): ) Full Product k-point Gauss ɛ = O(Z (2r )/m ) Z = O (k m Sparse Grid Smolyak ɛ = O(Z r (log Z) (m )(r+) ) Z = O ( 2 k k! mk)
17 Perturbation 7 Assume that stochastic input (e.g. κ(x, θ) = κ(x) + κ(x, θ)) is small variation around mean. Do Taylor expansion of solution w.r.t stochastic input quantity. Equate terms of equal order. Approximate mean and covariance of solution. Limited to small stochastic variations, only gives mean and covariance, nothing else; other statistics very difficult and unreliable. Derivatives often computed via dual / adjoint approach.
18 Response Surface 8 Assume a simple functional form for v(θ) u(θ). Usually low order polynomial. Determine by sampling i.e. computing the solution at some point θ m and (least squares) fitting. Use of normal, existing software. Once response surcase v(θ) has been determined, evaluate statistics as in Direct Integration approach with v in case of u, but now evaluation of solution much quicker through response surface. As described works only well for low dimension m of stochastic space. In some sense a Galerkin solution is also a response surface.
19 Stochastic Galerkin 9 Recipe: Stochastic ansatz and projection in stochastic dimensions u(θ) = β H β (θ)u (β) =: H(θ)u, H(θ) = (...H β (θ)...), u = (...u (β)...) T. Goal: Compute coefficients u (β). Through stochastic (intrusive) Galerkin-Methods, γ : E ((f(θ) A(θ)[H(θ)u])H γ (θ)) =, #dof space #dof stoch eqns. 2. Non-intrusive / uncoupled, direct projection, u (β) is also a statistic, integral computed through numerical (e.g. MC, QMC, or Smolyak) quadrature: β : u (β) = E (u(θ)h β (θ)), many problems of size #dof space.
20 Results of Galerkin Method 2 err. 4 mean Galerkin scheme m = 6 k = x y Prob{u(x) > 3.25},.8.6 MCsimulations.4.2 of response surface x y u (α) = E (u(ω)h α ) for α = (,,,,, ) x y Error 4 in u (α) Galerkin scheme x y
21 Galerkin-Methods for the Linear Case 2 γ satisfy: [ β R N(x) E (κ(x, θ)h β (θ)h γ (θ)) N(x) T dx ] u (β) = E (f(θ)h γ (θ)) }{{} =:f (γ) More efficient representation through expansion of κ in KLE and PCE and analytic computation of expectations. κ(x, θ) = κ j ξ j (θ)g j (x) j= r κ j ξ j (α) H α (θ)g j (x), α J 2k,m j= where J k,m = {α J α k, ı > m α ı = }.
22 Resulting Equations 22 Insertion of expansion of κ ( γ satisfy): ξ j (α) E (H α H β H γ ) N(x)κ }{{} j g j (x) N(x) T dx α β j }{{} =: (α) β,γ K j u (β) = f (γ) K j is stiffness matrix of a FEM discretisation for the material κ j g j (x). Use deterministic FEM program in black-box-fashion. Need:. Possibility to change material-parameter 2. Function returning mass-matrix R N(x)N(x)T dx for numerical KLE-decomposition. 3. Functions, which give residuum and (possibly) Jacobi-matrix K j. 4. Function, which will solve for a realisation, u(θ) = K (θ)f(θ)
23 Approximation Theory 23 Stability of discrete approximation under truncated KLE and PCE (M. & Keese, 23). Matrix stays uniformly positive definite. Convergence follows from Céas Lemma. Convergence rates under stochastic regularity (Benth & Gjerde, 998) in stochastic Hilbert spaces. Error estimation via dual weighted residuals possible. Stochastic Hilbert spaces start with formal PCE: R(θ) = α J R(α) H α (θ). Define for ρ and p norm (with (2N) β := j N (2j)β j ): R 2 ρ,p = α R (α) 2 (α!) +ρ (2N) pα.
24 Stochastic Hilbert Spaces Convergence Rates 24 Define for ρ, p : (S) ρ,p = {R(θ) = α J R (α) H α (θ) : R ρ,p < }. These are Hilbert spaces, the duals are denoted by (S) ρ, p, and L 2 (Ω) = (S),. Let P k,m be orthogonal projection from (S) ρ,p into the finite dimensional subspace spanned by {H α }, α J k,m. Theorem[Benth & Gjerde, Keese]: Let p R, q > p + r, where r.53 solves r = 2 r (r ), and let ρ. Then for any R (S) ρ,p : R P k,m (R) ρ,p R ρ,q c(m, k, q p), ( k+, where c(m, k, r) = c (r)m r r + c 2 (r) 2 (r )) and r c (r) = (2 r (r ) r), c 2 (r) = 2 r (r )c (r). E.g., if r = 2, then c(m, k, r) 2 = /2(m + (/2) k ).
25 Tensor Product Structure 25 The equations have the structure of a tensor product K u = ξ j (α) (α) K j u = f j α ξ j (α) (α) β,γ K j (α) β 2,γ K j (α) β,γ 2 K j (α) β 2,γ 2 K j u(β ). = α j..... u (β N) f (γ ).. f (γ N) #dof space #dof stoch linear equations. Exploit parallelism in the multiplication: parallel operator-sum, distributed block-vector Block-matrix efficiently stored and used in tensor-representation.
26 Sparsity Structure 26 Non-zero blocks for increasing degree of H α
27 Matrix Expansion Errors 27 H α is orthogonal on polynomials of degree < α (α) βγ = E (H αh β H γ ) = for α > β + γ, hence sum over α is finite. One can choose J such that expansion J ξ j (α) (α) K j u j= α is arbitrarily close to = N(x) E (κ(x, θ)h β (θ)h γ (θ)) N(x) T dx u (β), β R and hence global matrix uniformly positive definite. This stability argument is only valid for the (intrusive) Galerkin method.
28 Properties of Global Equations 28 K u = j α ξ j (α) (α) K j u = f Each K j is symmetric, and each (α). Block-matrix K is symmetric. SPDE ist positive definite. Appropriate expansion of κ K is uniformly positive definite stability. Solving the Equations: Never assemble block-matrix explicitly. Use K only as multiplication. Use Krylov Method (here CG) with pre-conditioner.
29 Block-Diagonaler Pre-Conditioner 29 Let K = K = stiffness-matrix for average material κ(x). Use deterministic solver as pre-conditioner: P = K = I K... K Good pre-conditioner, when variance of κ not too large. Otherwise use P = block-diag(k). This may again be done with existing deterministic solver. Block-diagonal P is well suited for parallelisation.
30 Parallelising the Matrix-Vector Product 3 γ : (K u) (γ) = J j N β n α α ξ (α) j (α) β,γ K j u β K j = deterministic solver. This may be a (lower-level) parallel program to do K j u β. Parallelise operator-sum in j several instances of deterministic solver in parallel. Distribute u and f Parallelise sum in β. Sum in α may also be done in parallel, but usually not essential.
31 Parallelisation 3 γ : (K u) (γ) = j β α ξ (α) j (α) β,γ K j u β Obviously Parallel in γ. Block-vectors u and f distributed. May be replicated, in order to reduce communication. Matrices K j distributed over processors. May be replicated, in order to reduce parallel communication. use more processors than number of K j. Several processor-groups, where each uses a subset of the K j and stores a subset of u and f.
32 Speedup Measurements I 32 On Cray T3E Constant problem size, N space = 75, N stoch = 54, (9 KLE-terms, degree 3 polynomial chaos), total number of equations, 55, 5. 2 operators K j Scaling up to > 5 6 equations. Distributed operator, distributed block-vector, (best use of memory): Efficiency Relative time used Number of processors Number of processors
33 Speedup Measurements II 33 Replicated operator, distributed block-vector (compromise): Efficiency Number of processors Operators distributed, replicated block-vector (parallelise for speed): Efficiency Number of processors
34 Speedup Measurements III 34 Execution time, problem size scaled with number of processors Time [s] Number of processors
35 Outlook 35 Stochastic problems at very beginning (like FEM in the 6 s). Deterministic programs can be used in this framework. Nonlinear (and instationary) problems possible (but much more work). Development of framework for stochastic coupling and parallelisation. Hierarchical parallelisation well possible, with different objectives: Memory Saver. Time Saver. Combination of both.
36 At Last Thank you for your attention! URL:
Hierarchical 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 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 informationPartial Differential Equations with Stochastic Coefficients
Partial Differential Equations with Stochastic Coefficients Hermann G. Matthies gemeinsam mit Andreas Keese Institut für Wissenschaftliches Rechnen Technische Universität Braunschweig wire@tu-bs.de http://www.wire.tu-bs.de
More informationStochastic Finite Elements: Computational Approaches to SPDEs
Stochastic Finite Elements: Computational Approaches to SPDEs Hermann G. Matthies Andreas Keese, Elmar Zander, Alexander Litvinenko, Bojana Rosić Technische Universität Braunschweig wire@tu-bs.de http://www.wire.tu-bs.de
More informationGalerkin Methods for Linear and Nonlinear Elliptic Stochastic Partial Differential Equations
ScientifiComputing Galerkin Methods for Linear and Nonlinear Elliptic Stochastic Partial Differential Equations Hermann G. Matthies, Andreas Keese Institute of Scientific Computing Technical University
More informationSampling and Low-Rank Tensor Approximations
Sampling and Low-Rank Tensor Approximations Hermann G. Matthies Alexander Litvinenko, Tarek A. El-Moshely +, Brunswick, Germany + MIT, Cambridge, MA, USA wire@tu-bs.de http://www.wire.tu-bs.de $Id: 2_Sydney-MCQMC.tex,v.3
More informationParametric Problems, Stochastics, and Identification
Parametric Problems, Stochastics, and Identification Hermann G. Matthies a B. Rosić ab, O. Pajonk ac, A. Litvinenko a a, b University of Kragujevac c SPT Group, Hamburg wire@tu-bs.de http://www.wire.tu-bs.de
More informationNon-Intrusive Solution of Stochastic and Parametric Equations
Non-Intrusive Solution of Stochastic and Parametric Equations Hermann G. Matthies a Loïc Giraldi b, Alexander Litvinenko c, Dishi Liu d, and Anthony Nouy b a,, Brunswick, Germany b École Centrale de Nantes,
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 informationOptimisation under Uncertainty with Stochastic PDEs for the History Matching Problem in Reservoir Engineering
Optimisation under Uncertainty with Stochastic PDEs for the History Matching Problem in Reservoir Engineering Hermann G. Matthies Technische Universität Braunschweig wire@tu-bs.de http://www.wire.tu-bs.de
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 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 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 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 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 informationNON-LINEAR APPROXIMATION OF BAYESIAN UPDATE
tifica NON-LINEAR APPROXIMATION OF BAYESIAN UPDATE Alexander Litvinenko 1, Hermann G. Matthies 2, Elmar Zander 2 http://sri-uq.kaust.edu.sa/ 1 Extreme Computing Research Center, KAUST, 2 Institute of Scientific
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 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 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 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 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 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 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 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 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 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 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 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 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 informationPlatzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen PARAMetric UNCertainties, Budapest STOCHASTIC PROCESSES AND FIELDS Noémi Friedman Institut für Wissenschaftliches Rechnen, wire@tu-bs.de
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 informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
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 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 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 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 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 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 informationProper Generalized Decomposition for Linear and Non-Linear Stochastic Models
Proper Generalized Decomposition for Linear and Non-Linear Stochastic Models Olivier Le Maître 1 Lorenzo Tamellini 2 and Anthony Nouy 3 1 LIMSI-CNRS, Orsay, France 2 MOX, Politecnico Milano, Italy 3 GeM,
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 informationMonte Carlo Methods for Uncertainty Quantification
Monte Carlo Methods for Uncertainty Quantification Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Mike Giles (Oxford) Monte Carlo methods 1 / 23 Lecture outline
More informationNumerical Solution I
Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions
More informationA Reduced Basis Approach for Variational Problems with Stochastic Parameters: Application to Heat Conduction with Variable Robin Coefficient
A Reduced Basis Approach for Variational Problems with Stochastic Parameters: Application to Heat Conduction with Variable Robin Coefficient Sébastien Boyaval a,, Claude Le Bris a, Yvon Maday b, Ngoc Cuong
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 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 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 informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
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 informationApproximating Infinity-Dimensional Stochastic Darcy s Equations without Uniform Ellipticity
Approximating Infinity-Dimensional Stochastic Darcy s Equations without Uniform Ellipticity Marcus Sarkis Jointly work with Juan Galvis USC09 Marcus Sarkis (WPI) SPDE without Uniform Ellipticity USC09
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 informationACM/CMS 107 Linear Analysis & Applications Fall 2017 Assignment 2: PDEs and Finite Element Methods Due: 7th November 2017
ACM/CMS 17 Linear Analysis & Applications Fall 217 Assignment 2: PDEs and Finite Element Methods Due: 7th November 217 For this assignment the following MATLAB code will be required: Introduction http://wwwmdunloporg/cms17/assignment2zip
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
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 informationParadigms of Probabilistic Modelling
Paradigms of Probabilistic Modelling Hermann G. Matthies Brunswick, Germany wire@tu-bs.de http://www.wire.tu-bs.de abstract RV-measure.tex,v 4.5 2017/07/06 01:56:46 hgm Exp Overview 2 1. Motivation challenges
More informationProbabilistic Collocation Method for Uncertainty Analysis of Soil Infiltration in Flood Modelling
Probabilistic Collocation Method for Uncertainty Analysis of Soil Infiltration in Flood Modelling Y. Huang 1,2, and X.S. Qin 1,2* 1 School of Civil & Environmental Engineering, Nanyang Technological University,
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 informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationUncertainty Quantification and related areas - An overview
Uncertainty Quantification and related areas - An overview Alexander Litvinenko, talk given at Department of Mathematical Sciences, Durham University Bayesian Computational Statistics & Modeling, KAUST
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 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 informationEstimating functional uncertainty using polynomial chaos and adjoint equations
0. Estimating functional uncertainty using polynomial chaos and adjoint equations February 24, 2011 1 Florida State University, Tallahassee, Florida, Usa 2 Moscow Institute of Physics and Technology, Moscow,
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 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 Solvers for the Euler Equations
43rd AIAA Aerospace Sciences Meeting and Exhibit 1-13 January 5, Reno, Nevada 5-873 Stochastic Solvers for the Euler Equations G. Lin, C.-H. Su and G.E. Karniadakis Division of Applied Mathematics Brown
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 informationSecond Order Elliptic PDE
Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic
More informationStochastic Galerkin Method for Elliptic Spdes: A White Noise Approach
Worcester Polytechnic Institute igital WPI Mathematical Sciences Faculty Publications epartment of Mathematical Sciences 7-1-26 Stochastic Galerkin Method for Elliptic Spdes: A White Noise Approach Luis
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 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 Framework for Statistical Timing Analysis using Non-Linear Delay and Slew Models
A Framework for Statistical Timing Analysis using Non-Linear Delay and Slew Models Sarvesh Bhardwaj, Praveen Ghanta, Sarma Vrudhula Department of Electrical Engineering, Department of Computer Science
More informationChapter 2 Spectral Expansions
Chapter 2 Spectral Expansions In this chapter, we discuss fundamental and practical aspects of spectral expansions of random model data and of model solutions. We focus on a specific class of random process
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 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 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 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 information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
More informationStochastic Spectral Methods for Uncertainty Quantification
Stochastic Spectral Methods for Uncertainty Quantification Olivier Le Maître 1,2,3, Omar Knio 1,2 1- Duke University, Durham, North Carolina 2- KAUST, Saudi-Arabia 3- LIMSI CNRS, Orsay, France 38th Conf.
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 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 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 informationAdaptive Stroud Stochastic Collocation Method for Flow in Random Porous Media via Karhunen-Loève Expansion
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 4, No. 1, pp. 12-123 Commun. Comput. Phys. July 28 Adaptive Stroud Stochastic Collocation Method for Flow in Random Porous Media via Karhunen-Loève Epansion
More informationPolynomial Chaos Expansion of random coefficients and the solution of stochastic partial differential equations in the Tensor Train format
Polynomial Chaos Expansion of random coefficients and the solution of stochastic partial differential equations in the Tensor Train format arxiv:1503.03210v1 [math.na] 11 Mar 2015 Sergey Dolgov Boris N.
More informationProbabilistic Structural Dynamics: Parametric vs. Nonparametric Approach
Probabilistic Structural Dynamics: Parametric vs. Nonparametric Approach S Adhikari School of Engineering, Swansea University, Swansea, UK Email: S.Adhikari@swansea.ac.uk URL: http://engweb.swan.ac.uk/
More informationLecture 1: Introduction to low-rank tensor representation/approximation. Center for Uncertainty Quantification. Alexander Litvinenko
tifica Lecture 1: Introduction to low-rank tensor representation/approximation Alexander Litvinenko http://sri-uq.kaust.edu.sa/ KAUST Figure : KAUST campus, 5 years old, approx. 7000 people (include 1400
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 information256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.
56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =
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 informationFast evaluation of mixed derivatives and calculation of optimal weights for integration. Hernan Leovey
Fast evaluation of mixed derivatives and calculation of optimal weights for integration Humboldt Universität zu Berlin 02.14.2012 MCQMC2012 Tenth International Conference on Monte Carlo and Quasi Monte
More informationA very short introduction to the Finite Element Method
A very short introduction to the Finite Element Method Till Mathis Wagner Technical University of Munich JASS 2004, St Petersburg May 4, 2004 1 Introduction This is a short introduction to the finite element
More informationFrom Completing the Squares and Orthogonal Projection to Finite Element Methods
From Completing the Squares and Orthogonal Projection to Finite Element Methods Mo MU Background In scientific computing, it is important to start with an appropriate model in order to design effective
More informationFFT-based Galerkin method for homogenization of periodic media
FFT-based Galerkin method for homogenization of periodic media Jaroslav Vondřejc 1,2 Jan Zeman 2,3 Ivo Marek 2 Nachiketa Mishra 4 1 University of West Bohemia in Pilsen, Faculty of Applied Sciences Czech
More informationApplication of QMC methods to PDEs with random coefficients
Application of QMC methods to PDEs with random coefficients a survey of analysis and implementation Frances Kuo f.kuo@unsw.edu.au University of New South Wales, Sydney, Australia joint work with Ivan Graham
More informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
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 informationMatrix assembly by low rank tensor approximation
Matrix assembly by low rank tensor approximation Felix Scholz 13.02.2017 References Angelos Mantzaflaris, Bert Juettler, Boris Khoromskij, and Ulrich Langer. Matrix generation in isogeometric analysis
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationStrain and stress computations in stochastic finite element. methods
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2007; 00:1 6 [Version: 2002/09/18 v2.02] Strain and stress computations in stochastic finite element methods Debraj
More informationA Stochastic Finite Element Method with a Deviatoric-volumetric Split for the Stochastic Linear Isotropic Elasticity Tensor
TECHNISCHE MECHANIK,,, (), 79 submitted: April, A Stochastic Finite Element Method with a Deviatoric-volumetric Split for the Stochastic Linear Isotropic Elasticity Tensor A. Dridger, I. Caylak,. Mahnken
More informationLehrstuhl Informatik V. Lehrstuhl Informatik V. 1. solve weak form of PDE to reduce regularity properties. Lehrstuhl Informatik V
Part I: Introduction to Finite Element Methods Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Necel Winter 4/5 The Model Problem FEM Main Ingredients Wea Forms and Wea
More information