Combining phase-space and time-dependent reduced density matrix approach to describe the dynamics of interacting fermions