Aerobreakup of drops is a fundamental two-phase flow problem that is essential to many spray applications. A parametric numerical study was performed by varying the gas stream velocity, focusing on the regime of moderate Weber numbers, in which the drop deforms to a forward bag. When the bag is unstable, it inflates and disintegrates into small droplets. Detailed numerical simulations were conducted using the volume-of-fluid method on an adaptive octree mesh to investigate the aerobreakup dynamics. Grid-refinement studies show that converged three-dimensional simulation results for drop deformation and bag formation are achieved by the refinement level equivalent to 512 cells across the initial drop diameter. To resolve the thin liquid sheet when the bag inflates, the mesh is refined further to 2048 cells across the initial drop diameter. The simulation results for the drop length and radius were validated against previous experiments, and good agreement was achieved. The high-resolution results of drop morphological evolution were used to identify the different phases in the aerobreakup process, and to characterize the distinct flow features and dominant mechanisms in each phase. In the early time, the drop deformation and velocity are independent of the Weber number, and a new internal-flow deformation model, which respects this asymptotic limit, has been developed. The pressure and velocity fields around the drop were shown to better understand the internal flow and interfacial instability that dictate the drop deformation. Finally, the impact of drop deformation on the drop dynamics was discussed.
more »
« less
Data-driven modeling of the aerodynamic deformation and drag for a freely moving drop in the sub-critical Weber number regime
Accurate prediction of the dynamics and deformation of freely moving drops is crucial for numerous droplet applications. When the Weber number is finite but below a critical value, the drop deviates from its spherical shape and deforms as it is accelerated by the gas stream. Since aerodynamic drag on the drop depends on its shape oscillation, accurately modeling the drop shape evolution is essential for predicting the drop's velocity and position. In this study, 2D axisymmetric interface-resolved simulations were performed to provide a comprehensive dataset for developing a data-driven model. Parametric simulations were conducted by systematically varying the drop diameter and free-stream velocity, achieving wide ranges of Weber and Reynolds numbers. The instantaneous drop shapes obtained in simulations are characterized by spherical harmonics. Temporal data of the drag and modal coefficients are collected from the simulation data to train a {Nonlinear Auto-Regressive models with eXogenous inputs} (NARX) neural network model. The overall model consists of two multi-layer perceptron networks, which predict the modal coefficients and the drop drag, respectively. The drop shape can be reconstructed with the predicted modal coefficients. The model predictions are validated against the simulation data in the testing set, showing excellent agreement for the evolutions of both the drop shape and drag.
more »
« less
- Award ID(s):
- 2321396
- PAR ID:
- 10510665
- Publisher / Repository:
- Elsevier
- Date Published:
- Journal Name:
- International Journal of Multiphase Flow
- Volume:
- 177
- ISSN:
- 0301-9322
- Page Range / eLocation ID:
- 104859
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Aerodynamic breakup of vaporizing drops is commonly seen in many spray applications. While it is well known that vaporization can modulate interfacial instabilities, the impact of vaporization on drop aerobreakup is poorly understood. Detailed interface-resolved simulations were performed to systematically study the effect of vaporization, characterized by the Stefan number, on the drop breakup and acceleration for different Weber numbers and density ratios. It is observed that the resulting asymmetric vaporization rates and strengths of Stefan flow on the windward and leeward sides of the drop hinder bag development and prevent drop breakup. The critical Weber number thus generally increases with the Stefan number. The modulation of the boundary layer also contributes to a significant increase of drag coefficient. Numerical experiments were performed to affirm that the drop volume reduction plays a negligible role and the Stefan flow is the dominant reason for the breakup suppression and drag enhancement observed.more » « less
-
null (Ed.)Aerobreakup of liquid drops are important to many droplet applications, such as fuel injection. When a liquid drop is subjected to a gas stream of high velocity, the drop can deform and break into small droplets. The drop aerobreakup is controlled by multiple dimensionless parameters. The Weber number (We) has been commonly used to characterize the different breakup regimes. While the effects of Weber and Ohnesorge numbers on the aerobreakup of a drop in unbounded domain have been extensively studied, the effect of the Reynolds number (Re) based on gas properties are less understood and will be investigated by 2D axis-symmetric and 3D detailed numerical simulations in the present paper. Attention will be focused on the moderate We regime, where the drop mostly breaks in the bag mode. In many previous studies for millimeter drops, Re is too large to be relevant. However, for applications where drops are small and the relative velocity is high, Re can be quite small when the drop breaks. Parametric simulations of Re and We are performed to systematically investigate the effect of Re on the drop aerobreakup dynamics. The simulations are performed using the Basilisk solver, where the mass-momentum consistent VOF method is used to capture the interfacial dynamics on an adaptive mesh. The reduced Re is found to induce significant changes in the drop acceleration, deformation, bag morphology, and the bag breakup dynamics, which in turn lead to significant variation in the size and spatial distributions of the children droplets formed.more » « less
-
The drag force on a spherical intruder in dense granular shear flows is studied using discrete element method simulations. Three regimes of the intruder dynamics are observed depending on the magnitude of the drag force (or the corresponding intruder velocity) and the flow inertial number: a fluctuation-dominated regime for small drag forces; a viscous regime for intermediate drag forces; and an inertial (cavity formation) regime for large drag forces. The transition from the viscous regime (linear force-velocity relation) to the inertial regime (quadratic force-velocity relation) depends further on the inertial number. Despite these distinct intruder dynamics, we find a quantitative similarity between the intruder drag in granular shear flows and the Stokesian drag on a sphere in a viscous fluid for intruder Reynolds numbers spanning five orders of magnitude. Beyond this first-order description, a modified Stokes drag model is developed that accounts for the secondary dependence of the drag coefficient on the inertial number and the intruder size and density ratios. When the drag model is coupled with a segregation force model for intruders in dense granular flows, it is possible to predict the velocity of gravity-driven segregation of an intruder particle in shear flow simulations.more » « less
-
We experimentally investigate the rise velocity of finite-sized bubbles in turbulence with a high energy dissipation rate of $$\unicode[STIX]{x1D716}\gtrsim 0.5~\text{m}^{2}~\text{s}^{-3}$$ . In contrast to a 30–40 % reduction in rise velocity previously reported in weak turbulence (the Weber number ( $We$ ) is much smaller than the Eötvös number ( $Eo$ ); $$We\ll 1more » « less
An official website of the United States government

