This paper introduces an alternative method for performing flash calculations in reacting systems, which has been developed using the theory of variable transformation of Ung and Doherty (1995). We use a simultaneous equation solving approach to perform two phase equilibrium calculations in chemically reactive systems. Thenumerical performance of proposed method is tested using several cases of study and our results show that it is efficient and generally can converge to the equilibrium compositions
