Scalable Algorithms for Optimal Control of Systems Governed by PDEs Under Uncertainty
|
|
- Ronald Perkins
- 5 years ago
- Views:
Transcription
1 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 University 2 Institute for Computational Engineering and Sciences The University of Texas at Austin 3 School of Natural Sciences University of California, Merced 3 Courant Institute of Mathematical Sciences New York University Advances in Uncertainty Quantification Methods, Algorithms and Applications Annual Meeting of KAUST SRI-UQ Center King Abdullah University of Science and Technology, Saudi Arabia January 5 10, 2016
2 From data to decisions under uncertainty Uncertain parameter m Mathematical Model A(m, u) = f Quantity of Interest (QoI) q(m) Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
3 From data to decisions under uncertainty Uncertain parameter m prior posterior Bayesian inversion π post(m y) π like(y m)π 0(m) Experimental data y πlike(y m) = πnoise(bu y) Mathematical Model A(m, u) = f Quantity of Interest (QoI) q(m) Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
4 From data to decisions under uncertainty Uncertain parameter m prior posterior Bayesian inversion π post(m y) π like(y m)π 0(m) Experimental data y πlike(y m) = πnoise(bu y) Mathematical Model A(m, u) = f experiments Quantity of Interest (QoI) q(m) Design of experiments min Ψ [ π post(m y; ξ) ] π(y)dy ξ ξ : experimental design Ψ: design objective/criterion Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
5 From data to decisions under uncertainty Uncertain parameter m prior posterior Bayesian inversion π post(m y) π like(y m)π 0(m) Experimental data y πlike(y m) = πnoise(bu y) Mathematical Model A(m, u) = f( ; z) z := control function experiments Quantity of Interest (QoI) q(z, m) Design of experiments min Ψ [ π post(m y; ξ) ] π(y)dy ξ ξ : experimental design Ψ: design objective/criterion Optimal control/design under uncertainty e.g. min q(z, m)µ(dm) z Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
6 Example: Groundwater contaminant remediation Source: Reed Maxwell, CSM Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
7 Example: Groundwater contaminant remediation Inverse problem Infer (uncertain) soil permeability from (uncertain) measurements of pressure head at wells and from a (uncertain) model of subsurface flow and transport Prediction (or forward) problem Predict (uncertain) evolution of contaminant concentration at municipal wells from (uncertain) permeability and (uncertain) subsurface flow/transport model Optimal experimental design problem Where should new observation wells be placed so that permeability is inferred with the least uncertainty? Optimal design problem Where should new remediation wells be placed so that (uncertain) contaminant concentrations at municipal wells are minimized? Optimal control problem What should the rates of extraction/injection at remediation wells be so that (uncertain) contaminant concentrations at municipal wells are minimized? Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
8 Optimal control of systems governed by PDEs with uncertain parameter fields PDE-constrained control objective Control of injection wells in porous medium flow (SPE10 permeability data) q = q(u(z, m)) where A(u, m) = f(z) q: control objective A: forward operator m: uncertain parameter field z: control function Problem: given the uncertainty model for m, find z that optimizes q(u(z, m)) Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
9 Optimization under uncertainty (OUU) H : parameter space, infinite-dimensional separable Hilbert space q(z, m): control objective functional m H : uncertain model parameter field, Optimization under uncertainty (OUU): min z q(z, m) z: control function Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
10 Optimization under uncertainty (OUU) H : parameter space, infinite-dimensional separable Hilbert space q(z, m): control objective functional m H : uncertain model parameter field, z: control function Risk-neutral optimization under uncertainty (OUU): min E m {q(z, m)} z E m {q(z, m)} = H q(z, m) µ(dm) Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
11 Optimization under uncertainty (OUU) H : parameter space, infinite-dimensional separable Hilbert space q(z, m): control objective functional m H : uncertain model parameter field, Risk-averse optimization under uncertainty (OUU): min E m {q(z, m)} + β var m {q(z, m)} z E m {q(z, m)} = q(z, m) µ(dm) H z: control function var m {q(z, m)} = E m {q(z, m) 2 } E m {q(z, m)} 2 Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
12 Optimization under uncertainty (OUU) H : parameter space, infinite-dimensional separable Hilbert space q(z, m): control objective functional m H : uncertain model parameter field, Risk-averse optimization under uncertainty (OUU): min E m {q(z, m)} + β var m {q(z, m)} z E m {q(z, m)} = q(z, m) µ(dm) H z: control function var m {q(z, m)} = E m {q(z, m) 2 } E m {q(z, m)} 2 Main challenges: Integration over infinite/high-dimensional parameter space Evaluation of q requires PDE solves Standard Monte Carlo approach (Sample Average Approximation) is prohibitive Numerous (n mc) samples required, each requires PDE solve Resulting PDE-constrained optimization problem has n mc PDE constraints Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
13 Some existing approaches for PDE-constrained OUU Methods based on stochastic collocation, sparse/adaptive sampling, POD,... Schulz & Schillings, Problem formulations and treatment of uncertainties in aerodynamic design, AIAA J, Borzì & von Winckel, Multigrid methods and sparse-grid collocation techniques for parabolic optimal control problems with random coefficients, SISC, Borzì, Schillings, & von Winckel, On the treatment of distributed uncertainties in PDE-constrained optimization, GAMM-Mitt Borzì & von Winckel, A POD framework to determine robust controls in PDE optimization, Computing and Visualization in Science, Gunzburger & Ming, Optimal control of stochastic flow over a backward-facing step using reduced-order modeling, SISC Gunzburger, Lee, & Lee, Error estimates of stochastic optimal Neumann boundary control problems, SINUM, Kunoth & Schwab, Analytic Regularity and GPC Approximation for Control Problems Constrained by Linear Parametric Elliptic and Parabolic PDEs, SICON, Tiesler, Kirby, Xiu, & Preusser, Stochastic collocation for optimal control problems with stochastic PDE constraints, SICON, Kouri, Heinkenschloss, Ridzal, & Van Bloemen Waanders, A trust-region algorithm with adaptive stochastic collocation for PDE optimization under uncertainty, SISC, Chen, Quarteroni, & Rozza, Stochastic optimal Robin boundary control problems of advection-dominated elliptic equations, SINUM, Kouri, A multilevel stochastic collocation algorithm for optimization of PDEs with uncertain coefficients, JUQ, Chen & Quarteroni, Weighted reduced basis method for stochastic optimal control problems with elliptic PDE constraint, JUQ, Chen, Quarteroni, & Rozza, Multilevel and weighted reduced basis method for stochastic optimal control problems constrained by Stokes equations, Num. Math Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
14 Control of injection wells in a porous medium flow m = mean of log permeability field q = target pressure at production wells state PDE: single phase flow in a porous medium n c (e m u) = z i f i (x) with Dirchlet lateral & Neumann top/bottom BCs uncertain parameter: log permeability field m: control variables: z i, mass source at injection wells; f i, mollified Dirac deltas control objective: q(z, m) := 1 2 Qu(z, m) q 2, q: target pressure dimensions: n s = n m = 3242, n c = 20, n q = 12 i=1 Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
15 Porous medium with random permeability field Distribution law of m: µ = N ( m, C) (Gaussian measure on Hilbert space H ) Take covariance operator as square of inverse of Poisson-like operator: C = ( κ + αi) 2 κ, α > 0 C is positive, self-adjoint, of trace-class; µ well-defined on H (Stuart 10) κ α correlation length; the larger α, the smaller the variance Random draws for κ = , α = 4 Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
16 OUU with linearized parameter-to-objective map Risk-averse optimal control problem (including cost of controls) min z E m {q(z, m)} + β var m {q(z, m)} + γ z 2 Linear approximation to parameter-to-objective map g m (z, ) := dq(z, ) dm q lin (z, m) = q(z, m) + g m (z, m), m m is the gradient with respect to m The moments of the linearized objective: E m {q lin (z, )} = q(z, m), var m {q lin (z, )} = g m (z, m), C[g m (z, m)] ( ) q lin (z, ) N q(z, m), g m (z, m), C[g m (z, m)] Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
17 Risk-averse optimal control problem with linearized parameter-to-objective map State-and-adjoint-PDE constrained optimization problem (quartic in z): min z Z J (z) := 1 2 Qu q 2 + β 2 gm( m), C[gm( m)] + γ 2 z 2 with g m( m) =e m u p, where n c (e m u) = z if i i=1 (e m p) = Q (Qu q) state equation adjoint equation Lagrangian of the risk-averse optimal control problem with q lin : L (z, u, p, u, p ) = 1 2 Qu q 2 + β e m u p, C[e m u p] + γ 2 2 z 2 + e m u, u n c z i f i, u i=1 + e m p, p + Q (Qu q), p Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
18 Gradient computation for risk averse optimal control State problem for risk-averse optimal control problem with q lin : e m u, ũ n c = z i f i, ũ i=1 e m p, p = Q (Qu q), p for all test functions p and ũ Adjoint problem for risk-averse optimal control problem with q lin : e m p, p = β e m u p, C[e m u p] e m u, ũ = Q (Qu q), ũ β e m ũ p, C[e m u p] Q Qp, ũ for all test functions p and ũ J Gradient: = γz j f j, u, j = 1,..., n c z j Cost of objective = 2 PDE solves; cost of gradient = 2 PDE solves Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
19 Risk-averse optimal control with linearized objective distribution Θ(z 0, ) Θlin(z 0, ) Θquad(z 0, ) initial (suboptimal) control z 0 distrib. of exact & approx objectives at z 0 Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
20 Risk-averse optimal control with linearized objective distribution Θ(z 0, ) Θlin(z 0, ) Θquad(z 0, ) initial (suboptimal) control z 0 distrib. of exact & approx objectives at z Θ(z opt lin, ) Θlin(z opt lin, ) Θquad(z opt lin, ) distribution optimal control z opt lin based on q lin distrib. of exact & approx objectives at z opt lin Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
21 Quadratic approximation to parameter-to-objective map Quadratic approximation to the parameter-to-control-objective map: q quad (z, m) = q(z, m) + g m (z, m), m m H m(z, m)(m m), m m g m : gradient of parameter-to-objective map H m : Hessian of parameter-to-objective map Observations: Quadratic approximation does not lead to a Gaussian control objective However, can derive analytic formulas for the moments of q quad in the infinite-dimensional Hilbert space setting Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
22 Analytic expressions for mean and variance with q quad Mean: E m {q quad (z, )} = q(z, m)+ 1 2 tr[ Hm (z, m) ] Variance: var m {q quad (z, )} = g m (z, m), C[g m (z, m)] tr[ Hm (z, m) 2] where H m = C 1/2 H m C 1/2 is the covariance-preconditioned Hessian Risk averse optimal control objective with q quad : J(z) = q(z, m) tr[ Hm (z, m) ] + β { g m (z, m), C[g m (z, m)] tr[ Hm (z, m) 2]} Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
23 Randomized trace estimator Randomized trace estimation: tr( H m ) 1 n tr tr( H 2 m) 1 n tr n tr n tr Hm ξ j, ξ j = 1 n tr H m ζ j, C[H m ζ j ] n tr H m ζ j, ζ j where ζ j = C 1/2 ξ j In computations, we use draws ζ j N (0, C) =: ν Straightforward to show: H m ζ, ζ ν(dζ) = tr( H m ), H H H m ζ, C[H m ζ] ν(dζ) = tr( H 2 m) Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
24 Risk-averse optimal control with quadraticized objective 1 min z Z 2 Qu q n tr 2n tr ζ j, η j + β 2 { g m ( m), C[g m ( m)] + 1 n tr 2n tr } C 1/2 η j 2 with g m ( m) =e m u p η j = e m (ζ j u p + υ j p + u ρ j ) }{{} H mζ j j {1,..., n tr} where (e m u) = N i=1 z if i (e m p) = Q (Qu q) (e m υ j ) = (ζ j e m u) (e m ρ j ) = Q Qυ j + (ζ j e m p) Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
25 Risk-averse optimal control with quadraticized objective 1 min z Z 2 Qu q n tr 2n tr ζ j, η j + β 2 { g m ( m), C[g m ( m)] + 1 n tr 2n tr } C 1/2 η j 2 with g m ( m) =e m u p η j = e m (ζ j u p + υ j p + u ρ j ) }{{} H mζ j j {1,..., n tr} where (e m u) = N i=1 z if i (e m p) = Q (Qu q) (e m υ j ) = (ζ j e m u) (e m ρ j ) = Q Qυ j + (ζ j e m p) Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
26 Risk-averse optimal control with quadraticized objective 1 min z Z 2 Qu q n tr 2n tr ζ j, η j + β 2 { g m ( m), C[g m ( m)] + 1 n tr 2n tr } C 1/2 η j 2 with g m ( m) =e m u p η j = e m (ζ j u p + υ j p + u ρ j ) }{{} H mζ j j {1,..., n tr} where (e m u) = N i=1 z if i (e m p) = Q (Qu q) (e m υ j ) = (ζ j e m u) (e m ρ j ) = Q Qυ j + (ζ j e m p) Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
27 Risk-averse optimal control with quadraticized objective 1 min z Z 2 Qu q n tr 2n tr ζ j, η j + β 2 { g m ( m), C[g m ( m)] + 1 n tr 2n tr } C 1/2 η j 2 with g m ( m) =e m u p η j = e m (ζ j u p + υ j p + u ρ j ) }{{} H mζ j j {1,..., n tr} where (e m u) = N i=1 z if i (e m p) = Q (Qu q) (e m υ j ) = (ζ j e m u) (e m ρ j ) = Q Qυ j + (ζ j e m p) Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
28 Risk-averse optimal control with quadraticized objective 1 min z Z 2 Qu q n tr 2n tr ζ j, η j + β 2 { g m ( m), C[g m ( m)] + 1 n tr 2n tr } C 1/2 η j 2 with g m ( m) =e m u p η j = e m (ζ j u p + υ j p + u ρ j ) }{{} H mζ j j {1,..., n tr} where (e m u) = N i=1 z if i (e m p) = Q (Qu q) (e m υ j ) = (ζ j e m u) (e m ρ j ) = Q Qυ j + (ζ j e m p) Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
29 Lagrangian for risk-averse optimal control with q quad L (z, u, p,{υ j} ntr, {ρ j} ntr, u, p, {υ j } ntr, {ρ j } ntr ) = 1 2 Qu q n tr n tr ζj, [ e m (ζ j u p + υ j p + u ρ j) ] + β e m u p, C[e m u p] 2 + β n tr C 1/2 [ e m (ζ j u p + υ j p + u ρ ] 2 j) 4n tr + e m u, u N i=1 zi fi, u + e m p, p + Q (Qu q), p n tr [ e + m υ j, υ j + ζje m u, υ ] j n tr [ e + m ρ j, ρ j + Q Qυ j, ρ j + ζje m p, ρ ] j Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
30 Adjoint & gradient for risk-averse optimal control w/q quad Adjoint problem for q quad approximation (e m ρ j ) = b (j) 1 j {1,..., n tr } (e m υ j ) + Q Qρ j = b (j) 2 j {1,..., n tr } n tr (e m p ) (ζ j e m ρ j ) = b 3 n tr (e m u ) + Q Qp (ζ j e m υ j ) = b 4 Gradient for q quad approximation L z j = γz j f j, u, j = 1,..., n c Cost of objective = 2 + 2n tr PDE solves; cost of gradient = 2 + 2n tr PDE solves Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
31 Risk-averse optimal control with quadraticized objective distribution Θ(z 0, ) Θlin(z 0, ) Θquad(z 0, ) initial (suboptimal) control z 0 distrib. of exact & approx objectives at z 0 Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
32 Risk-averse optimal control with quadraticized objective distribution Θ(z 0, ) Θlin(z 0, ) Θquad(z 0, ) initial (suboptimal) control z 0 distrib. of exact & approx objectives at z Θ(z opt quad, ) Θquad(z opt quad, ) distribution optimal control z opt quad based on q quad distrib. of exact & approx objectives at z opt quad Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
33 Effect of number of trace estimator vectors on distribution of control objective evaluated at optimal controls distribution n tr = 5 n tr = 20 n tr = 40 n tr = value of control objective Optimal controls z opt quad computed for each value of trace estimator using quadratic approximation of control objective, q quad Each curve based on 10,000 samples of distribution of q(z opt quad, m) (control objective evaluated at optimal control) Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
34 Comparison of distribution of control objective for optimal controls based on linearized and quadraticized objective distribution sample mean sample size sample variance sample size Θ(z opt quad, ) Θ(z opt lin, ) value of control objective Comparison of the distributions of q(z opt lin, m) with q(zopt quad, m) β = 1, γ = 10 5 and n tr = 40 trace estimation vectors KDE results are based on 10,000 samples Inserts show Monte Carlo sample convergence for mean and variance Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
35 Summary Some concluding remarks: Optimal control of PDEs with infinite-dimensional parameters Use of local approximations to parameter-to-objective map Quadratic approximation is effective in capturing the distribution of the control objective for the target problem Computational complexity of objective & gradient evaluation is independent of problem dimension (n m = 3242) objective + gradient cost = 125 PDE solves objective/gradient cost for Monte Carlo/SAA: 10,000 PDE solves differences even more striking for nonlinear PDE state problems Possible extensions: Higher moments (e.g., for forward UQ) Alternative risk measures Higher order approximations, e.g., third order Taylor expansion with proper tensor contraction of third order derivative tensor A. Alexanderian, N. Petra, G. Stadler, and O. Ghattas, Mean-variance risk-averse optimal control of systems governed by PDEs with random coefficient fields using second order approximations, in submission. Omar Ghattas (ICES, UT Austin) Optimal control under uncertainty UQAW / 23
Scalable algorithms for optimal experimental design for infinite-dimensional nonlinear Bayesian inverse problems
Scalable algorithms for optimal experimental design for infinite-dimensional nonlinear Bayesian inverse problems Alen Alexanderian (Math/NC State), Omar Ghattas (ICES/UT-Austin), Noémi Petra (Applied Math/UC
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 informationMathematisches Forschungsinstitut Oberwolfach. Numerical Methods for PDE Constrained Optimization with Uncertain Data
Mathematisches Forschungsinstitut Oberwolfach Report No. 04/2013 DOI: 10.4171/OWR/2013/04 Numerical Methods for PDE Constrained Optimization with Uncertain Data Organised by Matthias Heinkenschloss, Houston
More informationFast algorithms for the inverse medium problem. George Biros University of Pennsylvania
Fast algorithms for the inverse medium problem George Biros University of Pennsylvania Acknowledgments S. Adavani, H. Sundar, S. Rahul (grad students) C. Davatzikos, D. Shen, H. Litt (heart project) Akcelic,
More informationSparse Quadrature Algorithms for Bayesian Inverse Problems
Sparse Quadrature Algorithms for Bayesian Inverse Problems Claudia Schillings, Christoph Schwab Pro*Doc Retreat Disentis 2013 Numerical Analysis and Scientific Computing Disentis - 15 August, 2013 research
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 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 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 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 informationMulti-Objective Stochastic Optimization of Magnetic Fields. David S Bindel, Cornell University
Multi-Objective Stochastic Optimization of Magnetic Fields David S Bindel, Cornell University Current State of the Art (STELLOPT) Optimizer Calculate! " (physics + engineering targets) Adjust plasma boundary
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 informationTrust-Region SQP Methods with Inexact Linear System Solves for Large-Scale Optimization
Trust-Region SQP Methods with Inexact Linear System Solves for Large-Scale Optimization Denis Ridzal Department of Computational and Applied Mathematics Rice University, Houston, Texas dridzal@caam.rice.edu
More informationAt-A-Glance SIAM Events Mobile App
At-A-Glance SIAM Events 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.you can also visit www.tripbuildermedia.com/apps/siamevents
More informationDimension-Independent likelihood-informed (DILI) MCMC
Dimension-Independent likelihood-informed (DILI) MCMC Tiangang Cui, Kody Law 2, Youssef Marzouk Massachusetts Institute of Technology 2 Oak Ridge National Laboratory 2 August 25 TC, KL, YM DILI MCMC USC
More informationSampling and low-rank tensor approximation of the response surface
Sampling and low-rank tensor approximation of the response surface tifica Alexander Litvinenko 1,2 (joint work with Hermann G. Matthies 3 ) 1 Group of Raul Tempone, SRI UQ, and 2 Group of David Keyes,
More informationTowards Large-Scale Computational Science and Engineering with Quantifiable Uncertainty
Towards Large-Scale Computational Science and Engineering with Quantifiable Uncertainty Tan Bui-Thanh Computational Engineering and Optimization (CEO) Group Department of Aerospace Engineering and Engineering
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 informationLecture 3: Adaptive Construction of Response Surface Approximations for Bayesian Inference
Lecture 3: Adaptive Construction of Response Surface Approximations for Bayesian Inference Serge Prudhomme Département de mathématiques et de génie industriel Ecole Polytechnique de Montréal SRI Center
More informationPractical unbiased Monte Carlo for Uncertainty Quantification
Practical unbiased Monte Carlo for Uncertainty Quantification Sergios Agapiou Department of Statistics, University of Warwick MiR@W day: Uncertainty in Complex Computer Models, 2nd February 2015, University
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 informationMeasure-Theoretic parameter estimation and prediction for contaminant transport and coastal ocean modeling
Measure-Theoretic parameter estimation and prediction for contaminant transport and coastal ocean modeling Steven Mattis and Lindley Graham RMSWUQ 2015 July 17, 2015 Collaborators 2/45 The University of
More informationNumerical Solution I
Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions
More 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 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 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 informationIntegration of Sequential Quadratic Programming and Domain Decomposition Methods for Nonlinear Optimal Control Problems
Integration of Sequential Quadratic Programming and Domain Decomposition Methods for Nonlinear Optimal Control Problems Matthias Heinkenschloss 1 and Denis Ridzal 2 1 Department of Computational and Applied
More informationOn a Data Assimilation Method coupling Kalman Filtering, MCRE Concept and PGD Model Reduction for Real-Time Updating of Structural Mechanics Model
On a Data Assimilation Method coupling, MCRE Concept and PGD Model Reduction for Real-Time Updating of Structural Mechanics Model 2016 SIAM Conference on Uncertainty Quantification Basile Marchand 1, Ludovic
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 informationUT-LANL DOE/ASCR ISICLES Project PARADICE: Uncertainty Quantification for Large-Scale Ice Sheet Modeling and Simulation.
UT-LANL DOE/ASCR ISICLES Project PARADICE: Uncertainty Quantification for Large-Scale Ice Sheet Modeling and Simulation Omar Ghattas Representing the rest of the UT-Austin/LANL team: Don Blankenship Tan
More informationComparison of Clenshaw-Curtis and Leja quasi-optimal sparse grids for the approximation of random PDEs
MATHICSE Mathematics Institute of Computational Science and Engineering School of Basic Sciences - Section of Mathematics MATHICSE Technical Report Nr. 41.2014 September 2014 Comparison of Clenshaw-Curtis
More informationSuboptimal Open-loop Control Using POD. Stefan Volkwein
Institute for Mathematics and Scientific Computing University of Graz, Austria PhD program in Mathematics for Technology Catania, May 22, 2007 Motivation Optimal control of evolution problems: min J(y,
More informationRisk Averse Shape Optimization
Risk Averse Shape Optimization Sergio Conti 2 Martin Pach 1 Martin Rumpf 2 Rüdiger Schultz 1 1 Department of Mathematics University Duisburg-Essen 2 Rheinische Friedrich-Wilhelms-Universität Bonn Workshop
More informationPredictive Engineering and Computational Sciences. Local Sensitivity Derivative Enhanced Monte Carlo Methods. Roy H. Stogner, Vikram Garg
PECOS Predictive Engineering and Computational Sciences Local Sensitivity Derivative Enhanced Monte Carlo Methods Roy H. Stogner, Vikram Garg Institute for Computational Engineering and Sciences The University
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 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 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 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 informationUNCERTAINTY QUANTIFICATION IN LIQUID COMPOSITE MOULDING PROCESSES
UNCERTAINTY QUANTIFICATION IN LIQUID COMPOSITE MOULDING PROCESSES Bart Verleye 1, Dirk Nuyens 1, Andrew Walbran 3,2, and Ming Gan 2 1 Dept. Computer Science, KU Leuven, Celestijnenlaan 200A, B-3001 Leuven,
More informationJ~ J.~ Dennis D. Cox Professor of Statistics
An Approach for the Adaptive Solution of Optimization MAY, 2012 HOUSTON, TEXAS Beatrice M. Riviere Associate Professor of Computational and Applied Mathematics J~ J.~ Dennis D. Cox Professor of Statistics
More informationUncertainty quantification for inverse problems with a weak wave-equation constraint
Uncertainty quantification for inverse problems with a weak wave-equation constraint Zhilong Fang*, Curt Da Silva*, Rachel Kuske** and Felix J. Herrmann* *Seismic Laboratory for Imaging and Modeling (SLIM),
More informationQuantifying Stochastic Model Errors via Robust Optimization
Quantifying Stochastic Model Errors via Robust Optimization IPAM Workshop on Uncertainty Quantification for Multiscale Stochastic Systems and Applications Jan 19, 2016 Henry Lam Industrial & Operations
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 informationOptimization Tools in an Uncertain Environment
Optimization Tools in an Uncertain Environment Michael C. Ferris University of Wisconsin, Madison Uncertainty Workshop, Chicago: July 21, 2008 Michael Ferris (University of Wisconsin) Stochastic optimization
More informationProbing the covariance matrix
Probing the covariance matrix Kenneth M. Hanson Los Alamos National Laboratory (ret.) BIE Users Group Meeting, September 24, 2013 This presentation available at http://kmh-lanl.hansonhub.com/ LA-UR-06-5241
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 informationBayesian Inverse problem, Data assimilation and Localization
Bayesian Inverse problem, Data assimilation and Localization Xin T Tong National University of Singapore ICIP, Singapore 2018 X.Tong Localization 1 / 37 Content What is Bayesian inverse problem? What is
More informationParameter estimation by implicit sampling
Parameter estimation by implicit sampling Matthias Morzfeld 1,2,, Xuemin Tu 3, Jon Wilkening 1,2 and Alexandre J. Chorin 1,2 1 Department of Mathematics, University of California, Berkeley, CA 94720, USA.
More informationLog Gaussian Cox Processes. Chi Group Meeting February 23, 2016
Log Gaussian Cox Processes Chi Group Meeting February 23, 2016 Outline Typical motivating application Introduction to LGCP model Brief overview of inference Applications in my work just getting started
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 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 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 informationM E M O R A N D U M. Faculty Senate approved November 1, 2018
M E M O R A N D U M Faculty Senate approved November 1, 2018 TO: FROM: Deans and Chairs Becky Bitter, Sr. Assistant Registrar DATE: October 23, 2018 SUBJECT: Minor Change Bulletin No. 5 The courses listed
More informationVerification and Validation in Computational Engineering and Sciences
Verification and Validation in Computational Engineering and Sciences Serge Prudhomme Département de mathématiques et de génie industriel Ecole Polytechnique de Montréal, Canada SRI Center for Uncertainty
More informationReduced-Order Multiobjective Optimal Control of Semilinear Parabolic Problems. Laura Iapichino Stefan Trenz Stefan Volkwein
Universität Konstanz Reduced-Order Multiobjective Optimal Control of Semilinear Parabolic Problems Laura Iapichino Stefan Trenz Stefan Volkwein Konstanzer Schriften in Mathematik Nr. 347, Dezember 2015
More informationNumerical Analysis of Elliptic PDEs with Random Coefficients
Numerical Analysis of Elliptic PDEs with Random Coefficients (Lecture I) Robert Scheichl Department of Mathematical Sciences University of Bath Workshop on PDEs with Random Coefficients Weierstrass Institute,
More informationGroundwater Resources Management under. Uncertainty Harald Kunstmann*, Wolfgang Kinzelbach* & Gerrit van Tender**
Groundwater Resources Management under Uncertainty Harald Kunstmann*, Wolfgang Kinzelbach* & Gerrit van Tender**, C/f &OP3 Zwnc/?, Email: kunstmann@ihw. baum. ethz. ch Abstract Groundwater resources and
More informationPartial Differential Equations
Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This
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 informationBayesian parameter estimation in predictive engineering
Bayesian parameter estimation in predictive engineering Damon McDougall Institute for Computational Engineering and Sciences, UT Austin 14th August 2014 1/27 Motivation Understand physical phenomena Observations
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 informationGaussian Process Approximations of Stochastic Differential Equations
Gaussian Process Approximations of Stochastic Differential Equations Cédric Archambeau Centre for Computational Statistics and Machine Learning University College London c.archambeau@cs.ucl.ac.uk CSML
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 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 informationA Recursive Trust-Region Method for Non-Convex Constrained Minimization
A Recursive Trust-Region Method for Non-Convex Constrained Minimization Christian Groß 1 and Rolf Krause 1 Institute for Numerical Simulation, University of Bonn. {gross,krause}@ins.uni-bonn.de 1 Introduction
More informationIntroduction to Aspects of Multiscale Modeling as Applied to Porous Media
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part IV Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and
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 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 informationFinite Difference Methods (FDMs) 1
Finite Difference Methods (FDMs) 1 1 st - order Approxima9on Recall Taylor series expansion: Forward difference: Backward difference: Central difference: 2 nd - order Approxima9on Forward difference: Backward
More informationGlobal, Robust, Multi-Objective Optimization of Stellarators
Global, Robust, Multi-Objective Optimization of Stellarators David Bindel 29 March 2019 Department of Computer Science Cornell University 1 Quality Criteria What Makes a Good Stellarator? Half-Module (
More informationDownloaded 05/15/18 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 39, No. 5, pp. A2365 A2393 A DATA SCALABLE AUGMENTED LAGRANGIAN KKT PRECONDITIONER FOR LARGE-SCALE INVERSE PROBLEMS NICK ALGER, UMBERTO VILLA, TAN BUI-THANH, AND OMAR GHATTAS
More informationUncertainty quantification for Wavefield Reconstruction Inversion
Uncertainty quantification for Wavefield Reconstruction Inversion Zhilong Fang *, Chia Ying Lee, Curt Da Silva *, Felix J. Herrmann *, and Rachel Kuske * Seismic Laboratory for Imaging and Modeling (SLIM),
More informationResolving the White Noise Paradox in the Regularisation of Inverse Problems
1 / 32 Resolving the White Noise Paradox in the Regularisation of Inverse Problems Hanne Kekkonen joint work with Matti Lassas and Samuli Siltanen Department of Mathematics and Statistics University of
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 informationMethods of Data Assimilation and Comparisons for Lagrangian Data
Methods of Data Assimilation and Comparisons for Lagrangian Data Chris Jones, Warwick and UNC-CH Kayo Ide, UCLA Andrew Stuart, Jochen Voss, Warwick Guillaume Vernieres, UNC-CH Amarjit Budiraja, UNC-CH
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 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 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 informationOutline Lecture 2 2(32)
Outline Lecture (3), Lecture Linear Regression and Classification it is our firm belief that an understanding of linear models is essential for understanding nonlinear ones Thomas Schön Division of Automatic
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 information"ZERO-POINT" IN THE EVALUATION OF MARTENS HARDNESS UNCERTAINTY
"ZERO-POINT" IN THE EVALUATION OF MARTENS HARDNESS UNCERTAINTY Professor Giulio Barbato, PhD Student Gabriele Brondino, Researcher Maurizio Galetto, Professor Grazia Vicario Politecnico di Torino Abstract
More informationEnsemble Kalman filter assimilation of transient groundwater flow data via stochastic moment equations
Ensemble Kalman filter assimilation of transient groundwater flow data via stochastic moment equations Alberto Guadagnini (1,), Marco Panzeri (1), Monica Riva (1,), Shlomo P. Neuman () (1) Department of
More informationICES REPORT March Tan Bui-Thanh And Mark Andrew Girolami
ICES REPORT 4- March 4 Solving Large-Scale Pde-Constrained Bayesian Inverse Problems With Riemann Manifold Hamiltonian Monte Carlo by Tan Bui-Thanh And Mark Andrew Girolami The Institute for Computational
More informationSolution of Stochastic Nonlinear PDEs Using Wiener-Hermite Expansion of High Orders
Solution of Stochastic Nonlinear PDEs Using Wiener-Hermite Expansion of High Orders Dr. Mohamed El-Beltagy 1,2 Joint Wor with Late Prof. Magdy El-Tawil 2 1 Effat University, Engineering College, Electrical
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationMultilevel Monte Carlo for two phase flow and transport in random heterogeneous porous media
Eidgenössische Technische Hochschule Zürich Ecole polytechnique fédérale de Zurich Politecnico federale di Zurigo Swiss Federal Institute of Technology Zurich Multilevel Monte Carlo for two phase flow
More informationMULTISCALE FINITE ELEMENT METHODS FOR STOCHASTIC POROUS MEDIA FLOW EQUATIONS AND APPLICATION TO UNCERTAINTY QUANTIFICATION
MULTISCALE FINITE ELEMENT METHODS FOR STOCHASTIC POROUS MEDIA FLOW EQUATIONS AND APPLICATION TO UNCERTAINTY QUANTIFICATION P. DOSTERT, Y. EFENDIEV, AND T.Y. HOU Abstract. In this paper, we study multiscale
More 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 informationEfficient Variational Inference in Large-Scale Bayesian Compressed Sensing
Efficient Variational Inference in Large-Scale Bayesian Compressed Sensing George Papandreou and Alan Yuille Department of Statistics University of California, Los Angeles ICCV Workshop on Information
More informationCollocation methods for uncertainty quantification in PDE models with random data
Collocation methods for uncertainty quantification in PDE models with random data Fabio Nobile CSQI-MATHICSE, EPFL, Switzerland Acknowlegments: L. Tamellini, G. Migliorati (EPFL), R. Tempone, (KAUST),
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 informationOrganization. I MCMC discussion. I project talks. I Lecture.
Organization I MCMC discussion I project talks. I Lecture. Content I Uncertainty Propagation Overview I Forward-Backward with an Ensemble I Model Reduction (Intro) Uncertainty Propagation in Causal Systems
More informationNo Information Sharing in Oligopoly: The Case of Price Competition with Cost Uncertainty
No Information Sharing in Oligopoly: The Case of Price Competition with Cost Uncertainty Stephan O. Hornig and Manfred Stadler* Abstract We show that concealing cost information is a dominant strategy
More informationautomating the analysis and design of large-scale optimization algorithms
automating the analysis and design of large-scale optimization algorithms Laurent Lessard University of Wisconsin Madison Joint work with Ben Recht, Andy Packard, Bin Hu, Bryan Van Scoy, Saman Cyrus October,
More informationStochastic Analogues to Deterministic Optimizers
Stochastic Analogues to Deterministic Optimizers ISMP 2018 Bordeaux, France Vivak Patel Presented by: Mihai Anitescu July 6, 2018 1 Apology I apologize for not being here to give this talk myself. I injured
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 information. α β γ δ ε ζ η θ ι κ λ μ Aμ ν(x) ξ ο π ρ ς σ τ υ φ χ ψ ω. Α Β Γ Δ Ε Ζ Η Θ Ι Κ Λ Μ Ν Ξ Ο Π Ρ Σ Τ Υ Φ Χ Ψ Ω. Friday April 1 ± ǁ
. α β γ δ ε ζ η θ ι κ λ μ Aμ ν(x) ξ ο π ρ ς σ τ υ φ χ ψ ω. Α Β Γ Δ Ε Ζ Η Θ Ι Κ Λ Μ Ν Ξ Ο Π Ρ Σ Τ Υ Φ Χ Ψ Ω Friday April 1 ± ǁ 1 Chapter 5. Photons: Covariant Theory 5.1. The classical fields 5.2. Covariant
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 informationAdaptive surrogate modeling for response surface approximations with application to bayesian inference
Prudhomme and Bryant Adv. Model. and Simul. in Eng. Sci. 2015 2:22 DOI 10.1186/s40323-015-0045-5 RESEARCH ARTICLE Open Access Adaptive surrogate modeling for response surface approximations with application
More informationAdaptive and multilevel methods for parameter identification in partial differential equations
Adaptive and multilevel methods for parameter identification in partial differential equations Barbara Kaltenbacher, University of Stuttgart joint work with Hend Ben Ameur, Université de Tunis Anke Griesbaum,
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 information