Proper orthogonal decomposition-based reduced-order stabilized mixed finite volume element extrapolating model for the nonstationary incompressible Boussinesq equations