Computation Fluid Dynamics (CFD) simulation has become a routine design tool for i) predicting accurately the thermal performances of electronics set ups and devices such as cooling system and ii) optimizing configurations. Although CFD simulations using discretization methods such as finite volume or finite element can be performed at different scales, from component/board levels to larger system, these classical discretization techniques can prove to be too costly and time consuming, especially in the case of optimization purposes where similar systems, with different design parameters have to be solved sequentially. The design parameters can be of geometric nature or related to the boundary conditions. This motivates our interest on model reduction and particularly on reduced basis methods. As is well documented in the literature, the offline/online implementation of the standard RB method (a Galerkin approach within the reduced basis space) requires to modify the original CFD calculation code, which for a commercial one may be problematic even impossible. For this reason, we have proposed in a previous paper, with an application to a simple scalar convection diffusion problem, an alternative non-intrusive reduced basis approach (NIRB) based on a two-grid finite element discretization. Here also the process is two stages: offline, the construction of the reduced basis is performed on a fine mesh; online a new configuration is simulated using a coarse mesh. While such a coarse solution, can be computed quickly enough to be used in a rapid decision process, it is generally not accurate enough for practical use. In order to retrieve accuracy, we first project every such coarse solution into the reduced space, and then further improve them via a rectification technique. The purpose of this paper is to generalize the approach to a CFD configuration.