Modeling Heterogeneous Two-Phase Reacting Flow