Adjoint solver development using TAPENADE

Area: CFD

An adjoint solver is developed for a second order accurate 2-D finite volume code on unstructured triangular grids which solves the inviscid, compressible flow equations. The figure shows the adjoint Euler solution corresponding to density equation and lift coefficient for NACA0012 airfoil at 3 deg. angle of attack and Mach = 0.8.
The adjoint subroutines are generated using the AD tool TAPENADE. Following Giles, Ghate, and Duta, automatic differentiation is applied at the level of elementary subroutines like the flux routines. This leads to a very efficient adjoint solver in terms of memory requirements. The full source code of the flow and adjoint solver is available for download via the URL given below.