Risk Averse Shape Optimization
|
|
- Paula Warren
- 5 years ago
- Views:
Transcription
1 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 on PDE Constrained Optimization of Certain and Uncertain Processes 2009 Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 1 / 28
2 Outline 1 Introduction 2 Level set method Level set formulation Shape gradient Topological Derivative Optimization Algorithm 3 Risk Averse Functionals Uncertainty Expected Excess Excess Probability 4 Numerical Results Examples Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 2 / 28
3 Conceptual sketch D Γ0 O ΓD ΓN Figure: General setting in 2D Optimization Task min J(O) O ad = {O D : O Lipschitz-continuous } O O ad Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 3 / 28
4 Linear elasticity model The displacement u is given by the equation system PDE div(ae(u)) = f in O, u = 0 on Γ D, (Ae(u))n = g on Γ N Elastic body O R 3 O = Γ N Γ D, Γ D Volume forces f in O Neumann forces g on Γ N where e(u) = 1 2 ( u + ut ) is the linearized strain tensor and Hooke s law Aξ = 2µξ + λ(trξ)id, for any symmetric matrix ξ Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 4 / 28
5 Shape optimization problem Compliance J(O) = O f u dx + g u ds Γ N Least square error compared to target displacement ( J(O) = u u 0 2 dx O ) 1 2 Shape optimization problem min J(O) + lv(o) with l R, l > 0 O O ad O ad = {O D : O Lipschitz-continuous } Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 5 / 28
6 Outline 1 Introduction 2 Level set method Level set formulation Shape gradient Topological Derivative Optimization Algorithm 3 Risk Averse Functionals Uncertainty Expected Excess Excess Probability 4 Numerical Results Examples Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 6 / 28
7 Level set formulation Implicit description of the domain O via a level set function φ φ(x) = 0 <=> x O φ(x) < 0 <=> x O φ(x) > 0 <=> x O Figure: Levelset description in 2D Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 7 / 28
8 Outline 1 Introduction 2 Level set method Level set formulation Shape gradient Topological Derivative Optimization Algorithm 3 Risk Averse Functionals Uncertainty Expected Excess Excess Probability 4 Numerical Results Examples Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 8 / 28
9 Shape gradient We consider variations O t = (Id + t V)(O), t > 0 of a smooth elastic domain O for a smooth vector field V defined on the working domain D. The shape derivative of J(O) at O in direction V is defined as the Fréchet derivative of the mapping t J(O t ), i.e. J J(O t ) = J(O) + O, V + o( V ) Ot = (Id + t V) (O) Ot O cf. [Sokolowski, Zolesio 92], [Delfour, Zolesio ] Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 9 / 28
10 Shape gradient As a classical result of the shape sensitivity analysis the shape derivative takes the form < J ( [ ] ) (g u) O, V > = 2 + hg u + f u Aɛ(u) : ɛ(u) V n dν n Γ N + (Aɛ(u) : ɛ(u)) V n dν Γ D Here h denotes the mean curvature of O and n the outer normal. Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 10 / 28
11 Shape gradient in level set formulation When the domain O is implicitly deformed by varying the level set function φ φ t = φ + tψ the level set equation t φ + φ v n = 0 n = φ φ allows to define < J J, ψ >:=< φ O, ψ n φ > cf. [Osher, Sethian 88], [Burger, Osher 04] Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 11 / 28
12 Shape gradient in level set formulation We take into account a regularized gradient descent, based on the metric G(θ, ζ) = which is related to a Gaussian filter with width σ. The shape gradient is the solution of equation D θζ + σ2 θ ζ dx 2 G(grad φ J, θ) = < J, θ > θ H1,2 0 O (D) Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 12 / 28
13 Outline 1 Introduction 2 Level set method Level set formulation Shape gradient Topological Derivative Optimization Algorithm 3 Risk Averse Functionals Uncertainty Expected Excess Excess Probability 4 Numerical Results Examples Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 13 / 28
14 Topological Derivative Asymptotic behavier for infinitesimal small hole T (x) = lim ρ 0 J(O \ B ρ (x)) J(O) B ρ (x) Toppological derivative for the compliance D topo J(x) = π (λ + 2µ) 2µ (λ + µ) {4µAe(u i) : e(u i ) + (λ µ)tr Ae(u i )tr e(u i )} Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 14 / 28
15 Outline 1 Introduction 2 Level set method Level set formulation Shape gradient Topological Derivative Optimization Algorithm 3 Risk Averse Functionals Uncertainty Expected Excess Excess Probability 4 Numerical Results Examples Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 15 / 28
16 Optimization algorithm time continuous regularized gradient descent: with time discrete relaxation : φ(t) = grad φ J(φ) G(φ k+1 φ k, θ) = τ < J, θ > θ H1,2 0 O (D) additional ingrediens of the algorithm : multigrid method for the primal and the dual problem (d = 3) preconditioned CG (d = 2) cascadic optimization (from coarse to fine grid resolution) morphological smoothing when switching the grid resolution (σ = 2.5h or 4.5h) topological changes are performed every 10 steps Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 16 / 28
17 Outline 1 Introduction 2 Level set method Level set formulation Shape gradient Topological Derivative Optimization Algorithm 3 Risk Averse Functionals Uncertainty Expected Excess Excess Probability 4 Numerical Results Examples Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 17 / 28
18 Structure of Random Forces Assume that ω follows a discrete distribution with scenarios ω σ and probabilities π σ with S σ=1 π σ = 1 and basis loads (f k, g m ) spanning the load space: by linearity : f (ω) = K α k f k, g(ω) = k=1 ū(o, ω) = K α k u k f + k=1 M β m f m m=1 M β m u m g m=1 solves A(O, ū(o, ω σ ), ϕ) = l(o, ϕ, ω σ ) Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 18 / 28
19 Outline 1 Introduction 2 Level set method Level set formulation Shape gradient Topological Derivative Optimization Algorithm 3 Risk Averse Functionals Uncertainty Expected Excess Excess Probability 4 Numerical Results Examples Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 19 / 28
20 Expected Excess Considering the following problem { S } min π σ max{j(o, ω σ ) η, 0} O U ad σ=1 and the smooth approximation of the maximun function a2 + a a2 + ɛ + a max{a, 0} = =: Max(a) ɛ > we get the differentiable Expected Excess functional. expected excess S min{ee(o) := π σ Max(J(O, ω σ )) : O U} σ=1 Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 20 / 28
21 Outline 1 Introduction 2 Level set method Level set formulation Shape gradient Topological Derivative Optimization Algorithm 3 Risk Averse Functionals Uncertainty Expected Excess Excess Probability 4 Numerical Results Examples Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 21 / 28
22 Excess Probability Now consider min{p(j(o, ω) > η) : O U ad } The random stage ω follows a discrete distribution with finitely many scenarios ω σ and probabilities π σ according to S σ=1 π σ = 1 and we get min{p(j(o, ω σ ) > η) = S π σ H(J(O, ω σ ) η) where H(x) is supposed to be the Heaviside function. Again we use a smooth approximation H(x) tanh(kx) = 1 1+e 2k excess probability min{ep(o) := σ=1 S π σ H(J(O, ω σ )) : O U} σ=1 Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 22 / 28
23 Outline 1 Introduction 2 Level set method Level set formulation Shape gradient Topological Derivative Optimization Algorithm 3 Risk Averse Functionals Uncertainty Expected Excess Excess Probability 4 Numerical Results Examples Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 23 / 28
24 Test Setting % 95% Figure: Nonsymmetric stochastic loading Figure: The results for EV, EE and EP; the threshold η is set to 0.4 Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 24 / 28
25 Expected Exccess Figure: A sequence of results for the optimization with respect to the expected excess for η = 0.1, 0.2, 0.3, 0.4, 0.6, 0.8, 1.0., 1.5. Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 25 / 28
26 Cantilever load 1 load expected value expected value Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 26 / 28
27 Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 27 / 28
28 Thank you! Martin Pach (University Duisburg-Essen) Risk Averse Shape Optimization Trier09 28 / 28
A MULTISCALE APPROACH IN TOPOLOGY OPTIMIZATION
1 A MULTISCALE APPROACH IN TOPOLOGY OPTIMIZATION Grégoire ALLAIRE CMAP, Ecole Polytechnique The most recent results were obtained in collaboration with F. de Gournay, F. Jouve, O. Pantz, A.-M. Toader.
More informationSTRUCTURAL OPTIMIZATION BY THE LEVEL SET METHOD
Shape optimization 1 G. Allaire STRUCTURAL OPTIMIZATION BY THE LEVEL SET METHOD G. ALLAIRE CMAP, Ecole Polytechnique with F. JOUVE and A.-M. TOADER 1. Introduction 2. Setting of the problem 3. Shape differentiation
More informationShape Optimization under Uncertainty from a Stochastic Programming Point of View
Shape ptimization under Uncertainty from a Stochastic Programming Point of View Vom Fachbereich Mathematik der Universität Duisburg-Essen (Campus Duisburg) zur Erlangung des akademischen Grades eines Dr.
More informationRisk-averse structural topology optimization under random fields using stochastic expansion methods
Accepted Manuscript Risk-averse structural topology optimization under random fields using stochastic expansion methods Jesús Martínez-Frutos, David Herrero-Pérez, Mathieu Kessler, Francisco Periago PII:
More informationRisk Averse Elastic Shape Optimization with Parametrized Fine Scale Geometry
Risk Averse Elastic Shape Optimization with Parametrized Fine Scale Geometry Benedict Geihe Martin Lenz Martin Rumpf Ru diger Schultz March 27, 2012 Abstract Shape optimization of the fine scale geometry
More informationOptimal thickness of a cylindrical shell under dynamical loading
Optimal thickness of a cylindrical shell under dynamical loading Paul Ziemann Institute of Mathematics and Computer Science, E.-M.-A. University Greifswald, Germany e-mail paul.ziemann@uni-greifswald.de
More informationModelling of interfaces and free boundaries
University of Regensburg Regensburg, March 2009 Outline 1 Introduction 2 Obstacle problems 3 Stefan problem 4 Shape optimization Introduction What is a free boundary problem? Solve a partial differential
More informationShape optimisation with the level set method for contact problems in linearised elasticity
Shape optimisation with the level set method for contact problems in linearised elasticity Aymeric Maury, Grégoire Allaire, François Jouve To cite this version: Aymeric Maury, Grégoire Allaire, François
More informationarxiv: v1 [math.oc] 9 Dec 2013
arxiv:1312.2356v1 [math.oc] 9 Dec 2013 Multi-material phase field approach to structural topology optimization Luise Blank, M. Hassan Farshbaf-Shaker, Harald Garcke, Christoph Rupprecht, Vanessa Styles
More informationLevel Set Solution of an Inverse Electromagnetic Casting Problem using Topological Analysis
Level Set Solution of an Inverse Electromagnetic Casting Problem using Topological Analysis A. Canelas 1, A.A. Novotny 2 and J.R.Roche 3 1 Instituto de Estructuras y Transporte, Facultad de Ingeniería,
More informationA LINEARIZED APPROACH TO WORST-CASE DESIGN IN SHAPE OPTIMIZATION
1 A LINEARIZED APPROACH TO WORST-CASE DESIGN IN SHAPE OPTIMIZATION Grégoire ALLAIRE CMAP, Ecole Polytechnique Results obtained in collaboration with Ch. Dapogny (Rutgers University, formerly LJLL UPMC
More informationNumerical Methods of Applied Mathematics -- II Spring 2009
MIT OpenCourseWare http://ocw.mit.edu 18.336 Numerical Methods of Applied Mathematics -- II Spring 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationDetermination of thin elastic inclusions from boundary measurements.
Determination of thin elastic inclusions from boundary measurements. Elena Beretta in collaboration with E. Francini, S. Vessella, E. Kim and J. Lee September 7, 2010 E. Beretta (Università di Roma La
More informationA Post-treatement of the homogenization method in shape optimization
8 th World Congress on Structural and Multidisciplinary Optimization June 1-5, 2009, Lisbon, Portugal A Post-treatement of the homogenization method in shape optimization Olivier Pantz, Karim Trabelsi
More informationInformation geometry of mirror descent
Information geometry of mirror descent Geometric Science of Information Anthea Monod Department of Statistical Science Duke University Information Initiative at Duke G. Raskutti (UW Madison) and S. Mukherjee
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 informationSingular Perturbations of Stochastic Control Problems with Unbounded Fast Variables
Singular Perturbations of Stochastic Control Problems with Unbounded Fast Variables Joao Meireles joint work with Martino Bardi and Guy Barles University of Padua, Italy Workshop "New Perspectives in Optimal
More informationA Line search Multigrid Method for Large-Scale Nonlinear Optimization
A Line search Multigrid Method for Large-Scale Nonlinear Optimization Zaiwen Wen Donald Goldfarb Department of Industrial Engineering and Operations Research Columbia University 2008 Siam Conference on
More informationAdaptive Finite Element Methods Lecture Notes Winter Term 2017/18. R. Verfürth. Fakultät für Mathematik, Ruhr-Universität Bochum
Adaptive Finite Element Methods Lecture Notes Winter Term 2017/18 R. Verfürth Fakultät für Mathematik, Ruhr-Universität Bochum Contents Chapter I. Introduction 7 I.1. Motivation 7 I.2. Sobolev and finite
More informationOn the characterization of drilling rotation in the 6 parameter resultant shell theory
On the characterization of drilling rotation in the 6 parameter resultant shell theory Mircea Birsan and Patrizio Neff Chair for Nonlinear Analysis and Modelling Faculty of Mathematics, University Duisburg-Essen,
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 informationPH.D. PRELIMINARY EXAMINATION MATHEMATICS
UNIVERSITY OF CALIFORNIA, BERKELEY Dept. of Civil and Environmental Engineering FALL SEMESTER 2014 Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem
More informationCMS winter meeting 2008, Ottawa. The heat kernel on connected sums
CMS winter meeting 2008, Ottawa The heat kernel on connected sums Laurent Saloff-Coste (Cornell University) December 8 2008 p(t, x, y) = The heat kernel ) 1 x y 2 exp ( (4πt) n/2 4t The heat kernel is:
More informationNumerical methods in stochastic and two-scale shape optimization
Numerical methods in stochastic and two-scale shape optimization Dissertation zur Erlangung des Doktorgrades (Dr. rer. nat.) der Mathematisch-Naturwissenschaftlichen Fakultät der Rheinischen Friedrich-Wilhelms-Universität
More informationKramers formula for chemical reactions in the context of Wasserstein gradient flows. Michael Herrmann. Mathematical Institute, University of Oxford
eport no. OxPDE-/8 Kramers formula for chemical reactions in the context of Wasserstein gradient flows by Michael Herrmann Mathematical Institute, University of Oxford & Barbara Niethammer Mathematical
More informationAdaptive evolution : a population approach Benoît Perthame
Adaptive evolution : a population approach Benoît Perthame Adaptive dynamic : selection principle d dt n(x, t) = n(x, t)r( x, ϱ(t) ), ϱ(t) = R d n(x, t)dx. given x, n(x) = ϱ δ(x x), R( x, ϱ) = 0, ϱ( x).
More informationA sharp diffuse interface tracking method for approximating evolving interfaces
A sharp diffuse interface tracking method for approximating evolving interfaces Vanessa Styles and Charlie Elliott University of Sussex Overview Introduction Phase field models Double well and double obstacle
More informationA few words about the MTW tensor
A few words about the Ma-Trudinger-Wang tensor Université Nice - Sophia Antipolis & Institut Universitaire de France Salah Baouendi Memorial Conference (Tunis, March 2014) The Ma-Trudinger-Wang tensor
More informationKey words. Shape Derivatives, Topological Derivatives, Level Set Methods, Shape Optimization, Structural Design.
INCORPORATING TOPOLOGICAL DERIVATIVES INTO SHAPE DERIVATIVES BASED LEVEL SET METHODS LIN HE, CHIU-YEN KAO, AND STANLEY OSHER Abstract. Shape derivatives and topological derivatives have been incorporated
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 informationRome - May 12th Université Paris-Diderot - Laboratoire Jacques-Louis Lions. Mean field games equations with quadratic
Université Paris-Diderot - Laboratoire Jacques-Louis Lions Rome - May 12th 2011 Hamiltonian MFG Hamiltonian on the domain [0, T ] Ω, Ω standing for (0, 1) d : (HJB) (K) t u + σ2 2 u + 1 2 u 2 = f (x, m)
More informationStability of optimization problems with stochastic dominance constraints
Stability of optimization problems with stochastic dominance constraints D. Dentcheva and W. Römisch Stevens Institute of Technology, Hoboken Humboldt-University Berlin www.math.hu-berlin.de/~romisch SIAM
More informationHeaviside projection based aggregation in stress constrained topology optimization
Heaviside projection based aggregation in stress constrained topology optimization Cunfu Wang Xiaoping Qian Department of Mechanical Engineering, University of Wisconsin-Madison, 1513 University Avenue,
More informationScalable Algorithms for Optimal Control of Systems Governed by PDEs Under Uncertainty
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
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 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 informationAllen Cahn Equation in Two Spatial Dimension
Allen Cahn Equation in Two Spatial Dimension Yoichiro Mori April 25, 216 Consider the Allen Cahn equation in two spatial dimension: ɛ u = ɛ2 u + fu) 1) where ɛ > is a small parameter and fu) is of cubic
More informationMulti-phase structural optimization of multi-layered composites and viscoelastic materials
Multi-phase structural optimization of multi-layered composites and viscoelastic materials G Delgado, G. Allaire, M Hamdaoui To cite this version: G Delgado, G. Allaire, M Hamdaoui. Multi-phase structural
More informationPROBLEM OF CRACK UNDER QUASIBRITTLE FRACTURE V.A. KOVTUNENKO
PROBLEM OF CRACK UNDER QUASIBRITTLE FRACTURE V.A. KOVTUNENKO Overview: 1. Motivation 1.1. Evolutionary problem of crack propagation 1.2. Stationary problem of crack equilibrium 1.3. Interaction (contact+cohesion)
More informationComparison of Models for Finite Plasticity
Comparison of Models for Finite Plasticity A numerical study Patrizio Neff and Christian Wieners California Institute of Technology (Universität Darmstadt) Universität Augsburg (Universität Heidelberg)
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 informationStacking sequence and shape optimization of laminated composite plates via a level-set method
Stacking sequence and shape optimization of laminated composite plates via a level-set method G. Allaire 1, G. Delgado 2, March 26, 2015 1 CMAP, UMR CNRS 7641, Ecole Polytechnique, 91128 Palaiseau. 2 IRT
More informationR. M. Brown. 29 March 2008 / Regional AMS meeting in Baton Rouge. Department of Mathematics University of Kentucky. The mixed problem.
mixed R. M. Department of Mathematics University of Kentucky 29 March 2008 / Regional AMS meeting in Baton Rouge Outline mixed 1 mixed 2 3 4 mixed We consider the mixed boundary value Lu = 0 u = f D u
More informationInexact Newton Methods and Nonlinear Constrained Optimization
Inexact Newton Methods and Nonlinear Constrained Optimization Frank E. Curtis EPSRC Symposium Capstone Conference Warwick Mathematics Institute July 2, 2009 Outline PDE-Constrained Optimization Newton
More informationThe semi-geostrophic equations - a model for large-scale atmospheric flows
The semi-geostrophic equations - a model for large-scale atmospheric flows Beatrice Pelloni, University of Reading with M. Cullen (Met Office), D. Gilbert, T. Kuna INI - MFE Dec 2013 Introduction - Motivation
More informationStochastic Optimization with Inequality Constraints Using Simultaneous Perturbations and Penalty Functions
International Journal of Control Vol. 00, No. 00, January 2007, 1 10 Stochastic Optimization with Inequality Constraints Using Simultaneous Perturbations and Penalty Functions I-JENG WANG and JAMES C.
More informationDarwin and higher order approximations to Maxwell s equations in R 3. Sebastian Bauer Universität Duisburg-Essen
Darwin and higher order approximations to Maxwell s equations in R 3 Sebastian Bauer Universität Duisburg-Essen in close collaboration with the Maxwell group around Dirk Pauly Universität Duisburg-Essen
More informationComputing High Frequency Waves By the Level Set Method
Computing High Frequency Waves By the Level Set Method Hailiang Liu Department of Mathematics Iowa State University Collaborators: Li-Tien Cheng (UCSD), Stanley Osher (UCLA) Shi Jin (UW-Madison), Richard
More informationAn unconstrained multiphase thresholding approach for image segmentation
An unconstrained multiphase thresholding approach for image segmentation Benjamin Berkels Institut für Numerische Simulation, Rheinische Friedrich-Wilhelms-Universität Bonn, Nussallee 15, 53115 Bonn, Germany
More informationFoliations of hyperbolic space by constant mean curvature surfaces sharing ideal boundary
Foliations of hyperbolic space by constant mean curvature surfaces sharing ideal boundary David Chopp and John A. Velling December 1, 2003 Abstract Let γ be a Jordan curve in S 2, considered as the ideal
More informationTHE TOPOLOGICAL DESIGN OF MATERIALS WITH SPECIFIED THERMAL EXPANSION USING A LEVEL SET-BASED PARAMETERIZATION METHOD
11th. World Congress on Computational Mechanics (WCCM XI) 5th. European Conference on Computational Mechanics (ECCM V) 6th. European Conference on Computational Fluid Dynamics (ECFD VI) July 20-25, 2014,
More informationMATH 425, FINAL EXAM SOLUTIONS
MATH 425, FINAL EXAM SOLUTIONS Each exercise is worth 50 points. Exercise. a The operator L is defined on smooth functions of (x, y by: Is the operator L linear? Prove your answer. L (u := arctan(xy u
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 informationDesign of optimal RF pulses for NMR as a discrete-valued control problem
Design of optimal RF pulses for NMR as a discrete-valued control problem Christian Clason Faculty of Mathematics, Universität Duisburg-Essen joint work with Carla Tameling (Göttingen) and Benedikt Wirth
More informationA phase field model for the coupling between Navier-Stokes and e
A phase field model for the coupling between Navier-Stokes and electrokinetic equations Instituto de Matemáticas, CSIC Collaborators: C. Eck, G. Grün, F. Klingbeil (Erlangen Univertsität), O. Vantzos (Bonn)
More informationHolographic Geometries from Tensor Network States
Holographic Geometries from Tensor Network States J. Molina-Vilaplana 1 1 Universidad Politécnica de Cartagena Perspectives on Quantum Many-Body Entanglement, Mainz, Sep 2013 1 Introduction & Motivation
More informationMultilevel Preconditioning and Adaptive Sparse Solution of Inverse Problems
Multilevel and Adaptive Sparse of Inverse Problems Fachbereich Mathematik und Informatik Philipps Universität Marburg Workshop Sparsity and Computation, Bonn, 7. 11.6.2010 (joint work with M. Fornasier
More informationNumerical Approximation of Phase Field Models
Numerical Approximation of Phase Field Models Lecture 2: Allen Cahn and Cahn Hilliard Equations with Smooth Potentials Robert Nürnberg Department of Mathematics Imperial College London TUM Summer School
More informationAn Inexact Newton Method for Nonlinear Constrained Optimization
An Inexact Newton Method for Nonlinear Constrained Optimization Frank E. Curtis Numerical Analysis Seminar, January 23, 2009 Outline Motivation and background Algorithm development and theoretical results
More informationProper Orthogonal Decomposition in PDE-Constrained Optimization
Proper Orthogonal Decomposition in PDE-Constrained Optimization K. Kunisch Department of Mathematics and Computational Science University of Graz, Austria jointly with S. Volkwein Dynamic Programming Principle
More informationA Sobolev trust-region method for numerical solution of the Ginz
A Sobolev trust-region method for numerical solution of the Ginzburg-Landau equations Robert J. Renka Parimah Kazemi Department of Computer Science & Engineering University of North Texas June 6, 2012
More informationQUALIFYING EXAMINATION Harvard University Department of Mathematics Tuesday September 21, 2004 (Day 1)
QUALIFYING EXAMINATION Harvard University Department of Mathematics Tuesday September 21, 2004 (Day 1) Each of the six questions is worth 10 points. 1) Let H be a (real or complex) Hilbert space. We say
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationIntrinsic finite element modeling of a linear membrane shell problem
arxiv:3.39v [math.na] 5 Mar Intrinsic finite element modeling of a linear membrane shell problem Peter Hansbo Mats G. Larson Abstract A Galerkin finite element method for the membrane elasticity problem
More informationA Localized Linearized ROF Model for Surface Denoising
1 2 3 4 A Localized Linearized ROF Model for Surface Denoising Shingyu Leung August 7, 2008 5 Abstract 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 1 Introduction CT/MRI scan becomes a very
More informationStochastic Optimization Algorithms Beyond SG
Stochastic Optimization Algorithms Beyond SG Frank E. Curtis 1, Lehigh University involving joint work with Léon Bottou, Facebook AI Research Jorge Nocedal, Northwestern University Optimization Methods
More informationInformation geometry for bivariate distribution control
Information geometry for bivariate distribution control C.T.J.Dodson + Hong Wang Mathematics + Control Systems Centre, University of Manchester Institute of Science and Technology Optimal control of stochastic
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 informationNonlinear stabilization via a linear observability
via a linear observability Kaïs Ammari Department of Mathematics University of Monastir Joint work with Fathia Alabau-Boussouira Collocated feedback stabilization Outline 1 Introduction and main result
More informationAn optimal control problem for a parabolic PDE with control constraints
An optimal control problem for a parabolic PDE with control constraints PhD Summer School on Reduced Basis Methods, Ulm Martin Gubisch University of Konstanz October 7 Martin Gubisch (University of Konstanz)
More informationA Multigrid Method for Two Dimensional Maxwell Interface Problems
A Multigrid Method for Two Dimensional Maxwell Interface Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University USA JSA 2013 Outline A
More informationStochastic Gradient Descent with Variance Reduction
Stochastic Gradient Descent with Variance Reduction Rie Johnson, Tong Zhang Presenter: Jiawen Yao March 17, 2015 Rie Johnson, Tong Zhang Presenter: JiawenStochastic Yao Gradient Descent with Variance Reduction
More informationStatistical Measures of Uncertainty in Inverse Problems
Statistical Measures of Uncertainty in Inverse Problems Workshop on Uncertainty in Inverse Problems Institute for Mathematics and Its Applications Minneapolis, MN 19-26 April 2002 P.B. Stark Department
More informationDecision principles derived from risk measures
Decision principles derived from risk measures Marc Goovaerts Marc.Goovaerts@econ.kuleuven.ac.be Katholieke Universiteit Leuven Decision principles derived from risk measures - Marc Goovaerts p. 1/17 Back
More informationAn Introduction to Variational Inequalities
An Introduction to Variational Inequalities Stefan Rosenberger Supervisor: Prof. DI. Dr. techn. Karl Kunisch Institut für Mathematik und wissenschaftliches Rechnen Universität Graz January 26, 2012 Stefan
More informationAccelerated Dual Gradient-Based Methods for Total Variation Image Denoising/Deblurring Problems (and other Inverse Problems)
Accelerated Dual Gradient-Based Methods for Total Variation Image Denoising/Deblurring Problems (and other Inverse Problems) Donghwan Kim and Jeffrey A. Fessler EECS Department, University of Michigan
More informationOn Fréchet algebras with the dominating norm property
On Fréchet algebras with the dominating norm property Tomasz Ciaś Faculty of Mathematics and Computer Science Adam Mickiewicz University in Poznań Poland Banach Algebras and Applications Oulu, July 3 11,
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 informationLECTURE 4: GEOMETRIC PROBLEMS
LECTURE 4: GEOMETRIC PROBLEMS Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA Tutorial: Numerical Methods for FBPs Free Boundary Problems and Related
More informationUNCERTAINTY FUNCTIONAL DIFFERENTIAL EQUATIONS FOR FINANCE
Surveys in Mathematics and its Applications ISSN 1842-6298 (electronic), 1843-7265 (print) Volume 5 (2010), 275 284 UNCERTAINTY FUNCTIONAL DIFFERENTIAL EQUATIONS FOR FINANCE Iuliana Carmen Bărbăcioru Abstract.
More informationTHE LONGITUDINAL KAM-COCYCLE OF A MAGNETIC FLOW
THE LONGITUDINAL KAM-COCYCLE OF A MAGNETIC FLOW GABRIEL P. PATERNAIN Abstract. Let M be a closed oriented surface of negative Gaussian curvature and let Ω be a non-exact 2-form. Let λ be a small positive
More informationStochastic Gradient Descent in Continuous Time
Stochastic Gradient Descent in Continuous Time Justin Sirignano University of Illinois at Urbana Champaign with Konstantinos Spiliopoulos (Boston University) 1 / 27 We consider a diffusion X t X = R m
More informationMachine Learning. Support Vector Machines. Fabio Vandin November 20, 2017
Machine Learning Support Vector Machines Fabio Vandin November 20, 2017 1 Classification and Margin Consider a classification problem with two classes: instance set X = R d label set Y = { 1, 1}. Training
More information6.1 Formulation of the basic equations of torsion of prismatic bars (St. Venant) Figure 6.1: Torsion of a prismatic bar
Module 6 Torsion Learning Objectives 6.1 Formulation of the basic equations of torsion of prismatic bars (St. Venant) Readings: Sadd 9.3, Timoshenko Chapter 11 e e 1 e 3 Figure 6.1: Torsion of a prismatic
More informationRigidity and Non-rigidity Results on the Sphere
Rigidity and Non-rigidity Results on the Sphere Fengbo Hang Xiaodong Wang Department of Mathematics Michigan State University Oct., 00 1 Introduction It is a simple consequence of the maximum principle
More informationEinstein-Hilbert action on Connes-Landi noncommutative manifolds
Einstein-Hilbert action on Connes-Landi noncommutative manifolds Yang Liu MPIM, Bonn Analysis, Noncommutative Geometry, Operator Algebras Workshop June 2017 Motivations and History Motivation: Explore
More informationNeighboring feasible trajectories in infinite dimension
Neighboring feasible trajectories in infinite dimension Marco Mazzola Université Pierre et Marie Curie (Paris 6) H. Frankowska and E. M. Marchini Control of State Constrained Dynamical Systems Padova,
More informationAn Inexact Newton Method for Optimization
New York University Brown Applied Mathematics Seminar, February 10, 2009 Brief biography New York State College of William and Mary (B.S.) Northwestern University (M.S. & Ph.D.) Courant Institute (Postdoc)
More informationNumerical Methods for Large-Scale Nonlinear Systems
Numerical Methods for Large-Scale Nonlinear Systems Handouts by Ronald H.W. Hoppe following the monograph P. Deuflhard Newton Methods for Nonlinear Problems Springer, Berlin-Heidelberg-New York, 2004 Num.
More informationProblem Set Number 01, MIT (Winter-Spring 2018)
Problem Set Number 01, 18.377 MIT (Winter-Spring 2018) Rodolfo R. Rosales (MIT, Math. Dept., room 2-337, Cambridge, MA 02139) February 28, 2018 Due Thursday, March 8, 2018. Turn it in (by 3PM) at the Math.
More informationConvex Optimization in Computer Vision:
Convex Optimization in Computer Vision: Segmentation and Multiview 3D Reconstruction Yiyong Feng and Daniel P. Palomar The Hong Kong University of Science and Technology (HKUST) ELEC 5470 - Convex Optimization
More informationShape Optimization for Brake Squeal
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Shape Optimization for Brake Squeal Kohei Shintani 1 and Hideyuki Azegami 1 1 Nagoya University,
More informationGradient Descent and Implementation Solving the Euler-Lagrange Equations in Practice
1 Lecture Notes, HCI, 4.1.211 Chapter 2 Gradient Descent and Implementation Solving the Euler-Lagrange Equations in Practice Bastian Goldlücke Computer Vision Group Technical University of Munich 2 Bastian
More informationTHE SURFACE DIFFUSION FLOW WITH ELASTICITY IN THE PLANE
THE SURFACE DIFFUSION FLOW WITH ELASTICITY IN THE PLANE NICOLA FUSCO, VESA JULIN, AND MASSIMILIANO MORINI Abstract. In this paper we prove short-time existence of a smooth solution in the plane to the
More informationFundamentals of Linear Elasticity
Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy
More informationA new approach for Kirchhoff-Love plates and shells
A new approach for Kirchhoff-Love plates and shells Walter Zulehner Institute of Computational Mathematics JKU Linz, Austria AANMPDE 10 October 2-6, 2017, Paleochora, Crete, Greece Walter Zulehner (JKU
More informationPhysician Performance Assessment / Spatial Inference of Pollutant Concentrations
Physician Performance Assessment / Spatial Inference of Pollutant Concentrations Dawn Woodard Operations Research & Information Engineering Cornell University Johns Hopkins Dept. of Biostatistics, April
More informationStochastic Elastic-Plastic Finite Element Method for Performance Risk Simulations
Stochastic Elastic-Plastic Finite Element Method for Performance Risk Simulations Boris Jeremić 1 Kallol Sett 2 1 University of California, Davis 2 University of Akron, Ohio ICASP Zürich, Switzerland August
More informationWELL-POSEDNESS FOR HYPERBOLIC PROBLEMS (0.2)
WELL-POSEDNESS FOR HYPERBOLIC PROBLEMS We will use the familiar Hilbert spaces H = L 2 (Ω) and V = H 1 (Ω). We consider the Cauchy problem (.1) c u = ( 2 t c )u = f L 2 ((, T ) Ω) on [, T ] Ω u() = u H
More informationWaveform inversion and time-reversal imaging in attenuative TI media
Waveform inversion and time-reversal imaging in attenuative TI media Tong Bai 1, Tieyuan Zhu 2, Ilya Tsvankin 1, Xinming Wu 3 1. Colorado School of Mines 2. Penn State University 3. University of Texas
More information