In the present work, the risk of infectious disease transmission is evaluated based on a statistical analysis of respiratory droplet trajectory distribution. An analytical model recently developed by the authors allows the prediction of the trajectory and evaporation rate of exhaled droplets. The model is used to collect data from a sampling set of more than twenty thousand droplets distributed over a range of diameters from 0.1 μm to 1 mm for different respiratory scenarios. The analytical tool implements the governing equations of droplet transport, evaporation, energy balance, and chemical composition. It also features a two-dimensional unsteady empirical model of respiratory cloud including momentum dissipation and buoyancy. A discrete random walk approach to simulate the droplet turbulent dispersion, and the randomization of the droplet release within the exhalation period and the mouth cross section area complete the model enabling statistical analyses to be rightly performed. With the due boundary conditions, different types of respiratory events can be modeled easily. With additional information on the exhaled droplet size distribution and viral content, spatial maps of virus concentration are derived and associated with the risk of infectious disease transmission being able to discriminate between various transmission routes such as fomite, airborne, or direct inhalation. Different scenarios are presented including mouth breathing, nose breathing, speaking, coughing, and sneezing. The fluid dynamic behavior of respiratory droplets is explored on a size basis, and the role of ventilation discussed. Risk evaluation provides useful information for a knowledgeable discussion on the prevention needs and means from case to case.
I. INTRODUCTION
A primary mean of host-to-host infectious disease transmission is given by pathogen-laden respiratory droplets exhaled while breathing, speaking, coughing, or sneezing. Depending on their size, spanning several orders of magnitude, and the ambient thermohygrometric conditions, droplets can either fall to the ground as they evaporate, or dry out leaving an airborne nucleus made of salts and proteins. Any virus transported by the droplet may survive up to several hours if airborne and up to several days when deposited on surfaces. During this lapse of time it can be transmitted to the airways of a susceptible host either through contact with an infected surface (e.g., fomite route) or through inhalation (e.g., airborne route) fostering the pathogenesis.
The already relatively rich and varied literature on the subject experienced a sudden boost during the most recent SARS-CoV-2 pandemic as many studies tackled the issue from different perspectives.
From a fluid mechanics point of view, a fundamental role over infectious disease transmission is played by the transport and evaporation mechanisms governing the droplet trajectory and ultimately its fate. Trajectory evaluation cannot ignore the presence of the unsteady and intermittent humid air breath cloud within which the droplet is emitted, transported, and protected from evaporation. Further, turbulence induces random velocity fluctuations perturbing the trajectory and favoring droplet dispersion.
Due to the complexity of the physics involved several studies approached the problem numerically using Computational Fluid Dynamics (CFD) following an Eulerian and/or Lagrangian approach with Discrete Phase Model (DPM) simulations. Most works focus on the analysis of specific indoor scenarios in the presence of a given ventilation strategy. For instance, in Ref. 1, the exhaled droplet dispersion due to speaking in a ventilated office room is investigated using a Lagrangian approach: the impact of different ventilation strategies on droplet dispersion is discussed. In a similar fashion, the role of ventilation in hospital rooms is analyzed in Ref. 2. In Ref. 3, droplet dispersion due to breathing or talking, or following a cough or a sneeze on a bus is dealt with: larger droplets are modeled using a Lagrangian approach, smaller ones with an Eulerian approach. In Ref. 4, the focus switches to a cough in a supermarket aisle whereas different CFD tools are compared over the same scenario. Although these studies are interesting for their practicality, they are very case specific and their results cannot be generalized easily as things could change suddenly, for instance, assuming a different arrangement of the ventilation grilles or a different relative position between an infected and a susceptible individual. Some other work remain more general using CFD to focus on the exhalations from an isolated individual in a quiescent environment. In Ref. 5, the effect of head motion and pressure variation at the mouth during a sneeze is discussed finding how the exhaled breath cloud is mostly tilted upward favoring the droplet dispersion over larger distances. In Ref. 6, instead, the role of ambient relative humidity on the evaporation of droplet from a sneeze is analyzed. The effectiveness of face masks in absorbing the larger droplets and limiting the diffusion of smaller droplets is assessed numerically in Refs. 7 and 8, and experimentally in Ref. 9.
Leaving the CFD behind, other works in the literature tackle the problem using cheaper simulation means such as codes for the dynamic modeling of the air quality indices in buildings. In Ref. 10, a Monte Carlo method is used to predict the risk of contagion in a classroom for different ventilation control strategies, while in Ref. 11, multi-objective optimization has been applied in a similar indoor environment to find trade-off solutions between energy consumption and risk of infection.
Another less computationally intensive approach found in the literature regards the creation of models where the governing equations can be solved analytically on a droplet-to-droplet basis and evaluted numerically to better address the varying boundary conditions seen by the droplet during its flight. Important aspects of these models are those already mentioned above: the evaluation of the droplet transport and evaporation, the modeling of the breath cloud, and the role of turbulent dispersion. In Ref. 12, the breath cloud is modeled as a steady non-buoyant jet while turbulent dispersion is not accounted for. The focus of the analysis is on the role of relative humidity on droplet evaporation. Similar modeling assumptions are found in Ref. 13, where a multiple shell model is proposed to account for temperature and concentration gradients inside the droplet. This allows for better predictions of the final stage of the evaporation process. In Refs. 14 and 15, the breath cloud buoyancy is added to the model through the use of empirical formulas for steady jets first proposed in Ref. 16, while the breath cloud model is further expanded by including the role of the jet potential core and assuming a self-similar non-uniform air velocity profile in the jet. In Refs. 17 and 18, turbulent dispersion is modeled with a discrete random walk approach where the air velocity fluctuations in the breath cloud are stochastic functions of the turbulent kinetic energy and of the turbulent dissipation rate. Spatial fields for these quantities are thus needed, and are usually taken from empirical formulas for steady jets found in the literature.19 The introduction of the random walk model allows stochastic analyses of the droplet spatial distribution to be performed. Only recently the steady-state jet hypothesis was abandoned considering that respiratory events are composed of a jet phase where momentum is injected in the domain followed by a puff phase once the exhalation terminates, and these two phases have a different space vs time evolution. This was investigated experimentally in Ref. 20 and translated into a numerical model in Ref. 21 starting from momentum conservation hypothesis under the assumptions of steady jet phase and uniform spatial velocity in the cloud.
In a recent work by the authors,22 the breath cloud model was further extended to address any type of unsteady and intermittent velocity input at the mouth, being thus able to render also periodic respiratory events such as breathing, and any type of inlet velocity profile, given that the profile is assumed self-similar across the expanding breath cloud in space. Although the analysis stems from momentum conservation considerations, momentum dissipation in the cloud is modeled with the introduction of a momentum decay characteristic time. Finally, the droplet ejection is randomized in time, within the time span of the exhalation jet phase, and in space, within the mouth/nose cross section, thus adding to the statistical relevance of the stochastic analyses that can be performed with such a model. The randomization in time becomes significant as the breath cloud is no longer modeled as a steady jet.
While fluid mechanics can provide useful information on the droplet trajectory dispersion, which is a quantity certainly correlated with the risk of disease transmission, this is not enough to complete the picture. Additional elements are needed for a tentative quantification.
First of all, an estimate of the exhaled droplet number and size distribution is required. A relatively large number of papers addressed the issue focusing on different types of respiratory events since the pioneering work of Duguid.23 To cite a few: speaking and coughing are analyzed in Ref. 24, speaking, whispering, mouth breathing, and nose breathing in Ref. 25, breathing, speaking, and singing in Ref. 26. Droplet sizes from a fraction of micrometer up to a millimeter are found, depending on the type of respiratory event, but on their quantity and distribution there is no agreement in the literature. A limit may reside in the experimental methods and in the instruments employed not being able to capture the size correctly over such a large range of scales. In Ref. 27, many of these studies are compared finding that the results can be reasonably overlapped if properly normalized. This leads to the assumption of a power-law distribution of exponent 2 as best compromise solution for the droplet size in any type of respiratory event.21
Second, an estimate of the virus concentration in the saliva is needed to evaluate the potential virus shedding rate. For the case of SARS-CoV-2, an average concentration of was measured in Ref. 28 even though this value can change significantly as the disease progresses. Values as high as were measured in hospitalized patients.29
As for the infectious dose likely to cause the disease, it is believed to be around ,30 a value similar to that known for influenza. Furthermore, the infectious titer decay of SARS-CoV-2 virus has been measured in 1.1, 5.6, and half-life in air, on stainless steel, and on plastic surfaces, respectively.31
Information on the quantity of exhaled droplets and on the concentration of virus in the saliva allows the virion shedding rate from an infected individual to be estimated.32 Together with knowledge about the titer decay and the fraction of airborne droplets, this allows to estimate the virus accumulation rate in a closed environment,30 eventually including the diluting effect of ventilation.33 This is usually done under the assumption of well-mixed indoor air using the Wells-Riley equation.34 Finally, with information on the infectious dose, and the exposure time and breathing rate of a susceptible individual the risk of infection can be estimated. This is conveniently done introducing the concept of quantum generation rate,35 i.e., normalizing the virus shedding rate by the infectious dose. The well-mixed hypothesis essentially restricts the analysis to the case of airborne transmission, while other transmission routes also exist.
The World Health Organization, in fact, distinguishes between five routes:36 direct contact with an infected individual, indirect contact or fomite (i.e., contact with an infected surface), direct inhalation (i.e., inhalation of exhaled droplets over short distances), indirect inhalation or airborne (i.e., inhalation of suspended droplet nuclei accumulated in an environment), orofaecal (through infected feces contaminating air, water sources, or food).
In the present work, the analytical model developed by the authors and mentioned above is used to perform a series of statistical analyses on the droplet trajectory distribution for a large number of droplets and different types of respiratory events (namely, mouth breathing, nose breathing, speaking, coughing, and sneezing). These events are typified with given boundary conditions derived from the literature. Considering that no “standard” respiratory event actually exist, a sensitivity analysis37 on the impact of the various simulation parameters on the droplet trajectories is also presented. Following the simulation campaign, with additional information on the quantity and size of droplets exhaled in a real event, and on the virus concentration in the saliva, a virus content is associated with every simulated droplet. From the analysis of the droplet trajectories and velocity histories, an estimate of the virus spatial concentration in the neighborhood of an infected individual is derived. This allows to build virus concentration maps. By discriminating between the droplets falling to the ground and those remaining airborne, also a map of virus ground accumulation can be given, together with an estimate of the indirect inhalation transmission risk following from Wells-Riley model. By introducing an inhalation model, the spatial concentration map is translated into an inhalation map predicting the viral load inhaled by a susceptible individual as a function of its relative position with respect to the infected exhaler. The inhaled viral load is then associated with the probability of developing the disease. The analytical model is thus exploited to quantify the risk of infection related to the different transmission routes associated with the exhaled droplet trajectories: namely, indirect contact, direct inhalation, and indirect inhalation.
II. ANALYTICAL MODEL
The analytical model is briefly recalled here, for an extensive discussion the reader is referred to the previous work by the authors in Ref. 22. The model includes equations for the droplet transport, evaporation, energy balance, and nonvolatile fraction. An exhalation model including buoyancy and turbulent dispersion effects completes the work providing the velocity and the thermohygrometric boundary conditions as functions of space and time.
A. Droplet transport
B. Droplet evaporation
C. Droplet energy balance
D. Droplet nonvolatile fraction
E. Exhalation model
The breath cloud is a warm humid air stream gradually mixing with ambient air by entrainment. As such, its velocity and thermohygrometric properties also change gradually. Within this cloud, respiratory droplets are exhaled. As long as the droplet finds itself within the breath cloud boundaries, an exhalation model is necessary to provide the correct and ever changing ambient conditions to the analytical model outlined above.
F. Turbulent dispersion and droplet randomization
The random walk approach allows statistical analyses on the droplet trajectory distribution to be performed. This feature is further enhanced by considering how a droplet can be exhaled at any point within the mouth cross section and at any time within the exhalation duration, and this affects its trajectory. A randomization of the droplet initial position and time has been implemented where the probability distribution of these two quantities follows the shape of the spatial velocity profile and of the mass flow rate impulse in time, respectively, meaning that a droplet is most likely exhaled where and when the velocity at the mouth is higher. The initial droplet velocity is assumed aligned with the local exhaled flow and having an initial velocity chosen as a given fraction of the local air flow velocity.
G. Inhalation model
The analytical model presented in Sec. II was validated in the previous work from the authors,22 to which the reader is referred for additional information, against the very sparse literature available. An excellent agreement was found with the experimental works on droplet evaporation.41,46,47 Droplet transport and the resulting trajectory could be successfully validated against simpler analytical works14,48 based on the emission of isolated droplets with no turbulent dispersion. To the author knowledge, in fact, no experimental work on this topic exists in the literature. The turbulent dispersion model implemented, instead, is based on a simple and classical approach widely accepted in the numerical community.
III. EXHALATION DATA
In addition to the analytical model summarized above, further information is needed to quantify the potential risk of an infectious disease transmission. This regards the exhaled droplet size distribution and quantity, the virus concentration in the saliva, the infectious dose, and the virus titer decay.
In the following, is chosen as suggested in, Ref. 21, and is arbitrary assumed due to the lack of clear quantitative information in the literature. This choice enhances the virus airborne fraction in the results that follow so that predictions of the inhalation risk are likely in favor of safety. The virus concentration in the saliva is set to , this being the average concentration found in Ref. 28 for SARS-CoV-2, while the average infectious dose is assumed following Ref. 30. By definition, the average infectious dose represents the number of inhaled viruses sufficient to induce the disease with a probability of in a susceptible individual. Considering that droplet size in the range from 0.1 μm to have been reported for respiratory events,27 from Eq. (31) follows. The value of the constant B, instead, depends on the type of respiratory event. In Ref. 27, a review on the range of droplets size and concentration measured in several works for breathing, speaking, coughing, and sneezing is presented. Despite the results of the various works compared may even differ by orders of magnitude, by averaging the results logarithmically and assuming the exhaled air volume associated with each respiratory event, the coefficient B can be estimated using Eq. (27). The result of this operation is resumed in Table I. Data in the table are to be intended per respiratory event. For breathing and speaking, an exhaled flow of is assumed with breathing period of ,44 while for coughing and sneezing exhaled volumes of and are assumed, respectively.43 The large range mentioned above for the droplet size holds for coughing and sneezing, while during milder respiratory events such as breathing and speaking larger droplets are not expelled. An upper threshold of 10 μm for breathing, and 100 μm for speaking is chosen in the current work in line with some previous study.27,52
Type of event . | Exhaled air volume . | Period . | Range of diameters . | B . |
---|---|---|---|---|
Breathing | 4 s | μm | ||
Speaking | 4 s | μm | ||
Coughing | μm | |||
Sneezing | μm |
Type of event . | Exhaled air volume . | Period . | Range of diameters . | B . |
---|---|---|---|---|
Breathing | 4 s | μm | ||
Speaking | 4 s | μm | ||
Coughing | μm | |||
Sneezing | μm |
Table II resumes the results of Eqs. (27)–(30) given the assumptions above. The results are reported for the different respiratory events and are split by droplet size order of magnitude. In the case of periodic event hourly rates are given, while for coughing and sneezing, data are to be intended per event. It can be seen how moving from one order of magnitude to the next the number of droplets is reduced by a factor of 10, the overall droplet volume grows by a factor of 100, and the number of quanta emitted increases by a factor of 10. This follows from the values chosen for m and n coefficients. By considering that in standard temperature and relative humidity conditions the critical diameter (i.e., the diameter below which a droplet remains airborne) is slightly smaller than 100 μm, it can be deduced that with this choice in the case of coughing and sneezing the airborne fraction includes of the droplets exhaled, summing up to of the overall droplet volume, and of the quanta emitted. The values of the total liquid volume exhaled reported in the table are in line with those given in mass and mentioned in Ref. 53 for coughing ( ) and in Ref. 3 for sneezing ( ) events.
Type of event . | Range of diameters . | Droplet number ( ) . | Droplet cumulative volume ( ) . | Virus number ( ) . | Number of quanta ( ) . |
---|---|---|---|---|---|
Breathing | μm | ||||
μm | |||||
Total | |||||
Speaking | μm | ||||
μm | |||||
μm | |||||
Total | |||||
Coughing | μm | ||||
μm | |||||
μm | |||||
μm | |||||
Total | |||||
Sneezing | μm | ||||
μm | |||||
μm | |||||
μm | |||||
Total |
Type of event . | Range of diameters . | Droplet number ( ) . | Droplet cumulative volume ( ) . | Virus number ( ) . | Number of quanta ( ) . |
---|---|---|---|---|---|
Breathing | μm | ||||
μm | |||||
Total | |||||
Speaking | μm | ||||
μm | |||||
μm | |||||
Total | |||||
Coughing | μm | ||||
μm | |||||
μm | |||||
μm | |||||
Total | |||||
Sneezing | μm | ||||
μm | |||||
μm | |||||
μm | |||||
Total |
IV. SIMULATION SETUP
The analytical model outlined in Sec. II has been used to perform statistical analyses of the exhaled droplet trajectory distribution for the cases of mouth breathing, nose breathing, speaking, coughing, and sneezing. The simulation setup for what concerns the ambient, the droplet initial conditions and chemical composition, and the common breath cloud features is fixed and the modeling choices are resumed in Table III. Normal temperature and pressure conditions with 50% relative humidity are assumed. Two ventilation schemes are modeled: one with no ventilation and quiescent ambient air, the other with an average air velocity of blowing from behind on the infected individual. The breath cloud is assumed at and 95% relative humidity, while the entrainment coefficient has been shown in the literature to be for any respiratory event.43,44 Exhaled droplets are assumed to be released from height, being in thermal equilibrium with the breath cloud, and being composed of a 1% volume solid fraction made of salts and proteins.
Ambient | |
Pressure | |
Temperature | 20 °C |
Relative humidity | 50% |
Air velocity | and |
Air direction | horizontal |
Breath | |
Temperature | 34 °C |
Relative humidity | 95% |
Entrainment coefficient | 0.11 |
Droplet | |
Initial temperature | 34 °C |
Mouth/nose height | |
Initial velocity and direction | Aligned with and equal to the local breath velocity |
Initial solid volume fraction | 1% |
Solid fraction density | |
Solid fraction specific heat | |
Solid fraction molar mass |
Ambient | |
Pressure | |
Temperature | 20 °C |
Relative humidity | 50% |
Air velocity | and |
Air direction | horizontal |
Breath | |
Temperature | 34 °C |
Relative humidity | 95% |
Entrainment coefficient | 0.11 |
Droplet | |
Initial temperature | 34 °C |
Mouth/nose height | |
Initial velocity and direction | Aligned with and equal to the local breath velocity |
Initial solid volume fraction | 1% |
Solid fraction density | |
Solid fraction specific heat | |
Solid fraction molar mass |
Additional boundary conditions are needed to define the breath cloud behavior and these obviously depend on the type of respiratory event to be modeled. Their value is resumed in Table IV. In the literature, there is a limited agreement on the outflow velocities associated with different respiratory events, and it is rarely specified whether the values reported are peak or average ones. The values in Table IV derive from an attempt to critically analyze and average all these data. For breathing, for instance, a value of is reported in Ref. 54, while values between 0.5 and are found in Ref. 49 for light to heavy breathing. For coughing peak velocity is reported in Ref. 55, while average velocity in Ref. 15, with peaks of mentioned in Ref. 53. For sneezing, most common values are between (Ref. 56) and for the peak velocity,55 whereas a peak volume flow rate of is given in Ref. 57. Time durations of – are reported for coughing and sneezing43,55 with peak velocities found after from the beginning of the exhalation. For breathing, 15 cycles per minute are usually assumed,58 with the velocity peaking at .44 Mouth and nose diameters in Table IV are derived to match the exhaled air volume and the event time span with the average velocity expected. Since, in the present work, the breathing function in Eq. (17) is assumed, is always attained.
Type of event . | Opening diameter ( ) . | Breath direction (θ) . | Exhaled air volume ( ) . | Breathing function times ( ) . | Opening cross section . | Avg. volume flow rate ( ) . | Average and peak velocity ( ) . | Dissipation time ( ) . |
---|---|---|---|---|---|---|---|---|
Mouth breathing | −5° | |||||||
Nose breathinga | −60° | |||||||
Speaking | −5° | |||||||
Coughing | ||||||||
Sneezing |
Type of event . | Opening diameter ( ) . | Breath direction (θ) . | Exhaled air volume ( ) . | Breathing function times ( ) . | Opening cross section . | Avg. volume flow rate ( ) . | Average and peak velocity ( ) . | Dissipation time ( ) . |
---|---|---|---|---|---|---|---|---|
Mouth breathing | −5° | |||||||
Nose breathinga | −60° | |||||||
Speaking | −5° | |||||||
Coughing | ||||||||
Sneezing |
Data per nostril, the two streams are separated by 20° circumferentially on the horizontal.
According to Refs. 49 and 59, the direction of the exhaled breath is slightly tilted downward in the case of mouth breathing, while the angle with respect to the horizontal for nose breathing is (Refs. 44 and 60) and the streams coming from the nostrils are apart. A cough is usually directed downward by ,43 while during sneezing the direction of the exhalation changes in time due to the head motion,5 but on average is directed upward by . The breath cloud momentum dissipation characteristic time is set to for breathing and speaking, for coughing, and for sneezing: from preliminary CFD analyses of an exhaling individual in an empty room these appear to be reasonable choices providing a good match with regard to the extents of the breath cloud in time and space for the different respiratory events.
Although the rationale behind the modeling choices resumed in the tables and presented so far is rather reasonable, it must be highlighted that most of these data are subject to large margins of uncertainty for a number of reasons ranging from the difficulties in carrying out accurate measurements, to the lack of experimental data in the literature, to their inherent variability with respect to many factors. More simply, a standard breath, cough, or sneeze is just an abstraction that yet is needed to quantitatively try to investigate different scenarios in terms of infectious disease transmission risk.
V. POST-PROCESSING
A quanta load is associated with each of the 972 diameter ranges in the same way as shown in Table II for wider ranges, and to each simulated droplet a load is assigned. In this way, the viral distribution in each simulation campaign is modeled coherently whatever the size distribution and number of simulated droplets. It is therefore sufficient that the number of runs is large enough for the trajectory distribution to be statistically relevant and cover the space in a sufficiently uniformly way.
The same Eqs. (42)–(47) hold also for estimating the concentration on the ground, expressed in , given that the airborne quanta term is substituted by the amount deposited on the ground, the room volume is substituted by the surface area, and the characteristic time is computed omitting the ventilation term. To be mentioned that virus half-life on the ground depends on the surface type and its value is generally much higher than that found in air. Further, it must be considered that quanta deposition on the ground is not uniform, i.e., the well-mixed hypothesis technically does not hold. Nonetheless, the equations can still be applied to small surfaces whose concentration may well be assumed locally uniform: in this case the terms and Q/S must be intended as local viral deposition rate and local initial concentration, respectively.
In view of airborne infection risk assessment, the average volumetric concentrations over the time span of interest for periodic [Eq. (45)] and isolated [Eq. (47)] respiratory events must be multiplied by the volume inhaled [Eq. (35)] and the result fed into Eq. (37) to compute the quanta inhaled by a susceptible individual and the associated risk of infection, respectively. In virtue of the well-mixed hypothesis, the result of this operation consists in a uniform concentration or in a certain amount of inhaled quanta to be added to the quantities already computed for the direct inhalation transmission route.
With regard to fomite infection risk assessment, instead, the instantaneous spatial concentrations for periodic [Eq. (43)] and isolated [Eq. (46)] respiratory events, actually retrace Eqs. (39) and (38) enriching them with time-dependent terms that allow to address the role of virus titer decay on the surface concentration.
In the simulations performed, the droplet trajectory is tracked up to when the droplet reaches the ground or, in the case of airborne particles, up to when evaporation terminates, thermal equilibrium with the environment is reached, and a terminal velocity below a given threshold of is attained outside the area of influence of the breath cloud. When this occurs, it is assumed that the airborne droplet nuclei will be transported by the smallest air current, and in the long run will spread uniformly in the room. Information on the droplet trajectory is used to compute the viral inhalation associated with direct inhalation transmission route, information on the droplet deposition is used to evaluate the potential viral inhalation associated with fomite transmission route, and information on the airborne droplets is used to quantify the viral inhalation associated with airborne transmission route.
In the following, the airborne viral concentration is evaluated for a rather severe scenario of t = 1 h exposure starting from a virus free environment at time zero, in an average size ( ) poorly ventilated (1 air change per hour) office, assuming a virus half-life in air of as found in Ref. 31 for SARS-CoV-2. With these data, the room characteristic time computed with Eq. (44) results of . The viral concentration at the ground, instead, is evaluated for the same scenario at time , given a virus half-life on the surface of (i.e., ), which is an intermediate value among those mentioned in Ref. 31 for surfaces.
For post-processing purpose, the space in front of the infected individual is meshed with hexahedral cells of a given size (namely, ). In each cell quanta concentration or viral exposure associated with direct inhalation, transmission route is computed using Eqs. (33) and (34), thus providing practical viral concentration/exposure spatial maps. By further applying the inhalation model in Sec. II G to each cell center, similar viral inhalation maps are extracted from Eqs. (35) and (36). From these maps, applying Eq. (37), the risk for a susceptible individual of contracting an infectious disease when exposed to the exhalations from an infected individual can be quantified. The uniform airborne concentration after one hour exposure is then evaluated for the office scenario just mentioned from Eqs. (45) and (47), and its contribution added to the maps. The concentration in Eq. (47) is made dimensionally consistent with the viral exposure in Eq. (34) multiplying it by the exposure time.
For visualization needs, since these maps are three-dimensional, the spatial maps that follow are presented on a vertical x–z plane aligned with the direction of the exhalation where the concentrations in the spanwise direction y are cumulated and averaged over the breath cloud width at the given longitudinal distance from the mouth.
VI. RESULTS
The results of the simulation campaigns described in Sec. IV, including five types of respiratory event and two room ventilations, are presented in this section. The results are given in terms of droplet trajectories, viral concentration/exposure, inhalation, and risk of infection maps on the x–z vertical plane, and in terms of droplet deposition and viral concentration/deposition rate on the x–y ground plane. Animations of the instantaneous viral concentration in time can be found in the multimedia material available online and associated with the present work.
A sensitivity analysis on the effect of the main simulation parameters on the trajectories of droplet of different size is given in the end. This allows to appreciate how changes in the ambient conditions or in the intensity of the respiratory event may affect the results that, out of necessity, are given for a discrete, although large, number of scenarios.
A. Breathing
From Table II, it is noted that the virus exhalations associated with a breathing infected individual is relatively small and has been estimated in . The droplets exhaled are small enough (see Table I) to be fully transported by the breath cloud, so that the airborne fraction equals 100%. Following Eq. (45), this results in an average airborne viral concentration in the room of after one hour exposure, which translates into a minimum of inhaled by a susceptible individual breathing normally, and corresponding to a risk of infection of 0.14%. The difference between mouth and nose breathing scenarios is limited to the direction and the velocity of the exhalation as summarized in Table IV.
Figure 2 (Multimedia view) resumes the results obtained for the mouth breathing scenario without room ventilation. Figure 2(a) shows the exhaled droplet trajectories that intertwine and overlap within the breath cloud as a result of turbulent dispersion and the different initial conditions. Note that smaller droplet trajectories are drawn on top of larger droplet ones so that large droplet data are mostly covered due to the large number of trajectories at stake. In this scenario, the droplets are so small that the airborne stopping criteria in the simulations are met almost immediately once the droplet exits from the breath cloud boundaries. This may either occur due to the fading of the breath cloud momentum with the distance or to turbulent dispersion pushing the droplet away from the cloud prematurely. The effect of buoyancy, although marginal, is apparent from the curved trajectories. Applying Eq. (33) to each cell, trajectory information is transformed into the viral concentration map in Fig. 2(b). Of course viral concentration is higher close to the mouth of the infected individual and dilutes with the distance. A non-negligible concentration is perceived within the breath cloud extension up to a distance of even for the seemingly harmless mouth breathing scenario. Applying the inhalation model in Sec. II G to the concentration data, the viral inhalation rate in Fig. 2(c) is found, and further applying Eq. (37) to the inhalation rate the risk map in Fig. 2(d) follows. It is worth mentioning how to read these figures: in Fig. 2(c), for instance, each cell is colored according to the predicted viral inhalation rate a susceptible individual breathing normally would be subject to while staying with his mouth or nose at that particular location in space and finding himself in front of (i.e., aligned along the y direction with) a mouth breathing infected individual. The figure accounts both for direct inhalation and airborne transmission routes for the generic condition mentioned in Sec. V. The viral concentration associated with the latter is identified by the lower scale of the color bar and being orders of magnitude lower than that related to the peaks of the direct inhalation route it does not confound its reading. The same considerations hold for all the three maps in Figs. 2(b)–2(d) and for analogous maps presented in the following. The three maps look similar as they represent further elaborations of the same set of data given by the viral content of the single droplets and their trajectories as resumed in Fig. 2(a). The plume appears thicker in Fig. 2(c) since the inhalation rate, following the model in Sec. II G, is essentially a conditional sum of the concentrations over a radius hemisphere.
Figure 3 (Multimedia view) shows the same set of results for the case of nose breathing. Compared to the previous case, the difference in the results is marginal, if not for the different direction assumed by the trajectories. The results of the ventilated breathing scenarios are omitted for brevity. The marginal impact of ventilation on the breath cloud extension is such that the resulting maps are largely similar to those already shown.
B. Speaking
According to Table IV, the boundary conditions adopted for the breath cloud modeling in speaking and mouth breathing scenarios are the same. This choice was made considering that, on average, the breathing period and the overall volume of air exhaled while speaking are indeed similar to those characterizing breathing. Nonetheless, the continuous movement of the lips while speaking will induce irregular air outflow velocity patterns having sharp peaks on plosive sounds.58 As these patterns much depend on the type of sounds emitted and cannot be generalized, the choice fell on modeling the average behavior, knowing that in reality droplet dispersion will likely be more pronounced than predicted. The main difference between the mouth breathing and speaking models adopted here is that the latter is characterized by a double coefficient B in the droplet size PDF (see Table I) that translates in double exhalations for a given diameter range. Further, the droplet size range of interest extends up to 100 μm, resulting in a twentyfold quanta exhalation rate of compared to breathing, as summarized in Table II. While the results relative to droplets up to 10 μm in size are the same already seen for mouth breathing scenario, larger droplets are big enough to escape from the breath cloud due to gravity, eventually reaching the ground. After the simulations, a critical diameter of μm is found for both ventilation scenarios, so that the quanta airborne fraction equals 92% of the total ( ). This translates into an average airborne viral concentration in the room of after one hour exposure and a minimum of inhaled by a susceptible individual breathing normally. This corresponds to a risk of infection of 2.6%: a value much higher compared to that of the breathing scenario.
Figure 4 (Multimedia view) shows the results of the speaking scenario without room ventilation. Larger droplets are still small enough to have little inertia so that their horizontal velocity is quickly zeroed by drag once they exit from the breath cloud area of influence. As a consequence, their trajectories are mostly vertical rather than parabolic [see Fig. 4(a)]. The larger the droplets and the sooner they exit from the breath cloud, the shorter the horizontal distance traveled. In this scenario, for instance, supercritical droplets are deposited on the ground no further than from the infected individual, as shown in Fig. 4(e), and their dispersion in the y direction is very limited. Smaller droplets can be transported further, and if their size is μm they are able to escape from the breath cloud by gravity at some point along their journey. Overall, the viral concentration associated with these droplets can be relatively high as shown in Fig. 4(b), but the horizontal distance traveled by them is small enough not to represent a serious threat to the safety of a susceptible individual maintaining a reasonable social distance [see Figs. 4(c) and 4(d)]. It is rather of concern the quantity of airborne viruses that may accumulate with time in a poorly ventilated environment and against which there is nothing social distancing can do. The deposition rate map in Fig. 4(f) is derived from the droplet deposition data in Fig. 4(e). The map is characterized by several local peaks due to the relatively limited number of droplets touching the ground (only 189 in this scenario) compared to the relatively large number of cells with which the ground is discretized. To obtain a more uniform and representative map, data in each cell is iteratively smoothed by averaging its value with that of the neighboring cells. The result of this process is shown in Fig. 4(f). To be noted that the random-walk turbulent dispersion model adopted randomly deviates the droplet trajectory so that they may fall either in the y positive or negative half-plane indifferently. Considering the problem symmetry, all deposition coordinates are projected in the y-positive half-plane in Fig. 4(e) and the droplet viral content is split equally among the two half-planes to obtain the symmetrical plot in Fig. 4(f). To be also noted that the pixelated look that this and other maps might have is a consequence of the discretization used to create them. This effect is particularly evident where the extension of the map is limited in size. Nonetheless, finer discretizations would have led to less uniform and less significant maps, unless the number of simulations was increased exponentially.
The deposition rates in the figure are sufficient to justify a significant risk of fomite transmission in case an infected surface is touched and its viral load brought in contact with the airways of a susceptible individual. In total, the viral load deposited on the ground amounts to and from the simulations the infected area amounts to . From Eq. (43), after one hour of exposure the average viral concentration on the ground amounts to , and the risk of infection if a single square centimeter is touched and its viral load fully inhaled equals 6.2%. Further quantitative estimates can be derived easily from Eq. (37) as a function of the deposition time span, and thus of the actual viral concentration on the surface, and the extension of the contact. It is worth mentioning that the virus half-life on surfaces is generally much higher compared to air, and according to the data provided in Ref. 31 for SARS-CoV-2 the characteristic time associated with the virus titer decay on a surface can be as high as . This means a surface may still represent a non-negligible threat in view of fomite transmission even a few days after infection.
Trajectories are much different in the ventilated scenario shown in Fig. 5 (Multimedia view). Droplets are carried by the air stream, and having different terminal velocities depending on their size, they stratify in the vertical direction creating a compact front advancing in time as highlighted in Fig. 5(a) and in the multimedia material provided. As a result, the risk of infection [Fig. 5(b)] remains considerable also at larger distances from the infected individual, and at any height. A clear stratification is also apparent in the droplet deposition plot [Fig. 5(c)] where slightly smaller droplets are carried much further by the air flow before being deposited on the ground. The average deposition rate in Fig. 5(d) is lower compared to the previous non-ventilated case due to the larger area ( from the simulations) over which the same amount of virus is spread. Note that no droplet is deposited within a distance of from the infected individual since the falling time of larger droplets is and the air velocity . These results highlight the pros and the cons of ventilation. Ventilation, in fact, is generally associated with higher air change rates in a room, and as such it is most welcome in that it helps diluting the virus concentration lowering the risk of airborne transmission. On the other hand, if a susceptible individual accidentally finds himself downstream of an infected one, the chances of transmission by direct inhalation are largely increased. For what concerns fomite transmission, in the presence of ventilation the chance of inhaling viruses is higher due to the larger droplet dispersion on surfaces, but the risk of disease transmission is lower thanks to the greater virus dilution.
C. Coughing
Coughing, together with sneezing, is an isolated respiratory event characterized by the rapid exhalation of a large amount of droplets. This is reflected both in the values of the coefficient B in Table I, that for coughing is more than 300 times larger compared to speaking, and in the fivefold exhalation velocity reported in Table IV. Overall, the number of exhaled quanta associated with a single coughing event has been estimated in 172.8 in Table II. After the simulations, a critical diameter of μm is found for both ventilation scenarios, a value slightly lower than the one found for speaking. This is justified by the direction of the exhalation which is tilted downward and thus favors droplet deposition. The quanta airborne fraction equals 9% of the total and amounts to per event. This suggests that, for what concerns the airborne viral load, one cough equals min speaking. This quantity translates into an average airborne viral concentration in the room of in the hour following the respiratory event according to Eq. (47), inhaled by a susceptible individual breathing normally, and a risk of infection of 1.5% after 1 h exposure.
Figures 6 and 7 (Multimedia views) show the results of the coughing scenarios. The latter including room ventilation, the former without. Figures 6(a) and 7(a) show the droplet trajectories. Qualitatively, it is possible to discriminate between four types of droplet behavior: small droplets ( μm) whose trajectory is primarily due to their transport within the breath cloud and are affected by buoyancy, sub-critical intermediate droplets that remain airborne and whose trajectory is mostly governed by the local ventilation flow, supercritical intermediate droplets reaching the ground before evaporating and whose trajectory is affected by the local air flow to a lesser extent, and large droplets ( μm) whose trajectory tends to be parabolic and is not affected by air transport. Droplet dispersion on the ground follows a conic trace, as shown in Figs. 6(e) and 7(c), which is characterized by a slightly larger opening angle compared to that of the breath cloud ( vs ). This shows how the impact of turbulent dispersion is marginal on supercritical droplets. Thanks to the larger velocities with whom droplets are exhaled and transported within the breath cloud, the horizontal distance reached is larger compared to previous scenarios, and with that also the risk of infection at a distance becomes meaningful. Small and intermediate size droplets form a rather compact front characterized by high viral exposure that fades with the distance [Figs. 6(b)] and is perceived up to distance and more in the case of ventilated scenario. Large droplets, mostly falling to the ground at a distance between 1 and from the infected individual, further extend the risk area [Figs. 6(d) and 7(b)]. The role of ventilation on the transport of intermediate droplets is apparent from Figs. 7(c) and 7(d) compared to their non-ventilated counterparts.
Despite the amount of quanta deposited on the ground after a cough is forty times larger ( ) than that estimated after 1 h talk, much larger is also the area over which droplets are spread ( ), so that after one hour from the coughing event the average viral concentration on the ground, estimated with Eq. (46), is just half ( ) and so is the risk of infection after touching and inhaling the viral content of one square centimeter (3.4%). To be considered that these are just average data, while the actual viral concentration on the floor varies significantly from place to place and peaks can be six times larger than that or more. Figure 6(f) provides a more detailed picture of this. In the ventilated scenario, the infected area is slightly reduced as small droplets are transported further and none fall within of the infected individual [see Figs. 6(e) and 7(c)].
D. Sneezing
Data in Table IV show that sneezing modeling is characterized by triple exhalation velocity compared to coughing and its direction is tilted upward. This means droplets are expected to cover a much wider range of distances in their trajectories. Furthermore, the coefficient B in Table I is 10 times larger, and so is the amount of virus exhaled that reaches per event (see Table II). A critical diameter of μm has been found in the simulations, meaning that the quanta airborne fraction is 10% of the total and equals per event. The average airborne viral concentration in the room in the hour following the respiratory event is as high as , the quantity of airborne virus inhaled amounts to , and the risk of infection after one hour exposure spikes at 15.5%.
Figure 8 (Multimedia view) shows the results of the sneezing scenario without ventilation. What catches the eyes looking at the trajectories in Fig. 8(a) is the large horizontal extent covering a distance up to 6 m from the infected individual. Overall, results are qualitatively similar to those already shown for coughing where the different behavior of droplets of different size can be recognized, and the droplet dispersion in the y direction shown in Fig. 8(e) is of the same entity. The main differences are in the direction of the breath cloud, and in the tenfold quanta exhalation that, although dispersed over a much larger volume of space, translates into an increased viral exposure even at a large distance [Fig. 8(b)] and a risk of infection above 90% over most of the area covered by the droplet flow [Fig. 8(d)]. As a result, according to these values a susceptible individual eventually hit by the exhalations from a sneeze has a high chance of directly inhaling enough viruses to develop the infectious disease. These numbers may seem too large, and maybe they are considering the level of uncertainty on the information about the viral load in droplets and about the infectious dose that can be found in the sparse data available in the literature. Nonetheless, it is certainly true that in this scenario a distance of 6 m could at best be considered barely sufficient to avoid the risk of direct contagion.
Results of the ventilated scenario are omitted for brevity in that they only differ in the distance traveled by intermediate droplet in a similar fashion to what already seen in the coughing scenario.
To be noted how the maps in Fig. 8 appear less smooth compared to those obtained in previous scenarios, and the trajectories associated with the single droplets can be perceived in the maps, in particular for what concerns large droplets. Even though the number of droplet simulated is very high, in fact, this number is not enough to cover uniformly such a large portion of space. With the simulation of additional droplets the maps could be made smoother and a more gradual fading of the risk with distance would appear. The considerations made on the safety distance to be kept with respect to a sneezing infected individual, nonetheless, would not change significantly. The map that suffers the most from the limited number of droplets simulated is the one showing the viral concentration at the ground in Fig. 8(f) due to the small ratio between the number of droplets reaching the ground (5161 in this scenario) and the number of pixels with whom the ground is discretized. In order to obtain a more continuous and plausible map in the post-processing of Fig. 8(f) the cell size has been increased to . The infected area reaches an extension of , with an average viral concentration at the ground 1 h after exposure of , and a risk of infection of 4.2% in case the viral content of one square centimeter is inhaled. These values are similar to those already found for coughing since even though the viral load is tenfold, also the infected surface grows nearly proportionally.
E. Airborne and fomite transmission
The results shown so far in Sec. VI are very informative in that they provide an estimate for different infectious disease transmission routes.
As mentioned, maps on the x–z vertical plane are associated with direct inhalation transmission route for the various scenarios analyzed and, to the author's knowledge, they are a first attempt in the literature to provide spatial-dependent quantitative estimates for this route. The estimates given rely on the analytical model summarized in Sec. II. The scenarios investigated may either refer to a periodic respiratory event (e.g., breathing or speaking) where viruses are exhaled periodically in time generating a virus concentration field which is steady on average, or to an isolated respiratory event (e.g., coughing or sneezing) where viruses are exhaled over a short lapse of time.
Along with this, also airborne transmission needs to be accounted for. Following the Wells–Riley model and assuming a well-mixed environment, airborne transmission risk is easily computed. Nonetheless, it must be considered that in this case the risk of transmission depends on the airborne virus concentration in the environment, which is a quantity that changes with continuity in time. In particular, the airborne spatial concentration grows with time in the case of periodic event due to the progressive accumulation of airborne droplet nuclei, while it decays with time in the case of isolated event due to virus titer decay and the progressive dilution by ventilation. In order to match these results, that are clearly dependent on the exposure time, with those of the maps focusing on direct transmission where the concentration is either steady or the viral exposure limited in time, the average concentration associated with airborne transmission has been computed for the various scenarios for a given time lapse going from zero to 1 h, as mentioned in Sec. V. The resulting quantities have been added to the maps and are represented by the background value reported at the bottom of each color bar.
Even though airborne concentration is small compared to the peaks found in the maps, its relevance is connected to the fact that airborne transmission may occur over large exposure times and despite the social distance that may be kept. It is, thus, worth focusing on how the risk of airborne transmission does change with the exposure time along with periodic infectious exhalations or following an isolated respiratory event. This aspect is investigated in the current section.
The considerations just made also hold for fomite transmission since the average virus concentration on the surface depends on both the viral load discharged on it (either periodically or in a single event) and on the virus titer decay on the surface, i.e., on the time passed since the event. Further, the risk of infection depends on the extension of the surface touched and whose viral load is accidentally inhaled. In this case, of course, the well-mixed assumption only gives an estimate of the average risk, knowing that the concentration is not uniform in space as shown in the x–y ground maps provided. The ground concentration maps in Secs. VI B–VI D show the picture at a given moment in time, chosen as either 1 h from the beginning of a periodic event, or after an isolated event. By accepting the well-mixed hypothesis, an estimate of how the average fomite transmission risk changes in time can be given.
The governing equations are those reported in Sec. V [Eqs. (42)–(47)]. These have been used to extract plots of the time-dependent virus concentration and the risk of infection associated with airborne and fomite transmission routes for the scenarios investigated in the current work. The results are shown in Fig. 9. Table V resumes the number of quanta exhaled associated with each respiratory event subdividing the data into airborne and deposited quantities as found in the simulations. For the evaluation of airborne concentration, the poorly ventilated office scenario described in Sec. V is assumed. Similarly, for the evaluation of ground concentration the infected surface areas found in the simulations for the non-ventilated cases are taken.
Respiratory Event . | Airborne quanta . | Deposited quanta . | Infected surface . |
---|---|---|---|
Breathing | |||
Speaking | |||
Coughing | 15.02 | 157.77 | |
Sneezing | 168.22 | 1559.66 |
Respiratory Event . | Airborne quanta . | Deposited quanta . | Infected surface . |
---|---|---|---|
Breathing | |||
Speaking | |||
Coughing | 15.02 | 157.77 | |
Sneezing | 168.22 | 1559.66 |
Figure 9(a) shows the airborne average viral concentrations a susceptible individual would be subject to while staying in an infected room for a certain amount of time. To be noted that, as the characteristic time is the same for all the events, the concentrations associated with periodic events show the same trend and scale linearly with the airborne quanta, and also those associated with isolated events decay with the same rate. Data in Fig. 9(a) are processed through Eqs. (35) and (37) to provide the risk plot in Fig. 9(b). While the airborne transmission risk associated with breathing or to a single cough is pretty much negligible, it is shown how sneezing and prolonged speaking in poorly ventilated environments may spark high risks of infection.
For the same reason noted above, also the decay rate of the viral concentration on the ground is the same for every respiratory event and amounts to approximately one order of magnitude every twenty hours. Figure 9(c) shows the trends of the ground viral concentration in time for different respiratory events. These data are then recasted in terms of conditional risk of infection using Eq. (37) and assuming the viral load of a infected surface is inhaled. The two figures show that a small and comparable risk of fomite infection may follow from both coughing and sneezing: even though the viral load is high, it is spread over large surfaces. The risk associated with prolonged speaking instead may become relevant due to the limited area over which the virus is spread, at least under the hypothesis that the speaker does not move while speaking. Nonetheless, it must be considered that fomite infection risk may be largely limited adopting a proper behavior (e.g., with frequent hand washing and avoiding bringing the hands close to the airways), whereas the effectiveness of face masks against airborne transmission is only marginal.
F. Sensitivity analysis
The boundary conditions used in the previous scenarios are those resumed in Tables III and IV. In view of infection risk assessment, it is worth investigating to what extent the trajectories of the droplets are affected by the change in the simulation parameters. To do so, the horizontal distance traveled by a droplet before reaching the ground or before becoming airborne has been chosen as representative metric. Additional simulations are performed for each non-ventilated scenario in which the randomization elements described in Sec. II F are turned off, and droplets of different size are released from the breath cloud center ( ) at the velocity peak time . In this way, the horizontal distance found is close to maximum. Simulations are then repeated by changing one parameter at a time over a range of , where applicable, with respect to the values in the tables. All parameters are included in the analysis with the exception of the ambient ventilation that was set to zero in this analysis, and the times of the breathing function.
It is found that all parameters related to the droplet solid fraction and thermophysical properties are irrelevant. This is an expected results in that the behavior of large droplets is not affected by a solid fraction that is anyhow small, while airborne droplets are mostly transported by the breath cloud whatever their size and solid content. All quantities directly related to the local air density and the droplet evaporation rate (i.e., ambient pressure, temperature, and relative humidity) play a minor role only on the trajectories of large droplets, in particular for what concerns the sneezing scenario. Lower air densities in the ambient and slower evaporation rates, in fact, are associated with lower drag forces and larger droplet inertia, and translate into larger distances traveled. Nonetheless, small droplets have too little inertia to be affected significantly by the small change in density. The same occurs with the height at which the droplets are released: only large droplets characterized by nearly parabolic trajectories extend the range of the distance traveled when released from higher altitudes while smaller droplets are not affected by this parameter.
The most relevant impact on the droplet trajectories comes from parameters that are either directly related to the exhalation velocity (e.g., mouth diameter, exhaled volume, or initial velocity with respect to the cloud velocity) or to the rate at which the breath cloud velocity decreases with the distance (e.g., entrainment coefficient and momentum dissipation characteristic time). The latter mostly affect small droplets, the former large droplet trajectories. Furthermore, the breath cloud direction of course plays an important role on all droplet sizes.
VII. CONCLUSIONS
An analytical model recently developed by the authors22 and simulating the transport and the evaporation process of respiratory droplets has been extensively used to replicate different scenarios and evaluate the risk of infectious disease transmission from an infected individual. Beside the governing equations, the analytical tool implements a two-dimensional unsteady model of the buoyant breath cloud, turbulent dispersion through a random walk approach, and randomized droplet exhalation. These stochastic features allowed statistical analyses on the droplet trajectory distribution to be performed based on the simulation of a large number of droplets. With additional information on the number and the size distribution of the exhaled droplets and on the viral content of the saliva, a viral load has been associated with each droplet. By analyzing the time evolution of the droplet trajectories, these data have been translated into practical viral concentration/exposure spatial maps of the neighborhood of the infected individual. A simple inhalation model was proposed to reformulate these maps in terms of quantity of viruses potentially inhaled by a susceptible individual and ultimately in terms of the resulting risk of infection. As such, the maps presented relate to the risk by direct inhalation route for an individual staying in the surroundings of an infected person for a certain amount of time. In addition, discriminating between the droplets remaining airborne and those falling to the ground by gravity, the risks associated with airborne and fomite transmission routes were quantified for a reference indoor environment using the Wells–Riley model. Virus accumulation maps on the ground were also derived from the simulation data.
The scenarios investigated include different periodic or isolated respiratory events such as mouth breathing, nose breathing, speaking, coughing, and sneezing. Simulations are repeated for the cases of quiescent and ventilated environment. A sensitivity analysis quantifying the effect of the variability of the scenario boundary conditions on droplet trajectories is also included. To the author's knowledge this is the first work where a similar analytical model is applied to systematically investigate different real-life scenarios quantifying the associated risks of transmission and being able to discriminate between the different transmission routes. The analytical model has the advantage of being cheap in terms of computational resources required and is very precise in evaluating the droplet trajectories compared to numerical methods that rely on a discretization in space and time that is necessarily approximate due to the scales of the problem at stake. On the contrary, the level of detail with whom the unsteady breath cloud air flow can be modeled using an analytical model, although elaborated, is necessarily lower compared to what can be achieved with expensive 3D unsteady numerical methods such as CFD. Concerning this, it is expected that in a real-life scenario smaller droplets will likely be subject to increased dispersion compared to what predicted in the present work. Another limit of the present study, at least for what concerns the post-processing of the results, is in the fact that in the literature there is little agreement on the quantity and the size distribution of the droplets exhaled during the different types of respiratory event. Further, the viral content in the saliva varies during the disease and also with the type of disease, and the infectious dose that is sufficient to transmit the illness is often known only approximately. For sure viral concentration, exposure, and inhalation maps will scale linearly with the total viral load delivered with the respiratory event, while assuming different droplet size distributions may weaken or strenghten the relative importance of certain areas in the maps.
All these uncertainties may affect the results from a quantitative point of view, meaning that the values presented in the maps might need to be shifted but not to be altered in their shape and with respect to the range of distances reached by the droplets. For the sake of presenting plausible results for potential risk scenarios, these uncertain quantities have been chosen at best, after a thorough analysis of the literature available on the topic for what concerns the recent SARS-CoV-2 pandemic.
The maps presented have been shown to be very informative, enabling the investigation of respiratory droplets and infectious disease transmission to a great degree of detail from a fluid mechanics perspective. Of course, any further medical or clinical observation is beyond the scope of the present work. The results allowed a comprehensive analysis on the behavior of droplets of different size, ranging from small ones mostly transported by the breath cloud, to intermediate size droplets that are easily carried by the local air flow associated with ventilation, and to large ones whose trajectory depends primarily on the initial velocity with whom they are exhaled. Ventilation has been shown to stratify intermediate size droplets creating compact advancing fronts of high risk areas that could be dangerous for a susceptible individual accidentally finding himself downstream of an infected one. This is apparent in the speaking scenario where otherwise the droplets would not be transported far from where they have been generated. The social distancing rule we all became acquainted with during the last pandemic have been shown to be sufficient only to prevent direct inhalation from a breathing or a speaking individual given that no strong ventilation is present. The large viral load and the large distances traveled by respiratory droplets associated with stronger respiratory events such as coughing and sneezing would call for an increased social distancing that in most environments is simply not possible to be kept. Moreover, these types of events also imply non-negligible risks of airborne and fomite infection in indoor settings. As a consequence, in these scenarios other means of prevention aimed at impeding large quantities of respiratory droplets from diffusing freely in the ambient, such as the use of face masks or paying attention to coughing or sneezeing covering the mouth with the sleeve, become of the utmost importance. Prolonged speaking can also induce large concentrations of virus to gradually accumulate in a small poorly ventilated indoor environment thus fostering airborne transmission, while normal breathing is associated with minor, although non fully negligible, levels of risk.
AUTHOR DECLARATIONS
Conflict of Interest
The authors have no conflicts to disclose.
Author Contributions
Marco Cavazzuti: Conceptualization (lead); Data curation (lead); Formal analysis (lead); Funding acquisition (lead); Investigation (lead); Methodology (lead); Project administration (equal); Resources (equal); Software (lead); Supervision (equal); Validation (lead); Visualization (lead); Writing – original draft (lead); Writing – review & editing (lead). Paolo Tartarini: Funding acquisition (supporting); Project administration (equal); Resources (equal); Supervision (equal).
DATA AVAILABILITY
The data that support the findings of this study are available from the corresponding author upon reasonable request.
NOMENCLATURE
- b
-
Distance from centerline (m)
-
Droplet size PDF (1/m)
-
Drag coefficient
-
Viral concentration (1/m3; 1/m2)
-
Specific heat at constant pressure (J/kg K)
-
Drag force (N)
- D
-
Diameter (m)
-
Mass diffusivity (m2/s)
-
Viral deposition rate (1/m2 s)
-
Viral exposure (s/m3)
- g
-
Gravitational acceleration (m2/s)
- H
-
Height (m)
-
Latent heat of vaporization (J/kg)
-
Viral inhalation
-
Average infectious dose
-
Infectious dose PDF
- K
-
Evaporation rate (m2/s2)
- k
-
Turbulent kinetic energy (m2/s2)
- M
-
Molar mass (kg/kmol)
-
Mass flow rate (kg/s)
- N
-
Number
-
Nusselt number
- P
-
Momentum (Ns)
- p
-
Pressure (Pa)
- Q
-
Quanta
- R
-
Specific gas constant (J/kg K)
- r
-
Radius or radial distance (m)
-
Risk of infection
- S
-
Surface (m2)
- Sh
-
Sherwood number
- s
-
Position (m)
- T
-
Temperature (K)
- t
-
Time (s)
- u
-
Velocity (m/s)
-
Velocity fluctuation (m/s)
- V
-
Volume (m3)
-
Weight force (N)
- y
-
Molar fraction
- z
-
Vertical direction (m)
Greek symbols
- β
-
Entrainment coefficient
- γ
-
Cloud velocity scaling factor
- ε
-
Turbulent dissipation rate (m2/s3)
- ζ
-
Centerline dimensionless distance
-
Normal distributed random vector
- θ
-
Breath cloud direction (rad)
- λ
-
Thermal conductivity (W/mK)
- ξ
-
Radial dimensionless distance
- ρ
-
Density (kg/m3)
- τ
-
Characteristic time(s)
-
Volume fraction
- ω
-
Mass fraction
Subscripts
- a
-
Dry air
- ab
-
Airborne
- ac
-
Air change
- avg
-
Average
- c
-
Cell
- D
-
Exhaled droplets
- d
-
Drag
- ds
-
Dissipation
- e
-
Evaporation
- ep
-
Exposure
- ex
-
Exhalation
- gr
-
Ground
- hl
-
Half-life
- in
-
Inhalation
- lq
-
Liquid fraction
- m
-
Humid air mixture
- max
-
Maximum
- min
-
Minimum
- mt
-
Mouth
- p
-
Particle or droplet
- pk
-
Peak
- pr
-
Period
- r
-
Relative
- rm
-
Room
- s
-
Particle or droplet surface
- sl
-
Solid fraction
- st
-
Steady
- sv
-
Saliva
- t
-
Terminal
- V
-
Exhaled viruses
- v
-
Water vapor
-
Ambient condition
- 0
-
Initial condition
NOMENCLATURE
- b
-
Distance from centerline (m)
-
Droplet size PDF (1/m)
-
Drag coefficient
-
Viral concentration (1/m3; 1/m2)
-
Specific heat at constant pressure (J/kg K)
-
Drag force (N)
- D
-
Diameter (m)
-
Mass diffusivity (m2/s)
-
Viral deposition rate (1/m2 s)
-
Viral exposure (s/m3)
- g
-
Gravitational acceleration (m2/s)
- H
-
Height (m)
-
Latent heat of vaporization (J/kg)
-
Viral inhalation
-
Average infectious dose
-
Infectious dose PDF
- K
-
Evaporation rate (m2/s2)
- k
-
Turbulent kinetic energy (m2/s2)
- M
-
Molar mass (kg/kmol)
-
Mass flow rate (kg/s)
- N
-
Number
-
Nusselt number
- P
-
Momentum (Ns)
- p
-
Pressure (Pa)
- Q
-
Quanta
- R
-
Specific gas constant (J/kg K)
- r
-
Radius or radial distance (m)
-
Risk of infection
- S
-
Surface (m2)
- Sh
-
Sherwood number
- s
-
Position (m)
- T
-
Temperature (K)
- t
-
Time (s)
- u
-
Velocity (m/s)
-
Velocity fluctuation (m/s)
- V
-
Volume (m3)
-
Weight force (N)
- y
-
Molar fraction
- z
-
Vertical direction (m)
Greek symbols
- β
-
Entrainment coefficient
- γ
-
Cloud velocity scaling factor
- ε
-
Turbulent dissipation rate (m2/s3)
- ζ
-
Centerline dimensionless distance
-
Normal distributed random vector
- θ
-
Breath cloud direction (rad)
- λ
-
Thermal conductivity (W/mK)
- ξ
-
Radial dimensionless distance
- ρ
-
Density (kg/m3)
- τ
-
Characteristic time(s)
-
Volume fraction
- ω
-
Mass fraction
Subscripts
- a
-
Dry air
- ab
-
Airborne
- ac
-
Air change
- avg
-
Average
- c
-
Cell
- D
-
Exhaled droplets
- d
-
Drag
- ds
-
Dissipation
- e
-
Evaporation
- ep
-
Exposure
- ex
-
Exhalation
- gr
-
Ground
- hl
-
Half-life
- in
-
Inhalation
- lq
-
Liquid fraction
- m
-
Humid air mixture
- max
-
Maximum
- min
-
Minimum
- mt
-
Mouth
- p
-
Particle or droplet
- pk
-
Peak
- pr
-
Period
- r
-
Relative
- rm
-
Room
- s
-
Particle or droplet surface
- sl
-
Solid fraction
- st
-
Steady
- sv
-
Saliva
- t
-
Terminal
- V
-
Exhaled viruses
- v
-
Water vapor
-
Ambient condition
- 0
-
Initial condition