This work introduces a numerical multiphase model for multi-component mixtures, utilizing tabulated data for physical and transport properties across a spectrum of conditions from near-vacuum pressures to supercritical states. The property data are derived using Perturbed Chain Statistical Associating Fluid Theory (PC-SAFT), vapor-liquid equilibrium (VLE) calculations, entropy scaling methodologies, and Group Contribution (GC) methods. These techniques accurately reflect the thermodynamic behaviors of real fluids, avoiding the empirical estimation of Equation of State (EoS) input parameters. Implemented in OpenFOAM, the fluid dynamics solver is designed to address the three-dimensional Navier-Stokes equations for multi-component mixtures. The methodology integrates operator splitting to manage hyperbolic and parabolic steps distinctively. Hyperbolic terms are solved using the HLLC (Harten-Lax-van Leer-Contact) solver with temporal integration performed via a third-order Strong-Stability-Preserving Runge–Kutta (SSP-RK3) method. Viscous stress tensor contributions in the momentum equation are handled through an implicit velocity correction equation, while parabolic terms in the energy equation are explicitly solved. The simulation efficiency is further enhanced by adaptive Local Time Stepping and the Immersed Boundary (IB) method, which addresses interactions between the fluid and solid boundaries. Turbulence is resolved using the Wall Adaptive Large Eddy (WALE) model. Applied to high-pressure diesel fuel spray injections into non-reacting (nitrogen) gas environments, the model has been validated against Engine Combustion Network (ECN) data for the Spray-C configuration, featuring a fully cavitating multi-hole orifice. Results demonstrate that the model achieves accurate predictions across a broad range of tested conditions without the need for tuning or calibration parameters.