In the Non-equilibrium DPD-MD simulation, we obtained the trajectory of a system as a function of time in an NVE ensemble. The system traverses through multiple microstates before it finally reaches an equilibrium. While an equilibrium DPD-MD simulation was performed with a widely used Nose-Hover thermostat, this does not have the capability to ascertain the path of the system that reaches the final state when the system is coupled with the hydrodynamic forces. We used DPD-MD simulations with a non-isothermal algorithm to understand the mechanism of emulsification. This NVE ensemble fails to converge within the scope of the first-fluctuation dissipation theorem. Hence these simulations were performed with a modified DPD algorithm (using Second fluctuation dissipation theorem allowing the kinetic temperature to evolve). A velocity verlet method is also adopted to these simulations in order to prevent the self-penetration of the molecules within each other.