Necessity of ventilation for mitigating virus transmission quantified simply - arXiv.org
←
→
Page content transcription
If your browser does not render page correctly, please read the page content below
Necessity of ventilation for mitigating virus transmission quantified simply Eric G. Blackmana,b,∗, Gourab Ghoshala a Department of Physics and Astronomy, University of Rochester, Rochester, NY, 14627, USA b Laboratoryfor Laser Energetics, University of Rochester, Rochester NY, 14623, USA Summary arXiv:2006.11651v3 [physics.soc-ph] 18 Aug 2020 Background To mitigate the SARS-CoV-2 pandemic, officials have employed social distancing and stay- at-home measures, with increased attention to room ventilation emerging only more recently. Effective dis- tancing practices for open spaces can be ineffective for poorly ventilated spaces, both of which are commonly filled with turbulent air. This is typical for indoor spaces that use mixing ventilation. While turbulence initially reduces the risk of infection near a virion-source, it eventually increases the exposure risk for all occupants in a space without ventilation. To complement detailed models aimed at precision, minimalist frameworks are useful to facilitate order of magnitude estimates for how much ventilation provides safety, particularly when circumstances require practical decisions with limited options. Method Applying basic principles of transport and diffusion, we estimate the time-scale for virions injected into a room of turbulent air to infect an occupant, distinguishing cases of low vs. high initial virion mass loads and virion-destroying vs. virion-reflecting walls. We consider the effect of an open window as a proxy for ventilation. Findings When the airflow is dominated by isotropic turbulence, the minimum area needed to ensure safety depends only on the ratio of total viral load to threshold load for infection. Interpretation The minimalist estimates here convey simply that the equivalent of ventilation by modest sized open window in classrooms and workplaces significantly improves safety. Keywords: elsarticle.cls, LATEX, Elsevier, template 2010 MSC: 00-01, 99-00 Introduction The SARS CoV-2 virus, first reported in 2019 [1, 2] has since spread to at least 213 countries and territories leading to an unprecedented global pandemic [3]. The lack of therapeutics and vaccines have led public health officials to employ non-pharmaceutical interventions focused on physical distancing measures and masks, but comparatively much less on ventilation. With the resumption of visits to offices, bars, restaurants, salons and universities, where people are expected to be in close proximity in for prolonged periods of time, ventilation and HVAC filtering are safety precautions that must be prioritized [4, 5]. Airborne viral particles like SARS-CoV-2, which cause Covid-19, get trapped on moisture droplets that carry the virions [6]. Interior air is commonly, if not unavoidably, turbulent. This keeps droplets airbone much longer than their free-fall time [7]. The droplets form a spectrum of sizes, and those below 5µm can remain airborne in aerosols for at least 3 hours [8]. Since virions are attached to these aerosol droplets, basic principles of turbulent diffusion and transport or aerosols [9, 10, 11, 12, 13, 14, 15, 16] become directly applicable to guiding HVAC issues and the efficacy of masks [17]. Forced central air/heating in HVAC systems that use mixing ventilation without sufficient replenishment of fresh air exacerbates the danger of airborne virions. The rapid spread of aerosol droplets in room of turbulent air is exemplified by spraying ∗ Correspondingauthor Email addresses: blackman@pas.rochester.edu (Eric G. Blackman), gghoshal@pas.rochester.edu (Gourab Ghoshal) Preprint submitted to Journal of LATEX Templates August 20, 2020
injected mass of virions M mass of virions to infect one person Mc # of people injected load can infect N = M/Mc rate of virions encountering a face Ṁh density of virions uniformly spread over room ρ room radius R eddy size l turbulent eddy air speed vl eddy turnover time tl = l/vl eddy diffusivity νT = vl l/3 diffusion speed over scale x ≥ l vdif (x) = νT /x mean air flow vp face width h face area Ah = h2 required time for infection tc window area W time for load to diffuse via window tsaf Table 1: list of variables a scented aerosol and measuring how quickly a person on the other side of the room can detect it. The benefits of physical distancing are diminished by prolonged exposure to viral filled turbulent air in a closed room because virions are transported throughout by turbulent diffusion of the host droplets. Without a tightly sealed mask and proper eye protection, accumulated indoor-exposure is likely. Virions can also be transported by HVAC systems between rooms. People staying home may be exposed to the virus in poorly ventilated apartment buildings with forced circulating air. Staying at home in isolated houses is not equivalent to staying home within a poorly ventilated apartment building. Evidence for such non-local transport of viruses has been found in restaurants in China [18, 19], a call-center in South Korea [20], and a choir-setting in Washington state [21]. The study of airbone disease transmission and ventilation has a empirical and computational history [22, 23, 24, 25]. Much about the transmission modes of SARS CoV-2 remains unknown, such as the viral load required to cause infection [26]. Nevertheless, it is crucial to use basic principles that we do know to inform policy choices to reduce the risk of transmission [4]. As such, here we use employ basic concepts of turbulent transport to provide minimalist order-of-magnitude estimates to quantify the efficacy of ventilation. Specifically, we estimate the time scale for an individual to be infected in an enclosed space of arbitrary size, subjected to an injection of virions. We distinguish between virion-absorbing vs. virion-reflecting walls and cases in which the injected mass of virions is sufficient vs. insufficient to infect a person over one diffusion crossing time from the virial injector to the room boundary. We consider the role of an open window as a proxy for ventilation, and estimate the typical cross-section needed for safety. We also discuss the role of a systematic drift velocity in the room. Since the threshold viral load for severity of infection is unknown, despite the different loads produced by talking, breathing and other modes [27, 28, 29], we present our results simply in terms of a dimensionless quantity–the virion mass required to infect a single individual. Our results demonstrate simply the importance of ventilation and filtering, complementing, but in agreement with, more detailed computational [7, 18, 30] and empirical efforts [31]. 2
Research in Context Evidence before this study Coughing, sneezing and breathing release moisture droplets onto which SARS-CoV-2 bind. Droplets of size below 5µm in aerosols remain viable and airborne with a half life of three hours [8]. Interior air is commonly turbulent and diffuses aerosols throughout a room and building over much shorter time scales, spreading the infectious virions over distances much larger than the 6ft recommended physical distancing separation. Added value of this study We quantify the mitigating role of ventilation using simple from simple considerations of fluid transport and diffusion, by comparing the load of virions encountering a persons face in a room with and without a window. Specifically, this becomes a determination of the size of the window or filter needed to maintain safety. The results help guide practical decisions about safety and interior ventilation in spaces where people spend extended periods, such as classrooms and office buildings. Implications of all the available evidence Proper ventilation mitigates the airborne spread of SARS-CoV-2 and needs to be viewed along with masks and physical distancing as the third pillar of prevention, We show its efficacy quantitatively in a simple broadly accessible way. Results While progressively smaller droplets take longer and longer to fall out of the air [31] after the initial load (cough, breathing, or sneeze), we are mainly concerned with particles that remain suspended in the air after it is well mixed, namely these ≤ 5µ that can remain airborne for 3hr [8]. Such particles are equivalent carriers of virions for present purposes. We assume that a combined mass M of such viron-loaded aerosol droplets (VLADs) are injected in a room of radius R, and that the room has steady turbulent air of local eddy air speed vl and eddy scale l
Small viral injection mass, no ventilation The likely more commonly important case is when the initial injection of viral mass is small enough so that the room is fully mixed with VLADs before any individual is infected. First, we consider two sub-cases without ventilation: (i) completely absorbing walls (say if the walls are infused with anti-viral material such as copper fibers) and (ii) completely reflecting walls. For the former, sub-case virions are removed upon contact with the wall and no one in the room is infected. For the latter sub-case, the virion-carrying droplets remain airborne. Then, given that the accumulated mass of virions on a single person is Ṁh ' ρvdif Ah , (3) where 3M ρ= , (4) 4πR3 and vdif is the diffusion speed, which is simply the turbulent diffusion coefficient νT divided by the net displacement from the source to point of measurement. Given that we are in the well-mixed regime, the relevant displacement is just the eddy scale, since any new virion-mass accumulates via neighboring eddies. Therefore, 1 vdif = vl , (5) 3 Combining Eqns. (3), (4), and (5) into Eq. (1) yields, for the critical infection time, tl M/Mc −1 R 3 l −1 h −2 tc = 13.11 0.75 s 50 5m 0.75m 0.2m min, (6) where we have used fiducial values to facilitate estimates. Note, that while Mc is unknown, in all cases it enters as as fraction of the input viral mass. In the upper panel of Fig. 1 we plot Eq. (6) as a function of the room radius R and injected load N = M/Mc . By definition, to prevent infection an occupant must spend a duration t ≤ tc in the room. Role of an open window To quantify the effect of ventilation, we consider window of open area W > R2 /νT can be estimated as ML (t) = ṀW t. Dividing this by M and using equations (4) and (7) gives the fraction of VLAD mass lost from the room as 1 vl t W ML (t)/M ' . (8) 4π R R2 Safety would be achieved once ML ≥ M (1 − Mc /M ), corresponding to ≥ 99% replacement of the room air when Mc le0.01M . (If using a ventilation system measured in conventional ventilation units of air changes per hour (ACH) [? ], for Mc /M < 0.01 this would require > 10ACH to achieve this replacement in 1 hour since each ACH replaces only 63.2% of the air.) Solving Eq. (8) for the associated t = tsaf gives, 4πR3 tsaf = . (9) vl W 4
Now we set Eq. (6) equal to Eq. (9) to get the minimum open window area Wc ensuring tsaf < tc , so that enough virions are removed before anyone gets infected. This gives h 2 M/Mc 2 Wc ≥ 2 0.2m 50 m . (10) This is independent of the room size because the same viral mass in a larger room reduces the viral density. Therefore the flux of mass incident on a face, as well as on the surface of a window is reduced by the same factor. In the lower panel of Fig. 1b we show tc /tsaf as a function of W and N . The intersection of the curved surface with the plane is equation (10). Given that ρ is uniform in the room for the fully mixed case by diffusive action of turbulence, the ratio of mass flux through a window to that encountering a person’s face is approximately W/h2 . Increasing this ratio enhances the probability for viral mass encountering the window. Our example shows that the strategy for safety is to increase this ratio such that the leftover virial density in the room is insufficient to produce a critical infectious load for individuals during their period of stay. What is the role of masks in the context of the above estimates? They reduce the surface area of the face by at least a factor of 2 since droplets must slide around the smaller open area around the edges of the mask. A factor of ≥ 2 reduction in h2 has the equivalent effect of reducing the required size of the window area for safety by that same factor as seen from Eq. 10. A generalization of interest is the case in which there is a systematic bulk flow in the from one side to the other, for example from a pressure gradient or fan. In this case, the ratio of the flux through the window to that onto a face depends on the angles of the bulk flow to the window, θw and to the face, θf . If the bulk flow is aimed toward the window, then only the latter angle enters and this ratio of fluxes becomes (vT +vP )W (vT −vP cos θf )h2 where cos θf is the angle between the flow and the normal to the face. If |vP cos θf | >> vT , then this ratio becomes simply −W/(h2 cos θf ) highlighting that simply having ones back to the flow is the safer orientation, as the flux ratio is negative when the flow passes from back to front. In practice detailed modeling of flow around the head must be considered to assess the influence of boundary layers. Discussion Turbulence is hard to avoid in interior spaces and when the viral mass injected into a turbulent room is sufficiently high, a person can be infected over one turbulent diffusion time. as per equation (2). More likely are spaces such as school classrooms or work places where occupants remain in the room for extended periods after which virions are fully mixed are only exposed to a critical load over longer times. For reflecting walls and no open windows, the time scale for infection is given by Eq. (6) and illustrated in the upper panel of Fig. 1. Spending more time than this in such an unventilated room would lead to infection, independent of physical distancing. To provide simple practical information and intuition the role of ventilation, we estimated the minimum size of an open window needed to mitigate infection for a one-time viral load for the case that would otherwise cause infection on the time scale of equation (6). The critical window size is given by Eq. (10) and depends on one unknown, the critical VLAD mass for infection Mc . The lower panel of Fig. 1 exemplifies the effect, showing that a window area ≥ 2m2 is enough to prevent occupants being infected for a viral load that could potentially infect 50 people. Open windows are therefore extremely helpful. These conclusions are consistent with, but complementary to empirical approaches such as those of of [31]. Masks reduce the area of the face and thus decrease the required open window area to obtain the same protection. In addition to windows, retrofitting HVAC systems with UVC (ultra violet c) or other anti-viral filters will certainly also reduce exposure, particularly from air passing between rooms in building complexes and depending on the efficiency of the filter, there is a correspondence in efficacy to a window area of a particular size. In recognize that the efficacy of ventilation can be thouht of as air passing through a window or filter of a specific area, note that interior walls themselves provide a useful large surface area. If they were imbued with anti-viral materials and constructed to mitigate boundary layer effects, rooms could potentially be safe after one turbulent diffusion time following any VLAD injection. Such measures are also effective for other 5
Figure 1: Critical time for infection and condition for safety. Upper panel: the time tc for a non-ventilated room with reflecting walls according to Eq. (6) plotted as a function of the number of people that can be infected by the injected load and room radius R. Infection is prevented by spending a time t ≤ tc . Lower panel: mitigation by introducing ventilation via an open window of area W . The planar surface is the critical surface tc /tsaf = 1 above which the time scale for infection tc exceeds the time scale tsaf for which enough virons have diffused through the window. Points on the curved surface above the plane indicate safe occupancy duration. The line of intersection between the surfaces shows the window area needed as a function of the viral mass injected into the room according to Eq. (10). airborne illnesses such as the flu and offer the economic benefits of reducing sick days, even non-pandemic circumstances. Precision is not required for our key results and message, but its further pursuit involves more complex numerical modeling. These can include turbulent eddy and droplet size spectra with size-dependent airborne survival times; different droplet-eddy coupling times as a function of droplet and eddy sizes and humidity; repeated, continuous, or time-dependent injection of VLADs from different room-locations and sources [28, 29]; temperature gradients [32]; room geometry, inhomogeneous turbulence; and wall boundary layer effects [30]. Acknowledgments EB acknowldges support from NSF Grant AST-1813298 and Department of Energy grants DE-SC0020432 and DE-SC0020434. GG acknowledges support from NSF Grant IIS-2029095. We thank W.J. Forrest and J.A. Tarduno for related discussions. None of these funding sources had any involvement or influence in any aspect of the study. Author Contributions Authors Blackman and Ghoshal jointly conceived the idea for the study. Author Blackman carried out the basic calculations and drafted the initial version of the manuscript based on joint discussions. Ghoshal further contributed to writing the manuscript, and the interpretation and elucidation of the results, Corresponding author Blackman confirms full access to everything in this study and had final responsibility for the decision to submit for publication. 6
Conflicts of Interest We have no conflicts of interest to report. References References [1] F. Wu, S. Zhao, B. Yu, Y.-M. Chen, W. Wang, Z.-G. Song, Y. Hu, Z.-W. Tao, J.-H. Tian, Y.-Y. Pei, M.-L. Yuan, Y.-L. Zhang, F.-H. Dai, Y. Liu, Q.-M. Wang, J.-J. Zheng, L. Xu, E. C. Holmes, Y.-Z. Zhang, A new coronavirus associated with human respiratory disease in china, Nature 579 (7798) (2020) 265–269. doi:10.1038/s41586-020-2008-3. [2] P. Zhou, X.-L. Yang, X.-G. Wang, B. Hu, L. Zhang, W. Zhang, H.-R. Si, Y. Zhu, B. Li, C.-L. Huang, H.-D. Chen, J. Chen, Y. Luo, H. Guo, R.-D. Jiang, M.-Q. Liu, Y. Chen, X.-R. Shen, X. Wang, X.-S. Zheng, K. Zhao, Q.-J. Chen, F. Deng, L.-L. Liu, B. Yan, F.-X. Zhan, Y.-Y. Wang, G.-F. Xiao, Z.-L. Shi, A pneumonia outbreak associated with a new coronavirus of probable bat origin, Nature 579 (7798) (2020) 270–273. doi:10.1038/s41586-020-2012-7. [3] G. W. H. Organization, Who coronavirus disease (covid-19) dashboard (2020). URL https://covid19.who.int [4] M. Hadei, P. K. Hopke, A. Jonidi, A. Shahsavani, A letter about the airborne transmission of sars-cov-2 based on the current evidence, Aerosol and Air Quality Research 20 (4) (2020) 1–4. [5] J. L. Santarpia, D. N. Rivera, V. L. Herrera, M. J. Morwitzer, H. M. Creager, G. W. Santarpia, K. K. Crown, D. M. Brett-Major, E. R. Schnaubelt, M. J. Broadhurst, J. V. Lawler, S. P. Reid, J. J. Lowe, Aerosol and surface transmission potential of sars-cov-2, medRxiv 10.1101/2020.03.23.20039446. [6] Q. Li, X. Guan, P. Wu, X. Wang, L. Zhou, Y. Tong, R. Ren, K. S. M. Leung, E. H. Y. Lau, J. Y. Wong, X. Xing, N. Xiang, Y. Wu, C. Li, Q. Chen, D. Li, T. Liu, J. Zhao, M. Liu, W. Tu, C. Chen, L. Jin, R. Yang, Q. Wang, S. Zhou, R. Wang, H. Liu, Y. Luo, Y. Liu, G. Shao, H. Li, Z. Tao, Y. Yang, Z. Deng, B. Liu, Z. Ma, Y. Zhang, G. Shi, T. T. Y. Lam, J. T. Wu, G. F. Gao, B. J. Cowling, B. Yang, G. M. Leung, Z. Feng, Early transmission dynamics in wuhan, china, of novel coronavirus– infected pneumonia, New England Journal of Medicine 382 (13) (2020) 1199–1207. doi:10.1056/NEJMoa2001316. [7] L. Bourouiba, Turbulent Gas Clouds and Respiratory Pathogen Emissions: Potential Implications for Reducing Trans- mission of COVID-19, JAMA 323 (18) (2020) 1837–1838. doi:10.1001/jama.2020.4756. [8] N. van Doremalen, T. Bushmaker, D. H. Morris, M. G. Holbrook, A. Gamble, B. N. Williamson, A. Tamin, J. L. Harcourt, N. J. Thornburg, S. I. Gerber, J. O. Lloyd-Smith, E. de Wit, V. J. Munster, Aerosol and surface stability of sars-cov-2 as compared with sars-cov-1, New England Journal of Medicine 382 (16) (2020) 1564–1567. arXiv:https: //doi.org/10.1056/NEJMc2004973, doi:10.1056/NEJMc2004973. URL https://doi.org/10.1056/NEJMc2004973 [9] G. K. Batchelor, The Theory of Homogeneous Turbulence, Cambridge University Press, 1953. [10] G. Csanady, G. Csanady, Turbulent Diffusion in the Environment, Geophysics and Astrophysics Monographs, Springer Netherlands, 1973. [11] M. R. Maxey, The gravitational settling of aerosol particles in homogeneous turbulence and random flow fields, Journal of Fluid Mechanics 174 (1987) 441–465. doi:10.1017/S0022112087000193. [12] W. D. McComb, The physics of fluid turbulence, Oxford University Press, 1990. [13] A. K. Blackadar, Turbulence and diffusion in the atmosphere, Springer Nature, 1997. [14] A. Monin, A. Yaglom, J. Lumley, Statistical Fluid Mechanics: Mechanics of Turbulence, no. v. 1 in Dover books on physics, Dover Publications, 2007. [15] A. Monin, A. Yaglom, J. Lumley, Statistical Fluid Mechanics: Mechanics of Turbulence, no. v. 2 in Dover books on physics, Dover Publications, 2007. [16] K. R. Sreenivasan, Turbulent mixing: A perspective, Procs. Nat. Acad. Sci. U.S.A. 116 (37) (2019) 18175–18183. doi: 10.1073/pnas.1800463115. [17] K. A. Prather, C. C. Wang, R. T. Schooley, Reducing transmission of sars-cov-2, Sciencedoi:10.1126/science.abc6197. [18] Y. Li, H. Qian, J. Hang, X. Chen, L. Hong, P. Liang, J. Li, S. Xiao, J. Wei, L. Liu, M. Kang, Evidence for probable aerosol transmission of sars-cov-2 in a poorly ventilated restaurant, medRxiv 10.1101/2020.04.16.20067728. [19] J. Lu, J. Gu, K. Li, C. Xu, W. Su, Z. Lai, D. Zhou, C. Yu, B. Xu, Z. Yang, Covid-19 outbreak associated with air conditioning in restaurant, guangzhou, china, 2020, Emerg Infect Dis. 26 (7). [20] S. Y. Park, Y.-M. Kim, S. Yi, S. Lee, B.-J. Na, C. B. Kim, J. il Kim, H. S. Kim, Y. B. Kim, Y. Park, I. S. Huh, H. K. Kim, H. J. Yoon, H. Jang, K. Kim, Y. Chang, I. Kim, H. Lee, J. Gwack, S. S. Kim, M. Kim, S. Kweon, Y. J. Choe, O. Park, Y. J. Park, E. K. Jeong, Coronavirus disease outbreak in call center, south korea, Emerg Infect Dis. 26 (8). [21] L. Hamner, P. Dubbel, I. Capron, A. Ross, A. Jordan, J. Lee, J. Lynn, A. Ball, S. Narwal, S. Russell, D. Patrick, H. Leibrand, High sars-cov-2 attack rate following exposure at a choir practice — skagit county, washington, march 2020, MMWR Morb Mortal Wkly Rep 2020 69 (2020) 606–610. [22] E. C. Riley, G. Murphy, R. L. Riley, Airborne spread of measles in a suburban elementary school, Am. J. Epidemiol. 107 (5) (1978) 421–432. [23] C. Béghein, Y. Jiang, Q. Y. Chen, Using large eddy simulation to study particle motions in a room, Indoor Air 15 (4) (2005) 281–290. [24] A. A. Aliabadi, S. N. Rogak, K. H. Bartlett, S. I. Green, Preventing airborne disease transmission: review of methods for ventilation design in health care facilities, Adv Prev Med 2011 (2011) 124064. 7
[25] T. Smieszek, G. Lazzari, M. Salath?, Assessing the Dynamics and Control of Droplet- and Aerosol-Transmitted Influenza Using an Indoor Positioning System, Sci Rep 9 (1) (2019) 2185. [26] K. K.-W. To, O. T.-Y. Tsang, W.-S. Leung, A. R. Tam, T.-C. Wu, D. C. Lung, C. C.-Y. Yip, J.-P. Cai, J. M.-C. Chan, T. S.-H. Chik, D. P.-L. Lau, C. Y.-C. Choi, L.-L. Chen, W.-M. Chan, K.-H. Chan, J. D. Ip, A. C.-K. Ng, R. W.-S. Poon, C.-T. Luo, V. C.-C. Cheng, J. F.-W. Chan, I. F.-N. Hung, Z. Chen, H. Chen, K.-Y. Yuen, Temporal profiles of viral load in posterior oropharyngeal saliva samples and serum antibody responses during infection by sars-cov-2: an observational cohort study, The Lancet Infectious Diseases 20 (5) (2020) 565–574. doi:10.1016/S1473-3099(20)30196-1. [27] [28] G. Buonanno, L. Stabile, L. Morawska, Estimation of airborne viral emission: quanta emission rate of sars-cov-2 for infection risk assessment, medRxivdoi:10.1101/2020.04.12.20062828. [29] V. Stadnytskyi, C. E. Bax, A. Bax, P. Anfinrud, The airborne lifetime of small speech droplets and their potential importance in sars-cov-2 transmission, Procs. Nat. Acad. Sci. U.S.A. 117 (22) (2020) 11875–11877. doi:10.1073/pnas. 2006874117. [30] S. Lejeune, T. Gilet, L. Bourouiba, Edge effect: Liquid sheet and droplets formed by drop impact close to an edge, Phys. Rev. Fluids 3 (2018) 083601. doi:10.1103/PhysRevFluids.3.083601. [31] G. A. Somsen, C. van Rijn, S. Kooij, R. A. Bem, D. Bonn, Small droplet aerosols in poorly ventilated spaces and sars-cov-2 transmission, The Lancet Respiratory Medicine 8 (7) (2020) 658 – 659. doi:https://doi.org/10.1016/S2213-2600(20) 30245-9. URL http://www.sciencedirect.com/science/article/pii/S2213260020302459 [32] T. Elperin, N. Kleeorin, I. Rogachevskii, Turbulent Thermal Diffusion of Small Inertial Particles, Phys. Rev. Lett. 76 (2) (1996) 224–227. doi:10.1103/PhysRevLett.76.224. 8
You can also read