An adjoint method for shape optimization in unsteady viscous flows

A new method for shape optimization for unsteady viscous flows is presented. It is based on the continuous adjoint approach using a time accurate method and is capable of handling both inverse and direct objective functions. The objective function is minimized or maximized subject to the satisfaction of flow equations. The shape of the body is parametrized via a Non-Uniform Rational B-Splines (NURBS) curve and is updated by using the gradients obtained from solving the flow and adjoint equations. A finite element method based on streamline-upwind Petrov/Galerkin (SUPG) and pressure stabilized Petrov/Galerkin (PSPG) stabilization techniques is used to solve both the flow and adjoint equations. The method has been implemented and tested for the design of airfoils, based on enhancing its time-averaged aerodynamic coefficients. Interesting shapes are obtained, especially when the objective is to produce high performance airfoils. The effect of the extent of the window of time integration of flow and adjoint equations on the design process is studied. It is found that when the window of time integration is insufficient, the gradients are most likely to be erroneous.

[1]  Rainald Löhner,et al.  An adjoint‐based design methodology for CFD problems , 2004 .

[2]  Antony Jameson,et al.  Aerodynamic design via control theory , 1988, J. Sci. Comput..

[3]  David J. Lucia,et al.  Rocket Nozzle Flow Control Using a Reduced-Order Computational Fluid Dynamics Model , 2002 .

[4]  Michael B. Giles,et al.  The harmonic adjoint approach to unsteady turbomachinery design , 2002 .

[5]  Giuseppe Pontrelli,et al.  Flow of a viscoelastic fluid between two rotating circular cylinders subject to suction or injection , 1997 .

[6]  S. Obayashi,et al.  Aerodynamic inverse optimization with genetic algorithms , 1996, Proceedings of the IEEE International Conference on Industrial Technology (ICIT'96).

[7]  A. Jameson,et al.  Aerodynamic Shape Optimization of Complex Aircraft Configurations via an Adjoint Formulation , 1996 .

[8]  Siva Nadarajah,et al.  Multi-objective aerodynamic shape optimization for unsteady viscous flows , 2010 .

[9]  Gregory M. Hulbert,et al.  New Methods in Transient Analysis , 1992 .

[10]  M. D. Salas,et al.  Airfoil Design and Optimization by the One-Shot Method , 1995 .

[11]  O. C. Zienkiewicz,et al.  A simple error estimator and adaptive procedure for practical engineerng analysis , 1987 .

[12]  Niles A. Pierce,et al.  An Introduction to the Adjoint Approach to Design , 2000 .

[13]  Juan J. Alonso,et al.  Application of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations , 2002 .

[14]  M. Heinkenschloss,et al.  Optimal control of unsteady compressible viscous flows , 2002 .

[15]  P. Vaidyanathan Generalizations of the sampling theorem: Seven decades after Nyquist , 2001 .

[16]  Kenneth C. Hall,et al.  Linearized Euler predictions of unsteady aerodynamic loads in cascades , 1993 .

[17]  A. J. Jerri The Shannon sampling theorem—Its various extensions and applications: A tutorial review , 1977, Proceedings of the IEEE.

[18]  Gerald Farin,et al.  Curves and surfaces for computer aided geometric design , 1990 .

[19]  Jacques Periaux,et al.  Active Control and Drag Optimization for Flow Past a Circular Cylinder , 2000 .

[20]  J. C. Newman,et al.  Computationally efficient, numerically exact design space derivatives via the complex Taylor's series expansion method , 2003 .

[21]  Paul-Louis George,et al.  ASPECTS OF 2-D DELAUNAY MESH GENERATION , 1997 .

[22]  J. Nocedal,et al.  A Limited Memory Algorithm for Bound Constrained Optimization , 1995, SIAM J. Sci. Comput..

[23]  Max Gunzburger,et al.  SENSITIVITIES, ADJOINTS AND FLOW OPTIMIZATION , 1999 .

[24]  M. Rumpfkeil,et al.  The optimal control of unsteady flows with a discrete adjoint method , 2010 .

[25]  Ionel M. Navon,et al.  Suppression of vortex shedding for flow around a circular cylinder using optimal control , 2002 .

[26]  Dimitri J. Mavriplis,et al.  Unsteady Discrete Adjoint Formulation for Two-Dimensional Flow Problems with Deforming Meshes , 2008 .

[27]  Raymond M. Hicks,et al.  Wing design by numerical optimization , 1977 .

[28]  Antony Jameson,et al.  Optimum Shape Design for Unsteady Flows with Time-Accurate Continuous and Discrete Adjoint Methods , 2007 .

[29]  M. Kawahara,et al.  Shape Optimization of Body Located in Incompressible Navier--Stokes Flow Based on Optimal Control Theory , 2000 .

[30]  Gene Hou,et al.  Shape-sensitivity analysis and design optimization of linear, thermoelastic solids , 1992 .

[31]  Gene Hou,et al.  First- and Second-Order Aerodynamic Sensitivity Derivatives via Automatic Differentiation with Incremental Iterative Methods , 1996 .

[32]  S. Mittal,et al.  Computation of unsteady incompressible flows with the stabilized finite element methods: Space-time formulations, iterative strategies and massively parallel implementations , 1992 .

[33]  Antony Jameson,et al.  OPTIMUM SHAPE DESIGN FOR UNSTEADY FLOWS USING TIME ACCURATE AND NON-LINEAR FREQUENCY DOMAIN METHODS , 2003 .

[34]  Wr Graham,et al.  OPTIMAL CONTROL OF VORTEX SHEDDING USING LOW-ORDER MODELS. PART I-OPEN-LOOP MODEL DEVELOPMENT , 1999 .

[35]  S. Mittal,et al.  A stabilized finite element method for shape optimization in low Reynolds number flows , 2007 .

[36]  Olivier Soucy,et al.  SONIC BOOM REDUCTION VIA REMOTE INVERSE ADJOINT APPROACH , 2007 .

[37]  F. Guibault,et al.  Optimized Nonuniform Rational B-Spline Geometrical Representation for Aerodynamic Design of Wings , 2001 .

[38]  O. Pironneau,et al.  SHAPE OPTIMIZATION IN FLUID MECHANICS , 2004 .

[39]  S. Ravindran Control of flow separation over a forward-facing step by model reduction , 2002 .

[40]  A. J. Jerri Correction to "The Shannon sampling theorem—Its various extensions and applications: A tutorial review" , 1979 .

[41]  Rainald Löhner,et al.  CFD SHAPE OPTIMIZATION USING AN INCOMPLETE-GRADIENT ADJOINT FORMULATION , 2001 .

[42]  W. K. Anderson,et al.  Aerodynamic design optimization on unstructured grids with a continuous adjoint formulation , 1997 .

[43]  Laurent Cordier,et al.  Control of the cylinder wake in the laminar regime by Trust-Region methods and POD Reduced Order Models , 2005, Proceedings of the 44th IEEE Conference on Decision and Control.

[44]  Sanjay Mittal,et al.  Optimal airfoil shapes for low Reynolds number flows , 2009 .

[45]  Bijan Mohammadi,et al.  Optimization of aerodynamic and acoustic performances of supersonic civil transports , 2004 .

[46]  M. Heinkenschloss,et al.  Shape optimization in steady blood flow: A numerical study of non-Newtonian effects , 2005, Computer methods in biomechanics and biomedical engineering.

[47]  A Jameson,et al.  Computational Aerodynamics for Aircraft Design , 1989, Science.

[48]  S. Mittal,et al.  Incompressible flow computations with stabilized bilinear and linear equal-order-interpolation velocity-pressure elements , 1992 .

[49]  David W. Zingg,et al.  A General Framework for the Optimal Control of Unsteady Flows with Applications , 2007 .