Adaptive Collocation with Kernel Density Estimation
|
|
- Morgan Norman
- 5 years ago
- Views:
Transcription
1 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 University of Maryland at College Park with
2 Examples of 1 Background and Goals Examples of 5 with
3 Examples of Background and Goals Example: The Diffusion Equation The Joint Density Function Collocation Methods Example: the diffusion equation (a(x, ω) u) = f on D R d with suitable boundary conditions Uncertainty / randomness a = a(x, ω) is a random field For each fixed x D, a(x, ω) is a random variable Assumptions Coercivity: 0 < α 1 a α 2 < Problem is well posed Spatial correlation of random field: For x, x 1, x 2 D Mean µ(x) = E(a(x, )) Variance σ(x) 2 = E(a(x, ) 2 ) µ 2 Covariance function c(x 1, x 2 ) = E([a(x 1, ) µ(x 1 )] [a(x 2, ) µ(x 2 )]) is finite with
4 Examples of Example: The Diffusion Equation The Joint Density Function Collocation Methods Aims: Efficient computation of statistical properties of solution u(x, ω) Moments mean µ = Ω u(x, ω) dp(ω) = u(x, ξ) ρ(ξ)dξ Γ(Ω) variance Ω u(x, ω)2 dp(ω) µ 2 = Γ(Ω) u(x, ξ)2 ρ(ξ)dξ µ 2 Cumulative distribution functions, e.g., P(u(x, ω)) > α Issues addressed in this study: Costs of collocation methods Lack of knowledge of density function ρ(ξ) with
5 Examples of The Joint Density Function Example: The Diffusion Equation The Joint Density Function Collocation Methods Possible sources of diffusion coefficient Obtained from samples of a(x, ) Assume covariance function for a(x, ) is of standard type E.g. c(x 1, x 2 ) = exp ( x 1 x 2 1 /l) c(x 1, x 2 ) = exp ( ( x 1 x 2 2 /l) 2) Use to generate (truncated) Karhunen-Loève expansion a(x, ξ) = a 0 + σ m r=1 λr a r (x) ξ r Generate covariance matrix from samples, then as above Model itself has form of KL-expansion In all these scenarios The joint density function ρ(ξ) is generally not known with
6 Examples of Collocation Methods Example: The Diffusion Equation The Joint Density Function Collocation Methods Monte-Carlo (sampling) method: find u HE 1 (D) s.t. a(x, ξ k ) u vdx for all v HE 1 0 (D) D for a collection of samples {ξ k } L 2 (Γ) Collocation (Xiu, Hesthaven, Babuška, Nobile, Tempone, Webster) Choose {ξ k } in a special way (sparse grids), then construct discrete solution u (hp) (x, ξ) S E h T (p) to interpolate {u h (x, ξ k )} Advantages (vs. stochastic Galerkin): decouples algebraic system (like MC) applies in a straightforward way to nonlinear random terms vs. Monte Carlo: faster convergence for moderate number of parameters with
7 Examples of But Collocation may still be costly Example: The Diffusion Equation The Joint Density Function Collocation Methods Compare with Stochastic Galerkin (E., Miller, Phipps, Tuminaro): For comparable accuracy: # stochastic dof (collocation) 2 p (# stochastic dof (Galerkin)) { total polynomial degree for Galerkin Here p = level (depth) of sparse grid for collocation m = 4 uniform density Accuracy p = 1 p = 2 p = 3 p = 4 p = 5 Degrees of freedom p = 6 p = 5 p = 4 m = 5 uniform density Time p = 3 Error p = 2 p = 1 with
8 Examples of Adaptive Sparse Grid Collocation Influence of Density Function 1 Background and Goals 2 Adaptive Sparse Grid Collocation Influence of Density Function 3 4 Examples of 5 with
9 Examples of Adaptive Sparse Grid Collocation Adaptive Sparse Grid Collocation Influence of Density Function Collocation based on interpolation Consider function u(ξ) at left piecewise linear interpolant u I (ξ) ξ = child, node at next level uniform grid w(ξ ) = interpolation error at ξ Ma & Zabaras: If u i (ξ) is represented using hierarchical basis (left), then w(ξ ) = coefficient of basis function of interpolant on refined grid = hierarchical surplus Strategy: Refine grid using child ξ iff w(ξ ) > tolerance with
10 Examples of Adaptive Sparse Grid Collocation Adaptive Sparse Grid Collocation Influence of Density Function Collocation based on interpolation Consider function u(ξ) at left piecewise linear interpolant u I (ξ) ξ = child, node at next level uniform grid w(ξ ) = interpolation error at ξ Ma & Zabaras: If u i (ξ) is represented using hierarchical basis (left), then w(ξ ) = coefficient of basis function of interpolant on refined grid = hierarchical surplus Strategy: Refine grid using child ξ iff w(ξ ) > tolerance with
11 Examples of Adaptive Sparse Grid Collocation Influence of Density Function To Use this Idea with Diffusion Equation (a(x, ξ) u) = f Levels Start with ξ (0), compute solution u(x, ξ (0) ) Determines interpolant [A (0) u](x, ξ) Identify children {ξ (1) j } of ξ (0) Compute solutions {u(x, ξ (1) j )}, add {ξ (1) j } to set of collocation points according to adaptive strategy with tolerance test wj k > τ Determines interpolant [A (1) u](x, ξ) Repeat: identify children of level-1 points, compute solutions, etc. Result: Collocation solution A (k) u Approximate moments, distributions of u using A (k) u with
12 Examples of Adaptive Sparse Grid Collocation Influence of Density Function Algorithm: Adaptive interpolation with hierarchical basis functions Set A 0 (u)(ξ) = 0 Set k = 1 Set θadaptive 1 = θ1 repeat θ k+1 adaptive = for ξ k j w k j θadaptive k = u(ξj k ) A k 1 (u)(ξj k ) Hierarchical surplus if wj k > τ then θ k+1 adaptive = θk+1 adaptive child(ξ k endif endfor Set A k (u)(ξ) = k k = k + 1 until max( w k 1 j ) < τ i=1 j w i j ψi j (ξ) j ) Augment interpolant refined grid nodes only where interpolation error is large Interpolant = sum(levels) all level-k basis functions with
13 Examples of Adaptive Sparse Grid Collocation Influence of Density Function Will only add children of a collocation point if w k j > τ For an oscillatory function (below): points in refinement would tend to lie in oscillatory region u(ξ) = { ξ sin(1/ξ) ξ 0 0 ξ = 0 ξ with
14 Examples of Influence of Density Function Adaptive Sparse Grid Collocation Influence of Density Function What if density function is small on the oscillatory region? There may be no need to refine in this region u(ξ) = { ξ sin(1/ξ) ξ 0 0 ξ = 0 ξ Change criterion for refinement: refine if wj k ρ(ξj i ) is large with
15 Examples of Informally: only add children of a point if (i) the functon is changing in the region of that point, and (ii) the region is statistically significant What if we don t know ρ? Approximate ρ from the data using kernel density estimation Given N samples {ξ (i) } drawn from the distribution of ξ ˆρ(ξ) = 1 ( ) N Nh M k=1 K ξ ξ (i) h The kernel K satisfies R M K(ξ)dξ = 1 R M K(ξ)ξdξ = 0 R M K(ξ) ξ 2 dξ = k < K(ξ) 0 Comment: From here on: {ξ (i) } nearly always refers to samples with
16 Examples of KDE approximates the density by placing a bump over each observation Parameter h (bandwidth) controls smoothness Oversmoothed h too large Undersmoothed h too small with
17 Examples of In our examples: for bandwidth h, use maximum likelihood cross-validation: maximize CV (h) 1 N log(ˆρ i (ξ (i) )) N where ˆρ i (ξ) = 1 Nh M i=1 N k=1,k i ( ) ξ ξ (k) K h For K(ξ), use Epanechnikov kernel ( ) M 3 M K(ξ) = (1 ξi 2 )1 { 1 ξi 1} 4 i=1 Compactly supported, helps ensure coercivity for diffusion problem KDE material from Silverman with
18 Examples of One-Parameter Example Two-Parameter Example 1 Background and Goals Examples of One-Parameter Example Two-Parameter Example 5 with
19 Examples of One-Parameter Example Two-Parameter Example One-Parameter Example from Above Adaptive criterion and KDE: refine when w k j ˆρ(ξj i ) > τ: Collocation points are added only in statistically significant portions of the parameter space, determined from samples ξ with
20 Examples of Two-Parameter Example One-Parameter Example Two-Parameter Example Two dependent random variables, two statistically significant regions u(ξ 1, ξ 2 ) = ξ 1 ξ e (1 r) e (1 ˆr)2.1 r(ξ 1, ξ 2 ) = ξ ξ2 2 ˆr(ξ 1, ξ 2 ) = (ξ 1 5) 2 + (ξ 2 5) 2 ξ 1 log-normal with marginal density function ρ 1 (ξ 1 ) = ξ 2 = ξ 1 + η, η uniform (0, 1) (log ξ 1 1 ) 2 ξ 1 2π(.42 ) e 2(.4 2 ) with
21 Examples of One-Parameter Example Two-Parameter Example Plot of this function: Large near ξ = 1 and ξ 5 = 1 with
22 Examples of One-Parameter Example Two-Parameter Example Experiment Take 10,000 samples of (ξ 1, ξ 2 ) Shown at right, above Generate kernel density estimate ˆρ(ξ) by the method above Contour plot of ˆρ, shown at right below Observation: More samples near (0, 0) than (5, 5) = that large component plays more of a role in KDE with
23 Examples of One-Parameter Example Two-Parameter Example Experiment, continued So far: samples have been generated and used to define ˆρ Next: perform adaptive collocation with different values of τ Collocation points for τ = 10 3, τ = 10 4 shown at right N.B. Needed 5 global refinements to see large parts of u with
24 Examples of One-Parameter Example Two-Parameter Example Experiment, continued Use N = 10K samples {ξ (i) } N i=1 to compute sample mean E(u) µ MC 1 N N i=1 u(ξ(i) ) Compute sample variance in analogous manner From these: can get Monte-Carlo error with 95% confidence bound E(u) µ MC = Alterative: estimate mean using the collocation solution A(u) E(u) µ coll 1 N N i=1 (Au)(ξ(i) ) evaluated at Monte Carlo samples Comparison τ µ MC µ coll # coll. pts. (73) (101) (188) (3984) with
25 Examples of Error Analysis Example: Stochastic PDE 1 Background and Goals Examples of 5 Error Analysis Example: Stochastic PDE with
26 Error Analysis Background and Goals Examples of Error Analysis Example: Stochastic PDE Consider Chebyshev inequality error for Monte Carlo sample mean: That is: For a > 0, P ( E(u) µ MC a) Var(u) Take, say, P =.05, get a Var(u).05N N a 2. With 95% confidence, error of MC sample mean is bounded by Var(u).05N Collocation: Estimates the mean by µ coll = 1 N N i=1 (Au)(ξ(i) ) Error: E(u) µ coll E(u) µ MC + µ MC µ coll First term: MC error as above Second term: 1 N N i=1 (u(ξ(i) ) (Au)(ξ (i) )) Mean interpolation error with
27 Examples of Example: Stochastic PDE Error Analysis Example: Stochastic PDE Test problem: d dx a M(x, ξ) d dx u(x, ξ) = 1 x (0, 1) u(0, ξ) = u(1, ξ) = 0 a M = µ + M/2 1 k=0 λ k (ξ 2k cos(2πkx) + ξ 2k+1 sin(2πkx)) µ = 3, λ k = exp( k) ξ k uniformly distributed on [0, 1] Experiment: As above Generate N samples of ξ Use them to generate estimate ˆρ(ξ) Use ˆρ to generate adaptive collocation solution Au Compare with: Use the same N samples to perform MC simulation Differences from above: M > 2 random variables Function is solution to PDE with
28 Examples of Error Analysis Example: Stochastic PDE Representative Results: M = 10 parameters N τ (76) (1026) (1655) (5026) (8111) (92) (1170) (1189) (5773) (9404) (59) (1216) (1989) (5996) (9664) (93) (1120) (2041) (6095) (9787) (93) (1187) (2127) (6050) (9942) Monte-Carlo Mean interpolation error error 1 N N i=1 u(x, ξ(i) ) A(u)(x, ξ (i) ) l 2 (D) Parens: Number of DE solves needed for collocation Green: This number is smaller than N and yields smaller error with
29 Examples of Error Analysis Example: Stochastic PDE For different parameter sizes: For a smaller number, M = 5 Adaptive collocation is more effective N τ (28) (212) (301) (813) (1169) (28) (211) (315) (777) (1210) (33) (200) (297) (762) (1207) (33) (172) (286) (745) (1104) (33) (180) (294) (780) (1107) Monte-Carlo Mean interpolation error error 1 N N i=1 u(x, ξ(i) ) A(u)(x, ξ (i) ) l 2 (D) with
30 Examples of Error Analysis Example: Stochastic PDE For larger M = 20: Collocation is still effective N τ (41) (878) (1299) (4126) (6958) (41) (1045) (1738) (5545) (9106) (119) (1618) (2622) (8580) (14012) (156) (2459) (4169) (13389) (22276) (193) (3108) (4991) (15963) (26081) Monte-Carlo Mean interpolation error error 1 N N i=1 u(x, ξ(i) ) A(u)(x, ξ (i) ) l 2 (D) with
31 Examples of Concluding Remarks Adaptive collocation with kernel density estimation Enables adaptive construction of collocation points Avoids need to know density function explicitly Appears to be competitive with, potentially superior to, Monte Carlo methods with
Stochastic 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 informationSolving the steady state diffusion equation with uncertainty Final Presentation
Solving the steady state diffusion equation with uncertainty Final Presentation Virginia Forstall vhfors@gmail.com Advisor: Howard Elman elman@cs.umd.edu Department of Computer Science May 6, 2012 Problem
More informationASSESSMENT OF COLLOCATION AND GALERKIN APPROACHES TO LINEAR DIFFUSION EQUATIONS WITH RANDOM DATA
International Journal for Uncertainty Quantification, 11):19 33, 2011 ASSESSMENT OF COLLOCATION AND GALERKIN APPROACHES TO LINEAR DIFFUSION EQUATIONS WITH RANDOM DATA Howard C. Elman, 1, Christopher W.
More 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 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 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 information. Frobenius-Perron Operator ACC Workshop on Uncertainty Analysis & Estimation. Raktim Bhattacharya
.. Frobenius-Perron Operator 2014 ACC Workshop on Uncertainty Analysis & Estimation Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. uq.tamu.edu
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 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 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 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 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 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 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 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 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 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 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 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 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 informationRandom Matrix Eigenvalue Problems in Probabilistic Structural Mechanics
Random Matrix Eigenvalue Problems in Probabilistic Structural Mechanics S Adhikari Department of Aerospace Engineering, University of Bristol, Bristol, U.K. URL: http://www.aer.bris.ac.uk/contact/academic/adhikari/home.html
More informationAN ADAPTIVE REDUCED BASIS COLLOCATION METHOD BASED ON PCM ANOVA DECOMPOSITION FOR ANISOTROPIC STOCHASTIC PDES
AN ADAPTIVE REDUCED BASIS COLLOCATION METHOD BASED ON PCM ANOVA DECOMPOSITION FOR ANISOTROPIC STOCHASTIC PDES Heyrim Cho & Howard Elman 2 Department of Mathematics, University of Maryland, College Park,
More informationarxiv: v1 [math.na] 14 Sep 2017
Stochastic collocation approach with adaptive mesh refinement for parametric uncertainty analysis arxiv:1709.04584v1 [math.na] 14 Sep 2017 Anindya Bhaduri a, Yanyan He 1b, Michael D. Shields a, Lori Graham-Brady
More 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 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 informationDownloaded 01/28/13 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 27, No. 3, pp. 1118 1139 c 2005 Society for Industrial and Applied Mathematics HIGH-ORDER COLLOCATION METHODS FOR DIFFERENTIAL EQUATIONS WITH RANDOM INPUTS DONGBIN XIU AND JAN
More informationRandom Eigenvalue Problems Revisited
Random Eigenvalue Problems Revisited S Adhikari Department of Aerospace Engineering, University of Bristol, Bristol, U.K. Email: S.Adhikari@bristol.ac.uk URL: http://www.aer.bris.ac.uk/contact/academic/adhikari/home.html
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 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 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 sparse grid stochastic collocation method for elliptic partial differential equations with random input data
A sparse grid stochastic collocation method for elliptic partial differential equations with random input data F. Nobile R. Tempone C. G. Webster June 22, 2006 Abstract This work proposes and analyzes
More 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 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 informationPractical Bayesian Optimization of Machine Learning. Learning Algorithms
Practical Bayesian Optimization of Machine Learning Algorithms CS 294 University of California, Berkeley Tuesday, April 20, 2016 Motivation Machine Learning Algorithms (MLA s) have hyperparameters that
More informationFully Bayesian Deep Gaussian Processes for Uncertainty Quantification
Fully Bayesian Deep Gaussian Processes for Uncertainty Quantification N. Zabaras 1 S. Atkinson 1 Center for Informatics and Computational Science Department of Aerospace and Mechanical Engineering University
More informationMultilevel stochastic collocations with dimensionality reduction
Multilevel stochastic collocations with dimensionality reduction Ionut Farcas TUM, Chair of Scientific Computing in Computer Science (I5) 27.01.2017 Outline 1 Motivation 2 Theoretical background Uncertainty
More informationA sparse grid stochastic collocation method for partial differential equations with random input data
A sparse grid stochastic collocation method for partial differential equations with random input data F. Nobile R. Tempone C. G. Webster May 24, 2007 Abstract This work proposes and analyzes a Smolyak-type
More 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 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 informationTHE CHOICE OF AUXILIARY DENSITY FUNCTION IN STOCHASTIC COLLOCATION
IN STOCHASTIC COLLOCATION KELLIN RUMSEY & NEVIN MARTIN 12/06/2016 1 IN STOCHASTIC COLLOCATION 2 1. Introduction Stochastic collocation has been used in a wide range of application spaces as a way to efficiently
More 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 informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationFinite difference method for elliptic problems: I
Finite difference method for elliptic problems: I Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
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 informationBath Institute For Complex Systems
BICS Bath Institute for Complex Systems Finite Element Error Analysis of Elliptic PDEs with Random Coefficients and its Application to Multilevel Monte Carlo Methods Julia Charrier, Robert Scheichl and
More informationHigh Frequency Scattering by Convex Polygons Stephen Langdon
Bath, October 28th 2005 1 High Frequency Scattering by Convex Polygons Stephen Langdon University of Reading, UK Joint work with: Simon Chandler-Wilde Steve Arden Funded by: Leverhulme Trust University
More informationWELL POSEDNESS OF PROBLEMS I
Finite Element Method 85 WELL POSEDNESS OF PROBLEMS I Consider the following generic problem Lu = f, where L : X Y, u X, f Y and X, Y are two Banach spaces We say that the above problem is well-posed (according
More informationCheng Soon Ong & Christian Walder. Canberra February June 2018
Cheng Soon Ong & Christian Walder Research Group and College of Engineering and Computer Science Canberra February June 2018 Outlines Overview Introduction Linear Algebra Probability Linear Regression
More informationHigh Frequency Scattering by Convex Polygons Stephen Langdon
High Frequency Scattering by Convex Polygons Stephen Langdon University of Reading, UK Joint work with: Simon Chandler-Wilde, Steve Arden Funded by: Leverhulme Trust University of Maryland, September 2005
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 Quantification and hypocoercivity based sensitivity analysis for multiscale kinetic equations with random inputs.
Uncertainty Quantification and hypocoercivity based sensitivity analysis for multiscale kinetic equations with random inputs Shi Jin University of Wisconsin-Madison, USA Shanghai Jiao Tong University,
More informationLow-rank techniques applied to moment equations for the stochastic Darcy problem with lognormal permeability
Low-rank techniques applied to moment equations for the stochastic arcy problem with lognormal permeability Francesca Bonizzoni 1,2 and Fabio Nobile 1 1 CSQI-MATHICSE, EPFL, Switzerland 2 MOX, Politecnico
More 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 informationDimensionality reduction and polynomial chaos acceleration of Bayesian inference in inverse problems
Dimensionality reduction and polynomial chaos acceleration of Bayesian inference in inverse problems The MIT Faculty has made this article openly available. Please share how this access benefits you. Your
More informationAN EFFICIENT COMPUTATIONAL FRAMEWORK FOR UNCERTAINTY QUANTIFICATION IN MULTISCALE SYSTEMS
AN EFFICIENT COMPUTATIONAL FRAMEWORK FOR UNCERTAINTY QUANTIFICATION IN MULTISCALE SYSTEMS A Dissertation Presented to the Faculty of the Graduate School of Cornell University in Partial Fulfillment of
More 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 informationUNIVERSITY OF CALIFORNIA, SAN DIEGO. An Empirical Chaos Expansion Method for Uncertainty Quantification
UNIVERSITY OF CALIFORNIA, SAN DIEGO An Empirical Chaos Expansion Method for Uncertainty Quantification A Dissertation submitted in partial satisfaction of the requirements for the degree Doctor of Philosophy
More informationOptimal Prediction for Radiative Transfer: A New Perspective on Moment Closure
Optimal Prediction for Radiative Transfer: A New Perspective on Moment Closure Benjamin Seibold MIT Applied Mathematics Mar 02 nd, 2009 Collaborator Martin Frank (TU Kaiserslautern) Partial Support NSF
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationCBS Constants and Their Role in Error Estimation for Stochastic Galerkin Finite Element Methods. Crowder, Adam J and Powell, Catherine E
CBS Constants and Their Role in Error Estimation for Stochastic Galerkin Finite Element Methods Crowder, Adam J and Powell, Catherine E 2017 MIMS EPrint: 2017.18 Manchester Institute for Mathematical Sciences
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 informationUncertainty Quantification for multiscale kinetic equations with random inputs. Shi Jin. University of Wisconsin-Madison, USA
Uncertainty Quantification for multiscale kinetic equations with random inputs Shi Jin University of Wisconsin-Madison, USA Where do kinetic equations sit in physics Kinetic equations with applications
More informationStochastic collocation and mixed finite elements for flow in porous media
Stochastic collocation and mixed finite elements for flow in porous media Benjamin Ganis a, Hector Klie b, Mary F. Wheeler b, Tim Wildey b, Ivan Yotov a, and Dongxiao Zhang c a Department of Mathematics,
More 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 informationDecay profiles of a linear artificial viscosity system
Decay profiles of a linear artificial viscosity system Gene Wayne, Ryan Goh and Roland Welter Boston University rwelter@bu.edu July 2, 2018 This research was supported by funding from the NSF. Roland Welter
More informationSimultaneous Robust Design and Tolerancing of Compressor Blades
Simultaneous Robust Design and Tolerancing of Compressor Blades Eric Dow Aerospace Computational Design Laboratory Department of Aeronautics and Astronautics Massachusetts Institute of Technology GTL Seminar
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 Particle Methods for Rarefied Gases
CCES Seminar WS 2/3 Stochastic Particle Methods for Rarefied Gases Julian Köllermeier RWTH Aachen University Supervisor: Prof. Dr. Manuel Torrilhon Center for Computational Engineering Science Mathematics
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 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 informationLiquid-Rocket Transverse Triggered Combustion Instability: Deterministic and Stochastic Analyses
Liquid-Rocket Transverse Triggered Combustion Instability: Deterministic and Stochastic Analyses by W. A. Sirignano Mechanical and Aerospace Engineering University of California, Irvine Collaborators:
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 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 informationSupplemental Material for KERNEL-BASED INFERENCE IN TIME-VARYING COEFFICIENT COINTEGRATING REGRESSION. September 2017
Supplemental Material for KERNEL-BASED INFERENCE IN TIME-VARYING COEFFICIENT COINTEGRATING REGRESSION By Degui Li, Peter C. B. Phillips, and Jiti Gao September 017 COWLES FOUNDATION DISCUSSION PAPER NO.
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 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 informationA multiscale preconditioner for stochastic mortar mixed finite elements
A multiscale preconditioner for stochastic mortar mixed finite elements Mary F. Wheeler b, Tim Wildey b, and Ivan Yotov a a Department of Mathematics, University of Pittsburgh, Pittsburgh, PA 15260 b The
More informationA Kernel-based Collocation Method for Elliptic Partial Differential Equations with Random Coefficients
A Kernel-based Collocation Method for Elliptic Partial Differential Equations with Random Coefficients Gregory E. Fasshauer and Qi Ye Abstract This paper is an extension of previous work where we laid
More informationSparse Grids. Léopold Cambier. February 17, ICME, Stanford University
Sparse Grids & "A Dynamically Adaptive Sparse Grid Method for Quasi-Optimal Interpolation of Multidimensional Analytic Functions" from MK Stoyanov, CG Webster Léopold Cambier ICME, Stanford University
More 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 MULTILEVEL STOCHASTIC COLLOCATION METHOD FOR PARTIAL DIFFERENTIAL EQUATIONS WITH RANDOM INPUT DATA
A MULTILEVEL STOCHASTIC COLLOCATION METHOD FOR PARTIAL DIFFERENTIAL EQUATIONS WITH RANDOM INPUT DATA A. L. TECKENTRUP, P. JANTSCH, C. G. WEBSTER, AND M. GUNZBURGER Abstract. Stochastic collocation methods
More informationReliability Monitoring Using Log Gaussian Process Regression
COPYRIGHT 013, M. Modarres Reliability Monitoring Using Log Gaussian Process Regression Martin Wayne Mohammad Modarres PSA 013 Center for Risk and Reliability University of Maryland Department of Mechanical
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 informationHierarchical Matrices. Jon Cockayne April 18, 2017
Hierarchical Matrices Jon Cockayne April 18, 2017 1 Sources Introduction to Hierarchical Matrices with Applications [Börm et al., 2003] 2 Sources Introduction to Hierarchical Matrices with Applications
More informationMinimal Element Interpolation in Functions of High-Dimension
Minimal Element Interpolation in Functions of High-Dimension J. D. Jakeman, A. Narayan, D. Xiu Department of Mathematics Purdue University West Lafayette, Indiana Random Model as a function We can consider
More informationStochastic Capacitance Extraction Considering Process Variation with Spatial Correlation
Stochastic Capacitance Extraction Considering Process Variation with Spatial Correlation Tuck Chan EE Department UCLA, CA 90095 tuckie@ucla.edu Fang Gong EE Department UCLA, CA 90095 gongfang@ucla.edu
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 informationEstimation of uncertainties using the Guide to the expression of uncertainty (GUM)
Estimation of uncertainties using the Guide to the expression of uncertainty (GUM) Alexandr Malusek Division of Radiological Sciences Department of Medical and Health Sciences Linköping University 2014-04-15
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 informationA Primal-Dual Weak Galerkin Finite Element Method for Second Order Elliptic Equations in Non-Divergence Form
A Primal-Dual Weak Galerkin Finite Element Method for Second Order Elliptic Equations in Non-Divergence Form Chunmei Wang Visiting Assistant Professor School of Mathematics Georgia Institute of Technology
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 informationKullback-Leibler Designs
Kullback-Leibler Designs Astrid JOURDAN Jessica FRANCO Contents Contents Introduction Kullback-Leibler divergence Estimation by a Monte-Carlo method Design comparison Conclusion 2 Introduction Computer
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 informationModel Specification Testing in Nonparametric and Semiparametric Time Series Econometrics. Jiti Gao
Model Specification Testing in Nonparametric and Semiparametric Time Series Econometrics Jiti Gao Department of Statistics School of Mathematics and Statistics The University of Western Australia Crawley
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 informationA data-driven stochastic collocation approach for uncertainty quantification in MEMS
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 21; 83:575 597 Published online 9 March 21 in Wiley InterScience (www.interscience.wiley.com)..2844 A data-driven stochastic
More informationThe Bayesian approach to inverse problems
The Bayesian approach to inverse problems Youssef Marzouk Department of Aeronautics and Astronautics Center for Computational Engineering Massachusetts Institute of Technology ymarz@mit.edu, http://uqgroup.mit.edu
More informationAlberto Bressan. Department of Mathematics, Penn State University
Non-cooperative Differential Games A Homotopy Approach Alberto Bressan Department of Mathematics, Penn State University 1 Differential Games d dt x(t) = G(x(t), u 1(t), u 2 (t)), x(0) = y, u i (t) U i
More informationNonparametric Regression. Badr Missaoui
Badr Missaoui Outline Kernel and local polynomial regression. Penalized regression. We are given n pairs of observations (X 1, Y 1 ),...,(X n, Y n ) where Y i = r(x i ) + ε i, i = 1,..., n and r(x) = E(Y
More informationMath 660-Lecture 15: Finite element spaces (I)
Math 660-Lecture 15: Finite element spaces (I) (Chapter 3, 4.2, 4.3) Before we introduce the concrete spaces, let s first of all introduce the following important lemma. Theorem 1. Let V h consists of
More information