Numerical Solution of PDEs: Bounds for Functional Outputs and Certificates
|
|
- Jonathan Elliott
- 5 years ago
- Views:
Transcription
1 Numerical Solution of PDEs: Bounds for Functional Outputs and Certificates J. Peraire Massachusetts Institute of Technology, USA ASC Workshop, August 2005, Stanford University
2 People MIT: A.T. Patera, A.M. Budge, H. Ciria, and J. Wong UCS: J. Bonet UPC: N. Pares, A. Huerta NUS: Z.C. Xuan ASC Workshop, August
3 Acknowledgments Singapore-MIT Alliance Sandia National Laboratories DARPA/AFOSR ASC Workshop, August
4 Sleipner Platform Failure Sank in August 1991, causing an event registering 3.0 on the Richter scale and leaving nothing but a pile of debris at a depth of 220m Sinking traced to a failure of a concrete tricell FEM performed with NASTRAN underestimated shear stresses by 47% More precise simulation of underdesigned component predicted failure at 62m Actually sank at 65m ASC Workshop, August
5 Basic Question How do we know if the answer computed with a FE code is correct 1? given that: the solution may not be well behaved we may not have similar solutions to compare we may not have access to the source code the code may no longer exist!! 1 i.e. consistent with the mathematical model ASC Workshop, August
6 Basic Question How do we know if the answer computed with a FE code is correct? Provide a Certificate ASC Workshop, August
7 Certificates What is a Certificate? A data set that documents a given claim Can be used to rigorously proof correctness Simple to exercise Stand alone - access to the code used to compute it not required The stronger the claim the longer the certificate (usually) ASC Workshop, August
8 Certificates Current Paradigm COMPUTATIONAL CODE SOLUTION UNCERTAINTY... - ALGORITHM - SOFTWARE ASC Workshop, August
9 Certificates Proposed Paradigm COMPUTATIONAL CODE SOLUTION + CERTIFICATE + A-POSTERIORI BOUNDS ASC Workshop, August
10 Certificates Proposed Paradigm COMPUTATIONAL CODE SOLUTION + CERTIFICATE CERTAINTY GUARANTEED!! + A-POSTERIORI BOUNDS VERIFICATION TEST ASC Workshop, August
11 Certificates Disclaimer Physical Phenomenon Modelling Uncertainty Continuous Mathematical Model Discretization Uncertainty Discrete Mathematical Model Software Uncertainty Prediction ASC Workshop, August
12 Certificates Disclaimer Physical Phenomenon Modelling Uncertainty Continuous Mathematical Model Prediction with certified error bounds ASC Workshop, August
13 Certificates Examples Polynomial Bounds Given a polynomial F (x), x IR n Claim : F (x) γ, x Certificate : Polynomials f 1 (x),..., f m (x) s.t. or ( F (x) γ = n i=1 f 2 i m i=1 (x)) (F (x) γ) = m f 2 i (x) (SOS) i=n+1 f 2 i (x) ASC Workshop, August
14 Certificates Examples Bounds for solutions of IVP... Given ẋ = f(x, t), x(0) = x 0, (f(x, t) polynomial) Claim : x(t ) γ Certificate : Polynomial function B(x, t) s.t. B t (x, t) + B x (x, t)f(x, t) 0, x, t B(x T, T ) > B(x 0, 0), x T γ Parrilo, Doyle,... ASC Workshop, August
15 Certificates Examples...Bounds for solutions of IVP... B t (x, t) + B x (x, t)f(x, t) 0, x, t B(x T, T ) > B(x 0, 0), x T γ ASC Workshop, August
16 Certificates Examples...Bounds for solutions of IVP... B t (x, t) + B x (x, t)f(x, t) 0, x, t B(x T, T ) > B(x 0, 0), x T γ ASC Workshop, August
17 Certificates Examples...Bounds for solutions of IVP Given: ẋ = px 3 Barrier function B(x,t) 3 Separating level set x(0) [0.85, 0.95] x p [0.05, 0.2] x t t? x(2) [2.0, 2.5] x(2) / [2.0, 2.5] ASC Workshop, August
18 Objective Compute Certificates for Bounds of Outputs of PDE s Work with quantities of interest Work with equations of interest Guarantee certainty even for low cost Cost effective ASC Workshop, August
19 Objective Examples Elasticity Non-regular solution (Plane Stress) OUTPUTS l(u) = l(u) = Γ O u ds Γ 5 t x ds ASC Workshop, August
20 What can we do Today? Linear Functional Outputs for: - Convection-Diffusion-Reaction Equation (high Pe) - Linear Elasticity Equations - Stokes Equations (DG) Collapse Loads in Limit Analysis (SOCP) Energy Release Rates in Linear Elasticity ASC Workshop, August
21 Outline Problem Description Function Minimization/Duality in IR n Method Overview (for CG) 1.- Bounds for Energy 2.- Bounds for Arbitrary Outputs 3.- Bounds for Arbitrary Equations 4.- Domain Decomposition (Hybridization) Method Summary and Examples Extension to a non-linear Convex Problem: Limit Analysis ASC Workshop, August
22 Problem Description Let u(x) X, x IR d, be the solution of a PDE A u = f. e.g. A 2, 2 + U, etc. We are typically interested in outputs of the form s = l(u) IR, e.g. l(v) v(x 0 ), l(v) = v x dx,... ASC Workshop, August
23 Problem Description u(x) is not computable ( dimensional) In practice, we compute approximation ū(x), such that u ū = C( 0) (as cost increases ). - For a given ū, C is unknown, and, any output approximation s = l(ū), is uncertain. Existing error estimates are either, - certain but uncomputable, or, - computable but uncertain. ASC Workshop, August
24 Problem Description Approach Compute Strict upper and lower bounds for functional outputs of the Exact solutions of PDE s... and give Certificates ASC Workshop, August
25 Function Minimization f : Y IR Unconstrained e.g. Y IR n s = min v Y f(v) f(v ) ASC Workshop, August
26 Function Minimization Unconstrained Upper Bound s = min v Y f(v) f(v ) f( v) }{{} upper bound EASY! s +, v Y... we expect a reasonable bound if v v ASC Workshop, August
27 Function Minimization Unconstrained Upper Bound: Example... Linear-Quadratic Problem f(v) = 1 2 vt Av v T b, A IR n IR n SPD, b, v IR n Stationarity condition Av = b s = min v Y f(v) = f(v ) = 1 2 v T Av ASC Workshop, August
28 Function Minimization Unconstrained...Upper Bound: Example Let v = v + ɛ s + = f( v) = f(v ) ɛt Aɛ Error is O(ɛ 2 ) ASC Workshop, August
29 Function Minimization Unconstrained Lower Bound s = min v Y f(v) f(v ) It appears that we need to know v explicitly Not generally available?? ASC Workshop, August
30 Function Minimization f : Y IR, g : Y Λ Equality Constrained objective constraints e.g. Y IR n, Λ IR m s = min v Y g(v) = 0 f(v) ASC Workshop, August
31 Function Minimization Equality Constrained Lagrangian L(v, λ) : Y Λ IR L(v, λ) = f(v) + λ T g(v) It follows that λ : Lagrange multipliers s = min v Y max λ Λ L(v, λ) ASC Workshop, August
32 Function Minimization Equality Constrained Duality min v Y max λ Λ L(v, λ) max λ Λ min v Y L(v, λ) Proof: Let Then, L(v, λ (v)) max λ Λ L(v (λ), λ) min v Y L(v, λ) L(v, λ) L(v, λ (v)) L(v, λ) L(v (λ), λ), v, λ. ASC Workshop, August
33 Function Minimization Equality Constrained Lower Bound Since by duality, s = min v Y Then, s max λ Λ max λ Λ L(v, λ) max λ Λ G(λ) G( λ) }{{} lower bound min L(v, λ) } v Y {{} L(v (λ), λ) G(λ) Dual Function s, λ Λ... reasonable bounds to be expected when λ λ ASC Workshop, August
34 Function Minimization Unconstrained Lower Bound... s = min v Y f(v) Assume f(v) F (v, g(v)) and consider s = min v Y, w Λ g(v) w = 0 F (v, w) ASC Workshop, August
35 Function Minimization Unconstrained...Lower Bound Lagrangian L(v, w, λ) = F (v, w) + λ T (g(v) w) Dual function G(λ) = min L(v, w, λ) = v Y,w Λ L(v (λ), w (λ), λ) s = min v Y f(v) max λ Λ G(λ) G( λ) s, λ ASC Workshop, August
36 Function Minimization Unconstrained Lower Bound: Example... Linear-Quadratic Problem f(v) = v T K T Kv + v T b, K IR n IR n, b, v IR n f(v) F (v, g(v)) = g T g + v T b, g(v) = Kv IR n L(v, w, λ) = w T w + b T v + λ T (Kv w) ASC Workshop, August
37 Function Minimization Unconstrained...Lower Bound: Example Dual Function G(λ) = min v Y,w Λ L(v, w, λ) = { 1 2 λt λ if K T λ = b otherwise s = min v Y f(v) max G(λ) λ Λ 1 2 λ T λ s, λ s.t. Kλ = b ASC Workshop, August
38 Method Overview 1.- Energy s = J(u) Poisson s Equation: Find u X() 2 u = f(x), x, (+ b.c. s) Energy functional: J(v) : X IR J(v) = v v dx 2 fv dx ASC Workshop, August
39 Method Overview 1.- Energy s = J(u) Minimization Minimization formulation min v X J(v) = J(u) = uf dx ASC Workshop, August
40 Method Overview 1.- Energy s = J(u) Upper Bound Upper bound s + J(u h ), u h X h X ( trivial ) ASC Workshop, August
41 Method Overview 1.- Energy s = J(u) Lower Bound... Lower bound s ( harder ) Construct dual problem (J(u) =) J c (p) = max q Q f J c (q), ASC Workshop, August
42 Method Overview 1.- Energy s = J(u)...Lower Bound... s = min v X = min max v X q Q max min q Q v X = max q Q f Q f = {q Q ( v v 2vf) dx (q = v) ( q q + 2q v 2vf) dx ( q q + 2q v 2vf) dx q q dx q v dx = q=f fv dx, v X }{{} } ASC Workshop, August
43 Method Overview 1.- Energy s = J(u)...Lower Bound... s = min v X = min max v X q Q max min q Q v X = max q Q f Q f = {q Q J(v) J c (q) ( q q + 2q v 2vf) dx ( q q + 2q v 2vf) dx q v dx = q=f fv dx, v X }{{} } ASC Workshop, August
44 Method Overview 1.- Energy s = J(u)...Lower Bound... or, in a different way... q q dx 2 q q dx 2 q v dx + fv dx + }{{}}{{} (q v) 2 dx 0, v X, q Q v v dx 0, v X, q Q v v dx 0, v X, q Q f J c (q) + J(v) 0, v X, q Q f Q f = {q Q q v dx = fv dx, v X} ( q = f) J(v) J c (q), v X, q Q f ASC Workshop, August
45 Method Overview 1.- Energy s = J(u)...Lower Bound... Duality ASC Workshop, August
46 Method Overview 1.- Energy s = J(u)...Lower Bound... Then, s J c (p h ), p h (Q f ) h Q f. ASC Workshop, August
47 Method Overview 1.- Energy s = J(u)...Lower Bound Idea : We can exchange an infinite dimensional minimization problem by a finite dimensional feasibility problem while retaining the bounding property ASC Workshop, August
48 Method Overview 1.- Energy s = J(u) Lower Bound - Summary Claim : s = J(u) = Given 2 u = f(x) uf dx s Certificate : Any p h (Q f ) h Q f s.t. s J c (p h ) Recall: Q f = {q Q q v dx = fv dx, v X} ( q = f) ASC Workshop, August
49 Method Overview 2.- General Outputs s = l(u) Find s = l(u), where u X() (l(v) = or, 2 u = f(x), x, (+ b.c. s) ( u v fv) dx = 0, Modified Energy : E(v) : X IR E(v) v v dx v X f O v dx) fv dx E(u) = 0 ASC Workshop, August
50 Method Overview 2.- General Outputs s = l(u) Lagrangian s = l(u) = min v X ( v ψ fψ) dx = 0, ψ X l(v) +E(v) Lagrangian : L(v, ψ) : X X IR L(v, ψ) = E(v) + l(v) + ( v ψ fψ) dx s = l(u) = min v max ψ L(v, ψ) ASC Workshop, August
51 Method Overview 2.- General Outputs s = l(u) Lower Bound... Weak duality + Relaxation s = l(u) = min v max ψ L(v, ψ) max ψ min v L(v, ψ) min v L(v, ψ), ψ X ASC Workshop, August
52 Method Overview L(v, ψ) = 2.- General Outputs s = l(u) v v dx + l(v) +...Lower Bound... fv dx ( v ψ f ψ) dx For a given ψ, L(v, ψ), contains quadratic and linear terms in v identical to J(v) (for an appropriate f ψ). L(v, ψ) = v v dx 2 f ψv dx f ψ dx ASC Workshop, August
53 Method Overview 2.- General Outputs s = l(u)...lower Bound Idea : Write output as a constrained minimization problem. Relax constraint to obtain an energy-like minimization problem. Obtain lower bound by finding a feasible solution of the dual problem. ASC Workshop, August
54 Method Overview 2.- General Outputs s = l(u) Upper Bound Define l (v) = l(v) and compute, s l (u) Idea: s + s l (u) = l(u) Upper Bound for l(v) Lower Bound for l(v) ASC Workshop, August
55 Method Overview 2.- General Outputs s = l(u) Summary Given 2 u = f(x) Claim : s + s = l(u) s Certificate : ψ Xh X, p + h (Q f +) h Q f +, p h (Q f ) h Q f ASC Workshop, August
56 Method Overview 3.- Non-symmetric equations or, 2 u + U u = f(x), x, (+ b.c. s) ( u v + (U u)v fv) dx = 0, Modified Energy : E(v) : X IR E(v) v v dx v X fv dx E(u) = 0 ASC Workshop, August
57 Method Overview 3.- Non-symmetric equations Lagrangian... s = l(u) = min l(v) +E(v) v X ( v ψ+(u v)ψ fψ) dx = 0, ψ X Lagrangian : L(v, ψ) : X X IR L(v, ψ) = E(v) + l(v) + ( v ψ + (U v)ψ fψ) dx s = l(u) = min v max ψ L(v, ψ)... ASC Workshop, August
58 Method Overview 3.- Non-symmetric equations...lagrangian Idea : Non-symmetric terms do not contribute to the energy and only enter in the Lagrangian linearly. After relaxation, minimization problem retains convex structure. ASC Workshop, August
59 Method Overview 4.- Domain Decomposition Recall that a lower bound for s = J(u), is given by s = J c (q), q Q f Q f = {q Q q v dx = fv dx, v X} i.e. find q Q s.t. q = f, in... not trivial ASC Workshop, August
60 Method Overview 4.- Domain Decomposition v X() continuous ˆv ˆX() discontinuous X ˆX ASC Workshop, August
61 Method Overview 4.- Domain Decomposition v X() continuous ˆv ˆX() discontinuous Idea : solve local minimization problems use piecewise polynomial certificates ASC Workshop, August
62 Method Overview 4.- Domain Decomposition Continuity Form Jump bilinear form b : ˆX() Λ(Γ) IR b(ˆv, λ) = γ γ [ˆv] γ λ γ ds, X {ˆv ˆX b(ˆv, λ) = 0, λ Λ} ASC Workshop, August
63 Method Overview 4.- Domain Decomposition Relaxation... J(u) = min v X J(v) = min ˆv ˆX b(ˆv, λ) = 0, λ Λ J(ˆv) = min ˆv ˆX max λ Λ J(ˆv) + b(ˆv, λ) min ˆv ˆX J(ˆv) + b(ˆv, λ), λ Λ (equil.) } {{ } J λ(ˆv) Solve local minimization problems ASC Workshop, August
64 Method Overview 4.- Domain Decomposition...Relaxation Feasibility problem over each triangle q = f, in T e q n = λ, on T e has an explicit solution provided f and λ are polynomial Ladeveze... ASC Workshop, August
65 Method Overview Summary 1. Primal problem: u h X h Au h = f 2. Dual problem: ψ Xh A ψ = f O, (l(v) = f O v dx) ASC Workshop, August
66 Method Overview Summary 3. Domain decomposition (Equilibration) λ Global Solution Equilibrated Solution ASC Workshop, August
67 Method Overview Summary 4. Obtain lower bounds for local minimization problems s + s... and piecewise polynomial certificates 5. It can be shown that the bound gap can be written as s + s = with e 0 T e T H e... Adaptivity ASC Workshop, August
68 Examples Convection-Diffusion f =10 3 ν 2 u + U u = f f O =1 s = l(u) = f O u dx ASC Workshop, August
69 Examples Convection-Diffusion u: ψ: -7.1x x x10-04 Solution Adjoint ASC Workshop, August
70 Examples Convection-Diffusion Adaptive Solution gap = s = ± Gap: Elements Uniform refinement would require 6356 elements ASC Workshop, August
71 Examples Stokes Flows 1 1/2 1 1/2 σ + (p + p) = f u = 0 U X 2 X 1 w l(u) = p = 1 x 1 cylinder t x ds ASC Workshop, August
72 Examples Stokes Flows DG Results 1 Pressure h S S + 1/ / / Special attention paid to incompressibility constraint!! ASC Workshop, August
73 Examples Linear Elasticity Energy Release Rates... p Total Potential Energy Π(v) = 1 a(v, v) (f, v) g, v 2 Displacement solution u minimizes Π(v) 5 Crack tip 60 Π(u) = 1 2 a(u, u) = 1 2 u 30 Energy Release Rate J(u) p 20 δπ(u) = J (u) δl... l crack length ASC Workshop, August
74 Examples Linear Elasticity...Energy Release Rates... Mixed mode crack problem (Plane Strain, ν = 0.3) p Crac k tip O Crac k tip p ASC Workshop, August
75 Examples Linear Elasticity...Energy Release Rates... ASC Workshop, August
76 Examples Linear Elasticity...Energy Release Rates Mesh size H H/2 H/4 H/8 H/16 J (u H ) η χ e J J ASC Workshop, August
77 Nonlinear Extension Limit Analysis Compute Bounds on the Collapse Load under the assumption of rigid-plastic material behavior Linear vs. Rigid-Plastic ASC Workshop, August
78 Nonlinear Extension Limit Analysis Limit Analysis vs Non-linear Analysis CONS Limited Physics PROS Existence of solution Uniqueness of solution Computable Certifiable... ASC Workshop, August
79 Nonlinear Extension a(σ, v) = F (v) = fv dx + σ : ε(v) dx gv ds X F = {v X F (v) = 1} Σ = {σ f(σ) σ Y } ε(v) = { 0 if f(σ) < σy κ f σ if f(σ) = σ Y Limit Analysis ϕ = max σ Σ Formulation max σ Σ a(σ, v) = ϕf (v), v X = min max v X F σ Σ = max σ Σ a(σ, v) min a(σ, v) v X F ϕ a(σ, v) Upper Bound min a( σ, v) Lower Bound v X F ASC Workshop, August
80 Nonlinear Extension Limit Analysis Outline By choosing appropriate piecewise polynomial interpolations for v and σ we can obtain strict upper and lower bounds on ϕ Discrete minimization/maximization problems are convex (SOCP) and solved (globally) with an IPM ϕ + ϕ can be decomposed into elemental contributions Adaptivity ASC Workshop, August
81 Nonlinear Extension Limit Analysis Examples... Cantilever Beam in Plane Stress 1.6 g = 0 1 x x ASC Workshop, August
82 Nonlinear Extension Limit Analysis...Examples... UNIFORM ADAPTED ASC Workshop, August
83 Nonlinear Extension Limit Analysis...Examples... Uniform Mesh Number Number Low. Bound Upp. Bound Bound Low. Bound Upp. Bound of refin. of elem. λ LB h λ UB h Gap h Error (%) Error (%) Adaptive Mesh Number Number Low. Bound Upp. Bound Bound Low. Bound Upp. Bound of refin. of elem. λ LB h λ UB h Gap h Error (%) Error (%) ASC Workshop, August
84 Conclusions Uniform bounds on Relevant engineering outputs (linear functionals) of Exact weak solutions of linear PDEs, with a Stand-alone certificate of precision, including Non-symmetric operators, using Standard FE solutions and purely local subproblems. ASC Workshop, August
85 Conclusions Certificates allow to Standardize the use of more accurate and safer mathematical models (e.g. construction codes) Eliminate costlier-than-necessary computations Allow for true black boxes that can be used by non-experts in numerical analysis Document computations Address software error issues... ASC Workshop, August
86 Current Work Exploit Discontinuous Galerkin Discretizations Time dependent parabolic problems µ PDE s Non-coercive operators with positivity constraints on the solution Deformation theory of plasticity Recent papers can be found at: ASC Workshop, August
Computing Upper and Lower Bounds for the J integral in Two-Dimensional Linear Elasticity
Computing Upper and Lower Bounds for the J integral in Two-Dimensional Linear Elasticity Z.C. Xuan a 1 N. Parés b J. Peraire c 2 a Department of Computer Science, Tianjin University of Technology and Education,
More informationThe computation of bounds for linear-functional outputs of weak solutions to the two-dimensional elasticity equations
The computation of bounds for linear-functional outputs of weak solutions to the two-dimensional elasticity equations N. Parés a J. Bonet b A. Huerta a J. Peraire c 1 a Laboratori de Càlcul Numèric, Universitat
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 informationMixed Finite Element Methods. Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016
Mixed Finite Element Methods Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016 Linear elasticity Given the load f : Ω R n, find the displacement u : Ω R n and the
More information1 The Stokes System. ρ + (ρv) = ρ g(x), and the conservation of momentum has the form. ρ v (λ 1 + µ 1 ) ( v) µ 1 v + p = ρ f(x) in Ω.
1 The Stokes System The motion of a (possibly compressible) homogeneous fluid is described by its density ρ(x, t), pressure p(x, t) and velocity v(x, t). Assume that the fluid is barotropic, i.e., the
More informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
More informationApproximation of fluid-structure interaction problems with Lagrange multiplier
Approximation of fluid-structure interaction problems with Lagrange multiplier Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 30, 2016 Outline
More informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
More informationSpace-time Finite Element Methods for Parabolic Evolution Problems
Space-time Finite Element Methods for Parabolic Evolution Problems with Variable Coefficients Ulrich Langer, Martin Neumüller, Andreas Schafelner Johannes Kepler University, Linz Doctoral Program Computational
More informationCLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE
CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that
More informationMIXED FINITE ELEMENTS FOR PLATES. Ricardo G. Durán Universidad de Buenos Aires
MIXED FINITE ELEMENTS FOR PLATES Ricardo G. Durán Universidad de Buenos Aires - Necessity of 2D models. - Reissner-Mindlin Equations. - Finite Element Approximations. - Locking. - Mixed interpolation or
More informationOptimization for Machine Learning
Optimization for Machine Learning (Problems; Algorithms - A) SUVRIT SRA Massachusetts Institute of Technology PKU Summer School on Data Science (July 2017) Course materials http://suvrit.de/teaching.html
More informationBACKGROUNDS. Two Models of Deformable Body. Distinct Element Method (DEM)
BACKGROUNDS Two Models of Deformable Body continuum rigid-body spring deformation expressed in terms of field variables assembly of rigid-bodies connected by spring Distinct Element Method (DEM) simple
More informationLecture Note 5: Semidefinite Programming for Stability Analysis
ECE7850: Hybrid Systems:Theory and Applications Lecture Note 5: Semidefinite Programming for Stability Analysis Wei Zhang Assistant Professor Department of Electrical and Computer Engineering Ohio State
More informationUsing Extended Interval Algebra in Discrete Mechanics Fulvio Tonon University of Texas at Austin, USA
Using Extended Interval Algebra in Discrete Mechanics Fulvio Tonon University of Texas at, USA Traditional approach to modeling physical systems 1. Model a physical phenomenon using a differential equation
More informationFinite Elements for Elastic Shell Models in
Elastic s in Advisor: Matthias Heinkenschloss Computational and Applied Mathematics Rice University 13 April 2007 Outline Elasticity in Differential Geometry of Shell Geometry and Equations The Plate Model
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationLagrangian Duality Theory
Lagrangian Duality Theory Yinyu Ye Department of Management Science and Engineering Stanford University Stanford, CA 94305, U.S.A. http://www.stanford.edu/ yyye Chapter 14.1-4 1 Recall Primal and Dual
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 informationWELL POSEDNESS OF PROBLEMS I
Finite Element Method 85 WELL POSEDNESS OF PROBLEMS I Consider the following generic problem Lu = f, where L : X Y, u X, f Y and X, Y are two Banach spaces We say that the above problem is well-posed (according
More informationSolving large Semidefinite Programs - Part 1 and 2
Solving large Semidefinite Programs - Part 1 and 2 Franz Rendl http://www.math.uni-klu.ac.at Alpen-Adria-Universität Klagenfurt Austria F. Rendl, Singapore workshop 2006 p.1/34 Overview Limits of Interior
More informationBoundary Value Problems and Iterative Methods for Linear Systems
Boundary Value Problems and Iterative Methods for Linear Systems 1. Equilibrium Problems 1.1. Abstract setting We want to find a displacement u V. Here V is a complete vector space with a norm v V. In
More informationGustafsson, Tom; Stenberg, Rolf; Videman, Juha A posteriori analysis of classical plate elements
Powered by TCPDF (www.tcpdf.org) This is an electronic reprint of the original article. This reprint may differ from the original in pagination and typographic detail. Gustafsson, Tom; Stenberg, Rolf;
More informationKarush-Kuhn-Tucker Conditions. Lecturer: Ryan Tibshirani Convex Optimization /36-725
Karush-Kuhn-Tucker Conditions Lecturer: Ryan Tibshirani Convex Optimization 10-725/36-725 1 Given a minimization problem Last time: duality min x subject to f(x) h i (x) 0, i = 1,... m l j (x) = 0, j =
More informationConstitutive models. Constitutive model: determines P in terms of deformation
Constitutive models Constitutive model: determines P in terms of deformation Elastic material: P depends only on current F Hyperelastic material: work is independent of path strain energy density function
More information1 Discretizing BVP with Finite Element Methods.
1 Discretizing BVP with Finite Element Methods In this section, we will discuss a process for solving boundary value problems numerically, the Finite Element Method (FEM) We note that such method is a
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 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 informationAM 205: lecture 19. Last time: Conditions for optimality, Newton s method for optimization Today: survey of optimization methods
AM 205: lecture 19 Last time: Conditions for optimality, Newton s method for optimization Today: survey of optimization methods Quasi-Newton Methods General form of quasi-newton methods: x k+1 = x k α
More informationA primer on Numerical methods for elasticity
A primer on Numerical methods for elasticity Douglas N. Arnold, University of Minnesota Complex materials: Mathematical models and numerical methods Oslo, June 10 12, 2015 One has to resort to the indignity
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
More informationConvex Optimization. Lecture 12 - Equality Constrained Optimization. Instructor: Yuanzhang Xiao. Fall University of Hawaii at Manoa
Convex Optimization Lecture 12 - Equality Constrained Optimization Instructor: Yuanzhang Xiao University of Hawaii at Manoa Fall 2017 1 / 19 Today s Lecture 1 Basic Concepts 2 for Equality Constrained
More informationLecture 18: Optimization Programming
Fall, 2016 Outline Unconstrained Optimization 1 Unconstrained Optimization 2 Equality-constrained Optimization Inequality-constrained Optimization Mixture-constrained Optimization 3 Quadratic Programming
More informationEE 227A: Convex Optimization and Applications October 14, 2008
EE 227A: Convex Optimization and Applications October 14, 2008 Lecture 13: SDP Duality Lecturer: Laurent El Ghaoui Reading assignment: Chapter 5 of BV. 13.1 Direct approach 13.1.1 Primal problem Consider
More informationFEniCS Course. Lecture 0: Introduction to FEM. Contributors Anders Logg, Kent-Andre Mardal
FEniCS Course Lecture 0: Introduction to FEM Contributors Anders Logg, Kent-Andre Mardal 1 / 46 What is FEM? The finite element method is a framework and a recipe for discretization of mathematical problems
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationLecture 13: Constrained optimization
2010-12-03 Basic ideas A nonlinearly constrained problem must somehow be converted relaxed into a problem which we can solve (a linear/quadratic or unconstrained problem) We solve a sequence of such problems
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Variational Problems of the Dirichlet BVP of the Poisson Equation 1 For the homogeneous
More informationNew DPG techniques for designing numerical schemes
New DPG techniques for designing numerical schemes Jay Gopalakrishnan University of Florida Collaborator: Leszek Demkowicz October 2009 Massachusetts Institute of Technology, Boston Thanks: NSF Jay Gopalakrishnan
More information4. Algebra and Duality
4-1 Algebra and Duality P. Parrilo and S. Lall, CDC 2003 2003.12.07.01 4. Algebra and Duality Example: non-convex polynomial optimization Weak duality and duality gap The dual is not intrinsic The cone
More informationNumerical Methods for Partial Differential Equations
Numerical Methods for Partial Differential Equations Eric de Sturler University of Illinois at Urbana-Champaign Read section 8. to see where equations of type (au x ) x = f show up and their (exact) solution
More informationA posteriori error estimation for elliptic problems
A posteriori error estimation for elliptic problems Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in
More informationSpline Element Method for Partial Differential Equations
for Partial Differential Equations Department of Mathematical Sciences Northern Illinois University 2009 Multivariate Splines Summer School, Summer 2009 Outline 1 Why multivariate splines for PDEs? Motivation
More informationThe Plane Stress Problem
The Plane Stress Problem Martin Kronbichler Applied Scientific Computing (Tillämpad beräkningsvetenskap) February 2, 2010 Martin Kronbichler (TDB) The Plane Stress Problem February 2, 2010 1 / 24 Outline
More informationNonparametric density estimation for elliptic problems with random perturbations
Nonparametric density estimation for elliptic problems with random perturbations, DqF Workshop, Stockholm, Sweden, 28--2 p. /2 Nonparametric density estimation for elliptic problems with random perturbations
More informationTime-dependent Dirichlet Boundary Conditions in Finite Element Discretizations
Time-dependent Dirichlet Boundary Conditions in Finite Element Discretizations Peter Benner and Jan Heiland November 5, 2015 Seminar Talk at Uni Konstanz Introduction Motivation A controlled physical processes
More informationarxiv: v1 [math.na] 9 Oct 2018
PRIMAL-DUAL REDUCED BASIS METHODS FOR CONVEX MINIMIZATION VARIATIONAL PROBLEMS: ROBUST TRUE SOLUTION A POSTERIORI ERROR CERTIFICATION AND ADAPTIVE GREEDY ALGORITHMS arxiv:1810.04073v1 [math.na] 9 Oct 2018
More informationKINEMATICS OF CONTINUA
KINEMATICS OF CONTINUA Introduction Deformation of a continuum Configurations of a continuum Deformation mapping Descriptions of motion Material time derivative Velocity and acceleration Transformation
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationSpace time finite element methods in biomedical applications
Space time finite element methods in biomedical applications Olaf Steinbach Institut für Angewandte Mathematik, TU Graz http://www.applied.math.tugraz.at SFB Mathematical Optimization and Applications
More informationIterative Methods for Linear Systems
Iterative Methods for Linear Systems 1. Introduction: Direct solvers versus iterative solvers In many applications we have to solve a linear system Ax = b with A R n n and b R n given. If n is large the
More informationFrom Completing the Squares and Orthogonal Projection to Finite Element Methods
From Completing the Squares and Orthogonal Projection to Finite Element Methods Mo MU Background In scientific computing, it is important to start with an appropriate model in order to design effective
More informationAdaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation
Adaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA www.math.umd.edu/ rhn 7th
More informationLagrange Relaxation and Duality
Lagrange Relaxation and Duality As we have already known, constrained optimization problems are harder to solve than unconstrained problems. By relaxation we can solve a more difficult problem by a simpler
More informationLocking phenomena in Computational Mechanics: nearly incompressible materials and plate problems
Locking phenomena in Computational Mechanics: nearly incompressible materials and plate problems C. Lovadina Dipartimento di Matematica Univ. di Pavia IMATI-CNR, Pavia Bologna, September, the 18th 2006
More informationSolving Dual Problems
Lecture 20 Solving Dual Problems We consider a constrained problem where, in addition to the constraint set X, there are also inequality and linear equality constraints. Specifically the minimization problem
More informationApplied/Numerical Analysis Qualifying Exam
Applied/Numerical Analysis Qualifying Exam August 9, 212 Cover Sheet Applied Analysis Part Policy on misprints: The qualifying exam committee tries to proofread exams as carefully as possible. Nevertheless,
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Sobolev Embedding Theorems Embedding Operators and the Sobolev Embedding Theorem
More informationA posteriori error control for the binary Mumford Shah model
A posteriori error control for the binary Mumford Shah model Benjamin Berkels 1, Alexander Effland 2, Martin Rumpf 2 1 AICES Graduate School, RWTH Aachen University, Germany 2 Institute for Numerical Simulation,
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationTHE STOKES SYSTEM R.E. SHOWALTER
THE STOKES SYSTEM R.E. SHOWALTER Contents 1. Stokes System 1 Stokes System 2 2. The Weak Solution of the Stokes System 3 3. The Strong Solution 4 4. The Normal Trace 6 5. The Mixed Problem 7 6. The Navier-Stokes
More informationA Posteriori Error Estimation Techniques for Finite Element Methods. Zhiqiang Cai Purdue University
A Posteriori Error Estimation Techniques for Finite Element Methods Zhiqiang Cai Purdue University Department of Mathematics, Purdue University Slide 1, March 16, 2017 Books Ainsworth & Oden, A posteriori
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 informationThe Mortar Boundary Element Method
The Mortar Boundary Element Method A Thesis submitted for the degree of Doctor of Philosophy by Martin Healey School of Information Systems, Computing and Mathematics Brunel University March 2010 Abstract
More informationSpace time finite and boundary element methods
Space time finite and boundary element methods Olaf Steinbach Institut für Numerische Mathematik, TU Graz http://www.numerik.math.tu-graz.ac.at based on joint work with M. Neumüller, H. Yang, M. Fleischhacker,
More informationREGULAR LAGRANGE MULTIPLIERS FOR CONTROL PROBLEMS WITH MIXED POINTWISE CONTROL-STATE CONSTRAINTS
REGULAR LAGRANGE MULTIPLIERS FOR CONTROL PROBLEMS WITH MIXED POINTWISE CONTROL-STATE CONSTRAINTS fredi tröltzsch 1 Abstract. A class of quadratic optimization problems in Hilbert spaces is considered,
More informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
More information14. Duality. ˆ Upper and lower bounds. ˆ General duality. ˆ Constraint qualifications. ˆ Counterexample. ˆ Complementary slackness.
CS/ECE/ISyE 524 Introduction to Optimization Spring 2016 17 14. Duality ˆ Upper and lower bounds ˆ General duality ˆ Constraint qualifications ˆ Counterexample ˆ Complementary slackness ˆ Examples ˆ Sensitivity
More informationDuality. Geoff Gordon & Ryan Tibshirani Optimization /
Duality Geoff Gordon & Ryan Tibshirani Optimization 10-725 / 36-725 1 Duality in linear programs Suppose we want to find lower bound on the optimal value in our convex problem, B min x C f(x) E.g., consider
More informationquantitative information on the error caused by using the solution of the linear problem to describe the response of the elastic material on a corner
Quantitative Justication of Linearization in Nonlinear Hencky Material Problems 1 Weimin Han and Hong-ci Huang 3 Abstract. The classical linear elasticity theory is based on the assumption that the size
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
More informationLagrange Duality. Daniel P. Palomar. Hong Kong University of Science and Technology (HKUST)
Lagrange Duality Daniel P. Palomar Hong Kong University of Science and Technology (HKUST) ELEC5470 - Convex Optimization Fall 2017-18, HKUST, Hong Kong Outline of Lecture Lagrangian Dual function Dual
More informationHybridized Discontinuous Galerkin Methods
Hybridized Discontinuous Galerkin Methods Theory and Christian Waluga Joint work with Herbert Egger (Uni Graz) 1st DUNE User Meeting, Stuttgart Christian Waluga (AICES) HDG Methods October 6-8, 2010 1
More informationSparse Optimization Lecture: Dual Certificate in l 1 Minimization
Sparse Optimization Lecture: Dual Certificate in l 1 Minimization Instructor: Wotao Yin July 2013 Note scriber: Zheng Sun Those who complete this lecture will know what is a dual certificate for l 1 minimization
More informationConvex Optimization M2
Convex Optimization M2 Lecture 3 A. d Aspremont. Convex Optimization M2. 1/49 Duality A. d Aspremont. Convex Optimization M2. 2/49 DMs DM par email: dm.daspremont@gmail.com A. d Aspremont. Convex Optimization
More informationA posteriori error estimation in the FEM
A posteriori error estimation in the FEM Plan 1. Introduction 2. Goal-oriented error estimates 3. Residual error estimates 3.1 Explicit 3.2 Subdomain error estimate 3.3 Self-equilibrated residuals 3.4
More informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg
More informationLecture: Duality.
Lecture: Duality http://bicmr.pku.edu.cn/~wenzw/opt-2016-fall.html Acknowledgement: this slides is based on Prof. Lieven Vandenberghe s lecture notes Introduction 2/35 Lagrange dual problem weak and strong
More informationLocal pointwise a posteriori gradient error bounds for the Stokes equations. Stig Larsson. Heraklion, September 19, 2011 Joint work with A.
Local pointwise a posteriori gradient error bounds for the Stokes equations Stig Larsson Department of Mathematical Sciences Chalmers University of Technology and University of Gothenburg Heraklion, September
More informationA High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations
Motivation Numerical methods Numerical tests Conclusions A High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations Xiaofeng Cai Department of Mathematics
More informationSpace-time sparse discretization of linear parabolic equations
Space-time sparse discretization of linear parabolic equations Roman Andreev August 2, 200 Seminar for Applied Mathematics, ETH Zürich, Switzerland Support by SNF Grant No. PDFMP2-27034/ Part of PhD thesis
More informationConvex Optimization Boyd & Vandenberghe. 5. Duality
5. Duality Convex Optimization Boyd & Vandenberghe Lagrange dual problem weak and strong duality geometric interpretation optimality conditions perturbation and sensitivity analysis examples generalized
More information6-1 The Positivstellensatz P. Parrilo and S. Lall, ECC
6-1 The Positivstellensatz P. Parrilo and S. Lall, ECC 2003 2003.09.02.10 6. The Positivstellensatz Basic semialgebraic sets Semialgebraic sets Tarski-Seidenberg and quantifier elimination Feasibility
More informationBilinear Quadrilateral (Q4): CQUAD4 in GENESIS
Bilinear Quadrilateral (Q4): CQUAD4 in GENESIS The Q4 element has four nodes and eight nodal dof. The shape can be any quadrilateral; we ll concentrate on a rectangle now. The displacement field in terms
More informationSwitching, sparse and averaged control
Switching, sparse and averaged control Enrique Zuazua Ikerbasque & BCAM Basque Center for Applied Mathematics Bilbao - Basque Country- Spain zuazua@bcamath.org http://www.bcamath.org/zuazua/ WG-BCAM, February
More informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More information1.The anisotropic plate model
Pré-Publicações do Departamento de Matemática Universidade de Coimbra Preprint Number 6 OPTIMAL CONTROL OF PIEZOELECTRIC ANISOTROPIC PLATES ISABEL NARRA FIGUEIREDO AND GEORG STADLER ABSTRACT: This paper
More informationRegularity Theory a Fourth Order PDE with Delta Right Hand Side
Regularity Theory a Fourth Order PDE with Delta Right Hand Side Graham Hobbs Applied PDEs Seminar, 29th October 2013 Contents Problem and Weak Formulation Example - The Biharmonic Problem Regularity Theory
More informationEE/AA 578, Univ of Washington, Fall Duality
7. Duality EE/AA 578, Univ of Washington, Fall 2016 Lagrange dual problem weak and strong duality geometric interpretation optimality conditions perturbation and sensitivity analysis examples generalized
More informationLECTURE 3: DISCRETE GRADIENT FLOWS
LECTURE 3: DISCRETE GRADIENT FLOWS Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA Tutorial: Numerical Methods for FBPs Free Boundary Problems and
More informationParallel Discontinuous Galerkin Method
Parallel Discontinuous Galerkin Method Yin Ki, NG The Chinese University of Hong Kong Aug 5, 2015 Mentors: Dr. Ohannes Karakashian, Dr. Kwai Wong Overview Project Goal Implement parallelization on Discontinuous
More informationFormulation of the displacement-based Finite Element Method and General Convergence Results
Formulation of the displacement-based Finite Element Method and General Convergence Results z Basics of Elasticity Theory strain e: measure of relative distortions u r r' y for small displacements : x
More informationAM 205: lecture 19. Last time: Conditions for optimality Today: Newton s method for optimization, survey of optimization methods
AM 205: lecture 19 Last time: Conditions for optimality Today: Newton s method for optimization, survey of optimization methods Optimality Conditions: Equality Constrained Case As another example of equality
More information5. Duality. Lagrangian
5. Duality Convex Optimization Boyd & Vandenberghe Lagrange dual problem weak and strong duality geometric interpretation optimality conditions perturbation and sensitivity analysis examples generalized
More informationOn second order sufficient optimality conditions for quasilinear elliptic boundary control problems
On second order sufficient optimality conditions for quasilinear elliptic boundary control problems Vili Dhamo Technische Universität Berlin Joint work with Eduardo Casas Workshop on PDE Constrained Optimization
More informationAlexei F. Cheviakov. University of Saskatchewan, Saskatoon, Canada. INPL seminar June 09, 2011
Direct Method of Construction of Conservation Laws for Nonlinear Differential Equations, its Relation with Noether s Theorem, Applications, and Symbolic Software Alexei F. Cheviakov University of Saskatchewan,
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More information1. Let a(x) > 0, and assume that u and u h are the solutions of the Dirichlet problem:
Mathematics Chalmers & GU TMA37/MMG800: Partial Differential Equations, 011 08 4; kl 8.30-13.30. Telephone: Ida Säfström: 0703-088304 Calculators, formula notes and other subject related material are not
More informationLagrange duality. The Lagrangian. We consider an optimization program of the form
Lagrange duality Another way to arrive at the KKT conditions, and one which gives us some insight on solving constrained optimization problems, is through the Lagrange dual. The dual is a maximization
More informationOn the local stability of semidefinite relaxations
On the local stability of semidefinite relaxations Diego Cifuentes Department of Mathematics Massachusetts Institute of Technology Joint work with Sameer Agarwal (Google), Pablo Parrilo (MIT), Rekha Thomas
More information