J. Fluid Mech. (2014), vol. 745, pp. 537–563.
doi:10.1017/jfm.2014.88

c Cambridge University Press 2014

537

Violent expiratory events: on coughing andsneezing
Lydia Bourouiba1,2, †, Eline Dehandschoewercker3 and John W. M. Bush1
1 Department of Mathematics, Massachusetts Institute of Technology, Cambridge, MA 02130, USA
2 Department of Civil and Environmental Engineering, Massachusetts Institute of Technology, Cambridge,

MA 02130, USA

3 PMMH - ESPCI, O207 10, rue Vauquelin, 75005 Paris, France(Received 4 January 2013; revised 29 December 2013; accepted 7 February 2014)

Violent respiratory events such as coughs and sneezes play a key role in transferring
respiratory diseases between infectious and susceptible individuals. We present theresults of a combined experimental and theoretical investigation of the fluid dynamics
of such violent expiratory events. Direct observation of sneezing and coughing events
reveals that such flows are multiphase turbulent buoyant clouds with suspended
droplets of various sizes. Our observations guide the development of an accompanying
theoretical model of pathogen-bearing droplets interacting with a turbulent buoyant
momentum puff. We develop in turn discrete and continuous models of droplet fallout
from the cloud in order to predict the range of pathogens. According to the discrete
fallout model droplets remain suspended in the cloud until their settling speed matches
that of the decelerating cloud. A continuous fallout model is developed by adapting
models of sedimentation from turbulent fluids. The predictions of our theoreticalmodels are tested against data gathered from a series of analogue experiments inwhich a particle-laden cloud is ejected into a relatively dense ambient. Our study
highlights the importance of the multiphase nature of respiratory clouds, specificallythe suspension of the smallest drops by circulation within the cloud, in extending therange of respiratory pathogens.
Key words: biological fluid dynamics, multiphase and particle-laden flows, turbulent flows

1. Introduction

Emerging or reemerging viral and bacterial infectious diseases (e.g. H1N1, H5N1,
SARS, tuberculosis) have increasingly high human and economic consequences(IMF/World Bank 2006). Confined environments, such as airplanes, hospitals andschools, serve as mixers where pathogens can stay suspended and spread fromhost to host. Hence, understanding the dynamics of pathogens indoors is criticalto improving the modelling and control of epidemics (Settles 2006; Tang et al.
2006; Weber & Stilianakis 2008). Nevertheless, the transmission mechanisms of eventhe most common respiratory diseases remain poorly understood. Three modes oftransmission are discussed in the medical literature. Self-inoculation may arise through† Email address for correspondence: lbouro@mit.edu

538

L. Bourouiba, E. Dehandschoewercker and J. W. M. Bushdirect contact with the mucus (or other bodily fluids) of an infectious subject. Large
droplet transmission may arise through the spraying of infected droplets directly onto
the conjunctiva or mucus of a susceptible host via coughing or sneezing. Airborne
transmission may arise through inhalation of relatively small infected droplets or thepathogen-bearing solid residues of size <5–10 µm, referred to as droplet nuclei, thatcan form from the small droplets via evaporation (e.g. Nicas, Nazaroff, & Hubbard
2005; Tellier 2006). The first two modes are termed direct short-range routes ofpathogen transmission, both requiring the close proximity of individuals, while thethird is an indirect long-range route of transmission.
The literature remains undecided as to the relative importance of long-range airborne
transmission and short-range large droplet transmission for many diseases, including
SARS, avian influenza and whooping cough. For example, Tellier (2006, 2009)
and Clark & de Calcina-Goff (2009) vigorously argue the importance of airborne
transmission for avian influenza viruses. For SARS, a posteriori studies revealed
that the transmission involved a combination of two modes, short-range large droplet
transmission via direct contact and airborne transmission via transport of droplet
nuclei in the ventilation system of healthcare facilities (e.g. Wong et al. 2004; Yu
et al. 2004). The basic mechanisms of direct and indirect pathogen transfer between
people remain poorly understood (e.g. Tang et al. 2006; Weber & Stilianakis 2008;
Bourouiba & Bush 2012). Experimental studies have been aimed at measuring the size
distribution of droplets emitted by various respiratory functions, such as coughing andsneezing (e.g. Wells 1934, 1955; Duguid 1945; Loudon & Roberts 1967; Papineni &Rosenthal 1997; Yang et al. 2007; Morawska et al. 2009b; Johnson et al. 2011). The
discrepancy between reported data on droplet size spectra poses a basic impediment
to the mathematical modelling of respiratory disease transmission. Fluid dynamics
studies of the cloud resulting from violent expirations have been presented by Settles
and coworkers (e.g. Settles 2006; Tang, Path & Settles 2008; Tang et al. 2009),
who combined schlieren optical method and particle image velocimetry in order toreveal the velocity field within the initial phase of a cough cloud (Jonassen, Settles &Tronosky 2006). Although the dynamics of the gas phase were detailed, the droplet
dynamics were not given particular attention.
Violent expirations release multiphase turbulent flows that are generally composedof buoyant hot moist air and suspended droplets of various sizes. These droplets
contain components such as pathogens and minerals that can form droplet nuclei after
evaporation of the drop’s liquid phase. At low Reynolds number, the Stokes settling
speed US = gd2 /(18µ)(ρd − ρ) of a droplet of diameter d and density ρd in an ambient
gas phase of density ρ is proportional to its surface area which necessarily decreases
with time due to evaporation. It remains a common belief that most contamination of
new individuals occurs via direct contact with large droplets emitted by an infectious
individual. Wells (1934, 1955) challenged this view and was the first to examine
the role of airborne transmission of respiratory disease transmission. He compared
the time for complete evaporation to the settling time of various droplets ranging in
diameter from 1 to 1000 µm. He reported that drops with diameter d > 100 µm settle
to the ground in less than 1 s, without significant evaporation (Wells 1934). Whileevaporation rates will in general depend on ambient temperature and humidity, drops
with d < 100 µm will typically become droplet nuclei before settling (Wells 1955).
Droplets with d < 5–10 µm rapidly evolve into droplet nuclei with settling speeds
less than 3 mm s−1 , and so may be readily suspended and advected by the cloud of
air emitted by the cougher, or resuspended by any ambient flow as may arise, forexample, through air conditioning. Thus, suspended droplet nuclei are expected to be
critical elements in long-range airborne transmission.

Coughing and sneezing539

The size distribution of respiratory droplets has been the subject of an increasing
number of studies (e.g. Duguid 1946), with recent focus on the improvement of the
precision of measurement of the small submicrometre range (e.g. Morawska 2006;
Johnson et al. 2011; Zayas et al. 2012). An example of droplet size distribution forcoughs is shown in figure 17, which indicates a peak drop size of ∼15 µm, theassociated settling speed being 6.5 mm s−1 in ambient air typical of winter indoor
conditions. In addition to having potential for long-range transmission, it is generally
thought that smaller drops and particles better penetrate the respiratory tract, reaching
deeper target tissues within the lungs. Indeed, experiments of inoculation of animals
with respiratory pathogens and solid particles (e.g. Sonkin 1951; Wells 1955) suggest
that, given the same viral load, inhalation of the atomized solution results in higher
infection and death rates than direct intranasal inoculation.
2. Observations

Puffs, jets, thermals and plumes are turbulent bodies of fluid ejected from a
localized source either continuously (in the case of plume and jets) or discontinuously(in the case of puffs and thermals). Thermals and plumes are buoyancy driven, whilepuffs and jets are momentum driven (Scorer 1978; Turner 1979). As they proceedfrom the source, these flows retain roughly the same shape and their lateral extent
increases linearly with distance from the source. This observation of self-similarity
led to the entrainment hypothesis (Morton, Taylor & Turner 1956) which allows one
to render such turbulent flows analytically tractable by assuming that the details of thevelocity fluctuations can be averaged out. According to the entrainment hypothesis,
the speed w of the ambient fluid drawn into the cloud is proportional to the average
cloud speed U, so the cloud radius r increases linearly with the distance s fromthe source: w = αU and r = αs (Morton et al. 1956) (figure 1a). The values ofthe entrainment coefficient α are observed to be roughly constant although withconsiderable variability for each type of flow, be it thermal or plume, puff or jet.
This turbulent plume theory has proven to be sufficient to capture key features ofmany geophysical and environmental flows, including the emission of industrial
effluent, the emission of volcanic hydrothermal plumes or thermals and the ejectionsof bubble plumes in lakes (e.g. Scorer 1978; Cardoso & Woods 1993; Gonnermann& Manga 2007; Woods 2010). Here, we adapt such an approach in order to describe
the dynamics of violent expiratory events. The fate of the suspended droplets isdetermined by drawing an analogy with multiphase particle-laden clouds as arises
in a variety of geophysical systems such as riverine outflows or industrial waste
discharges (e.g. Scorer 1978; Socolofsky, Crounse, & Adams 2002; Bush, Thurber, &Blanchette 2003; Hunt et al. 2007).
We here present the results of a combined experimental and theoretical study
aimed at elucidating the dynamics of pathogen-bearing droplet dispersion by violent
expiratory events. In § 2, we present the results of our visualization of real humancoughs and sneezes. An analogue experimental study of multiphase puffs is presented
in § 3, and allows us to refine our physical picture of expiratory clouds. The resulting
theoretical model is developed in § 4. The application of our theoretical model toclinical data is presented in § 5.
We began by visualizing the dynamics of coughs and sneezes produced by adults.
A schematic illustration of the experimental set-up used in our study is presented infigure 2. The head of the subject was positioned in front of a black felt backdrop inorder to create a black-field effect. The area was surrounded with lights oriented to540

L. Bourouiba, E. Dehandschoewercker and J. W. M. Bush
(a)

(b)
w

s
U

r

U

g

y
s
x

r

F IGURE 1. (a) Schematic of the spreading and settling of the particles from a cloud
emitted horizontally. Note the self-similarity of the cloud as it grows by entrainment of thesurrounding fluid. (b) A schematic illustration of the trajectory of the multiphase buoyant
puff of radius r(t) at position s(t) from the source. Gravity deflects the cloud motion
from the horizontal, causing it to rise. Note that in this example, the puff is initially shot
downward at an angle θ0 6= 0.
yg
LightsxBlack felt

High speed camera

F IGURE 2. Experimental set-up for the visualization of the violent respiratory events using
high-speed cameras.

maximize light scattering by the emitted droplets. High-speed imaging of various
violent expirations was recorded at 1000–4000 frames per second (f.p.s.) witha Phantom high-speed video. A monochrome camera was used to reduce lightsensitivity. Note that despite the high density of ejecta apparent in figures 3 and 4,
considerably more droplets and mist were observed with the naked eye than with thecamera. Nevertheless, our visualizations were sufficient for our purposes. In certain
experiments, a smoke generator enabled tracking of the emitted gas phase. Coughs
and sneezes were found to consist of a turbulent cloud of buoyant gas with suspended
droplets. The largest droplets follow a ballistic trajectory relatively unaffected by theflow in the gas phase, while the smaller drops were suspended to varying degrees
within the turbulent gas cloud, thereby having their range extended.

541

Coughing and sneezing
(a)

(b)

(c)

(d )
70 cm

(f)

70 cm

Smoke
generator

50 cm

(e)

F IGURE 3. High-speed images of a cough recorded at 1000 frames per second (f.p.s.)
reveal the dynamics of the expelled gas and liquid phases. The sequence is displayed
for the times (a) 0.006 s, (b) 0.01 s, (c) 0.029 s and (d) 0.106 s. (e) Large droplets
are ejected and their trajectories shown in this streak image. (f ) A typical cough airflow
is visualized using a smoke generator and recorded at 2000 f.p.s. Note that (e) and (f )are the superposition of the instantaneous images of the droplets and cloud trajectories
throughout the sequence.

Figure 3(a–d) shows the visualization of a cough, recorded at 1000 f.p.s. Thissequence shows the time evolution of the emitted spray until 106 ms after the onsetof visible droplet emissions. The early images (figure 3a–c) indicate the cone-likeshape of the ejecta cloud and the high density of droplets emitted by the cougher.
The largest droplets are tracked visually, resulting in the streak image shown infigure 3(e). The shape of the trajectories suggests a ballistic path for the larger
droplets whose dynamics is not greatly affected by the cough cloud. We also usedsmoke to enhance our visualization of the cough cloud. Figure 3(e,f ) allow us tocompare the trajectories of the emitted droplets with the path of the air cloud, as
visualized with smoke. The air flow emitted was recorded at 2000 f.p.s. The smoke
cloud was initially directed downward at an angle 24 ± 7◦ from the horizontal andassociated with an entrainment coefficient of 0.24 ± 0.02 in the initial jet phaseclose to the mouth and 0.132 ± 0.06 in the second puff phase. The airflow was
expelled for approximately 300 ms as highlighted by the use of smoke in figure 3(f ),where we superpose the video sequence from the onset of the cough to the end of

542

L. Bourouiba, E. Dehandschoewercker and J. W. M. Bush

(a)

(b)

(c)
(e)

(d )
(f)

70 cm
70 cm

(g)
(h)

g
y

V0

c

x
B0

I0

F IGURE 4. (Colour online) High-speed images recorded at 1000 f.p.s. of a sneeze attimes (a) 0.006 s, (b) 0.029 s, (c) 0.106 s, (d) 0.161 s, (e) 0.222 s and (f ) 0.341 s.
(g) Trajectories of the largest droplets are revealed through a streak image. (h) A
schematic of the initial cloud with characteristic initial momentum I0 , buoyancy B0 andvolume V0 emitted into quiescent ambient air.

the recording. The duration of coughs was estimated to be 200–500 ms in Gupta,
Lin & Chen (2009), with peak flows arising at 57–110 ms for female subjects and
57–96 ms for males and opening angle close to the mouth leading to approximate
entrainment coefficient of 0.21. Our observations are thus consistent with the ranges
previously reported. The expelled volume of the female subject in figure 3(f ) was
estimated to be 0.6 ± 0.1 l by recording the volumes of water displaced by a series

543

Coughing and sneezing
(a)

35 cm

(b)

35 cm

F IGURE 5. (Colour online) (a) The effect of buoyancy on the sneeze cloud is apparent
in its upward curvature. (b) The large droplet trajectories are shown in the streak image
recorded at 2000 f.p.s.

of coughs directed into a bath via a tube. By way of comparison, measurements onlarge groups of patients reported that the volume expired during coughs was in therange 0.25–1.25 l for females and 0.4–1.6 l for males (McCool 2006; Gupta et al.
2009; Park et al. 2010). The mouth opening area in this figure was estimated to be
3.4 cm2 during the cough, which is also consistent with values reported in previous
studies, 4 ± 0.95 cm2 for males and 3.37 ± 1.4 cm2 for females (Gupta et al. 2009).
Figure 4 shows analogous visualizations of a sneeze recorded at 1000 f.p.s. The top
series of images shows the evolution of the cloud from 0.006 to 0.34 s from the onsetof the sneeze. The cloud is evidently more dense than that generated by the cough(figure 3); moreover, the range of visible droplets is larger. Figure 4(g) shows the
trajectories of the largest droplets emitted in the sneeze, whose dynamics are markedly
different from those of the smaller drops transported within the cloud. Sneezes werealso observed to be shorter than coughs for the same individual, the observed duration
being approximately 200 ms. Our visualizations capture the dynamics of the cloud
up to approximately 1 m from the source. Using the observed ejection speed andgeometry, we estimate a typical Reynolds number of 104 for coughing and 4 × 104
for sneezing. The estimated entrainment coefficients for the sneeze shown in figure 5
are 0.13 ± 0.02 in the initial jet phase and 0.055 ± 0.01 in the subsequent puff phase.
Figure 5 shows the comparison between the trajectory of the exhaled sneeze cloud
and the larger droplets recorded at 2000 f.p.s. The cloud trajectory was particularly
helpful to visualize the effect of buoyancy. The values of momentum and buoyancyof the cloud vary with distance from the source as ambient air is entrained and thedroplets settle out of the cloud. We note that the apparent kink in the trajectoryof the cloud, a common feature in the sneezes recorded, occurs at a time where
the average momentum and buoyancy of the cloud become comparable, as will be
discussed further in section § 4.
Figure 6 shows an image sequence of a sneeze recorded at 2000 f.p.s. The
circulation within the cloud is readily apparent on the video and highlighted inthis figure with a pair of arrows. Note the change of orientation of the arrows, which
again illustrates the influence of buoyancy in deflecting the cloud upward. One canclearly see the trace of the mist droplets that remain suspended and recirculated within
the sneeze cloud even beyond 70 cm from the source. The residual mist suspended
in the cloud remains apparent until the end of the video. We also observe the settling

544

L. Bourouiba, E. Dehandschoewercker and J. W. M. Bush

20 cm

0.008 s0.034 s0.06 s0.13 sF IGURE 6. An image sequence of a sneeze cloud indicating the interaction between thegas and fluid phases. The arrows illustrate the apparent circulation within the cloud. Thechange in direction of the cloud is also shown by the shift in the direction of the arrows.
The total release time of this sneeze was 200 ms. Time since initial release is indicatedin each frame. The smallest mist droplets remain suspended until the end of the video,
over 1.35 s.
y

g
w(s)
Real
sources

Drop-laden
cloud

r
U(s)

o
Virtualorigin

x
xfall

xmax

Fallout
from cloud

Deposition:
Range of contamination

F IGURE 7. Physical picture of the evolution of a cough or sneeze cloud emitted
horizontally. Negatively buoyant suspended droplets are circulated within the cloud until
they settle out (at a horizontal distance xfall ) and reach a final position xmax , corresponding
to the range of contamination. The virtual origin is shown behind the sneezer source andis computed by extrapolating the measured trajectory of the cloud. The cloud climbs in
response to its increasing buoyancy as the negatively buoyant droplets fall out.

of small, long-lived droplets from the cloud at speeds as low as 0.06 m s−1 . Thissettling speed and time to complete evaporation are consistent with pure water droplets
smaller than 50 µm diameter that have a settling speed of less than 0.068 m s−1 andtime to complete evaporation of less than 6 s (Holterman 2003). Smaller settling
speeds were also observed, but difficult to estimate owing to their irregular motion.
Figure 7 illustrates the physical picture that emerges from a synthesis of our flow
visualizations, prior work on cloud visualization using schlieren (e.g. Tang et al.
2009) and droplet sampling (e.g. Morawska et al. 2009b). Violent respiratory events
release droplet-bearing turbulent momentum puffs of moist, buoyant air ejected in a
direction that depends on the individual. As the buoyant momentum puff advances, it
entrains ambient air, its width increasing and its speed decreasing with distance from

Coughing and sneezing545

the source (Turner 1979) as illustrated in figure 1(a,b). While large droplets follow
a ballistic trajectory largely unaffected by the turbulent cloud, smaller droplets may
remain suspended in the cloud and are circulated therein until they settle out. In § 4.2
we first describe the fate of suspended droplets using a discrete fallout model, where
drops remain suspended in the cloud until their settling speeds match the speed of
the cloud. A more sophisticated fallout description is then developed that accounts
for a continuous fallout of the droplets along the entirety of the cloud’s trajectory.
Such dynamics are illustrated for horizontally directed flows in figure 7.
3. Analogue experiments

As an experimental analogue to droplet-laden coughs and sneezes, we examined the
dynamics of turbulent, dyed, particle-laden fluid injected into a quiescent water tankfrom a piston (figure 8). The cloud buoyancy resulted from the use of salty waterin the tank and fresh water in the payload (see table 1). The water tank of height
60 cm, length 3.25 m and width 85 cm was divided into two parts, a release sectionof approximately 75 cm in length that contained the source piston, and an illuminatedobservation section of approximate length 2.5 m. A wall was used to separate thetwo regions in order to minimize disturbances to the observation area caused by
the piston release. The piston was placed 23 cm above the base of the tank, which
minimized wall and surface effects on the puff trajectory, such as attraction to theboundary via the Coanda effect. The piston, with internal radius 2.6 cm and opening
diameter of 0.9 cm, was filled and its payload (88 cm3 ) released over a characteristic
time of 375 ± 50 ms. This was fine-tuned to minimize the release time while still
maintaining a puff subsequent cloud dynamics. During release, the Reynolds numberof the fluid payload was approximately 5000. The trajectory of the resulting turbulent
puff cloud and the deposition pattern of its particles were recorded with a Nikon SLR
D90 camera at 24 f.p.s. A backdrop of black felt was used to minimize reflections
from the tank wall.
Three dozen experiments were carried out, five of which are detailed here (I–V see
table 2). One experiment (I) involved a buoyant puff cloud free of beads, while theother experiments (II–V) involved bead-laden buoyant clouds. The settling speed ofthe particles was measured in fresh water. The difference in density between the fresh
dyed water contained in the piston and the clear water of the tank ranged from 5.4 to
10 mg cm3 . The tank fluid density was measured using a densitometer (DMA 35N,
Anton Paar) at several different locations in the tank in order to ensure that it was
well-mixed. Care was taken to ensure that residual disturbances caused by the filling
process dissipated before release. Note that data acquisition ceased when the buoyant
puff approached the bottom or the top of the water tank. In addition, some of theclouds grew asymmetrically; which can be attributed to residual motion originatingfrom the release mechanism, for example, the formation of an asymmetric vortex ringat the front of the cloud (e.g. Scorer 1957). Such experiments were discarded fromconsideration.
As the cloud progresses, the ambient fluid is entrained and mixes with the cloud,leading to an increase in size and hence a deceleration of the puff (Turner 1979).
The entrainment of ambient fluid modifies the total mass and volume, but the cloudroughly retains its shape. The self-similar growth of the turbulent puff leads to r = αs,
where r is the half-width of the cloud and s the distance from the source (seefigure 1a). In reality, details of source conditions displace the effective origin of thecloud relative to the piston (Hunt & Kaye 2001). The virtual origin can be computedby extrapolating the measured r against s plot, as illustrated in figure 7.

546

L. Bourouiba, E. Dehandschoewercker and J. W. M. BushPartical-laden payload
2.6 cm
0.9 cm
7.5 cmLighting

CameraPistonObservation area

60 cm

1.5 m

23 cm

3.25 cm

85 cm

CameraF IGURE 8. A schematic illustration of our analogue experiment, conducted in a tankcontaining salty water. Buoyant, dyed, particle-laden payloads were released using the
rapid triggering of the piston. The release and observation regions were separated by awall to minimize disturbances to the particle-laden puff by the piston release.

Nd (mm)
ρp (g cm−3 )ums (cm s−1 )uts (cm s−1 )Repφρtank (g cm−3 )ρf (g cm−3 )v0 (m s−1 )B0 (kg m s−2 )I0 (kg m s−1 )

I

II

III

IV

V

0
—
—
—
—
—
0
1.00290.999 75
1.007
2.72 × 10−3
8.86 × 10−2

252
2.86 ± 0.04
2.4
35.5 ± 0.5
32.3
9343.51 × 10−2
1.00500.999 55
1.007
−3.77 × 10−2
9.29 × 10−2935
1.36 ± 0.04
2.4
22 ± 0.5
19.3
270
1.4 × 10−2
1.00500.999 55
1.174
−1.22 × 10−2
1.05 × 10−14200
0.78 ± 0.08
2.5
11.3 ± 0.5
11.8
921.19 × 10−2
1.00530.99951.118−1.04 × 10−2
1 × 10−1

1258
0.46 ± 0.04
2.5
6.6 ± 0.25
6.65
30
7.3 × 10−4
1.00960.999 55
1.102
7.73 × 10−3
9.70 × 10−2

TABLE 1. Parameters used in experiments I to V: number of beads N, bead diameter d,bead density ρp , measured settling speed ums , theoretical settling speed uts , and associatedparticle Reynolds number Rep = uts d/ν, where ν is a water kinematic viscosity. Here, φ isthe initial volume fraction of particles in the payload of volume 88 cm3 . We also report
the density of the tank water ρtank and fluid released ρf , as well as initial cloud velocity
v0 , momentum I0 and buoyancy B0 . The dynamic viscosity of water is µ = ρν = 1.002 ×
10−3 kg m−1 s−1 .

While these clouds are buoyant, they are initially driven by momentum, as is thecase for coughs and sneezes. Two phases of the cloud evolution were apparent in all
experiments. The first phase is dominated by a jet-like dynamics, corresponding to0.166 ± 0.01
0.10 ± 0.01
0.13 ± 0.022
0.10 ± 0.007
0.134 ± 0.003

0.028 ± 0.006
0.034 ± 0.003
0.018 ± 0.003
0.02 ± 0.003
0.021 ± 0.0029

α20.089 < 0.097 < 0.105
0.063 < 0.07 < 0.075
0.061 < 0.074 < 0.0856
0.057 < 0.061 < 0.066
0.075 < 0.077 < 0.08

αavg
25
35
35
25
19

xT (cm)—22 ± 237 ± 657 ± 463 ± 2

t
xfall
(cm)

—
29 ± 341 ± 660 ± 465 ± 2

t
xmax
(cm)

—
31
37
63
66

o
xmax
(cm)

TABLE 2. Entrainment coefficients measured in the two phases of the cloud dynamics in each experiment: in the jet phase α1 , in the puff
phase α2 , and their average αavg . Here, xT is the observed distance of transition from the jet to puff phases extracted from the cloud trajectorytt
and size data (e.g. figure 9). Here, xfall
is the predicted fallout distance of particles from the cloud. Here, xmax
is the predicted maximumodistance of deposition of particles (or range). Finally, xmax is the observed maximum distance of deposition of the particles at the bottomof the tank. All theoretical predictions are computed using αavg .

I
II
III
IV
V

α1

Coughing and sneezing
547

548

L. Bourouiba, E. Dehandschoewercker and J. W. M. Bush(a) 120
Experiment II
100s (cm)

80
60
40
20
0

5

10

15

20

25

30

35

t (s)
(b) 8
7

Experiment II

6r (cm)

5
4
3
2
1
0

10

20

30

40

50

60

70

80

90

100

s (cm)F IGURE 9. (Colour online) (a) Relationship between the distance travelled by the cloud’s
geometric centre s and the time t for experiment II. (b) Relationship between distances travelled and the mean radius of the puff cloud r for experiment II. Note that theobserved transition in slope is common to all experiments I–V. The slopes are used todetermine entrainment coefficients in the first jet phase α1 and in the second puff phase
α2 , as reported in table 2. Characteristic error bars are shown at different stages.

the high-speed release of the payload. The second phase is dominated by a puff-likedynamics, which is characterized by the self-similar growth of the puff cloud. Fora jet, the conservation of momentum flux M0 ∼ ρr2 ( ds/ dt)2 combined with the selfsimilar growth of the cloud r ∼ αs leads to ds/ dt ∼ 1/s; hence, s ∼ t1/2 . For a puff,
the conservation of momentum I0 ∼ ρr3 ds/ dt combined with the self-similar growthof the cloud r ∼ αs leads to ds/ dt ∼ 1/s3 ; hence, s ∼ t1/4 . Distinct jet and puff phases
were apparent in all our experiments.
Figure 9 illustrates the dependence of the position s on time t and the corresponding
dependence of r on s. The change of slope in the relationship between r and s
indicates the transition between the jet and puff regimes and the associated change inentrainment coefficient α. Similar changes of slope and scalings were observed for all
experiments. The entrainment coefficient of each phase was extracted and reported as
α1 and α2 in table 2. Here, xT , the distance of transition from the jet to puff phases,
is also given. The range of coefficients measured is 0.09 < α1 < 0.18 for the jet phaseand 0.015 < α2 < 0.037 for the puff phase (see table 2). Reported values of α forthermals ranged from 0.2 to 0.34 with a mean of 0.25 (Scorer 1957) or from 0.13 to
0.53 (Richards 1961). Richards (1965) showed that non-buoyant puffs and thermals

549

Coughing and sneezing
(a)

(b)

(c)

(d )

(e)

(f)

F IGURE 10. Illustration of the tracking of the front and geometric centre of the turbulent
cloud in experiment I as it moves away from the source. The successive images are taken
at times (a) 0.08 s, (b) 0.3 s, (c) 0.46 s, (d) 1.54 s, (e) 2.25 s and (f ) 3.5 s after release.
Scale bar, 10 cm.

are similar, with α ranging from 0.15 to 0.58. Richards (1965) noted that there is alarger variability in the entrainment coefficient of puffs than thermals.
Figure 11 shows the trajectory of the buoyant puff cloud released without beads(experiment I). The inset shows the superposition of images tracing the shape of thecloud as it grows away from the source. The front edge and shape of the cloud weretracked on the video frames. As was observed in both our experiments and previous
studies of thermals, the cloud was best described as an ellipsoid (Scorer 1957). The
tracking of the position of the geometric centre of the cloud (figure 10) was used toconstruct the time evolution of its position s. Note that in the inset and snapshots offigures 10–11, the upturn in the cloud is evidence of the influence of buoyancy on
the cloud motion.
Figure 12 shows the trajectories of the puff cloud when beads are contained in
the payload. Note that the addition of beads modifies the dynamics of the cloud,decreasing the entrainment coefficient and so increasing the particle range. Quite
generally, turbulence affects the mixing and coalescence of suspended particles anddroplets in various multiphase flows (Ghosh et al. 2005). The best fit entrainment
coefficient for the multiphase cloud was lower than for homogeneous clouds as isconsistent with the recent observations of Lai et al. (2013). Indeed, figure 12 showsthat the trajectories of the clouds II–V are best captured by the average of the jetand puff phase entrainment coefficients αavg . The average values all fall close to 0.07,
the value reported by Lai et al. (2013) in their study of a two-phase sediment cloud.
The pattern of bead deposition is shown in figure 13. Such images were analysedoto extract the maximum distance of deposition, their range xmax.
4. Theoretical model

We proceed by developing a theoretical model based on the physical picture ofviolent expiratory events as discrete emissions of droplet-laden multiphase turbulent
clouds (figure 7). This model is used to rationalize the role of the multiphase cloudin extending the observed range of the beads reported in our analogue experiments
and to yield insight into real coughs and sneezes.
Recall that the expiratory puff is treated as an isolated mass of turbulent fluid
moving through an undisturbed environment, driven initially primarily by its
momentum rather than its buoyancy; hence, it is initially closer to a multiphase

550

L. Bourouiba, E. Dehandschoewercker and J. W. M. Bush
0.20
Experiment
0.15
21 cm
0.10

y (m)
0.05

0

– 0.05

0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

x (m)

F IGURE 11. Trajectory of the turbulent buoyant puff cloud for experiment I, where nobeads are included. The insets show the superposition of the images of the shape ofthe cloud upon release (at the end of the jet phase), followed by the superposition of
the images retracing the full evolution of the cloud. The recorded trajectory is comparedwith the theoretical predictions (based on (4.7)–(4.11)) obtained using the jet entrainment
coefficient α1 = 0.17 (dashed line) as the upper bound and the average entrainment
coefficient αavg = 0.08 (dot-dashed line) as the lower bound. Both coefficients were
measured from our experiments (see table 2).

puff than a thermal. However, buoyancy cannot be entirely neglected and ultimately
influences the cloud trajectory as is evident in figure 5. Hence, we model the cloud
as a buoyant multiphase puff with initial buoyancy B0 and initial momentum I0 .
For a one-phase cloud without suspended particles, the cloud buoyancy is definedas B = Vc (ρa − ρc )g, where Vc and ρc are the volume and density of the cloud,respectively. Here, ρa is the density of the ambient air and g the acceleration due
to gravity (see figure 4h). If the cloud has initial momentum I0 and buoyancy B0 ,one expects a transition from momentum to buoyancy dominated dynamics at a timet∗ = I0 /B0 . In our analogue experiments, particles generally fall out prior to such atransition; however, small droplets may remain suspended in real cough and sneeze
clouds throughout their evolution.
We model the multiphase flow with initial densityρc (0) =NX(ρp − ρf )φn (0) + ρf ,

(4.1)

n=1

where φn (0) is the volume fraction of suspended material, specifically, the dropletsor particles of diameter dn , i.e. Nn vn /V, where vn = π dn3 /6 is the drop volume, V thecloud volume, and Nn is the number of drops or particles suspended in the cloudof diameter dn . The cloud moves forward, entraining ambient fluid of density ρaand acquiring a volume V(t) = V(0) + Va (t) provided the drops remain suspended
in the cloud. Note that the initial volume fraction of the fluid phase in real coughs
and sneezes reported previously is relatively small, e.g. 10−5 (Duguid 1946), or551

Coughing and sneezing
(a)

0.30

(b) 0.25

Experiment

0.25
0.20

y (m)

0.15

0.15
0.10

Experiment II

0.10

0

0
–0.05

– 0.05
0

0.30

y (m)

0.2

0.4

0.6

0.8

1.0

1.2

0.25

0.20

0.20

0.15

0.15
Experiment IV

0.10

0.05

0

0
0.2

0.4

0.6

0.8

1.0

1.2

– 0.05

0.2

0.4

0.6

0.8

1.0

Experiment

Experiment V

0.10

0.05

0

0

(d) 0.30

Experiment

0.25

–0.05

Experiment III

0.05

0.05

(c)

Experiment

0.20

0

0.2

0.4

x (m)

0.6

0.8

1.0

x (m)

F IGURE 12. Trajectories of the buoyant multiphase puff clouds of experiments II–V. The
theoretical cloud trajectories based on (4.7)–(4.11) obtained with the range of entrainment
coefficients displayed in both the figures and in table 2 are shown against the observed
experimental trajectories. The number N of particles used for each experiment is also
displayed, their properties being listed in table 1.

even less, e.g. 10−7 (Johnson et al. 2011). The low volume fraction (1 %) allows
us to neglect the modifications to the droplet speeds required at higher particle
concentration (Blanchette 2003).
The cloud density evolves according to


V(0)
V(0)
ρ(0) + ρa 1 −,(4.2)
ρc (t) =V(t)V(t)
or alternativelyNX
V(0)ρc (t) =φn (t)(ρp − ρf ),(ρf − ρa ) + ρa +
V
n=1

(4.3)

where the volume fraction of suspended material of diameter dn at time t is φn (t) =Nn (t)vn /V. We consider the case where all of the suspended material has the samedensity ρp . Within this framework, droplets or particles remain suspended within the
puff until their settling speed usn exceeds the mean circulation speed within the puff.
As the puff speed necessarily decreases with distance from the source, so does thesize of the suspended material: the largest drops fall out of the cloud directly afterejection, while the smallest drops and droplet nuclei can remain suspended within the
puff for a significant distance.

552

L. Bourouiba, E. Dehandschoewercker and J. W. M. Bush
(a)

2 cm

(b)

2 cm

F IGURE 13. Deposition pattern of beads in experiment III. We show (a) the original
image and (b) the treated image from which the particle counting and estimation of
maximum particle range is obtained. Here, the maximum observed deposition range isoxmax= 37 cm, as listed in table 2.

4.1. Trajectory
The trajectory of the puff is governed by the evolution of its momentum I andbuoyancy B. The buoyancy force is vertical, causing the cloud to follow a curvilinear
trajectory denoted by (s, θ) as illustrated in figure 1(b). The equations of motion forthe cloud are decomposed into horizontal x and vertical y components. The horizontal
and vertical momentum; Ix = |I| cos θ and Iy = |I| sin θ, evolve according tod|I| cos θ= 0,dt

d|I| sin θ= B(t),
dt

(4.4)

where
|I| = ρc ηr3dsds= ρc ηα 3 s3 .
dt
dt

(4.5)

Here, B(t) = Vc (t)(ρa − ρc (t))g is the cloud buoyancy, ds/ dt = w is the mean speedof the cloud along its curvilinear trajectory, V = ηr3 = ηα 3 s3 is the cloud volume(see figure 1) and |I|2 = Ix2 + Iy2 . Equation (4.4) can be solved given initial conditions,
Ix (0) = I0 cos θ0 , Iy (0) = I0 sin θ0 , B(0), V(0) = V0 , s(0) = s0 , |I|(0) = I0 and θ = θ0 .
Note that if the cloud is spherical, η = 4π /3; however, it has been observed in both
our experiments and previous studies of thermals that the cloud is better described as
an ellipsoid with radius r and height kr, leading to η = k4π /3 with k = 9/(4π ) (Scorer
1957). Hence, the volume of the cloud can be expressed as V = 3α 3 s3 . This form isused throughout the remainder of the paper and is consistent with the observations
and developments of Scorer (1957).
The total mass of the cloud increases through entrainment according to d(ρc V)/ dt =ηρa α 3 s2 ( ds/ dt). If ρa is constant, then necessarily dB/ dt = 0: entraining fluid with no553

Coughing and sneezing

buoyancy does not change the cloud buoyancy. The conservation of initial buoyancy
B(0) = B0 , which is valid prior to particle fallout leads to the following simplificationof (4.4):
Ix = |I| cos θ = I0 cos θ0

Iy = |I| sin θ = B0 t + I0 sin θ0 .

(4.6)

In this case, |I|2 = Ix2 + Iy2 and (4.5) leads toB2 t2 + I02 + 2I0 B0 sin θ0 t
ds
= 0
dt
ηα 3 ρc s3

1/2
.

(4.7)

,

(4.8)

Using (4.3), equation (4.7) becomes
B2 t2 + I02 + 2I0 B0 sin θ0 t
ds
= 0
dt
Zp + ηα 3 ρa s3

1/2

where
Zp = V0 (ρf − ρa ) + Vp (ρp − ρf )with Vp =NX

Nn vp .

(4.9)

n=1

The governing equation for θ(t) can also be obtained by manipulation of (4.6):
tan θ =

B0
t + tan θ0 .
I0 cos θ0(4.10)

After some algebra, the trajectory of the cloud (s, θ ) can be written in the form:





tan θ
1
B0 t
I02
4
+ ln tan θ +and θ = arctan,[Zp s + ηρf s /4] =2B0 cos θcos θI0(4.11)
with initial conditions (s0 , θ0 ) = (0, 0) at t0 = 0. Note that this result breaks down forvertical trajectories, that is, when θ = ±π/2.
4.2. Fallout
As the cloud entrains fluid and decelerates, its mean velocity eventually becomescomparable to the settling speed of the largest suspended droplets. We firstdescribe the particle fallout as being discrete, assuming that droplets fallout
when their settling speed matches the mean speed of the cloud ds/ dt, whichcorresponds to its characteristic internal circulation speed. The Stokes settling
speed us = gd2 /(18µ)(ρp − ρa ) was previously used to describe the settling ofdroplets emitted by violent expirations (e.g. Wells 1934, 1955). However, a moregeneral formulation is required to describe the large size range of droplets anddroplet nuclei. Specifically, inertial effects are important for large droplets, while theshortcomings of the non-slip hypothesis becomes important for particles smaller than
1 µm (Hinds 1999; Crowder et al. 2001; Lauga, Brenner, & Stone 2005). Porosityof the residual droplet nuclei in real coughs and sneezes can also substantially affecttheir settling speeds (e.g. Kajihara 1971; Edwards et al. 1998; Crowder et al. 2001).
The dependence of the settling speed us on drop size is characterized in the Appendix.

554

L. Bourouiba, E. Dehandschoewercker and J. W. M. BushUrw
yfalloutrc

F IGURE 14. (Colour online) Schematic of the spreading and settling of the particles froma cloud emitted vertically, which is roughly analogous to flows generated by patientssneezing or coughing upward while lying down in a hospital setting. A particle may
remain suspended in the cloud until its settling speed Us exceeds the speed of thedecelerating cloud U. The fallout height yfallout and range of contamination rc are bothindicated and would be given by the smallest positive root of (4.13).

4.2.1. Discrete fallout
For discrete particle fallout, the fallout criterion of droplets or particles of diameterd is thus Us = ds/ dt. In the simplified case of a vertical buoyant puff (see figure 14),with the initial condition θ0 = π/2, the governing equations becomedy
B0 t + I0
ds
=
=
.
dt
dt Zp + ρa ηα 3 y3

(4.12)

Here, the fallout criterion Us = dy/ dt leads to the following polynomial equation for
the fallout height y−(ηUs ρa α 3 )2 6 ρa ηα 3 4 Us2 Zp ηα 3 ρa 3
y +
y −
y
2B0
4
B0
+ Zp y + [I02 /(2B0 ) − (Us Zp )2 /(2B0 )] = 0.

(4.13)The smallest positive root determines the height and time of fallout for dropletswith settling speed Us . For the horizontal configuration (figure 7) the derivation is
analogous; however, the fallout criterion leads to a transcendental equation rather thana polynomial.
We first assume that, following fallout from the cloud, the particles of diameterd are in free fall with initial position (sfall (tfall ), θfall (tfall )) or (xfall , yfall ) and exit
velocity (Us cos θfall , Us sin θfall ). While the fallout time and position are important, thehorizontal range of deposition (xrange ) is of primary interest for disease transmission
by emitted pathogens. We report in table 2 the theoretically predicted fallout distance
t
t
xfall
and the maximum distance of deposition (or range) xmax
. The latter is computedtusing the free fall equations with the initial conditions prescribed by fallout at xfall
,
t
t
t
t
specifically, u = u(xfall , yfall , tfall ). The predicted range xmax is compared with thatoobserved xmax
in our analogue experiments.
We have thus far assumed that the particles or droplets remain in the cloud until
the mean cloud velocity and so its internal speed of recirculation matches the settling
speed of the particles. This implicitly assumes that the recirculation within the cloud isacting to suspend particles. We note that such is not the case at the edge of the cloud555

Coughing and sneezing(a)Growing particle-laden
turbulent puff

(b)
A(t)

r(t)

s

r(t)
Height hConcentration
Area A

Concentration

Us

Us

F IGURE 15. (a) Schematic illustration of sedimentation from a turbulent ambient as
described by Martin & Nokes (1988). Particles with settling speed less than thecharacteristic turbulent speed remain suspended unless they settle into the lower boundary
layer. The result is a suspended particle concentration that decays exponentially over atimescale h/Us . (b) Schematic of our model of continuous fallout as described by thesystem (4.15)–(4.20). Particles now settle out of the cloud continuously prior to the final
fallout time, tfallout , as predicted by the discrete fallout model.

where the recirculation is largely horizontal. Consequently, the droplets can leave thecloud before reaching the fallout length, as is consistent with our observations of a
relatively continuous deposition of particles in our analogue experiments (figure 13).
We thus modify our discrete fallout model through considering the change in buoyancy
and volume of the turbulent puff associated with continuous particle settling.
4.2.2. Continuous fallout
Martin & Nokes (1988) proposed a simple model to describe particle settling
through a turbulent fluid in a bounded container. They assume that the fluid is well
mixed and turbulent and that the particle settling speed is substantially less thanthe characteristic turbulent speed. The particles thus remain suspended unless they
settle into the lower boundary layer, in which case they sediment out. Their model
predicts an exponential decrease of the particle concentration with a timescale ofh/Us , where h is the height of the bounded container and Us is the settling speedof the suspended particles (see figure 15). Blanchette (2003) applied this model todescribe the sedimentation of bidisperse particles from riverine gravity currents, amodel further developed by Clarke, Bush & Bush (2009). Following such ideas,
we propose a model of continuous particle fallout, in which particles settle out
continuously from the cloud prior to final fallout.
Based on our observations (see also Lee et al. 2013), we assume that the settling ofparticles occurs through the lower half surface of the turbulent puff cloud with area
A(t)/2. This gives the rate of change of the number of particles within the cloud:dN= −Us A(t)C(t)/2.
dt(4.14)

Assuming that the particles are homogeneously mixed within the cloud, the particle
concentration C = N(t)/V(t); hence, we obtain
dN
3Us N
=−
,
dt
2αs

(4.15)

where we used our previous expressions for the volume and surface area of the cloud,V(t) = ηα 3 s(t)3 and A(t) = 3ηα 2 s(t)2 . Unlike in the model of Martin & Nokes (1988),

556

L. Bourouiba, E. Dehandschoewercker and J. W. M. Bushthe volume and surface area of the cloud is changing due to entrainment. We cannevertheless solve for the deposition pattern and its effect on the trajectory of thecloud by complementing (4.15) with the following system of coupled differential
equations:dsdNdV= 3ηα 3 s2
− vp
dt
dt
dt


dρc
1
dN
ds
=
vp
(ρc − ρp ) + 3ηα 3 s2 (ρa − ρc )
dt
V
dt
dt
dB
dV
d(Vρc )
= gρa
−g
dt
dt
dt
dIy
= B
dt
q
ds
= Iy2 + (I0 cos θ0 )2 /(Vρc ),dt(4.16)(4.17)(4.18)(4.19)(4.20)

where the cloud density ρc is defined in (4.2) and depends explicitly on the particle
concentration. We can then use the system (4.15)–(4.20) to derive the trajectory of thet
cloud prior to tfallout and predict the associated particle range xmax
. The latter quantities
can be obtained for either of our two fallout models, with fallout happening either
instantaneously at tfall or continuously prior to tfall .
The predicted behaviour of the particle-bearing multiphase buoyant cloud is now
compared with the experimental observations. The trajectories of the clouds infigures 11 and 12 are compared with our theoretical prediction obtained using therange of entrainment coefficients recorded in table 2. The comparison between
predictions of the discrete and continuous fallout models is shown in figure 16. Thettheoretically predicted fallout distance xfall
and the maximum distance of depositiont(or range) of the particles xmax are given in table 2 and compared with the observedomaximum distance of deposition xmax
(e.g. figure 13) for each experiment also intable 2. Note that, except for the larger beads used in experiment II, beads fell outotof the cloud in the puff phase. The comparisons of xmax
and xmax
show a satisfactory
agreement between experiments and theory for the bead deposition range. We also seethat the multiphase cloud substantially increases the range of the beads. Indeed, using
a simple ballistic calculation (in which drag is neglected), one obtains upper bounds
for the maximum deposition distance of 21, 25, 24 and 24 cm from the source forthe experiments II, III, IV and V, respectively. Accounting for drag reduces thesevalues to 4, 2, 1 and 0.5 cm, respectively. The multiphase cloud dynamics extends
these ranges to 31, 37, 63 and 66 cm, respectively (see table 2), i.e. by a factor ofup to 130 for the smallest beads. Our theoretical model of a turbulent multiphase
cloud with continuous fallout rationalizes this observed increase of range, which ismost pronounced for smaller particles more readily suspended by the cloud.
5. Application to clinical data

Our combined experimental and theoretical study has shown how particles can besuspended in turbulent clouds, their range thus being extended. We now consider the
implications of these findings for real coughs and sneezes. Cough clouds initiallycontain both buoyant gas and negatively buoyant droplets. We proceed by using ourmodel to describe their dynamics using the initial parameters of coughs recordedin the clinical literature. Initial velocity, opening angle, moisture, volume and other557

Coughing and sneezing
0.30
Continuous fallout
theory0.25
0.20

Experiment II

0.15

y (m)

0.10

Fallout

0.05
Discrete fallout
theory

0
–0.05
–0.10

0

0.2

0.4

0.6

0.8

1.0

1.2

x (m)

F IGURE 16. Comparison of observed cloud trajectory (experiment II) and trajectories
predicted on the basis of our discrete and continuous fallout models with entrainment
coefficient αavg = 0.07. The observed point of final bead fallout is indicated. Note that the
dip in the trajectory predicted by the discrete fallout model is due to its underestimating
the cloud buoyancy through neglect of the continuous fallout of negatively buoyant
particles.
1600

Number of droplets

1400
1200
1000
800
600
400
200
0
2

4

8

16 24 32 40 50 75 100 125 150 200 250 500 1000 2000

Droplet diameter ( m)F IGURE 17. (Colour online) Histogram of droplet size in coughs (reproduced from
Duguid 1946).

parameters recorded for clinical coughs are taken from various sources (see table 3).
The droplet size distribution considered is that reported by Duguid (1946) (figure 17).
Figure 18(a,b) show the trajectories of two multiphase buoyant clouds predictedby the theoretical models developed in § 4. Figure 18(a) shows the trajectory andfallout distance of a multiphase cough cloud composed of large droplets of diameter
700 µm, while 18(b) illustrates the same for a cloud with droplets of diameter

558

L. Bourouiba, E. Dehandschoewercker and J. W. M. Bush34 ◦ C with 23 ◦ C ambient
Höppe (1981)
30–35 ◦ C in mouth and nose
Morawska et al. (2009a)
RHf
85 %McCutchan & Taylor (1981)
85–100 %Morawska et al. (2009a)
ρf air exhaled
0.98 kg m−3
Picard et al. (2008)
ρd droplets
993 kg m−3 in ambient at 34 ◦ C
Nicas et al. (2005)
V
1.2 l
Gupta et al. (2009)
Estimated over 100 coughsuMean of 11.2 m s−1
Gupta et al. (2009), Tang et al. (2009)
α0.2116
Tang et al. (2009)
From opening angle θ = 23.9 ± 3.4◦
Tf

TABLE 3. Characteristics of a typical cough. Parameters include the temperature of thefluid exhaled Tf , its relative humidity RHf and the resulting density ρf computed using
the evaporation model of Picard et al. (2008). The ambient air density ρa = 1.172 kg m−3
for typical winter conditions (temperature 23 ◦ C, relative humidity 19.1 %). The dynamic
viscosity of air is 1.9 × 10−5 kg s−1 m−1 . The density of the droplets emitted ρd , the total
volume of a cough fluid V, its mean initial velocity u, and its entrainment coefficient are
estimated based on the opening angle of the cough cloud α. The initial volume fractionof droplets used is φ = 1.89 × 10−7 , with initial momentum I0 = 1.31 × 10−2 kg m s−1 andbuoyancy B0 = 2.3 × 10−3 kg m s−2 .

30 µm. Note that the height reached for the cloud containing 30 µm droplets iscomparable to ceiling heights (4–6 m) of most facilities where ventilation systems
reside. Figure 18(c) shows the mean speed of the cough cloud decaying with time(and distance) from the source. The settling speeds of droplets of size 30–700 µm
are also compared with the mean speed of the cloud, yielding the fallout time ofeach droplet group. We assume that upon leaving the cloud, the droplets or particles
have the velocity and position of the cloud at fallout. The resulting distance at which
the droplets reach the floor (estimated at 2 m below the cougher’s mouth) can thus
be computed.
The fallout and contamination ranges are compared in figure 18(d). We see againthat the cloud dynamics plays an important role in extending the range of all droplets,
particularly the smallest ones. Indeed, without the cloud entrainment, the trajectory ofthe droplets of diameter d 6 100 µm expelled at cough cloud velocity would only
be up to a dozen centimetres away from the mouth. This range is extended to metresaway from the cougher by the cloud entrainment. For example, figure 18(b) shows thatdroplets of diameter 30 µm can have a horizontal range of up 2.5 m away from thecougher due to cloud dynamics. Moreover, when relatively large particles are ejectedin the forward direction from the cloud, their range can be extended significantly.
Indeed, their higher exit velocities relative to smaller drops that fall out later can leadto a non-monotonic variation of the contamination range as a function of the droplet
diameter, as indicated in figure 18(d).
6. Summary and conclusion

We have examined the factors that influence the transport of droplet-borne
pathogens in violent respiratory expirations. We have presented experiments of humancoughs and sneezes, which demonstrate that violent respiratory flows are turbulent,
multiphase puffs. We observed that the initial emitted turbulent fluid entrains ambient559

Coughing and sneezing
Height from cougher y (m)(a) 5.0

(b)

8
7

4.0

6
Fallout

5

3.0

4
2.0

3

Fallout

2

1.0

1
0

0.5

1.0

1.5

2.0

0

2.5

Distance from cougher x (m)

Cloud
speed

1.5

2.0

2.5

3.0

1.8

Distance (m)

0.1

1.0

(d ) 2.2(c) 101

0.5

Distance from cougher x (m)Us (d )
Settling

1.4

Contamination range

1.0
0.6
Fallout length

10–3 10–2 10–1 100

101

102

103

0.2
0

100 200 300 400 500 600 700

Time (s)F IGURE 18. Application of our model to realistic cough clouds containing monodisperse
droplets. Trajectories of the clouds containing only particles of diameter (a) d = 700 µmor (b) d = 30 µm. (c) The cloud speed is compared with the settling speed of the dropletsof various sizes, thus indicating their fallout time. (d) The non-monotonic dependence ofthe range on drop diameter arises because larger suspended droplets fall out earlier in thelife of the cloud, hence with higher exit velocities than smaller drops that fall out later.

fluid in a self-similar manner, leading to the increase of its size and decrease ofits mean speed with distance from the source. We have proposed both discrete andcontinuous fallout models that draw from previous models of multiphase plume(Morton et al. 1956; Bush et al. 2003) and sedimentation dynamics (Martin & Nokes
1988). In the discrete fallout model, droplets remain suspended until their settling
speed exceeds the mean speed of the cloud, at which point they fall out. In thecontinuous fallout model, droplets leave the bottom part of the cloud continuously
at a rate prescribed by the size of the cloud and the settling speed of the droplets.
Analogue experiments were used to test the range of validity of our theoretical model.
Both the multiphase cloud trajectory and the deposition range of the analogue droplets
were compared against experimental results and showed satisfactory match. Thus, weapplied the theoretical model to an example of a clinical data set focusing on coughparameters in winter indoor conditions.
Our key findings are as follows. The turbulent multiphase cloud plays a criticalrole in extending the range of the majority of pathogen-bearing drops that accompany
human coughs and sneezes. Smaller droplets (less than 50 µm diameter) can remainsuspended in the cloud long enough for the cough to reach heights where ventilation

560

L. Bourouiba, E. Dehandschoewercker and J. W. M. Bushsystems can be contaminated (4–6 m). A droplet of diameter d = 10 µm evaporates in
0.027 s, during which it would fall a distance of approximately 0.08 mm at a settling
speed of approximately 3 mm s−1 . It would thus clearly remain suspended in a coughor sneeze cloud metres away from the cougher. Indeed, the cloud model presented
in this study predicts that the range of droplets of diameter d 6 100 µm would be
extended by a factor of 5 to more than 200 as d decreases from 100 to 10 µm.
Moreover, following fallout from the cloud, such small droplets or their associateddroplet nuclei can be resuspended by ambient air currents that usually have speeds ofthe order of ∼1 cm s−1 (e.g. Melikov & Kaczmarczyk 2012).
The ambient conditions also influence the buoyancy of the cloud and so the rangeof contamination of its suspended droplets. Changing the buoyancy from summer towinter indoor conditions can result in a variation of the range of deposition of theorder of metres for the relatively large droplets (diameter d > 50 µm) to dozens ofmetres for the smallest droplets and droplet nuclei (diameter d < 10 µm). We notealso that droplet evaporation can enhance the cloud buoyancy, increasing its vertical
momentum and thus the chances of the cloud reaching the ceiling and contaminating
the ventilation system.
Finally, we note that problems arising in disease transmission are numerous,and constitute a rich class of problems to which the tools of fluid dynamics can be
profitably brought to bear. The use of multiphase turbulent flow theory has enabled us
to revise the existing physical picture of respiratory disease transmission. Specifically,
we have developed a theoretical model that can be used with clinical data to yield
improved estimates of the range of airborne respiratory disease transmission. A morethorough integration of our theoretical model with clinical data is the subject of ourcontinuing work.
Acknowledgements

The authors thank Alex Techet at the Experimental and Hydrodynamics Laboratory
at MIT for the use of the water tank where the analogue experiments were performed.
The authors thank the MIT Edgerton Center for giving access to some of the highspeed cameras used to visualize the flows. The authors gratefully acknowledges thefinancial support of the National Science Foundation (Grant DMS-1022356) and LB
likewise the Natural Sciences and Engineering Research Council of Canada.
Appendix

For relatively large droplets, three regimes of drag are used: the turbulent regime(Re > 1000), the transition regime (0.5 < Re 6 1000) and the Stokes regime (Re < 0.5).
For the largest drops in the turbulent regime Re > 1000, Ust = 8ag/(3 × 0.44ρa )(ρp −ρa ). In the intermediate regime, Ust = µ/(2aρa ) exp[−3.07 + 0.9935J − 0.0178J 2 ], with
J = log[4(ρp − ρa )ρa (2a)3 g/(3µ2 )] (see e.g. Hinds 1999, p. 57). The Stokes settling
speed is Ust = 2ga2 /(9µ)(ρp − ρa ).
In standard conditions, the mean free path in air is λ ≈ 100 nm and its ratiowith the size of a suspended particle d defines the Knudsen number Kn = λ/d. Slip
becomes important for 0.1 . Kn , corresponding to particles with diameters d < 1 µm,
while it is negligible for particles with d > 10 µm (Lauga et al. 2005; Crowder et al.
2001; Hinds 1999; Kajihara 1971; Cunningham 1910). In between these two sizes, instandard air conditions, additional drag is felt by the droplet due to the breakdown of
the no-slip hypothesis and a correction factor to Stokes law needs to be introduced.

Coughing and sneezing

561

REFERENCESB LANCHETTE , F. A. 2003 Sedimentation in a stratified ambient. PhD thesis, MIT, Cambridge.
B OUROUIBA , L. & B USH , J. W. M. 2012 Drops and bubbles. In Handbook of Environmental Fluid
Dynamics Volume One: Overview and Fundamentals Taylor & Francis.
B USH , J. W. M., T HURBER , B. A. & B LANCHETTE , F. 2003 Particle clouds in homogeneous andstratified environments. J. Fluid Mech. 489, 29–54.
C ARDOSO , S. S. S. & W OODS , A. W. 1993 Mixing by a turbulent plume in a confined stratified
region. J. Fluid Mech. 250, 277–305.
C LARK , R. P. & DE C ALCINA -G OFF , M. L. 2009 Some aspects of the airborne transmission of
infection. J. R. Soc. Interface 6, S767–S782.
C LARKE , J. C., B USH , A. B. G. & B USH , J. W. M. 2009 Freshwater discharge, sediment transport,
and modeled climate impacts of the final drainage of glacial lake Agassiz. J. Climate 22,
2161–2180.
C ROWDER , T. M., ROSATI , J. A., S CHROETER , J. D., H ICKEY, A. J. & M ARTONEN , T. B. 2001
Fundamental effects of particle morphology on lung delivery: predictions of Stokes’ law andthe particular relevance to dry powder inhaler formulation and development. Pharm. Res. 19,
239–245.
C UNNINGHAM , E. 1910 On the velocity of steady fall of spherical particles through fluid medium.
Proc. R. Soc. Lond. A 83 (563), 357–365.
D UGUID , J. P. 1945 The numbers and the sites of origin of the droplets expelled during expiratory
activities. Edinburgh Med. J. LII (II), 385–401.
D UGUID , J. P. 1946 The size and the duration of air-carriage of respiratory droplets and droplet-nuclei.
J. Hyg. 44 (6), 471–479.
E DWARDS , D. A., B EN -J EBRIA , A. & L ANGER , R. 1998 Recent advances in pulmonary drug delivery
using large, porous inhaled particles. J. Appl. Physiol. 85, 379–385.
G HOSH , S., DÁVILA , J., H UNT, J., S RDIC , A., F ERNANDO , H. & J ONAS , P. 2005 How turbulence
enhances coalescence of settling particles with applications to rain in clouds. Proc. R. Soc.
A: Math. Phys. Engng 461 (2062), 3059–3088.
G ONNERMANN , H. M. & M ANGA , M. 2007 The fluid mechanics inside a volcano. 321–356.
G UPTA , J. K., L IN , C. H. & C HEN , Q. 2009 Flow dynamics and characterization of a cough. Indoor
Air 19, 517–525.
H INDS , W. C. 1999 Aerosol Technology: Properties, Behavior, and Measurement of Airborne Particles.
John Wiley & Sons, Inc..
H OLTERMAN , H. J. 2003 Kinetics and evaporation of water drops in air. Tech. Rep. IMAG Report
2003, Wageningen, Netherlands.
H ÖPPE , P. 1981 Temperatures of expired air under varying climatic conditions. Intl J. Biometeorol.
25, 127–132.
H UNT, J. C. R., D ELFOS , R., E AMES , I. & P ERKINS , R. J. 2007 Vortices, complex flows andinertial particles. Flow Turbul. Combust. 79, 207–234.
H UNT, G. R. & K AYE , N. G. 2001 Virtual origin correction for lazy turbulent plumes. J. Fluid
Mech. 435, 377–396.
IMF/W ORLD B ANK , 2006 Press briefing: Avian flu. In Proc. IMF/World Bank Annual Meetings
Suntec.
J OHNSON , G., M ORAWSKA , L., R ISTOVSKI , Z., H ARGREAVES , M., M ENGERSEN , K., C HAO , C.,
WAN , M., L I , Y., X IE , X., K ATOSHEVSKI , D. & C ORBETT, S. 2011 Modality of human
expired aerosol size distributions. J. Aerosol Sci. 42, 839–851.
J ONASSEN , D. R., S ETTLES , G. S. & T RONOSKY, M. D. 2006 Schlieren PIV for turbulent flows.
Opt. Lasers Engng 44, 190–207.
K AJIHARA , M. 1971 Settling velocity and porosity of suspended particle. J. Oceanogr. Soc. Japan27, 158–162.
L AI , A., Z HAO , B., L AW, A. -K. & A DAMS , E. 2013 Two-phase modeling of sediment clouds.
Environ. Fluid Mech. 1–29.
L AUGA , E., B RENNER , M. P. & S TONE , H. A. 2005 Microfluidics: the no-slip boundary condition.
Handbook of Experimental Fluid Dynamics. Springer.

562

L. Bourouiba, E. Dehandschoewercker and J. W. M. Bush

L EE , W. Y., L I , A. C. Y., L EE , J. H. W. & A SCE , F. 2013 Structure of a horizontal sediment-laden
momentum jet. J. Hydraul. Engng 139, 124–140.
L OUDON , R. & ROBERTS , M. 1967 Relation between the airborne diameters of respiratory droplets
and the diameter of the stains left after recovery. Nature 213, 95–96.
M ARTIN , D. & N OKES , R. 1988 Crystal settling in a vigorously convecting magma chamber. Nature
332, 534–536.
M C C OOL , F. D. 2006 Global physiology and pathophysiology of cough. Chest 129 (1 Suppl),
48S–53S.
M C C UTCHAN , J. W. & TAYLOR , C. L. 1981 Respiratory heat exchange with varying temperatureand humidity of inspired air. J. Appl. Physiol. 4, 121–135.
M ELIKOV, A. K. & K ACZMARCZYK , J. 2012 Air movement and perceived air quality. Build. Environ.
47, 400–409.
M ORAWSKA , L. 2006 Droplet fate in indoor environments, or can we prevent the spread of infection?.
Indoor Air 16, 335–347.
M ORAWSKA , L., J OHNSON , G., R ISTOVSKI , Z., H ARGREAVES , M., M ENGERSEN , K., C ORBETT,
S., C HAO , C., L I , Y. & K ATOSHEVSKI , D. 2009a Size distribution and sites of origin ofdroplets expelled from the human respiratory tract during expiratory activities. Aerosol Sci.
40, 256–259.
M ORAWSKA , L., J OHNSON , G., R ISTOVSKI , Z., H ARGREAVES , M., M ENGERSEN , K., C ORBETT,
S., C HAO , C., L I , Y. & K ATOSHEVSKI , D. 2009b Size distribution and sites of origin ofdroplets expelled from the human respiratory tract during expiratory activities. J. Aerosol Sci.
40, 256–269.
M ORTON , B. R., TAYLOR , G. I. & T URNER , J. S. 1956 Turbulent gravitational convection frommaintained and instantaneous sources. Proc. R. Soc. Lond. 234, 1–23.
N ICAS , M., NAZAROFF , W. W. & H UBBARD , A. 2005 Toward understanding the risk of secondary
airborne infection: emission of respirable pathogens. J. Occup. Environ. Hyg. 2, 143–154.
PAPINENI , R. S. & ROSENTHAL , F. S. 1997 The size distribution of droplets in the exhaled breathof healthy human subjects. J. Aerosol Med. 10, 105–116.
PARK , J., K ANG , S., L EE , S. & K IM , W. C. D. 2010 How respiratory muscle strength correlateswith cough capacity in patients with respiratory muscle weakness. Yonsei Med. J. 51, 392–397.
P ICARD , A., DAVIS , R. S., G LÄSER , M. & F UJII , K. 2008 Revised formula for the density of moist
air (CIPM-2007). Metrologia 45, 149–155.
R ICHARDS , J. M. 1961 Experiments on the penetration of an interface by buoyant thermals. J. Fluid
Mech. 11, 369–384.
R ICHARDS , J. M. 1965 Puff motions in unstratified surroundings. J. Fluid Mech. 21, 97–106.
S CORER , R. S. 1957 Experiments on convection of isolated masses of buoyant fluid. J. Fluid Mech.
2, 583–594.
S CORER , R. S. 1978 Environmental Aerodynamics. Ellis Horwood.
S ETTLES , G. S. 2006 Fluid mechanics and homeland security. Annu. Rev. Fluid Mech. 38, 87–110.
S OCOLOFSKY, S. A., C ROUNSE , B. C. & A DAMS , E. E. Environmental fluid mechanics: theories
and applications. In Multi-Phase Plumes in Uniform, Stratified, and Flowing Environments
American Society of Civil Engineers.
S ONKIN , L. S. 1951 The role of particle size in experimental airborne infection. Am. J. Hyg. 53,
337–354.
TANG , J. W., L I , Y., E AMES , I., C HAN , P. & R IDGWAY, G. 2006 Factors involved in the aerosol
transmission of infection and control of ventilation in healthcare premises. J. Hosp. Infect.
64, 100–114.
TANG , J. W., L IEBNER , T. J., C RAVEN , B. A. & S ETTLES , G. S. 2009 A schlieren optical studyof the human cough with and without wearing masks for aerosol infection control. J. R. Soc.
Interface 6, S727–S736.
TANG , J. W., PATH , F. R. C. & S ETTLES , G. S. 2008 Coughing and aerosols. New Engl. J. Med.
359, e19.
T ELLIER , R. 2006 Review of aerosol transmission of influenza A virus. Emerg. Infect. Dis. 12,
1657–1662.

Coughing and sneezing

563

T ELLIER , R. 2009 Aerosol transmission of influenza A virus: a review of new studies. J. R. Soc.
Interface 6, S783–S790.
T URNER , J. S. 1979 Buoyancy Effects in Fluids. Cambridge University Press.
W EBER , T. P. & S TILIANAKIS , N. I. 2008 Inactivation of influenza A virus in the environment andmodes of transmission: a critical review. J. Infect. 57, 361–373.
W ELLS , W. F. 1934 On air-born infection. Study II. Droplet and droplet nuclei. Am. J. Hyg. 20,
611–618.
W ELLS , W. F. 1955 Airborne Contagion and Air Hygiene: An Ecological Study of Droplet Infection.
Harvard University Press.
W ONG , T. W., L EE , C. K., TAM , W., L AU , J. T., Y U , T. S., L UI , S. F., C HAN , P. K., L I , Y.,
B RESEE , J. S., S UNG , J. J. & PARASHAR , U. D. 2004 Cluster of SARS among medical
students exposed to single patient—Hong Kong. Emerg. Infect. Dis. 10, 269–276.
W OODS , A. W. 2010 Turbulent plumes in nature. Annu. Rev. Fluid Mech. 42, 391–412.
YANG , S., L EE , G. W. M., C HEN , C. -M., W U , C. -C. & Y U , K. -P. 2007 The size and concentrationof droplets generated by coughing in human subjects. J. Aerosol Med. 20, 484–494.
Y U , I. T., L I , Y., W ONG , T. W., TAM , W., C HAN , A. T., L EE , J. H., L EUNG , D. Y. & H O , T.
2004 Evidence of airborne transmission of the severe acute respiratory syndrome virus. New
Engl. J. Med. 350, 1731–1739.
Z AYAS , G., C HIANG , M., W ONG , E., M AC D ONALD , F., L ANGE , C., S ENTHILSELVAN , A. &K ING , M. 2012 Cough aerosol in healthy participants: fundamental knowledge to optimize
droplet-spread infectious respiratory disease management. BMC Pulmonary Med. 12 (1), 11
doi:10.1186/1471-2466-12-11.