A new Σ-Y Eulerian model construction has been built around an IMEX-RK3 algorithm.
A density-based algorithm is used to resolve inconsistencies of segregated solutions.
The new model and the existing one have been compared on diesel sprays simulations.
New model improves the accuracy, stability and parallel efficiency of simulations.
New Runge-Kutta construction has long-term value for its potential.