Verification of a High-Order FEM-based CFD Code using the Method of Manufactured Solutions

A high-order computational fluid dynamics (CFD) code capable of solving compressible turbulent flow problems was developed. The CFD code employs the Flowfield Dependent Variation (FDV) scheme implemented in a Finite Element Method (FEM) framework. The FDV scheme is basically derived from the Lax-Wendroff Scheme (LWS) involving the replacement of LWS’s explicit time derivatives with a weighted combination of explicit and implicit time derivatives. The code utilizes linear, quadratic and cubic isoparametric quadrilateral and hexahedral Lagrange finite elements with corresponding piecewise shape functions that have formal spatial accuracy of second-order, third-order and fourth-order, respectively. In this paper, the results of observed order-of-accuracy of the implemented FDV FEM-based CFD code involving grid and polynomial order refinements on uniform Cartesian grids are reported. The Method of Manufactured Solutions (MMS) is applied to governing 2-D Euler and Navier-Stokes equations for flow cases spanning both subsonic and supersonic flow regimes. Global discretization error analyses using discrete 𝐿2 norm show that the spatial order-of-accuracy of the FDV FEM-based CFD code converges to the shape function polynomial order plus one, in excellent agreement with theory. Uniquely, this procedure establishes the wider applicability of MMS in verifying the spatial accuracy of a high-order CFD code.