Sampling and low-rank tensor approximation of the response surface
|
|
- Marjory Powell
- 5 years ago
- Views:
Transcription
1 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, Extreme Computing Research Center KAUST, 3 TU Braunschweig, Germany
2 KAUST Figure : KAUST campus, 5 years old, approx people (include 1400 kids), 100 nations. 2 / 43
3 Stochastic Numerics Group at KAUST 3 / 43
4 4 / 43
5 Points of my study and of my work. 5 / 43
6 The structure of the talk 1. Results of MUNA (input: α, Ma, output CL, CD, cp, cf) 2. Low-rank response surface 3. Numerics (input: α, Ma, output CL, CD, density, pressure,...) 4. Use the response surface to get a good start value 5. Random geometry of the airfoil 6. Comparison of 5 methods 6 / 43
7 Results of MUNA project ( ) are published 390 pages, our results on pp / 43
8 Sparse Gauss-Hermite grid, 13 points. Assume that RVs α and Ma are Gaussian with mean st. dev. σ/mean σ α Ma Then uncertainties in the solution lift CL and drag CD are CL CD / 43
9 Example 1: 500 MC realisations of cp in dependence on α i and Ma i 9 / 43
10 Example 2: 500 MC realisations of cf in dependence on α i and Ma i 10 / 43
11 Example 3: 5% and 95% quantiles for cp from 500 MC realisations. 11 / 43
12 Example 4: 5% and 95% quantiles for cf from 500 MC realisations. 12 / 43
13 Compute difference between the mean solution and the deterministic solution (α = 2.74, Ma = 0.73). The difference is large! It is comparable with the lowest physical value! (a) deterministic turbulent viscosity (2e 5, 3.6e 4) (b) deterministic density (0.45, 1.25) 13 / 43
14 Response surface Let RF q(x, θ), θ = (θ 1,..., θ M,...) is approximated: q(x, θ) = β J q β (x)h β (θ) (1) q β (x) = 1 H β (θ)q(x, θ) P(dθ) 1 n q H β (θ i )q(x, θ i )w i, β! Θ β! where n q - number of quadrature points. Using low-rank format, obtain tion Logo q β (x) Lock-up = 1 β! [q(x, θ 1),..., q(x, θ nq )] [H β (θ 1 )w 1,..., H β (θ nq )w nq ] T i=1 14 / 43
15 PCE coefficients Denote c β := 1 β! [H β(θ 1 )w 1,..., H β (θ nq )w nq ] T R nq and approximate the set of realisations in low-rank format: [q(x, θ 1 ),..., q(x, θ nq )] AB T. The matrix of all PCE coefficients will be R N J [...q β (x)...] AB T [...c β...], β J. (2) Put all together, obtain low-rank representation of RS q(x, θ) = q β (x)h β (θ) = [...q β (x)...]h T (θ), (3) β J where H(θ) = (..., H β (θ),...). 15 / 43
16 Application of low-rank response surface Now, having discretised RS q(x, θ) q(θ) = AB T [...c β...]h T (θ) (4) Sample RV θ 10 6 times and then use the obtained sample to compute errorbars, quantiles, cumulative density function. 16 / 43
17 Use in MC sampling solution storage Inflow and air-foil shape uncertain. Data compression achieved by updated SVD: Made from 600 MC Simulations, SVD is updated every 10 samples. n = 260, 000 Z = 600 Updated SVD: Relative errors, memory requirements: rank k pressure turb. kin. energy memory [MB] e-2 4.0e e-2 5.9e e-3 1.5e Dense matrix R costs 1250 MB storage. 17 / 43
18 Conjunction of two low-rank matrices W k = [q(x, θ 1 ),..., q(x, θ Z )], W = [q(x, θ Z+1 ),..., q(x, θ Z+m )] Suppose W k = AB T R n Z is given Suppose W R n m contains new m solution vectors Compute C R n k and D R m k such that W CD T. Build Anew := [A C] R n 2k and B T new = blockdiag[b T D T ] R 2k (Z+m) Rank-k truncation of Wnew = AnewB T new costs O((n + Z + m)k 2 + k 3 ) 18 / 43
19 Mean and variance in the rank-k format u := 1 Z Cost is O(k(Z + n)). Z u i = 1 Z i=1 Z A b i = Ab. (5) i=1 C = 1 Z 1 W cw T c 1 Z 1 U kσ k Σ T k UT k. (6) Cost is O(k 2 (Z + n)). Lemma: Let W W k 2 ε, and u k be a rank-k approximation of the mean u. Then a) u u k ε Z, b) C C k 1 Z 1 ε2. 19 / 43
20 Tensor product structure Story does not end here as one may choose S = k S k, approximated by S B = K k=1 S B k, with S Bk S k. Solution represented( as a tensor of grade K + 1 K ) in W B,N = k=1 S B k U N. For higher grade tensor product structure, more reduction is possible, but that is a story for another talk, here we stay with K = / 43
21 Numerics Numerical experiments Compare MC with PCE-137, PCE-281 RAE-2822, Case 1, α = 1.93 and Ma = 0.676: (mean and variance) No difference, 137 is enough. Compare MC with PCE-29, PCE-201 tion Logo RAE-2822, Lock-up Case 9, α = 2.79 and Ma = 0.73: (mean and variance) 21 / 43
22 Case 1, density (a) ρ PCE137 ρ MC F = relative error: (b) ρ PCE281 ρ MC F = relative error: (c) var(ρ PCE137 ) var(ρ MC ) F = 0.25 (d) var(ρ PCE281 ) var(ρ MC ) F = / 43
23 Case 1, Pressure (a) p PCE137 p MC F = relative error: (b) p PCE281 p MC F = relative error: (c) var(p PCE137 ) var(p MC ) F = 0.47 (d) var(p PCE281 ) var(p MC ) F = / 43
24 Case 9, density (a) ρ PCE29 ρ MC F = 7.0 relative error: (b) ρ PCE201 ρ MC F = 7.2 relative error: (c) var(ρ PCE29 ) var(ρ MC ) F = (d) var(ρ PCE201 ) var(ρ MC ) F = / 43
25 Case 9, Pressure (a) p PCE29 p MC F = relative error: (b) p PCE201 p MC F = relative error: (c) var(p PCE29 ) var(p MC ) F = 1.3 (d) var(p PCE201 ) var(p MC ) F = / 43
26 Relative error, density mean, trans-sonic flow, Z = 2600 samples 26 / 43
27 Relative error, density variance, trans-sonic flow, Z = / 43
28 Figure : Comparison of the mean pressures computed with PCE and with MC. (Left) p := p PCE137 p MC, Case 1 without shock, (Right) p := p PCE201 p MC, Case 9 with shock. 28 / 43
29 Sampling of response surface with residual Collocation points are θ l, l = 1..Z. Algorithm: 1. Compute RS q(x, θ) AB T [...c β...]h T (θ) from l points. 2. (l = l + 1), evaluate RS in θ l+1, obtain q(x, θ l+1 ). 3. Compute residual r(q(x, θ)). Only if r is large, solve expensive determ. problem. 4. Update A, B T, [...c β...] and go to (2). If we are lucky, we solve the determ. problem only few times, otherwise we must solve the determ. problem Z times for all tion Logo θ 1,...,θ Lock-up Z. 29 / 43
30 Errors, Case 1 Figure : (left) Relative errors in the Frobenius and the maximum norms for pressure and density. (right) 10 points (α, Ma) were chosen in the neigbourhood of α = 1.93 and Ma = / 43
31 Density, mean and variance The mean density and variance of the density. Case 9, RAE / 43
32 Density with preconditioning failed Density evaluated from two different PCE-based response surfaces (of order p = 2 and p = 4). Failed. 32 / 43
33 Uncertain geometry of the airfoil Uncertain geometry of the airfoil. Comparison of 5 methods. 33 / 43
34 Uncertain geometry of the airfoil [1]. of airfoil geometry-induced aerodynamic uncertainties - comparison of approaches, D. Liu, A. Litvinenko,C. Schillings, V. Schulz, AIAA, 2015 Input: ψ (x, ω) - Gaussian RF, R(x, ω) = s(x) arccos (1 2Φ (ψ (x, ω))), G(ω) = G + R(x, ω) n(x), perturbed airfoil 34 / 43
35 Uncertainty quantification - comparison of methods To compare the efficiency of methods in estimating the statistics. Methods: gradient-employing gradient-assisted point-collocation polynomial chaos (GAPC) gradient-assisted radial basis function (GARBF) gradient-enhanced Kriging (GEK) not gradient-employing quasi-monte Carlo quadrature (QMC), with low discrepancy sequence polynomial chaos (PC), with coefficients estimated by sparse grid quadrature Criteria: computation cost for a certain accuracy in the statistics. Cost: measured in elapse time-penalized sample number M M = 2N for gradient-employing method M = N for others Accuracy: judged by comparing with a reference statistics obtained by a QMC of N = / 43
36 Three realizations of perturbed geometry Original 1st realization 2nd realization 3rd realization 0.02 y x Target statistics: means µ L and µ D, standard deviations σ L and σ D, exceedance probabilities P L,κ = Pro{C L µ L κ σ L } and P D,κ = Pro{C D µ D + κ σ D } with κ = 2, / 43
37 On estimating the mean of C L On estimating the stdv of C L 1e 03 1e 04 QMC PC SGH GEK GEPC GERBF 1e 03 1e 04 QMC PC SGH GEK GEPC GERBF Error Error 1e 05 1e 05 1e M M 1e 01 On estimating P L,2 On estimating P L,3 1e 02 1e 02 Error QMC 1e 03 PC SGH GEK GEPC GERBF tion Logo 1e 04 Lock-up 3ς 1 Error 1e 03 1e 04 QMC PC SGH GEK GEPC GERBF 3ς M M Figure : Absolute errors 37 / 43
38 On estimating the mean of C D On estimating the stdv of C D 1e 04 QMC PC SGH GEK GEPC GERBF 1e 04 QMC PC SGH GEK GEPC GERBF Error 1e 05 Error 1e 05 1e M 1e M On estimating P D,2 On estimating P D,3 1e 02 1e 01 1e 02 QMC PC SGH GEK GEPC GERBF 3ς 1 Error 1e 03 QMC PC SGH GEK GEPC 1e 04 GERBF 3ς 1 Error 1e 03 1e M M Figure : Absolute errors 38 / 43
39 Figure : Estimated pdf (in dash line) of C L by QMC, GEK, GEPC and GERBF at M = / 43 pdf of C L by QMC, M=33 pdf of C L by GEK, M= reference pdf pdf by QMC reference pdf pdf by GEK pdf pdf C L C L pdf of C L by GEPC, M=33 pdf of C L by GERBF, M= reference pdf pdf by GEPC reference pdf pdf by GERBF pdf pdf C L C L
40 Figure : Estimated pdf (in dash line) of C D by QMC, GEK, GEPC and GERBF at M = / pdf of C D by QMC, M=33 reference pdf pdf by QMC pdf of C D by GEK, M=33 reference pdf pdf by GEK pdf 400 pdf C L x 10 3 pdf of C D by GEPC, M=33 reference pdf pdf by GEPC C L x pdf of C D by GERBF, M=33 reference pdf pdf by GERBF pdf pdf C L x C L x 10 3
41 Uncertainty quantification - conclusion Conclusion: Gradient-employing surrogate methods outperform the others GEK seems more efficient than other gradient-employing methods, especially when M is small, and when estimating the far-end exceedance probability and the pdf. Since PC only has two data point, its performance in this comparison may not indicate its asymptotic capacity. 41 / 43
42 Acknowledgement 1. Nathalie Rauschmayr for pictures 2. Project MUNA, German Luftfahrtforschungsprogramm funded by the Ministry of Economics (BMWA). 3. Elmar Zander for sglib. Type in your terminal git clone git://github.com/ezander/sglib.git initialize all variables, run startup.m You will find: generalised PCE, sparse grids, (Q)MC, stochastic Galerkin, linear solvers, KLE, covariance matrices, statistics, quadratures (multivariate Chebyshev, Laguerre, Lagrange, Hermite ) etc There are: many examples, many test, rich demos To 42 / 43
43 Literature 1. PCE of random coefficients and the solution of stochastic partial differential equations in the Tensor Train format, S. Dolgov, B. N. Khoromskij, A. Litvinenko, H. G. Matthies, 2015/3/11, arxiv: Efficient analysis of high dimensional data in tensor formats, M. Espig, W. Hackbusch, A. Litvinenko, H.G. Matthies, E. Zander Sparse Grids and Applications, 31-56, 40, Application of hierarchical matrices for computing the Karhunen-Loeve expansion, B.N. Khoromskij, A. Litvinenko, H.G. Matthies, Computing 84 (1-2), 49-67, 31, Efficient low-rank approximation of the stochastic Galerkin matrix in tensor formats, M. Espig, W. Hackbusch, A. Litvinenko, H.G. Matthies, P. Waehnert, Comp. & Math. with Appl. 67 (4), , Numerical Methods for Uncertainty and Bayesian Update in Aerodynamics, A. Litvinenko, H. G. Matthies, Book Management and Minimisation of Uncertainties and Errors in Numerical Aerodynamics pp , / 43
NON-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 informationQuantification of airfoil geometry-induced. aerodynamic uncertainties - comparison of. approaches
Quantification of airfoil geometry-induced aerodynamic uncertainties - comparison of approaches Dishi Liu 1 Lilienthalplatz 7, 3818 Braunschweig, Germany Alexander Litvinenko 2 47 KAUST, M.B. 4436, Thuwal,
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 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 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 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 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 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 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 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 informationHierarchical matrix techniques for maximum likelihood covariance estimation
Hierarchical matrix techniques for maximum likelihood covariance estimation tifica Alexander Litvinenko, Extreme Computing Research Center and Uncertainty Center, KAUST (joint work with M. Genton, Y. Sun
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 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 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 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 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 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 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 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 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 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 informationPolynomial chaos expansions for structural reliability analysis
DEPARTMENT OF CIVIL, ENVIRONMENTAL AND GEOMATIC ENGINEERING CHAIR OF RISK, SAFETY & UNCERTAINTY QUANTIFICATION Polynomial chaos expansions for structural reliability analysis B. Sudret & S. Marelli Incl.
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 informationOverview. Bayesian assimilation of experimental data into simulation (for Goland wing flutter) Why not uncertainty quantification?
Delft University of Technology Overview Bayesian assimilation of experimental data into simulation (for Goland wing flutter), Simao Marques 1. Why not uncertainty quantification? 2. Why uncertainty quantification?
More 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 informationUncertainty Quantification in Computational Models
Uncertainty Quantification in Computational Models Habib N. Najm Sandia National Laboratories, Livermore, CA, USA Workshop on Understanding Climate Change from Data (UCC11) University of Minnesota, Minneapolis,
More informationUncertainty Management and Quantification in Industrial Analysis and Design
Uncertainty Management and Quantification in Industrial Analysis and Design www.numeca.com Charles Hirsch Professor, em. Vrije Universiteit Brussel President, NUMECA International The Role of Uncertainties
More informationA regularized least-squares method for sparse low-rank approximation of multivariate functions
Workshop Numerical methods for high-dimensional problems April 18, 2014 A regularized least-squares method for sparse low-rank approximation of multivariate functions Mathilde Chevreuil joint work with
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 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 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 informationCurriculum Vitae. Alexander Litvenenko
Curriculum Vitae Alexander Litvinenko https://ecrc.kaust.edu.sa/pages/litvinenko.aspx http://stochastic_numerics.kaust.edu.sa/pages/litvinenkoalexander.aspx Name: Nationality: Languages: Alexander Litvenenko
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 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 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 informationPC EXPANSION FOR GLOBAL SENSITIVITY ANALYSIS OF NON-SMOOTH FUNCTIONALS OF UNCERTAIN STOCHASTIC DIFFERENTIAL EQUATIONS SOLUTIONS
PC EXPANSION FOR GLOBAL SENSITIVITY ANALYSIS OF NON-SMOOTH FUNCTIONALS OF UNCERTAIN STOCHASTIC DIFFERENTIAL EQUATIONS SOLUTIONS M. Navarro, O.P. Le Maître,2, O.M. Knio,3 mariaisabel.navarrojimenez@kaust.edu.sa
More informationQuadrature for Uncertainty Analysis Stochastic Collocation. What does quadrature have to do with uncertainty?
Quadrature for Uncertainty Analysis Stochastic Collocation What does quadrature have to do with uncertainty? Quadrature for Uncertainty Analysis Stochastic Collocation What does quadrature have to do with
More 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 informationSteps in Uncertainty Quantification
Steps in Uncertainty Quantification Challenge: How do we do uncertainty quantification for computationally expensive models? Example: We have a computational budget of 5 model evaluations. Bayesian inference
More 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 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 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 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 informationAt A Glance. UQ16 Mobile App.
At A Glance UQ16 Mobile App Scan the QR code with any QR reader and download the TripBuilder EventMobile app to your iphone, ipad, itouch or Android mobile device. To access the app or the HTML 5 version,
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 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 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 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 informationExperiences with Model Reduction and Interpolation
Experiences with Model Reduction and Interpolation Paul Constantine Stanford University, Sandia National Laboratories Qiqi Wang (MIT) David Gleich (Purdue) Emory University March 7, 2012 Joe Ruthruff (SNL)
More informationSequential Importance Sampling for Rare Event Estimation with Computer Experiments
Sequential Importance Sampling for Rare Event Estimation with Computer Experiments Brian Williams and Rick Picard LA-UR-12-22467 Statistical Sciences Group, Los Alamos National Laboratory Abstract Importance
More informationRobust Optimal Control using Polynomial Chaos and Adjoints for Systems with Uncertain Inputs
2th AIAA Computational Fluid Dynamics Conference 27-3 June 211, Honolulu, Hawaii AIAA 211-369 Robust Optimal Control using Polynomial Chaos and Adjoints for Systems with Uncertain Inputs Sriram General
More informationA Posteriori Adaptive Low-Rank Approximation of Probabilistic Models
A Posteriori Adaptive Low-Rank Approximation of Probabilistic Models Rainer Niekamp and Martin Krosche. Institute for Scientific Computing TU Braunschweig ILAS: 22.08.2011 A Posteriori Adaptive Low-Rank
More informationsuch that an objective like the drag is minimized subject to important aerodynamic restrictions
OPTIMAL AERODYNAMIC DESIGN UNDER SHAPE UNCERTAINTIES VOLKER SCHULZ AND CLAUDIA SCHILLINGS Abstract. The unavoidable presence of uncertainties poses several difficulties to the numerical treatment of optimization
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 informationInverse Problems in a Bayesian Setting
Inverse Problems in a Bayesian Setting Hermann G. Matthies, Elmar Zander, Bojana V. Rosić, Alexander Litvinenko, Oliver Pajonk arxiv:1511.00524v1 [math.pr] 2 Nov 2015 Institute of Scientific Computing
More informationA Non-Intrusive Polynomial Chaos Method For Uncertainty Propagation in CFD Simulations
An Extended Abstract submitted for the 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada January 26 Preferred Session Topic: Uncertainty quantification and stochastic methods for CFD A Non-Intrusive
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 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 informationQuantifying conformation fluctuation induced uncertainty in bio-molecular systems
Quantifying conformation fluctuation induced uncertainty in bio-molecular systems Guang Lin, Dept. of Mathematics & School of Mechanical Engineering, Purdue University Collaborative work with Huan Lei,
More informationUncertainty Quantification in Discrete Fracture Network Models
Uncertainty Quantification in Discrete Fracture Network Models Claudio Canuto Dipartimento di Scienze Matematiche, Politecnico di Torino claudio.canuto@polito.it Joint work with Stefano Berrone, Sandra
More 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 informationAccuracy, Precision and Efficiency in Sparse Grids
John, Information Technology Department, Virginia Tech.... http://people.sc.fsu.edu/ jburkardt/presentations/ sandia 2009.pdf... Computer Science Research Institute, Sandia National Laboratory, 23 July
More 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 informationEvaluation of Non-Intrusive Approaches for Wiener-Askey Generalized Polynomial Chaos
49th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference 6t 7 - April 28, Schaumburg, IL AIAA 28-892 Evaluation of Non-Intrusive Approaches for Wiener-Askey Generalized
More informationScalable Algorithms for Optimal Control of Systems Governed by PDEs Under Uncertainty
Scalable Algorithms for Optimal Control of Systems Governed by PDEs Under Uncertainty Alen Alexanderian 1, Omar Ghattas 2, Noémi Petra 3, Georg Stadler 4 1 Department of Mathematics North Carolina State
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 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 informationNonlinear stochastic Galerkin and collocation methods: application to a ferromagnetic cylinder rotating at high speed
Nonlinear stochastic Galerkin and collocation methods: application to a ferromagnetic cylinder rotating at high speed Eveline Rosseel Herbert De Gersem Stefan Vandewalle Report TW 541, July 29 Katholieke
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 informationA new Hierarchical Bayes approach to ensemble-variational data assimilation
A new Hierarchical Bayes approach to ensemble-variational data assimilation Michael Tsyrulnikov and Alexander Rakitko HydroMetCenter of Russia College Park, 20 Oct 2014 Michael Tsyrulnikov and Alexander
More informationDynamic System Identification using HDMR-Bayesian Technique
Dynamic System Identification using HDMR-Bayesian Technique *Shereena O A 1) and Dr. B N Rao 2) 1), 2) Department of Civil Engineering, IIT Madras, Chennai 600036, Tamil Nadu, India 1) ce14d020@smail.iitm.ac.in
More informationA comparison of polynomial chaos and Gaussian process emulation for uncertainty quantification in computer experiments
University of Exeter Department of Mathematics A comparison of polynomial chaos and Gaussian process emulation for uncertainty quantification in computer experiments Nathan Edward Owen May 2017 Supervised
More informationUncertainty Quantification of Radionuclide Release Models using Non-Intrusive Polynomial Chaos. Casper Hoogwerf
Uncertainty Quantification of Radionuclide Release Models using Non-Intrusive Polynomial Chaos. Casper Hoogwerf 1 Foreword This report presents the final thesis of the Master of Science programme in Applied
More informationAccepted Manuscript. SAMBA: Sparse approximation of moment-based arbitrary polynomial chaos. R. Ahlfeld, B. Belkouchi, F.
Accepted Manuscript SAMBA: Sparse approximation of moment-based arbitrary polynomial chaos R. Ahlfeld, B. Belkouchi, F. Montomoli PII: S0021-9991(16)30151-6 DOI: http://dx.doi.org/10.1016/j.jcp.2016.05.014
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 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 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 informationGradient Enhanced Universal Kriging Model for Uncertainty Propagation in Nuclear Engineering
Gradient Enhanced Universal Kriging Model for Uncertainty Propagation in Nuclear Engineering Brian A. Lockwood 1 and Mihai Anitescu 2 1 Department of Mechanical Engineering University of Wyoming 2 Mathematics
More informationModel Calibration under Uncertainty: Matching Distribution Information
Model Calibration under Uncertainty: Matching Distribution Information Laura P. Swiler, Brian M. Adams, and Michael S. Eldred September 11, 008 AIAA Multidisciplinary Analysis and Optimization Conference
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 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 informationEnabling Advanced Automation Tools to manage Trajectory Prediction Uncertainty
Engineering, Test & Technology Boeing Research & Technology Enabling Advanced Automation Tools to manage Trajectory Prediction Uncertainty ART 12 - Automation Enrique Casado (BR&T-E) enrique.casado@boeing.com
More information4. Numerical Quadrature. Where analytical abilities end... continued
4. Numerical Quadrature Where analytical abilities end... continued Where analytical abilities end... continued, November 30, 22 1 4.3. Extrapolation Increasing the Order Using Linear Combinations Once
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 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 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 informationWhat Should We Do with Radioactive Waste?
What Should We Do with Radioactive Waste? Andrew Cliffe School of Mathematical Sciences University of Nottingham Simulation of Flow in Porous Media and Applications in Waste Management and CO 2 Sequestration
More informationGiovanni Migliorati. MATHICSE-CSQI, École Polytechnique Fédérale de Lausanne
Analysis of the stability and accuracy of multivariate polynomial approximation by discrete least squares with evaluations in random or low-discrepancy point sets Giovanni Migliorati MATHICSE-CSQI, École
More informationApplication and validation of polynomial chaos methods to quantify uncertainties in simulating the Gulf of Mexico circulation using HYCOM.
Application and validation of polynomial chaos methods to quantify uncertainties in simulating the Gulf of Mexico circulation using HYCOM. Mohamed Iskandarani Matthieu Le Hénaff Carlisle Thacker University
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 informationSpectral Representation of Random Processes
Spectral Representation of Random Processes Example: Represent u(t,x,q) by! u K (t, x, Q) = u k (t, x) k(q) where k(q) are orthogonal polynomials. Single Random Variable:! Let k (Q) be orthogonal with
More informationAnalytical Non-Linear Uncertainty Propagation: Theory And Implementation
Analytical Non-Linear Uncertainty Propagation: Theory And Implementation Kohei Fujimoto 1), Daniel J. Scheeres 2), and K. Terry Alfriend 1) May 13, 2013 1) Texas A&M University 2) The University of Colorado
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 informationFast Numerical Methods for Stochastic Computations: A Review
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 5, No. 2-4, pp. 242-272 Commun. Comput. Phys. February 2009 REVIEW ARTICLE Fast Numerical Methods for Stochastic Computations: A Review Dongbin Xiu Department
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 informationUncertainty Quantification and Validation Using RAVEN. A. Alfonsi, C. Rabiti. Risk-Informed Safety Margin Characterization. https://lwrs.inl.
Risk-Informed Safety Margin Characterization Uncertainty Quantification and Validation Using RAVEN https://lwrs.inl.gov A. Alfonsi, C. Rabiti North Carolina State University, Raleigh 06/28/2017 Assumptions
More informationSequential approaches to reliability estimation and optimization based on kriging
Sequential approaches to reliability estimation and optimization based on kriging Rodolphe Le Riche1,2 and Olivier Roustant2 1 CNRS ; 2 Ecole des Mines de Saint-Etienne JSO 2012, ONERA Palaiseau 1 Intro
More informationHyperbolic Polynomial Chaos Expansion (HPCE) and its Application to Statistical Analysis of Nonlinear Circuits
Hyperbolic Polynomial Chaos Expansion HPCE and its Application to Statistical Analysis of Nonlinear Circuits Majid Ahadi, Aditi Krishna Prasad, Sourajeet Roy High Speed System Simulations Laboratory Department
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 informationLong Time Propagation of Stochasticity by Dynamical Polynomial Chaos Expansions. Hasan Cagan Ozen
Long Time Propagation of Stochasticity by Dynamical Polynomial Chaos Expansions Hasan Cagan Ozen Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy in the Graduate
More information