Applications of the Discrete Adjoint Method in Computational Fluid Dynamics

Abstract

The discrete adjoint method allows efficient evaluation of the derivative of a function I(s) with respect to parameters s in situations where I depends on s indirectly, via an intermediate variable w(s), which is computationally expensive to evaluate. In this thesis two applications of this method in the context of computational fluid dynamics are considered. The first is shape optimisation, where the discrete adjoint approach is employed to compute the derivatives with respect to shape parameters for a performance functional depending on the solution of a mathematical flow model which has the form of a discretised system of partial differential equations. In this context particular emphasis is given to efficient solution strategies for the linear systems arising in the discretisation of the flow models. Numerical results for two example problems are presented, demonstrating the utility of the approach. The second application, in adaptive mesh design, allows efficient evaluation of the derivatives of an a posteriori error estimate with respect to the positions of the nodes in a finite element mesh. This novel approach makes additional information available which may be utilised to guide the automatic design of adaptive meshes. Special emphasis is given to problems with anisotropic solution features, for which adaptive anisotropic mesh refinement can deliver significant performance improvements over existing adaptive hrefinement approaches. Two adaptive solution algorithms are presented and compared to existing approaches by applying them to a reaction-diffusion model problem