Anisotropic grid adaptation for functional outputs: application to two-dimensional viscous flows

An anisotropic, unstructured grid adaptive method is presented for improving the accuracy of functional outputs of viscous, compressible flow simulations for general discretizations. The procedure merges output error control with Hessian-based anisotropic grid adaptation. An adjoint formulation is used to relate the estimated functional error to the local residual errors of both the primal and adjoint solutions. This relationship allows local error contributions to be used as indicators in a grid adaptive method designed to produce specially tuned grids for accurately estimating the chosen functional. Element stretching and orientation information is obtained from interpolation error estimates for linear triangular finite elements. The proposed adaptive method is implemented using a standard second-order upwind finite volume discretization, although the procedure is applicable to other types of discretizations such as the finite element method. A series of airfoil test cases, including separated, high-lift flows, are presented to demonstrate the approach; the functionals considered are the lift and drag coefficients. The proposed adaptive method is shown to be superior in terms of reliability and output accuracy relative to pure Hessian-based adaptation.