Conclusions
We have developed an efficient code for performing numerical simulations of fluid flow carrying evaporating droplets.
We have maintained a flexible code that can been easily parallelized for performing direct numerical simulations at higher Reynolds numbers and can be converted for performing large eddy simulations by the inclusion of subgrid scale models.
Our experiences illustrate how the available computational resources influence not only the selection of the numerical method, but also the physical model used.