The standard reference is _Mixed Finite Element Methods and Applications_ by Boffi-Brezzi-Fortin. The combination $P_1-P_0$ will not, in general, satisfy the discrete _inf-sup_ condition and therefore you must either stabilize the problem or use the so-called nonconforming $P_1-P_0$ approximation where the velocity is continuous only at the triangle edge midpoints (i.e., Crouzeix-Raviart shape functions).