TOWARD THE SOLUTION OF KKT SYSTEMS IN AERODYNAMICAL SHAPE OPTIMISATION
|
|
- Marsha Hart
- 5 years ago
- Views:
Transcription
1 CIRM, Marseille, 7-10 juillet TOWARD THE SOLUTION OF KKT SYSTEMS IN AERODYNAMICAL SHAPE OPTIMISATION Alain Dervieux, François Courty, INRIA, Sophia-Antipolis alain.dervieux@sophia.inria.fr 1
2 CIRM, Marseille, 7-10 juillet OVERVIEW Automated Differentiation One Shot SQP Multilevel preconditioner Application to sonic boom 2
3 CIRM, Marseille, 7-10 juillet The minimisation problem Find x opt = (Y opt,u opt ) in IR N Y IR N u Shape such that J(x opt ) = min Ψ(x) = 0 J(x) Lift, Boom j (u) = ( ) J u (Y,u) < Ψ(Y,u) = 0 ( ) (State: Euler) ( ) Ψ J Y (Y,u) Π = Y (Y,u) ( ) Ψ u (Y,u) Π, 1 > = 0 (Adjoint state) (Stationarity) 3
4 CIRM, Marseille, 7-10 juillet I. Progress in Automated Differentiation TROPICS: An INRIA project (team+research program) focalised on AD and applications : - a new generation tool: TAPENADE - the analysis of its application to optimal control Maurizio Araya, François Courty, Laurent Hascoet, Valérie Pascual, Alain Dervieux 4
5 CIRM, Marseille, 7-10 juillet Progress in inverse mode AD f = f p f p 1 f 1 f (x)ẋ = (f p f p 1 f p 2 f 1 (x)). (f p 1 f p 2 f 1 (x))..... (f 2 f 1 (x)). (f 1(x))ẋ. f t (x).y = (f 1 t (x)). (f 2 t f 1 (x))..... (f p 1 t f p 2 f 1 (x)). (f p t f p 1 f p 2 f 1 (x)). y. 5
6 CIRM, Marseille, 7-10 juillet DO i=1,n... A <- B... ENDDO Case of a loop 6
7 CIRM, Marseille, 7-10 juillet DO i=1,n... A <- B SAVE A(i)=A... ENDDO DO i=n,1,-1... RESTORE A(i) B <- function(a(i),ā(i))... ENDDO Case of a loop, end d 7
8 CIRM, Marseille, 7-10 juillet Inverse differentiation of a fixed point algorithm do iter = 1, update state variable if(residual.lt.eps)stop loop enddo 8
9 CIRM, Marseille, 7-10 juillet Inverse differentiation of a fixed point algorithm do iter = 1, update ST AT E ST ORE ST AT E(iter) if(residual.lt.eps)stop loop enddo do enddo iter = last iter,1 READ ST AT E(iter) update ADJOINT ST AT E - exact but subject to discontinuities - subject to instabilities, if update is nonlinear - very large storage requirement. 9
10 CIRM, Marseille, 7-10 juillet Inverse differentiation of a fixed point algorithm - Keep the differentiated algorithm, but replace the inverse-order read data by the direct-order ones : piggy back method of Griewank-Faure, - Keep the differentiated algorithm, but replace the inverse-order read data by the converged state: Hovland-Mohammadi-Bishof, etc. - Take the converged state, apply an adapted algorithm (preconditioned fixed point). Griewank-Faure, Piggyback Differentiation for Optimization, to be published in Proceedings of the Workshop on PDE Constrained Optimization in Santa FE, July 2001, Springer Lecture Notes on Computational Science and Engineering, Hovland-Mohammadi-Bishof, Automatic Differentiation of Navier-Stokes computations, Argonne National Laboratory MCS-P , 1997 Courty-Koobus-Dervieux-Hascoet 10
11 CIRM, Marseille, 7-10 juillet DO i=1,nelement... A <- B... ENDDO Inverse differentiation of an assembly loop An assembly loop can be executed in any order of the index. 11
12 CIRM, Marseille, 7-10 juillet Inverse differentiation of an assembly loop DO i=1,nnelement... A <- B SAVE A(i)=A... ENDDO Enormous memory storage DO i=nnelement,1,-1... RESTORE A(i) B <- function(a(i),ā(i))... ENDDO 12
13 CIRM, Marseille, 7-10 juillet DO i=1,n... A <- B SAVE A=A... Inverse differentiation by iterationwise adjoining... RESTORE A B <- function(a,ā)... ENDDO Hascoët, The Data-Dependence Graph of Adjoint Programs, INRIA-RR Courty-Koobus-Dervieux-Hascoet 13
14 CIRM, Marseille, 7-10 juillet Return to minimisation problem Find x opt = (Y opt,u opt ) in IR N Y IR N u Shape such that J(x opt ) = min Ψ(x) = 0 J(x) Lift, Bang j (u) = J u (Y,u) < ( Ψ u (Y,u) ( Ψ Y (Y,u) Ψ(Y,u) = 0 ) (State: Euler) Π = J Y (Y,u) (Adjoint state) ) Π, 1 > = 0 (Stationarity) 14
15 CIRM, Marseille, 7-10 juillet II. SQP / One shot optimisation algorithms A classical Sequential Quadratic Programming algorithm: For k=1,... In k : -Compute state update v (Newton). -Compute adjoint Π and gradient. -Compute Hessian 2 xl(x,π) (BFGS). -Compute step h for optimisation subproblem. Let d = v + h. Linear search in k : d ρd Trust region heuristics: Compare decrements of Merit function and quadratic model. Update Trust Region ball k. Next k. Byrd, Third SIAM Conference on Optimization, Houston, TX,1987, Omojokun, PhD Thesis, 1991 Byrd-Gilbert-Nocedal, Math. Programming, 2000 Nocedal-Wright, Numerical Optimization, Springer, 1999 N.B.: Still Three linearised state system at each iteration. 15
16 CIRM, Marseille, 7-10 juillet One Shot Methods S : a fixed (small) number of iterations of a given algorithm Iteration on state: Y k = Y k 1 SΨ( Ψ Iteration on adjoint Π k = Π k 1 S Π Y Iteration on control: u = u ρg(y k,u,π k ) ( J Y )) Ta asan, One Shot methods for optimal control of distributed parameter systems I: finite dimensional control, Report 91-2, Nasa report , ICASE, Marco-Beux, Multilevel optimization : application to shape optimum design with a One Shot method INRIA research report, no 2068,
17 CIRM, Marseille, 7-10 juillet Orientation In the seventies, gradient methods : exact solution of state and co-state, convergence control based on descent direction, Eighties: SQP, one Newton step for state, convergence control based on descent direction. Next step: One-Shot + descent direction 17
18 CIRM, Marseille, 7-10 juillet Assumption 1 : A few assumptions... Ψ(Y SΨ(Y,u),u) 0.99 Ψ(Y,u) Ψ (Π Sr Π ) J Y Y 0.99 Ψ Π J Y Y Assumption 2 : After k 0 applications of S, the gradient candidate, g(y k,u,π k ) is close enough to gradient for the satisfaction of Wolfe condition: k 0, k k 0, ρ k / 0 < ρ k ρ max, d k = M 1 g(y k,u,π k ), { J(Y k,u + ρ k d k ) J(Y k,u) 0.95ρ k < g(y k,u,π k ), d k > where Y k = Y (u + ρ k d k ) with k applications of S 18
19 CIRM, Marseille, 7-10 juillet SQP-like One-Shot algorithm Pour k=0,... *State partial restoration Y = Y + δy ( = Y SΨ Ψ ( )) J Adjoint partial restoration Π = Π + S Π Y Y Compute pseudo gradient g(y k,u,π k ) Linear search for descending under Wolfe conditions (restore state). If not successful, go to * (Optional) Trust Region heuristics Next k. 19
20 CIRM, Marseille, 7-10 juillet Application to a model problem Min Y (u) Y target 2, tel que Y = u on [0,1], Y (0) = Y (1) = 0 P 1 Finite Elements. Dimension of Y, u and adjoint Π is N = 500. Classical SQP algorithm One Shot SQP 20
21 CIRM, Marseille, 7-10 juillet Comparison of the efficiency of the two algorithms Algo Opt. iterations Conj. Grad. iterations Byrd-Omojokun, N= One shot, N= Byrd-Omojokun, N= One shot, N= Tab. 1 Efficiency of the two algorithms from the total number of conjugate gradient iterations for setting all residuals to F. Courty, A. Dervieux, A SQP-like one shot algorithm for Optimal Control problems, submitted to Mathematical Programming 21
22 CIRM, Marseille, 7-10 juillet Application to sonic boom reduction j(γ) = α 1 (C D C target D ) 2 + α 2 (C L C target L ) 2 + α 3 Ω B ( p) 2 dv α 1 = 1.0 α 2 = 10.0 α 3 = Fig. 1 Aircraft shape and gradient 22
23 CIRM, Marseille, 7-10 juillet Application to sonic boom reduction (cont d) Pressure in a plan under the aircraft (10 iterations) Pressure in two cuts, parallel to symmetry axis 23
24 CIRM, Marseille, 7-10 juillet Wing optimisation ONERA M6. initial and optimised wing profiles. ONERA M6. Pressure in a plane under wing. 24
25 CIRM, Marseille, 7-10 juillet Efficiency of the one-shot approach N Y = = N u = 780 one-shot SQP optimisation is equivalent to 10 complete state+adjoint solutions. 200 iterations: 14h 40 min 25
26 CIRM, Marseille, 7-10 juillet Efficiency of the one-shot approach: finer meshes N Y = = N u = iterations cost 6h CPU instead of 44h, - the one-shot optimisation CPU is equivalent to 15 state+adjoint - of course, with the finer mesh, convergence is slower. 26
27 CIRM, Marseille, 7-10 juillet III. Which preconditioners for Shape Optimal Design? u n+1 = u n ρg n 27
28 CIRM, Marseille, 7-10 juillet III. Which preconditioners for Shape Optimal Design? u n+1 = u n ρ C g n
29 CIRM, Marseille, 7-10 juillet III. Which preconditioners for Shape Optimal Design? u n+1 = u n ρ C g n Damp high frequencies: Reuther-Jameson, Mohammadi-Pironneau: - Laplace-Beltrami C = (Id µ ) 1 Identify the degree of the Fourier transform of Hessian: Arian-Ta asan - If degree = 2: Laplace-Beltrami C = (Id µ ) 1 - If degree = +1: C =Neumann-Dirichlet Pseudo-differential operator Reuther-Jameson, Aerodynamic Shape Optimization of Wing and Wing-Body Configurations Using Control Theory, AIAA Paper , 1995 Mohammadi-Pironneau, Applied shape optimization for fluids, Clarendon Press - Oxford, 2001 Arian-Ta asan, Analysis of the Hessian for aerodynamic optimization: Inviscid flow., Comput. Fluids, 28, 7, p ,
30 CIRM, Marseille, 7-10 juillet Shape Optimal Design for the Dirichlet problem Find γ, the optimal displacement of boundary in normal direction: γ = ArgMin J(Z); Z = 1 in Ω γ ; Z Ωγ = 0 Hadamard s variational formula: Unbounded operator in C k. γ γ ρ Z Π n n Regularity loss is 1. Palmerio-Dervieux, CRAS,A,280, , Murat-Simon, Sur le contrôle d un domaine géométrique. Publication 76015, Lab. d Analyse Numérique, Univ. Paris 6, Pironneau, Optimal shape design for elliptic systems, Springer-Verlag,
31 CIRM, Marseille, 7-10 juillet Hadamard formula for the Euler system (formal) j (γ 0,θ) = (F (W ) Π + G(W ) Π + H(W ) Π Ω γ0 x y z )( n V )θ d Ω γ0 + ( p Π + p Π)( n V )θ d Ω γ0 Ω γ0 + (Π 5 p + W 1 Π 1 + W 2 Π 2 + W 3 Π 3 + W 4 Π 4 + W 5 φ 5 ) q Ω γ0 γ δγ d Ω γ 0 - spatial derivatives of pressure p - spatial derivatives of boundary parameter γ Formally, one degree of smoothness is lost Beux-Dervieux, Exact-gradient shape optimization of a 2D Euler flow, Finite Elements in Analysis and Design, 12, p ,
32 CIRM, Marseille, 7-10 juillet Additive multilevel preconditioner (Au,v) = (f,v) v V, f given in V. V 1 V k V Let Q k be the projector: u V, v V k, (Q k u,v) = (u,v). Cg = + k ( ) k 1 2 a (Q k g Q k 1)g Property: C is continuous from H s 0(Ω) to H s+a 0 (Ω) (s, a in adhoc intervals) Bramble, Pasciak, and Xu, Parallel multilevel preconditioners, Math. Comput., 1990, 55, 191,p 1-22 Kunoth, PhD thesis, University of Berlin, Courty, Dervieux, Multilevel Functional Preconditioning for Shape Optimization, submitted to Numerische Mathematik. 31
33 CIRM, Marseille, 7-10 juillet Application to a model problem u opt = ArgMin 1 2 ( u 2 fu) dv
34 CIRM, Marseille, 7-10 juillet Application to a model problem u opt = ArgMin 1 2 ( u 2 fu) dv Verification of the optimal value of a, mesh independant convergence 32
35 CIRM, Marseille, 7-10 juillet Extension to unstructured : agglomeration 4 cells grouped in a single macro-cell Lallemand, Steve, Dervieux, Unstructured multiggridding by volume agglomeration : current status, Computer Fluids, 21, 3, , 1992 Marco, Dervieux, Multilevel parametrization for aerodynamical optimization of 3D shapes, Finite Elements in Analysis and Design, 26, p ,
36 CIRM, Marseille, 7-10 juillet Smooth transfer operators P : prolongation, from coarse to fine. P : restriction, from fine to coarse. L : A smoothing operator, L its transpose: (Lu) i = (1 θ)u i + θ Meas(j) u j j V(i) {i} j V(i) {i} Meas(j) where j belongs to the set V(i) of neighboring cells of i, Meas(j) is cell j measure. 34
37 CIRM, Marseille, 7-10 juillet Building the preconditioner in the unstructured case 1. Build the different levels by agglomeration, 2. Define operators L m,p m,pm,l m, 3. Define: Q k = L m P m PmL m 4.Compute: Cg = n k 1 m k ( ) k 1 2 a (Q k g Q k 1)g 35
38 CIRM, Marseille, 7-10 juillet Application to sonic boom (coarse) Optimal shape design by SQP/conjugate gradient: Gradient convergence Experimental evaluation of the optimal parameter a opt 36
39 CIRM, Marseille, 7-10 juillet Application to sonic boom (FINE) Optimal shape design by SQP/conjugate gradient: Gradient convergence Experimental evaluation of the optimal parameter a opt 37
40 CIRM, Marseille, 7-10 juillet To conclude... - The previous results are recent. It remains to combine One-Shot and Multi-Level. - The Multi-Level approach is being extended to non-symmetric problems. The whole optimality system would then be preconditioned by this mean. - Optimality systems can also be constructed for optimal meshing, and we shall try to solve them with the above methods. 38
MULTILEVEL FUNCTIONAL PRECONDITIONING FOR SHAPE OPTIMISATION. F. Courty, A. Dervieux
MULTILEVEL FUNCTIONAL PRECONDITIONING FOR SHAPE OPTIMISATION F. Courty, A. Dervieux INRIA, Centre de Sophia-Antipolis, BP 93, 2004 Route des Lucioles, 06902 Sophia-Antipolis Cedex, France Francois.Courty@inria.fr,
More informationAutomatic Differentiation for Optimum Design, Applied to Sonic Boom Reduction
Automatic Differentiation for Optimum Design, Applied to Sonic Boom Reduction Laurent Hascoët, Mariano Vázquez, Alain Dervieux Laurent.Hascoet@sophia.inria.fr Tropics Project, INRIA Sophia-Antipolis AD
More informationSAROD Conference, Hyderabad, december 2005 CONTINUOUS MESH ADAPTATION MODELS FOR CFD
1 SAROD Conference, Hyderabad, december 2005 CONTINUOUS MESH ADAPTATION MODELS FOR CFD Bruno Koobus, 1,2 Laurent Hascoët, 1 Frédéric Alauzet, 3 Adrien Loseille, 3 Youssef Mesri, 1 Alain Dervieux 1 1 INRIA
More informationSensitivity analysis by adjoint Automatic Differentiation and Application
Sensitivity analysis by adjoint Automatic Differentiation and Application Anca Belme Massimiliano Martinelli Laurent Hascoët Valérie Pascual Alain Dervieux INRIA Sophia-Antipolis Project TROPICS Alain.Dervieux@inria.fr
More informationCertification of Derivatives Computed by Automatic Differentiation
Certification of Derivatives Computed by Automatic Differentiation Mauricio Araya Polo & Laurent Hascoët Project TROPCS WSEAS, Cancún, México, May 13, 2005 1 Plan ntroduction (Background) Automatic Differentiation
More informationHigher-Order Methods
Higher-Order Methods Stephen J. Wright 1 2 Computer Sciences Department, University of Wisconsin-Madison. PCMI, July 2016 Stephen Wright (UW-Madison) Higher-Order Methods PCMI, July 2016 1 / 25 Smooth
More informationNumerical Methods for PDE-Constrained Optimization
Numerical Methods for PDE-Constrained Optimization Richard H. Byrd 1 Frank E. Curtis 2 Jorge Nocedal 2 1 University of Colorado at Boulder 2 Northwestern University Courant Institute of Mathematical Sciences,
More informationTrust-Region SQP Methods with Inexact Linear System Solves for Large-Scale Optimization
Trust-Region SQP Methods with Inexact Linear System Solves for Large-Scale Optimization Denis Ridzal Department of Computational and Applied Mathematics Rice University, Houston, Texas dridzal@caam.rice.edu
More informationMS&E 318 (CME 338) Large-Scale Numerical Optimization
Stanford University, Management Science & Engineering (and ICME) MS&E 318 (CME 338) Large-Scale Numerical Optimization 1 Origins Instructor: Michael Saunders Spring 2015 Notes 9: Augmented Lagrangian Methods
More informationLecture 15: SQP methods for equality constrained optimization
Lecture 15: SQP methods for equality constrained optimization Coralia Cartis, Mathematical Institute, University of Oxford C6.2/B2: Continuous Optimization Lecture 15: SQP methods for equality constrained
More informationOPER 627: Nonlinear Optimization Lecture 14: Mid-term Review
OPER 627: Nonlinear Optimization Lecture 14: Mid-term Review Department of Statistical Sciences and Operations Research Virginia Commonwealth University Oct 16, 2013 (Lecture 14) Nonlinear Optimization
More informationNumerical Optimization Algorithms
Numerical Optimization Algorithms 1. Overview. Calculus of Variations 3. Linearized Supersonic Flow 4. Steepest Descent 5. Smoothed Steepest Descent Overview 1 Two Main Categories of Optimization Algorithms
More informationIterative Solvers in the Finite Element Solution of Transient Heat Conduction
Iterative Solvers in the Finite Element Solution of Transient Heat Conduction Mile R. Vuji~i} PhD student Steve G.R. Brown Senior Lecturer Materials Research Centre School of Engineering University of
More informationAn Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems
An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems P.-O. Persson and J. Peraire Massachusetts Institute of Technology 2006 AIAA Aerospace Sciences Meeting, Reno, Nevada January 9,
More informationToward Practical Aerodynamic Design Through Numerical Optimization
Toward Practical Aerodynamic Design Through Numerical Optimization David W. Zingg, and Laura Billing, Institute for Aerospace Studies, University of Toronto 4925 Dufferin St., Toronto, Ontario M3H 5T6,
More informationEfficient Hessian Calculation using Automatic Differentiation
25th AIAA Applied Aerodynamics Conference, 25-28 June 2007, Miami Efficient Hessian Calculation using Automatic Differentiation Devendra P. Ghate and Michael B. Giles Oxford University Computing Laboratory,
More informationLecture 14: Newton s Method
10-725/36-725: Conve Optimization Fall 2016 Lecturer: Javier Pena Lecture 14: Newton s ethod Scribes: Varun Joshi, Xuan Li Note: LaTeX template courtesy of UC Berkeley EECS dept. Disclaimer: These notes
More informationREPORTS IN INFORMATICS
REPORTS IN INFORMATICS ISSN 0333-3590 A class of Methods Combining L-BFGS and Truncated Newton Lennart Frimannslund Trond Steihaug REPORT NO 319 April 2006 Department of Informatics UNIVERSITY OF BERGEN
More informationOne-Shot Domain Decomposition Methods for Shape Optimization Problems 1 Introduction Page 565
1 One-Shot Domain Decomposition Methods for Shape 2 Optimization Problems 3 Rongliang Chen 1 and Xiao-Chuan Cai 2 4 1 College of Mathematics and Econometrics, Hunan University, Changsha, Hunan 410082,
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More 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 informationOptimizing Runge-Kutta smoothers for unsteady flow problems
Optimizing Runge-Kutta smoothers for unsteady flow problems Philipp Birken 1 November 24, 2011 1 Institute of Mathematics, University of Kassel, Heinrich-Plett-Str. 40, D-34132 Kassel, Germany. email:
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 informationSuboptimal Open-loop Control Using POD. Stefan Volkwein
Institute for Mathematics and Scientific Computing University of Graz, Austria PhD program in Mathematics for Technology Catania, May 22, 2007 Motivation Optimal control of evolution problems: min J(y,
More informationIntegration of Sequential Quadratic Programming and Domain Decomposition Methods for Nonlinear Optimal Control Problems
Integration of Sequential Quadratic Programming and Domain Decomposition Methods for Nonlinear Optimal Control Problems Matthias Heinkenschloss 1 and Denis Ridzal 2 1 Department of Computational and Applied
More 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 informationNumerical Optimization Professor Horst Cerjak, Horst Bischof, Thomas Pock Mat Vis-Gra SS09
Numerical Optimization 1 Working Horse in Computer Vision Variational Methods Shape Analysis Machine Learning Markov Random Fields Geometry Common denominator: optimization problems 2 Overview of Methods
More informationNumerical Optimization of Partial Differential Equations
Numerical Optimization of Partial Differential Equations Part I: basic optimization concepts in R n Bartosz Protas Department of Mathematics & Statistics McMaster University, Hamilton, Ontario, Canada
More informationNumerical Optimization: Basic Concepts and Algorithms
May 27th 2015 Numerical Optimization: Basic Concepts and Algorithms R. Duvigneau R. Duvigneau - Numerical Optimization: Basic Concepts and Algorithms 1 Outline Some basic concepts in optimization Some
More informationMultidisciplinary System Design Optimization (MSDO)
Multidisciplinary System Design Optimization (MSDO) Numerical Optimization II Lecture 8 Karen Willcox 1 Massachusetts Institute of Technology - Prof. de Weck and Prof. Willcox Today s Topics Sequential
More informationAlgorithms for Constrained Optimization
1 / 42 Algorithms for Constrained Optimization ME598/494 Lecture Max Yi Ren Department of Mechanical Engineering, Arizona State University April 19, 2015 2 / 42 Outline 1. Convergence 2. Sequential quadratic
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 informationImplicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
More informationOn Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities
On Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities Heiko Berninger, Ralf Kornhuber, and Oliver Sander FU Berlin, FB Mathematik und Informatik (http://www.math.fu-berlin.de/rd/we-02/numerik/)
More informationMS&E 318 (CME 338) Large-Scale Numerical Optimization
Stanford University, Management Science & Engineering (and ICME) MS&E 318 (CME 338) Large-Scale Numerical Optimization Instructor: Michael Saunders Spring 2015 Notes 11: NPSOL and SNOPT SQP Methods 1 Overview
More informationPDE-Constrained and Nonsmooth Optimization
Frank E. Curtis October 1, 2009 Outline PDE-Constrained Optimization Introduction Newton s method Inexactness Results Summary and future work Nonsmooth Optimization Sequential quadratic programming (SQP)
More informationCONTINUOUS MESH ADAPTATION MODELS FOR CFD
Computational Fluid Dynamics JOURNAL CONTINUOUS MESH ADAPTATION MODELS FOR CFD Youssef Mesri, 1 Frédéric Alauzet, 3 Adrien Loseille, 3 Laurent Hascoët, 1 Bruno Koobus, 1,2 Alain Dervieux 1 Abstract Continuous
More informationConstrained optimization: direct methods (cont.)
Constrained optimization: direct methods (cont.) Jussi Hakanen Post-doctoral researcher jussi.hakanen@jyu.fi Direct methods Also known as methods of feasible directions Idea in a point x h, generate a
More informationNumerical Optimization of Partial Differential Equations
Numerical Optimization of Partial Differential Equations Bartosz Protas Department of Mathematics & Statistics McMaster University, Hamilton, Ontario, Canada URL: http://www.math.mcmaster.ca/bprotas Rencontres
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 informationCompressible Navier-Stokes (Euler) Solver based on Deal.II Library
Compressible Navier-Stokes (Euler) Solver based on Deal.II Library Lei Qiao Northwestern Polytechnical University Xi an, China Texas A&M University College Station, Texas Fifth deal.ii Users and Developers
More informationIntro to Automatic Differentiation
Intro to Automatic Differentiation Thomas Kaminski (http://fastopt.com) Thanks to: Simon Blessing (FastOpt), Ralf Giering (FastOpt), Nadine Gobron (JRC), Wolfgang Knorr (QUEST), Thomas Lavergne (Met.No),
More informationA MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY
A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,
More informationSF2822 Applied Nonlinear Optimization. Preparatory question. Lecture 9: Sequential quadratic programming. Anders Forsgren
SF2822 Applied Nonlinear Optimization Lecture 9: Sequential quadratic programming Anders Forsgren SF2822 Applied Nonlinear Optimization, KTH / 24 Lecture 9, 207/208 Preparatory question. Try to solve theory
More informationShort title: Total FETI. Corresponding author: Zdenek Dostal, VŠB-Technical University of Ostrava, 17 listopadu 15, CZ Ostrava, Czech Republic
Short title: Total FETI Corresponding author: Zdenek Dostal, VŠB-Technical University of Ostrava, 17 listopadu 15, CZ-70833 Ostrava, Czech Republic mail: zdenek.dostal@vsb.cz fax +420 596 919 597 phone
More informationImproving L-BFGS Initialization for Trust-Region Methods in Deep Learning
Improving L-BFGS Initialization for Trust-Region Methods in Deep Learning Jacob Rafati http://rafati.net jrafatiheravi@ucmerced.edu Ph.D. Candidate, Electrical Engineering and Computer Science University
More informationCopyright 2003 by the American Institute of Aeronautics and Astronautics Inc. All rights reserved
1 ISABE29-1225 COMPUTATION OF SECOND-ORDER SENSITIVITIES OF THE TOTAL PRESSURE LOSSES FUNCTIONAL IN CASCADE FLOWS T. Zervogiannis, D.I. Papadimitriou and K.C. Giannakoglou National Technical University
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 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 informationA multilevel, level-set method for optimizing eigenvalues in shape design problems
A multilevel, level-set method for optimizing eigenvalues in shape design problems E. Haber July 22, 2003 Abstract In this paper we consider optimal design problems that involve shape optimization. The
More informationOn Lagrange multipliers of trust-region subproblems
On Lagrange multipliers of trust-region subproblems Ladislav Lukšan, Ctirad Matonoha, Jan Vlček Institute of Computer Science AS CR, Prague Programy a algoritmy numerické matematiky 14 1.- 6. června 2008
More informationAn Investigation of the Attainable Efficiency of Flight at Mach One or Just Beyond
An Investigation of the Attainable Efficiency of Flight at Mach One or Just Beyond Antony Jameson Department of Aeronautics and Astronautics AIAA Aerospace Sciences Meeting, Reno, NV AIAA Paper 2007-0037
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 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 informationReduced-Hessian Methods for Constrained Optimization
Reduced-Hessian Methods for Constrained Optimization Philip E. Gill University of California, San Diego Joint work with: Michael Ferry & Elizabeth Wong 11th US & Mexico Workshop on Optimization and its
More informationNewton s Method. Javier Peña Convex Optimization /36-725
Newton s Method Javier Peña Convex Optimization 10-725/36-725 1 Last time: dual correspondences Given a function f : R n R, we define its conjugate f : R n R, f ( (y) = max y T x f(x) ) x Properties and
More informationModule III: Partial differential equations and optimization
Module III: Partial differential equations and optimization Martin Berggren Department of Information Technology Uppsala University Optimization for differential equations Content Martin Berggren (UU)
More informationIs My CFD Mesh Adequate? A Quantitative Answer
Is My CFD Mesh Adequate? A Quantitative Answer Krzysztof J. Fidkowski Gas Dynamics Research Colloqium Aerospace Engineering Department University of Michigan January 26, 2011 K.J. Fidkowski (UM) GDRC 2011
More informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
More informationChapter 6. Finite Element Method. Literature: (tiny selection from an enormous number of publications)
Chapter 6 Finite Element Method Literature: (tiny selection from an enormous number of publications) K.J. Bathe, Finite Element procedures, 2nd edition, Pearson 2014 (1043 pages, comprehensive). Available
More informationSeminal papers in nonlinear optimization
Seminal papers in nonlinear optimization Nick Gould, CSED, RAL, Chilton, OX11 0QX, England (n.gould@rl.ac.uk) December 7, 2006 The following papers are classics in the field. Although many of them cover
More informationSolving Distributed Optimal Control Problems for the Unsteady Burgers Equation in COMSOL Multiphysics
Excerpt from the Proceedings of the COMSOL Conference 2009 Milan Solving Distributed Optimal Control Problems for the Unsteady Burgers Equation in COMSOL Multiphysics Fikriye Yılmaz 1, Bülent Karasözen
More informationMultigrid and Domain Decomposition Methods for Electrostatics Problems
Multigrid and Domain Decomposition Methods for Electrostatics Problems Michael Holst and Faisal Saied Abstract. We consider multigrid and domain decomposition methods for the numerical solution of electrostatics
More informationAn Active Set Strategy for Solving Optimization Problems with up to 200,000,000 Nonlinear Constraints
An Active Set Strategy for Solving Optimization Problems with up to 200,000,000 Nonlinear Constraints Klaus Schittkowski Department of Computer Science, University of Bayreuth 95440 Bayreuth, Germany e-mail:
More informationDerivatives for Time-Spectral Computational Fluid Dynamics using an Automatic Differentiation Adjoint
Derivatives for Time-Spectral Computational Fluid Dynamics using an Automatic Differentiation Adjoint Charles A. Mader University of Toronto Institute for Aerospace Studies Toronto, Ontario, Canada Joaquim
More informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
More informationMultigrid and Iterative Strategies for Optimal Control Problems
Multigrid and Iterative Strategies for Optimal Control Problems John Pearson 1, Stefan Takacs 1 1 Mathematical Institute, 24 29 St. Giles, Oxford, OX1 3LB e-mail: john.pearson@worc.ox.ac.uk, takacs@maths.ox.ac.uk
More informationA STUDY OF MULTIGRID SMOOTHERS USED IN COMPRESSIBLE CFD BASED ON THE CONVECTION DIFFUSION EQUATION
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
More informationInfeasibility Detection and an Inexact Active-Set Method for Large-Scale Nonlinear Optimization
Infeasibility Detection and an Inexact Active-Set Method for Large-Scale Nonlinear Optimization Frank E. Curtis, Lehigh University involving joint work with James V. Burke, University of Washington Daniel
More informationA Trust-region-based Sequential Quadratic Programming Algorithm
Downloaded from orbit.dtu.dk on: Oct 19, 2018 A Trust-region-based Sequential Quadratic Programming Algorithm Henriksen, Lars Christian; Poulsen, Niels Kjølstad Publication date: 2010 Document Version
More informationMultigrid Methods for Elliptic Obstacle Problems on 2D Bisection Grids
Multigrid Methods for Elliptic Obstacle Problems on 2D Bisection Grids Long Chen 1, Ricardo H. Nochetto 2, and Chen-Song Zhang 3 1 Department of Mathematics, University of California at Irvine. chenlong@math.uci.edu
More informationMultigrid Based Optimal Shape and Topology Design in Magnetostatics
Multigrid Based Optimal Shape and Topology Design in Magnetostatics Dalibor Lukáš 1 Department of Applied Mathematics, VŠB Technical University of Ostrava, 17. listopadu 15, 708 33 Ostrava Poruba, Czech
More informationParameter Identification in Partial Differential Equations
Parameter Identification in Partial Differential Equations Differentiation of data Not strictly a parameter identification problem, but good motivation. Appears often as a subproblem. Given noisy observation
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 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 informationNumerical Analysis of Electromagnetic Fields
Pei-bai Zhou Numerical Analysis of Electromagnetic Fields With 157 Figures Springer-Verlag Berlin Heidelberg New York London Paris Tokyo Hong Kong Barcelona Budapest Contents Part 1 Universal Concepts
More informationPreconditioned conjugate gradient algorithms with column scaling
Proceedings of the 47th IEEE Conference on Decision and Control Cancun, Mexico, Dec. 9-11, 28 Preconditioned conjugate gradient algorithms with column scaling R. Pytla Institute of Automatic Control and
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
More informationIntroduction. New Nonsmooth Trust Region Method for Unconstraint Locally Lipschitz Optimization Problems
New Nonsmooth Trust Region Method for Unconstraint Locally Lipschitz Optimization Problems Z. Akbari 1, R. Yousefpour 2, M. R. Peyghami 3 1 Department of Mathematics, K.N. Toosi University of Technology,
More informationNewton s Method. Ryan Tibshirani Convex Optimization /36-725
Newton s Method Ryan Tibshirani Convex Optimization 10-725/36-725 1 Last time: dual correspondences Given a function f : R n R, we define its conjugate f : R n R, Properties and examples: f (y) = max x
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 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 informationOn Lagrange multipliers of trust region subproblems
On Lagrange multipliers of trust region subproblems Ladislav Lukšan, Ctirad Matonoha, Jan Vlček Institute of Computer Science AS CR, Prague Applied Linear Algebra April 28-30, 2008 Novi Sad, Serbia Outline
More informationApplication of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations
Application of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations Matthew McMullen and Antony Jameson and Juan J. Alonso Dept. of Aeronautics & Astronautics Stanford University
More informationA SHORT NOTE COMPARING MULTIGRID AND DOMAIN DECOMPOSITION FOR PROTEIN MODELING EQUATIONS
A SHORT NOTE COMPARING MULTIGRID AND DOMAIN DECOMPOSITION FOR PROTEIN MODELING EQUATIONS MICHAEL HOLST AND FAISAL SAIED Abstract. We consider multigrid and domain decomposition methods for the numerical
More informationAn Introduction to Algebraic Multigrid (AMG) Algorithms Derrick Cerwinsky and Craig C. Douglas 1/84
An Introduction to Algebraic Multigrid (AMG) Algorithms Derrick Cerwinsky and Craig C. Douglas 1/84 Introduction Almost all numerical methods for solving PDEs will at some point be reduced to solving A
More informationGradient Computations for Optimal Design of Turbine Blades
Gradient Computations for Optimal Design of Turbine Blades K. Arens 1, P.Rentrop 1, and S.O.Stoll 2 1 TU München, Zentrum Mathematik arens@ma.tum.de 2 TU Karlsruhe, IWRMM Summary. The optimal profile of
More informationProgramming, numerics and optimization
Programming, numerics and optimization Lecture C-3: Unconstrained optimization II Łukasz Jankowski ljank@ippt.pan.pl Institute of Fundamental Technological Research Room 4.32, Phone +22.8261281 ext. 428
More informationAdjoint code development and optimization using automatic differentiation (AD)
Adjoint code development and optimization using automatic differentiation (AD) Praveen. C Computational and Theoretical Fluid Dynamics Division National Aerospace Laboratories Bangalore - 560 037 CTFD
More informationminimize x subject to (x 2)(x 4) u,
Math 6366/6367: Optimization and Variational Methods Sample Preliminary Exam Questions 1. Suppose that f : [, L] R is a C 2 -function with f () on (, L) and that you have explicit formulae for
More informationStability-Constrained Aerodynamic Shape Optimization of a Flying Wing Configuration
13th AIAA/ISSMO Multidisciplinary Analysis Optimization Conference 13-15 September 2010, Fort Worth, Texas AIAA 2010-9199 13th AIAA/ISSMO Multidisciplinary Analysis and Optimization Conference September
More informationA Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems
A Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems Etereldes Gonçalves 1, Tarek P. Mathew 1, Markus Sarkis 1,2, and Christian E. Schaerer 1 1 Instituto de Matemática Pura
More informationDetermination of Feasible Directions by Successive Quadratic Programming and Zoutendijk Algorithms: A Comparative Study
International Journal of Mathematics And Its Applications Vol.2 No.4 (2014), pp.47-56. ISSN: 2347-1557(online) Determination of Feasible Directions by Successive Quadratic Programming and Zoutendijk Algorithms:
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 informationNew constructions of domain decomposition methods for systems of PDEs
New constructions of domain decomposition methods for systems of PDEs Nouvelles constructions de méthodes de décomposition de domaine pour des systèmes d équations aux dérivées partielles V. Dolean?? F.
More informationThe Lifted Newton Method and Its Use in Optimization
The Lifted Newton Method and Its Use in Optimization Moritz Diehl Optimization in Engineering Center (OPTEC), K.U. Leuven, Belgium joint work with Jan Albersmeyer (U. Heidelberg) ENSIACET, Toulouse, February
More informationTotal Pressure Losses Minimization in Turbomachinery Cascades, using a New Continuous Adjoint Formulation
Total Pressure Losses Minimization in Turbomachinery Cascades, using a New Continuous Adjoint Formulation D.I. Papadimitriou, K.C. Giannakoglou Laboratory of Thermal Turbomachines, National Technical University
More informationDirichlet-Neumann and Neumann-Neumann Methods
Dirichlet-Neumann and Neumann-Neumann Methods Felix Kwok Hong Kong Baptist University Introductory Domain Decomposition Short Course DD25, Memorial University of Newfoundland July 22, 2018 Outline Methods
More informationKey words. minimization, nonlinear optimization, large-scale optimization, constrained optimization, trust region methods, quasi-newton methods
SIAM J. OPTIM. c 1998 Society for Industrial and Applied Mathematics Vol. 8, No. 3, pp. 682 706, August 1998 004 ON THE IMPLEMENTATION OF AN ALGORITHM FOR LARGE-SCALE EQUALITY CONSTRAINED OPTIMIZATION
More information