\frac{\partial (\rho Y _i)}{\partial t} + \nabla \cdot ( \rho \vec v Y _i) = - \nabla \cdot \vec J _i + R _i + S _i .
Ri is the net rate of production of species i by chemical reaction and Si is the rate of creation by addition from the dispersed phase and the user defined source.