Analysis of finite-volume discrete adjoint fields for two-dimensional compressible Euler flows

This work deals with a number of questions relative to the discrete and continuous adjoint fields associated to Euler equations and classical aerodynamic functions. The consistency of the discrete adjoint equations with the corresponding continuous adjoint partial differential equation is one of them. It is has been established or at least discussed only for a handful of numerical schemes and a contribution of this article is to give the adjoint consistency conditions for the 2D Jameson-Schmidt-Turkel scheme in cell-centred finite-volume formulation. The consistency issue is also studied here from a new heuristic point of view by discretizing the continuous adjoint equation for the discrete flow and adjoint fields. Both points of view prove to provide useful information. Besides, it has been often noted that discrete or continuous inviscid lift and drag adjoint exhibit numerical divergence close to the wall and stagnation streamline for a wide range of subsonic and transonic flow conditions. This is analyzed here using the physical source term perturbation method introduced in reference [Giles and Pierce, AIAA Paper 97-1850, 1997]. With this point of view, the fourth physical source term of appears to be the only one responsible for this behavior. It is also demonstrated that the numerical divergence of the adjoint variables corresponds to the response of the flow to the convected increment of stagnation pressure and diminution of entropy created at the source and the resulting change in lift and drag.

[1]  Stefan Ulbrich,et al.  A Sensitivity and Adjoint Calculus for Discontinuous Solutions of Hyperbolic Conservation Laws with Source Terms , 2002, SIAM J. Control. Optim..

[2]  G. R. Shubin,et al.  A comparison of optimization-based approaches for a model computational aerodynamics design problem , 1992 .

[3]  Jason E. Hicken,et al.  Dual consistency and functional accuracy: a finite-difference perspective , 2014, J. Comput. Phys..

[4]  Frédéric Alauzet,et al.  An a priori anisotropic goal-oriented error estimate for viscous compressible flow and application to mesh adaptation , 2019, J. Comput. Phys..

[5]  Arthur Stück,et al.  Dual-consistency study for Green-Gauss gradient schemes in an unstructured Navier-Stokes method , 2017, J. Comput. Phys..

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

[7]  Ralf Hartmann,et al.  Adjoint Consistency Analysis of Discontinuous Galerkin Discretizations , 2007, SIAM J. Numer. Anal..

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

[9]  J. Lions,et al.  Contrôle optimal de systèmes gouvernés par des équations aux dérivées partielles , 1968 .

[10]  J. Eric,et al.  Aerodynamic Design Optimization on Unstructured Meshes Using the Navier-Stokes Equations , 1998 .

[11]  A. Majda,et al.  The stability of multidimensional shock fronts , 1983 .

[12]  Enrique Zuazua,et al.  Systematic Continuous Adjoint Approach to Viscous Aerodynamic Design on Unstructured Grids , 2007 .

[13]  Enrique Zuazua,et al.  2-D Euler Shape Design on Nonregular Flows Using Adjoint Rankine-Hugoniot Relations , 2008 .

[14]  John C. Vassberg,et al.  In Pursuit of Grid Convergence for Two-Dimensional Euler Solutions , 2009 .

[15]  Michael B. Giles,et al.  Analytic adjoint solutions for the quasi-one-dimensional Euler equations , 2001, Journal of Fluid Mechanics.

[16]  Adrian Sandu,et al.  On the properties of discrete adjoints of numerical methods for the advection equation , 2008 .

[17]  S. Hiernaux,et al.  Aerodynamic optimization using Navier-Stokes equations and optimal control theory , 1999 .

[18]  Maxime Nguyen-Dinh,et al.  Qualification des simulations numériques par adaptation anisotropique de maillages , 2014 .

[19]  Carlos Lozano,et al.  Watch Your Adjoints! Lack of Mesh Convergence in Inviscid Adjoint Solutions , 2019, AIAA Journal.

[20]  Marcus Meyer,et al.  Stabilisation of discrete steady adjoint solvers , 2015, J. Comput. Phys..

[21]  On the Properties of Solutions of the 2D Adjoint Euler Equations , 2018, Computational Methods in Applied Sciences.

[22]  R. Dwight,et al.  Effect of Approximations of the Discrete Adjoint on Gradient-Based Optimization , 2006 .

[23]  D. Zingg,et al.  Newton-Krylov Algorithm for Aerodynamic Design Using the Navier-Stokes Equations , 2002 .

[24]  Antony Jameson,et al.  Optimum aerodynamic design using the Navier-Stokes equations , 1997 .

[25]  P. Courtier,et al.  Variational Assimilation of Meteorological Observations With the Adjoint Vorticity Equation. I: Theory , 2007 .

[26]  Michael B. Giles,et al.  Improved- lift and drag estimates using adjoint Euler equations , 1999 .

[27]  C. Lozano.  A note on the dual consistency of the discrete adjoint quasi-one-dimensional Euler equations with cell-centered and cell-vertex central discretizations , 2016 .

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

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

[30]  Andrew J. Majda,et al.  The validity of the modified equation for nonlinear shock waves , 1985 .

[31]  O. Pironneau.  On optimum design in fluid mechanics , 1974 .

[32]  P. Luchini,et al.  Adjoint Equations in Stability Analysis , 2014, 2404.17304.

[33]  David L. Darmofal,et al.  Adaptive precision methodology for flow optimization via discretization and iteration error control , 2004 .

[34]  On the adjoint solution of the quasi‐1D Euler equations: the effect of boundary conditions and the numerical flux function , 2005 .

[35]  Carlos Lozano,et al.  Singular and Discontinuous Solutions of the Adjoint Euler Equations , 2018, AIAA Journal.

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

[37]  C. Hirsch,et al.  Numerical Computation of Internal and External Flows. By C. HIRSCH. Wiley. Vol. 1, Fundamentals of Numerical Discretization. 1988. 515 pp. £60. Vol. 2, Computational Methods for Inviscid and Viscous Flows. 1990, 691 pp. £65. , 1991, Journal of Fluid Mechanics.

[38]  Arthur Stück,et al.  An adjoint view on flux consistency and strong wall boundary conditions to the Navier-Stokes equations , 2015, J. Comput. Phys..

[39]  Jacques Peter,et al.  Goal oriented mesh adaptation using total derivative of aerodynamic functions with respect to mesh coordinates – With applications to Euler flows , 2012 .

[40]  Richard P. Dwight,et al.  Heuristic a posteriori estimation of error due to dissipation in finite volume schemes and application to mesh adaptation , 2008, J. Comput. Phys..

[41]  T. J. Mitchell,et al.  Bayesian design and analysis of computer experiments: Use of derivatives in surface prediction , 1993 .

[42]  A. Jameson,et al.  Numerical solution of the Euler equations by finite volume methods using Runge Kutta time stepping schemes , 1981 .

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

[44]  Denis Sipp,et al.  Stability, Receptivity, and Sensitivity Analyses of Buffeting Transonic Flow over a Profile , 2015 .

[45]  Rolf Rannacher,et al.  An optimal control approach to a posteriori error estimation in finite element methods , 2001, Acta Numerica.

[46]  Lei Shi,et al.  Adjoint-based error estimation and mesh adaptation for the correction procedure via reconstruction method , 2015, J. Comput. Phys..

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

[48]  R. Dwight,et al.  Numerical sensitivity analysis for aerodynamic optimization: A survey of approaches , 2010 .

[49]  J. Essers,et al.  An Optimal Control Theory Based Algorithm to Solve 2D Aerodynamic Shape Optimisation Problems for Inviscid and Viscous Flows , 2000 .