The continuous adjoint approach to the k–ω SST turbulence model with applications in shape optimization

In this article, the gradient of aerodynamic objective functions with respect to design variables, in problems governed by the incompressible Navier–Stokes equations coupled with the k–ω SST turbulence model, is computed using the continuous adjoint method, for the first time. Shape optimization problems for minimizing drag, in external aerodynamics (flows around isolated airfoils), or viscous losses in internal aerodynamics (duct flows) are considered. Sensitivity derivatives computed with the proposed adjoint method are compared to those computed with finite differences or a continuous adjoint variant based on the frequently used assumption of frozen turbulence; the latter proves the need for differentiating the turbulence model. Geometries produced by optimization runs performed with sensitivities computed by the proposed method and the ‘frozen turbulence’ assumption are also compared to quantify the gain from formulating and solving the adjoint to the turbulence model equations.

[1]  D. Spalding,et al.  A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows , 1972 .

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

[3]  F. Menter Two-equation eddy-viscosity turbulence models for engineering applications , 1994 .

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

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

[6]  M. Giles,et al.  Adjoint equations in CFD: duality, boundary conditions and solution behaviour , 1997 .

[7]  A. Jameson,et al.  Optimum Aerodynamic Design Using the Navier–Stokes Equations , 1997 .

[8]  W. K. Anderson,et al.  Airfoil Design on Unstructured Grids for Turbulent Flows , 1999 .

[9]  Andreas Griewank,et al.  Algorithm 799: revolve: an implementation of checkpointing for the reverse or adjoint mode of computational differentiation , 2000, TOMS.

[10]  D K Smith,et al.  Numerical Optimization , 2001, J. Oper. Res. Soc..

[11]  F. Menter,et al.  Ten Years of Industrial Experience with the SST Turbulence Model , 2003 .

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

[13]  Byung Joon Lee,et al.  Automated design methodology of turbulent internal flow using discrete adjoint formulation , 2007 .

[14]  Andrea Walther,et al.  Automatic Differentiation of an Entire Design Chain for Aerodynamic Shape Optimization , 2007 .

[15]  Kyriakos C. Giannakoglou,et al.  A continuous adjoint method with objective function derivatives based on boundary integrals, for inviscid and viscous flows , 2007 .

[16]  C. Othmer A continuous adjoint formulation for the computation of topological and surface sensitivities of ducted flows , 2008 .

[17]  Kyriakos C. Giannakoglou,et al.  Continuous adjoint approach to the Spalart–Allmaras turbulence model for incompressible flows , 2009 .

[18]  Kyriakos C. Giannakoglou,et al.  Adjoint wall functions: A new concept for use in aerodynamic shape optimization , 2010, J. Comput. Phys..

[19]  Ralf Hartmann,et al.  Adjoint-based error estimation and adaptive mesh refinement for the RANS and k-ω turbulence model equations , 2011, J. Comput. Phys..

[20]  Enrique Zuazua,et al.  Continuous adjoint approach for the Spalart−Allmaras model in aerodynamic optimization , 2012 .

[21]  Ionel M. Navon,et al.  Impact of non‐smooth observation operators on variational and sequential data assimilation for a limited‐area shallow‐water equation model , 2012 .

[22]  F. Thiele,et al.  Discrete Adjoint based Sensitivity Analysis for Optimal Flow Control of a 3D High-Lift Configuration , 2013 .

[23]  Uwe Naumann,et al.  A Discrete Adjoint Model for OpenFOAM , 2013, ICCS.

[24]  Juan J. Alonso,et al.  A hybrid adjoint approach applied to turbulent flow simulations , 2013 .

[25]  Kyriakos C. Giannakoglou,et al.  The continuous adjoint approach to the k–ε turbulence model for shape optimization and optimal active control of turbulent flows , 2015 .

[26]  Johanna Weiss,et al.  Optimal Shape Design For Elliptic Systems , 2016 .