The dispersion of spherical droplets in source–sink flows and their relevance to the COVID-19 pandemic

In this paper, we investigate the dynamics of spherical droplets in the presence of a source–sink pair flow field. The dynamics of the droplets is governed by the Maxey–Riley equation with the Basset–Boussinesq history term neglected. We find that, in the absence of gravity, there are two distinct behaviors for the droplets: small droplets cannot go further than a specific distance, which we determine analytically, from the source before getting pulled into the sink. Larger droplets can travel further from the source before getting pulled into the sink by virtue of their larger inertia, and their maximum traveled distance is determined analytically. We investigate the effects of gravity, and we find that there are three distinct droplet behaviors categorized by their relative sizes: small, intermediate-sized, and large. Counterintuitively, we find that the droplets with a minimum horizontal range are neither small nor large, but of intermediate size. Furthermore, we show that in conditions of regular human respiration, these intermediate-sized droplets range in size from a few μm to a few hundred μm. The result that such droplets have a very short range could have important implications for the interpretation of existing data on droplet dispersion.


I. INTRODUCTION
The transport of inertial particles in fluid flows occurs in many problems arising in engineering and biology, such as the build-up of microplastics in the ocean 1 and respiratory virus transmission through tract droplets. [2][3][4] The Maxey-Riley equation 5 describes the motion of a finite-sized spherical particle in an ambient fluid flow. The equation is a representation of Newton's second law in which the forces acting on the particle include a Stokesian drag force, an added mass force, a gravity force, the force due to the undisturbed flow, and a Basset-Boussinesq history term. The equation takes the form of a second-order, implicit integro-differential equation with a singular kernel and with a forcing term that is singular at the starting time. 6 The equation has been applied to model the dynamics of aerosol comprising particles of various density ratios, 7 the feeding mechanism of jellyfish, 8,9 and the dynamics of inertial particles in vortical flows. [10][11][12] The equation has also been applied to droplet-laden flows with a phase change at sub-Kolmogorov scales. 13 The Basset-Boussinesq term accounts for the drag due to the production of vorticity as the particle is accelerated from rest. It is difficult to include this term numerically and is often omitted on the assumption that particles move in a quasistatic manner. 14 This assumption breaks down in bubbly and slurry flows, where the Basset-Boussinesq term accounts for a quarter of the forces on the particle 14 when the density ratio R = 2ρ ′ f /(ρ ′ f + 2ρ ′ p ) is greater than 2/3, where ρ ′ f is the fluid density and ρ ′ p is the particle density. Recent advances 15 have shown that the full Maxey-Riley equation can be represented as a forced, time-dependent Robin boundary condition of the 1D diffusion equation. Here, the authors found that a particle settling under gravity relaxes to its terminal velocity according to t −1/2 ; however, if the Basset-Boussinesq term is neglected, it relaxes exponentially quickly. 16 In this paper, we examine the transport of inertial particles in source-sink flows. 17 Such a flow could represent the trajectories of water droplets emitted from coughing, sneezing, [2][3][4] or breathing and in the presence of extraction, such as an air-conditioning unit or air current. 18 Our simplified mathematical model yields to analytic treatment in certain limits of large and small droplets. This enables us to provide important physical insight into this complex problem, but we remark that the effects such as drag non-linearity 19 and turbulent dispersion 20 are not taken into account. Since the dynamics of settling droplets is significantly affected by their size, it is important to understand the impact that the emitted droplet size has on the destination of such a droplet in a source-sink flow. In particular, since droplets are vectors for infectious diseases such as COVID- 19, it is imperative that we understand the droplet dynamics in such flows to mitigate the spread of the disease. This paper is organized as follows: in Sec. II, the mathematical model is presented and non-dimensionalized. The results are presented in Sec. III for small (Sec. III A) and intermediate-sized (Sec. III B) droplets in the absence of gravity. Gravitational effects are considered for small droplets in Sec. III C and for intermediatesized droplets in Sec. III D. In Sec. IV, we present applications of our results for human breathing without (Sec. IV B) and with (Sec. IV C) the inclusion of extraction. Finally, we discuss our findings in Sec. V.

II. MATHEMATICAL MODEL
Consider a source producing air of density ρ ′ air and viscosity ν ′ air , with volume flux of Q ′ 1 , containing spherical liquid droplets of density ρ ′ drop , which are emitted with a characteristic velocity U ′ . Let us represent the 3D velocity field u ′ source (x ′ ) at a position x ′ of the emitted air as a point source of strength Q ′ 1 , centered at the origin in the Cartesian coordinates, 17 We include an extraction unit as a point sink of strength Q ′ 2 located at a position x ′ 0 as follows: The resulting airflow is given by the linear superposition of these two flows, The natural timescale of the problem emerges as T ′ = |x ′ 0 |/U ′ . We non-dimensionalize (3) according to which gives the non-dimensionalized expression for the airflow velocity with where X(t) is the position of the droplet at time t, v(t) =Ẋ(t) is its velocity, the dot indicates the time derivative, and with a ′ being the droplet radius, g ′ being the acceleration due to gravity vector, Re = U ′ |x ′ 0 |/ν ′ air is the Reynolds number, and St is the particle Stokes number. Note here that the Faxén correction terms 5 have not been omitted: they are identically zero since Δu = 0.
The approximate ratio of the Basset history drag to Stokes drag is O(St 1/2 ), which, for the range of St we are interested in, is generally much less than one. In the remainder of this paper, we neglect the Basset history term since we anticipate that its magnitude is negligible compared to the Stokes drag term for the parameters of interest to us, and the resulting equations arė subject to the initial conditions v(0) = u(X(0), 0), where X(0) lie on a circle surrounding the origin of radius |X(0)|. In (5), taking the limit hence, we can ensure that the non-dimensional initial velocity has unit magnitude by requiring |X(0)| = √ Λ.

A. Computational considerations
The resulting equations (8) are a set of three coupled secondorder non-linear ordinary differential equations (ODEs) for the position vector X(t). The algebra involved in computing the material derivative in (8) is straightforward, but cumbersome, and it is omitted here. This set of equations does not admit analytical solutions, in general, and so it must be solved numerically.
We solved the equations by expressing them as a system of six first-order ODEs using the MATLAB ® ode15s algorithm, a variable-step, variable-order solver based on the numerical differentiation formulas. 21 This was performed on a laptop equipped with

Physics of Fluids
In (10), for small droplets (St ≪ R) emitted from the source, the balance is between the first term on the left-hand side and the righthand side so that the velocity rapidly adjusts to the background flow v ≈ u| X .
We are interested in whether droplets move away from or toward the sink. To this end, we look for trajectories for which If we take x 0 = [1, 0, 0], then the trajectory that emerges from the source and travels in the direction of the negative x-axis is the one that gets the greatest distance away from the sink. Hence, let us consider this inequality in the first component, and along the line y = 0, We are interested in where the flow field changes direction, since this indicates the maximum distance the droplets emitted at the source can travel before moving toward the sink. To this end, let us choose a point x = −λ along y = 0 and z = 0; then, this inequality tells us that This inequality can hold only if γ > 1. This makes sense, since the flow is directed toward the sink only if the sink is stronger than the source. Figure 1 shows the trajectories for small droplets (St ≪ R) in the presence of a source-sink pair: the source is located at the origin (green disk) and the sink is located at x = 1 along the x-axis (red disk). For γ = 1 [ Fig. 1(a)], we have equal strength and droplets can take large excursions from the source before returning to the sink. As γ increases, the trajectories emanating from the source occupy an increasingly compacted region [Figs. 1(b)-1(d)]. We can use this inequality above to define a region such that small droplets do not get further than a distance |λ| before traveling toward the sink. The circle with radius |λ| is shown in Fig. 1 (dashed curve). Observe that, as one gets increasingly close to the source (λ → 0), the inequality tends to meaning that in order to maintain trajectories moving away from a given test point, the sink strength needs to increase quadratically with the distance of the test point to the source.
B. Intermediate-sized droplets in the absence of gravity which represents a balance between inertia and drag forces. Provided γ > 1, and in the absence of gravity, in the long-term, the droplet will always migrate toward the sink. However, in the case of intermediate-sized droplets, the maximum distance traveled by the droplet before it moves toward the sink is given by |v(0)|/(R/St).
Since the initial velocity of the droplet is chosen to be the same as the surrounding fluid, then we can write the maximum distance as |u(X(0), 0)|/(R/St). As explained above [see (9)], in our nondimensionalization, our characteristic velocity U ′ was chosen to be that of the outlet. Hence, in this non-dimensionalization, |u(X(0), 0)| = 1. Figure 2 shows the trajectories of intermediate-sized droplets for γ = 5 in the absence of gravity. The striking feature of the plot is the shift from a regime where the maximal extent of the trajectories as predicted by (14) is no longer valid and must be replaced with a circle of radius St/R. In Fig. 2(a), St/R = 0.1 so that the droplets are slowed down rapidly before following the fluid flow. In Fig. 2(b), St/R = 1, meaning that the droplets are slowed down over the area covered by the unit circle, before being brought to the sink as ideal tracers. Finally, in Fig. 2(c), St/R = 2 so that the droplets travel a non-dimensional distance of 2 before being slowed down enough to be pulled into the sink.
Hence, we find that, in the absence of gravity, we can have two very different behaviors depending on whether we have small droplets St ≪ R or intermediate-sized droplets St ≥ R. Small droplets cannot get further than a distance ( √ γ + 1)/(γ − 1) from the source before traveling toward the sink, but intermediate-sized droplets are not restricted by this and can travel further than this, provided

C. The effect of gravity on small droplets
As the droplets move from the source to the sink, gravity attempts to pull them vertically downwards. Over the timescale of the problem, i.e., the average time it takes for a droplet to travel from the source to sink, gravity may or may not have an appreciable effect. Intuitively, one would imagine that smaller droplets are influenced more by the airflow than gravity: for stronger sinks, the effect of gravity is comparatively less. Intuitively, one would also expect that this holds true, provided that the source and sink are not too far away. The gravitational vector is non-dimensionalized according to The dashed circle indicates the predicted maximal distance that a droplet can travel in this regime, calculated using the inequality (14). The source is indicated by a green filled circle, and the sink is indicated by a red filled circle.

Physics of Fluids
U ′2 /|x ′ 0 | as shown in (7) so it depends on the initial speed and the distance between the source and sink.
For St ≪ R < 2/3, and in the absence of gravity, there are three fixed points: the source, the sink, and a saddle point located at x = −|λ| along the x-axis (Fig. 1). When gravitational effects are included, the fixed point at x = −|λ| moves clockwise about the origin as the effect of gravity is increased [see Fig. 3(a)]. A fourth fixed point (saddle) is created far from the source-sink pair, which gradually moves toward the sink [Figs. 3(b) and 3(c)] as the effect of gravity is increased. In Fig. 3(d), the separatrices (indicated as the red dashed-dotted curves) show that there is a wedge of trajectories that escape the pull of the sink. As might be expected, these trajectories are those that point directly away from the sink. Our results show that even for small droplets, gravity can be important if either the sink is far away or if the ejection speed is too low.

D. The effect of gravity on intermediate-sized droplets
Small droplets are deflected by gravity but generally feel the pull of the sink. Whether or not they are pulled in is determined by the interaction of gravity, the angle of their trajectory, and γ. As the droplets become larger, gravitational effects dominate and the sink becomes ineffective. In Fig. 4, we show how the droplet trajectories behave as St is increased. Figure 4(a) shows the familiar situation where the droplets are so small that gravity does not appreciably affect their trajectory over the characteristic lengthscale.
As St is increased, Fig. 4(b) shows that there are a range of trajectories with ejection angles α (defined with respect to the positive sense of the x-axis) around the source, which are deflected downwards away from the sink. This is consistent with Sec. III C. In these plots, R = 0.001, Λ = 0.0001, and |g| = 0. The dashed circle indicates the predicted maximal distance that a droplet can travel in this regime, calculated using the inequality (14). The dashed-dotted circle indicates the maximal distance predicted by the inertia-drag balance, giving radius equal to St/R. The source is indicated by a green filled circle and the sink is indicated by a red filled circle.

Physics of Fluids
However, at a critical St ≈ 2.5 × 10 −6 , each trajectory is deflected downwards by gravity [ Fig. 4(c)]. In this case, the maximum horizontal distance traveled by the droplets is very small. Interestingly, this trend is not monotonic. Further increasing St, the trajectories adopt a ballistic trajectory [ Fig. 4(d)]. Such droplets can move in very close proximity to the sink but are not pulled into it [ Fig. 4(d)].

IV. EXAMPLES OF APPLICATION
A. Background on respiratory virus transmission One of the possible applications of this paper is to underpin more sophisticated analytical or numerical models to study the transmission of respiratory viruses. In medical applications, it is common practice to categorize the emitted fluid particles as larger droplets from 5 μm to 1 mm in diameter, which have a ballistic trajectory, and aerosol that remains airborne. 23 Droplets smaller than 5 μm and the desiccated droplet nuclei are known as aerosol, which can remain airborne for several hours. [24][25][26] Respiratory viruses are transmitted from the virus-laden fluid particles to the recipient through (1) aerosol inhalation, (2) droplet deposition on the recipient's mouth, nose, or conjunctiva, or (3) droplet deposition on a surface and successive transmission through physical contact. 27 The SARS-CoV-2 virus, for example, has a diameter of 70 nm-90 nm, 28 and it is carried by droplets and aerosol. 26,29 The model proposed in this paper can provide new insights into the aerosol transmission, i.e., through those particles whose Stokes number is not sufficiently large to have a ballistic trajectory. The relative importance of aerosol (1) and droplet [ (2) and (3)] virus transmission is not always known, and it is yet to be established for SARS-CoV-2. 30 Counterintuitively, it has been argued that aerosol could be more dangerous than larger droplets. 31 Smaller droplets (≤5 μm) suspended in aerosol might carry a higher concentration of virus than larger droplets (>5 μm). 30,32,33 The largest droplets are less likely to penetrate deeply in the respiratory system and might be deactivated by the effective first structural and defense barrier of the mucosa. 34 Conversely, aerosolized virus half-life exceeds 1 h 26 and can be transported airborne through inhalation deep into the lungs, [35][36][37][38] avoiding the defense mechanisms of the upper respiratory system. Furthermore, aerosol inoculation has been shown to cause more severe symptoms than droplets administered by intranasal inoculation, and the dose of influenza required for inoculation by the aerosol route is 2-3 orders of magnitude lower than the dose required by intranasal inoculation. 2,33,38 To apply our model to aerosol dispersion, we consider the particles ejected by a person talking. A person ejects about tens of fluid particles per second with diameters between 39 0.1 μm to 1 mm and with a speed of the order 40 of 1 m s −1 . Because this is the most frequent source of aerosol, this accounts for most of the aerosol inhaled by other people. 41,42 Coughing leads to the ejection of 100-1000 fluid particles per second with a speed around 10 m s −1 , while sneezing generates 1000-10 000 fluid particles per second with a speed of up to 43 20 m s −1 . The values presented in this paragraph should be taken as indicative because there is a significant variability between different experimental studies. 2,33,36,38,[44][45][46][47][48][49][50][51][52][53][54][55] Some of the physics that is not considered in this work is the particle-particle interaction and evaporation. In fact, fluid particles are ejected through a jet that transports particles in the range of 2 μm-150 μm, 45,56,57 i.e., the aerosol, while the largest droplets have a ballistic trajectory independent of the surrounding flow. 2,57,58 The jet can be either laminar or turbulent when breathing and speaking, while coughing and sneezing always results in a turbulent jet with a diameter-based Reynolds number higher 2 than 10 4 . Once ejected, the air jet extends along a straight trajectory; its diameter increases linearly with the traveled distance, while the mean velocity linearly decreases, and the turbulent statistics remain constant (i.e., the jet is self similar 59 ). Once the largest particles with a ballistic trajectory have left the air jet, the jet bends upwards due to the buoyancy force caused by the temperature and thus density difference. 2 Smaller size particles (≤100 μm) are transported by the jet while they evaporate. Once a droplet exits the jet, it falls at its settling speed. For a particle with a diameter of 50 μm and 10 μm, the settling speed is less than 0.06 m s −1 and 0.03 m s −1 , respectively. The smallest of these two droplets is likely to land in the form of a desiccated nucleus. In fact, while a droplet with a diameter of 50 μm evaporates in about 6 s, a 10 μm droplet evaporates in less than 2,60 0.1 s, although their survivability also depends on the ambient temperature and relative humidity. 20 Once these droplets leave the jet, they can still be transported by ambient air currents, which have speeds typically in excess 61 of 0.01 m s −1 . These currents are modeled by the sink-source flow field discussed in this paper.

Physics of Fluids
A key issue that is discussed in this study is the extent to which the cloud of droplets and aerosol are displaced into the neighboring environment, as this is associated with virus transmission risk. Previous studies estimated that the overall horizontal range of the droplets generated while breathing and coughing before they land on the ground is around 1 m-2 m. [56][57][58] These studies led to the Centers for Disease Control and Prevention (CDC) 62 and World Health Organization (WHO) 63 social distancing guidelines. Nonetheless, the complex physics involved, which includes knowledge of the particle size distribution, their speed of evaporation, the viral charge of droplets of different size, the diffusivity of the virus-laden particles, etc., makes it difficult to assess what is the effective dispersion of the virus-laden fluid particles into the environment once ejected. It was found that the largest droplets generated by sneezing can reach a distance as far as 8 m, 2,3,57 while aerosol dispersion is highly dependent on the temperature, humidity, and air currents. For these reasons, this paper does not aim to provide definitive measures for the aerosol displacements but contributes to building a body of evidence around this complex question.

B. Predicted droplet dispersion
Currently, there is a large amount of disagreement in the reported spectra of droplet sizes in respiratory events. 2 The analysis is complicated by various factors including the evaporation of the droplets as they travel from the source, which, in turn, is influenced by ambient humidity and temperature. Recent mathematical modeling of droplet emission during talking have categorized droplets into one of the three groups: 64 small (<75 μm), intermediate (75 μm-400 μm), and large (>400 μm). Small droplets approximately follow the air and can travel a great distance by weakly feeling the effects of gravity. Large droplets can also travel a large distance due to their inertia. However, the intermediate-sized droplets feel strongly both gravity and drag, and their trajectory is a complex interaction of   66,67 b Taken from wind tunnel experiments. 55 these effects. Similar trends were observed in computational fluid dynamics simulations of previous authors. 65 In this section, we examine the problem from a much simplified perspective: we ignore evaporation entirely. We model the situation as a point source emitting droplets of various sizes in the presence of gravitational forces and compute the maximum horizontal distance traveled by these droplets. In this case, Q ′ 2 = 0 l min −1 , and other quantities such as jet speed, direction, and spread are taken from recent experimental studies of the authors: 66 these quantities are summarized in Table I.

Physics of Fluids
We find that for both heavy and quiet breathing, the maximum distance traveled by droplets L ′ (and the corresponding flight time τ ′ ) depends strongly on the droplet diameter (see Fig. 5). As expected, small droplets can travel many meters; however, we see that there is an intermediate range of droplet diameters where the horizontal distance is minimized. For quiet breathing, this minimum occurs between 69 μm < d < 77 μm, while for heavy breathing, this minimum occurs between 50 μm < d < 55 μm. This multi-modal behavior is reminiscent of that in previous experimental studies that measured the size distributions of droplets in various respiratory events such as talking and coughing 54,55 and sneezing. 43 The multimodal behavior observed in experiments is attributed to the different generation modes: bronchiolar, laryngeal, and oral. In our simplified model, we do not have any assumption on the biological origin of the droplet: the existence of the minimum is a characteristic of the droplets themselves and cannot be used as an indicator of the underlying droplet size distribution. The time it takes τ ′ decreases monotonically with increasing droplet diameter, as shown in Figs. 5(b) and 5(c).

ARTICLE scitation.org/journal/phf
In order to unpick the physics, observe that the drag force scales with the diameter of the droplet, but the weight of the droplet scales with the diameter cubed; hence, for large droplets, the drag force is negligible in comparison with the inertia of the droplet. As shown before, the droplets are slowed down exponentially in the horizontal direction and are accelerated in the vertical direction by gravity, giving the maximum horizontal range of the droplet (when nominally Y = −1), For large droplets (St ≫ R), we can then estimate that the maximum meaning that the trajectories are ballistic, and we expect that, for St ≫ R, the maximum distance becomes independent of St, in agreement with the observation that large droplets' trajectories are independent of the surrounding flow. 2,57,58 For small droplets St ≪ R, the drag decreases linearly with decreasing droplet diameter, but the weight rapidly decreases cubically with decreasing diameter. Hence, small droplets follow the airflow faithfully with little influence from gravity. Such droplets can get great distances before falling, as shown in the left-hand side of Fig. 5(a).
In the case of small droplets, the horizontal component of the droplet's trajectory follows the airflow like a tracer, and the droplet falls at its Stokesian settling velocity. Upon inspection, we find that the maximum horizontal distance L (when nominally Y = −1) tends to the following asymptote as St → 0: We can therefore estimate that droplets for which L > 1 or, equivalently, (i.e., the droplets travel farther in the horizontal direction than the vertical direction) weakly feel gravity.
In between these two extreme cases, the drag force on the droplet is the same order of magnitude as the gravitational force. By balancing these two effects, we can approximate the upper bound of St where the droplets become ballistic, Such droplets are not light enough to get carried any great distance by the ambient airflow but do not have large enough inertia to become ballistic. Hence, we have the following designations: (I) small droplets with St satisfying St < 2RΛ Reverting to dimensional quantities, we have the following range of intermediate-sized droplets:  Plugging in the numbers from Table I, we have the approximate range for quiet breathing and for heavy breathing. Our upper bound is in good agreement with previous categorizations of droplets, 64 although our lower bound seems to be smaller than those found by previous authors.

ARTICLE
scitation.org/journal/phf for breathing, and the typical speed of a jet in normal breathing conditions is of the order of U ′ = 1 m s −1 . In violent respiratory events, such as sneezing or coughing, these values could be significantly higher. 2 Finally, the extraction unit is located a distance of |x ′ 0 | = 0.2 m from the person. These quantities are summarized in Table II. Based on these numbers, the non-dimensional parameters that govern the trajectory of the droplet are determined to be R = 0.001 15, Re = 12, 225, Λ = 0.000 22, and |g| = 1.96, and the Stokes number ranges approximately from 10 −9 to 10 −1 . The parameter γ relates the flux of the extraction unit to the flux of a human's breath, and its effect will be examined. In particular, if we suppose that the envisaged extraction unit has a volume flux approximately equal to that of a standard vacuum cleaner (2832 l min −1 ), then we can approximate that γ ≈ 436. In Fig. 7, we show the efficacy of such extraction for a range of St. Extraction is very effective at low St; however, for St > 8.5 × 10 −5 , such extraction is ineffective. This upper bound of the Stokes number corresponds to water droplets of diameter 71 μm. Droplets larger than this will not be collected by extraction. In the droplet classification of Sec. IV B, the effective range of extraction corresponds to nonballistic droplets.

V. DISCUSSION AND CONCLUSION
In this paper, we have presented a simplified mathematical model for droplet dispersion from a source and in the presence an aerosol extractor. In the absence of gravity, and for St ≪ R, droplets behave as ideal tracers, and the maximum distance that they can travel before being extracted is a function of γ only. In this case, there are two (source, and sink if γ = 1) or three (source, sink, and saddle if γ > 1) fixed points. The fixed points in this study are colinear, and the position of the saddle depends on γ alone, for any given distance between the source and sink. For moderate St, the droplets' inertia carry them far away from the source until they are slowed down by drag forces and pulled into the sink. In this case, the maximum distance that the droplets can travel is given by R/St.
When gravity effects are taken into account, the saddle point for St ≪ R is no longer colinear but moves on an arc, clockwise about the source, and a fourth fixed point (saddle) emerges approximately below the sink fixed point. For fixed γ, this fixed point moves closer to the source as the magnitude of gravity is increased. In this case, there is a set of trajectories that are pulled away from the sink by gravity. For moderate St, gravity plays an increasingly important role, and there is a critical value of gravity that pulls all trajectories vertically downwards away from the source. For yet larger St, the trajectories adopt a ballistic trajectory, with even those that travel close to the sink not being pulled in.
COVID-19 has brought increased awareness of the risks of aerosol generating procedures (AGPs) across all fields of medicine, highlighting the need for a deeper understanding of droplet dispersion and categorization during respiration and AGPs. Clinicians recognize that our historical approaches to protection during AGPs are no longer adequate and that many additional precautions are necessary. In order to develop the most effective solutions, a critical first step is understanding the behavior of droplets generated during AGPs. This paper allows us to predict this behavior and inform our understanding of "at risk" zones in the vicinity of an AGP. In particular, we performed simulations relevant to human respiration, as well as simulations to inform the development of an aerosol extractor for use in clinical settings. These simulations can help to guide recommendations on maximum safe distances between the source and sink.
Additionally, these models provide a better understanding of the behavior of individual droplets of various sizes, which may be present in a wide range of aerosols contaminated with viruses or other pathogens. This may help clinicians to make better informed decisions regarding safety while performing AGPs and in their choices of the type of PPE they wear. Finally, these models provide a basis on which aerosol and droplet contamination from a wide range of surgical, medical, dental, and veterinary AGPs can be modeled while taking into account airflows in confined clinical spaces. In this case, we found that, for St ≤ 8.5 × 10 −5 , all of the aerosol is extracted and that gravity has a minimal effect; this St corresponds to droplets with approximate diameter equal to 0.07 mm. Droplets larger than this are affected by gravity, and for St = 10 −2 , corresponding to droplets equal to 0.78 mm, none of the droplets are extracted. Such large droplets would be typically captured by personal protective equipment (PPE), such as FFP1 masks that have pore sizes typically smaller than 1 μm.
We determined the maximum range of droplets ejected from the source in the absence of a sink and found that the range is minimized for intermediate-sized droplets. We find that, in human respiration, this pertains to droplets within the observed range of the ejected droplets. This could have implications for the interpretation for data obained from experiments on biological subjects, in particular, those that attribute observed bi-and tri-modal droplet dispersion to biological functions. Our studies suggest that the bi-modal nature of the curve is a function of the droplet's Stokes number and not necessarily linked to a specific biological function.
In our model, we neglected the Basset history term in the Maxey-Riley equation. The Basset history term is of significant importance for bubbly flows, where it can account for a quarter of the instantaneous force on a bubble. 14 Generally speaking, for R ≪ 2/3, this term can be safely ignored for small and intermediatesized droplets. Recent studies have also shown that neglecting it in the modeling of raindrop growth leads to a substantial overestimate of the growth rate of the droplet. Hence, for the solutions that become ballistic, we expect that such trajectories would be influenced by the Basset history term that should be included. To do this efficiently, there is a very promising method developed recently. 15 Since this is not the focus of our study (such droplets can be captured by other forms of PPE), we do not perform such a study here.
If the aerosol route of transmission is confirmed to be important by the World Health Organization, 20,69 we will need to reconsider guidelines on social distancing, ventilation systems, and shared spaces. To ensure that we put in place the correct mitigating measures, for example, face coverings, we need to have a better understanding of the different droplet behaviors and their different dispersion mechanisms depending on their size. This paper contributes to this debate by providing a new framework for categorizing droplets depending on their dispersion mechanism.