A vector penalty-projection approach for the time-dependent incompressible magnetohydrodynamics flows