skip to main content

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Thursday, February 13 until 2:00 AM ET on Friday, February 14 due to maintenance. We apologize for the inconvenience.


Title: ADVECTION AND DEPOSITION OF MICRODROPLETS IN STAGNATION POINT FLOW

We investigate trajectories of microscale evaporating droplets in a stagnation point flow near a wall of a respiratory airway. The configuration is motivated by the problem of advection and deposition of microscale droplets of respiratory fluids in human airways during transmission of infectious diseases, such as tuberculosis and COVID-19. Laminar boundary layer equations are solved to describe the airflow while the equations of motion of the droplet include contributions from gravity, aerodynamic drag, and Saffman force. Evaporation is accounted for at both the droplet surfaceand the wall of the respiratory airway and is shown to delay droplet deposition as compared to the predictions of isothermal models. Evaporation at the airway wall has a stronger effect on droplet trajectories than evaporation at the droplet surface, leading to droplets being advected away by the flow and thus avoiding deposition in the stagnation point flow region.

 
more » « less
Award ID(s):
2009741
PAR ID:
10557380
Author(s) / Creator(s):
;
Publisher / Repository:
Begell House
Date Published:
Journal Name:
Interfacial Phenomena and Heat Transfer
Volume:
11
Issue:
4
ISSN:
2169-2785
Page Range / eLocation ID:
103 to 112
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We derived equations and closed-form solutions of transit time for a viscous droplet squeezing through a small circular pore with a finite length at microscale under constant pressures. Our analyses were motivated by the vital processes of biological cells squeezing through small pores in blood vessels and sinusoids and droplets squeezing through pores in microfluidics. First, we derived ordinary differential equations (ODEs) of a droplet squeezing through a circular pore by combining Sampson flow, Poiseuille flow, and Young–Laplace equations and took into account the lubrication layer between the droplet and the pore wall. Second, for droplets wetting the wall with small surface tension, we derived the closed-form solutions of transit time. For droplets with finite surface tension, we solved the original ODEs numerically to predict the transit time. After validations against experiments and finite element simulations, we studied the effects of pressure, viscosity, pore/droplet dimensions, and surface tension on the transit time. We found that the transit time is inversely linearly proportional to pressure when the surface tension is low compared to the critical surface tension for preventing the droplet to pass and becomes nonlinear when it approaches the critical tension. Remarkably, we showed that when a fixed percentage of surface tension to critical tension is applied, the transit time is always inversely linearly proportional to pressure, and the dependence of transit time on surface tension is nonmonotonic. Our results provided a quick way of quantitative calculations of transit time for designing droplet microfluidics and understanding cells passing through constrictions.

     
    more » « less
  2. The nasopharynx, at the back of the nose, constitutes the dominant initial viral infection trigger zone along the upper respiratory tract. However, as per the standard recommended usage protocol (“Current Use”, or CU) for intranasal sprays, the nozzle should enter the nose almost vertically, resulting in sub-optimal nasopharyngeal drug deposition. Through the Large Eddy Simulation technique, this study has replicated airflow under standard breathing conditions with 15 and 30 L/min inhalation rates, passing through medical scan-based anatomically accurate human airway cavities. The small-scale airflow fluctuations were resolved through use of a sub-grid scale Kinetic Energy Transport Model. Intranasally sprayed droplet trajectories for different spray axis placement and orientation conditions were subsequently tracked via Lagrangian-based inert discrete phase simulations against the ambient inhaled airflow field. Finally, this study verified the computational projections for the upper airway drug deposition trends against representative physical experiments on sprayed delivery performed in a 3D-printed anatomic replica. The model-based exercise has revealed a new “Improved Use” (or, IU) spray usage protocol for viral infections. It entails pointing the spray bottle at a shallower angle (with an almost horizontal placement at the nostril), aiming slightly toward the cheeks. From the conically injected spray droplet simulations, we have summarily derived the following inferences: (a) droplets sized between 7–17  μ m are relatively more efficient at directly reaching the nasopharynx via inhaled transport; and (b) with realistic droplet size distributions, as found in current over-the-counter spray products, the targeted drug delivery through the IU protocol outperforms CU by a remarkable 2 orders-of-magnitude. 
    more » « less
  3. The phase transition from subcritical to supercritical conditions, referred to as transcritical behavior, significantly impacts the evaporation and fuel–air mixing in high-pressure liquid-fuel propulsion systems. Transcritical behavior is characterized as a transition from classical two-phase evaporation to a single-phase gas-like diffusion regime as surface tension and latent heat of vaporization reduce. However, the interfacial behavior represented by the surface tension coefficient and evaporation rate during this transition which are crucial inputs for Computational Fluid Dynamics (CFD) simulations of practical transcritical fuel spray is still missing. This study aims at developing new evaporation rate and surface tension models for transcritical n-dodecane droplets using molecular dynamics (MD) simulations irrespective of the droplet size. As MD simulations are primarily limited to the nanoscale, the new models can bridge the gap between MD and continuum simulations and enable the direct application of these findings to microscopic droplets. A new characteristic timescale, i.e., “undroplet time,” is defined which marks the transition from classical two-phase evaporation to single-phase gas-like diffusion behavior. The undroplet time indicates the onset of droplet core disintegration and penetration of nitrogen molecules into the droplet, which occurs after the vanishment of the surface tension. By normalizing the time with respect to the undroplet time, the rate of surface tension decay, evaporation rate, and the rate of droplet mass depletion become independent of the droplet size. Calculation of pairwise correlation coefficients for the entire MD results shows that both surface tension coefficient and evaporation rate are strongly correlated with the background temperature, while pressure and droplet size play a less significant role past the critical point. Therefore, new models for surface tension coefficient and evaporation rate spanning from sub- to supercritical conditions are developed as a function of background pressure and temperature, which can be used in continuum simulations. The identified phase change behavior based on the undroplet time shows a good agreement with the phase change regime maps obtained using microscale experiments and nanoscale MD predictions. 
    more » « less
  4. Abstract

    The significance of respiratory droplet transmission in spreading respiratory diseases such as COVID-19 has been identified by researchers. Although one cough or sneeze generates a large number of respiratory droplets, they are usually infrequent. In comparison, speaking and singing generate fewer droplets, but occur much more often, highlighting their potential as a vector for airborne transmission. However, the flow dynamics of speech and the transmission of speech droplets have not been fully investigated. To shed light on this topic, two-dimensional geometries of a vocal tract for a labiodental fricative [f] were generated based on real-time MRI of a subject during pronouncing [f]. In these models, two different curvatures were considered for the tip tongue shape and the lower lip to highlight the effects of the articulator geometries on transmission dynamics. The commercial ANSYS-Fluent CFD software was used to solve the complex expiratory speech airflow trajectories. Simultaneously, the discrete phase model of the software was used to track submicron and large size respiratory droplets exhaled during [f] utterance. The simulations were performed for high, normal, and low lung pressures to explore the influence of loud, normal, and soft utterances, respectively, on the airflow dynamics. The presented results demonstrate the variability of the airflow and droplet propagation as a function of the vocal tract geometrical characteristics and loudness.

     
    more » « less
  5. The dynamics of an evaporating droplet in an unsteady flow is of practical interest in many industrial applications and natural processes. To investigate the transport and evaporation dynamics of such droplets, we present a numerical study of an isolated droplet in an oscillating gas-phase flow. The study uses a one-way coupled two-phase flow model to assess the effect of the amplitude and the frequency of a sinusoidal external flow field on the lifetime of a multicomponent droplet containing a non-volatile solute dissolved in a volatile solvent. The results show that the evaporation process becomes faster with an increase in the amplitude or the frequency of the gas-phase oscillation. The liquid-phase transport inside the droplet also is influenced by the unsteadiness of the external gas-phase flow. A scaling analysis based on the response of the droplet under the oscillating drag force is subsequently carried out to unify the observed evaporation dynamics in the simulations under various conditions. The analysis quantifies the enhancement in the droplet velocity and Reynolds number as a function of the gas-phase oscillation parameters and predicts the effects on the evaporation rate. 
    more » « less