Transport and evaporation of virus-containing droplets exhaled by men and women in typical cough events.

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. © Springer Nature B.V. 2022.

Introduction

In light of the on-going COVID-19 pandemic, the understanding of how virus-containing droplets exhaled during respiratory events (e.g., when coughing, sneezing or talking) travel in the surrounding environment has gained considerable importance [1-3]. An intense debate has raised within the scientific community on the possibility that airborne transmission could play a relevant role in the spreading of SARS-CoV-2 and similar infections [4-6], resulting in a renovated research effort in the analysis of this complex and cross-disciplinary problem which includes experimental tests [7, 8] as well as investigations of computational nature [9-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-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 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-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 Sect. 2, we introduce the governing equations and how these are solved numerically. Then, in Sect. 3 we present the main results of our analysis. Finally, in Sect. 4 we draw some conclusions and perspectives.

Problem formulation

Airflow dynamics

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 equationswhere and p are the velocity and pressure fields, respectively, is the kinematic viscosity and 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 , we evolve the supersaturation field , which is governed by a simple advection-diffusion equation [20]where is the water vapor diffusivity. Equation (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 °C and a temperature of the exhaled air of 30 °C [21]. The previous set of equations are solved numerically within a domain of length , width and height . The fluid is initially at rest, i.e. , and at the ambient supersaturation (i.e., considering the environment at a relative humidity of ). The exhaled air is assumed to be fully saturated [21] (i.e. ) and is injected through a round opening of area mimicking the mouth, at a distance from the ground of . 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 with a peak velocity of for the male and of for the female subject. The resulting Reynolds numbers (based on the peak velocity and on the mouth average radius) are about and , respectively. For the other domain boundaries, we prescribe the no-slip condition at the bottom () and left () boundaries and the free-slip condition at the top boundary (), while applying the Dirichlet condition . For both the velocity and supersaturation field, we impose a convective outlet boundary condition at the right boundary (). Finally, periodic boundary conditions are enforced in the spanwise direction (i.e., and ).

Fig. 2

(Left) Time-varying inlet airflow velocity of cough according to Ref. [16]. The blue and orange lines are used to distinguish the velocity profiles of a male and female subject, respectively. (Right) Time history of the mean streamwise velocity component v. The lines show the expected scaling while the symbols are the results of our simulations. (Color figure online)

Droplet dynamics

The exhaled droplets are modelled as an ensemble of N spherical particles of radius dispersed within the airflow. Since the droplet volume fraction for coughing is always smaller than (see e.g. [19, 22]), 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 , and where and are the position and velocity vectors of the i-th droplet, respectively, and is the gravitational acceleration. The dynamics is affected by a Brownian white-noise process 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 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 , we assume the droplet to be composed by a dry nucleus with density surrounded by a salty water layer with density . The dry nucleus is composed by a soluble phase (NaCl) and a insoluble phase (mucus), giving an overall density which can be expressed aswhere is the mass fraction of soluble material (NaCl) with respect to the total dry nucleus and and the density of the insoluble (mucus) and soluble (NaCl) parts. Thus, the density of the droplet can be computed aswhere the radius of the (dry) solid part of the droplet when NaCl is totally crystallized is given bybeing the mass fraction of dry nucleus with respect to the total droplet. Finally, the droplet radius evolves according to a condensation model that has been successfully employed in the analysis of rain formation processes [20, 24–26], i.e.where is the droplet condensational growth rate defined asbeing the saturation vapor pressure and A and B two model parameterswhere 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.

Table 1

Physical/chemical parameters representative of expiratory events and adopted in the present investigation

Density of soluble aerosol part (NaCl)	ρ_s	2.2×10^3kg/m^3
Density of insoluble aerosol part (mucus)	ρ_u	1.5×10^3kg/m^3
Density of dry nucleus	ρ_N	1.97×10^3kg/m^3
Mass fraction of soluble material (NaCl) w.r.t. the total dry nucleus	ϵ_m	0.75
Mass fraction of dry nucleus w.r.t. the total droplet	C	1%
Specific gas constant of water vapor	R_v	4.6×10^2J/(kgK)
Heat conductivity of dry air	k_a	2.6×10^-2W/Km
Latent heat for evaporation of liquid water	L_w	2.3×10^6J/kg
Saturation vapor pressure	e_sat	0.616kPa
Droplet condensational growth rate	C_R	1.5×10^-10m^2/s
Surface tension between moist air and salty water	σ_w	7.6×10^-2J/m^2
Molar mass of NaCl	M_s	5.9×10^-2kg/mol
Molar mass of water	M_w	1.8×10^-2kg/mol
Number of ions into which a salt molecule dissociates	n_s	2
Osmotic coefficient of salt in solution	Φ_s	1.2

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 with the falling between 1 and . Droplets are set initially at rest and randomly distributed within a sphere of radius 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].

Numerical discretisation

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-32], see also https://groups.oist.jp/cffu/code. In our simulations, the numerical domain is discretised with uniform grid spacing 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].

Results

Fig. 1 shows a typical side view of our results obtained after 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.

Fig. 1

Side view of the instantaneous relative humidity field (color coded) and exhaled droplet positions (blue and orange spheres, not in scale) after from the starting of cough, for both the male (top) and female (bottom) subject case. (Color figure online)

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 (time-dependent) 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 and the female one .

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 approximately of the viral load has reached the ground; this large percentage is however due to a very small number of droplets, around , as shown in the right panel of Fig. 3, thus indicating that of droplets remain airborne after (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 , reducing to 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].

Fig. 3

Top: Cumulative viral load (left) and normalized droplet number (right) settling to the ground as a function of time (both quantities are expressed as percentage). Bottom: Probability density function of the distance from the mouth when droplets reach the ground. In the figures, the blue and orange colors distinguish the results for male and female, respectively

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.

Fig. 4

(Left) as a function of time experienced by two representative droplets. The solid and dashed lines are used to distinguish a small and large droplet. (Right) Probability density function of the droplet evaporation time , i.e. the time for the droplet to shrink to its final radius. The observation time is

The delay in the dynamics for the female subject in turns affects the evaporation process. We thus quantify this feature by measuring the time 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 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 , with values ranging between and , for women it increases to , with the tail of the distributions even reaching .

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 from the mouth. In the male-subject simulation, considerable viral load reaches the distance covering an area of around 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 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 while in the male case it approaches .

Fig. 5

(Left) Cumulative viral load per unit area reaching a distance of from the mouth after . (Right) Trajectory of the viral load center of mass (computed considering only the airborne droplets): horizontal position (solid line) and vertical position (dash-dotted line). The colored lines indicate the results from the simulation while the black dotted ones are extrapolations over longer times

Conclusions

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 one- to two-meter-distance are not sufficiently conservative to protect individuals from long-range 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.