Numerical investigation of the aerodynamic breakup of diesel droplets under various gas pressures

[EN] Abstract 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.Financial support from the MSCA-ITN-ETN of the European Union’s H2020 programme, under REA grant agreement n. 675676 is acknowledged.Stefanitsis, D.; Malgarinos, I.; Strotos, G.; Nikolopoulos, N.; Kakaras, E.; Gavaises, M. (2017). Numerical investigation of the aerodynamic breakup of diesel droplets under various gas pressures. En Ilass Europe. 28th european conference on Liquid Atomization and Spray Systems. Editorial Universitat Politècnica de València. 1052-1059. https://doi.org/10.4995/ILASS2017.2017.4690OCS1052105

Coupling a local adaptive grid refinement technique with an interface sharpening scheme for the simulation of two-phase flow and free-surface flows using VOF methodology

This study presents the implementation of an interface sharpening scheme on the basis of the Volume of Fluid (VOF) method, as well as its application in a number of theoretical and real cases usually modelled in literature. More specifically, the solution of an additional sharpening equation along with the standard VOF model equations is proposed, offering the advantage of “restraining” interface numerical diffusion, while also keeping a quite smooth induced velocity field around the interface. This sharpening equation is solved right after volume fraction advection; however a novel method for its coupling with the momentum equation has been applied in order to save computational time. The advantages of the proposed sharpening scheme lie on the facts that a) it is mass conservative thus its application does not have a negative impact on one of the most important benefits of VOF method and b) it can be used in coarser grids as now the suppression of the numerical diffusion is grid independent. The coupling of the solved equation with an adaptive local grid refinement technique is used for further decrease of computational time, while keeping high levels of accuracy at the area of maximum interest (interface). The numerical algorithm is initially tested against two theoretical benchmark cases for interface tracking methodologies followed by its validation for the case of a free-falling water droplet accelerated by gravity, as well as the normal liquid droplet impingement onto a flat substrate. Results indicate that the coupling of the interface sharpening equation with the HRIC discretization scheme used for volume fraction flux term, not only decreases the interface numerical diffusion, but also allows the induced velocity field to be less perturbed owed to spurious velocities across the liquid–gas interface. With the use of the proposed algorithmic flow path, coarser grids can replace finer ones at the slight expense of accuracy

Numerical investigation of the aerodynamic breakup of a parallel moving droplet cluster

The present work examines numerically the aerodynamic breakup of a cluster of Diesel droplets moving in parallel with respect to the gas flow. Two- and three-dimensional simulations of the incompressible Navier-Stokes equations together with the VOF method are performed for Weber (We) numbers in the range of 5 up to 60 and non-dimensional distance between the droplets (H/D0) ranging from 1.25 to 20. The numerical results indicate that the proximity of droplets affects their breakup for distances H/D0≤5. For low droplet proximity distances (H/D0≤2.5), the droplets experience the so-called shuttlecock breakup mode, which has been also identified for droplets in tandem formations in a previous authors’ work and is characterized by an oblique peripheral stretching of the droplet. With decreasing H/D0 the breakup initiation time decreases, while the drag coefficient increases relative to that of isolated droplets. When the distance between the droplets is low enough (H/D0<1.5), this can result in critical We number, i.e. minimum We number leading to breakup, lower than that of an isolated droplet at the same conditions

Predicting droplet deformation and breakup for moderate Weber numbers

The present work examines numerically the deformation and breakup of free falling droplets subjected to a continuous cross flow. The model is based on the solution of the Navier–Stokes equations coupled with the Volume of Fluid (VOF) methodology utilized for tracking the droplet-air interface; an adaptive local grid refinement is implemented in order to decrease the required computational cost. Neglecting initially the effect of the vertical droplet motion, a 2D axisymmetric approximation is adopted to shed light on influential numerical parameters. Following that, 3D simulations are performed which include inertial, surface and gravitational forces. The model performance is assessed by comparing the results against published experimental data for the bag breakup and the sheet thinning breakup regimes. Furthermore, a parametric study reveals the model capabilities for a wider range of Weber numbers. It is proved that the model is capable of capturing qualitatively the breakup process, while the numerical parameters that best predict the experimental data are identified

Numerical investigation of aerodynamic droplet breakup in a high temperature gas environment

The Navier-Stokes equations, energy and vapor transport equations coupled with the VOF methodology and a vaporization rate model are numerically solved to predict aerodynamic droplet breakup in a high temperature gas environment. The numerical model accounts for variable properties and uses an adaptive local grid refinement technique on the gas-liquid interface to enhance the accuracy of the computations. The parameters examined include Weber (We) numbers in the range 15-90 and gas phase temperatures in the range 400-1000 K for a volatile n-heptane droplet. Initially isothermal flow conditions are examined in order to assess the effect of Weber (We) and Reynolds (Re) number. The latter was altered by varying the gas phase properties in the aforementioned temperature range. It is verified that the We number is the controlling parameter, while the Re number affects the droplet breakup at low We number conditions. The inclusion of droplet heating and evaporation mechanisms has revealed that heating effects have generally a small impact on the phenomenon due to its short duration except for low We number cases. Droplet deformation enhances heat transfer and droplet evaporation. An improved 0-D model is proposed, able to predict the droplet heating and vaporization of highly deformed droplets

Search for dark matter produced in association with bottom or top quarks in √s = 13 TeV pp collisions with the ATLAS detector

A search for weakly interacting massive particle dark matter produced in association with bottom or top quarks is presented. Final states containing third-generation quarks and miss- ing transverse momentum are considered. The analysis uses 36.1 fb−1 of proton–proton collision data recorded by the ATLAS experiment at √s = 13 TeV in 2015 and 2016. No significant excess of events above the estimated backgrounds is observed. The results are in- terpreted in the framework of simplified models of spin-0 dark-matter mediators. For colour- neutral spin-0 mediators produced in association with top quarks and decaying into a pair of dark-matter particles, mediator masses below 50 GeV are excluded assuming a dark-matter candidate mass of 1 GeV and unitary couplings. For scalar and pseudoscalar mediators produced in association with bottom quarks, the search sets limits on the production cross- section of 300 times the predicted rate for mediators with masses between 10 and 50 GeV and assuming a dark-matter mass of 1 GeV and unitary coupling. Constraints on colour- charged scalar simplified models are also presented. Assuming a dark-matter particle mass of 35 GeV, mediator particles with mass below 1.1 TeV are excluded for couplings yielding a dark-matter relic density consistent with measurements