Tube bundles in steam boilers of nuclear power plants and nuclear on-board stokehold are known to be exposed to high levels of vibrations. This coupled fluid-structure problem is very complex to numerically set up, because of its three-dimensional characteristics and because of the large number of degrees of freedom involved. A complete numerical resolution of such a problem is currently not viable, all the more so as a precise understanding of this system behaviour needs a large amount of data, obtained by very expensive calculations. We propose here to apply the now classical reduced order method called Proper Orthogonal Decomposition to a case of 2D flow around a tube bundle. Such a case is simpler than a complete steam generator tube bundle; however, it allows observing the POD projection behaviour in order to project its application on a more realistic case. The choice of POD leads to reduced calculation times and could eventually allow parametrical investigations thanks to a low data quantity. But, it implies several challenges inherent to the fluid-structure characteristic of the problem. Previous works on the dynamic analysis of steam generator tube bundles already provided interesting results in the case of quiescent fluid [J.F. Sigrist, D. Broc; Dynamic Analysis of a Steam Generator Tube Bundle with Fluid-Structure Interaction; Pressure Vessel and Piping, July 27–31, 2008, Chicago]. Within the framework of the present study, the implementation of POD in academic cases (one-dimensional equations, 2D-single tube configuration) is presented. Then, firsts POD modes for a 2D tube bundle configuration is considered; the corresponding reduced model obtained thanks to a Galerkin projection on POD modes is finally presented. The fixed case is first studied; future work will concern the fluid-structure interaction problem. Present study recalls the efficiency of the reduced model to reproduce similar problems from a unique data set for various configurations as well as the efficiency of the reduction for simple cases. Results on the velocity flow-field obtained thanks to the reduced-order model computation are encouraging for future works of fluid-structure interaction and 3D cases.