Mixing under high pressure conditions plays a central role in several engineering applications, such as direct-injection engines and liquid rocket engines. Numerical flow simulations have become a complementary tool to study the mixing process under these conditions but require complex thermodynamic modeling as well as validation with accurate experimental data. For this reason, we use experiments of supercritical single-phase jet mixing from the literature, where the mixing is quantified by the mixture speed of sound, as a reference for our work. We here focus on the thermodynamic modeling of multi-component flows under high pressure conditions and the analytical calculation of the mixture speed of sound. Our thermodynamic model is based on cubic equations of state extended for multi-components. Using an extension of OpenFOAM, we perform large-eddy simulations of hexane and pentane injections and compare our results with the experimentally measured mixture speed of sound at specific positions. The simulation results show the same characteristic trends, indicating that the mixing effects are well reproduced in the simulations. Additionally, the effect of the sub-grid scale modeling is assessed by comparing results using different models (Smagorinsky, Vreman, and Wall-Adapting Local Eddy-viscosity). The comprehensive simulation data presented here, in combination with the experimental data, provide a benchmark for numerical simulations of jet mixing in high pressure conditions.