Unsteady Discrete Adjoint Formulation for Two-Dimensional Flow Problems with Deforming Meshes

A method to apply the discrete adjoint for computing sensitivity derivatives in two-dimensional unsteady flow problems is presented. The approach is to first develop a forward or tangent linearization of the nonlinear flow problem in which each individual component building up the complete flow solution is differentiated against the design variables using the chain rule. The reverse or adjoint linearization is then constructed by transposing and reversing the order of multiplication of the forward problem. The developed algorithm is very general in that it applies directly to the arbitrary Lagrangian-Eulerian form of the governing equations and includes the effect of deforming meshes in unsteady flows. It is shown that an unsteady adjoint formulation is essentially a single backward integration in time and that the cost of constructing the final sensitivity vector is close to that of solving the unsteady flow problem integrated forward in time. It is also shown that the unsteady adjoint formulation can be applied to time-integration schemes of different orders of accuracy with minimal changes to the base formulation. The developed technique is then applied to three optimization examples, the first in which the shape of a pitching airfoil is morphed to match a target time-dependent load profile, the second in which the shape is optimized to match a target time-dependent pressure profile, and the last in which the time-dependent drag profile is minimized without any loss in lift.

[1]  Antony Jameson,et al.  Optimal Control of Unsteady Flows Using Time Accurate and Non-Linear Frequency Domain Methods , 2003 .

[2]  Dimitri J. Mavriplis,et al.  Navier-stokes applications to high-lift airfoil analysis , 1995 .

[3]  A. Jameson,et al.  A COMPARISON OF THE CONTINUOUS AND DISCRETE ADJOINT APPROACH TO AUTOMATIC AERODYNAMIC OPTIMIZATION , 2000 .

[4]  Antony Jameson,et al.  Computational Fluid Dynamics for Aerodynamic Design: Its Current and Future Impact , 2001 .

[5]  P. Roe.  Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[6]  Dimitri J. Mavriplis,et al.  Higher-order Time Integration Schemes for Aeroelastic Applications on Unstructured Meshes , 2006 .

[7]  D. Mavriplis.  Multigrid Solution of the Discrete Adjoint for Optimization Problems on Unstructured Meshes , 2006 .

[8]  Antony Jameson,et al.  OPTIMAL CONTROL OF UNSTEADY FLOWS USING A TIME ACCURATE METHOD , 2002 .

[9]  Antony Jameson,et al.  NON-LINEAR FREQUENCY DOMAIN BASED OPTIMUM SHAPE DESIGN FOR UNSTEADY THREE-DIMENSIONAL FLOWS , 2006 .

[10]  M. Giles,et al.  Adjoint Code Developments Using the Exact Discrete Approach , 2001 .

[11]  D. Mavriplis.  Discrete Adjoint-Based Approach for Optimization Problems on Three-Dimensional Unstructured Meshes , 2007 .

[12]  Dimitri J. Mavriplis,et al.  MULTIGRID TECHNIQUES FOR UNSTRUCTURED MESHES , 1995 .

[13]  D. Darmofal,et al.  An implicit, exact dual adjoint solution method for turbulent flows on unstructured grids , 2004 .

[14]  Rainald Löhner,et al.  An Adjoint-Based Design Methodology for CFD Optimization Problems , 2003 .

[15]  A. Jameson,et al.  Aerodynamic shape optimization techniques based on control theory , 1998 .

[16]  R. Errico.  What is an adjoint model , 1997 .

[17]  D. Venditti,et al.  Anisotropic grid adaptation for functional outputs: application to two-dimensional viscous flows , 2003 .

[18]  W. K. Anderson,et al.  Recent improvements in aerodynamic design optimization on unstructured meshes , 2001 .

[19]  Oktay Baysal,et al.  LIMIT-CYCLE SHAPE -OPTIMIZATION USING TIME-DEPENDE.NT TRANSONIC EQUATION , 1999 .

[20]  M. Giles,et al.  Algorithm Developments for Discrete Adjoint Methods , 2003 .

[21]  Luis Santos,et al.  Aerodynamic shape optimization using the adjoint method , 2007 .

[22]  J. Batina.  Unsteady Euler airfoil solutions using unstructured dynamic meshes , 1989 .

[23]  J. Peraire,et al.  Practical Three-Dimensional Aerodynamic Design and Optimization Using Unstructured Meshes , 1997 .

[24]  R. Dwight,et al.  Eect of Various Approximations of the Discrete Adjoint on Gradient-Based Optimization , 2006 .

[25]  David G. Wilson,et al.  2. Constrained Optimization , 2005 .

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

[27]  E. Nielsen,et al.  Efficient Construction of Discrete Adjoint Operators on Unstructured Grids Using Complex Variables , 2005 .

[28]  Charbel Farhat,et al.  Design and analysis of ALE schemes with provable second-order time-accuracy for inviscid and viscous flow simulations , 2003 .

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

[30]  D. Venditti,et al.  Grid adaptation for functional outputs: application to two-dimensional inviscid flows , 2002 .

[31]  David L. Darmofal,et al.  An Exact Dual Adjoint Solution Method for Turbulent Flows on Unstructured Grids , 2003 .

[32]  Dimitri J. Mavriplis,et al.  Construction of the discrete geometric conservation law for high-order time-accurate simulations on dynamic meshes , 2006, J. Comput. Phys..

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