RocketGioud, T.Odier, N.Schmitt, T. engines may face various thermodynamic conditions, including subcritical regimes where liquid and gaseous phases coexist in a reactive, compressible environment. In this context, a multi-fluid approach, which assumes equilibrium of temperature, pressure, velocity, and Gibbs potentials, is considered to simulate these flows. Surface tension forces, which play a crucial role in atomization phenomena, are accounted for in this work. The method is evaluated on a reactive configuration representative of rocket engine conditions (test bench investigated at Technische Universität München (TUM)), presenting very promising results and highlighting the benefits of simulating the liquid jet compared to a method that models it through the injection of Lagrangian particles.