key: cord-0518445-cwkngfbf
authors: Olivieri, Stefano; Cavaiola, Mattia; Mazzino, Andrea; Rosti, Marco Edoardo
title: Transport and evaporation of virus-containing droplets exhaled by men and women in typical cough events
date: 2022-01-11
journal: nan
DOI: nan
sha: 45f8e543d896adcd2af1ef88aa6b96f0b3491eda
doc_id: 518445
cord_uid: cwkngfbf

The spreading of the virus-containing droplets exhaled during respiratory events, e.g., cough, is an issue of paramount importance for the prevention of many infections such as COVID-19. According to the scientific literature, remarkable differences can be ascribed to several parameters that govern such complex and multiphysical problem. Among these, a particular influence appears associated with the different airflows typical of male and female subjects. Focusing on a typical cough event, we investigate this aspect by means of highly-resolved direct numerical simulations of the turbulent airflow in combination with a comprehensive Lagrangian particle tracking model for the droplet motion and evaporation. We observe and quantify major differences between the case of male and female subjects, both in terms of the droplet final reach and evaporation time. Our results can be associated with the different characteristics in the released airflow and thus confirm the influence of the subject gender (or other physical properties providing different exhalation profiles) on both short-range and long-range airborne transmission.

problem which includes experimental tests [7, 8] as well as investigations of computational nature [9, 10, 11, 12, 13] .

In the latter framework, in order to simulate expiratory events in an accurate way, firstly we have to select proper and representative initial conditions. This means to characterise, on one hand, the airflow that is released from the mouth and, on the other, the initial size distribution of the saliva droplets that are exhaled at the same time. Both information can be obtained from the currently available data [14, 15, 16, 17, 18] , measured by performing experiments over diverse human subjects with the goal of providing a standard and representative airflow velocity profile and droplet size distribution. Looking at such evidence, however, remarkable differences are typically found in terms of the exhaled airflow between male and female subjects [14, 16, 18] . E.g., Ref. [18] reported that the initial coughing velocity is about 40% larger for men and observed similar differences for talking, while Ref. [16] performed a detailed analysis of the time-varying flow rate of cough highlighting the same quantitative discrepancy. Note that, on the other hand, such characteristic difference is not observed concerning the initial droplet size distribution [14, 15, 17] .

Motivated by such intriguing issue, here we present a numerical study on the dispersion of droplets released when coughing with emphasis on the characteristics associated with the gender of the exhaling subject. In our previous work [11, 12] we have considered such kind of expiratory event in light of its importance for the effective transmission of virus-laden droplets, focusing on the accurate description of the physical process governing the transport and evaporation of droplets [11] as well as the influence of the droplet size distribution and environmental conditions on the prediction of the final reach by direct (or short-range) transmission [12] . We now turn our attention into another aspect that can have a relevant role in the problem: the variability with respect to the human subject.

A note should be added with respect to the objective of the present work. Here, our attention is devoted to the droplet exhalation, transport and evaporation by coughing, which is generally referred to as direct (or short-range) transmission, as opposed to indirect (or long-range) transmission associated with aerosol generated especially while breathing or speaking. Recent findings suggest that indirect transmission is the primary driver of transmission in the case of SARS-CoV-2 [4, 5, 6] . Remarkably (and as it will shown in this work), along with a prevailing contribution in terms of direct transmission, to some extent coughing (and sneezing) will also act as long-range transmission [19] . Furthermore, it has to be noted that the present analysis is carried out under some simplified conditions (as detailed in the following).

The remainder of the paper is structured as follows. In Sec. 2, we introduce the governing equations and how these are solved numerically. Then, in Sec. 3 we present the main results of our analysis. Finally, in Sec. 4 we draw some conclusions and perspectives.

We consider the airflow generated during a violent expiratory event such as coughing into a quiescent environment [16] . The flow is governed by the incompressible Navier-Stokes Fig. 1 Side view of the instantaneous relative humidity field (color coded) and exhaled droplet positions (blue and orange spheres, not in scale) after 7.6 s from the starting of cough, for both the male (top) and female (bottom) subject case.

where u u u and p are the velocity and pressure fields, respectively, ν a = 1.8 × 10 −5 m 2 /s is the kinematic viscosity and ρ a = 1.18 kg/m 3 is the density of the air. The airflow emitted from the mouth is more humid than the ambient condition; instead of focusing on the evolution of the relative humidity field RH, we evolve the supersaturation field s = RH − 1, which is governed by a simple advection-diffusion equation [20] 

where D v = 2.5 × 10 −5 m 2 /s is the water vapor diffusivity. Eq. (3) is valid when the saturated vapor pressure is constant, an assumption that holds true as long as the ambient and the exhaled air temperature do not differ substantially. Note that in this simplified framework, the supersaturation field behaves passively without any dependence from the fluid temperature. In our simulations, we assume the ambient temperature to be 25 o C and a temperature of the exhaled air of 30 o C [21] . The previous set of equations are solved numerically within a domain of length L x = 4 m, width L y = 1.25 m and height L z = 2.5 m. The fluid is initially at rest, i.e. u u u(x x x, 0) = 0 0 0, and at the ambient supersaturation s(x x x, 0) = s a = RH a − 1 = −0.4 (i.e., considering the environment at a relative humidity of 60%). The exhaled air is assumed to be fully saturated [21] (i.e. s mouth = 0) and is injected through a round opening of area A mouth = 4.5 cm 2 mimicking the mouth, at a distance from the ground of z mouth = 1.6 m. The injected airflow is prescribed according to the experimental measurements reported by Ref. [16] and shown in the left panel of Fig. 2 . The duration of the exhalation is around 0.5 s with a peak velocity of 13 m/s for the male and of 7 m/s for the female subject. The resulting Reynolds numbers (based on the peak velocity and on the mouth average radius) are about 9 × 10 3 and 5 × 10 3 , respectively. For the other domain boundaries, we prescribe the no-slip condition at the bottom (z = 0) and left (x = 0) boundaries and the free-slip condition at the top boundary (z = L z ), while applying the Dirichlet condition s = s a . For both the velocity and supersaturation field, we impose a convective outlet boundary condition at the right boundary (x = L x ). Finally, periodic boundary conditions are enforced in the spanwise direction (i.e., y = 0 and y = L y ).

The exhaled droplets are modelled as an ensemble of N spherical particles of radius R i dispersed within the airflow. Since the droplet volume fraction for coughing is always smaller than 10 −5 (see e.g. [22, 19] ), we can safely neglect the backreaction to the flow and humidity fields. The dynamics of each droplet is thus governed by the well-known set of equations [23] Ẋ

for i = 1, . . . , N, and where X X X i and U U U i are the position and velocity vectors of the i-th droplet, respectively, and g g g is the gravitational acceleration. The dynamics is affected by a Brownian white-noise process η η η i and by the Stokes time

For the sake of simplicity, no Reynolds-based correction is used for the viscous drag; it was checked a posteriori that such choice is fully justified during the puff phase, when the Reynolds number experienced by the droplets is always sufficiently small, reaching a value up to O(10) only in the initial jet phase when the droplets are expelled by the mouth. In order to define the density of the i-th droplet ρ D i , we assume the droplet to be composed by a dry nucleus with density ρ N surrounded by a salty water layer with density ρ w = 9.97 × 10 2 kg/m 3 . The dry nucleus is composed by a soluble phase (NaCl) and a insoluble phase (mucus), giving an overall density which can be expressed as

where ε m is the mass fraction of soluble material (NaCl) with respect to the total dry nucleus and ρ u and ρ s the density of the insoluble (mucus) and soluble (NaCl) parts. Thus, the density ρ D i of the droplet can be computed as

where the radius of the (dry) solid part of the droplet when NaCl is totally crystallized is given by

being C the mass fraction of dry nucleus with respect to the total droplet. Finally, the droplet radius R i evolves according to a condensation model that has been successfully employed in the analysis of rain formation processes [24, 20, 25, 26] , i.e.

where C R is the droplet condensational growth rate defined as

being e sat the saturation vapor pressure and A and B two model parameters e sat = 6.1078 × 10 2 e (17.27 T /(T +237.3)) Pa,

where ε v = ε m (ρ N /ρ s ) is the volume fraction of dry nucleus with respect to the total droplet. Note that in Eqs. (11) , (10) and (12) the (ambient) temperature T is expressed in degrees Celsius. The list of physical and chemical parameters involved in the model is completed by Table 1 which provides the values selected for the present investigation.

Concerning the initial droplet size distribution, we assume here the one from Ref. [27] , still considered as a reference on the subject. Accordingly, we consider initial droplet radii approximately ranging from 1 to 1000 µm with the 95% falling between 1 and 50 µm. Droplets are set initially at rest and randomly distributed within a sphere of radius 1 cm located inside the circular pipe from which the exhaled airflow is released. The size distribution is the same for the male and female case, this choice being justified since for the droplet size distribution no significant differences between male and female subjects are reported in the current literature [14, 15, 17 ].

The system of equations is solved numerically using the DNS code named Fujin. The equations are discretised in space with the (second-order) central finite-difference method and in time with the (second-order) Adams-Bashfort scheme, except for the droplet dynamic equations (4), (5) and (9) which are advanced in time using the explicit Euler scheme with smaller substeps to deal with the numerical constraint required by very small droplets. A fast and efficient FFT-based approach is used to solve the resulting Poisson equation for pressure. The solver has been extensively validated in a variety of problems [28, 29, 30, 31, 32] , see also https://groups.oist.jp/cffu/code. In our simulations, the numerical domain is discretised with uniform grid spacing ∆ x = 3.5 mm in all directions, resulting in a total number of approximately 0.3 billion grid points. The convergence of the results with respect to the numerical and parametric settings was verified against grid refinement, number of sampled droplets and droplet release time [11, 12] .

3 Results Fig. 1 shows a typical side view of our results obtained after 7s from the respiratory event. The figure clearly shows that the resulting droplet dynamics can be macroscopically classified in two distinct behaviors: i) some droplets leave the humid and turbulent air puff released from the mouth and travel vertically within the still ambient fluid, eventually settling on the ground at a relatively short distance from their emission point, showing a predominantly ballistic motion; ii) other droplets travel for long distance within the humid air forming a cloud of so-called airborne droplets. Although the same general classification can be done for the male and female case, in the latter the distance traveled by both sets of droplets is substantially reduced. To assess the risk of virus transmission, we define the viral load as the ratio between the initial volume of a subset of droplets (e.g., settling or remaining airborne) and the total initial volume of all exhaled droplets.

The exhalation process occurring in violent air expulsions is composed by two different regimes. In the early evolution stage, air is injected into the ambient from the mouth and the resulting flow is a jet, whose dynamics is determined by the conservation of the (timedependent) momentum flux. In the late evolution stage, the cloud stops to receive momentum from the source and freely evolves in the ambient as a turbulent puff. In the puff stage, the momentum of the cloud is constant; in this case, under the hypothesis of self-similarity and following the results by Kovasznay et al. [33] we obtain the following decay law for the mean velocity of the puff:

The right panel in Fig. 2 shows the mean velocity resulting from our simulations along with the predicted scaling. Differences can be clearly noticed between the male and female case in the initial jet stage (for which we have different inlet velocity profiles), while both cases agree very well with the theoretical decay law from around half second. The different initial jet phase results in a different coefficient for the scaling law, with the male profile following the law 0.45 t −3/4 and the female one 0.19 t −3/4 . Having characterised the main features of the emitted airflow, we now move to the analysis of the droplet transmission mechanisms, focusing at first on the settling droplets. The left panel in Fig. 3 shows as a function of time the viral load of settling droplets. We can observe that in the first few seconds after exhalation, the puff rapidly loses viral load carried by large droplets to the ground. For both the male and female subjects, after around 5 s ap- proximately 99% of the viral load has reached the ground; this large percentage is however due to a very small number of droplets, around 5%, as shown in the right panel of Fig. 3 , thus indicating that 95% of droplets remain airborne after 60 s (note that, at this time, all droplets are still within the computational domain). The distance reached by the settling droplets is reported in the bottom of Fig. 3 where we show the probability density function of the distance travelled by the large droplets when they reach ground. For men, the maximum distance at the ground reached by the large droplets is 1.5 m, reducing to 1 m for women. It is worth noticing that for these cases the distance is within the social-distancing limits suggested by the World Health Organisation. Note however that, when varying the environmental condition and/or the initial size distribution, such agreement can be lost, showing that the current recommendations may substantially underestimate the safety margin [12] .

Let us now focus on the smaller droplets that experience airborne or, more specifically, indirect (or long-range) transmission. Droplets travel into a non-saturated field and thus evaporate; in Fig. 4 -left we report the temporal evolution of the supersaturation field s sampled along the Lagrangian trajectory of a typical small and large droplets. For the large droplet (dashed lines), we observe that the supersaturation field remains similar to the one inside the mouth for long time, indicating an initial slow dynamics due to the droplet inertia, followed by a rapid decay towards the ambient value as soon as the droplet leaves the turbulent puff. On the contrary, the supersaturation field felt by the small droplet (solid lines) shows the so-called "plateaux-and-cliffs" behaviour, where the scalar field displays dramatic fluctuations occurring in small regions (called cliffs) separating larger areas where the scalar is well mixed (called plateaux). Because of this, small droplets tend to remain long in the large well-mixed regions where they can equilibrate with the local value of supersaturation. The resulting dynamical process is thus made of equilibrium phases alternating with phases of rapid evaporation. Focusing on the influence of the subject gender, the same behaviour can be observed for both men (blue) and women (orange), with the latter showing an overall delay in the dynamics due to the reduced velocity of the exhaled airflow.

The delay in the dynamics for the female subject in turns affects the evaporation process. We thus quantify this feature by measuring the time τ evap needed to reach the final equilibrium size for each airborne droplet, as shown in the right panel of Fig. 4 ; in particular, we report the probability density function of τ evap for both the male and female case. In general we find broad probability density functions, which is the fingerprint of turbulent fluctuations. However, while for men the mean evaporation time is around 0.4 s, with values ranging between 0.1 s and 1 s, for women it increases to 0.6 s, with the tail of the distributions even reaching 1.5 s.

As a further consequence, the observed delay in the evaporation significantly affects the droplet motion. Indeed, we find different predictions for the viral load carried by airborne droplets for the two subjects. In particular, the left panel of Fig. 5 shows the cumulative viral load per unit area reaching the distance of 2 m from the mouth. In the male-subject simulation, considerable viral load reaches the distance covering an area of around 1 m in size, while in the female-subject case much less viral load reaches the same distance. To fully quantify the distance travelled by airborne droplets, we track the position of the center of mass of the cloud made exclusively of such subset. For the first 60 s we can directly compute this quantity from our numerical results while for longer times we extrapolate the trajectory up to the location where the center of mass eventually reaches the ground [12] . It clearly appears that, in the absence of external airflow, small airborne droplets can travel several meters: for both men and women, the cloud reaches the ground in about 20 min; due to the different horizontal speed, however, in the female case the center of mass stops at about 1.5 m while in the male case it approaches 3 m.

The present work investigates the influence of the gender of the emitting subject in the dynamics of violent expiratory events in order to better characterise the transport and evaporation process undergone by the exhaled saliva droplets, a topic of paramount relevance for devising improved safety recommendations to face the spreading of airborne infections such as the recent coronavirus outbreak. To this aim, we focus on the typical cough of both a male and female subject case, and choose a representative condition where the temperature difference between the ambient and exhaled air is such that the supersaturation field can be modelled as a passive scalar. Under such assumptions, a numerical investigation is therefore carried out by means of highly-resolved direct numerical simulations complemented by a Lagrangian model to evolve the droplets released in the flow.

Among our main findings we observe that, although the physical process is found to be essentially the same from a qualitative viewpoint, substantial quantitative differences occur between men and women in terms of the droplet final reach and evaporation time. The horizontal distance travelled by droplets is found to be generally larger for men, and so is for the cumulative viral load reaching a 2-meter distance. For women the droplet evaporation time is larger due to a slower dynamics while residing within the turbulent puff. Overall, the results can be associated to the different characteristic airflow, with a typically stronger expulsion for men.

Our findings further suggest that the current safety guidelines recommending from oneto two-meter-distance are not sufficiently conservative to protect individuals from longrange airborne transmission, either considering typical expiratory events for male or female subjects. The gender of the subject appears as another parameter (along with, e.g., environmental conditions) having an influence on the spreading of virus-containing droplets. However, it should be pointed out that the gender is typically correlated with other physical parameters, such as the weight of the subject. It cannot therefore be excluded at all that the latter has a more direct effect in the observed dynamics. Similarly, there are other parameters, such as the age or height of the subject, that can be relevant and should be properly considered to expand the existing knowledge on the topic. Finally, the present work focuses solely on cough and the extension of our findings to other expiratory events, such as sneeze or talk, is not straightforward due to the peculiar features in the airflow generated by different kinds of exhalation [18] , thus representing an open issue for future investigations.

Airborne or droplet precautions for health workers treating COVID-19?

Biological fluid dynamics of airborne COVID-19 infection

The flow physics of COVID-19

It is Time to Address Airborne Transmission of COVID-19

Mounting evidence suggests coronavirus is airborne -but health advice has not caught up

The coronavirus pandemic and aerosols: Does COVID-19 transmit via expiratory particles?

Quantity, size distribution, and characteristics of cough-generated aerosol produced by patients with an upper respiratory tract infection

An experimental framework to capture the flow dynamics of droplets expelled by a sneeze

Host-to-host airborne transmission as a multiphase flow problem for science-based social distance guidelines

Extended lifetime of respiratory droplets in a turbulent vapor puff and its implications on airborne disease transmission

Turbulence role in the fate of virus-containing droplets in violent expiratory events

Fluid dynamics of COVID-19 airborne infection suggests urgent data for a scientific design of social distancing

Short-range exposure to airborne virus transmission and current guidelines

The size and concentration of droplets generated by coughing in human subjects

Exhaled droplets due to talking and coughing

Flow dynamics and characterization of a cough

Modality of human expired aerosol size distributions

Study on the initial velocity distribution of exhaled air from coughing and speaking

Violent expiratory events: on coughing and sneezing

Droplet condensation in turbulent flows

Size distribution and sites of origin of droplets expelled from the human respiratory tract during expiratory activities

Settling velocity and concentration distribution of heavy particles in homogeneous isotropic turbulence

Equation of motion for a small rigid sphere in a nonuniform flow

Microphysics of clouds and precipitation

The equivalent size of cloud condensation nuclei

Droplet feedback on vapor in a warm cloud

The size and the duration of air-carriage of respiratory droplets and droplet-nuclei. Epidemiology & Infection

Numerical simulation of turbulent channel flow over a viscous hyper-elastic wall

Flowing fibers as a proxy of turbulence statistics

Droplets in homogeneous shear turbulence

Increase of turbulent drag by polymers in particle suspensions

Dispersed fibers change the classical energy budget of turbulence via nonlocal transfer

Unsteady turbulent puffs

The authors acknowledge the computational time provided by HPCI on the Oakbridge-CX cluster under the grant hp200157 of the "HPCI Urgent Call for Fighting against COVID-19" and the computer time provided by the Scientific Computing section of Research Support Division at OIST. The research was supported by the Okinawa Institute of Science and Technology Graduate University (OIST) with subsidy funding from the Cabinet Office, Government of Japan. A.M. thanks the financial support from the Compagnia di San Paolo, project MINIERA no. I34I20000380007. Useful discussions with Agnese Seminara are warmly acknowledged.