This paper aims at presenting the DropletSMOKE++ solver, a comprehensive multidimensional computational framework for the evaporation of fuel droplets, under the influence of a gravity field and an external fluid flow. The Volume Of Fluid (VOF) methodology is adopted to dynamically track the interface, coupled with the solution of energy and species equations. The evaporation rate is directly evaluated based on the vapor concentration gradient at the phase boundary, with no need of semi-empirical evaporation sub-models.The strong surface tension forces often prevent to model small droplets evaporation, because of the presence of parasitic currents. In this work we by-pass the problem, eliminating surface tension and introducing a centripetal force toward the center of the droplet. This expedient represents a major novelty of this work, which allows to numerically hang a droplet on a fiber in normal gravity conditions without modeling surface tension. Parasitic currents are completely suppressed, allowing to accurately model the evaporation process whatever the droplet size. DropletSMOKE++ shows an excellent agreement with the experimental data in a wide range of operating conditions, for various fuels and initial droplet diameters, both in natural and forced convection. The comparison with the same cases modeled in microgravity conditions highlights the impact of an external fluid flow on the evaporation mechanism, especially at high pressures.Non-ideal thermodynamics for phase-equilibrium is included to correctly capture evaporation rates at high pressures, otherwise not well predicted by an ideal gas assumption. Finally, the presence of flow circulation in the liquid phase is discussed, as well as its influence on the internal temperature field. DropletSMOKE++ will be released as an open-source code, open to contributions from the sci-

No abstract

No abstract

No abstract

In this paper a numerical methodology for surface tension modeling is presented, with an emphasis on the implementation in the OpenFOAM R framework. The methodology relies on a combination of (i) a well-balanced approach based on the Ghost Fluid Method (GFM), including the jump of density and pressure directly in the numerical discretization of the pressure equation, and (ii) Height Functions to evaluate the interface curvature, implemented, to the authors' knowledge, for the first time in OpenFOAM R . The method is able to significantly reduce spurious currents (almost to machine accuracy) for a stationary droplet, showing second order convergence both for the curvature and the interface shape. Accurate results are also obtained for additional test cases such as translating droplets, capillary oscillations and rising bubbles, for which numerical results are comparable to what obtained by other numerical codes in the same conditions. Finally, the Height Functions method is extended to include the treatment of contact angles, both for sessile droplets and droplets suspended under the effect of gravity, showing a very good agreement with the theoretical prediction.The code works in parallel mode and details on the actual implementation in OpenFOAM R are included to facilitate the reproducibility of the results.