skip to main content


Title: A Dual Scale Model for Reconstructing Sub-Filter Shear-Induced Instabilities Using a Vortex Sheet Method
A method to predict sub-filter shear-induced velocities on a liquid-gas phase interface for use in a dual scale LES model is presented. The method reconstructs the sub-filter velocity field in the vicinity of the interface by introducing a vortex sheet at the interface. The vortex sheet is transported by an unsplit geometric volume and surface area advection scheme with a Piece- wise Linear Interface Construction (PLIC) representation of the material interface. At each step and desired location the shear-induced velocities can be calculated by integrating the vortex sheet and other relevant quantities over the liquid-gas surface with the sub-grid velocity recon- struction limited to a small number of cells near the phase interface. The vortex sheet method is tested and compared against prior literature.  more » « less
Award ID(s):
1803657
NSF-PAR ID:
10281849
Author(s) / Creator(s):
;
Date Published:
Journal Name:
ICLASS 2021, 15th Triennial International Conference on Liquid Atomization and Spray Systems
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    A method to predict sub-filter shear-induced velocities on a liquid-gas phase interface for use in a dual scale LES model is presented. The method reconstructs the sub-filter velocity field in the vicinity of the interface by introducing a vortex sheet at the interface. The vortex sheet is transported by an unsplit geometric volume and surface area advection scheme with a Piecewise Linear Interface Construction (PLIC) representation of the material interface. At each step and desired location the shear-induced velocities can be calculated by integrating the vortex sheet and other relevant quantities over the liquid-gas interface with the sub-grid velocity reconstruction limited to a small number of cells near the phase interface. The vortex sheet method is tested and compared against prior literature. 
    more » « less
  2. A method to predict sub-filter shear-induced velocities on a liquid-gas phase interface for use in a dual scale LES model is presented and compared against prior work on Vortex Sheet methods. The method reconstructs the sub-filter velocity field in the vicinity of the interface by employing a vortex sheet at the interface location. The vortex sheet is transported by an unsplit geometric volume and surface area advection scheme with a Piecewise Linear Interface Construction (PLIC) representation of the material interface. At each step, the vorticity field is constructed by evaluating a volume integral of the vortex sheet and a numerical spreading parameter near the liquid-gas interface. A Poisson equation can then be constructed and solved for the vector potential; the self-induced velocities due to the vortex sheet are subsequently evaluated from the vector potential. The described vortex sheet method is tested and compared against prior literature. 
    more » « less
  3. null (Ed.)
    A method to compute sub-filter velocities due to shear induced instabilities on a liquid-gas interface for use in a dual scale LES-DNS model is presented. The method reconstructs the sub-filter velocity field as the sum of a prescribed base velocity profile and a perturbation velocity field determined by the Orr-Sommerfeld equations. The base velocity profile is approximated as an error function appropriately scaled with flow parameters, and the perturbation velocity field is computed by solving the Orr-Sommerfeld equations with appropriate boundary and interface conditions. The perturbation velocities of the Orr-Sommerfeld equations are expanded into Chebyshev polynomials to create a linear eigenvalue problem as outlined by Schmid and Henningson (2001). Finally the eigenvalue problem is solved using a standard linear algebra package and used to evaluate the perturbation velocities. The Chebyshev method is tested under a variety of flow parameters and initial interface disturbances. Results are presented and compared against prior literature and asymptotic solutions. 
    more » « less
  4. Understanding the development and breakup of interfacial waves in a two-phase mixing layer between the gas and liquid streams is paramount to atomization. Due to the velocity difference between the two streams, the shear on the interface triggers a longitudinal instability, which develops to interfacial waves that propagate downstream. As the interfacial waves grow spatially, transverse modulations arise, turning the interfacial waves from quasi-two-dimensional to fully three-dimensional. The inlet gas turbulence intensity has a strong impact on the interfacial instability. Therefore, parametric direct numerical simulations are performed in the present study to systematically investigate the effect of the inlet gas turbulence on the formation, development and breakup of the interfacial waves. The open-source multiphase flow solver, PARIS, is used for the simulations and the mass–momentum consistent volume-of-fluid method is used to capture the sharp gas–liquid interfaces. Two computational domain widths are considered and the wide domain will allow a detailed study of the transverse development of the interfacial waves. The dominant frequency and spatial growth rate of the longitudinal instability are found to increase with the inlet gas turbulence intensity. The dominant transverse wavenumber, determined by the Rayleigh–Taylor instability, scales with the longitudinal frequency, so it also increases with the inlet gas turbulence intensity. The holes formed in the liquid sheet are important to the disintegration of the interfacial waves. The hole formation is influenced by the inlet gas turbulence. As a result, the sheet breakup dynamics and the statistics of the droplets formed also change accordingly. 
    more » « less
  5. The evaporation of droplets on surfaces is a ubiquitous phenomenon essential in nature and industrial applications ranging from thermal management of electronics to self-assembly-based fabrication. In this study, water droplet evaporation on a thin quartz substrate is analyzed using an unsteady two-step arbitrary Lagrangian-Eulerian (ALE) moving mesh model, wherein the evaporation process is simulated during the constant contact radius (CCR) and contact angle (CCA) modes. The numerical model considers mass transfer in the gas domain, flow in the liquid and gas domains, and heat transfer in the solid, liquid, and gas domains. Besides, the model also accounts for interfacial force balance, including thermocapillary stresses, to obtain the instantaneous droplet shape. Experiments involving droplet evaporation on unheated quartz substrates agree with model predictions of contact radius, contact angle, and droplet volume. Model results indicating temperature and velocity distribution across an evaporating water droplet show that the lowest temperatures are at the liquid-gas interface, and a single vortex exists for the predominant duration of the droplet's lifetime. The temperature of the unheated substrate is also significantly reduced due to evaporative cooling. The interfacial evaporation flux distribution, which depends on heat transfer across the droplet and advection in the surrounding medium, shows the highest values near the three-phase contact line. In addition, the model also predicts evaporation dynamics when the substrate is heated and exposed to different advection conditions. Generally, higher evaporation rates result from higher substrate heating and advection rates. However, substrate heating and advection in the surrounding gas have minimal effects on the relative durations of CCR and CCA modes for a given receding contact angle. Specifically, in this case, a 40× increase in substrate heating rate or 7.5× increase in gas velocity can only change these relative durations by 3%. This study also highlights the importance of surface wettability, which affects evaporation dynamics for all the conditions explored by the numerical model. 
    more » « less