Search for: All records

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. 

The interfacial instability in a two-phase mixing layers between parallel gas and liquid streams is important to two-phase atomization. Depending on the inflow conditions and fluid properties, interfacial instability can be convective or absolute. The goal of the present study is to investigate the impact of gas viscosity on the interfacial instability. Both interface-resolved simulations and linear stability analysis (LSA) have been conducted. In LSA, the Orr–Sommerfeld equation is solved to analyze the spatio-temporal viscous modes. When the gas viscosity decreases, the Reynold number (Re) increases accordingly. The LSA demonstrates that when Re is higher than a critical threshold, the instability transitions from the absolute to the convective (A/C) regimes. Such a Re-induced A/C transition is also observed in the numerical simulations, though the critical Re observed in simulations is significantly lower than that predicted by LSA. The LSA results indicate that the temporal growth rate decreases with Re. When the growth rate reaches zero, the A/C transition will occur. The Re-induced A/C transition is observed in both confined and unconfined mixing layers and also in cases with low and high gas-to-liquid density ratios. In the transition from typical absolute and convective regimes, a weak absolute regime is identified in the simulations, for which the spectrograms show both the absolute and convective modes. The dominant frequency in the weak absolute regime can be influenced by the perturbation introduced at the inlet. The simulation results also show that the wave propagation speed can vary in space. In the absolute instability regime, the wave propagation speed agrees well with the absolute mode celerity near the inlet and increases to the Dimotakis speed further downstream. 

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. 

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. 

Flow-blurring atomization is an innovative twin-fluid atomization approach that has demonstrated superior effectiveness in producing fine sprays compared to traditional airblast atomization methods. In flow-blurring atomizers, the high-speed gas flow is directed perpendicular to the liquid jet. Under specific geometric and physical conditions, the gas penetrates back into the liquid nozzle, resulting in a highly unsteady bubbly two-phase mixing zone. Despite the remarkable atomization performance of flow-blurring atomizers, the underlying dynamics of the two-phase flows and breakup mechanisms within the liquid nozzle remain poorly understood, primarily due to the challenges in experimental measurements of flow details. In this study, detailed interface-resolved numerical simulations are conducted to investigate the two-phase flows generated by a planar flow-blurring atomizer. By varying key dimensionless parameters, including the dynamic-pressure ratio, density ratio, and Weber number, over wide ranges, we aim to comprehensively characterize their effects on the two-phaseflow regimes and breakup dynamics.