In this paper, a proper orthogonal decomposition based reduced-order model is presented for parametrized multiphysics computations. Our application physics is Doppler feedback in a simplified model of the molten salt fast reactor concept. The reduced model is created using the method of snapshots where the offline training set is obtained by exercising a full-order model created with the OpenFOAM based multiphysics solver, GeN-Foam. The steady state models solve the multi-group diffusion k-eigenvalue equations with moving precursors together with the energy equation. A fixed velocity field is assumed throughout the computations, hence the momentum and continuity equations are not solved. The discrete empirical interpolation method is used for the efficient coupling of the ROM solvers, while the input parameter space is surveyed using the improved distributed latin hypercube sampling algorithm.