Multigrid and stochastic sparse-grids techniques for PDE control problems with random coefficients
|
|
- Martina Stone
- 5 years ago
- Views:
Transcription
1 Multigrid and stochastic sparse-grids techniques for PDE control problems with random coefficients Università degli Studi del Sannio Dipartimento e Facoltà di Ingegneria, Benevento, Italia
2 Random fields A random field is a generalization of a stochastic process such that the underlying parameter can be a multidimensional vector space or a manifold. Values in a random field are spatially correlated. c(y, u, δ, σ) = 0, with δ(ω), σ(ω) where ω O.
3 Random parameter fields The triple (O, A, P) denotes a probability space where O is the space of elementary events, A is the sigma-algebra of subsets of O, and P the probability measure on O. We assume that δ = δ(t, ω), σ = σ(x, ω) > 0 where t (0, T ), x Ω, and ω O. Examples of finite-dimensional random fields and δ(t, ω) = δ 0 + W (ω) sin(2πt/t ) σ(x 1, x 2, ω) = σ 0 + exp{ [Y 1 (ω) cos(πx 2 ) + Y 3 (ω) sin(πx 2 )] e 1/8 where W, Y j [ 1, 1], j = 1, 2, 3, 4. +[Y 2 (ω) cos(πx 1 ) + Y 4 (ω) sin(πx 1 )] e 1/8 }
4 Random field modeling The values of the stochastic functions δ(t, ω) and σ(x, ω) are timely and spatially correlated in a way characterized by a covariance structure. The mean and the covariance of σ(x, ω) are nown respectively as σ 0 (x) = E(σ)(x) := σ(x, ω)dp(ω) and C σ (x, x ) = O O (σ(x, ω) σ 0 (x))(σ(x, ω) σ 0 (x ))dp(ω). C σ defines the ernel of a compact and self-adjoint operator.
5 Karhunen Loève random field representation The Karhunen Loève (KL) expansion of the random field σ(x, ω) is based on a spectral decomposition of the covariance ernel of the stochastic process. Ω C σ (x, x ) z j (x )dx = λ j z j (x), x Ω Define the following uncorrelated random variables Y j (ω) = 1 (σ(x, ω) σ 0 (x)) z j (x)dx, j = 1, 2,... λj Ω with zero mean and unit variance. The truncated KL expansion of σ(x, ω) is given by N σ N (x, ω) = σ 0 (x) + λj z j (x) Y j (ω) We assume the eigenvalues decay sufficiently fast. j=1
6 Optimal control of PDE with random data A model of the stochastic dynamical system A description of the control mechanism A criterion that models the purpose of the control (deterministic) and the cost of its action We have a constrained minimization problem minimize J(y, u) under the constraint c(y, u, ω) = 0 ω O
7 Optimality system For an event ω O, we have a deterministic control problem min u Uad J(y, u) := h(y) + ν g(u) s.t. c(y, u, ω) = 0 J : Y U R The existence of cy 1 enables a distinction between y, the state variable, and u U ad U, the optimization variable in the admissible set. So we have the mapping u J(y(u), u) in the form u IFT y(u) J(y(u), u) =: Ĵ(u) The solution of this optimization problem is characterized by the following optimality system c(y, u, ω) = 0 c y (y, u, ω) p = h (y) (ν g (u) + cu p, v u) 0 for all v U ad We have Ĵ(u) = ν g (u) + c u p(u).
8 Steady models with random inputs An explosive reaction-diffusion problem { y + δ e y = f in Ω y = g on Ω, δ is a reaction coefficient that is considered subject to randomness and is described by a spatially varying random field. Two realizations of δ(x 1, x 2, ω):
9 A singular optimal control problem We discuss the following min u L2 (Ω) J(y, u) y + δ(ω) e y = f + u in Ω, ω O y = 0 on Σ = Ω Control required to attain a desired target configuration y d (x) A cost functional of tracing type can be considered J(y, u) = 1 2 y y d 2 L 2 (Ω) + ν 2 u 2 L 2 (Ω)
10 Optimality system for the steady model Corresponding to ω O, the optimal solution is characterized by y + δ(ω) e y = u in Ω, p + δ(ω) e y p + (y y d ) = 0 in Ω, νu p = 0 in Ω, y = 0, p = 0 on Ω. Solving an optimal control problem with random coefficients means to determine the -dim space of controls for all ω O. This space represents the solution of the stochastic control problem.
11 Time-dependent models with random inputs t y + G δ (y) + σ y = u in Q = Ω (0, T ) y = y 0 in Ω {t = 0} y = 0 on Σ = Ω (0, T ) The nonlinear term G δ (y) models the reaction inetics for the state y and δ is the reaction parameter. Here σ is the diffusion coefficient. Two realizations of σ(x 1, x 2, ω):
12 Reaction-diffusion process controlled through source terms We discuss the following min u L2 (Q) J(y, u) t y + G δ(ω) (y) + σ(ω) y = u in Q = Ω (0, T ), ω O y = y 0 in Ω {t = 0} y = 0 on Σ = Ω (0, T ) Control required to trac a desired trajectory y d (x, t) reach a desired terminal state y T (x) A cost functional of tracing type can be considered J(y, u) = α 2 y y d 2 L 2 (Q) + β 2 y(, T ) y T 2 L 2 (Ω) + ν 2 u 2 L 2 (Q)
13 Optimality system Corresponding to the event ω O, the optimal solution is characterized by t y + G δ(ω) (y) + σ(ω) y = u in Q t p + G δ(ω) (y)p + σ(ω) p + α(y y d) = 0 in Q νu p = 0 in Q y = 0, p = 0 on Σ With initial condition y(x, 0) = y 0 (x) for the state variable (evolving forward in time). And the terminal condition for the adjoint variable (evolving bacward in time) p(x, T ) = β(y(x, T ) y T (x)). By solving optimal control problems with random coefficients, we explore the -dim space of controls for all configurations of the parameters space. The space of controls represents the solution of the stochastic control problem.
14 The issue of a robust control An important tas is to define a unique control for the governing random PDE model that provides good tracing features for all configurations of the random coefficients. Consider the following problem min u U E(Ĵ)(u) := min u U O Ĵ ω (u) dp(ω) where Ĵω(u) = Ĵ(y ω(u), u). A modeling issue... E(J(σ, δ)) = J(E(σ), E(δ)) J 2 σ 2 s2 σ J 2 δ 2 s2 δ J 2 σ δ s σs δ + O(s 4 ) with standard deviations (s σ, s δ ). J(y(u E(σ), E(δ) ), u E(σ), E(δ) ) J(y(E(u)), E(u)).
15 Collocation approach in the probability space In this framewor, the PDE model is approximated by FEM, FVM or FDM discretization and the stochastic domain is approximated using multidimensional interpolating functions. Stochastic optimal solutions are sought in a tensor-product space P(Γ) [Y U] where P represents the span of tensor-product interpolation polynomials in the N + M-dimensional stochastic space We have the random fields σ(x, ω) = σ(x, Y 1 (ω),..., Y N (ω)) and δ(t, ω) = δ(t, W 1 (ω),..., W M (ω)). ρ denotes the joint probability density of (Y 1,..., Y N, W 1,..., W M ) with support Γ = N j=1 ΓY j M j=1 ΓW j. N = N + M.
16 Lagrange interpolation Consider a function f on the one-dim space Γ j = [ 1, 1] with m j nodal values (z 1,..., z mj ). Lagrange interpolation operator m j m j I (f )(z) = f (z i ) L j i (z), Lj i (z) = i=1 =1, i z z z i z. Assume p mj (f ) is the best approximation polynomial of degree m j 1 f p mj (f ) f I (f ) (1 + Λ) f p mj (f ), mj Λ = max z Γj i=1 Lj i (z) is the Lebesgue constant depending on the distribution of the nodes. The Chebyshev Gauss Lobatto (CGL) nodes are given by ( ) (i 1)π z i = cos, i = 1,..., m j, m j 1 m j = 2 j (m 1 = 1). Here j > 1 is the order index for the grid with nodes z j = (z j 1, zj 2,..., zj m j ).
17 Multivariate interpolation In the multidimensional case the determination of the best approximation polynomial is an unsolved problem. For f : R N R the practice is to consider tensor-product interpolation polynomials on a full tensor-product grid. m 1 I(f )(z) =... i 1=1 m N i N =1 ( ) f (z j1 i 1,..., z j N in ) L j1 i 1... L j N in (z), Thus the mean value of a function f : R N R is given by m 1 E(f ) =... i 1=1 m N i N =1 and the variance is Var(f ) = E([f E(f )] 2 ). ( ) f (z j1 i 1,..., z j N in ) L j1 i 1... L j N in (z) ρ(z) dz, Γ
18 Smolya sparse-grids The full tensor-product interpolation formulae require m 1... m N function evaluations and this effort grows exponentially with the number of dimensions: Curse of dimensionality. A Smolya scheme provides multivariate interpolation as linear combination of tensor-product formulas. A full tensor-product N -dimensional grid of order J with j 1 = j 2 = j N = J in each dimension is given by ẑ J = z j1 z j2 z j N The sparse grid z J of order J is composed of a strict subset of full grids. For example, for N = 2 and J = 4 we have {z 1 z 4 } {z 2 z 3 } {z 3 z 2 } {z 4 z 1 }
19 Sparse-grid interpolation and integration The sparse grids are the union of all possible reduced full grids with orders that sum up to J z J = j =J +N 1 z j1 z j2 z j N with vector index notation j = (j 1,.., j N ), where j = j j N. On the sparse-grids we have an interpolation formula and we have a sparse-grid quadrature.
20 Approximation error in the probabilistic and physical spaces The optimality systems are discretized by finite differences and bacward Euler scheme. The discretization error of the deterministic solution y h y c h 2, p h p c h 2, and u h u c h 2. The stochastic solution is affected by an error having three components: 1) The error due to truncation of the KL expansion, e N ; 2) The discretization error (above), e d ; 3) The Smolya sparse grid interpolation error due to the collocation, e S c χ g/(1+log(2n )). (χ is the total number of sparse-grids points.) The total error is e total = e N + e d + e S.
21 CS-Multigrid - Sparse-grids solution process For i = 1,..., χ 1. Define coefficient functions δ i, σ i. 2. Call CSMG to solve the ith optimality system to Ĵ(u) < tol. 3. Collect results for statistics. End Statistics: compute E( y ), Var ( y ), Sew ( y ), and E(u).
22 Collective Smoothing Multi-Grid -V(m 1, m 2 )-Cycle Set B 1 (w (0) 1 ) A 1 1 (e.g., iterating with S 1 starting with w (0) 1 ). For = 2,..., L define B in terms of B 1 as follows. 1. Set the starting approximation w (0). 2. Pre-smoothing. Define w (l) w (l) for l = 1,..., m 1, by = S (w (l 1), f ) 3. Coarse grid correction. Set w (m1+1) = w (m1) + I 1 (w 1 Î 1 w 1 = B 1 (Î 1 w (m 1 ) h ) 4. Post-smoothing. Define w (l) w (l) I 1 5. Set B (w (0) ) f = w (m1+m2+1). w (m1) ) where (f A (w (m 1 ) )) + A 1 (Î 1 w (m 1 ) i ). for l = m 1 + 2,, m 1 + m 2 + 1, by = S (w (l 1), f )
23 Multigrid components: Smoothing Coupling between state and control variables. For time dependent problems preserve opposite orientation of state and adjoint equations. t = 0 y 0 y t = T p p 0 Pointwise smoother: Collective Gauss-Seidel Newton Iteration (CGSN), Time-Splitted Collective Gauss-Seidel Iteration (TS-CGS) Blocwise smoother: Time-Line Collective Gauss-Seidel Iteration (TL-CGS)
24 Mean MG Convergence factors -2-4 Log nu gamma Log nu gamma Convergence factors of TL-CGS scheme (left) and TS-CGS scheme E(η) TL-CGS TS-CGS γ \ ν The mean E(η) for TL-CGS and TS-CGS multigrid schemes Università (ν 1 = degli ν 2 Studi = 1); del Sannio, Italy
25 Elliptic case: experimental setting We consider a random reaction coefficient given by δ(x 1, x 2, ω) = δ 0 + exp{ [Y 1 (ω) cos(πx 2 ) + Y 3 (ω) sin(πx 2 )] e 1/8 +[Y 2 (ω) cos(πx 1 ) + Y 4 (ω) sin(πx 1 )] e 1/8 } where δ 0 = 2 and Y j [ 1, 1], j = 1, 2, 3, 4. This field is characterized by a squared exponential covariance typical of Gaussian processes The desired target configuration is given by y d (x 1, x 2 ) = 1 sin(π x 1 ) sin(π x 2 ). This target is not attainable by any control. Multigrid The and mean stochastic field reaction sparse-grids coefficient techniques E(δ) for PDE control problems with random coefficients
26 Numerical results for the elliptic case Table: Results with the CSMG sparse-grids scheme; y = y y d. J = 2, χ = 41 ν N x N y CPU(s) E( y ) Var ( y ) Sew ( y ) E(ρ obs ) J = 3, χ = 137 ν N x N y CPU(s) E( y ) Var ( y ) Sew ( y ) E(ρ obs )
27 Parabolic case: experimental setting We consider the case of a nonmonotone nonlinear reaction term given by G δ (y) = δ e y that is used to model explosive combustion phenomena. The desired target trajectory is given by y d (x, t) = (1 + t) (x 1 x 2 1 ) (x 2 x 2 2 ) cos(4π t). This is an oscillating function whose amplitude increases linearly with time. We tae y T (x) = y d (x, T ). We use one pre- and one post-smoothing steps (ν 1 = ν 2 = 1) and h = 1/4 is the coarsest space mesh size. We tae Ω = (0, 1) (0, 1) and T = 1. The stopping criteria is r y + r p /ν <
28 TS-CGS MG convergence and tracing properties Results for α = 1, β = 0, with TS-CGS MG scheme. y = y y d. J = 2, χ = 61 ν N x N y N t CPU(s) E( y ) Var ( y ) Sew ( y ) E(ρ obs ) J = 3, χ = 241 ν N x N y N t CPU(s) E( y ) Var ( y ) Sew ( y ) E(ρ obs )
29 TL-CGS MG convergence and tracing properties Results for α = 1, β = 0 and J = 2, χ = 61, with TL-CGS MG scheme. y = y y d. ν N x N y N t CPU(s) E( y ) Var ( y ) Sew ( y ) E(ρ obs )
30 Searching for a robust control We consider the function E(u) defined as the mean function of the optimal controls corresponding to each point of the coefficient configuration space. Another possible way to define a unique control function is to consider the control computed using average coefficient values. The question is then how good are y E(u) and y E(σ), E(δ) for tracing the desired target for different optimization settings.
31 Elliptic case: results with robust control candidates Table: Results with a given control (left) and moments of tracing ability for controls corresponding to different realization of the reaction coefficient; N x N y = , χ = 137. ν E( y y d ) Var ( y y d ) E( y ) Var ( y ) y E(u) y ue(δ) y E(u) y ue(δ) y E(u) y ue(δ)
32 Parabolic case: results with robust control candidates α = 1, β = 0, and ν = 10 6 ; N x N y N t = , χ = 241. E( y y d ) Var ( y y d ) Sew ( y y d ) y E(u) y ue(σ), E(δ) α = 0, β = 1, and ν = 10 6 ; N x N y N t = , χ = 241. E( y y d ) Var ( y y d ) Sew ( y y d ) y E(u) y ue(σ), E(δ) The difference in the tracing property of the two strategies is minimal thus motivating forthcoming investigation efforts.
33 A model reduction approach Define a robust control as û = argmin u L 2 (Q) O J(y ω (u), u) dp(ω) To approximately solve this problem we consider argmin u L2 (Q) J(y ω (u), u) dp(ω) argmin u L2 (Q) O χ w j J(y ωj (u), u) Further, we construct a POD basis of the optimal controls. Let u j (x, t), j = 1,..., χ represent the controls corresponding to different configurations, the POD basis functions are obtained from T A j = u j (x, t)u (x, t)dxdt, A = V SV. Ω 0 The th POD basis function is φ (x, t) = 1 s m < χ and s = S. j=1 m V j u j (x, t), where j=1
34 Minimization in the POD space We use the representation u(x, t) = m c φ (x, t) =1 and find min J(u) := χ j=1 w j J(y ωj (u), u). The gradient w.r.t the POD coefficients is given by c J(u) = Then we can use c (l+1) where u (l) = m χ w j (νc p ωj (u)) = νc j=1 =1 c(l) χ w j p ωj (u) j=1 = c (l) α c J(u (l) ), = 1,..., m φ (x, t).
Multilevel methods for PDE control problems
Multilevel methods for PDE control problems Alfio Borzì Università degli Studi del Sannio Dipartimento e Facoltà di Ingegneria, Benevento, Italia Alfio Borzì (Università degli Studi del Sannio, Multilevel
More informationMultigrid methods for pde optimization
Multigrid methods for pde optimization Alfio Borzì Università degli Studi del Sannio Dipartimento e Facoltà di Ingegneria, Benevento, Italia Alfio Borzì (Università degli Studi del Sannio, Multigrid Italy)
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 informationEfficient Solvers for Stochastic Finite Element Saddle Point Problems
Efficient Solvers for Stochastic Finite Element Saddle Point Problems Catherine E. Powell c.powell@manchester.ac.uk School of Mathematics University of Manchester, UK Efficient Solvers for Stochastic Finite
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationBang bang control of elliptic and parabolic PDEs
1/26 Bang bang control of elliptic and parabolic PDEs Michael Hinze (joint work with Nikolaus von Daniels & Klaus Deckelnick) Jackson, July 24, 2018 2/26 Model problem (P) min J(y, u) = 1 u U ad,y Y ad
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 informationMultilevel stochastic collocations with dimensionality reduction
Multilevel stochastic collocations with dimensionality reduction Ionut Farcas TUM, Chair of Scientific Computing in Computer Science (I5) 27.01.2017 Outline 1 Motivation 2 Theoretical background Uncertainty
More informationNumerical Analysis Comprehensive Exam Questions
Numerical Analysis Comprehensive Exam Questions 1. Let f(x) = (x α) m g(x) where m is an integer and g(x) C (R), g(α). Write down the Newton s method for finding the root α of f(x), and study the order
More informationMultilevel accelerated quadrature for elliptic PDEs with random diffusion. Helmut Harbrecht Mathematisches Institut Universität Basel Switzerland
Multilevel accelerated quadrature for elliptic PDEs with random diffusion Mathematisches Institut Universität Basel Switzerland Overview Computation of the Karhunen-Loéve expansion Elliptic PDE with uniformly
More 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 informationKarhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques
Institut für Numerische Mathematik und Optimierung Karhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques Oliver Ernst Computational Methods with Applications Harrachov, CR,
More 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 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 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 informationSimultaneous Robust Design and Tolerancing of Compressor Blades
Simultaneous Robust Design and Tolerancing of Compressor Blades Eric Dow Aerospace Computational Design Laboratory Department of Aeronautics and Astronautics Massachusetts Institute of Technology GTL Seminar
More 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 informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 24: Preconditioning and Multigrid Solver Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 5 Preconditioning Motivation:
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 informationNUMERICAL METHODS FOR ENGINEERING APPLICATION
NUMERICAL METHODS FOR ENGINEERING APPLICATION Second Edition JOEL H. FERZIGER A Wiley-Interscience Publication JOHN WILEY & SONS, INC. New York / Chichester / Weinheim / Brisbane / Singapore / Toronto
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 informationSolving the stochastic steady-state diffusion problem using multigrid
IMA Journal of Numerical Analysis (2007) 27, 675 688 doi:10.1093/imanum/drm006 Advance Access publication on April 9, 2007 Solving the stochastic steady-state diffusion problem using multigrid HOWARD ELMAN
More 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 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 informationThe Conjugate Gradient Method
The Conjugate Gradient Method Classical Iterations We have a problem, We assume that the matrix comes from a discretization of a PDE. The best and most popular model problem is, The matrix will be as large
More 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 information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 9
Spring 2015 Lecture 9 REVIEW Lecture 8: Direct Methods for solving (linear) algebraic equations Gauss Elimination LU decomposition/factorization Error Analysis for Linear Systems and Condition Numbers
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 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 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 informationOutline. 1 Boundary Value Problems. 2 Numerical Methods for BVPs. Boundary Value Problems Numerical Methods for BVPs
Boundary Value Problems Numerical Methods for BVPs Outline Boundary Value Problems 2 Numerical Methods for BVPs Michael T. Heath Scientific Computing 2 / 45 Boundary Value Problems Numerical Methods for
More informationIndex. higher order methods, 52 nonlinear, 36 with variable coefficients, 34 Burgers equation, 234 BVP, see boundary value problems
Index A-conjugate directions, 83 A-stability, 171 A( )-stability, 171 absolute error, 243 absolute stability, 149 for systems of equations, 154 absorbing boundary conditions, 228 Adams Bashforth methods,
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 informationModel order reduction & PDE constrained optimization
1 Model order reduction & PDE constrained optimization Fachbereich Mathematik Optimierung und Approximation, Universität Hamburg Leuven, January 8, 01 Names Afanasiev, Attwell, Burns, Cordier,Fahl, Farrell,
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 informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationConstrained Minimization and Multigrid
Constrained Minimization and Multigrid C. Gräser (FU Berlin), R. Kornhuber (FU Berlin), and O. Sander (FU Berlin) Workshop on PDE Constrained Optimization Hamburg, March 27-29, 2008 Matheon Outline Successive
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 informationarxiv: v2 [math.na] 8 Apr 2017
A LOW-RANK MULTIGRID METHOD FOR THE STOCHASTIC STEADY-STATE DIFFUSION PROBLEM HOWARD C. ELMAN AND TENGFEI SU arxiv:1612.05496v2 [math.na] 8 Apr 2017 Abstract. We study a multigrid method for solving large
More informationStochastic 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 informationAdaptive algebraic multigrid methods in lattice computations
Adaptive algebraic multigrid methods in lattice computations Karsten Kahl Bergische Universität Wuppertal January 8, 2009 Acknowledgements Matthias Bolten, University of Wuppertal Achi Brandt, Weizmann
More informationAspects of Multigrid
Aspects of Multigrid Kees Oosterlee 1,2 1 Delft University of Technology, Delft. 2 CWI, Center for Mathematics and Computer Science, Amsterdam, SIAM Chapter Workshop Day, May 30th 2018 C.W.Oosterlee (CWI)
More informationSuboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids
Suboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids Jochen Garcke joint work with Axel Kröner, INRIA Saclay and CMAP, Ecole Polytechnique Ilja Kalmykov, Universität
More informationSchwarz Preconditioner for the Stochastic Finite Element Method
Schwarz Preconditioner for the Stochastic Finite Element Method Waad Subber 1 and Sébastien Loisel 2 Preprint submitted to DD22 conference 1 Introduction The intrusive polynomial chaos approach for uncertainty
More informationarxiv: v1 [math.oc] 21 Jan 2019
An Efficient Parallel-in-Time Method for Optimization with Parabolic PDEs Sebastian Götschel Michael L. Minion arxiv:1901.06850v1 [math.oc] 21 Jan 2019 January 17, 2019 Abstract To solve optimization problems
More informationA sparse grid stochastic collocation method for partial differential equations with random input data
A sparse grid stochastic collocation method for partial differential equations with random input data F. Nobile R. Tempone C. G. Webster May 24, 2007 Abstract This work proposes and analyzes a Smolyak-type
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
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 informationOn Multigrid for Phase Field
On Multigrid for Phase Field Carsten Gräser (FU Berlin), Ralf Kornhuber (FU Berlin), Rolf Krause (Uni Bonn), and Vanessa Styles (University of Sussex) Interphase 04 Rome, September, 13-16, 2004 Synopsis
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 informationPOD for Parametric PDEs and for Optimality Systems
POD for Parametric PDEs and for Optimality Systems M. Kahlbacher, K. Kunisch, H. Müller and S. Volkwein Institute for Mathematics and Scientific Computing University of Graz, Austria DMV-Jahrestagung 26,
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 informationFundamentals of Data Assimila1on
014 GSI Community Tutorial NCAR Foothills Campus, Boulder, CO July 14-16, 014 Fundamentals of Data Assimila1on Milija Zupanski Cooperative Institute for Research in the Atmosphere Colorado State University
More informationMultigrid absolute value preconditioning
Multigrid absolute value preconditioning Eugene Vecharynski 1 Andrew Knyazev 2 (speaker) 1 Department of Computer Science and Engineering University of Minnesota 2 Department of Mathematical and Statistical
More informationarxiv: v1 [math.na] 14 Sep 2017
Stochastic collocation approach with adaptive mesh refinement for parametric uncertainty analysis arxiv:1709.04584v1 [math.na] 14 Sep 2017 Anindya Bhaduri a, Yanyan He 1b, Michael D. Shields a, Lori Graham-Brady
More informationSolving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners
Solving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners Eugene Vecharynski 1 Andrew Knyazev 2 1 Department of Computer Science and Engineering University of Minnesota 2 Department
More informationKarhunen-Loeve Expansion and Optimal Low-Rank Model for Spatial Processes
TTU, October 26, 2012 p. 1/3 Karhunen-Loeve Expansion and Optimal Low-Rank Model for Spatial Processes Hao Zhang Department of Statistics Department of Forestry and Natural Resources Purdue University
More informationIntroduction. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA
Introduction to Finite and Spectral Element Methods Using MATLAB Second Edition C. Pozrikidis University of Massachusetts Amherst, USA (g) CRC Press Taylor & Francis Group Boca Raton London New York CRC
More informationSparse Grids. Léopold Cambier. February 17, ICME, Stanford University
Sparse Grids & "A Dynamically Adaptive Sparse Grid Method for Quasi-Optimal Interpolation of Multidimensional Analytic Functions" from MK Stoyanov, CG Webster Léopold Cambier ICME, Stanford University
More information1. Fast Iterative Solvers of SLE
1. Fast Iterative Solvers of crucial drawback of solvers discussed so far: they become slower if we discretize more accurate! now: look for possible remedies relaxation: explicit application of the multigrid
More informationA multigrid method for large scale inverse problems
A multigrid method for large scale inverse problems Eldad Haber Dept. of Computer Science, Dept. of Earth and Ocean Science University of British Columbia haber@cs.ubc.ca July 4, 2003 E.Haber: Multigrid
More informationNon-linear least squares
Non-linear least squares Concept of non-linear least squares We have extensively studied linear least squares or linear regression. We see that there is a unique regression line that can be determined
More informationSUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS
SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS A. RÖSCH AND R. SIMON Abstract. An optimal control problem for an elliptic equation
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More informationNonlinear stochastic Galerkin and collocation methods: application to a ferromagnetic cylinder rotating at high speed
Nonlinear stochastic Galerkin and collocation methods: application to a ferromagnetic cylinder rotating at high speed Eveline Rosseel Herbert De Gersem Stefan Vandewalle Report TW 541, July 29 Katholieke
More informationThe HJB-POD approach for infinite dimensional control problems
The HJB-POD approach for infinite dimensional control problems M. Falcone works in collaboration with A. Alla, D. Kalise and S. Volkwein Università di Roma La Sapienza OCERTO Workshop Cortona, June 22,
More informationLecture 1: Center for Uncertainty Quantification. Alexander Litvinenko. Computation of Karhunen-Loeve Expansion:
tifica Lecture 1: Computation of Karhunen-Loeve Expansion: Alexander Litvinenko http://sri-uq.kaust.edu.sa/ Stochastic PDEs We consider div(κ(x, ω) u) = f (x, ω) in G, u = 0 on G, with stochastic coefficients
More informationLesson 4: Stationary stochastic processes
Dipartimento di Ingegneria e Scienze dell Informazione e Matematica Università dell Aquila, umberto.triacca@univaq.it Stationary stochastic processes Stationarity is a rather intuitive concept, it means
More informationSimulating with uncertainty : the rough surface scattering problem
Simulating with uncertainty : the rough surface scattering problem Uday Khankhoje Assistant Professor, Electrical Engineering Indian Institute of Technology Madras Uday Khankhoje (EE, IITM) Simulating
More informationProper Orthogonal Decomposition. POD for PDE Constrained Optimization. Stefan Volkwein
Proper Orthogonal Decomposition for PDE Constrained Optimization Institute of Mathematics and Statistics, University of Constance Joined work with F. Diwoky, M. Hinze, D. Hömberg, M. Kahlbacher, E. Kammann,
More informationAlgebraic Multigrid as Solvers and as Preconditioner
Ò Algebraic Multigrid as Solvers and as Preconditioner Domenico Lahaye domenico.lahaye@cs.kuleuven.ac.be http://www.cs.kuleuven.ac.be/ domenico/ Department of Computer Science Katholieke Universiteit Leuven
More informationA sparse grid stochastic collocation method for elliptic partial differential equations with random input data
A sparse grid stochastic collocation method for elliptic partial differential equations with random input data F. Nobile R. Tempone C. G. Webster June 22, 2006 Abstract This work proposes and analyzes
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Dept. of
More informationMultigrid finite element methods on semi-structured triangular grids
XXI Congreso de Ecuaciones Diferenciales y Aplicaciones XI Congreso de Matemática Aplicada Ciudad Real, -5 septiembre 009 (pp. 8) Multigrid finite element methods on semi-structured triangular grids F.J.
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Higher Order
More informationITERATIVE METHODS FOR NONLINEAR ELLIPTIC EQUATIONS
ITERATIVE METHODS FOR NONLINEAR ELLIPTIC EQUATIONS LONG CHEN In this chapter we discuss iterative methods for solving the finite element discretization of semi-linear elliptic equations of the form: find
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 informationIntroduction. Chapter One
Chapter One Introduction The aim of this book is to describe and explain the beautiful mathematical relationships between matrices, moments, orthogonal polynomials, quadrature rules and the Lanczos and
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 informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationProbabilistic Structural Dynamics: Parametric vs. Nonparametric Approach
Probabilistic Structural Dynamics: Parametric vs. Nonparametric Approach S Adhikari School of Engineering, Swansea University, Swansea, UK Email: S.Adhikari@swansea.ac.uk URL: http://engweb.swan.ac.uk/
More informationChapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS
Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS 5.1 Introduction When a physical system depends on more than one variable a general
More informationReview. Numerical Methods Lecture 22. Prof. Jinbo Bi CSE, UConn
Review Taylor Series and Error Analysis Roots of Equations Linear Algebraic Equations Optimization Numerical Differentiation and Integration Ordinary Differential Equations Partial Differential Equations
More informationMultigrid Method ZHONG-CI SHI. Institute of Computational Mathematics Chinese Academy of Sciences, Beijing, China. Joint work: Xuejun Xu
Multigrid Method ZHONG-CI SHI Institute of Computational Mathematics Chinese Academy of Sciences, Beijing, China Joint work: Xuejun Xu Outline Introduction Standard Cascadic Multigrid method(cmg) Economical
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 informationKey words. PDE-constrained optimization, deferred correction, time-dependent PDE, coupled system
A RATIONAL DEFERRED CORRECTION APPROACH TO PDE-CONSTRAINED OPTIMIZATION STEFAN GÜTTEL AND JOHN W PEARSON Abstract The accurate and efficient solution of time-dependent PDE-constrained optimization problems
More informationPreface to the Second Edition. Preface to the First Edition
n page v Preface to the Second Edition Preface to the First Edition xiii xvii 1 Background in Linear Algebra 1 1.1 Matrices................................. 1 1.2 Square Matrices and Eigenvalues....................
More informationNonlinear Diffusion. Journal Club Presentation. Xiaowei Zhou
1 / 41 Journal Club Presentation Xiaowei Zhou Department of Electronic and Computer Engineering The Hong Kong University of Science and Technology 2009-12-11 2 / 41 Outline 1 Motivation Diffusion process
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 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 informationNumerical Methods for Engineers and Scientists
Numerical Methods for Engineers and Scientists Second Edition Revised and Expanded Joe D. Hoffman Department of Mechanical Engineering Purdue University West Lafayette, Indiana m MARCEL D E К К E R MARCEL
More informationIntroduction LECTURE 1
LECTURE 1 Introduction The source of all great mathematics is the special case, the concrete example. It is frequent in mathematics that every instance of a concept of seemingly great generality is in
More informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
More informationA Vector-Space Approach for Stochastic Finite Element Analysis
A Vector-Space Approach for Stochastic Finite Element Analysis S Adhikari 1 1 Swansea University, UK CST2010: Valencia, Spain Adhikari (Swansea) Vector-Space Approach for SFEM 14-17 September, 2010 1 /
More informationON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * 1. Introduction
Journal of Computational Mathematics, Vol.6, No.6, 008, 85 837. ON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * Tao Tang Department of Mathematics, Hong Kong Baptist
More informationEXPLORING THE LOCALLY LOW DIMENSIONAL STRUCTURE IN SOLVING RANDOM ELLIPTIC PDES
MULTISCALE MODEL. SIMUL. Vol. 5, No. 2, pp. 66 695 c 27 Society for Industrial and Applied Mathematics EXPLORING THE LOCALLY LOW DIMENSIONAL STRUCTURE IN SOLVING RANDOM ELLIPTIC PDES THOMAS Y. HOU, QIN
More informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More 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 informationwhere r n = dn+1 x(t)
Random Variables Overview Probability Random variables Transforms of pdfs Moments and cumulants Useful distributions Random vectors Linear transformations of random vectors The multivariate normal distribution
More information