Special issue in honor of Prof. David Y.H. Pui for his “50 Years of Contribution in Aerosol Science and Technology” (III)

Dong-Bin Kwak1, Hayden D. Fischer1, David Y.H. Pui This email address is being protected from spambots. You need JavaScript enabled to view it.1,2

1 Particle Technology Laboratory, Mechanical Engineering, University of Minnesota, Minneapolis 55455, USA
2 School of Science and Engineering, The Chinese University of Hong Kong, Shenzhen, Guangdong 518172, China


Received: November 21, 2022
Revised: January 9, 2023
Accepted: January 18, 2023

 Copyright The Author(s). This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are cited.


Download Citation: ||https://doi.org/10.4209/aaqr.220409  


Cite this article:

Kwak, D.B., Fischer, H.D., Pui, D.Y.H. (2023). Saliva Droplet Evaporation Experiment and Simple Correlation of Evaporation-falling Curve under Different Temperatures and RH. Aerosol Air Qual. Res. 23, 220409. https://doi.org/10.4209/aaqr.220409


HIGHLIGHTS

  • Ratio of saliva residue to initial droplet size is fixed to be 0.216.
  • Saliva correction factor for evaporation coefficient was determined to be 0.857.
  • Lifetime of saliva droplets is expected to be longer than water droplets.
  • Simple equation for calculating the evaporation-falling curve was presented.
 

ABSTRACT


In the present study, water and saliva droplet evaporation rates are experimentally evaluated in various environmental conditions: temperature and relative humidity. We found that the ratio of saliva residue size to saliva initial droplet size is 0.216. We also found that the saliva correction factor which is defined as a ratio of water evaporation coefficient to saliva evaporation coefficient was not affected by environmental conditions and was determined as 0.857. By considering the evaporation and gravitational settling effects, the saliva airborne lifetime, i.e., the evaporation-falling curve is numerically calculated and characterized. In addition to this numerical calculation results, we present the simple correlation equation for obtaining the evaporation-falling curve under different temperatures and relative humidity.


Keywords: SARS-CoV-2, Evaporation-falling curve, Saliva evaporation


1 INTRODUCTION


Social distancing and punitive measures to prevent the spread of COVID-19 through human-to-human transmission are having a huge impact on the environment, health, and economy (Sarkodie and Owusu, 2021). Despite the current growing vaccination rates, the continuous emergence of newly mutated viruses such as Beta, Delta, and Omicron variants makes it difficult to prevent virus transmissions. As of January 9, 2023, there were more than 657 million confirmed cases and 6.6 million reported deaths worldwide due to COVID-19 (WHO, 2023). In addition, the recent study (Liu et al., 2021) has shown that omicron variants can evade the immune protection conferred by vaccines and natural infections. Therefore, as the current serious COVID-19 pandemic is getting worse, not only social measures like social distancing and border restrictions but also academic research to mitigate the risk of SARS-CoV-2 virus transmission is still required.

Centers for Disease Control and Prevention (CDC) claims that there are three main routes for people to become infected with SARS-CoV-2 by infectious virus-laden respiratory droplets (CDC 2021): (1) inhalation of fine respiratory droplets or aerosol particles, i.e., airborne transmission; (2) deposition of respiratory droplets in the mouth, nose, or eyes, i.e., droplet transmission; (3) touching mucous membranes with dirty hands (direct contact with respiratory fluid or indirect touch viral contaminated surfaces), i.e., fomite transmission. Recently, to mitigate SARS-CoV-2 transmission, various methods like room ventilation (Dao and Kim, 2022; Faulkner et al., 2022; Cao et al., 2021; Chen et al., 2021), surface decontamination (Kim et al., 2021; Kwak et al., 2021a), and wearing facial masks (Fakir and Kim, 2022; Ou et al., 2020) are suggested. Even though there is a controversial discussion regarding the dominant SARS-CoV-2 transmission route, airborne or droplet transmission plays a vital role in the spread of diseases (Cheng et al., 2022; Lai et al., 2020; Wan et al., 2009). People can produce a considerable amount of respiratory droplets like saliva during exhalation, e.g., breathing, talking, singing, coughing, and sneezing (Edwards et al., 2021; Alsved et al., 2020; Asadi et al., 2019, 2020; Echternach et al., 2020; Stadnytskyi et al., 2020; Morawska et al., 2009; Xie et al., 2009; Yang et al., 2007; Duguid, 1946). For determining the risk of airborne or droplet transmission of SARS-CoV-2, the airborne lifetime of respiratory droplets should be estimated.

Various methods were employed for determining saliva droplet size distribution, even though there are significant discrepancies between previous studies. By using an aerodynamic particle sizer, Alsved et al. (2020) measured exhaled respiratory droplets during the singing and talking and revealed that most of the droplets are 0.5−0.7 µm. By using the solid impaction method, Duguid (1946) reported that 95% of respiratory droplets were produced between 2 µm and 100 µm when coughing or sneezing and Xie et al. (2009) reported that the average size of saliva droplets during coughing and sneezing was approximately 50−100 µm. By using an aerodynamic particle sizer and scanning mobility particle sizer, Yang et al. (2007) reported that respiratory droplets ranging from 0.62 to 15.9 µm with an average mode of 8.4 µm were produced by coughing. Han et al. (2013) investigated particle size distribution of the droplets exhaled by sneezing using a laser particle size analyzer and found two different distributions from twenty subjects: the geometric mean sizes for unimodal distribution, and bimodal distribution were 360.1 µm, and 74.4 µm, respectively. By using the interferometric Mie imaging technique, Chao et al. (2009) reported that geometric mean sizes of respiratory droplets were 13.5 µm for coughing, and 16 µm for talking.

Respiratory droplet residue is the final product of a respiratory droplet after complete evaporation produced by the atomization of human secretions along the airway (Liu et al., 2017; Hare, 1964). Numerous studies on the evaporation of water have previously been conducted. Al-Shammiri (2002) investigated the effects of salinity on water evaporation by measuring the evaporation rate of distilled water, sea water (34,000 ppm), and two brine solutions of 56,000 and 69,000 ppm. Furthermore, Zhang (2011) also explored the correlation between salinity and evaporation as well as the effects of varied relative humidity and temperature conditions on evaporation. Both confirmed that an increase in salinity will result in a decrease in the evaporation rate. Zhang (2011) also determined that altering relative humidity with a constant temperature will be inversely related while altering temperature with constant relative humidity and the water evaporation rate provides an inconclusive relationship. However, Liu et al. (2017) found through a cough jet that medium-sized droplets, around 60 µm, are more influenced by humidity than smaller and larger droplets. A study by Hisatake et al. (1993) explored the evaporation rate of water within a vessel, undergoing varied vessel diameter, relative humidity, temperature, and air velocity. By comparing the theory of water vapor evaporation rate through the boundary layer covering the water surface in a vessel to experimental results, Hisatake et al. (1993) found that the theory provides good estimates. Other variables of environmental conditions have also been analyzed.

In this study, saliva evaporation characteristics, which are very important for understanding the airborne transmission of SARS-CoV-2, were experimentally observed under various ambient conditions (20% ≤ RH ≤ 70%, 15°C ≤ T ≤ 30°C). To depict saliva in the air, a microgravity environment was implemented using an acoustic levitator. By comparing the evaporation characteristics of saliva obtained through the experiment with that of water, the data were quantified in two aspects: (1) residue size obtained after complete evaporation of saliva; (2) ratio of saliva evaporation coefficient. Using the results obtained through the experiment, the evaporation-falling curve was obtained through numerical analysis, and a simple correlation curve that can obtain the evaporation-falling curve without numerical analysis was presented. The simple correlation of the evaporation-falling curve is expected to be utilized to determine the risk of airborne virus transmission in various atmospheric conditions without complicated calculations.

 
2 EXPERIMENTAL METHOD


The experimental setup used to measure the change of droplet sizes during evaporation was illustrated in Fig. 1. The acoustic levitator (Tec5 AG, Oberursel, Germany) was used for floating the droplets in the air without contacting the surface. Droplets can be suspended in a standing wave which is generated from caused by multiple reflections between a piezoelectric transducer and a concave reflector, resulting from the balance of gravity and forces caused by axial acoustic pressure. In our tested condition, droplets with diameters of 100−1000 µm can be suspended. The acoustic levitator was placed in a controlled environmental chamber (Environmental Growth Chamber, GCW15, Chargin Falls, OH, USA). The evaporation process of acoustically levitated saliva or water was recorded by a microscope camera with a 300X industry zoom c-mount lens (HY-2307K, Hayear Electronics, Shenzhen, China). The temperature and relative humidity inside the environmental chamber changed from 15°C to 30°C and from 20% to 70%, respectively, so that a total of 28 conditions were tested.

Fig. 1. (a) Environmental chamber to control temperature and relative humidity and (b) experimental setup to measure changes in the size of acoustically levitated water or saliva droplets during evaporation.Fig. 1. (a) Environmental chamber to control temperature and relative humidity and (b) experimental setup to measure changes in the size of acoustically levitated water or saliva droplets during evaporation.

To obtain the size evolution during the droplet evaporation, the images were acquired every second. For handling lots of images obtained from the microscope camera, the images were processed with a certain rule by using ImageJ 1.53e software (Schneider et al., 2012), as shown in Fig. 2(a). The original images of Fig. 2 show a distinct oval shape in the middle of the levitated droplets due to the light refracted from the vertically installed LED light in Fig. 1(b). For finding the edges of droplets, the Sobel edge detector algorithm which measures 2D spatial gradient on images was implemented, resulting in distinguishing the droplets and background noise. For calculating the droplet size based on the droplet projected area, RGB color images need to be binary images. Through averaging each RGB value of images, 8-bit grayscale images were obtained.

 Fig. 2. (a) Flow chart of the image processing algorithm and (b) examples of obtained binary image and droplet size through the image processing algorithm. The scale bar indicates 100 µm.Fig. 2. (a) Flow chart of the image processing algorithm and (b) examples of obtained binary image and droplet size through the image processing algorithm. The scale bar indicates 100 µm.

By comparing 16 different auto thresholds in ImageJ software, the triangle method (Zack et al., 1977) was chosen for converting 8-bit grayscale images to binary images. Based on the outline of binary images, the images are filled with white color. The droplet sizes were determined by the projected area of binary images. Fig. 2(b) shows examples of obtained binary images and corresponding droplet diameters. Some images were not analyzed when blur images were obtained due to shaking droplets as presented in Fig. S2 in the Supplementary Information.

 
3 EVAPORATION MODELING


For smaller Knudsen numbers, the evaporation depends on the diffusion rate of molecules on the droplet surface. The rate of evaporation can be expressed as follows (Niimura and Hasegawa 2019; Frohn and Roth, 2000; Hinds, 1982):

 

where dp is droplet diameter, D is diffusion coefficient, M is the molecular weight, R is the gas constant of 8.3145 J mol–1 K1, ρp is the droplet density, pv, is the partial vapor pressure remote from the droplet surface, i.e. ambient partial vapor pressure, pv,s is the partial vapor pressure at the droplet surface, T is ambient temperature, Ts is droplet surface temperature, and ϕ is the correction factor due to experimental conditions. To compensate for the difference between water and saliva evaporation, the saliva correction factor (ϕsaliva) was introduced. In the present study, since the test substances were water and saliva (which contains water 99.5%) (Radha et al., 2016; Schipper et al., 2007), several physical parameters can be estimated (Huang, 2018; Kukkonen et al., 1989; Hinds, 1982):

 

and

 

Here, pamb is ambient pressure, pv,sat is saturation vapor pressure, and SR is the saturation ratio. The correction factor for the acoustic levitation experiment (ϕ) is suggested as a function of the Sherwood number (Schiffter and Lee, 2007; Tuckermann et al., 2002). It should be noticed that the correction factor for the acoustic levitation experiment proportionally increases with the intensity of the ultrasonic field (Tuckermann et al., 2002). For our conditions, the correction factor was determined to be 1.1552, as presented in Fig. S2 in the Supplementary Information. Through integrating Eq. (1) with respect to time, the rate of evaporation which is not influenced by the different initial droplet sizes can be expressed as follows:

 

To simplify the dimensionless Eq. (7), the evaporation coefficient (βv) was introduced. The evaporation coefficient (βv) and the simplified dimensionless Eq. (7) are as follows:

 

and

 

If the vapor pressures and temperatures slowly change during the evaporation experiment, the evaporation coefficient is nearly constant in time. As shown in Fig. 3, it can be seen that dimensionless Eq. (7) (or Eq. (9)) is not affected by different initial droplet sizes under the same temperature and relative humidity conditions, which characteristic is known as the d2-law. It should be noted that the slope in Fig. 3(b) indicates the corrected evaporation coefficient (βvϕ). Thus, the simplified dimensionless equation (i.e., Eq. (4)) allows a direct comparison of the water and saliva droplet evaporation process with different initial sizes.

Fig. 3. Results of water droplet evaporation experiment at T = 20°C and RH = 25%. (a) Droplet diameter as a function of time for five droplets with various initial droplet diameters. (b) Evolution of (dp2/dp,o2) as a function of 4t/dp,o2 with five different initial droplet diameters. T, temperature; RH, relative humidity; dp, droplet diameter; dp,o, initial droplet diameter; t, time.Fig. 3. Results of water droplet evaporation experiment at T = 20°C and RH = 25%. (a) Droplet diameter as a function of time for five droplets with various initial droplet diameters. (b) Evolution of (dp2/dp,o2) as a function of 4t/dp,o2 with five different initial droplet diameters. T, temperature; RH, relative humidity; dp, droplet diameter; dp,o, initial droplet diameter; t, time.

 
4 RESULTS AND DISCUSSION


Fig. 4 shows the evaporation process of acoustically levitated pure water and saliva droplets at a temperature of 20°C and a relative humidity of 50%. Because the evaporation process takes place at the droplet surface, both water and saliva droplets evaporate at the rate proportional to the surface area, resulting in a constant slope (i.e., evaporation coefficient) regardless of time, as shown in Fig. 4(b). However, two distinct evaporation characteristics of saliva compared to pure water were confirmed: (1) before about 600 seconds, the slope in Fig. 4(b) was slightly less steep than that of pure water; (2) after about 600 seconds, saliva did not evaporate any more, unlike pure water, and created a residue particle. Therefore, this study describes the characteristics of saliva residue in Section 4.1 and compares water and saliva evaporation coefficients in Section 4.2.

Fig. 4. Comparison of water and saliva droplet evaporation experiment results at T = 20°C and RH = 50%. (a) Droplet diameter as a function of time for water and saliva droplets. Inserted images are the example original images during the saliva droplet evaporation. (b) Evolution of (dp2/dp,o2) as a function of 4t/dp,o2 for water and saliva droplets. T, temperature; RH, relative humidity; dp, droplet diameter; dp,o, initial droplet diameter; t, time.Fig. 4. Comparison of water and saliva droplet evaporation experiment results at T = 20°C and RH = 50%. (a) Droplet diameter as a function of time for water and saliva droplets. Inserted images are the example original images during the saliva droplet evaporation. (b) Evolution of (dp2/dp,o2) as a function of 4t/dp,o2 for water and saliva droplets. T, temperature; RH, relative humidity; dp, droplet diameter; dp,o, initial droplet diameter; t, time.


4.1 Saliva Residue

Under a temperature between 15°C and 30°C and relative humidity between 20% and 70%, 32 saliva droplet samples with saliva initial droplet size ranging from 338 µm to 974 µm were tested for determining the ratio of saliva residue size to saliva initial droplet size, i.e., ψ = dp,residue/dp,o.

Nicas et al. (2005) predicted the saliva residue size to be about half that of the initial saliva droplet (ψ ≈ 0.5), using the saliva composition data by Effros et al. (2002). Liu et al. (2017) suggested that when the initial salt mass fraction was 0.9% and the initial solid volume ratio was 1.8%, the ratio of the saliva residue to the initial saliva droplet size was 0.32 (ψ = 0.32). However, it is worthwhile to mention that according to Radha et al. (2016)’s studies and Schipper et al. (2007)’s review, the saliva composition of solid constituents which includes inorganic solids (e.g., NaCl) and organic solids (e.g., mucin, serum) is only 0.5%. Therefore, both model (Liu et al., 2017; Nicas et al., 2005) needs to be further improved for predicting the ratio of saliva residue size to saliva initial droplet size.

Recently, the average ratios of saliva residue size to saliva initial droplet size measured by Lieber et al. (2021) and Stiti et al. (2022) were 0.198, and 0.235, respectively. The results of saliva residue diameter (dp,residue) are shown as a function of saliva initial droplet diameter (dp,o) in Fig. 5. The average value of present measurement data corresponded to a ratio (ψ) of 0.216, and the coefficient of determination (R2) was 0.988. When comparing the average ratios (ψ) with other studies, our measurement values are considered to be reliable. In addition, although the effect of temperature, relative humidity, and initial droplet size was examined, any relation between ratio (ψ) and other parameters (RH, T, and dp,o) was not found as presented in Fig. S3 in Supplementary Information. Therefore, regardless of other parameters (i.e., RH, T, and dp,o), the ratio of saliva residue size to saliva initial droplet size is fixed to be 0.216 (ψ = 0.216).

Fig. 5. Comparison of saliva initial droplet diameter (dp,saliva) and saliva residue diameter (dp,o). The blue open square symbol and solid line indicate the experimental results and linear fitting curve (ψ = 0.216, R2 = 0.988), respectively. The black dash-dotted line and dashed line indicate the averaged ratio (ψ) obtained by Lieber et al. (2021) and Stiti et al. (2022), respectively.Fig. 5. Comparison of saliva initial droplet diameter (dp,saliva) and saliva residue diameter (dp,o). The blue open square symbol and solid line indicate the experimental results and linear fitting curve (ψ = 0.216, R2 = 0.988), respectively. The black dash-dotted line and dashed line indicate the averaged ratio (ψ) obtained by Lieber et al. (2021) and Stiti et al. (2022), respectively.

 
4.2 Saliva Evaporation Coefficient

Fig. 6 shows the water and saliva evaporation coefficients as a function of relative humidity and temperature, respectively. As the relative humidity increases, the evaporation coefficient decreases almost linearly. This is because the smaller the amount of water vapor in the atmosphere, the more water vapor the atmosphere can receive from the droplet surface, so the lower the relative humidity, the faster the droplet evaporates, resulting in a higher evaporation coefficient value. On the other hand, as the temperature increases, the number of evaporating water molecules increases because the water molecules move more vigorously at high temperatures (i.e., high diffusion coefficient). Therefore, the evaporation coefficient increases as the temperature increases. In addition, Eq. (8) for obtaining the water evaporation coefficient shows that the water evaporation coefficient can be analytically obtained without much error from the experimental results, as shown in Fig. 6.

Fig. 6. Effect of (a) relative humidity at T = 20°C and (b) temperature at RH = 40% on the evaporation coefficient (βv). The red solid line which is denoted as water theory was calculated using Eq. (8).Fig. 6. Effect of (a) relative humidity at T = 20°C and (b) temperature at RH = 40% on the evaporation coefficient (βv). The red solid line which is denoted as water theory was calculated using Eq. (8).
 

As illustrated in Fig. 6, it can be seen that under the same conditions, the evaporation coefficient of saliva droplets was mostly measured lower than that of water droplets (i.e., the saliva droplet slightly slowly evaporates compared to the water droplet). This is attributed to the fact that the ions and proteins in saliva create an attraction, resulting in preventing water molecules from going into the atmosphere. The water evaporation coefficient (βv,water) and the saliva evaporation coefficient (βv,saliva) were compared to obtain the saliva correction factor (ϕsaliva), as illustrated in Fig. 7. The effect of temperature and relative humidity on the evaporation coefficient is presented in Fig. S4 in Supplementary Materials. Any relation between (saliva evaporation coefficient) and (RH, and T) was not found under the tested conditions (20% ≤ RH ≤ 70%, 15°C ≤ T ≤ 30°C).

Fig. 7. Comparison of water and saliva droplet evaporation coefficient. The red dashed line represents the linear fitting curve (ϕsaliva = 0.857, R2 = 0.984). ϕsaliva, saliva correction factor.Fig. 7. Comparison of water and saliva droplet evaporation coefficient. The red dashed line represents the linear fitting curve (ϕsaliva = 0.857, R2 = 0.984). ϕsaliva, saliva correction factor.
 

Because the saliva evaporation coefficient is assumed to be independent of RH and T for our tested conditions, the saliva correction factor can be assumed to be a ratio of water evaporation coefficient to saliva evaporation coefficient, i.e., ϕsaliva = βv,saliva/βv,water. The saliva correction factor was determined to be 0.857 with a high coefficient of determination (ϕsaliva = 0.857, R2 = 0.984). By using this factor, the saliva droplet evaporation in the air can be calculated.

 
4.3 Characteristics of Saliva Airborne Lifetime

For predicting airborne lifetime (tlife), the quiescent ambient conditions are assumed. The settling velocity (VT) in quiescent ambient conditions can be calculated depending on Reynolds number (Rep = ρVdp/µ). If particle or droplet size approaches the mean free path of the gas, slip phenomena occurs at the surface the particles. To account for the slip effect in Stoke’s law, the slip correction factor was (Cc) introduced. For particle motion with Rep > 1.0, settling velocity should be determined by equating the general form of Newton’s resistance to the force of gravity (Hinds, 1982). Therefore, the settling velocity can be formulated as follows:

 

Here, the correlation of the drag coefficient (CD) (Hinds, 1982) and slip correction factor (Cc) (Kwak et al., 2020, 2021b, 2023; Kim et al., 2005) can be calculated as

 

and

 

where Kn is Knudsen number (Kn = 2λ/dp). It is worthwhile to mention that the trial-and-error solution should be employed for obtaining settling velocity for large Reynolds numbers (Rep > 1.0). The detailed airborne lifetime calculation process and method are described in Fig. 8. The airborne lifetime without evaporation, and after perfectly evaporated can be calculated as follows:

 

and

 

where VT,o and tlife,o are the settling velocity and airborne lifetime without evaporation, i.e., dp(t) = dp,o, respectively. VT,r and tlife,r are the settling velocity and airborne lifetime after completely evaporated, i.e., dp(t) = ϕsalivadp,o, respectively. y* is the height above the ground. The time taken for complete evaporation (teva) (i.e., the time it takes to evaporate from the initial droplet size to the residue size) can be determined as follows:

 

To investigate the effect of saliva initial droplet size on airborne lifetime, the classical Wells’ evaporation-falling curve (Wells, 1934) was reproduced under the temperature of 25°C, relative humidity of 0%, and height of 2 m, as shown in Fig. 9. The initial saliva droplet diameters ranging from 1 µm to 500 µm were divided with a logarithmic constant interval of 0.02.

Fig. 8. Flow diagram for obtaining the saliva evaporation-falling curve.Fig. 8. Flow diagram for obtaining the saliva evaporation-falling curve.

Fig. 9. Results of calculated saliva airborne lifetime, i.e., evaporation-falling curve (red open symbol), saliva airborne lifetime as initial droplet size (purple dashed line), saliva airborne lifetime as residue size (green dash-dotted line), and time taken to completely evaporate (blue solid line) under the temperature of 25°C and relative humidity of 0%. Regime I, II, and III can be categorized by dominant effects: Regime I, the effect of residue size; Regime II, transition; Regime III, the effect of initial droplet size.Fig. 9. Results of calculated saliva airborne lifetime, i.e., evaporation-falling curve (red open symbol), saliva airborne lifetime as initial droplet size (purple dashed line), saliva airborne lifetime as residue size (green dash-dotted line), and time taken to completely evaporate (blue solid line) under the temperature of 25°C and relative humidity of 0%. Regime I, II, and III can be categorized by dominant effects: Regime I, the effect of residue size; Regime II, transition; Regime III, the effect of initial droplet size.

Fig. 9 can be divided into three regimes (e.g., Regime I, II, and III) according to dominant effects. In Regime I (effect of residue size), the initial droplet size is so small that the time taken for complete evaporation (teva) is much less than 10 seconds. Therefore, when the initial droplet diameter is smaller than 50−70 µm, the saliva airborne lifetime is almost equal to the time it takes to reach the ground by free-falling as the residue size (tlife,r, green dash-dotted line in Fig. 9), resulting in the airborne lifetime of more than 300 seconds. In Regime III (effect of initial droplet size), complete evaporation takes more than 10 seconds due to the large initial droplet size. The settling velocity is too large, so it is hardly affected by evaporation. Therefore, the airborne lifetime is approximately equal to the time it takes to reach the ground by free-falling as the initial droplet size (tlife,o, purple dashed line in Fig. 9). In Regime II (transition), the time it takes for complete evaporation and the time it takes to reach the ground by free-falling as the initial droplet size are not significantly different. Thus, a combined results of the two effects (residue size and initial droplet size) is shown. According to a van Doremalen (2000)’s study, the half-lives of SARS-CoV-2 in aerosols were approximately 4000 seconds, corresponding to the airborne lifetime of initial saliva droplets diameter of ~20 µm. This indicates that respiratory droplets with an initial diameter smaller than 20 µm will mostly contribute to airborne transmission due to their long airborne lifetime as well as long virus survival time. For respiratory droplets with an initial diameter larger than 100 µm, they will deposit on the ground or object surface due to a short airborne lifetime of fewer than 10 seconds, contributing mostly to fomite transmission. It is worthwhile to mention that once saliva droplets deposit on surfaces, the half-lives of SARS-CoV-2 dramatically increased (e.g., 1.09 hours in aerosols; 3.46 hours in cardboard; 5.63 hours in steel; 6.81 hours in plastic) (van Doremalen et al., 2020).

Fig. 10 shows the effect of temperature and relative humidity on the saliva airborne lifetime. As temperature increases or relative humidity decreases, the time taken to completely evaporate (blue line in Fig. 10) decreases. The time it takes to reach the ground by free-falling to either the residue size (tlife,r, green line in Fig. 10) or the initial droplet size (tlife,o, purple line in Fig. 10) is not significantly affected by ambient conditions (e.g., temperature and humidity). However, in Regime II (transition), the ambient conditions have a large influence on the saliva airborne lifetime. For example, for an initial saliva droplet size of 79.4 µm, when the temperature was increased from 10°C to 30°C, the saliva airborne lifetime at RH = 30% increased from 99 to 172 seconds. On the other hand, when the relative humidity was changed from 10% to 70%, the saliva airborne lifetime at T = 20°C was greatly reduced from 166 to 18 seconds. The reason that the saliva airborne lifetime rapidly changes in Regime II according to temperature and relative humidity is that it is more affected by the evaporation coefficient than in the other two Regimes.

Fig. 10. Effect of (a) temperature at RH = 30% and (b) relative humidity at T = 20°C on saliva airborne lifetime (red), and time taken to completely evaporate (blue). T, temperature; RH, relative humidity.Fig. 10. Effect of (a) temperature at RH = 30% and (b) relative humidity at T = 20°C on saliva airborne lifetime (red), and time taken to completely evaporate (blue). T, temperature; RH, relative humidity.

To predict the saliva airborne lifetime without complicate calculation which is presented in Fig. 8, the simple correlation equation was introduced. As shown in Fig. 8, the saliva airborne lifetime in Regimes I and III is almost the same as in tlife,r, and tlife,o, respectively, and in Regime II, the graph has an inverted S-shaped curve. Therefore, a correlation equation was created by slightly modifying the sigmoid function. The saliva airborne lifetime (i.e., evaporation-falling curve) can be predicted by using a simple correlation equation as follows:

 

where dp* is droplet diameter when the completely evaporated time (teva) and an airborne lifetime without evaporation (tlife,o) is the same, i.e., dp at teva = tlife,o. The correlation constants in Eq. (16) were determined as A = 30, and B = 0.05 for the height above the ground (y*) of 2 m. As shown in Fig. 11, there is almost no difference between correlation calculation results and numerical calculation results. More detailed comparison results between correlation and numerical calculation are presented in Fig. S5 in Supplementary Materials. The error is less than 0.7% for Regimes I and III and less than 2.2% for Regime II. It is worthwhile to mention that it takes about 1−3 hours to calculate by numerical method with a personal computer (2.8-GHz-CPU, Matlab 2018a), but it takes about a few seconds when calculating through the correlation equation.

Fig. 11. Correlation results of the saliva airborne lifetime (i.e., evaporation-falling curve) (red lines) under (a) various temperatures at RH = 30%, and (b) various relative humidity at T = 20°C. Open symbols represent the numerical calculation results. Grey lines represent saliva airborne lifetime as initial droplet size or residue size. T, temperature; RH, relative humidity.Fig. 11. Correlation results of the saliva airborne lifetime (i.e., evaporation-falling curve) (red lines) under (a) various temperatures at RH = 30%, and (b) various relative humidity at T = 20°C. Open symbols represent the numerical calculation results. Grey lines represent saliva airborne lifetime as initial droplet size or residue size. T, temperature; RH, relative humidity.

 
5 CONCLUSIONS


In the present study, we compared the characteristics of water and saliva evaporation phenomena by using a levitator. To carefully control and adjust the temperature and relative humidity conditions, all experiments were performed in the environmental chamber. The image processing algorithm was developed to obtain the size evolution during droplet evaporation. By introducing a dimensionless evaporation equation, the evaporation coefficients of water and saliva are defined. The experimentally obtained water evaporation results showed that the initial droplet sizes under the same environmental condition were less affected by the evaporation coefficients, i.e., the experimental results well follow the d2-law.

Under the evaporation experiment condition of a temperature (T) between 15°C and 30°C, relative humidity (RH) between 20% and 70%, and saliva initial droplet size (dp,o) ranging from 338 µm to 974 µm, we found the following things:

  • After the volatile liquids inside the saliva droplets are evaporated, the saliva residues are created. Regardless of other parameters (i.e., RH, T, and dp,o), the ratio of saliva residue size to saliva initial droplet size is fixed to be 0.216 (ψ = 0.216, R2 = 0.988).

  • Saliva droplet slowly evaporates compared to water droplet due to attraction by ions/proteins in saliva droplet. By comparing water and saliva droplet evaporation coefficients, the saliva correction factor was determined to be 0.857. (ϕsaliva = 0.857, R2 = 0.984).

  • Due to the slower evaporation rate and residue of saliva droplets, the lifetime of saliva droplets is expected to be longer than water droplets.

The detailed evaporation-falling curve calculation method is presented. For obtaining the saliva evaporation-falling curve, the authors recommend using the numerical calculation method which is presented in Fig. 8 if the readers want to calculate more accurate results, otherwise, we recommend for readers who want to see the general trend of the evaporation-falling curve according to various environmental conditions to use Eq. (16).

 
ACKNOWLEDGMENTS


The authors would like to express our deepest appreciation to Dr. Chao-Hsin Lin of Boeing Company for suggesting a great research topic. The authors would also like to extend our deepest gratitude to Dr. Bruce Forsyth of Boston Scientific for donating the acoustic levitator, and Dr. Seong Chan Kim and Douglas Brinkman for helping set up the levitation experiment in the controlled environmental chamber. The authors express their appreciation to the Center for Filtration Research: 3M Corporation, Applied Materials, Inc., BASF Corporation, Boeing Company, Corning Co., China Yancheng Environmental Protection Science and Technology City, Cummins Filtration Inc., Donaldson Company, Inc., Entegris, Inc., Ford Motor Company, Guangxi Wat Yuan Filtration System Co., Ltd, LG Electronics Inc., MSP Corporation, Parker Hannifin, Samsung Electronics Co., Ltd., Xinxiang Shengda Filtration Technology Co., Ltd., Shigematsu Works Co., Ltd., TSI Inc., W. L. Gore & Associates, Inc., and the affiliate member National Institute for Occupational Safety and Health (NIOSH). URL: http://www.me.umn.edu/cfr/.


REFERENCES


  1. Al-Shammiri, M. (2002). Evaporation rate as a function of water salinity. Desalination 150, 189–203. https://doi.org/10.1016/S0011-9164(02)00943-8

  2. Alsved, M., Matamis, A., Bohlin, R., Richter, M., Bengtsson, P.E., Fraenkel, C.J., Medstrand, P., Löndahl, J. (2020). Exhaled respiratory particles during singing and talking. Aerosol Sci. Technol. 54, 1245–1248. https://doi.org/10.1080/02786826.2020.1812502

  3. Asadi, S., Wexler, A.S., Cappa, C.D., Barreda, S., Bouvier, N.M., Ristenpart, W.D. (2020). Effect of voicing and articulation manner on aerosol particle emission during human speech. PLoS One 15, e0227699. https://doi.org/10.1371/journal.pone.0227699

  4. Asadi, S., Wexler, A.S., Cappa, C.D., Barreda, S., Bouvier, N.M., Ristenpart, W.D. (2019). Aerosol emission and superemission during human speech increase with voice loudness. Sci. Rep. 9, 2348. https://doi.org/10.1038/s41598-019-38808-z

  5. Cao, Q., Kuehn, T.H., Kim, S.C., Ou, Q., Pei, C., Pui, D.Y.H. (2021). An experimentally validated analytical model for aerosol number concentration reduction in classrooms. Aerosol Air Qual. Res. 21. 210038. https://doi.org/10.4209/aaqr.210038

  6. Centers for Disease Control and Prevention (CDC) (2021). Scientific Brief: SARS-CoV-2 Transmission. (Accessed 20 December 2021).

  7. Chao, C.Y.H., Wan, M.P., Morawska, L., Johnson, G.R., Ristovski, Z.D., Hargreaves, M., Mengersen, K., Corbett, S., Li, Y., Xie, X., Katoshevski, D. (2009). Characterization of expiration air jets and droplet size distributions immediately at the mouth opening. J. Aerosol Sci. 40, 122–133. https://doi.org/10.1016/j.jaerosci.2008.10.003

  8. Chen, W., Kwak, D.B., Anderson, J., Kanna, K., Pei, C., Cao, Q., Ou, Q., Kim, S.C., Kuehn, T.H., Pui, D.Y.H. (2021). Study on droplet dispersion influenced by ventilation and source configuration in classroom settings using low-cost sensor network. Aerosol Air Qual. Res. 21, 210232. https://doi.org/10.4209/aaqr.210232

  9. Cheng, P., Luo, K., Xiao, S., Yang, H., Hang, J., Ou, C., Cowling, B.J., Yen, H.-L., Hui, D.S.C., Hu, S., Li, Y. (2022). Predominant airborne transmission and insignificant fomite transmission of SARS-CoV-2 in a two-bus COVID-19 outbreak originating from the same pre-symptomatic index case. J. Hazard. Mater. 425. 128051. https://doi.org/10.1016/j.jhazmat.2021.128051

  10. Dao, H.T., Kim, K.S. (2022). Behavior of cough droplets emitted from Covid-19 patient in hospital isolation room with different ventilation configurations. Build. Environ. 209. 108649. https://doi.org/10.1016/j.buildenv.2021.108649

  11. Duguid, J.P. (1946). The size and the duration of air-carriage of respiratory droplets and droplet-nuclei. Epidemiol. Infect. 44, 471–479. https://doi.org/10.1017/S0022172400019288

  12. Echternach, M., Gantner, S., Peters, G., Westphalen, C., Benthaus, T., Jakubaß, B., Kuranova, L., Döllinger, M., Kniesburges, S. (2020). Impulse dispersion of aerosols during singing and speaking: A potential COVID-19 transmission pathway. Am. J. Respir. Crit. Care Med. 202, 1584–1587. https://doi.org/10.1164/rccm.202009-3438LE

  13. Edwards, D.A., Ausiello, D., Salzman, J., Devlin, T., Langer, R., Beddingfield, B.J., Fears, A.C., Doyle-Meyers, L.A., Redmann, R.K., Killeen, S.Z., Maness, N.J., Roy, C.J. (2021). Exhaled aerosol increases with COVID-19 infection, age, and obesity. PNAS 118, e2021830118. https://doi.org/​10.1073/pnas.2021830118

  14. Effros, R.M., Hoagland, K.W., Bosbous, M., Castillo, D., Foss, B., Dunning, M., Gare, M., Lin, W., Sun, F. (2002). Dilution of respiratory solutes in exhaled condensates. Am. J. Respir. Crit. Care Med. 165, 663–669. https://doi.org/10.1164/ajrccm.165.5.2101018

  15. Fakir, M.H., Kim, J.K. (2022). Prediction of individual thermal sensation from exhaled breath temperature using a smart face mask. Build. Environ. 207, 108507. https://doi.org/10.1016/j.​buildenv.2021.108507

  16. Faulkner, C.A., Castellini, J.E., Zuo, W., Lorenzetti, D.M., Sohn, M.D. (2022). Investigation of HVAC operation strategies for office buildings during COVID-19 pandemic. Build. Environ. 207, 108519. https://doi.org/10.1016/j.buildenv.2021.108519

  17. Frohn, A., Roth, N. (2000). Dynamics of droplets. Springer Science & Business Media.

  18. Han, Z.Y., Weng, W.G., Huang, Q.Y. (2013). Characterizations of particle size distribution of the droplets exhaled by sneeze. J. R. Soc. Interface 10, 20130560. https://doi.org/10.1098/rsif.​2013.0560

  19. Hare, R. (1964). The transmission of respiratory infections. Proc. R. Soc. Med. 57, 221–230. https://doi.org/10.1177/003591576405700329

  20. Hinds, W.C. (1982). Aerosol Technology: Properties, Behaviour, and Measurement of Airborne Particles. John Wiley & Sons.

  21. Hisatake, K., Tanaka, S., Aizawa, Y. (1993). Evaporation rate of water in a vessel. J. Appl. Phys. 73, 7395–7401. https://doi.org/10.1177/003591576405700329

  22. Huang, J. (2018). A simple accurate formula for calculating saturation vapor pressure of water and Ice. J. Appl. Meteorol. Climatol. 57, 1265–1272. https://doi.org/10.1175/JAMC-D-17-0334.1

  23. Kim, J.H., Mulholland, G.W., Kukuck, S.R., Pui, D.Y.H. (2005). Slip correction measurements of certified PSL nanoparticles using a nanometer differential mobility analyzer (Nano-DMA) for knudsen number from 0.5 to 83. J. Res. Natl. Inst. Stand. Technol. 110, 31–54. https://doi.org/​10.6028/jres.110.005

  24. Kim, S.C., Kwak, D.B., Kuehn, T., Pui, D.Y.H. (2021). Characterization of handheld disinfectant sprayers for effective surface decontamination to mitigate severe acute respiratory coronavirus virus 2 (SARS-CoV-2) transmission. Infect. Control Hosp. Epidemiol. 42, 901–903. https://doi.org/10.1017/ice.2020.1423

  25. Kukkonen, J., Vesala, T., Kulmala, M. (1989). The interdependence of evaporation and settling for airborne freely falling droplets. J. Aerosol Sci. 20, 749–763. https://doi.org/10.1016/0021-8502(89)90087-6

  26. Kwak, D.B., Kim, S.C., Lee, H., Pui, D.Y.H. (2020). Experimental study of nanoparticle transport and penetration efficiency on a sharp-bent tube (elbow connection). Int. J. Heat Mass Transf. 155, 119816. https://doi.org/10.1016/j.ijheatmasstransfer.2020.119816

  27. Kwak, D.B., Kim, S.C., Kuehn, T.H., Pui, D.Y.H. (2021a). Quantitative analysis of droplet deposition produced by an electrostatic sprayer on a classroom table by using fluorescent tracer. Build. Environ. 205, 108254. https://doi.org/10.1016/j.buildenv.2021.108254

  28. Kwak, D.B., Kim, S.C., Lee, H., Pui, D.Y.H. (2021b). Numerical investigation of nanoparticle deposition location and pattern on a sharp-bent tube wall. Int. J. Heat Mass Transf. 164, 120534. https://doi.org/10.1016/j.ijheatmasstransfer.2020.120534

  29. Kwak, D.B., Wang, J.S., Lee, H., Pui, D.Y.H. (2023). Numerical study of nanoparticle penetration characteristics in forked tubes using tracking particle identification. Powder Technol. 415, 118172. https://doi.org/10.1016/j.powtec.2022.118172

  30. Lai, C.C., Shih, T.P., Ko, W.C., Tang, H.J., Hsueh, P.R. (2020). Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and coronavirus disease-2019 (COVID-19): The epidemic and the challenges. Int. J. Antimicrob. Agents 55, 105924. https://doi.org/10.1016/j.ijantimicag.2020.​105924

  31. Lieber, C., Melekidis, S., Koch, R., Bauer, H.J. (2021). Insights into the evaporation characteristics of saliva droplets and aerosols: Levitation experiments and numerical modeling. J. Aerosol Sci. 154, 105760. https://doi.org/10.1016/j.jaerosci.2021.105760

  32. Liu, L., Wei, J., Li, Y., Ooi, A. (2017). Evaporation and dispersion of respiratory droplets from coughing. Indoor Air 27, 179–190. https://doi.org/10.1111/ina.12297

  33. Liu, L., Iketani, S., Guo, Y., Chan, J.F.W., Wang, M., Liu, L., Luo, Y., Chu, H., Huang, Y., Nair, M.S., Yu, J., Chik, K.K.H., Yuen, T.T.T., Yoon, C., To, K.K.W., Chen, H., Yin, M.T., Sobieszczyk, M.E., Huang, Y., Wang, H.H., et al. (2022). Striking antibody evasion manifested by the Omicron variant of SARS-CoV-2. Nature 602, 676–681. https://doi.org/10.1038/s41586-021-04388-0

  34. Morawska, L., Johnson, G.R., Ristovski, Z.D., Hargreaves, M., Mengersen, K., Corbett, S., Chao, C.Y.H., Li, Y., Katoshevski, D. (2009). Size distribution and sites of origin of droplets expelled from the human respiratory tract during expiratory activities. J. Aerosol Sci. 40, 256–269. https://doi.org/10.1016/j.jaerosci.2008.11.002

  35. Nicas, M., Nazaroff, W.W., Hubbard, A. (2005). Toward understanding the risk of secondary airborne infection: Emission of respirable pathogens. J. Occup. Environ. Hyg. 2, 143–154. https://doi.org/10.1080/15459620590918466

  36. Niimura, Y., Hasegawa, K. (2019). Evaporation of droplet in mid-air: Pure and binary droplets in single-axis acoustic levitator. PLoS One 14, e0212074. https://doi.org/10.1371/journal.pone.​0212074

  37. Ou, Q., Pei, C., Chan Kim, S., Abell, E., Pui, D.Y.H. (2020). Evaluation of decontamination methods for commercial and alternative respirator and mask materials – view from filtration aspect. J. Aerosol Sci. 150, 105609. https://doi.org/10.1016/j.jaerosci.2020.105609

  38. Radha, R., John Blesswin, A., Selva Mary, G. (2016). A simple innovative approach DNA-based saliva security system for user authentication. Indian J. Sci. Technol. 9, 1–6. https://doi.org/​10.17485/ijst/2016/v9i37/102064

  39. Sarkodie, S.A. and Owusu, P.A. (2021). Global assessment of environment, health and economic impact of the novel coronavirus (COVID-19). Environ. Dev. Sustain. 23, 5005–5015. https://doi.org/10.1007/s10668-020-00801-2

  40. Schiffter, H., Lee, G. (2007). Single-droplet evaporation kinetics and particle formation in an acoustic levitator. Part 1: Evaporation of water microdroplets assessed using boundary-layer and acoustic levitation theories. J. Pharm. Sci. 96, 2274–2283. https://doi.org/10.1002/jps.20860

  41. Schipper, R.G., Silletti, E., Vingerhoeds, M.H. (2007). Saliva as research material: Biochemical, physicochemical and practical aspects. Arch. Oral Biol. 52, 1114–1135. https://doi.org/10.1016/​j.archoralbio.2007.06.009

  42. Schneider, C.A., Rasband, W.S., Eliceiri, K.W. (2012). NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 9, 671–675. https://doi.org/10.1038/nmeth.2089

  43. Stadnytskyi, V., Bax, C.E., Bax, A., Anfinrud, P. (2020). The airborne lifetime of small speech droplets and their potential importance in SARS-CoV-2 transmission. PNAS 117, 11875–11877. https://doi.org/10.1073/pnas.2006874117

  44. Stiti, M., Castanet, G., Corber, A., Alden, M., Berrocal, E. (2022). Transition from saliva droplets to solid aerosols in the context of COVID-19 spreading. Environ. Res. 204, 112072. https://doi.org/10.1016/j.envres.2021.112072

  45. Tuckermann, R., Bauerecker, S., Neidhart, B. (2002). Evaporation rates of alkanes and alkanols from acoustically levitated drops. Anal. Bioanal. Chem. 372, 122–127. https://doi.org/​10.1007/s00216-001-1132-7

  46. van Doremalen, N., Bushmaker, T., Morris, D.H., Holbrook, M.G., Gamble, A., Williamson, B.N., Tamin, A., Harcourt, J.L., Thornburg, N.J., Gerber, S.I., Lloyd-Smith, J.O., de Wit, E., Munster, V.J. (2020). Aerosol and surface stability of SARS-CoV-2 as compared with SARS-CoV-1. N. Engl. J. Med. 382, 1564–1567. https://doi.org/10.1056/NEJMc2004973

  47. Wan, M.P., Sze To, G.N., Chao, C.Y.H., Fang, L., Melikov, A. (2009). Modeling the fate of expiratory aerosols and the associated infection risk in an aircraft cabin environment. Aerosol Sci. Technol. 43, 322–343. https://doi.org/10.1080/02786820802641461

  48. Wells, W.F. (1934). On air-borne infection. Study II. Droplets and droplet nuclei. Am. J. Epidemiol. 20, 611–618. https://doi.org/10.1093/oxfordjournals.aje.a118097

  49. World Health Organization (WHO) (2023). Coronavirus (COVID-19) Dashboard. (Accessed 9 January 2023).

  50. Xie, X., Li, Y., Sun, H., Liu, L. (2009). Exhaled droplets due to talking and coughing. J. R. Soc. Interface 6, S703–S714. https://doi.org/10.1098/rsif.2009.0388.focus

  51. Yang, S., Lee, G.W.M., Chen, C.M., Wu, C.C., Yu, K.P. (2007). The size and concentration of droplets generated by coughing in human subjects. J. Aerosol Med. 20, 484–494. https://doi.org/​10.1089/jam.2007.0610

  52. Zack, G.W., Rogers, W.E., Latt, S.A. (1977). Automatic measurement of sister chromatid exchange frequency. J. Histochem. Cytochem. 25, 741–753. https://doi.org/10.1177/25.7.70454

  53. Zhang, T. (2011). Study on surface tension and evaporation rate of human saliva, saline, and water droplets. MS thesis, West Virginia University. https://doi.org/10.33915/etd.2271


Share this article with your colleagues 

 

Subscribe to our Newsletter 

Aerosol and Air Quality Research has published over 2,000 peer-reviewed articles. Enter your email address to receive latest updates and research articles to your inbox every second week.

7.3
2022CiteScore
 
 
77st percentile
Powered by
Scopus
 
   SCImago Journal & Country Rank

2022 Impact Factor: 4.0
5-Year Impact Factor: 3.4

The Future Environment and Role of Multiple Air Pollutants

Aerosol and Air Quality Research partners with Publons

CLOCKSS system has permission to ingest, preserve, and serve this Archival Unit
CLOCKSS system has permission to ingest, preserve, and serve this Archival Unit

Aerosol and Air Quality Research (AAQR) is an independently-run non-profit journal that promotes submissions of high-quality research and strives to be one of the leading aerosol and air quality open-access journals in the world. We use cookies on this website to personalize content to improve your user experience and analyze our traffic. By using this site you agree to its use of cookies.