Add assert in chem. non-equilibrium that it only works for numPhases > 1