Airfoil shape optimization using adjoint method and automatic differentiation
|
|
- Warren Powers
- 5 years ago
- Views:
Transcription
1 Airfoil shape optimization using adjoint method and automatic differentiation Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore AeSI CFD Symposium August, 29 Praveen. C (TIFR-CAM) Shape Optimization AeSI, Aug 29 1 / 25
2 Objectives and controls Objective function I(β) = I(β, Q) mathematical representation of system performance Control variables β Parametric controls β R n Infinite dimensional controls β : X Y Shape β set of admissible shapes State variable Q: solution of an ODE or PDE R(β, Q) = = Q = Q(β) Praveen. C (TIFR-CAM) Shape Optimization AeSI, Aug 29 2 / 25
3 Mathematical formulation Constrained minimization problem min β I(β, Q) subject to R(β, Q) = Find δβ such that δi < (to first order) δi = I [ I I δβ + β Q δq = β + I ] Q δβ Q β }{{} G Steepest descent δβ = ɛg, ɛ > δi = ɛgg = ɛ G 2 < How to compute gradient G cheaply and accurately? Praveen. C (TIFR-CAM) Shape Optimization AeSI, Aug 29 3 / 25
4 Elements of shape optimization 1 Shape parameterization 2 Surface grid generation/deformation 3 Domain grid generation/deformation 4 Flow solution (Euler/Navier-Stokes solver) 5 Adjoint flow solution 6 Optimization method Shape parameters β Surface grid X s Volume grid X CFD solution Q I di dβ = di dx dx s dx dx s dβ Praveen. C (TIFR-CAM) Shape Optimization AeSI, Aug 29 4 / 25
5 Adjoint approach For shape optimization: I = I(X, Q) di dx = I X + I Q Q X Flow sensitivity Q X ; costly to evaluate Differentiate state equation R(X, Q) = R X + R Q Q X = Introducing an adjoint variable Ψ, we can write di dx = I X + I Q Q X + Ψ [ R X + R Q ] Q X Praveen. C (TIFR-CAM) Shape Optimization AeSI, Aug 29 5 / 25
6 Adjoint approach Collect terms involving the flow sensitivity di dx = I R + Ψ X X + Choose Ψ so that flow sensitivity vanishes [ I R + Ψ Q Q ] Q X I R + Ψ Q Q = or ( ) R Ψ + Q ( ) I = Q Gradient di dx = I R + Ψ X X Praveen. C (TIFR-CAM) Shape Optimization AeSI, Aug 29 6 / 25
7 Optimization steps β = X s = X Solve the flow equations to steady-state dq dt + R(X, Q) = = Q, I(X, Q) Solve adjoint equations to steady-state dψ dt + ( ) R Ψ + Q Compute gradient wrt grid X ( ) I = = Ψ Q di dx = I R + Ψ X X di dβ = di dx dx s dx dx s dβ = β β ɛ di dβ Praveen. C (TIFR-CAM) Shape Optimization AeSI, Aug 29 7 / 25
8 Continuous and discrete approaches Continuous approach (differentiate and discretize) PDE Adjoint PDE Discrete adjoint Discrete approach (discretize and differentiate) PDE Discrete PDE Discrete adjoint We use the discrete approach R(X, Q) = represent the finite volume equations which are algebraic equations Use ordinary calculus to differentiate Need to compute I Q, I X, ( ) R Ψ, Q ( ) R Ψ X Praveen. C (TIFR-CAM) Shape Optimization AeSI, Aug 29 8 / 25
9 Automatic differentiation Computer code available to compute I(X, Q), R(X, Q) Code is made of composition of elementary functions T = X, T r = F r (T r 1 ) Y = F (X) = F p F p 1... F 1 (T ) Use differentiation by parts formula Ẏ = F (X)Ẋ = F p(t p 1 )F p 1(T p 2 )... F 1(T )Ẋ Automated using AD tools Computer code P Automatic Differentiation New code Ṗ Praveen. C (TIFR-CAM) Shape Optimization AeSI, Aug 29 9 / 25
10 Reverse differentiation Reverse mode computes transpose: (X, Ȳ ) X X = [F (X)] Ȳ = [F 1(T )] [F 2(T 1 )]... [F p(t p 1 )] Ȳ Forward sweep and then reverse sweep T Func: T F 1 T 1 2 T p 1 F2... Grad: T p 1, Ȳ [F p] T Tp 2 [F p 1 ]T Tp 3... T F p [F 1 ] T Forward variables T j required in reverse order: store or recompute Reverse mode useful to compute ( ) I, Q ( ) I, X ( ) R Ψ, Q ( ) R Ψ X Praveen. C (TIFR-CAM) Shape Optimization AeSI, Aug 29 1 / 25
11 Differentiation: Example A simple example f = (xy + sin x + 4)(3y 2 + 6) Computer code, f = t 1 t 1 = x t 2 = y t 3 = t 1 t 2 t 4 = sin t 1 t 5 = t 3 + t 4 t 6 = t t 7 = t 2 2 t 8 = 3t 7 t 9 = t t 1 = t 6 t 9 Praveen. C (TIFR-CAM) Shape Optimization AeSI, Aug / 25
12 F77 code: costfunc.f s ubroutine c o s t f u n c ( x, y, f ) t1 = x t2 = y t3 = t1 t2 t4 = s i n ( t1 ) t5 = t3 + t4 t6 = t5 + 4 t7 = t2 2 t8 = 3. t7 t9 = t t1 = t6 t9 f = t1 end Praveen. C (TIFR-CAM) Shape Optimization AeSI, Aug / 25
13 Automatic Differentiation: Reverse mode SUBROUTINE COSTFUNC_B(x, xb, y, yb, f, fb) t1 = x t2 = y t3 = t1*t2 t4 = SIN(t1) t5 = t3 + t4 t6 = t5 + 4 t7 = t2**2 t8 = 3.*t7 t9 = t t1b = fb t6b = t9*t1b t9b = t6*t1b t8b = t9b t7b = 3.*t8b t5b = t6b t3b = t5b t2b = t1*t3b + 2*t2*t7b t4b = t5b t1b = t2*t3b + COS(t1)*t4b yb = t2b xb = t1b fb =. END Praveen. C (TIFR-CAM) Shape Optimization AeSI, Aug / 25
14 Implementation of AD CFD code written with many subroutines Subroutines differentiated individually Then assembled together to form adjoint solver Only non-linear portions differentiated with AD numerical flux (Roe) Limiters Linear portions differentiated manually Leads to an efficient code with less memory requirements Praveen. C (TIFR-CAM) Shape Optimization AeSI, Aug / 25
15 Shape parameterization Parameterize the deformations [ ] [ ] xs x () s = y s y s () + h(ξ) = [ nx n y m β k B k (ξ) k=1 ] h(ξ) Hicks-Henne bump functions n A B ξ h(ξ) B k (ξ) = sin p (πξ q k ), q k = log(.5) log(ξ k ) Move points along normal to reference line AB ξ Exact derivatives dxs dβ computed can be Praveen. C (TIFR-CAM) Shape Optimization AeSI, Aug / 25
16 Grid deformation Interpolate displacement of surface points to interior points using RBF f(x, y) = a + a 1 x + a 2 y + N b j r r j 2 log r r j j=1 where r = (x, y) Initial grid Deformed grid Results in smooth grids Exact derivatives dx dx s can be computed Praveen. C (TIFR-CAM) Shape Optimization AeSI, Aug / 25
17 NUWTUN flow solver Based on the ISAAC code of Joseph Morrison Finite volume scheme Structured, multi-block grids Roe flux MUSCL reconstruction with Koren limiter Implicit scheme Source code of NUWTUN available online Praveen. C (TIFR-CAM) Shape Optimization AeSI, Aug / 25
18 Convergence tests Residue Adjoint residual.1 1x1-6 1x1-7 1x1-8 1x1-9 Flow residual 1x1-1 1x1-11 1x1-12 1x1-13 1x1-14 1x Number of iterations Cl Number of iterations Cl Cd Cd Convergence characteristics for the flow and adjoint solutions, and convergence of lift and drag coefficients, for RAE2822 airfoil at M =.73 Praveen. C (TIFR-CAM) Shape Optimization AeSI, Aug / 25
19 Validation of adjoint gradients 8 AD FD 6 4 Dot-product test Limiters can cause non-differentiability Koren limiter: dependance on parameter Check adjoint derivatives against finite difference Gradient Gradient Hicks-Henne parameter AD FD Hicks-Henne parameter Praveen. C (TIFR-CAM) Shape Optimization AeSI, Aug / 25
20 Test cases NACA12: M =.8, α = 1.25 o I = C d C l C l C d RAE2822: M =.729, α = 2.31 o Penalty approach Constrained minimization I = C d + C d 1 C l C l min I = C d C d s.t. C l = C l Praveen. C (TIFR-CAM) Shape Optimization AeSI, Aug 29 2 / 25
21 NACA12: Maximize L/D Initial -Cp -.5 Initial conmin_frcg optpp_q_newton steep conmin_frcg optpp_q_newton steep x/c Method I 1C d C l N fun N grad Initial conmin frcg optpp q newton steep Praveen. C (TIFR-CAM) Shape Optimization AeSI, Aug / 25
22 RAE2822: Drag minimization, penalty approach Cp.5 y.2 Initial conmin_frcg optpp_q_newton steep Initial conmin_frcg optpp_q_newton steep x x Method I 1C d C l N fun N grad Initial conmin frcg optpp q newton steep Praveen. C (TIFR-CAM) Shape Optimization AeSI, Aug / 25
23 RAE2822: Lift-constrained drag minimization Cp.5.2 Initial conmin_mfd fsqp ipopt Initial conmin_mfd fsqp ipopt x/c Method I 1C d C l N fun N grad Initial conmin mfd fsqp ipopt Praveen. C (TIFR-CAM) Shape Optimization AeSI, Aug / 25
24 RAE2822: Lift-constrained drag minimization Cp.5.2 ipopt Initial -.5 ipopt Initial x/c Praveen. C (TIFR-CAM) Shape Optimization AeSI, Aug / 25
25 Sensitivity to perturbations.4 Initial Optimized 2 Initial Optimized.3 15 Drag coefficient.2 Lift/Drag Mach number Mach number (a) (b) Variation of (a) drag coefficient and (b) L/D with Mach number for RAE2822 airfoil and optimized airfoil Need for robust aerodynamic optimization Praveen. C (TIFR-CAM) Shape Optimization AeSI, Aug / 25
Adjoint approach to optimization
Adjoint approach to optimization Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in Health, Safety
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 informationNumerical Solution of Partial Differential Equations governing compressible flows
Numerical Solution of Partial Differential Equations governing compressible flows Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore
More informationA COUPLED-ADJOINT METHOD FOR HIGH-FIDELITY AERO-STRUCTURAL OPTIMIZATION
A COUPLED-ADJOINT METHOD FOR HIGH-FIDELITY AERO-STRUCTURAL OPTIMIZATION Joaquim Rafael Rost Ávila Martins Department of Aeronautics and Astronautics Stanford University Ph.D. Oral Examination, Stanford
More informationA Crash-Course on the Adjoint Method for Aerodynamic Shape Optimization
A Crash-Course on the Adjoint Method for Aerodynamic Shape Optimization Juan J. Alonso Department of Aeronautics & Astronautics Stanford University jjalonso@stanford.edu Lecture 19 AA200b Applied Aerodynamics
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 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 informationCOMPLETE CONFIGURATION AERO-STRUCTURAL OPTIMIZATION USING A COUPLED SENSITIVITY ANALYSIS METHOD
COMPLETE CONFIGURATION AERO-STRUCTURAL OPTIMIZATION USING A COUPLED SENSITIVITY ANALYSIS METHOD Joaquim R. R. A. Martins Juan J. Alonso James J. Reuther Department of Aeronautics and Astronautics Stanford
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 informationAdjoint-based aerodynamic shape optimization
IT Licentiate theses 2003-012 Adjoint-based aerodynamic shape optimization OLIVIER AMOIGNON UPPSALA UNIVERSITY Department of Information Technology Adjoint-based aerodynamic shape optimization BY OLIVIER
More informationAn Investigation of the Attainable Efficiency of Flight at Mach One or Just Beyond
45 th Aerospace Sciences Meeting and Exhibit, January 8 11, 2007, Reno, Nevada An Investigation of the Attainable Efficiency of Flight at Mach One or Just Beyond Antony Jameson Department of Aeronautics
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 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 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 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 informationFinite volume method on unstructured grids
Finite volume method on unstructured grids Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationOptimal control problems with PDE constraints
Optimal control problems with PDE constraints Maya Neytcheva CIM, October 2017 General framework Unconstrained optimization problems min f (q) q x R n (real vector) and f : R n R is a smooth function.
More informationNumerical Optimization
Constrained Optimization Computer Science and Automation Indian Institute of Science Bangalore 560 012, India. NPTEL Course on Constrained Optimization Constrained Optimization Problem: min h j (x) 0,
More informationAERO-STRUCTURAL WING DESIGN OPTIMIZATION USING HIGH-FIDELITY SENSITIVITY ANALYSIS
AERO-STRUCTURAL WING DESIGN OPTIMIZATION USING HIGH-FIDELITY SENSITIVITY ANALYSIS Joaquim R. R. A. Martins and Juan J. Alonso Department of Aeronautics and Astronautics Stanford University, Stanford, CA
More informationLagrangian Coherent Structures (LCS)
Lagrangian Coherent Structures (LCS) CDS 140b - Spring 2012 May 15, 2012 ofarrell@cds.caltech.edu A time-dependent dynamical system ẋ (t; t 0, x 0 )=v(x(t;,t 0, x 0 ),t) x(t 0 ; t 0, x 0 )=x 0 t 2 I R
More informationChallenges and Complexity of Aerodynamic Wing Design
Chapter 1 Challenges and Complexity of Aerodynamic Wing Design Kasidit Leoviriyakit and Antony Jameson Department of Aeronautics and Astronautics Stanford University, Stanford CA kasidit@stanford.edu and
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 informationEFFICIENT AERODYNAMIC OPTIMIZATION USING HIERARCHICAL KRIGING COMBINED WITH GRADIENT
EFFICIENT AERODYNAMIC OPTIMIZATION USING HIERARCHICAL KRIGING COMBINED WITH GRADIENT Chao SONG, Xudong YANG, Wenping SONG National Key Laboratory of Science and Technology on Aerodynamic Design and Research,
More informationDiscrete and continuous adjoint method for compressible CFD J. Peter ONERA
Discrete and continuous adjoint method for compressible CFD J. Peter ONERA J. Peter 1 1 ONERA DMFN October 2, 2014 J. Peter (ONERA DMFN) Adjoint method for compressible CFD October 2, 2014 1 / 60 Outline
More informationA Multi-Dimensional Limiter for Hybrid Grid
APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy
More informationApplications of adjoint based shape optimization to the design of low drag airplane wings, including wings to support natural laminar flow
Applications of adjoint based shape optimization to the design of low drag airplane wings, including wings to support natural laminar flow Antony Jameson and Kui Ou Aeronautics & Astronautics Department,
More information5. FVM discretization and Solution Procedure
5. FVM discretization and Solution Procedure 1. The fluid domain is divided into a finite number of control volumes (cells of a computational grid). 2. Integral form of the conservation equations are discretized
More information1.2 Derivation. d p f = d p f(x(p)) = x fd p x (= f x x p ). (1) Second, g x x p + g p = 0. d p f = f x g 1. The expression f x gx
PDE-constrained optimization and the adjoint method Andrew M. Bradley November 16, 21 PDE-constrained optimization and the adjoint method for solving these and related problems appear in a wide range of
More informationA Short Essay on Variational Calculus
A Short Essay on Variational Calculus Keonwook Kang, Chris Weinberger and Wei Cai Department of Mechanical Engineering, Stanford University Stanford, CA 94305-4040 May 3, 2006 Contents 1 Definition of
More informationAerodynamic Inverse Design and Shape Optimization via Control Theory
Aerodynamic Inverse Design and Shape Optimization via Control Theory 1 1 Thomas V. Jones Professor of Engineering Department of Aeronautics & Astronautics Stanford University SciTech 2015 January 5, 2015
More informationFar Field Noise Minimization Using an Adjoint Approach
Far Field Noise Minimization Using an Adjoint Approach Markus P. Rumpfkeil and David W. Zingg University of Toronto Institute for Aerospace Studies 4925 Dufferin Street, Toronto, Ontario, M3H 5T6, Canada
More informationDifferential Equation Types. Moritz Diehl
Differential Equation Types Moritz Diehl Overview Ordinary Differential Equations (ODE) Differential Algebraic Equations (DAE) Partial Differential Equations (PDE) Delay Differential Equations (DDE) Ordinary
More informationAdjoint based multi-objective shape optimization of a transonic airfoil under uncertainties
EngOpt 2016-5 th International Conference on Engineering Optimization Iguassu Falls, Brazil, 19-23 June 2016. Adjoint based multi-objective shape optimization of a transonic airfoil under uncertainties
More informationSTUDY ON THE IMPROVED KRIGING-BASED OPTIMIZATION DESIGN METHOD
7 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES STUDY ON THE IMPROVED KRIGING-BASED OPTIMIZATION DESIGN METHOD Song Wenping, Xu Ruifei, Han Zhonghua (National Key Laboratory of Science and Technology
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 informationLecture V. Numerical Optimization
Lecture V Numerical Optimization Gianluca Violante New York University Quantitative Macroeconomics G. Violante, Numerical Optimization p. 1 /19 Isomorphism I We describe minimization problems: to maximize
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 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 informationVector Calculus, Maths II
Section A Vector Calculus, Maths II REVISION (VECTORS) 1. Position vector of a point P(x, y, z) is given as + y and its magnitude by 2. The scalar components of a vector are its direction ratios, and represent
More informationCompatible Discretization Schemes
Compatible Discretization Schemes Marc Gerritsma 1 1 Collaborators Artur Palha, Guido Oud, Jasper Kreeft & Mick Bouman Marc Gerritsma Burgers Course CFD I, 25-29 January 2010 1 / 81 What if... Why mimetic
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) Key Concepts Current Semester 1 / 25 Introduction The purpose of this section is to define
More informationFEniCS Course Lecture 13: Introduction to dolfin-adjoint. Contributors Simon Funke Patrick Farrell Marie E. Rognes
FEniCS Course Lecture 13: Introduction to dolfin-adjoint Contributors Simon Funke Patrick Farrell Marie E. Rognes 1 / 25 dolfin adjoint automatic adjoint models for FEniCS P. E. Farrell, S. W. Funke, D.
More informationTwo-Body Problem. Central Potential. 1D Motion
Two-Body Problem. Central Potential. D Motion The simplest non-trivial dynamical problem is the problem of two particles. The equations of motion read. m r = F 2, () We already know that the center of
More informationActive Flux for Advection Diffusion
Active Flux for Advection Diffusion A Miracle in CFD Hiroaki Nishikawa National Institute of Aerospace! NIA CFD Seminar! August 25, 2015 In collaboration with the University of Michigan Supported by NASA
More informationFinal: Solutions Math 118A, Fall 2013
Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or
More informationAn introduction to Mathematical Theory of Control
An introduction to Mathematical Theory of Control Vasile Staicu University of Aveiro UNICA, May 2018 Vasile Staicu (University of Aveiro) An introduction to Mathematical Theory of Control UNICA, May 2018
More informationSeparation of Variables in Linear PDE: One-Dimensional Problems
Separation of Variables in Linear PDE: One-Dimensional Problems Now we apply the theory of Hilbert spaces to linear differential equations with partial derivatives (PDE). We start with a particular example,
More informationLINEAR AND NONLINEAR PROGRAMMING
LINEAR AND NONLINEAR PROGRAMMING Stephen G. Nash and Ariela Sofer George Mason University The McGraw-Hill Companies, Inc. New York St. Louis San Francisco Auckland Bogota Caracas Lisbon London Madrid Mexico
More informationFinite difference method for elliptic problems: I
Finite difference method for elliptic problems: I Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationIn this chapter, we study the calculus of vector fields.
16 VECTOR CALCULUS VECTOR CALCULUS In this chapter, we study the calculus of vector fields. These are functions that assign vectors to points in space. VECTOR CALCULUS We define: Line integrals which can
More informationEdwin van der Weide and Magnus Svärd. I. Background information for the SBP-SAT scheme
Edwin van der Weide and Magnus Svärd I. Background information for the SBP-SAT scheme As is well-known, stability of a numerical scheme is a key property for a robust and accurate numerical solution. Proving
More informationAn Introduction to Partial Differential Equations
An Introduction to Partial Differential Equations Ryan C. Trinity University Partial Differential Equations Lecture 1 Ordinary differential equations (ODEs) These are equations of the form where: F(x,y,y,y,y,...)
More informationReview and Unification of Methods for Computing Derivatives of Multidisciplinary Computational Models
his is a preprint of the following article, which is available from http://mdolabenginumichedu J R R A Martins and J Hwang Review and unification of methods for computing derivatives of multidisciplinary
More informationLogistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations
Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations S. Y. Ha and J. Park Department of Mathematical Sciences Seoul National University Sep 23, 2013 Contents 1 Logistic Map 2 Euler and
More informationActive sets, steepest descent, and smooth approximation of functions
Active sets, steepest descent, and smooth approximation of functions Dmitriy Drusvyatskiy School of ORIE, Cornell University Joint work with Alex D. Ioffe (Technion), Martin Larsson (EPFL), and Adrian
More informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationCFD Based Optimization of a Vertical Axis Wind Turbine
SCHOOL OF ENGINEERING INSTITUTE FOR ENERGY SYSTEMS CFD Based Optimization of a Vertical Axis Wind Turbine Emmanouil Falagkaris Supervisory Team Prof. David Ingram, Dr Ignazio Maria Viola, Dr Daniel Friedrich
More informationGradient-Based Optimum Aerodynamic Design Using Adjoint Methods
Gradient-Based Optimum Aerodynamic Design Using Adjoint Methods Lei Xie Dissertation submitted to the Faculty of the Virginia Polytechnic Institute and State University in partial fulfillment of the requirements
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 informationChapter 3 Numerical Methods
Chapter 3 Numerical Methods Part 2 3.2 Systems of Equations 3.3 Nonlinear and Constrained Optimization 1 Outline 3.2 Systems of Equations 3.3 Nonlinear and Constrained Optimization Summary 2 Outline 3.2
More information. Frobenius-Perron Operator ACC Workshop on Uncertainty Analysis & Estimation. Raktim Bhattacharya
.. Frobenius-Perron Operator 2014 ACC Workshop on Uncertainty Analysis & Estimation Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. uq.tamu.edu
More informationON ESTIMATION OF TEMPERATURE UNCERTAINTY USING THE SECOND ORDER ADJOINT PROBLEM
ON ESTIMATION OF TEMPERATURE UNCERTAINTY USING THE SECOND ORDER ADJOINT PROBLEM Aleksey K. Alekseev a, Michael I. Navon b,* a Department of Aerodynamics and Heat Transfer, RSC, ENERGIA, Korolev (Kaliningrad),
More informationFirst, Second, and Third Order Finite-Volume Schemes for Diffusion
First, Second, and Third Order Finite-Volume Schemes for Diffusion Hiro Nishikawa 51st AIAA Aerospace Sciences Meeting, January 10, 2013 Supported by ARO (PM: Dr. Frederick Ferguson), NASA, Software Cradle.
More informationMultiple Integrals and Vector Calculus (Oxford Physics) Synopsis and Problem Sets; Hilary 2015
Multiple Integrals and Vector Calculus (Oxford Physics) Ramin Golestanian Synopsis and Problem Sets; Hilary 215 The outline of the material, which will be covered in 14 lectures, is as follows: 1. Introduction
More informationAdditive Manufacturing Module 8
Additive Manufacturing Module 8 Spring 2015 Wenchao Zhou zhouw@uark.edu (479) 575-7250 The Department of Mechanical Engineering University of Arkansas, Fayetteville 1 Evaluating design https://www.youtube.com/watch?v=p
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS 1. Basic Terminology A differential equation is an equation that contains an unknown function together with one or more of its derivatives. 1 Examples: 1. y = 2x + cos x 2. dy dt
More informationEntropy stable schemes for compressible flows on unstructured meshes
Entropy stable schemes for compressible flows on unstructured meshes Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore deep@math.tifrbng.res.in http://math.tifrbng.res.in/
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 informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS Basic Terminology A differential equation is an equation that contains an unknown function together with one or more of its derivatives. 1 Examples: 1. y = 2x + cos x 2. dy dt =
More informationLesson Rigid Body Dynamics
Lesson 8 Rigid Body Dynamics Lesson 8 Outline Problem definition and motivations Dynamics of rigid bodies The equation of unconstrained motion (ODE) User and time control Demos / tools / libs Rigid Body
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 informationProper Orthogonal Decomposition Extensions and Their Applications in Steady Aerodynamics
Proper Orthogonal Decomposition Extensions and Their Applications in Steady Aerodynamics by Bui Thanh Tan B.Eng. in Aeronautics Ho Chi Minh City University of Technology, Vietnam, 2 SUBMITTED TO THE SMA
More informationReduced-Order Multiobjective Optimal Control of Semilinear Parabolic Problems. Laura Iapichino Stefan Trenz Stefan Volkwein
Universität Konstanz Reduced-Order Multiobjective Optimal Control of Semilinear Parabolic Problems Laura Iapichino Stefan Trenz Stefan Volkwein Konstanzer Schriften in Mathematik Nr. 347, Dezember 2015
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS Chapter 1 Introduction and Basic Terminology Most of the phenomena studied in the sciences and engineering involve processes that change with time. For example, it is well known
More informationA Study of Transonic Flow and Airfoils. Presented by: Huiliang Lui 30 th April 2007
A Study of Transonic Flow and Airfoils Presented by: Huiliang Lui 3 th April 7 Contents Background Aims Theory Conservation Laws Irrotational Flow Self-Similarity Characteristics Numerical Modeling Conclusion
More informationClassification of partial differential equations and their solution characteristics
9 TH INDO GERMAN WINTER ACADEMY 2010 Classification of partial differential equations and their solution characteristics By Ankita Bhutani IIT Roorkee Tutors: Prof. V. Buwa Prof. S. V. R. Rao Prof. U.
More informationScalable parallel approach for high-fidelity steady-state aeroelastic analysis and adjoint derivative computations
This is a preprint of the following article, which is available from http://mdolab.engin.umich.edu Gaetan K. W. Kenway, Graeme J. Kennedy, and Joaquim R. R. A. Martins, Scalable Parallel Approach for High-
More informationNumerical Optimal Control Overview. Moritz Diehl
Numerical Optimal Control Overview Moritz Diehl Simplified Optimal Control Problem in ODE path constraints h(x, u) 0 initial value x0 states x(t) terminal constraint r(x(t )) 0 controls u(t) 0 t T minimize
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS 1. Basic Terminology A differential equation is an equation that contains an unknown function together with one or more of its derivatives. 1 Examples: 1. y = 2x + cos x 2. dy dt
More informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
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 informationAP Calculus 2004 AB FRQ Solutions
AP Calculus 4 AB FRQ Solutions Louis A. Talman, Ph. D. Emeritus Professor of Mathematics Metropolitan State University of Denver July, 7 Problem. Part a The function F (t) = 8 + 4 sin(t/) gives the rate,
More informationTowards Reduced Order Modeling (ROM) for Gust Simulations
Towards Reduced Order Modeling (ROM) for Gust Simulations S. Görtz, M. Ripepi DLR, Institute of Aerodynamics and Flow Technology, Braunschweig, Germany Deutscher Luft und Raumfahrtkongress 2017 5. 7. September
More informationPOD/DEIM 4DVAR Data Assimilation of the Shallow Water Equation Model
nonlinear 4DVAR 4DVAR Data Assimilation of the Shallow Water Equation Model R. Ştefănescu and Ionel M. Department of Scientific Computing Florida State University Tallahassee, Florida May 23, 2013 (Florida
More informationAlberto Bressan. Department of Mathematics, Penn State University
Non-cooperative Differential Games A Homotopy Approach Alberto Bressan Department of Mathematics, Penn State University 1 Differential Games d dt x(t) = G(x(t), u 1(t), u 2 (t)), x(0) = y, u i (t) U i
More informationAn Introduction to the Discontinuous Galerkin Method
An Introduction to the Discontinuous Galerkin Method Krzysztof J. Fidkowski Aerospace Computational Design Lab Massachusetts Institute of Technology March 16, 2005 Computational Prototyping Group Seminar
More informationMathematics Qualifying Exam Study Material
Mathematics Qualifying Exam Study Material The candidate is expected to have a thorough understanding of engineering mathematics topics. These topics are listed below for clarification. Not all instructors
More informationCurved grid generation and DG computation for the DLR-F11 high lift configuration
ECCOMAS Congress 2016, Crete Island, Greece, 5-10 June 2016 Curved grid generation and DG computation for the DLR-F11 high lift configuration Ralf Hartmann 1, Harlan McMorris 2, Tobias Leicht 1 1 DLR (German
More informationPROJECTION METHODS FOR DYNAMIC MODELS
PROJECTION METHODS FOR DYNAMIC MODELS Kenneth L. Judd Hoover Institution and NBER June 28, 2006 Functional Problems Many problems involve solving for some unknown function Dynamic programming Consumption
More informationEfficient gradient computation of a misfit function for FWI using the adjoint method
Efficient gradient computation of a misfit function for FWI using the adjoint method B. Galuzzi, department of Earth Sciences A. Desio, University of Milan E. Stucchi, department of Earth Sciences A. Desio,
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationDrag Characteristics of a Low-Drag Low-Boom Supersonic Formation Flying Concept
Drag Characteristics of a Low-Drag Low-Boom Supersonic Formation Flying Concept Yuichiro Goto, Shigeru Obayashi and Yasuaki Kohama Tohoku University, Sendai, Japan In this paper, a new concept for low-drag,
More informationChapter 9b: Numerical Methods for Calculus and Differential Equations. Initial-Value Problems Euler Method Time-Step Independence MATLAB ODE Solvers
Chapter 9b: Numerical Methods for Calculus and Differential Equations Initial-Value Problems Euler Method Time-Step Independence MATLAB ODE Solvers Acceleration Initial-Value Problems Consider a skydiver
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 informationHp-Adaptivity on Anisotropic Meshes for Hybridized Discontinuous Galerkin Scheme
Hp-Adaptivity on Anisotropic Meshes for Hybridized Discontinuous Galerkin Scheme Aravind Balan, Michael Woopen and Georg May AICES Graduate School, RWTH Aachen University, Germany 22nd AIAA Computational
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS Basic Terminology A differential equation is an equation that contains an unknown function together with one or more of its derivatives. 1 Examples: 1. y = 2x + cos x 2. dy dt =
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 informationAnswers to Problem Set # 01, MIT (Winter-Spring 2018)
Answers to Problem Set # 01, 18.306 MIT (Winter-Spring 2018) Rodolfo R. Rosales (MIT, Math. Dept., room 2-337, Cambridge, MA 02139) February 28, 2018 Contents 1 Nonlinear solvable ODEs 2 1.1 Statement:
More informationInverse problems Total Variation Regularization Mark van Kraaij Casa seminar 23 May 2007 Technische Universiteit Eindh ove n University of Technology
Inverse problems Total Variation Regularization Mark van Kraaij Casa seminar 23 May 27 Introduction Fredholm first kind integral equation of convolution type in one space dimension: g(x) = 1 k(x x )f(x
More informationNumerical algorithms for one and two target optimal controls
Numerical algorithms for one and two target optimal controls Sung-Sik Kwon Department of Mathematics and Computer Science, North Carolina Central University 80 Fayetteville St. Durham, NC 27707 Email:
More information