International audienceJet atomization is at the core of many industrial applications such as in cryogenic combustion chambers. Since Direct Numerical Simulations (DNS) of these two-phase flows in real engines are still out of reach, reduced-order models must be built to develop predictive numerical tools. However great care must be taken on the choice of these models in order to reach sound mathematical properties and predictive simulations after a validation process.The contribution of this talk is three-fold. First, we present an Euler-Euler modelling strategy. It uses a novel hierarchy of diffuse interface models, with a proper description of various disequilibrium levels of the mixture inspired from [1,2] and a special attention devoted to the choice of these models to respect the entropy inequality. These diffuse interface models are then coupled to an element of the Kinetic-Based Moment Method (KBMM) for the dispersed flow [3].Secondly, to cope with the strong discontinuities encountered in jet atomization, a robust and accurate numerical method using multi-slope MUSCL technique is applied [4].Finally, relying on the previous two points, simulations of a jet atomization in a cryogenic combustion chamber in subcritical conditions have been conducted and results show thatphysical properties are recovered