Preparation and assembly of discrete adjoint CFD codes

A methodology for constructing the sensitivity of the incompressible Navier–Stokes equations is presented as the context for differentiating high-level Fortran source code using source-transformation automatic differentiation tools. The methodology aims to be scalable and to retain all the compile-time and run-time safety the language offers. The incompressible solver is presented as it is the standard kernel in industrial CFD software. To complement this paper, the software used for this work has been made available on www.gpde.net.