Air change rate effects on the airborne diseases spreading in Underground Metro wagons

Metrics details

Abstract

The effect of the rate of change of fresh air inside passengers’ wagons for Underground Metro on the spreading of airborne diseases like COVID-19 is investigated numerically. The study investigates two extreme scenarios for the location of the source of infection within the wagon with four different air change rates for each. The first scenario considers the source of infection at the closest point to the ventilation system while the other places the infection source at the farthest point from the wagon ventilation system. The effect of the wagon windows’ status (i.e. closed or open) is also studied. It is found that under all conditions, open windows are always favored to decrease the infection spreading potential. A higher air change rate also decreases the infection spreading up to a certain value, beyond which the effect is not noticeable. The location of the infection source was found to greatly affect the infection spreading as well. The paper gives recommendations on the minimum air change rate to keep the infection spreading potentials to a minimum considering different times the passengers stay in the wagon.

Introduction

The Middle East and North Africa (MENA) region is one of the hottest places on earth with temperatures up to 50 °C during summer (Krarti and Ihm 2016 ; Hijazi and Howieson 2018 ). Egypt, as part of the MENA region, has an average annual cooling degree-day of about 2500 for a base temperature of 22 °C and an average heating degree-day of about 300 for a base temperature of 18 °C (Hijazi and Howieson 2018 ). This high cooling and heating load demand necessitates the use of ventilation equipment to ensure that the proper temperature and humidity conditions are met (Gupta et al. 2009 ). Metro stations are no exception to this fact with stricter constrains on the used ventilation mechanism duo to its huge impact on the Indoor Air Quality (IAQ) (Kim et al. 2012 ; Kwon et al. 2015 ; Martins et al. 2015 ; Qiao et al. 2015 ).

Egypt is one of the most populous countries in the MENA region with an estimated population of 100 million (Ahmed et al. 2020 ). Cairo has an estimated population of 20 million as of 2016 making it the city with the highest population in the MENA region (Ahmed et al. 2020 ). It is estimated that 97% of the population in Egypt are using public transportation (Ahmed et al. 2020 ). The reason for this high percentage is the standard of living in addition to the incentives given by the government to promote the use of more energy efficient and environmentally friendly modes of transportation (Ahmed et al. 2020 ).

Metro is considered as one of the most plausible modes of transportation in Greater Cairo (GC) (Awad 2002 ; Eldeeb et al. 2018 ; Eldeeb et al. 2019 ; Ahmed et al. 2020 ). The Greater Cairo Underground Metro (GCUM) is the main mode of transportation used by residents and visitors of the GC area with a capacity of up to 50,000 passengers/hour for some lines during peak hours (Eldeeb et al. 2019 ).

Several studies have investigated the importance of IAQ for the wellbeing of people (Ghanizadeh and Godini 2018 ; Guo et al. n.d. ) while others considered the effect of the IAQ for metro stations on the transport of particulate matter and airborne diseases (Kim et al. 2012 ; Hijazi and Howieson 2018 ; Guo et al. n.d. ; Yan et al. 2015 ). As recently confirmed by the World Health Organization (WHO) and several studies, COVID-19 is one of the airborne diseases which poses a great danger and high infection probability for crowded places with mechanical ventilation systems like metro stations (Correia et al. 2020 ).

A recent review suggests the modification of the ventilation system as a solution to mitigate the effect of the airborne pathogens by increasing the fresh air flow rate and the use of better air filters (Eldeeb et al. 2019 ; Wen et al. 2020 ).

Despite the numerous studies about the IAQ of metro stations in Egypt and other countries, the studies are limited to the stations and platforms whether underground or surface stations. An investigation about the effect of the ventilation systems and the rate of fresh air change on the spreading of airborne pathogens inside the passengers’ wagons in Egypt is missing.

During the outbreak of the COVID-19 and after the lockdown measures have gradually relaxed, it is expected that many people will spend more time in ventilated public transportations or even buildings, so it is quite needed to discuss and study the flow field effects of the HVAC systems in public transportations mainly subways since it is one of the most crowded places.

Morawska et al. n.d. ) proposed few ways to avoid the buildup of the viral contamination by avoiding air recirculation and increasing the ventilation rates to avoid the airborne transmission of the virus, though it discussed mainly the indoor transmission, but it was very important to study the exceptional measures to decrease the infection possibilities in the public transportation too, including subways.

Also, regarding the ventilation, the same conclusions were obtained in (Kohanski et al. 2020 ; Marcone 2020 ) where it assured the importance of introducing more flow rates in the ventilation systems in addition to extra filtration for the sources of the flow.

Chun et al. ( 2014 ) introduced a simple mathematical model to predict the number of infected people depending on few simple parameters, like number of sensitive people present in the environment, number of infected people capable of transmitting the virus, and mainly regarding the ventilation systems like the room fresh airflow rate.

In (Gupta et al. 2009 ), a set of equations was developed that can be used to generate boundary conditions for predicting infectious virus transport by CFD due to coughing. The work focused on the flow rate and direction and air velocity that can be determined from the human size and mouth opening area. The inputs required are height, weight, and gender of a person.

In (Zhu et al. 2006 ), a numerical analysis is conducted to analyze the indoor flow field assuming coughing and respiration. They assumed coughing to be a steady phenomenon, with a velocity equal to the maximum initial velocity of the coughed airflow determined by PIV experiments as the boundary condition. Then, a parametric study is conducted to analyze the transport process for droplets with diameters of 30, 50, 100, 200, 300, and 500 micrometers using the Lagrangian equation to show how saliva droplets were expelled from the mouth by coughing, dispersed in the air, and finally attached to surfaces. Based on the results, the transport characteristics of saliva droplets produced by coughing are investigated in a calm indoor environment. The patient human model was assumed to cough consecutively, while the other human model was assumed to inhale constantly. Andrade et al. ( 2018 ) analyzed the risk of infection for practicing persons in a gymnasium considering different gyms, namely, ventilated with either split system or central system air conditioners. The study advised that ineffective ventilation increases the risk of infections such as influenza and tuberculosis. Abbasi and Samaei ( 2019 ) studied the effects of the incubation temperature on the density and composition of airborne fungi in an indoor and outdoor space of hospital. They concluded that the incubation temperature had a remarkable effect on airborne fungi.

In (Khare et al. 2015 ), a study on covering a cough is conducted experimentally. It is found that it can be useful in reducing the transmission of airborne infectious diseases. In order to do this, smoke is used to visualize the airflow exhaled by as much as 16 people. Their mouths were covered by several means, as cupped hands, fists, elbows, and tissues. Then that study developed simplified models for predicting the airflow based on smoke visualization data. It was found that covering a cough with a tissue, a cupped hand, or an elbow can significantly reduce the exhaled air horizontal velocity and cause the particles to move upward with the thermal plumes generated by a human body. In contrast with an uncovered cough, a covered cough or a cough with the head turned away may prevent direct exposure.

Mathematical modeling is used for the prediction of the spread of the virus whether on a global scale like national or even international scale using statistical models and also is used to predict the spread of the virus over building scales by tracing the path of the infection particles using differential equations. For a review on different models, see Kohanski et al. ( 2020 ).

Regarding the mathematical modeling of airborne diseases modelled as aerosols, one could think of purely Eulerian techniques (Hijazi and Howieson 2018 ) or mixed Eulerian–Lagrangian techniques (Kim et al. 2012 ; Guo et al. n.d. ) to be able to capture or track the particles causing the infection, respectively. Using OpenFOAM package, Pendar and Páscoa ( 2020 ) applied a fully coupled Eulerian–Lagrangian technique resulting in a deeper understanding of the saliva-disease-carrier droplet transmission mechanisms. Dbouka and Drikakis ( 2020 ) used a mathematical model considering the multiphase nature of the phenomenon and considered heat transfer effects to investigate transport, dispersion, and evaporation of saliva particles and they concluded that a further work is required to quantify the influence of the environment’s relative humidity (HR) and temperature. The effect of HR and the ambient wind on the value of the “physical-distancing” is investigated using Eulerian–Lagrangian formulation by Feng et al. ( 2020 ) using computational fluid-particle dynamic (CFPD) simulations with the capability of predicting droplet size changes due to condensation and evaporation.

The goal of this paper is to study the transport of the cough-generated aerosol from a patient in a metro wagon subjected to different scenarios of ventilation rates and patient position and make recommendations on the ventilation for reducing the risk of infection of COVID 19 and the similar airborne viruses.

The paper is structured as follows: in “The mathematical model”, the mathematical model of the problem is stated, followed by describing the geometry, mesh, and the boundary conditions of the cases in “Problem setup”. Results are discussed in “Results and discussion”, and conclusions are drawn in “Summary and conclusions”.

The mathematical model

The problem at hand is modeled as unsteady, compressible, turbulent, and non-isothermal ideal gas flow with dispersed incompressible constant diameter cough-generated aerosols. Fig. 1 shows a simplified schematic for the interaction between the gas and the cough-generated aerosols.

Fig. 1

and the turbulent heat flux is:

$$ {\mathbf{q}}_t=-{h}_t\mathbf{\nabla}T, $$

(11)

With the Boussinesq hypothesis, the turbulence problem is reduced to the determination of two variables, namely, the turbulent viscosity μt and the turbulent thermal conductivity ht which are determined from appropriate turbulence model specifically in this work k-ω model.

The evolution equation for the turbulent kinetic energy is (Couaillier 2000 ):

$$ {\rho}_a\left(\frac{\partial k}{\partial t}+{\mathbf{V}}_{\boldsymbol{a}}\cdotp \nabla k\right)={\boldsymbol{\uptau}}_t:\mathbf{D}+\nabla \cdotp \left[\left(\mu +{\mu}_t\right)\nabla k\right]-{\rho}_a\varepsilon -R, $$

(12)

and for the dissipation rate of k is:

$$ {\rho}_a\left(\frac{\partial \varepsilon }{\partial t}+{\mathbf{V}}_{\boldsymbol{a}}\cdotp \nabla \varepsilon \right)=1.57\frac{\varepsilon }{k}{\boldsymbol{\uptau}}_t:\mathbf{D}+\nabla \cdotp \left[\left(\mu +\frac{\mu_t}{1.3}\right)\nabla \varepsilon \right]-2c\frac{\varepsilon^2}{k}+E, $$

(13)

with

\( R=2\mu {\left(\nabla \sqrt{k}\right)}^2,E=\frac{2\mu {\mu}_t}{\rho }{\left(\frac{\partial^2{V}_t}{\partial {n}^2}\right)}^2 \) , \( c=1-0.3\mathit{\exp}\left(-{R_{et}}^2\right),{h}_t=\frac{c_{p{\mu}_t}}{{\mathit{\Pr}}_t} \)and \( {R}_{et}=\frac{\rho {k}^2}{\mu \varepsilon} \).

The turbulent viscosity is by (Couaillier 2000 ):

$$ {\mu}_t=0.09{R}_{et}\mathit{\exp}\left(\frac{-2.5}{1+0.02{R}_{et}}\right) $$

In the above equations, ρa is the air density, Va is air velocity vector, p is the thermodynamics pressure to be determined from the perfect gas equation of state, and m is the mass of the water droplet.

Problem setup

Geometry

The geometry of the used metro wagon is illustrated in Fig. 2 . The volume of the available space is 58 m3. It contains three HVAC openings, with area of 0.8 m2 each and a distance of 1.18 m between each other. The wagon has four windows, open with 15 degrees.

Fig. 2

Full size table

Results and discussion

In order to quantify the risk of COVID-19 propagation, a case study is simulated as follows. At time 0, one patient coughs, exhaling air-water mixture, while other people breathe normally, and the ventilation system continuously supplies fresh air. It is assumed that water is initially dispersed, and the infection is carried on these droplets. Hence, the water content in the whole metro wagon is used to identify the presence of the infection. Namely, the ratio of water volume divided by the whole volume of the wagon is measured in this study. The simulations are done using ANSYS CFX on 64 cores. The computational time for one simulation is 576 hours/core.

Three sets of numerical simulations were conducted. The first set simulates the patient to sit close to the ventilation suction window, which represents the best-case scenario, and the flow rate of the ventilation system changes from 16 to 40 CFM per passenger, corresponding to air changes per hour (ACPH) of 0.47 to 1.17 which were calculated using Equation ( 23 ).

$$ \mathrm{ACPH}=60\ast \mathrm{CFM}/\mathrm{volume} $$

(23)

The second set represents the case where the patient sits away from the ventilation suction window, which represents the worst-case scenario, and the ventilation flow rate changes also. The third and fourth cases study the window’s effect, by comparing cases where the wagon windows are closed with those where the windows are open, in both the best and worst cases. In these cases, the ventilation flow rate changes from 16 to 24 CFM. The simulated cases are summarized in Table 4 .

Table 4 Investigated cases

Full size image

Summary and conclusions

The effect of ventilation systems on the spread of airborne diseases like COVID-19 is numerically investigated using Eulerian–Lagrangian approach considering turbulence and thermal effects. The infection particles are modeled as being aerosol particles. From the simulation results, the following conclusions are clear:

1.

Opening the windows helps in changing the air inside the wagon. Quantitatively speaking, opening windows reduced the time required to reach 50% WVR by 40%, for the same ACPH.

2.

By increasing the ventilation inlet flow rate, the infection spreading decreases considerably. In the current simulation, increasing the ACPH from 0.47 to 0.7 reduced the WVR from 19% to 6% after 10 minutes.

3.

As the HVAC inlet flow rate increases beyond a certain value, the infection hazard does not show much difference. Changing the ACPH value from 0.7 to 0.94 reduced the WVR from 6% to 2% after 10 minutes.

4.

The position of the infected person relative to the ventilation windows has a major role in the infection spreading. The closer the infected person to the suction ventilation window is, the faster the cough generated aerosols leave the wagon. In the best-case scenario, after 10 minutes from the cough, the WVR values are 19% and 35% for the cases where the patient is close to the ventilation suction and inlet windows, respectively.

5.

Finally, this work recommends keeping the windows open and increase the volume flow rate above the average.

Data availability

Not applicable.

References

Abbasi F, Samaei MR (2019) The effect of temperature on airborne filamentous fungi in the indoor and outdoor space of a hospital. Environ Sci Pollut Res 26:16868–16876