A stable and conservative high-order solver for the Reynolds-Averaged

This paper describes the development of an highly efficient parallel multiblock structured code for aerodynamic applications. The goal of our research is to assess whether or not high-order energy stable schemes are more efficient for such problems. The spatial part of the Reynolds-Averaged Navier-Stokes equations are solved making use of high-order energy stable discretization techniques based on Summation By Parts (SBP) finite difference operators and Simultaneous Approximation Term (SAT) boundary treatment [1, 2, 3, 4]. The SBP/SAT schemes we employ are up to 5th order accurate. The solver is conservative, implicit and fully coupled with a modified version of the Spalart-Allmaras turbulence model[5]. Thanks to the energy stability property of the SBP/SAT schemes, a significantly reduced amount of artificial dissipation is needed compared to schemes which do not posses this (or a similar) property. As it will be shown in the results, this leads to an higher accuracy of the numerical solutions.