The present study investigates numerically the aerodynamic breakup of Diesel droplets for a wide range of ambient pressures encountered in engineering applications relevant to oil burners and internal combustion engines. The numerical model solves the Navier-Stokes equations coupled with the Volume of Fluid (VOF) methodology utilized for capturing the interface between the liquid and the surrounding gas. An adaptive local grid refinement technique is used to increase the accuracy of the numerical results around the interface. The Weber (We) numbers examined are in the range of 14 to 279 which correspond to bag, multimode and sheet-thinning breakup regimes. Model results are initially compared against published experimental data and show a good agreement in predicting the drop deformation and the different breakup modes. The predicted breakup initiation times for all cases lie within the theoretical limits given by empirical correlations based on the We number. Following the model validation, the effect of density ratio on the breakup process is examined by varying the gas density (or equivalently the ambient pressure), while the We number is kept almost constant equal to 270; ambient gas pressure varies from 1 up to 146bar and the corresponding density ratios (ε) range from 700 down to 5. Results indicate that the predicted breakup mode of sheet-thinning remains unchanged for changing the density ratio. Useful information about the instantaneous drag coefficient (Cd) and surface area as functions of the selected non-dimensional time is given. It is shown that the density ratio is affecting the drag coefficient, in agreement with previous numerical studies.
KeywordsDroplet breakup, Diesel, VOF, density ratio, breakup initiation time.
IntroductionDroplet motion, deformation and breakup are observed in a wide variety of engineering applications such as in the injection systems of oil burners and internal combustion engines. Droplet deformation and subsequent breakup are caused by aerodynamic forces exerted on the drop by the surrounding gas, while surface tension and viscosity of the drop hinder deformation and tend to restore it to a spherical shape. The non-dimensional numbers that account for these effects are the Weber (We), Ohnesorge (Oh) and Reynolds (Re) numbers as well as the density (ε) and viscosity (N) ratios of the two phases [1]; the timescale used to non-dimensionalise the time is the shear breakup timescale [2]. For low Oh numbers (Oh<0.1) the droplet breakup is mainly controlled by the We number and according to [1] four different breakup regimes are observed as the We number increases. For We numbers in the range of 11 up to 35 the bag breakup mode is encountered during which the drop deforms into a bag resembling shape. As the We number is further increased up to 80, the multimode regime starts to appear which is essentially a combination of the bag and sheet-thinning modes. In the sheet-thinning breakup mode (80