Stability Properties of POD-Galerkin Approximations for the Compressible Navier-Stokes Equations
Résumé
Fluid flows are very often governed by the dynamics of a small number of coherent structures, i.e., fluid features which keep their individuality during the evolution of the flow. The purpose of this paper is to study a low order simulation of the Navier-Stokes equations on the basis of the evolution of such coherent structures. One way to extract from flow simulations some basis functions which can be interpreted as coherent structures is by Proper Orthogonal Decomposition (POD). Then, by means of a Galerkin projection, it is possible to find the system of ODEs which approximates the problem in the finite dimensional space spanned by the POD basis functions. It is found that the low order modeling of relatively complex flow simulations, such as the laminar vortex shedding from an airfoil at incidence and the turbulent vortex shedding from a square cylinder, provides good qualitative results compared with reference computations. In this respect, it is shown that the accuracy of numerical schemes based on simple Galerkin projection is insufficient and numerical stabilization is needed. To conclude, we approach the issue of the optimal selection of the norm used in the POD for the compressible Navier-Stokes equations by several numerical tests.