Reformulate diffusion laws to match with mass reference velocities

As discussed previously: if Darcys law and Navier Stokes law give mass averaged velocities we need to adapt the diffusion laws to calculate the gradient based on mass fractions not mole fractions.