Shape sensitivity analysis for the compressible Navier–Stokes equations via discontinuous Galerkin methods

Abstract This paper describes the formulation of adjoint-based sensitivity analysis and optimization techniques for high-order discontinuous Galerkin discretizations applied to viscous compressible flow. The flow is modeled by the compressible Navier–Stokes equations and the discretization of the viscous flux terms is based on an explicit symmetric interior penalty method. The discrete adjoint equation arising from the sensitivity derivative calculation is formulated consistently with the analysis problem, including the treatment of boundary conditions. In this regard, the influence on the sensitivity derivatives resulting from the deformation of curved boundary elements must be taken into account. Several numerical examples are used to examine the order of accuracy (up to p = 4 ) achieved by the current DG discretizations, to verify the derived adjoint sensitivity formulations, and to demonstrate the effectiveness of the discrete adjoint algorithm in steady and unsteady design optimization for both two- and three-dimensional viscous design problems.

[1]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[2]  Dimitri J. Mavriplis,et al.  Multigrid algorithms for high-order discontinuous Galerkin discretizations of the compressible Navier-Stokes equations , 2009, J. Comput. Phys..

[3]  D. Mavriplis,et al.  Robust Computation of Turbulent Flows Using a Discontinuous Galerkin Method , 2012 .

[4]  R. Bouard,et al.  Experimental determination of the main features of the viscous flow in the wake of a circular cylinder in uniform translation. Part 1. Steady flow , 1977, Journal of Fluid Mechanics.

[5]  Dimitri J. Mavriplis,et al.  High-order discontinuous Galerkin methods using an hp-multigrid approach , 2006, J. Comput. Phys..

[6]  Dimitri J. Mavriplis,et al.  Discrete Adjoint Based Time-Step Adaptation and Error Reduction in Unsteady Flow Problems , 2007 .

[7]  Edward N. Tinoco,et al.  Summary of the Fourth AIAA CFD Drag Prediction Workshop , 2010 .

[8]  R. M. Hicks,et al.  Wing Design by Numerical Optimization , 1977 .

[9]  Ralf Hartmann,et al.  An optimal order interior penalty discontinuous Galerkin discretization of the compressible Navier-Stokes equations , 2008, J. Comput. Phys..

[10]  Zhi Yang,et al.  Unstructured Dynamic Meshes with Higher-order Time Integration Schemes for the Unsteady Navier-Stokes Equations , 2005 .

[11]  A. D. Hall,et al.  The PORT Mathematical Subroutine Library , 1978, TOMS.

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

[13]  David L. Darmofal,et al.  p-Multigrid solution of high-order discontinuous Galerkin discretizations of the compressible Navier-Stokes equations , 2005 .

[14]  Patrick Knupp,et al.  Code Verification by the Method of Manufactured Solutions , 2000 .

[15]  Dimitri J. Mavriplis,et al.  A Discrete Adjoint-Based Approach for Optimization Problems on Three-Dimensional Unstructured Meshes , 2006 .

[16]  Edward N. Tinoco,et al.  Summary of Data from the Second AIAA CFD Drag Prediction Workshop (Invited) , 2004 .

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

[18]  Li Wang,et al.  Adjoint-based h-p adaptive discontinuous Galerkin methods for the 2D compressible Euler equations , 2009, J. Comput. Phys..

[19]  Dimitri J. Mavriplis,et al.  An hp-Adaptive Discontinuous Galerkin Solver for Aerodynamic flows on Mixed-Element Meshes , 2011 .

[20]  Li Wang,et al.  Implicit solution of the unsteady Euler equations for high-order accurate discontinuous Galerkin discretizations , 2006, J. Comput. Phys..

[21]  J. V. D. Vegt,et al.  Slip flow boundary conditions in discontinuous Galerkin discretizations of the Euler equations of gas dynamics , 2002 .

[22]  Jaime Peraire,et al.  Practical 3D aerodynamic design and optimization using unstructured meshes , 1996 .

[23]  Dimitri J. Mavriplis,et al.  Adjoint Sensitivity Formulation for Discontinuous Galerkin Discretizations in Unsteady Inviscid Flow Problems , 2010 .

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

[25]  Anthony Skjellum,et al.  Using MPI - portable parallel programming with the message-parsing interface , 1994 .

[26]  Chi-Wang Shu,et al.  The Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws. IV. The multidimensional case , 1990 .

[27]  S. Rebay,et al.  Numerical evaluation of two discontinuous Galerkin methods for the compressible Navier–Stokes equations , 2002 .

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

[29]  Victor W. Sparrow,et al.  Two-dimensional focusing of sonic boom noise penetrating an air-water interface , 1997 .

[30]  Ralf Hartmann,et al.  Error estimation and anisotropic mesh refinement for 3d laminar aerodynamic flow simulations , 2010, J. Comput. Phys..

[31]  S. Rebay,et al.  High-Order Accurate Discontinuous Finite Element Solution of the 2D Euler Equations , 1997 .

[32]  Derek M. Causon,et al.  On the Choice of Wavespeeds for the HLLC Riemann Solver , 1997, SIAM J. Sci. Comput..