This paper discusses the progress made towards the development of a high-fidelity computational tool to simulate turbulent flow with shocks and particles or droplets. We have extended our recently developed high-order Eulerian-Lagrangian (EL) method to include the effects of droplet evaporation. The EL method approximates the Euler equations governing the gas dynamics with the improved high order weighted essentially non-oscillatory (WENO-Z) scheme, while individual particles are traced in the Lagrangian frame using high-order time integration schemes. A high-order ENO interpolation determines the carrier phase properties at the particle location. A high-order central weighing deposits the particle influence on the carrier phase. The droplet module is verified against the D2-law and for a uniform gas flow interaction with a rectangular cloud of droplets.