Accurate and robust CFD algorithms applied to 3D arbitrary polyhedral grids

A parallel, accurate, robust and grid-transparent CFD solver was developed to solve compressible flow on 3D arbitrary polyhedral grids. To improve spacial accurate, a constrained least-squares reconstruction method is developed. To accelerate convergence, a matrix free implicit method GMRES+LU-SGS, and a parallel method using OpenMP are presented on shared-memory parallel systems with help of a special grid reordering technique. Several typical test cases, including subsonic, transonic and hypersonic flows, prove that the resulting solver performs fast, accurate and robust.