A combined BDF-semismooth Newton approach for time-dependent Bingham flow