Due to the longer auto-ignition time with liquid fuels compared with hydrogen, understanding the interaction of shock waves with sprays and subsequent vapor mixing is significant for designing ramjets/scramjets with liquid fuel sprays. In this study, an Eulerian-Lagrangian framework is developed based on the OpenFOAM platform. In this solver, detailed multicomponent transport models for Eulerian gas-phase species properties are included. In addition, Lagrangian spray break-up, atomization, and evaporation models are added to simulate liquid phase, and an equilibrium wall function is added to model the near-wall properties. The newly developed solver is used to conduct large eddy simulations on nonreactive liquid jets in supersonic crossflow with liquid sprays. The liquid penetration length is compared with the experimental data, showing a very good agreement. Effects of evaporation and fuel properties (e.g., heat capacity and enthalpy of evaporation) on penetration length, temperature, Sauter mean diameter, and volumetric parcel flux are discussed in this study. It is shown that evaporation effects primarily show up in the temperature field. For n-heptane sprays, such impact could be conducted to other properties of the flow field like spray plume size, particle size distribution, and volumetric flux, which is caused by the smaller enthalpy of evaporation and heat capacity compared to water.