NONLOCALITY AND STOCHASTICITY TWO EMERGENT DIRECTIONS FOR APPLIED MATHEMATICS. Max Gunzburger
|
|
- Bryan Warren
- 5 years ago
- Views:
Transcription
1 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 s C cientific omputing Florida State University 1/2748
2 STOCHASTICITY
3 RECENT DEVELOPMENTS IN NUMERICAL METHODS FOR SPDES Lots about uncertainty in this workshop: Brian Adams Managing Complexity in Simulation-Based Uncertainty Quantification Dan Cacuci Experimentally Validated Best Estimate + Uncertainty Modeling of Complex Systems: The Cornerstone of the Emerging Field of Predictive Science Chris Jones Mathematical Challenges of Climate Research: Data Assimilation and Uncertainty
4 The ubiquity of uncertainty quantification (UQ) there is uncertainty everywhere uncertainty should be quantified, or else the world will come to an end but deterministic models have served us quite well in lots of situations but knowing about uncertainty might be able to avert disastrous situations, might save money, and might end up with better policies at least that is the hope (or dream) one has about UQ
5 The ubiquity of probability common definition of UQ uncertainty quantification = given statistics about the inputs of a system, determine statistics about the output of the system but what is usually meant uncertainty quantification = given information about the pdf of the inputs of a system, determine information about the pdf of the output of the system
6 but, because there are other ways to quantify uncertainty, what it should really be uncertainty quantification = given information about the uncertainty in the inputs of a system, determine information about the uncertainty in the output of the system to a statistician UQ statistics
7 a word about randomness computations are always deterministic pseudo random number generators are deterministic algorithms that s why computations are repeatable lots of probability and statistics theorems to not apply a word about rare events they are not always rare - just because an event happening has probability zero, does not mean it doesn t happen
8 The ubiquity of SPDEs everything in the universe can be described by stochastic PDEs what we mean is, perhaps, everything in the world can be described by PDEs with random inputs if this were only true!!!!! to a probabilist SPDE PDE with random inputs
9 But in this talk, we define stochasticity = probabilistic approach to UQ for systems governed by PDEs SPDE = PDE with random inputs within this talk, we give a short course on numerical methods for SPDEs in particular, we will lecture on recent developments but first, we need to talk about random inputs in particular, we will go over what is assumed (justly or unjustly) about random inputs before we get into the short course about recent developments assumptions marked with a * are assumed for old developments as well
10 The only thing certain about uncertainty is that it is itself uncertain This is certainly true about uncertain inputs especially those that are defined probabilistically
11 The ubiquity of independence* random parameters are assumed to be independent - makes life simple - may or may not be true in practice
12 The ubiquity of hypercubes* random parameters are assumed to be defined on a hypercube (of possible infinite extent) can t easily do constrained parameters, e.g., γ 2 + β * For all practical purposes this is also true for old developments if one has more than a few parameters
13 The ubiquity of Gauss or uniform* when in doubt, assume the inputs have a uniform or Gaussian PDF - if you know the mean or if you know a range (or impose one) for a random parameter, assume it is uniformly distributed - if you know the mean and variance, assume it s Guassian need for lots more model calibration determining statistical information about inputs
14 The ubiquity of making things positive* how can a coefficient in an elliptic equation be Gaussian? no problema assume it is e γ where γ is Gaussian (log-normal coefficient) why?
15 The ubiquity of correlation it is assumed that random fields are correlated colored noise this means one has to know something about how the random field is correlated one seldom knows how the random field is correlated, so, most often, it is assumed that correlation function = e (x x ) 2 where L = correlation length L 2 or e x x L A random field can be viewed as a function whose value at a point or at an instant of time is random
16 Now, we are ready for the short course on recent developments in numerical methods for SPDEs
17 A one-slide, very short course in recent developments in numerical methods for SPDEs SPDE = a partial differential equation with random inputs stochastic finite element method = discretize dependence of solution on spatial variables using a finite element method stochastic Galerkin method = discretize dependence of solution on random parameters using a Galerkin approach Karhunen-Loevy expansion = do an SVD of the correlation matrix to discretize a correlated random field into a finite set of uncorrelated parameters polynomial chaos method = approximate dependence of inputs and/or solutions on parameters using global orthogonal polynomials stochastic sampling method = sample random parameters at a selected set of sampling points, e.g., Monte Carlo sampling stochastic collocation method = sample random parameters at quadrature points for some quadrature rule sparse grid method = use a sparse grid quadrature rule, e.g., a Smolyak rule, to define the sampling points
18 What has to be assumed about the outputs for these recently developed methods to work? Remember, we already assumed a lot about the inputs for these recently developed methods to be applicable method working = is better than, e.g., Monte Carlo sampling
19 The ubiquity of the curse obtaining accurate statistics about the output of a system requires multiple simulations of the system or the simulation of a humongous coupled stochastic-spatial-temporal system recent developments in numerical methods for SPDEs are meant to mitigate the curse of dimensionality as the parameter dimension increases, the cost of obtaining accurate statistics grows quickly the idea is that by increasing accuracy, one can reduce costs e.g., the number of sample points can be reduced do polynomial chaos and sparse grid collocation methods really do this? in the best-case scenario, they do somewhat for moderate parameter dimension Very important fact: in principle, the accuracy of Monte Carlo sampling is independent of the parameter dimension
20 The ubiquity of smoothness the solution must be assumed to possess lots of derivatives with respect to the parameters quantities of interest are also assumed to possess lots of derivatives with respect to the parameters without the necessary smoothness, recently developed methods break down Often, the pointwise solution of the SPDE is not what one is interested in, but rather, a functional of that solution is of interest
21 THE UBIQUITY OF WHITE NOISE White noise random fields, are by far the most popular means for modeling uncertainty in inputs Gaussian white noise is by far the most popular white noise white noise is often invoked when one knows nothing about the nature of the uncertainty lots of beautiful math about white noise very popular in the modeling community The value of a white noise random field at a point or at an instant of time is identically distributed but is independent of and uncorrelated from its values at other points and other instants of time Recall that polynomial chaos and stochastic collocation only apply to correlated noise
22 Problemas numero uno: white noise has infinite energy the power spectrum of white noise is flat approximations of white noise such as 1 N t i=1 χ (ti 1,t i ]ω i where ω i, i = 1,..., N, are independently sampled from a given (usually Guassian) PDF, have unboundedly increasing energy as t 0 numero dos: white noise is not observed
23 Why do we get away with using white noise? one truncates white noise (band limit it to a finite spectrum) of course, then it is not longer white noise is is another model for noise at least for the additive noise case and for problems with smoothing solution operators, the solution is not white in fact, it can be much smoorter Is there an alternative to white noise? one that has finite energy one that is observed pink noise or more, generally, 1/f α noise
24 Noises having a 1/f α power spectrum α = 0 white noise α = 2 Brownian motion (or red noise) α = 1 pink noise finite energy for α > 1 ubiquitously observed in natural, social, man-made,... systems Is 1/f α noise the noise for you? it may not be what you want, but if you don t know anything about your noise, why not?
25 Realizations of approximate 1/f α noise for α = 0, 0.5, 1, 1.5, and 2
26 Corresponding approximate realizations of solution of 1D elliptic problem with 1/f α additive noise
27 NONLOCALITY
28 In many situations, atomistic and finer scale models invariably involve nonlocal interactions some interactions are of short range, e.g., effectively extending over a few neighboring atoms Leonard-Jones potential in molecular dynamics some interactions are of long range, e.g., effectively extending over many neighboring atoms Coulomb potential in molecular dynamics Classical continuum models, on the other hand, interactions are local in nature particles interact only through contact one of the truly seminal ideas in modeling the physical world (Cauchy)
29 Classical models have been extremely useful (this is an understatement!) however Classical models cannot accurately model or completely break down in many important situations plasticity, radiation, fracture,... Classical models, e.g., of elasticity, do not possess a length scale other than those characterizing the domain or that appear in spatially dependent coefficients, e.g., material properties think of the wave equation u tt u xx = f as a result, they cannot be used on their own for multiscale modeling gives rise to, e.g., AtC (atomistic-to-continuum) coupling models for which classical continuum models are coupled to atomistic models for example, classical elasticity coupled to molecular dynamics Thus, there is a huge industry devoted to fixing classical models
30 Higher-order gradient methods, e.g., au xx replaced by au xx bu xxxx now we have the length scale b/a now we have a multiscale model if b/a is small, then we reduce to the case u xx if b/a is large, then we reduce to the case u xxxx if b/a is neither small nor large, the we have u xx + a b u xxxx Nonlinear rheologies, e.g., nonlinear constitutive laws a length scale can be embedded into the constitutive law very popular in solid mechanics (e.g., plasticity), electromagnetism (e.g., ferromagnetics), fluids (polymers, mantle convection, or ice-sheet flows)
31 Nonlocal (integral) terms added to differential terms u xx replaced by u xx + k(x, y)u(y)dy a length scale can be easily be made to appear e.g., through an interaction radius for the kernel k(x, y) nonlinear integral terms also popular e.g., to model radiation effects well-used idea in time dependent problems e.g., to model problems with memory or delay also used for spatial modeling in all three types of fixes for classical models, still have spatial derivatives
32 Why we do not want spatial derivatives derivatives make it difficult to model singular behavior, e.g., crack propagation derivatives make it even more difficult to model nucleation of singularities, e.g., crack nucleation Why we do want nonlocality coupling (local) continuum models to (nonlocal) atomistic models is messy, e.g., ghost force effects especially at small scales, nonlocality is observed, e.g., in diffusion and solid mechanics experiments Why we do want multiscale having a model that accurately describes behaviors over a wide range of scales is the holy grail, or at least, the pot of gold at the end of the rainbow
33 Ideal model (perhaps) contains a length scale δ is a multiscale model when viewed on a scale of δ, it behaves like a nonlocal, discrete (e.g., atomistic) model when viewed on a scale much larger than δ, it behaves like a local, classical continuum model is free of spatial derivatives is itself a continuum model Such models exist for diffusion and mechanics lots of great math for nonlocal diffusion along with some simulations impressive simulations and some recent math for nonlocal mechanics
34 NONLOCAL DIFFUSION The parabolic heat equation modeling diffusion: u t (D(x) w(x) ) = b(x, t) in Ω R d The nonlocal equation modeling nonlocal diffusion: with u t + 2 ( u(x ) u(x) ) µ(x, x ) dx = b(x) in Ω R d B δ (x) µ(x, x ) = µ(x, x) 0 B δ (x) = {x Ω : x x δ}
35 suppose that and Ω Ω µ(x, x ) dx dx = 1 µ(x, x ) = µ(x, x) 0 so that µ(x, x ) can be interpreted as the joint probability density of moving between x and x then ( u(x ) u(x) ) µ(x, x ) dx B δ (x) = u(x )µ(x, x ) dx u(x) B δ (x) B δ (x) µ(x, x ) dx is the rate at which u enters x less the rate at which u departs x
36 the nonlocal diffusion equation can be derived from a nonlocal random walk you get drift if µ(x, x ) = µ(x, x) the nonlocal diffusion equation is an example of a differential Chapman-Kolmogorov equation
37 Constrained-value problems for the operator ( Lu = 2 u(x ) u(x) ) µ(x, x ) dx B δ (x) to obtain a solution of the nonlocal problem Lu(x) = b(x) data has to be specified on measurable sets thus, for example, we have the constrained-value problem { Lu(x) = b(x) for x Ω u(x) = g(x) where meas(ω o ) 0) and Ω Ω o = for x Ω o this is the nonlocal analog of the Dirichlet boundary-value problem for a second-order elliptic operator nonlocal analogs of Neumann problems can also be defined
38 What about spaces for well posedness of constrained-value problem? for appropriate (singular) kernels µ(x, x ) L : H s (Ω) H s (Ω) for 0 s 1 this is in contrast to second-order elliptic operators that map H 1 0(Ω) H 1 (Ω) this is also in contrast to fractional Laplacians that can map H s (Ω) H s (Ω) for s 1 2 note that for kernels such that s < 1 2, L is bounded when acting on functions with jump discontinuities note that the case s = 0 is allowable solving the constrained value problem results in no smoothing
39 Connections to differential operators if u(x) is smooth, then, L(u) second-order elliptic operator as δ 0 Nonlinear versions of nonlocal constrained-value problems are easily defined
40 NONLOCAL MECHANICS Generalizations of the nonlocal equation ρu tt = 2 µ(x, x ) ( u(x ) u(x) ) dx B δ (x) define nonlocal models for mechanics the peridynamics model for mechanics Everything said about the nonlocal diffusion operator holds for the operator 2 µ(x, x ) ( u(x ) u(x) ) dx B δ (x) and its peridynamic generalizations constrained-value problems analogous to boundary-value problems in classical elasticity can be defined reduced smoothness gains with respect to the data for solutions of constrainedvalue problems can again be obtained, including the L 2 L 2 case
41 for appropriately chosen kernels, the operator is bounded when acting on functions with jump discontinuities generalizations of the operators reduce to, as δ 0, to classical linear elasticity operators nonlinear versions can be defined Furthermore, peridynamics has been shown to reduce to molecular dynamics when viewed at the atomistic scale Furthermore, for appropriate kernels, displacements u having jump discontinuities are allowable cracks are allowable Generalizations are easily defined that also allow for the nucleation of cracks
42 Constrained-value problem in peridynamics may be discretized by a collocation method to obtain a discrete system that looks like the equations for a system of interacting particles in fact, peridynamics has been incorporated at Sandia into their molecular dynamics software package LAMMPS
43 Implications for finite element methods (based on weak formulations of constrainedvalue problems) discontinuous finite element spaces are conforming optimal accuracy can be obtained by using abrupt grid refinement, this remains true even if solutions contain jump discontinuities
44 a multiscale computational model is now realizable by merely varying the grid size h and choosing appropriate finite element spaces where nothing bad is happening, choose h δ standard continuous finite element spaces discrete equations are automatically the same as that for the finite element discretization of the corresponding PDE model where bad things are happening, choose h < δ discontinuous finite element spaces discrete equations automatically reduce to particle dynamics-type equations
45 Dispersion behavior
Fast 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 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 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 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 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 informationDinesh Kumar, Mehrdad Raisee and Chris Lacor
Dinesh Kumar, Mehrdad Raisee and Chris Lacor Fluid Mechanics and Thermodynamics Research Group Vrije Universiteit Brussel, BELGIUM dkumar@vub.ac.be; m_raisee@yahoo.com; chris.lacor@vub.ac.be October, 2014
More informationA NONLOCAL VECTOR CALCULUS WITH APPLICATION TO NONLOCAL BOUNDARY-VALUE PROBLEMS. Max Gunzburger
A NONLOCAL VECTOR CALCULUS WITH APPLICATION TO NONLOCAL BOUNDARY-VALUE PROBLEMS Max Gunzburger Department of Scientific Computing Florida State University Collaboration with: Richard Lehoucq Sandia National
More informationIntroduction and some preliminaries
1 Partial differential equations Introduction and some preliminaries A partial differential equation (PDE) is a relationship among partial derivatives of a function (or functions) of more than one variable.
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 informationUncertainty Quantification in Computational Science
DTU 2010 - Lecture I Uncertainty Quantification in Computational Science Jan S Hesthaven Brown University Jan.Hesthaven@Brown.edu Objective of lectures The main objective of these lectures are To offer
More informationA 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 informationA VECTOR CALCULUS AND FINITE ELEMENT METHODS FOR NONLOCAL DIFFUSION EQUATIONS. Max Gunzburger
A VECTOR CALCULUS AND FINITE ELEMENT METHODS FOR NONLOCAL DIFFUSION EQUATIONS Max Gunzburger Department of Scientific Computing Florida State University Collaboration with Richard Lehoucq Sandia National
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 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 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 informationConcepts in Global Sensitivity Analysis IMA UQ Short Course, June 23, 2015
Concepts in Global Sensitivity Analysis IMA UQ Short Course, June 23, 2015 A good reference is Global Sensitivity Analysis: The Primer. Saltelli, et al. (2008) Paul Constantine Colorado School of Mines
More informationSolving Boundary Value Problems (with Gaussians)
What is a boundary value problem? Solving Boundary Value Problems (with Gaussians) Definition A differential equation with constraints on the boundary Michael McCourt Division Argonne National Laboratory
More informationAn Introduction to Numerical Methods for Differential Equations. Janet Peterson
An Introduction to Numerical Methods for Differential Equations Janet Peterson Fall 2015 2 Chapter 1 Introduction Differential equations arise in many disciplines such as engineering, mathematics, sciences
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationPascal s Triangle on a Budget. Accuracy, Precision and Efficiency in Sparse Grids
Covering : Accuracy, Precision and Efficiency in Sparse Grids https://people.sc.fsu.edu/ jburkardt/presentations/auburn 2009.pdf... John Interdisciplinary Center for Applied Mathematics & Information Technology
More information13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs)
13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs) A prototypical problem we will discuss in detail is the 1D diffusion equation u t = Du xx < x < l, t > finite-length rod u(x,
More informationDirichlet s principle and well posedness of steady state solutions in peridynamics
Dirichlet s principle and well posedness of steady state solutions in peridynamics Petronela Radu Work supported by NSF - DMS award 0908435 January 19, 2011 The steady state peridynamic model Consider
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
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 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 informationContinuous and Discontinuous Finite Element Methods for a Peridynamics Model of Mechanics DRAFT
Continuous and Discontinuous Finite Element Methods for a Peridynamics Model of Mechanics Xi Chen and Max Gunzburger Department of Scientific Computing, Florida State University, Tallahassee FL 3306-40,
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 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 informationAN ALGORITHMIC INTRODUCTION TO NUMERICAL METHODS FOR PDES WITH RANDOM INPUTS DRAFT. Max Gunzburger
AN ALGORITHMIC INTRODUCTION TO NUMERICAL METHODS FOR PDES WITH RANDOM INPUTS Max Gunzburger Department of Scientific Computing Florida State University mgunzburger@fsu.edu with Clayton Webster John Burkardt
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 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 informationUtilizing Adjoint-Based Techniques to Improve the Accuracy and Reliability in Uncertainty Quantification
Utilizing Adjoint-Based Techniques to Improve the Accuracy and Reliability in Uncertainty Quantification Tim Wildey Sandia National Laboratories Center for Computing Research (CCR) Collaborators: E. Cyr,
More informationMath Tune-Up Louisiana State University August, Lectures on Partial Differential Equations and Hilbert Space
Math Tune-Up Louisiana State University August, 2008 Lectures on Partial Differential Equations and Hilbert Space 1. A linear partial differential equation of physics We begin by considering the simplest
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 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 informationStatistics for Data Analysis. Niklaus Berger. PSI Practical Course Physics Institute, University of Heidelberg
Statistics for Data Analysis PSI Practical Course 2014 Niklaus Berger Physics Institute, University of Heidelberg Overview You are going to perform a data analysis: Compare measured distributions to theoretical
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 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 informationNumerical Analysis and Methods for PDE I
Numerical Analysis and Methods for PDE I A. J. Meir Department of Mathematics and Statistics Auburn University US-Africa Advanced Study Institute on Analysis, Dynamical Systems, and Mathematical Modeling
More 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 informationCourse Description for Real Analysis, Math 156
Course Description for Real Analysis, Math 156 In this class, we will study elliptic PDE, Fourier analysis, and dispersive PDE. Here is a quick summary of the topics will study study. They re described
More informationLecture Introduction
Lecture 1 1.1 Introduction The theory of Partial Differential Equations (PDEs) is central to mathematics, both pure and applied. The main difference between the theory of PDEs and the theory of Ordinary
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 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 informationSequential Monte Carlo Samplers for Applications in High Dimensions
Sequential Monte Carlo Samplers for Applications in High Dimensions Alexandros Beskos National University of Singapore KAUST, 26th February 2014 Joint work with: Dan Crisan, Ajay Jasra, Nik Kantas, Alex
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 informationHairer /Gubinelli-Imkeller-Perkowski
Hairer /Gubinelli-Imkeller-Perkowski Φ 4 3 I The 3D dynamic Φ 4 -model driven by space-time white noise Let us study the following real-valued stochastic PDE on (0, ) T 3, where ξ is the space-time white
More informationDepartment of Electrical- and Information Technology. ETS061 Lecture 3, Verification, Validation and Input
ETS061 Lecture 3, Verification, Validation and Input Verification and Validation Real system Validation Model Verification Measurements Verification Break model down into smaller bits and test each bit
More informationDipartimento di Scienze Matematiche
Exploiting parallel computing in Discrete Fracture Network simulations: an inherently parallel optimization approach Stefano Berrone stefano.berrone@polito.it Team: Matìas Benedetto, Andrea Borio, Claudio
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 informationClass Meeting # 2: The Diffusion (aka Heat) Equation
MATH 8.52 COURSE NOTES - CLASS MEETING # 2 8.52 Introduction to PDEs, Fall 20 Professor: Jared Speck Class Meeting # 2: The Diffusion (aka Heat) Equation The heat equation for a function u(, x (.0.). Introduction
More informationMcGill University Department of Mathematics and Statistics. Ph.D. preliminary examination, PART A. PURE AND APPLIED MATHEMATICS Paper BETA
McGill University Department of Mathematics and Statistics Ph.D. preliminary examination, PART A PURE AND APPLIED MATHEMATICS Paper BETA 17 August, 2018 1:00 p.m. - 5:00 p.m. INSTRUCTIONS: (i) This paper
More informationPDEs, part 1: Introduction and elliptic PDEs
PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,
More informationChapter 9. Non-Parametric Density Function Estimation
9-1 Density Estimation Version 1.2 Chapter 9 Non-Parametric Density Function Estimation 9.1. Introduction We have discussed several estimation techniques: method of moments, maximum likelihood, and least
More information1 Introduction. 2 Diffusion equation and central limit theorem. The content of these notes is also covered by chapter 3 section B of [1].
1 Introduction The content of these notes is also covered by chapter 3 section B of [1]. Diffusion equation and central limit theorem Consider a sequence {ξ i } i=1 i.i.d. ξ i = d ξ with ξ : Ω { Dx, 0,
More informationGeometry and the Dirichlet Problem in Any Co-dimension
Geometry and the Dirichlet Problem in Any Co-dimension Max Engelstein (joint work with G. David (U. Paris Sud) and S. Mayboroda (U. Minn.)) Massachusetts Institute of Technology January 10, 2019 This research
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 informationINTRODUCTION TO PDEs
INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial
More informationPartial Differential Equations
Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives with respect to those variables. Most (but
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 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 informationA Brief Overview of Uncertainty Quantification and Error Estimation in Numerical Simulation
A Brief Overview of Uncertainty Quantification and Error Estimation in Numerical Simulation Tim Barth Exploration Systems Directorate NASA Ames Research Center Moffett Field, California 94035-1000 USA
More informationMacroscopic Failure Analysis Based on a Random Field Representations Generated from Material Microstructures
Macroscopic Failure Analysis Based on a Random Field Representations Generated from Material Microstructures Reza Abedi Mechanical, Aerospace & Biomedical Engineering University of Tennessee Knoxville
More informationLecture 3 Partial Differential Equations
Lecture 3 Partial Differential Equations Prof. Massimo Guidolin Prep Course in Investments August-September 2016 Plan of the lecture Motivation and generalities The heat equation and its applications in
More informationDiscretization of PDEs and Tools for the Parallel Solution of the Resulting Systems
Discretization of PDEs and Tools for the Parallel Solution of the Resulting Systems Stan Tomov Innovative Computing Laboratory Computer Science Department The University of Tennessee Wednesday April 4,
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 informationNonparametric Bayesian Methods (Gaussian Processes)
[70240413 Statistical Machine Learning, Spring, 2015] Nonparametric Bayesian Methods (Gaussian Processes) Jun Zhu dcszj@mail.tsinghua.edu.cn http://bigml.cs.tsinghua.edu.cn/~jun State Key Lab of Intelligent
More informationPartial differential equations
Partial differential equations Many problems in science involve the evolution of quantities not only in time but also in space (this is the most common situation)! We will call partial differential equation
More informationIntroduction to Random Diffusions
Introduction to Random Diffusions The main reason to study random diffusions is that this class of processes combines two key features of modern probability theory. On the one hand they are semi-martingales
More informationPreliminary Examination in Numerical Analysis
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August 7, 06, 0 am pm. Submit solutions to four (and no more) of the following six problems. Show all your work, and justify
More informationSystem Identification
System Identification Arun K. Tangirala Department of Chemical Engineering IIT Madras July 27, 2013 Module 3 Lecture 1 Arun K. Tangirala System Identification July 27, 2013 1 Objectives of this Module
More informationReal Analysis Prof. S.H. Kulkarni Department of Mathematics Indian Institute of Technology, Madras. Lecture - 13 Conditional Convergence
Real Analysis Prof. S.H. Kulkarni Department of Mathematics Indian Institute of Technology, Madras Lecture - 13 Conditional Convergence Now, there are a few things that are remaining in the discussion
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 informationHilbert Space Methods for Reduced-Rank Gaussian Process Regression
Hilbert Space Methods for Reduced-Rank Gaussian Process Regression Arno Solin and Simo Särkkä Aalto University, Finland Workshop on Gaussian Process Approximation Copenhagen, Denmark, May 2015 Solin &
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 informationParameter Selection Techniques and Surrogate Models
Parameter Selection Techniques and Surrogate Models Model Reduction: Will discuss two forms Parameter space reduction Surrogate models to reduce model complexity Input Representation Local Sensitivity
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 informationLearning gradients: prescriptive models
Department of Statistical Science Institute for Genome Sciences & Policy Department of Computer Science Duke University May 11, 2007 Relevant papers Learning Coordinate Covariances via Gradients. Sayan
More informationChapter 9. Non-Parametric Density Function Estimation
9-1 Density Estimation Version 1.1 Chapter 9 Non-Parametric Density Function Estimation 9.1. Introduction We have discussed several estimation techniques: method of moments, maximum likelihood, and least
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 informationRandom Vibrations & Failure Analysis Sayan Gupta Indian Institute of Technology Madras
Random Vibrations & Failure Analysis Sayan Gupta Indian Institute of Technology Madras Lecture 1: Introduction Course Objectives: The focus of this course is on gaining understanding on how to make an
More 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 informationHypocoercivity and Sensitivity Analysis in Kinetic Equations and Uncertainty Quantification October 2 nd 5 th
Hypocoercivity and Sensitivity Analysis in Kinetic Equations and Uncertainty Quantification October 2 nd 5 th Department of Mathematics, University of Wisconsin Madison Venue: van Vleck Hall 911 Monday,
More informationDiffusion and cellular-level simulation. CS/CME/BioE/Biophys/BMI 279 Nov. 7 and 9, 2017 Ron Dror
Diffusion and cellular-level simulation CS/CME/BioE/Biophys/BMI 279 Nov. 7 and 9, 2017 Ron Dror 1 Outline How do molecules move around in a cell? Diffusion as a random walk (particle-based perspective)
More informationc 1 v 1 + c 2 v 2 = 0 c 1 λ 1 v 1 + c 2 λ 1 v 2 = 0
LECTURE LECTURE 2 0. Distinct eigenvalues I haven t gotten around to stating the following important theorem: Theorem: A matrix with n distinct eigenvalues is diagonalizable. Proof (Sketch) Suppose n =
More informationNONLOCAL PROBLEMS WITH LOCAL DIRICHLET AND NEUMANN BOUNDARY CONDITIONS BURAK AKSOYLU AND FATIH CELIKER
NONLOCAL PROBLEMS WITH LOCAL DIRICHLET AND NEUMANN BOUNDARY CONDITIONS BURAK AKSOYLU AND FATIH CELIKER Department of Mathematics, Wayne State University, 656 W. Kirby, Detroit, MI 480, USA. Department
More informationComputational Techniques Prof. Sreenivas Jayanthi. Department of Chemical Engineering Indian institute of Technology, Madras
Computational Techniques Prof. Sreenivas Jayanthi. Department of Chemical Engineering Indian institute of Technology, Madras Module No. # 05 Lecture No. # 24 Gauss-Jordan method L U decomposition method
More informationStochastic Structural Dynamics Prof. Dr. C. S. Manohar Department of Civil Engineering Indian Institute of Science, Bangalore
Stochastic Structural Dynamics Prof. Dr. C. S. Manohar Department of Civil Engineering Indian Institute of Science, Bangalore Lecture No. # 33 Probabilistic methods in earthquake engineering-2 So, we have
More information5.4 Continuity: Preliminary Notions
5.4. CONTINUITY: PRELIMINARY NOTIONS 181 5.4 Continuity: Preliminary Notions 5.4.1 Definitions The American Heritage Dictionary of the English Language defines continuity as an uninterrupted succession,
More informationAn introduction to PDE-constrained optimization
An introduction to PDE-constrained optimization Wolfgang Bangerth Department of Mathematics Texas A&M University 1 Overview Why partial differential equations? Why optimization? Examples of PDE optimization
More informationGeneral Technical Remarks on PDE s and Boundary Conditions Kurt Bryan MA 436
General Technical Remarks on PDE s and Boundary Conditions Kurt Bryan MA 436 1 Introduction You may have noticed that when we analyzed the heat equation on a bar of length 1 and I talked about the equation
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 informationCS 542G: The Poisson Problem, Finite Differences
CS 542G: The Poisson Problem, Finite Differences Robert Bridson November 10, 2008 1 The Poisson Problem At the end last time, we noticed that the gravitational potential has a zero Laplacian except at
More informationDownloaded 08/15/18 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 40, No. 4, pp. A2152 A2173 c 2018 Society for Industrial and Applied Mathematics STOCHASTIC DOMAIN DECOMPOSITION VIA MOMENT MINIMIZATION DONGKUN ZHANG, HESSAM BABAEE, AND GEORGE
More informationUniversity of Regina. Lecture Notes. Michael Kozdron
University of Regina Statistics 252 Mathematical Statistics Lecture Notes Winter 2005 Michael Kozdron kozdron@math.uregina.ca www.math.uregina.ca/ kozdron Contents 1 The Basic Idea of Statistics: Estimating
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 informationRegularization by noise in infinite dimensions
Regularization by noise in infinite dimensions Franco Flandoli, University of Pisa King s College 2017 Franco Flandoli, University of Pisa () Regularization by noise King s College 2017 1 / 33 Plan of
More informationLast Update: April 7, 201 0
M ath E S W inter Last Update: April 7, Introduction to Partial Differential Equations Disclaimer: his lecture note tries to provide an alternative approach to the material in Sections.. 5 in the textbook.
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 informationPhysics 6303 Lecture 11 September 24, LAST TIME: Cylindrical coordinates, spherical coordinates, and Legendre s equation
Physics 6303 Lecture September 24, 208 LAST TIME: Cylindrical coordinates, spherical coordinates, and Legendre s equation, l l l l l l. Consider problems that are no axisymmetric; i.e., the potential depends
More information