Ventilation mode performance against airborne respiratory infections in small office spaces: limits and rational improvements for Covid-19

The Coronavirus disease (Covid-19), caused by the SARS-CoV-2 virus, has produced significant social and economic disruptions in different countries. Current evidence suggests a strong correlation between the infection and the cohabitation of indoor spaces. International organizations and experts consider that the airborne transmission through aerosols can occur in specific conditions and that inadequate ventilation increases the risk of infection. As a result, the increase in ventilation rates and air filtration efficiencies are recommended for public buildings in the context of Covid-19, with significant impacts on energy consumption, and a paradigm shift in the design of ventilation systems is necessary for this new context. Therefore, this study has assessed the comparative performance of the displacement ventilation and the mixed ventilation mode on reducing the risk of long-range airborne infection for the Covid-19 in a small office application. A coupled multizone-CFD (Computational Fluid Dynamics) software developed by the National Institute of Standards and Technology was used in this study to assess the relative performance of several design solutions related to different ventilation modes, filter efficiencies, and outdoor air flow rates. The results demonstrate that the displacement ventilation technique produces a better overall performance in reducing the SARS-CoV-2 airborne infection risk than the conventional mixed ventilation for all the studied cases.


Covid-19 pandemic context
The Coronavirus disease , caused by the SARS-CoV-2 virus, has spread throughout a wide geographic area [1,2], causing significant social and economic disruptions in different countries [1], as well as more than one million deaths to this date [2]. The restrictions on individual mobility, commercial activity, international travel, and school closures, related to the pandemics' control, produced an unprecedented impact in the US stock market [3] as well as in the world economy [4].

Covid-19 airborne transmission in indoor spaces
The World Health Organization (WHO) [5] and the US Centers for Disease Control and Prevention (CDC) [6] suggest that the virus is mainly transmitted from person to person through respiratory droplets produced by expiratory activities. Some studies have demonstrated that the virus may remain viable after aerosolization in the indoor environment [7][8][9]. Qian et al. [10] evaluated 328 Covid-19 outbreaks involving three or more infected persons in 120 cities in China. The results show a strong correlation between the infection and the cohabitation of indoor spaces [10]. The recent outbreak cases at the Skagit Valley choral [11] in the USA, the call center in South Korea [12], and the bus travel episode in China [13], with a large number of infected subjects, have cast a light on the role of super spreaders and airborne propagation. Morawska and Milton [14] produced a letter, supported by 239 scientists, appealing to the medical community and to the relevant bodies to recognize the potential for airborne spread of . The different vocabulary between the exposure science community and the medical infectious disease community regarding airborne transmission results in a communication difficulty that hampers collaborative efforts [15][16][17]. An arbitrary cutoff size of 5 μm is often used to separate the droplets' direct transmission and the airborne transmission [16,18,19]. Particles much larger than 5 µm may be transported through air currents for a significant distance (when compared to a typical indoor room scale) [15,16]. The cutoff size is affected by the droplet composition [18] and indoor airflow patterns [16][17][18], which introduces high uncertainty [18]. Small droplets (d ≤ 100 μm) could remain airborne for significant enough time (e.g., tens of seconds to several minutes) to be inhaled, depending on the indoor airflow patterns [17]. Specialized researchers on aerosol science and exposure science identify an urgent need to clarify the methodology to distinguish between aerosols and larger droplets (that leave the suspension in a low timescale due to gravity) using the threshold size of 100 µm, rather than the traditional size of 5 µm [16,20]. This characteristic size is more effective to separate their aerodynamic behavior, inhalation capability, and efficacy of intervention measures [16,20]. The American Society of Heating, Refrigeration and Air Conditioning Engineers (ASHRAE) suggests that the Transmission of SARS-CoV-2 through the air is sufficiently likely that airborne exposure to the virus should be controlled [21]. The WHO [5] and the CDC [6] consider that airborne transmission through aerosols can occur in specific conditions and that inadequate ventilation increases that risk.

Motivation, objectives, and justification
of the present study ASHRAE [22], the Federation of European Heating, Ventilation and Air-Conditioning Associations (REHVA) [23], the CDC [24], and the WHO [25] recommend measures related to the increase in ventilation rates and filtration efficiencies for public buildings in the context of Covid-19. MERV-13 filters, as rated by ANSI/ASHRAE 52.2 standard [26], are recommended by ASHRAE [22] and the CDC [24]. However, systems designed to meet the minimum requirements of a ventilation standard like ANSI/ASHRAE 62.1 standard [27] will have a MERV-8 filter [26], leading to a filter upgrade necessity in many cases. This filter upgrade may be hampered due to existing fan pressure constraints or may lead to an increase in energy consumption related to the higher filter pressure loss. Regarding ventilation, ASHRAE [22], REHVA [23], the CDC [24], and the WHO [25] recommend increasing outdoor airflow rate to the maximum allowable value (100% if possible). Providing 100% outdoor airflow rate to a commercial building produces an energy penalty that is seldom justified [28]. Therefore, the majority of the existing commercial buildings worldwide use the minimum affordable ventilation rate [28]. Increasing the ventilation rate in existing buildings may lead to indoor moisture control disruption in hot & humid climates and will increase energy consumption. A paradigm shift in the design of ventilation systems is necessary for this new context, with a focus on source control, advanced air distribution, and occupant-based microenvironment control [29]. ASHRAE suggests that research should address methods to quantify the relative airborne infection control performance and cost-effectiveness of specific engineering strategies [21]. Previous studies have demonstrated that the use of high values of ventilation rates has a limited efficiency on airborne contaminant control [30][31][32][33][34] and that an air distribution strategy that reduces the contaminant residence timescale at the room breathing zone provides better performance [32,[35][36][37]. Several studies have demonstrated that the displacement ventilation mode (DV) is a promising technique for enhancing the ventilation performance and reducing the risk of airborne infection in hospital settings [38] and commercial applications [39][40][41]. However, in all these studies, the contaminant concentration was modeled as null in the supply air, in a scenario that would resemble a 100% outdoor air ventilation system or a high-efficiency air filtration. Since recirculation and lower efficiency air filters are often applicable in commercial and office spaces [27], a complimentary evaluation of the DV ventilation performance with these features is necessary. Therefore, the objective of this study is to assess the comparative ventilation performance of the DV and the mixed ventilation (MV) system on reducing the risk of airborne infection for the Covid-19 in an office application for cases involving the use of air recirculation and different filter efficiencies. Given the current discussion by the scientific community regarding the appropriate particle size range for airborne transmission, this study has focused on controlling the fraction of the emitted respiratory droplets that may remain suspended in the air for a longstanding period, consistent with the long-range and long-time scales of airborne transmission in indoor settings. Exposure to the exhalation jet of an infector (short-range airborne or direct spray route) that occurs in close contact and leads to higher doses [18] was not considered in this study to focus on the ventilation role in limiting the longrange airborne transmission. This study is aligned with the ASHRAE proposition [21] and expands the evaluation of the DV technique on airborne infection control to commercial applications with air recirculation and different filtration efficiencies, contributing to a better understanding of the rational limits of ventilation control in the context of Covid-19 and other respiratory diseases.

General description of the method
Buonanno et al. [42,43] proposed a method to estimate the rate of generation of SARS-CoV-2 infectious quanta that can be useful for assessing the long-range airborne infection risk of different ventilation configurations. In this approach, the quantum is the dose of respiratory droplets' residues required to cause infection in 63% of the susceptible occupants of an enclosed space. The quanta emission rate is estimated by the method as a function of the viral load in the sputum, respiratory activity, and activity level of the infectious subject [42]. As many uncertainties persist at the moment regarding the new Coronavirus disease, the method presents many uncertainties and data extrapolation from the SARS-CoV-1 (also a coronavirus). Despite that, the authors sustain its relevance for allowing a quantitative assessment during the course of the pandemic, when limited data are available [42]. The method has been used to assess the efficiency of infection control measures related to ventilation and quarantine interventions [42,43] and to assess the quanta emission rate related to specific outbreak cases in retrospective studies [11,43]. Therefore, the method proposed by Buonanno et al. [42,43] was used in this study in a coupled multizone-CFD analysis to estimate the long-range airborne infection risk of a cohabitant in a model of a small office configuration with two occupants. The quanta emission rate, estimated by the method of Buonanno et al. [42,43], was used as a source input in a numerical model of the office's ventilation system to predict the resultant quanta concentration at the breathing zone of the cohabitant (n b ). CONTAM, a public domain software developed by the National Institute of Standards and Technology (NIST) [44], and its Computational Fluid Dynamics (CFD) module, developed by Srebric et al. [45], and improved by Wang et al. [46], were used in this study to model the office ventilation system and to assess the resultant quanta concentration of several design solutions related to ventilation modes, filter efficiencies, and outdoor air flow rates. As in the method of Buonanno et al. [42,43], the resultant quanta concentration at the breathing zone of the cohabitant was integrated over time using the Wells-Riley equation [47] to determine the infection risk (R) of each scenario: where IR is the inhalation rate [m 3 h −1 ] and is the exposure time. The infection risk of each occupant was assessed for each studied case, exchanging the position of the emitter and the susceptible. The following premises were adopted for the simulations: the mode of infection investigated was through the long-range airborne route only (direct and indirect contact risks were not investigated, to focus on ventilation performance); the simulation period was set as 8 h, related to a commercial shift of work; all occupants remain in the space during all the simulation period; inhalation rate was set as IR = 0.96 m 3 h −1 , consistent with an office activity [42]. "Appendix A" presents a flowchart that encompasses all the simulation steps. "Appendix B" presents the verification of the mathematical model of infection.

Description of the modeled office and its basic ventilation system
The configuration of the studied office is shown in Fig. 1 and is the same as the small office modeled by Chen et al. [39] for the ASHRAE research project RP-949. This setting was adopted given that: it resembles a typical configuration in the USA for small offices [39]; high-quality experimental data on airflow and contaminant distribution are available through the ASHRAE RP-949 benchmark [39] for CFD validation. The floor area is 18.8 m 2 with two occupants (#1 and #2), and the total heat dissipation from people, light fixtures, and equipment is 634 W.
The air-conditioning system was modeled for the base case with the characteristics presented in Table 1, according to the ASHRAE RP-949 benchmark features [39] and the ASHRAE 62.1 standard [27] procedure for ventilation calculations regarding the outdoor air flow rate for offices. Two different ventilation modes were investigated with the same airflow and thermal load rates: MV (mixed) and DV (displacement).

Estimative on quanta source strength for SARS-CoV-2
The approach proposed by Buonanno et al. [42,43] is based on the hypothesis that the emitted droplets have the same viral load as the sputum of the infectious subject (aerosol source).
The emitted viral load is therefore estimated by a mass balance, expressed in terms of quanta emission rate (ER q , quanta h −1 ), by the following equation [43]:   where c v is the viral load in the sputum (RNA copies mL −1 ), c i is the conversion factor (quanta RNA copies −1 ), V d is the droplets generation rate by the expiratory activity (mL h −1 ), and ∅ d is the droplet volume concentration at the exhaled flow (mL m −3 ).

Droplets generation rate
A wide range of droplets' sizes is generated by the expiratory activity, from the sub-micrometer order to the millimeter order [18,49]. Following their exhalation, dehydration due to liquid evaporation reduces the size of the droplets [18,49]. The final dehydrated equilibrium sizes depend on droplet composition and surrounding air relative humidity [18] and range from 20% [49,50] to 50% [51] of the initial sizes. Gravity acts on the deposition of the larger droplets, while the small droplets follow the airflow motion [18]. The threshold size between small and large droplets is understood to vary from 60 µm to 120 µm [18]. This study has focused on the fraction of the small droplets that may remain suspended in the air for a longstanding period, consistent with the time scales of indoor ventilation control. De Oliveira et al. [18] mention that recent studies [15,17] suggest that values around 20 µm should be used for ventilation calculations and that the approach of using a single value for all applications is inappropriate because the particle cutoff size is flow dependent. Barzant and Bush [49] provide a method to assess the 'critical particle size' for long-range airborne transmission in a given scenario, based on the well-mixed sedimentation and ventilation timescales. For the cases assessed in this study, with an effective ventilation rate ranging from λ c = 2.9 to λ c = 4.2, the critical particle diameter according to Barzant and Bush [49] ranges from 9 µm to 12 µm (for a 'well-mixed' scenario). The real airflow field in an indoor setting presents several flow structures that modify the particle transport predictions by the 'well-mixed' scenario, as jets, thermal plumes, recirculation, entrainment, etc. The transient simulations performed by Vuorinen et al. [17] incorporated these flow features through CFD and demonstrated that the transport pattern of massless particles and 10-µm particles by the indoor airflow was relatively similar, while the transport of 20-µm particles was noticeable different and relatively limited, when compared to the former. This limited transport is especially highlighted if the analysis is focusing on a long-timescale transport as in the case of this study. Therefore, the use of respiratory droplets' residues with a maximum diameter in the O(10 1 )μm is a reasonable assumption for the current analysis. Different studies have produced estimates on V d with considerable discrepancies, given the uncertainty and complexity of the associated experimental measurements. A range from O(1) to O(10) nL h −1 for speaking was used in different studies [18,52], based on past experimental data [53,54]. Mikszewski et al. [55] suggest the value of 4900 nL h −1 based on the recent measurements of Stadnytskyi et al. [50] for loudly speaking. Prentiss et al. [56] mention that these discrepancies may be related to the sensitivity of each experimental method for the measurement of particles in the submicron and supermicron ranges and suggest that the estimate of Stadnytskyi et al. [50] is more appropriate for their study. Buonanno et al. [43] have adopted a range from 1000 to 2500 nL h −1 for breathing and 1000 to 4500 nL h −1 for speaking, consistent with Stadnytskyi et al. [50].

Viral load and infectious dose
The viral load ranges over many orders of magnitude across patients [43,56], and for each patient, it changes during the disease progression [43,56]. There is a risk of pre-symptomatic transmission of SARS-CoV-2, as well as transmission from asymptomatic individuals [57], and the higher viral loads are expected on or just prior to the day of symptom onset [57]. A survey on studies [58][59][60][61][62][63][64] that reported the viral load in the sputum throughout the infection is presented in Table 2. The reported values from these studies present the trend of higher values in the early days of the infection, with posterior reduction on the course of the disease.
Although high viral loads have been measured in these studies, there is still little experimental evidence of the actual amount of viable virus in aerosol samples from indoor air. The failure on the detection of viable SARS-CoV-2 from aerosol samples is often reported [9]. Lednicky et al. [9] consider that usual air sample methods may contribute to this failure due to their harsh collection process that may inactivate the virus. Using a new air sampling technology (VIVAS), Lednicky et al. [9] were able to identify viable virions from aerosol samples. The amount of airborne virus detected was small [9], and the authors [9] address technical issues that may have contributed to this result, related to the collection complexity and the high air exchange rates of the sample location (hospital). Further research on the topic is necessary [9], given the complexity and the uncertainties related to this new disease. Some studies [18,42,43,49,56] have presented estimates on the number of necessary virions to cause an infection on a susceptible person (N 0 , RNA copies quanta −1 ), based on epidemiological data. The associated c i = 1∕N 0 from these studies [18,42,43,49,56] is assessed on the range from O(10 −1 ) to O(10 −3 ). Buonanno et al. propose the use of c i = 1.4 × 10 −3 (quanta RNA copie −1 ) in their most recent research [55], improving the predictive estimation approach reported in previous works [42,43].

Quanta generation rate
Buonanno et al. have provided estimates on the quanta emission rate using Eq. 2 and a Monte Carlo method to account for variations on input data [43,55]. Probability density functions were considered using normal [43] and lognormal [55] distributions for the viral loads [43,55], droplets emission rates [43], and infectious dose [43]. Table 3 presents the quanta generation rate estimated by these studies for two expiratory activities related to office work,.
The median values are related to the average viral load on the day of symptom onset reported in reference studies [58][59][60][61][62][63]. The average viral load produced a relatively low quanta emission rate spanning from O(10 −1 ) to O(10 1 ), regardless of the expiratory activity. An emission rate on a higher-order (10 2 ) may be obtained for the higher percentiles (≥ 90th) estimates of the probability density functions in the case of vocalization. These resultant values are consistent with the emission rates of SARS-CoV-2 analyzed by other studies for oral breathing and speaking in light activities. Barzant and Bush [49] have inferred a value of 14, 15, and 45 quanta h −1 for the cases of the Ningbo tour bus [13], the Diamond Princess cruise ship [65], and the Wuhan city outbreak [10], respectively. Therefore, this study selected two orders of emission rates for evaluation, representing the boundaries of a possible range in an office application: 1 and 100 quanta h −1 . The lower order was evaluated with the breathing aerosol size distribution, while the higher order was evaluated with the unmodulated vocalization aerosol size distribution, given the data in Table 3 and the former discussion.

The relative performance of higher filter efficiencies and ventilation scale-up on infection risk
Different filter efficiencies and outdoor air flow rates were also selected for evaluation with the MV and the DV modes to study the relative performance of different design arrangements on reducing the airborne infection risk of SARS-CoV-2 and to assess the quality of the numerical simulations. These different design arrangements were selected according to the recommendations of specialized organizations on filter and ventilation scale-up in the context of Covid-19 [22][23][24][25]. Therefore, the base cases presented in Table 1 were also evaluated with MERV-10 and MERV-13 filters (rated by ANSI/ASHRAE 52.2 standard [26]) and with 100% outdoor air exchange rate (no recirculation). It is expected that the infection risk that is predicted by the numerical simulations decreases with the filter efficiency upgrade and outdoor air flow rate improvement.

Coupled multizone-CFD simulation
A model of the office's ventilation system was constructed in the coupled multizone-CFD software (CONTAM), as depicted in Fig. 3.
In the coupled simulation, one of the zones was modeled in detail through the CFD technique (office), while the others use the "perfect-mixing" assumption (air-handling unit and ducts). In this method, the multizone module provides the boundary conditions to the CFD module. After solving the detailed flow field, the CFD module returns the updated variables to the multizone model. An iterative solution is performed, where information of airflow rates, concentrations, and pressures at the boundaries is exchanged between the two codes in each step solution until a convergence criterion is reached. The paper from Wang and Chen [66] provides detailed information on this coupled method. "Appendix C" presents the grid refinement studies and the validation of the CONTAM CFD code.

CFD airflow simulation
The detailed indoor airflow simulations were performed through the CFD module of CONTAM that uses a RANS (Reynolds Averaged Navier-Stokes) formulation with the Chen and Xu [67] zero-equation turbulence model. This turbulence model uses data from DNS (direct numerical simulation) for indoor airflow simulations and was selected since it has already proved to produce useful results in indoor airflow simulations [46,68]. The simulations were performed in non-isothermal steady-state conditions for 3-D incompressible flow (air). The heat output from occupants, lights, and equipment is the same as the benchmark by Chen et al. [39] for the small office case and was set in the CFD model accordingly. Given the coupled solution between the CONTAM and its CFD module, the numerical model of the office and its ventilation system, depicted in Fig. 3b, was set and calibrated to produce boundary conditions in the CFD zone that matches the operational conditions presented in Table 1 for the airflow rates and temperatures.
In this paper, equations' residuals were used to assess convergence. The solutions were considered converged when the scaled residuals were less than 10 −4 .

Contaminant transport simulation
This research aims to study the airborne transport of pathogen-laden droplets' residues, micrometer particles aerosolized from expiratory activities. For this class of particles, gravitational sedimentation and deposition may be neglected [69]. Some experimental studies [38,70] have demonstrated that a tracer gas method using SF 6 or N 2 O may be a suitable technique for studying the distribution of respiratory droplets' residual particles indoors. Olander et al. [71] consider that the tracer gas technique is suitable to describe the behavior of suspended particles with an aerodynamic diameter up to 5 μm. The CFD simulations performed by Vuorinen et al. [17] demonstrated that the transport pattern of massless particles and 10-µm particles by the indoor airflow was sufficiently similar to enable the passive scalar modeling for this class of particles. So, given the focus of this study in modeling the respiratory droplets' residues with a maximum diameter in the O(10 1 )μm (as discussed in Sect. 2.3.1), the air contaminant distribution was modeled using the transport of a passive scalar (quantum). This passive scalar was defined with the same properties of a tracer gas (SF 6 ) that has already been used to assess the transport of respiratory droplets' residues in indoor settings [38]. This simplified assumption was selected to include the "quanta" infection control metric in the coupled multizone-CFD simulations given that this study's purpose is to compare different scenarios' performance rather than precisely predict the exposure dose. Therefore, given the flow field characteristics and the quanta emission rate (source), the solution of the passive scalar transport equation predicts the resultant quanta concentration field. The CFD code uses the Reynolds' Analogy to account for the turbulent contaminant diffusivity through the turbulent viscosity and the turbulent Schmidt number (Sc t = 0.7).
In the tested scenarios, a quanta source was generated at the breathing region of the infectious subject. The source location was the same as the contaminant source of the ASHRAE RP-949 benchmark [39] for each infectious subject. The adopted CFD code does not provide resources to include a 'sink' term on the calculations, and therefore, the viral inactivation with time was not considered in this study. The relaxation time of droplets' residual particles is negligible [72], and therefore, the quanta source approach was consistent with a source injected with no relative velocity to the room airflow. The following method was adopted to account for the difference in filter efficiency regarding particle size and the passive scalar quanta approach: the total-quanta emission source was sub-divided into different passive scalars sources, and artificial filtration efficiency was numerically imposed for each passive scalar in the CON-TAM software, according to the representative respiratory droplets' residual size and air-handling filter rating. The method of Buonanno et al. [42] uses an approximation for the small droplets size distribution at the exhaled flow. This approximation is based on the size distribution of respiratory droplets' residues measured by Morawska et al. [53] for dehydrated particles up to 10 µm with four channels with midpoint diameters of D 1 = 0.8 µm, D 2 = 1.8 µm, D 3 = 3.5 µm, and D 4 = 5.5 µm. These midpoint diameters of each channel were used for the calculations in this study, as in the method of Buonanno et al. [42]. Morawska et al. [53] consider that the residence time of the respiratory droplets inside the probe used in their measurements is consistent with their estimated dehydration time, ensuring an equilibrium size for the particle distribution used in this study. Therefore, to account for the differences in filter efficiency regarding particle size, four passive scalars were included, with the same properties of a tracer gas (SF 6 ) but with different filter penetration efficiencies (related to 0.8-, 1.8-, 3.5-, and 5.5-µm droplets' residues size bins). The air-handling unit's filter efficiency was numerically imposed for each passive scalar in the multizone software, given the MERV rating and data by Kowalski et al. [73] related to the filter efficiency for each particle size. Table 4 presents the numerical filter efficiency for each passive scalar. The quanta emission rate in each size bin ER i is presented in Table 4 and was calculated based on the total emission rate ER q and the relative volume of aerosols, by Eq. 14 and data from Morawska et al. [53]: The resultant quanta concentration at the breathing zone of the cohabitant (n b ) is the sum of the quanta concentration of each size bin. The solutions were considered converged when the scaled residuals were less than 10 −6 for the passive scalar transport. Some differences in the exposure dose may arise for the cutoff at the 10 µm particle size. These differences are related to the possible exposure of larger droplets' residues that belong to the same size order but have slightly higher sizes (e.g., 20 µm). The approach used in this study to model the transport of the infectious quanta (as a passive scalar rather than the droplets particles) reduces these differences, as long as the quanta emission rate is consistent with the long-range airborne transmission.

Results and discussion
The results of the predicted dose and infection risk of occupants #1 and #2 with different ventilation modes, filter efficiencies, and outdoor air flow rates are depicted in Fig. 4 (occupant #1) and 5 (occupant #2). Figures 4a, c,  Fig. 5a, c present the results of the lower-order of emission rates (1 quanta h −1 ), while Figs. 4b, d, 5b, d present the results of the higher order (100 quanta h −1 ). The predicted doses and infection risks were lower for all the DV cases when compared to the same MV case. These results show a trend of reduction in the dose and infection risk according to the increase in filter efficiency and outdoor air flow rate, suggesting, in addition to the validation tests ("Appendix C"), consistent numerical simulation results. The results of Figs. 4 and 5 also imply that the infection risk is highly sensitive to the quanta generation rate, evidencing that a source control measure (e.g., the use of masks) may be a suitable and rational strategy for reducing the infection risk in a pandemic scenario. The higher-risk scenarios of this study were associated with an intensive vocalization expiratory activity, reinforcing that a social convention on the 'vocalization etiquette' would also be a rational strategy in indoor public spaces during the pandemic. The results also demonstrate that relatively low  Table 3), the base air-conditioning and ventilation configuration according to the ASHRAE standards may be able to produce a low airborne infection risk of a large proportion of the predicted scenarios involving infectious subjects. Improving the filter efficiency produced little benefits in reducing the infection risk of these cases of lower quanta emission rates. The upgrade on filter class from MERV-8 to MERV-13 resulted in a net increase in the clean air delivery rate (CADR) of 42.6% (from 37 L s −1 to 53 L s −1 ) for these cases. The resultant dose for all the cases is much lower than one quantum (Figs. 4a and 5a). The combined effect of these features in the exponential dose-response model was a low variation on the infection risk with the upgrade on filter efficiency (Figs. 4c and 5c).
However, the probability of a high-viral load subject or a 'superspreader' that produces a high emission rate (≥ 100 quanta h −1 ) may not be neglected [42,43]. For these cases, the accumulated dose for the MV mode is considerably higher than one quantum (Figs. 4b and 5b), reaching a plateau level where the infection risk is near 100%, regardless of the filter class (Figs. 4d and 5d). Increasing the filter efficiency presents limited effectiveness in reducing the infection risk of these cases due to the exponential characteristics of the dose-response infection model. Therefore, the upgrade to a higher-efficiency filter (MERV-13) provides little reduction in the long-range airborne infection risk with mixed ventilation (Figs. 4d and 5d) for the studied cases. The use of a displacement ventilation mode is more effective than the upgrade on filter efficiency and outdoor air flow rate in mixed ventilation mode for reducing that risk. The high ventilation effectiveness performance of the DV mode produced lower accumulated doses than the MV mode for all the cases. In the higher emission scenarios of 100 quanta h −1 , the dose was sufficiently low with the DV mode  Figs. 4b and 5b) to enable a higher sensitivity of reduction on infection risk due to filter upgrade by the dose-response model (Figs. 4d and 5d).
Regarding the relative position between infector and susceptible, the accumulated dose and infection risk of occupant #2 was higher than for occupant #1 in all the studied cases (Figs. 4 and 5). Occupant #2 is located in a downstream position relative to occupant #1 and the supply airflow. The associated contaminant transport and distribution were affected by these relative positions. These results demonstrate that the airflow distribution provided by the MV mode did not comply with the 'well-mixed' assumption, emphasizing the benefits of a detailed airflow calculation through the CFD technique for improving the analysis of ventilation design solutions.
Increasing the outdoor air flow rate (from 32 to 100% air exchange rate) yields approximately the same performance as the use of MERV-13 filters for the studied cases with air recirculation. It is expected that the use of HEPA filters produces a similar result to the 100% outdoor air exchange rate case, given the high HEPA filter efficiency (> 99.97% for 0.3 µm particles). The increase in outdoor air flow rates may disturb the moisture control in hot & humid climates and will increase energy consumption. On the other hand, the filter upgrade may be accomplished in several existing buildings, with strategies like the use of a lower air facevelocity at the filter (resulting in lower pressure drop). This strategy must be analyzed in a 'case-by-case' method, geared to match the existing fan available pressure to the requirements of the new filter pressure loss. The filter upgrade will increase the energy consumption, but improving the outdoor air will also produce an energy penalty related to the precooling or pre-heating necessity. A life-cycle cost analysis for each possible solution would provide the necessary consistency for selection. The simulations' results suggest that the displacement ventilation technique produces a better The streamlines included in Fig. 6 demonstrate that the DV mode reduces the residence timescale of a contaminant generated by the occupants when compared to the MV mode. The typical residence time for contaminants was assessed as 4 min for the DV mode and 13 min for the MV mode. These residence timescales are considerably lower than the virus inactivation half-life timescale (order of 1 h [7]), which supports the validity of the results despite not including the virus inactivation process due to software constraints. The buoyant transport provided by the thermal plumes in DV mode (Fig. 6a) and the location of the exhaust grille at the ceiling reduce the contaminants' residence time in the DV mode when compared to the MV mode. Therefore, the better overall performance on reducing the SARS-CoV-2 airborne infection risk of the DV mode in the studied cases is a result of its superior room air scavenging capability when compared with the MV mode. Figure 7 depicts the iso-surfaces of quanta concentration provided by the MV and DV mode for one of the studied cases.
The air distribution through a high-momentum jet above the occupied zone is characteristic of MV systems, and the related effects of entrainment, recirculation, and turbulence in air mixing result in the spread of the contaminant source throughout a larger proportion of the room (Fig. 7c and d) when compared to the DV mode ( Fig. 7a and b). The use of the displacement ventilation mode produces the contaminant's buoyant transport to the ceiling region and out of the breathing zone ( Fig. 7a and b), increasing the ventilation effectiveness, and it is a characteristic feature of the displacement ventilation mode [74]. This technique has already been associated with the enhancement of thermal comfort [75] and energy use reduction [76]. It may be a relevant solution in this new post-pandemic era, where the design of new buildings demands features that address a better infection control solution in future pandemic or outbreak scenarios. Specific room configuration might impact the results, and a CFD analysis should be performed for each specific configuration. This study contributes to providing a basic framework that can be applied for these future analyses. An experimental study on the airborne transport of a real infectious agent by the room ventilation system presents some difficulties and limitations regarding safety concerns. The use of a numerical simulation overcomes these limitations, and this study demonstrates that this technique may provide valuable data for studying the relative performance of different ventilation solutions in the future challenges for this new post-pandemic scenario. In this research, the airborne infection risk was evaluated using the method of Buonanno et al. [42,43] and the results of the contaminant concentration provided by the coupled multizone-CFD calculations. This approach links the capabilities of these two methods and expands their potential use for future analysis on ventilation performance and infection control.

Conclusions
This study has used a coupled multizone-CFD code and an infection risk model to assess the performance of different design arrangements on reducing the SARS-CoV-2 longrange airborne infection risk in a small office. The code was validated by the use of a high-quality benchmark. Comparing the results of the predicted dose and the infection risk of different solutions led to the following conclusions: (a) Relatively low infection risk was predicted for all the cases with a viral load of the infectious subject within the average values reported in specialized studies. This implies that a base air-conditioning and ventilation configuration according to the ASHRAE standards may be able to produce a low long-range airborne infection risk of a large proportion of the predicted scenarios involving infectious subjects in the cases assessed in this study (droplets' residual size up to 10 µm).
(b) Relatively high infection risk was predicted for all the cases with vocalization and a viral load of the infectious subject within the maximum values reported in specialized studies. The increase in outdoor air flow rates or filter efficiency in MV has little impact on reducing the risk of long-range airborne infection for these cases in this study (droplets' residual size up to 10 µm). Additional infection control measures are necessary for this scenario. (c) The displacement ventilation technique produces a better overall performance on reducing the SARS-CoV-2 long-range airborne infection risk than the conventional mixed ventilation for all the studied cases in this study (droplets' residual size up to 10 µm). The benefits of upgrading filter classes are highly dependent on the source term (quanta emission rate) and the method for modeling the infection risk. For the dose-response method used in this study, the use of a DV mode in the higher-emission cases produced a sufficiently low This study has some limitations. Other factors affect the SARS-CoV-2 infection risk mitigation, as the hands' hygiene, social distancing, the use of masks, and surface disinfection, among others. The main objective of this study was not to present a solution that mitigates the infection risk in an office but rather to promote a comparative study in the airborne infection control performance of different ventilation arrangements. The infection risk method presents many uncertainties, as many persist at the moment regarding the new Coronavirus disease. It uses data extrapolation from the SARS-CoV-1 (also a coronavirus), and high uncertainty is related to the infectious dose. The estimative on infectious droplets emission relies on a mass conservation criterion, and therefore, the droplets' viral load is proportional to the droplets' mass concentration. There are little data, at the time of this study, regarding the SARS-CoV-2 viral load distribution related to the droplets size distribution of commonly expiratory activities. However, there are available data on the relationship between viral load and droplets mass concentration for other respiratory viruses [77] that present this trend. The method uses experimental data of Morawska et al. [53] for an approximation regarding the size distribution of expiratory droplets' residues. Morawska et al. [53] consider that a sufficient residence time inside the probe used in their measurements is consistent with an equilibrium size [53]. However, Stadnytskyi et al. [50] mention that droplet dehydration might have been incomplete due to high humidity in the experiment of Morawska et al. [53]. These experimental studies are inherently complex, and the result of this former possibility is that the droplets' residues might reach smaller sizes than assessed in this study. However, these differences might not be substantial, given that the measurements of Morawska et al. [53] include some size reduction and that the evaporation time for the small droplets (the focus of this study) is considerably low. Some uncertainties arise from the method's approach using the 'quanta' metric to study the transport of infectious droplets residues and the simplifications related to the 'passive scalar' model. The buoyant flow acceleration by thermal plumes, typical for DV systems, may improve the transport of droplets with other sizes that were not assessed in this study (e.g., > 10 µm). The fraction of these droplets that follow the airflow path will be transported to the unoccupied zone and scavenged by the exhaust. However, another fraction might become airborne for a sufficient time to be inhaled by the co-occupant. In this case, depending on the emitter-receptor transport timescale, the infection route may not fulfill the long-timescale airborne route criterion selected in this present study. It is a complex phenomenon, and the 'passive scalar' approach limits an appropriate analysis. Future work could include a particle-tracking approach to study this important matter. The viral inactivation with time was modeled in the method of Buonanno et al. [42] as an additional removal rate of λ vi = 0.63 h −1 based on the virus airborne half-life assessment by van Doremalen et al. [7]. This approach may be suitable for 'well-mixed' calculations, but the extension to this study's CFD detailed calculations may not apply, and therefore, the viral inactivation with time was not considered. However, the resultant contaminant residence timescales were considerably lower than the virus inactivation half-life timescale (order of 1 h [7]), which supports the validity of the results despite not including the virus inactivation process.
Despite these limitations, the study provides relevant information on the relative performance of different design solutions, allowing a quantitative assessment during the pandemic of a new disease, when limited data are available. Experimental tests on the performance of air-conditioning systems on infection control using infectious agents have proven to be difficult, and more research is necessary on this topic. On the other hand, the use of numerical simulations may be used for parametric tests and to study the relative performance of different solutions, and the framework presented in this study may contribute to this task. Future progress on the research of this new disease may demand the revision of some parameters of the model. The basic framework of this study is a practical application of a numerical tool already established that can be easily adapted to update these parameters. This study has focused on small office applications, and a specific CFD analysis using the framework of this study may be used to assess specific configurations. Future work will assess the performance of other applications. Three sets of grids were tested, with representative element sizes of 0.2 m (coarse), 0.1 m (intermediate), and 0.05 m (fine), and progressive refinement. The simulation results for the SF 6 concentration were compared with the ASHRAE RP-949 [39] benchmark's experimental data for the measurements at the nine poles. The simulations used structured hexagonal grid meshes, manually built via the software graphical interface. The results of the grid refinement study are depicted in Fig. 10. The finer grid, with a representative element size of 0.05 m, was chosen due to little improvements from the intermediate mesh and an overall better agreement with the experimental data. This grid size (102 × 74 × 48) is close to the grid independent size (123 × 86 × 54) used by Wang and Zhai [78] in their simulations for the same case.
The same benchmark was also simulated with another CFD code (Fluent) that has been widely used in indoor airflow simulations [34,79,80] to provide a comparison parameter and a complimentary assessment on the quality of the results obtained by the CONTAM CFD code. The Fluent code was used with the RNG-κε turbulence model, given the superior performance among other eddy viscosity models for flow prediction in similar problems [80]. A non-structured mesh comprising tetrahedral elements with a representative element size of 0.05 m was built and used for the Fluent simulation. The results were also compared with the simulations provided by Chen et al. [39] for the same benchmark. Figure 11 depicts the results of dimensionless velocity and SF 6 profiles predictions at selected poles and a comparison with the experimental data from the ASHRAE RP 949 benchmark and the simulations provided by Chen et al. [39] for the same case.
The CFD simulations performed in this study for the indoor flow field prediction involve inherent complexity, given the asymmetric three-dimensional geometry and the specific characteristics of a non-isothermal transitional indoor airflow [81,82]. The Archimedes Number (Ar, relative magnitude of buoyancy forces to the inertial forces) for the simulated cases is Ar = 70, which increases the simulation complexity and reduces the accuracy of the simulated results. Some differences between the experimental and predicted data arise for the validation case (Fig. 11), due to: modeling simplifications; boundary conditions uncertainties; heat transfer and turbulence model limitations; discretization errors and numerical approximations. However, The CON-TAM CFD simulation of the validation case reproduces the vertical contaminant stratification trend measured by Chen Fig. 10 Results of the grid refinement tests a Plan view of the small office with the location of the measurement poles, heat sources and SF 6 source, according to the benchmark of Chen et al. [39]. Simulation results of SF 6 concentration for three sets of grids at poles 3 (b), 4 (c), 7 (d), and 8 (e) et al. [39] (Fig. 11), which is a characteristic feature of the displacement ventilation mode. It also produced a reasonable prediction on contaminant concentration at the breathing height and a similar performance when compared to the Fluent results and the simulations performed by Chen et al. [39] for the same case (Fig. 11). Therefore, the CONTAM CFD modeling was considered to provide acceptable predictions on ventilation performance for the office application, despite the inherent modeling complexity. Fig. 11 Results of the validation test with the CONTAM and Fluent codes and comparison with the simulations provided by Chen et al. [39] for the same benchmark (a) Plan view of the small office with the location of the measurement poles, heat sources and SF 6 source. Simulation results of dimensionless velocity for poles 2 (b) and 3 (c). Simulation results of dimensionless SF 6 concentration for poles 2 (d), 3(e), 4 (f), 6 (g), and 8 (h)