1 INTRODUCTION
The outbreak of COVID-19 disease at superspreading events intensifies the scientific discussion on the airborne infection transfer mechanism responsible for the global spread of this disease 1-4; most of all for combating strategies and then to the issue of rebooting economic activity. Limited information is available on the formation of virus-laden droplets/aerosols, dispersion, and evaporation of expiratory droplet cloud in the indoor environment 5, 6. Furthermore, there is a wide variation in the input parameters such as droplet size distribution, concentration, volume of expiratory fluid, release conditions, etc. Accordingly, a considerable amount of uncertainty is associated with modeling practices used for inhalation exposure assessment. Although there are few methods of exposure assessment available in the literature 7-9, it is necessary to address larger scientific questions which form the basis of these calculations. In this study, we wish to focus on the issue of residence time of cough droplet and droplet nuclei, ejected from the expiratory activity of infected subjects as localized puff to start with, in an indoor space. The residence times are not only important for exposure assessment, but are also important for assessing the performance of air cleaners, or for deciding the minimum level of Clean Air Delivery Rate (CADR) expected of them, as effective intervention technologies for minimizing the infection transfer risks in confined environments.
Generally, the mean residence time is derived from the reciprocal of the removal rate, which is a concept based on uniformly mixed chamber models10. This concept is best suited for fine aerosol particles of about 1 μm or less, predominantly occurring in the context of indoor pollution. On the other extreme is the models of residence time for falling-evaporating droplets, as described by Wells11 long ago and revisited with considerable mathematical precision by Xie et al.12, which are best suited to describe the fate of droplets of 20 μm or more. There is a wide gap in modeling the residence time for droplets in the intervening size range (say, 5–20 μm), which due to natural convective mixing in a room, will be undergoing dispersal in the air currents during their fall under gravity, and cannot be described by either of the above models. With increasing concern expressed over the yet to be understood role of aerosol route of infection spread in the COVID-19 scenario, the assessment of the impact of size range is of considerable urgency. The first motivating factor for this paper is to formulate a model that smoothly transits from the falling droplet to a uniformly mixed chamber model (Falling-to-Mixing-Plate-out) across the entire size range of polydisperse droplets released in a cough.
A second and more important motivation for this study stems from the following consideration. In a recent paper13, we pointed out that the size distribution metrics of corona virus-laden polydisperse droplets expirated from infected subjects, could be very different from the original droplet distributions, due to the well-known Poissonian fluctuations13-17 in virus incorporation propensities into droplets at the time of their formation. To focus on this viral load-dependent distinctiveness, a terminology “virusol” was proposed for the ensemble of virus-laden droplets. From both exposure evaluation and control technology point of view, the residence time of the virusolized fraction is a matter of concern, and not that of the entire droplet population. The Poissonian theory shifts focus on larger size fraction effectively rendering the residence times sensitive to the viral load in the ejecta of subjects. Examining this dependence is an important aspect of this study.
Size distribution of expiratory droplets is one of the basic input parameters required to determine the lifetime of the virus-laden droplets in the ambient environment. There are several experimental studies on the size distribution and concentration of expiratory droplets 7, 18-20. These studies show a large variation in the size distribution parameters of the droplets ejected from the infected subjects. Equally important is the data on viral load among patients and the variability across several findings have been summarized by Anand and Mayya13. The combined effect of size dispersity and viral load has a complex influence on the lifetime of the virus-laden droplets as we will see in this paper.
Computational fluid dynamics (CFD) is an increasingly adopted approach to address the transport of droplets along with their evaporative dynamic, from the point of emission to their settling on surfaces 21-24. In a strictly theoretical sense, this would be most accurate for modeling residence time distributions. However, a large number of parameters are required for carrying out simulation using fluid dynamics. Therefore, CFD approaches are highly input parameter intensive. For example, to handle turbulence in a room it would be necessary to specify all the processes responsible for creating turbulence. It should then include both mechanical and thermal sources, and this would make the CFD models quite useful in specific contexts, but lack generic predictive power. To have generality, CFD would require averaging over a large number of scenarios before arriving at the reasonable estimates of a representative droplet residence time. On the other hand, phenomenological model based on turbulent diffusion as the governing transport mechanism in addition to gravitational settling addresses the dispersive mechanisms directly through an assumed diffusion coefficient and hence will be more robust in providing generic estimates of the residence time distributions. Besides, treating droplet transport within the framework of an assumed bulk diffusion process obviates the need to examine the origin of its occurrence for the purposes of describing the migration dynamics of the droplets. Equally important, turbulent bulk diffusion models are amenable to analytical solutions for falling droplets which is useful for comparison purposes for providing insights into the role of individual processes.
In the phenomenological approach, the emitted droplet begins to disperse in bulk space as it settles under gravity simultaneously evaporating and undergoing wall removal processes. The removal rates are matched with the boundary conditions using Lai and Nazaroff model10. This phenomenological model integrates smoothly the wall removal as well as air exchange removal without having to concern about the points of air entry and exit for various realistic emission scenarios. We term it as a residence time distribution problem since the released droplets are not uniformly mixed in the indoor environment, and the inverse of the mean residence time is equivalent to the removal rate of uniformly mixed model. The dynamics of these evaporating and settling virusol puff is studied under different conditions by varying the important input parameters such as relative humidity, differential temperature between the droplet and ambient atmosphere, release height, bulk turbulent coefficient, friction velocity, initial droplet size distribution, multi-component nature of the droplets, and viral load in the biological fluid. The issue of the dependence of the mean residence time of virus-laden droplets (virusols) on the viral load forms a key component of our analysis.
2 MATERIALS AND METHODS
Let us consider an expiratory event in which droplets of initial concentration (N0) are released from an infected person at a height (z0) in a closed environment (Figure 1). The convective diffusive expansion of the puff containing airborne droplets is modeled using the mixing transition model, which combines homogeneous mixed room model with falling droplet model seamlessly. Many researchers 6, 25, 26 argue that the cut-off diameter for the aerosol definition in the present context should be a single discrete number (e.g., <5 μm or <20 μm). However, in the present study, a higher cut-off diameter (100 μm) is considered so that the transition to inhalable aerosol size from the evaporation of continuous size spectrum of ejected droplets can be modeled seamlessly.
2.1 Falling-to-mixing-plate-out model
2.2 Mean residence time (MRT)
The definition of residence time is more general compared to the lifetime definition used in single exponential decay expressions. In the latter case, the mean residence time refers to time required to decrease 1/eth of original concentration, equivalent to ~63% removal in the single exponential decay process. Equation (3) is a measure of multi-exponential decay due to various processes acting on the removal/survival of the expiratory droplets.
2.3 Vertical motion of evaporating and settling droplets
The droplets undergo evaporation in a homogeneous temperature Tamb with a constant ambient water vapor pressure at a given relative humidity, till they reach equilibrium conditions. The evaporation model neglects the effect of Stefan flow on heat and mass transfer between the droplet and the surrounding gas since this effect leads to less than 0.5% change in the model variables. The airborne droplets in the puff are assumed to be spherical particles consist of non-volatile components such as Na+, K+, Cl- ions, and lactate and glycoproteins (about 0.71% mole fraction) with water as a major constituent.
2.4 Virusol system
The virusol droplet number distribution is then obtained by multiplying the number-size distribution with at any instant of time t. It is observed that the variation of virus concentration in the biological fluid is very wide (102–1011 RNA copies/ml), and studies indicate that there is a possible linkage of viral load with the severity of disease29. Based on data, we approximately classify viral load into two categories, i.e., (i) mild-to-moderate cases – <106 RNA copies/ml and (ii) severe cases – >106 RNA copies/ml.
3 RESULTS AND DISCUSSION
The dispersion of a puff containing virus-laden droplets and aerosols exhaled during a cough expiratory event by an infected subject is analyzed using 1-D model described in the Materials and Methods section.
3.1 Input parameters and initial conditions
The puff is released at a height of 1.5 m (represents the elevation of infected subject’s standing position) in an indoor environment of volume 48 m3 (height of the room = 3 m) having an air-exchange rate of 1.0 h−1, representing a typical office environment. We consider the puff of cough droplets with a temperature of 35°C is released into an environment at the ambient temperature of 25°C at 50% RH, ideal condition for comfort working environment. A total number of 5000 droplets with lognormal number-size distribution (CMD = 14 μm and GSD = 2.6) is considered in this study. These expiratory droplets contain residues7, settle by gravity while evaporating in ambient air. Other important input parameters/constant and their range are given in Table 1. All the governing equations are numerically solved together in the Mathematica30.
3.2 Comparison of Falling-to-Mixing-Plate-out model results with falling droplet and uniformly mixed models
Early works of Wells11 and other studies12, 31, 32 highlight the dependence of lifetime of the droplets on ambient parameters such as temperature, relative humidity, air speed, and droplet properties such as its size, constituents, and temperature. The falling droplet evaporation model considers the droplet evolution processes in a localized space without considering the boundary wall deposition and ventilation effects. Also, some of these models do not consider detailed droplet composition, and hence, large variation in the final droplet diameters (~20%–50% of the initial values) which will lead to considerable error in the modeling of airborne droplet concentration. On the other hand, homogeneously mixed model10 assumes uniform mixing of expiratory droplets/aerosols in the indoor volume immediately. This model lacks the handling of spatial variation of the aerosol concentration that affects other accompanied processes. A study is carried out to compare these two models with the present formulation for non-evaporating rigid droplets released at a height of 0.3 m in an indoor room environment with a ventilation rate of 0.5 h−1. The results are presented in Figure 2. From Figure 2, it is quickly apparent that neither falling droplet model nor homogeneously mixed chamber model could adequately describe all the effects if we include coarse as well as submicron particles for assessing risks.
The present model (Equations (1–9-1–9)) results show smooth transition of residence time from the falling evaporating model to uniform mixing model as illustrated in Figure 2. These results confirm that the present approach seamlessly integrates all the important processes that governs the residence time of virus-laden droplets and aerosols in a typical indoor environment.
3.3 Temporal evolution of droplet diameter of individual droplets undergoing evaporation
The rate of evaporation of droplets and its composition determines the final size of airborne droplets/aerosols (virusols) which play a key role in governing their residence time in the indoor atmosphere. To study these effects, the virusols residence time is estimated at different RH, and compared with that of droplets without solute. The temporal evolution of single droplet diameter (14 μm is chosen since its CMD of the cough droplets considered in this study) at different RH and composition is shown in Figure 3A. At a lower RH (10%), the water vapor content in the atmosphere is less and hence the rate of evaporation is highest. So, the droplet quickly evaporates and reaches a steady size value (~6 μm) in a time period of 0.12 s and 0.3 s for an RH value of 10% and 50% respectively. At higher RH values (90%), the evaporation rate is least due to the saturated water vapor pressure in the indoor environment and hence, it takes more time (> 1 s) to reach the final size. Thus, the results (time scale) clearly show that it is important to integrate the evaporation process with other processes described in the governing equations.
Apart from the atmospheric conditions, the composition of the droplet also plays an important role in determining the rate of evaporation and final size of the droplet. The presence of salt (residue) content reduces the vapor pressure of water thereby leading to lower rate of evaporation and larger final size of the droplet. The comparison of variation of droplet size as a function of time for RH = 50% clearly indicates this fact as the final size reaches to 6.1 μm in the presence of solute while the pure water droplet evaporates completely in ~0.3 s (Figure 3A). The rate of evaporation of droplet thus depends on the volatile content of the droplet along with the ambient conditions.
Furthermore, the final size of all the droplets as a function of initial droplet size is presented in Figure 3B. The overall reduction in the final size is about ~55% for all particle sizes (i.e., droplet with a size of 100 µm reduces to ~45 µm) in the case of RH = 50% and for the given residue composition. This size reduction is more pronounced if solute content in the saliva is reduced as discussed above. It is to be noted that the initial droplet may be large enough to settle down under the force of gravity but simultaneous evaporation will reduce the size significantly (depending upon the atmospheric conditions), thereby reducing the gravitational effects and leading to higher residence time in the air.
3.4 Effect of AER on MRT
The effect of AER on virusol mean residence time is studied with two different relationships of AER vs D ((i) D as a function of AER given by23, and (ii) D is independent of AER), is shown in Figure 4. The bulk diffusion coefficient varies from ~10−3 m2/s to 1.1 × 10−2 m2/s for AER of 0 to 5 h−1 in a room volume of 48 m3. The virusol mean residence time decreases from 960 s to 410 s as the AER increases from 0 h−1 to 5 h−1, in the case of D ∝ AER. For the constant D case, the variation of mean residence time is from 1040 s (0 h−1) to 500 s (5 h−1). These results show that the virusol mean residence time decreases monotonically when the AER increases; however, the reduction is only 0.43 times when the AER is increased from 0 h−1 to 5 h−1. Hence, a large AER is required to reduce the virusol exposure significantly.
3.5 Effect of release height and ventilation rate on MRT
The virusol mean residence time of droplets as a function of initial droplet size is compared for different heights of release (z0 = 0.3 m or 1.5 m) during an expiratory event. The lower value of release height (0.3 m) is chosen to illustrate that the cough cloud does not settle immediately on the ground and buoyancy helps to lift the puff and droplets survives for reasonable time compared to higher release height in the indoors. The mean residence time of the virusols varies from few seconds to hours for droplets of different sizes from 0.4 to 100 µm (Figure 5). Higher size droplets settle faster due to gravitational effects and smaller size droplets have largest residence time in the air as expected. As the release height increases, the corresponding residence time of the droplets also increases for sizes >10 µm as the settling distance increases with height (~7 times increase in the residence time for 100 µm droplet between 0.3 m and 1.5 m release heights). In the case of smaller size particles (~2500 s for aerosols <10 µm for AER of 1.0 h−1), MRT difference is not significant with respect to z0 as the removal of these particles is predominantly due to diffusion and ventilation, and hence, residence times of smaller size virusols are independent of height of release. On the other hand, the study results show that the ventilation rate strongly affects the residence time of smaller size virusols compared to larger ones (Figure 5). At smaller sizes, the lower ventilation rate leads to highest virusol residence time (~6000 s for an AER of 0.001 h−1) compared to high ventilation rate (~700 s for an AER of 5 h−1) for particles <5 µm. Since AER is linearly related to bulk diffusion coefficient23, the combined effect of ventilation removal and convective diffusion determines the survival of airborne particles in the indoor atmosphere. Also, the effect of AER is negligible in the large particle size domain due to other competing processes such as evaporation and gravitational settling acting on them. The effect of friction velocity on residence time of the virusols is insignificant as the study shows that the variation of virusol residence time is within 4% if the friction velocity varies from 1 cm/s to 5 cm/s.
3.6 Effect of viral load on number-size distribution and MRT of virusol system
Viral concentration in the biological fluid and the process of incorporation of virus in the expiratory droplets/aerosols plays a crucial role in determining the total residence time of virus-laden airborne particles, which is addressed through the concept of virusols. Most of the studies consider virus load in the airborne droplet as linearly proportional to the droplet volume. However, the atomized droplets during the expiratory activity show that the presence of virus in the expiratory droplets/aerosols is a function of viral concentration and droplet size, and it is governed by the Poisson distribution as given by Equation (9).
The initial and final (at t = 100 s) number-size distribution of virusols are compared with initial complete droplet size spectrum in Figure 6. The complete droplet spectrum which has a CMD of 14 μm contains the droplets/aerosols starting from ~0.3 μm onwards with viral load proportional to their volume without considering the Poisson distribution of virions in the droplets. On the other hand, it is very clear from Figure 6 that the introduction of Poisson probability limits the lower particle size of virusols, and the probability of the droplet being virus-laden (virusol) increases with its size. For example, the initial virusol number-size distribution (at t = 0) shows that the lower cut-off virusol diameter is ~5 μm and ~1 μm for viral load of 108 and 1011 RNA copies/ml (corresponding to severe cases) respectively. This shows that there is very low probability that smaller droplets will contain any virus in the case of low viral load (mild-to-moderate cases) in biological fluid (e.g., saliva). It is interesting to note that there is a larger difference in the concentrations at the lower side of the size spectrum between complete droplets and virusols as compared to that at larger sizes. This is due to the fact that the value of P (Equation (9)) tends to unity for large droplets (size >40 μm). Comparison of number-size distribution at t = 100 s with other initial size spectra shows that most of the particles greater than ~70 μm are removed from the system mainly due to gravitational settling; also, the evaporation process plays an important role in shifting the droplet spectrum in the range of (22–70) μm significantly (Figure 6).
We further study the effect of virus concentration on the residence time and exposure for a wide range of concentration from 103 to 1013 RNA copies/ml, and the results are presented in Figure 7 for various RH. The virusol residence time varies from ~100 s (RH = 90% for viral load <106 RNA copies/ml) to ~1200 s (RH = 10% and viral load >1011 RNA copies/ml). The results further show that MRT of virusols in the severe cases (viral load >1011 RNA copies/ml) is eight times to that of mild-to-moderate cases (<106 RNA copies/ml) for RH = 50%, and this ratio increases with RH. As RH increases from 10% to 90%, the evaporation rate reduces which leads to larger final size of droplet; hence faster settling and smaller residence time. The study results show that the MRT of virusols increase with viral load and saturates to a constant value beyond 1011 RNA copies/mL for the given conditions, which clearly show larger dependence of virusol residence time on virus load.
The exposure to these virusols can be obtained from the mean residence time using the relationship, , which is defined as the number of virions breathed in during MRT of virusols (), where BR is the breathing rate (~10 lpm), is the room volume, and Cv is the total virions released during the expiratory event. For example, if we consider a viral load of 106 RNA copies/ml (RH = 50%), then the mean residence time is 150 s (Figure 7) and corresponding exposure is 0.23 virions considering room volume of 48 m3. This result shows that exposure to a single virion is possible only if the viral load in the biological fluid is more than 106 RNA copies/ml. However, if the viral load is 1010 RNA copies/ml (severe case), then the mean residence time is 1050 s (Figure 7) and the corresponding exposure is 1.6 × 104 virions. These study results show that higher residence time and survival fraction would lead to higher exposure of person contracting the virus, beyond a certain cut-off viral load in the biological fluid. The residence time of the droplets thus plays a significant role in estimating the viral exposure to an individual, and it depends on many parameters such as particle size, viral load, process of incorporation of viral copies in the droplet, height of release, turbulent diffusion coefficient, AER, and RH.
4 CONCLUSIONS
The present conjoined model, “Falling-to-Mixing-Plate-out”, smoothly transits from evaporating and falling droplet model11 to mixed reservoir/all surface-plate out model10 by considering all the important processes that govern the virusol residence time and exposure in an indoor environment. The residence time estimates from the conjoined model lie between that obtained by pure gravitational settling and classical surface-plate out model. The present 1-D model is good enough to quickly estimate virusol residence time and exposure to the individual as their results are closer to 3-D system, particularly for droplets >0.1 µm. The present work thus provides a comprehensive theoretical study from the droplet emission to inhalation exposure using residence time estimates.
Study results indicate that virusol system residence times could be far less than that for the complete droplet cloud system, and smaller droplets (<20 μm) are more likely to be blanks in mild-to-moderate cases. Virusol system will settle down faster than the complete droplet system and hence less likely to be airborne when the biological viral load is <108 RNA copies/ml. However, greater viral load in the biological fluid leads to longer residence time of virusol system. This has important implications from an air cleaning perspective aimed at COVID-19 risk mitigation in enclosed spaces. It may not be necessary to capture ultrafine particles as these would be generally harmless; one may relax filtration efficiency criteria limiting them to essentially larger particles, thereby opting for coarser filters and reducing pressure drops and increasing air circulation rates. This has far reaching implications in adopting intervention technologies in indoor spaces.
Also, there is a significant shift in the geometric mean diameter of virusols as compared to that of the complete size spectrum and it has practical implication to risk assessment. On the whole, we have effectively demonstrated that virusol size distribution and residence times are distinctly different from the ensemble of total droplets and will be a critically useful concept both for intervention technologies and lung deposition modeling. These implications are under investigation.
ACKNOWLEDGEMENTS
Authors (AS & JK) would like to acknowledge Shri Kapil Deo Singh and Dr M S Kulkarni of Bhabha Atomic Research Centre for their support and encouragement. Authors would like to acknowledge Ms Riya Dey for the preparation of the artwork (Figure 1).
CONFLICT OF INTEREST
The authors declare no competing interests.
AUTHOR CONTRIBUTIONS
Y.S.M. and A.S. initiated the study and suggested the research method; S.A. and S.B. collected and analyzed data for the research method; J.K., Y.S.M., and A.S. performed numerical calculations and prepared figures; A.S. and Y.S.M. wrote the manuscript; S.B. and J.K. helped in preparation of the manuscript. All the authors reviewed the manuscript.
"load" - Google News
June 12, 2021 at 02:00PM
https://ift.tt/3ztYWKz
Modeling the viral load dependence of residence times of virus‐laden droplets from COVID‐19‐infected subjects in indoor environments - Wiley
"load" - Google News
https://ift.tt/2SURvcJ
https://ift.tt/3bWWEYd
Bagikan Berita Ini
0 Response to "Modeling the viral load dependence of residence times of virus‐laden droplets from COVID‐19‐infected subjects in indoor environments - Wiley"
Post a Comment