The assumption of phase- and chemical reaction equilibrium will introduce index problems when the dynamics of such systems are modeled in a conventional manner. By applying a modeling method which guarantees the resulting model to be semi-explicit index one, a system with phase- and chemical reaction equilibrium is modeled. Both the modeling method and the resulting model are presented. Three methods for initialization of differential algebraic equations are discussed, and one of the methods is recommended for the actual model.