Despite great practical interest in how sprays emanate from fuel injectors, the near-nozzle region has remained a challenge for spray modelers. Recently, Eulerian models have shown promise in capturing the fast gas-liquid interactions in the near field. However, with the inclusion of compressibility, it can be difficult to maintain consistency between the hydrodynamic and thermodynamic variables. In order to resolve numerical inconsistencies that occur in segregated solutions of Eulerian spray model equations as well as to provide good scalability and stability, a new construction of a E-Y model is introduced. This construction is built around an IMEX-RK3 algorithm which offers accuracy and efficiency. The new algorithm is compared to an existing implementation for speed and is validated against experimental measurements of spray evolution in order to test the accuracy. The predictions of the new construction are slightly more accurate and, when tested on 256 processors, are 34 times faster. (C) 2016 Elsevier Ltd. All rights reserved.