A novel parameterized non-intrusive reduced order model (P-NIROM) based on proper orthogonal decomposition (POD) has been developed. This P-NIROM is a generic and efficient approach for model reduction of parameterized partial differential equations (P-PDEs). Over existing parameterized reduced order models (P-ROM) (most of them are based on the reduced basis method), it is non-intrusive and independent on partial differential equations and computational codes. During the training process, the Smolyak sparse grid method is used to select a set of parameters over a specific parameterized space (Ω p ∈ R P ). For each selected parameter, the reduced basis functions are generated from the snapshots derived from a run of the high fidelity model. More generally, the snapshots and basis function sets for any parameters over Ω p can be obtained using an interpolation method. The P-NIROM can then be constructed by using our recently developed technique [50,53] where either the Smolyak or radial basis function (RBF) methods are used to generate a set of hyper-surfaces representing the underlying dynamical system over the reduced space.The new P-NIROM technique has been applied to parameterized Navier-Stokes equations and implemented with an unstructured mesh finite element model. The capability of this P-NIROM has been illustrated numerically by two test cases: flow past a cylinder and lock exchange case. The prediction capabilities of the P-NIROM have been evaluated by varying the viscosity, initial and boundary conditions. The results show that this P-NIROM has captured the quasi-totality of the details of the flow with CPU speedup of three orders of magnitude. An error analysis for the P-NIROM has been carried out.