Energy Analysis of Control Measures for Reducing Aerosol Transmission of COVID ‐ 19 in the Tourism Sector of the “Costa Daurada” Spain

: In this paper, the use of HVAC systems and non ‐ HVAC control measures to reduce virus ‐ laden bioaerosol exposure in a highly occupied indoor space is investigated. A simulation tool was used to model the fate and transport of bioaerosols in an indoor space in the hotel industry (bar or pub) with three types of HVAC system (central air handling system (CAHS), dedicated outdoor air system (DOAS), and wall unit system (WUS)). Non ‐ HVAC control measures such as portable air cleaners (PAC) and local exhaust fans were considered. Occupant exposure was evaluated for 1 μ m bioaerosols, which transport SARS ‐ CoV ‐ 2, for 3 h/day of continuous source and exposure. The combined effects of ventilation (400 l/s of outdoor air), recirculated air filtration (90% efficacy), and a PAC with a capacity up to 900 m 3 /h mitigated the (normalized) integrated exposure of the occupant by 0.66 to 0.51 (CAHS) and 0.43 to 0.36 (DOAS). In the case of WUS, the normalized integrated exposure was reduced by up to 0.2 when the PAC with a capacity of up to 900 m 3 /h was used. The corresponding electricity consumed increased by 297.4 kWh/year (CAHS) and 482.7 kWh/year (DOAS), while for the WUS it increased by 197.1 kWh/year.


Introduction
The tourism sector accounted for about 12.4% of the gross domestic product (GDP) in 2019 [1] and is the linchpin of the Spanish economy, but because of the COVID-19 pandemic, it is estimated to have dropped by around 44% in 2020 [2]. Spain is also one of the leading destinations for international tourism, along with Italy, France, and the United States. Therefore, Spain is one of the few countries in the world suffering the consequences of this dual situation. The hospitality sector is a key part of international and domestic tourism in Spain and has undeniably been one of the hardest hit during the pandemic, with a dramatic fall in occupancy levels (e.g., overnight stays in hotels dropped by 95.1% in June 2020 in contrast to the same month in 2019 [3]). In general, the tourism sector is highly sensitive and vulnerable to risk situations caused by external factors such as the current pandemic. The vulnerability of the sector to COVID-19 has been quantified in the 50 provinces of Spain in a recent study [4].
The principal mode by people with COVID-19 is exposure to respiratory droplets carrying the virus (SARS-CoV-2). Respiratory droplets are produced by human exhalation (e.g., breathing, speaking, singing, coughing, sneezing) and span a wide range of sizes that can be divided into two main groups depending on how long they remain suspended in air [5][6][7][8][9][10]: larger droplets (typically >100 μm in size), some of which are visible, fall out of the air within seconds or minutes, while still close to the source; smaller droplets (considered 100 μm under stagnant ambient conditions [10]) and particles (formed when small droplets dry rapidly in the airstream) can remain suspended for several minutes or even hours and travel far from the source on air currents. However, the precise droplet size limit is dependent on the ambient air conditions and airflow patterns (e.g., still air or turbulent air) [10]. In the case of large droplets (e.g., 200 μm), the gravitational force (gravity settling) mainly determines the particle movement in air, while ambient air condition has a reduced effect in contrast to the small droplets. Therefore, the drying rate is usually insufficient for the large droplets to settle to the ground or on other nearby surrounding surfaces. Once respiratory droplets have been exhaled and as they move away from the source, their concentration decreases as they fall from the air (largest droplets first, and then the smallest), and as the remaining smaller droplets and particles become diluted in the volume of air they encounter.
Pathogens (e.g., SARS-CoV-2) expelled by the respiratory processes of infected people have three potential modes of transmission [11,12]: (i) infection by suspended particles (bioaerosols) via inhaling airborne respiratory droplets, (ii) droplets that travel directly to the mucosal surfaces of the susceptible host from the infected individual, and (iii) direct contact of respiratory droplets with contaminated surfaces (fomites). The importance of viral transmission via small airborne microdroplets and droplet nuclei (known as 'aerosols') has been widely discussed in the context of the COVID-19 pandemic [7,9,[13][14][15]. Airborne transmission, therefore, is one of three commonly accepted modes of viral transmission. Figure 1 illustrates the typical trajectories of particles in air for sizes between 0.01 μm and 100 μm. In addition to particle size, the behavior of these trajectories depends on factors such as initial velocity and flow field velocity. In Figure 1, the gravitational force ( ) and aerodynamic forces (lift force, , and drag/air resistance force, ) that act on the particle in air, which dictate its movement, are shown with particle diameter and corresponding dominant type of transmission mode. Airborne transmission spreads infection by exposing people to virus-containing respiratory droplets consisting of smaller droplets and particles that can remain suspended in the air over long distances (usually >1-2 m, which is recommended as social distancing) and time (typically hours) [18]. Several outbreaks of COVID-19 infection are thought to have been caused by airborne transmission, most of which have been in areas with poor ventilation [19,20]. Moreover, recent research findings also strongly support that SARS-CoV-2 can circulate through some heating, ventilation, and air conditioning (HVAC) systems [21,22].
According to the US Center for Disease Control and Prevention (CDC), the main infection control strategies in buildings or enclosed spaces can be classified as pathogen elimination, engineering and administration control, and personal protection [23]. Control strategies that can disinfect/purify the air of pathogens while ensuring the building's high energy performance have become crucial in recent years. The main strategies to mitigate and/or eliminate possible pathogen transmission are:  For airborne transmission: reduce the concentration of aerosolized pathogens (e.g., SARS-CoV-2) by diluting them with fresh air provided by the ventilation system, contaminated air filtration, and disinfection/sterilization of contaminated air. In addition, reduce the risk of inhalation by using masks and visors.  In case of pathogens spread by droplets: physical/social distancing, limits on the gathering of people, and restrictions on movements using quarantining measures.  For propagation of pathogens due to contact: handwashing, regular cleaning, and decontamination of surfaces that are hotspots for the spread of pathogens.
HVAC systems are often implemented in buildings and enclosed spaces to maintain a healthy and comfortable indoor environment. The type of HVAC systems and how they are operated can have a major impact on the transmission of respiratory infectious diseases by aerosols from an infected subject to uninfected occupants sharing the same indoor space in buildings [24][25][26]. Therefore, decreasing the exposure of uninfected occupants to virus-laden bioaerosols by using engineering control mechanisms is a key step in curtailing the spread of infectious diseases such as COVID-19, which is caused by SARS-CoV-2. Therefore, the type of HVAC system and how it is operated are of great importance for the quality of the building's indoor air, reducing the risk of infection and improving occupant wellbeing during the current pandemic and beyond.
Several methods can be used to control the spread of pathogen-laden aerosols inside buildings, for example, dilution, particulate filtration, and sterilization/disinfection [8,24]. Indoor air can be diluted using natural or mechanical ventilation and an in-room air cleaner. The recirculated air from the HVAC system is also filtered and/or sterilized using HVAC filters and/or ultraviolet germicidal irradiation (UVGI), respectively, before being supplied to the occupied space of the building. However, it is often recommended that air recirculation in a building's HVAC system be avoided to prevent the spread of pathogens (e.g., virus) inside [8]. Stand-alone portable air cleaning and disinfection devices can also be used to reduce (or eliminate) viruses from indoor space.
Diluting air in which pathogens such as SARS-CoV-2 might exist is seen as one of the more practical solutions for mitigating the risk of airborne transmission in indoor spaces. This can be done by applying the guidelines and recommendations of organizations and societies like the Federation of European Heating, Ventilation, and Air Conditioning Associations (REHVA) [23], ASHRAE [27], and others [28]. The main ventilation strategies recommended by REHVA and ASHRAE to mitigate the spread of COVID-19 in buildings are to increase the outside air fraction of the supply air, to open windows for natural ventilation, and to extend the operational hours of the HVAC system. However, present ventilation standards like ASHRAE 62.1-2019 do not consider infection control as their aim. Thus, the minimum ventilation rate for preventing airborne transmission of the virus is unknown [19]. Thereby, HVAC systems and other non-HVAC control measures should be further investigated so that air quality standards can be improved by diluting, disinfecting, and cleaning air with current technologies. This paper studies the relative reduction in virus-laden bioaerosol exposure in a highly occupied indoor space typical of the Spanish hospitality sector using (i) different types of HVAC system and (ii) non-HVAC control measures. Then, the energy consumption was estimated for each control measure considered based on Tarragona (Spain) weather conditions. This paper is organized as follows: first, the HVAC systems are described (Section 2), and details of the modeling methodology are presented (Section 3); then, Section 4 presents and discusses the results of the simulation and, finally, the paper concludes with a summary of the key findings.

Description of HVAC Systems
This study focuses on a highly occupied indoor space, typical of the Spanish hospitality industry (e.g., bar and pub) with different HVAC system types. The HVAC systems considered are a central air handling system (CAHS), a dedicated outdoor air system (DOAS), and a wall unit system (WUS), which are applied to serve the occupied space as illustrated in Figure 2. The CAHS shown in Figure 2a can be placed on the rooftop or adjacent to the building. It is a constant air volume (CAV) system that supplies the space with a mix of outdoor air and recirculated air at the required supply airflow rate, while a portion of return air from indoors is exhausted out of the building. The DOAS depicted in Figure 2b uses a separate outdoor air unit (OAU) which supplies the required outdoor air directly to the indoor space after conditioning it to the indoor set point temperature. The air within the space is recirculated, filtered, and conditioned using a room-level air-conditioning system (CS in Figure 2b). The WUS considered in the study, illustrated in Figure 2c, is a through-the-wall room air conditioner, which does not supply outdoor air for ventilation or filter the return room air.  In the case of the CAHS shown in Figure 2a, a constant air volume (CAV) HVAC configuration is used to deliver the cooling and heating requirements of the conditioned space with a water-cooled electrical chiller and natural gas-fired boiler. The CHW loop design setpoint temperature and ∆ are 8 °C and 6 °C, respectively, while the boiler HW loop design setpoint temperature and ∆ are 50 °C and 10 °C. The cooling water inlet temperature is set at 29.4 °C for the electrical chiller's condenser that uses an evaporative cooling tower. The cooling coil and heating coil (which use chilled water and hot water, respectively) setpoint temperatures are 16 °C and 22 °C. In addition, the zone cooling and heating supply air temperatures are set to 16 °C and 35 °C, respectively. An air-to-air heat-recovery-type plate heat exchanger with sensible and latent heat recovery was considered. It should be noted that this simulation did not consider economizer, humidifier, dehumidifier, and coil preheat.
The DOAS, Figure 2b, consists of an OAU with an air-to-air heat recovery heat exchanger and a thermal conditioning system (CS) which recirculates, conditions, and filters the air in the indoor space. The DOAS is based on a four-pipe fan coil unit (FCU) connected to the chilled water and hot water loops of a central chiller plant and natural gas boiler (similar to the above CAHS). The supply fan total efficiency and motor efficiency were set to 70% and 90%. The was set to 400 l/s on the basis of full occupancy and RITE's standard (provided in Table 1). The performance of the electrical chiller and gas boiler was modelled using an average COP of 4.8 and a boiler efficiency of 92%, respectively. The FCU cooling and heating coils' design setpoint temperatures were 24 °C and 20 °C, respectively. The recovery heat exchanger was modelled using a sensible and latent heat exchanger effectiveness of 70% and 65%, respectively.
The WUS, Figure 2c, was assumed to be an in-room air conditioner using a singlespeed direct expansion (DX) cooling/heating coil with design cooling supply air ∆T of 8 °C and a heating design supply air temperature of 50 °C. A typical cooling COP of about 3.0 was considered for the DX system (coil) in this study. There is no mechanical OA supplied to the space, so the corresponding cooling (or heating) demand was not estimated.
As an example of a highly occupied indoor space, a representative bar/pub in Tarragona (Spain) was selected. It was served by three types of HVAC system, and its dimensions are given in Table 1. The bar/pub is part of a bigger building which is, at the same time, a hotel. Figure 3 illustrations a 3D view of the supposed representative singlezone room considered in the simulation. The main façade faces north. The bar/pub also has a door (2 m × 2 m facing south on the front wall) and two windows facing east on the left side of the building. As the hotel is a conditioned space at similar temperature as the bar/pub, it is also assumed that the floor, ceiling, and right and back walls of the room are adiabatic. The thermal transmittance, U, of the components of the building envelope are set to the limit values of thermal transmittance allowed for the January climate zone of Tarragona (Spain) according to the Technical Building Code of Spain. The details of U values and other input parameters are listed in Table 1.   [38]; e maximum power to be installed per lighted surface according to use and average maintained horizontal illuminance (Em) based on the Spanish regulation; 1 (Em ≤ 600 lux); 2 (Em > 600 lux); f VEEI = value for the energy efficiency of the installation; g for category IDA 3 from RITE [36].

Modelling of Exposure to Virus-Laden Bioaerosols in an Indoor Space
In order to predict the transport and fate of an aerosolized virus within an indoor space, a simulation tool newly developed by the National Institute of Standards and Technology (NIST) of the United States, known as FaTIMA (Fate and Transport of Indoor Microbiological Aerosols), was used in this study [39][40][41]. This tool is based on a dynamic model of a single-zone room served by a mechanical ventilation system and uses mass conservation in the room to relate the aerosol sources. It also includes various removal mechanisms such as recirculated air filters, an in-room air cleaner, and a local exhaust fan besides removal due to natural processes like deposition (on surfaces) and decay. The aerosol concentration throughout the zone is uniform (i.e., a well-mixed approximation) [42,43]. However, more detailed methods can be used to study the influence of spatial variations on aerosols by using CFD simulations (e.g., in Refs. [44,45]). In this case, the governing mass balance for a single-zone volume of a room becomes [40]: where is the aerosol concentration in air (kg/m 3 ), is the aerosol penetration factor (0 ≤ ≤1) from outside air into the room, is the aerosol generation rate within the zone (kg/s), , is the total aerosol deposition on surfaces (kg/s), V is zone volume (m 3 ), is the volumetric airflow rate (m 3 /s), t is time (s), and is the in-room air cleaner (i.e., PAC) aerosol filtration efficiency (-).
The total surface area of deposited aerosols is estimated using: where , is the deposition area for surface i (m 2 ), is the number of surfaces in the zone (floor, walls, ceiling, and others, if any), and , is the aerosol deposition velocity for surface i (m/s). (Subscripts: ac = air cleaner in room, d = deposition, inf = infiltration, lx = local exhaust, oa = outdoor air, ret = return, sup = supply, tot = total and, z = zone.) The surface loading for each surface, , (kg/m 2 ), is related to the deposition rate using Equations (3) and (4).
where is the aerosol deposition rate (1/s). According to Equation (1), the rate of change in mass of aerosols within the zone's air volume is related to the aerosol source (i.e., outdoor air by infiltration , supply airflow , and aerosol generation within the zone ) and the removal rate from the zone (the return airflow , the local exhaust airflow , the air cleaner in the room , and the deposition on surfaces). The initial virus-laden aerosol concentration in both outside and indoor air is considered to be zero (i.e., 0.0) as only the effect of a contagious occupant sharing indoor space with other uninfected occupants is analyzed in this study.
The simulation tool model is based on CONTAM, which is a multi-zone indoor air quality and ventilation system computer program developed by NIST [46]. CONTAM computes on mass-based units [39,40] and inputs aerosol sources on a number basis (#), instead of a mass basis, since it converts mass into number-based variables using the specified aerosol diameter and density. Therefore, we analyze aerosols of a single size with an average aerodynamic diameter of 1 μm and a density of 1000 kg/m 3 (approximated by using water density because of the aqueous properties of exhaled respiratory aerosols). The influence of relative humidity and temperature on airborne aerosol transmission is not taken into consideration. All airflow rate calculations are performed on a mass basis since the effect of temperature and pressure on the density of air is neglected and takes the density of air at 20 °C and 101.325 kPa (i.e., 1.2041 kg/m 3 ).
The simulation outputs are obtained in terms of the time history of the aerosol concentration, the surface loading, and the integrated exposure that an occupant experiences in the zone [39,40]. The simulation period is set to be carried out for 24 h with a time step of 15 s. The integrated exposure, IE (kg•s/m 3 ) given by Equation (5), is the integral of the aerosol concentration in air of the space/room in which the occupant is exposed during the exposure duration, indicated by the specific start and end time (i.e., and , respectively) [40].
However, the occupant inhalation or body uptake rates, and the aerosol removal rate due to the presence of an occupant in the room, are not accounted for in this parameter. The IE, Equation (5), is used to analyze the effect of control measures on the relative reduction of virus-laden aerosol exposure of occupants under several scenarios. The normalized integrated exposure (NIE) is used to evaluate the relative reduction in virusladen aerosol exposure caused by the control measures. Hence, NIE for a given HVAC system, , is defined as the ratio of IE with and without control measures as given by Ng et al. [47]: where is the aerosol filtering efficiency (i.e., protection efficiency) of the mask worn by an exposed occupant, and the HVAC system can be CAHS, DOAS, or WUS.
The equivalent clean ACH (air changes per hour), given by Equation (7), is used to denote the "clean" airflow provided to the single-zone room as a combination of outdoor air supplied by the mechanical ventilation system (ACH ), outdoor air delivered via infiltration (ACH ), and the control measure ( ). The control measure can be a single measure or a combination of measures applied to provide the equivalent clean air (i.e., 0 ) to the zone, which can be HVAC and/or non-HVAC controls. The consists of enhanced filtration in the recirculated air using HVAC filters based on MERV rating, the addition of a PAC, or the use of a local exhaust fan. Aerosol deposition on interior surfaces is not taken into account in the . The estimated relative reduction in virus-laden bioaerosols exposure, in terms of obtained by using Equation (6), correspond to the equivalent clean ACH (and its value is calculated using Equation (7)) delivered by the applied control measures.
The outdoor air flowrate, , is the minimum requirement defined by the RITE. This depends on the building's use, which, in turn, defines the indoor air quality category and the minimum requirement of outdoor air per person (or floor area) [36]. Then, the ACH is calculated as the ratio of (in l/s) and zone volume ( in m 3 ): where is the outdoor air flowrate (l/s/person), and is the number of occupants in the zone (full occupancy), which may be obtained from the default occupancy based on the ASHRAE standard 62.192019 or local regulations for a specified space type [37]. The outdoor air infiltration rate, , through the building's envelope is obtained using a typical value for ACH , although infiltration is dependent on several factors (for instance, wind speed and direction, indoor/outdoor temperature difference, and HVAC system operation [48]) and is outside the scope of this study.
The supply airflow rate, , is determined from the peak cooling/heating loads estimated in the design conditions set for the highly occupied indoor space, which is modelled and simulated with the EnergyPlus TM simulation program [49]. Then, the outdoor air intake fraction, , of the ventilation system is calculated as the ratio of supply and outdoor air flow rates (i.e., , 1.0). The total outdoor air ACH, ACH , , can be calculated using [40]: where the balance airflow rate of the zone is given by Equation (10).
The recirculation airflow rate, , which is the amount of return air, if any, recirculated back to the zone by the air handling system is calculated as: According to recent guidelines for HVAC system operation in buildings during the COVID-19 pandemic, air should not be recirculated if at all possible [28]. If recirculation cannot be avoided, the recirculated air should be filtered before being mixed with the fresh outdoor air intake in the AHU. Thereby, different air filters, rated by MERV 4-16, are implemented in the FaTIMA model using the efficiency prediction model developed by Kowalski et al. [50]. However, even though it is validated with experimental data, there are some discrepancies between this filter efficiency model and the minimum requirements available in ASHRAE 52.2-2017 for the MERV rating [47]. Therefore, the overall/equivalent clean ACH due to filtration of the recirculated air is: where is the filter efficiency, which depends on the aerosol/particle size. Since the virus-laden aerosol concentration in the outdoor air is assumed to be zero, an outdoor air filter is not used in this study.
The clean ACH delivered by the air cleaner in a room, ACH , is calculated as: where is the aerosol filtering efficiency of the air cleaner obtained from the MERV rating charts for an aerosol diameter of 1 μm. and , are the actual fan airflow rate and the maximum fan airflow rate, respectively, of the air cleaner, and is the fraction of fan flow rates (i.e., , and 0 1), which makes it possible to use air cleaners with multiple airflow rate settings. Therefore, the clean air delivery rate, CADR, of the air cleaner, is calculated as: The removal of aerosols from the room air because of deposition on surfaces (i.e., floor, walls, ceilings, and others) is estimated by using the effective deposition rate, , defined as:  [40]. The value of determined by using Equation (15) does not include the aerosol deactivation rate. Since viruses may lose viability to cause infection during the airborne phase, the simulation model includes the biological decay rate of SARS-CoV-2 ( in h −1 ) estimated by using Equation (16) and a half-life ( ⁄ ) of 1.1 h [16,51].

Estimation of the Energy Consumption of the Control Measures
This subsection describes the methods used to estimate the annual energy consumption of each control measure, that is to say, mechanical ventilation using OA with several values of ( 0.0 1.0 ), recirculated air filtration with MERV rating (≥MERV 13), PACs, and local exhaust fans.
The energy consumption of ventilation using OA depends on whether it has to be cooled or heated to bring it to the desired level of comfort (i.e., temperature and relative humidity) and, therefore, is highly dependent on climate conditions (e.g., [52,53]). The influence of daily variation can be reduced by using an economizer, which automatically increases ventilation rates when the outdoor air conditions are favorable from an energy consumption viewpoint. However, such control strategies (e.g., free cooling at night) are not commonly implemented in existing buildings. Moreover, heat-recovery technologies can be implemented to minimize the energy impact of ventilation.
In this study, the energy consumption associated with ventilation and other heating/cooling loads in the indoor space is simulated using the EnergyPlus TM software , v9.3.0 [49], which provides a reliable estimate of the building's energy consumption based on the weather conditions of the location and the building use.
The amount of energy required to sustain the desired indoor conditions is associated with the cooling or heating demand of the indoor space. The cooling/heating demand is estimated by considering the gains and losses in building heat such as heat transmission via building envelope elements, infiltration air heat gains/losses, fresh OA ventilation heat gains/losses, and internal heat gains from occupants, solar radiation, artificial lighting, and electrical equipment. The thermal characteristics of the building envelope are regulated in Spain by the requirements of the Technical Building Code (CTE-DB-HE), which must be obeyed by both new and renovated buildings. The CTE-DB-HE has recently been modified to comply with the requirements for achieving nearly zero-energy buildings (NZEBs) [29]. The U-values of building elements (i.e., exterior wall, ground floor, flat roof, windows, and door) were chosen based on the limit values established by the current (after 2020) [29] and previous CTE-DB-HE (2006-2019) [30][31]34].
As far as the auxiliary components of the CAV HVAC system are concerned, the supply (and return) fan and motors were modelled using their typical efficiency values (i.e., 70% total fan efficiency and 90% electrical efficiency). The supply fan total efficiency and motor efficiency of the DOAS were also set to 70% and 90%.
The implementation of an HVAC filter in the recirculated air duct needs additional fan power because of the pressure drops raised due to the particulate filtration. The additional annual electrical energy consumption, in kWhel, can be reasonably estimated by considering the pressure drop of the recirculated airflow across the filter using [54,55]: where ∆ is the average pressure drop across the HVAC filter in Pa, ƞ and ƞ are the total mechanical efficiency (70%) and electric motor efficiency (65%), respectively, of the variable speed fan, and (in h) is the annual fan operational time. The fan was assumed to operate for an hour before the building was occupied and then operated throughout the period it was occupied.
For an in-room air cleaner (i.e., PAC), the electrical energy consumption ( in kWh/year) can be estimated from the device's rated power using Equation (18). Theoretical models can also be found in the literature for estimating and optimizing the overall particle removal rate for a given level of electricity consumption [56]. However, it is advisable to use power consumption data obtained from the device's manufacturer specifications at the maximum fan speed (and other fan speeds, i.e., flow rates, when variable speed operation devices are used) [43]. The PAC energy efficiency must be determined with a standard test, and it is defined as a clean airflow rate per unit of electrical power consumed (e.g., it will be assigned to one of the classes of Eurovent Certita Certification for Clean Air Efficiency, that is to say, between A-class (>13 m 3 /h/W) and Eclass (<2 m 3 /h/W)) [57,58]. According to a recent online survey on commercial PAC, the average electrical power consumed per CADR is about 0.15 W/(m 3 /h) [59].
where is the PAC specific energy consumption in W/(m 3 /h) and CADR the clean air delivery rate in m 3 /h of the air cleaner.
The electricity consumption of a local exhaust fan, , in kWh, can be estimated using the rated power of the fan motor indicated on its nameplate. Even though the amount of electricity consumed can be estimated more accurately when models calibrated using actual operating data are used, the nameplate power rating and design airflow rate of the fan also provide a reasonable estimation for this type of energy analysis.

Results and Discussion
Sections 4.1-4.3 present and discuss the simulation results obtained regarding the effect of HVAC and non-HVAC control on the reduction of exposure to virus-laden bioaerosols of uninfected occupants in a highly occupied indoor space (i.e., a typical bar or pub in Tarragona, Spain). The level of bioaerosol exposure on uninfected occupants sharing the same indoor space with an infected individual is investigated first in Section 4.1 in baseline conditions (Table 2) and with the other input modeling parameters listed in Table 1. Then, the effect of the control measures listed in Table 2 on the reduction of infections resulting from aerosol exposure is analyzed in an individual (Section 4.2) and combined form of implementation (Section 4.3). The energy consumption associated with each control measure is estimated and discussed.
The susceptibility of individuals exposed in indoor spaces to the presence of an asymptomatic or pre-symptomatic subject infected with SARS-CoV-2 has been examined in recent studies (e.g., [15,16,44]). The aerosol generation rate, , of an infected occupant was set to 686 #/min by considering the average adult exhalation rate and the number of exhaled aerosols (in through the nose and out through the mouth). The exhaled aerosols number is determined based on experiments carried out on healthy individuals (i.e., 98 #/min) [60]. This value does not assume that all the aerosols contain viruses or define a specific concentration of aerosol in occupied indoor space that is sufficient to cause COVID-19 [47]. In order to analyze the impact of each control measure, in this study, results are given in terms of relative reduction in infectious human aerosol exposure, so the main results are normalized values. Hence, the absolute value of the aerosol generation rate has little impact on the outcome of a study of this type and an arbitrary can also be used (e.g., 500 #/min used by Ng et al. [47]). Additionally, although the model can accommodate intermittent aerosol generation sources like coughing, such sources are not considered in this analysis within the occupancy duration (i.e., 3 h per day, 20:00-23:00). Table 2. Input modeling parameters to simulate the reduction of virus-laden bioaerosols exposure in a highly occupied indoor space (bar or pub) by using HVAC and non-HVAC control measures.  (Figure 2b). c WUS ( Figure  2c). d for particle (aerosol) size of between 0.3-1.0 μm. The MERV of HVAC filters is based on particle filtering efficiency values in [40,50]. n/a = not applicable (or not used). e based on values available in [47].

Description and Parameter
The supply airflow predicted by the EnergyPlus simulation, , from the internal and external heat gains of the fully occupied room (listed in Table 1) was about 620 l/s for the CAHS, 618 l/s for the DOAS, and 627 l/s for the WUS. According to Spain's CTE, the influence of the building's envelope , on the calculated was minor due to the high internal heat gains in the space (e.g., the CAHS increased by <1.0%, from 620 l/s to 625 l/s, when the was based on the CTE until 2019 instead of CTE-DB-HE1 2019). The simulated building data was obtained using typical meteorological year (TMY3) data from Tarragona (Spain).
The main drawback of the DOAS is that it is designed for a specific mechanical ventilation rate (e.g., 400 l/s in this case) and cannot readily accommodate considerably increased ventilation rates, while the WUS has no mechanical ventilation with outdoor air. Therefore, the increase in OAF was only applicable in the case of CAHS from about 0.65 up to one (i.e., 620 l/s). It is also assumed that the HVAC systems (i.e., CAHS, DOAS, and WUS) began to operate an hour before the building was occupied. Figure 4 shows the simulated thermal load (heating and cooling loads) and temperature profile (outdoor and indoor air) for the design days in summer and winter conditions and for the three HVAC systems studied. As the conditioned zone is highly occupied, the thermal loads are especially affected by occupancy density (i.e., of 100 #/100 m 2 ), resulting in a much higher requirement for cooling than heating. In this respect, heating requirements are mostly required on cold days to heat the outdoor airflow into the space.
In summer, the three HVAC systems have similar cooling requirements and profiles. However, in winter, their thermal loads and profiles are very different. In this regard, after the bar had been heated for the first hour after opening with CAHS, heating was still required during the following hours as the heating requirements of the outside air flow rate were higher than the internal cooling loads. On the other hand, the DOAS and the WUS were only required for heating during the first hour of opening, and, after that, cooling was required. For the WUS, no heating was required as no outside air was supplied. For the DOAS, no heating was required because the outside airflow did not need to be heated up as much as for the CAHS. Moreover, for the DOAS, the outside airflow rate decreased the cooling required in comparison with WUS.

Baseline Analysis
The design supply airflow rate of each HVAC type (i.e., CAHS, DOAS, and WUS) was based on the peak cooling loads. Table 3 shows the peak cooling and heating loads, and the annual thermal demands (i.e., cooling and heating) of the simulated indoor space under full occupancy (i.e., 50 individuals). Table 3. Simulated annual thermal demands (cooling and heating) and peak cooling and heating loads of a highly occupied indoor space served by three types of HVAC system (Figure 2 Figure 5a shows the normalized and integrated exposure (NIE and IE) of the bar/pub with a simulated volume of 150 m 3 for each type of HVAC system (Figure 2) based on the base-case conditions listed in Table 2. The exposure to aerosols of the room occupants, evaluated using IE and NIE, was estimated for 3 h/day of continuous exposure to a continuous bioaerosols source for the same duration. The NIE values reported in Figure  5a were obtained by normalizing the IE associated with each HVAC system type by using the lowest IE (which was for the CAHS). The base-case IE value of the WUS was over seven times higher than of the CAHS and DOAS, largely because OA was not supplied by the WUS ventilation system but was by CAHS and DOAS. The ACH of OA in the CAHS and DOAS were 9.97 h −1 and 9.90 h −1 , respectively, while for the WUS, it was only 0.3 h −1 (due to ACH of OA delivered by the infiltration). Therefore, the type of HVAC system clearly has an impact on the IE value, especially in the case of WUS (i.e., 2,024,900#•s/m 3 ), which has no mechanical OA supply system.
The IE reported in Figure 5a does not include the decay rate of the virus in aerosols (~0.63 h −1 ) and the aerosols deposited on surfaces (~0.067 h −1 ). The surface loading is detailed in Appendix A. In addition, the transient aerosol airborne concentration and the exposure within 24 h are shown in Appendix B.
The associated annual energy consumption (i.e., electricity and natural gas) of the three HVAC systems is shown in Figure 5b. This figure is the energy required by the HVAC system including the power required by the chiller and the auxiliaries and the natural gas required by the boiler. Finally, the heating and cooling required are provided by the cooling and heating coils described in Section 2.

Effect of Individual Control Measures
This section analyzes the effect of the control measures listed in Table 2 on the reduction of virus-laden aerosol exposure of uninfected occupants sharing the same indoor space with a COVID-19-infected individual (asymptomatic or pre-symptomatic).  Effect of the outdoor air fraction (OAF): Figure 6a shows the effect of increased mechanical ventilation using OA on the NIE for the CAHS versus OAF. The associated annual energy consumption (i.e., consumed electricity and natural gas) of this enhanced mechanical ventilation mitigation strategy, for several OAF values, is illustrated in Figure 6b. The equivalent clean ACH increased from 9.9 h −1 to 15.2 h −1 when the OAF was increased from ~0.65 to 1.00, at an outdoor airflow rate ranging from 400 l/s to 620 l/s. The OA ventilation rate per person increased from about 8.0 l/s/person to 12.4 l/s/person in the same range of OAF. The annual consumed electricity dropped by about 8.9% when the outdoor airflow rate was increased from the base case (i.e., 400 l/s) to 620 l/s, while natural gas consumption increased by more than twofold in the same range (Figure 6b). This is due to the 16.3% drop in annual cooling demand and the almost 1.9-fold rise in annual heating demand when the outdoor airflow rate increased from 400 l/s to 620 l/s (Figure 7). The peak cooling and heating demand of the fully occupied indoor space is also shown in Figure 7 at the different outdoor airflow rates (i.e., OAF) considered in this study.  Effect of recirculated air filtration: Figure 8a,b show the influence on the NIE of recirculated air filtered with HVAC filters with MERV ratings between MERV 13 and MERV 16 for the CAHS (Figure 2a) and DOAS (Figure 2b), respectively. The additional electricity consumed for recirculated air filtration (see Figure 8a,b) is estimated from the power consumption of the fan required to overcome the pressure drop caused by the HVAC filter, and the fan operational time (i.e., 4.0 h per day for 365 days per annum). The average pressure drops of the MERV-rated HVAC filters, reported in Azimi and Stephens [54], were used to estimate the power consumption of the fan. Since there is no HVAC filter for the return air stream of the WUS (Figure  2c), the influence of return/recirculated air filtration was not analyzed. In both CAHS and DOAS, the calculated NIE values were almost the same (0.68−0.70 and 0.43−0.45, respectively) when HVAC filters with MERV 13 and higher ratings were implemented for the aerosol size considered (i.e., a diameter of 1.0 μm). In addition, the NIE decreased more in the DOAS than in the CAHS as the recirculated airflow was higher. The equivalent clean ACH increased from 14.66 to 15.13 h −1 for the CAHS and 23.25 to 24.58 h −1 for the DOAS when the HVAC filter efficiency rating was raised from MERV 13 to MERV 16. In the CASH, the increase in electricity consumption when there were filters in the recirculation flow ranged from 3.4% to 4.5%, whereas in the DOAS it ranged from 9.0% to 11.7%.

Portable air cleaner (PAC):
The use of a PAC to reduce the virus-laden bioaerosol exposure of building occupants is an advantage when extra ventilation using outdoor air and/or recirculated air filtration is not possible, as in the case of WUS. Hence, the effect of PAC on the reduction of aerosol exposure was simulated with each type of HVAC system (i.e., CAHS, DOAS, and WUS, as shown in Figure 2) serving the indoor space. The effect of the PAC on NIE was simulated using six different CADR values, which are specified in Table 2 in terms of (i.e., 1.0-6.0 h −1 ). For the three HVAC types, Figure 9 compares the influence of PAC on the reduction of aerosol exposure to the base-case in terms of NIE. When the was increased from 1.0 h −1 to 6 h −1 , the maximum airflow rate via the PAC was approximately 42 l/s to 253 l/s with a filtering efficiency of 99% for an aerosol diameter of 1.0 μm. The electricity consumed by the PAC to deliver the clean ACH is also depicted in Figure  9. The influence of the PAC CADR on the NIE is similar in the CASH and DOAS (it decreases in both from 0.92, when the CADR is 150 m 3 /h, to 0.65, when the CADR is 900 m 3 /h), and much higher in the WUS (it decreases from 0.61, when the CADR is 150 m 3 /h, to 0.20, when the CADR is 900 m 3 /h). The increase in energy consumption with the PAC unit is small, always below 200 kWh/year.


Exhaust fan: The commercially available box fans can be installed in a room (e.g., in windows or doorways) to exhaust contaminated air out of the indoor space, thus reducing the exposure of occupants to infectious bioaerosols. The effect of the exhaust fan on NIE, which was simulated with different capacities (maximum airflow rates, Table 2), is shown in Figure 10. Electricity consumed by the exhaust fans can be estimated as it was for the PAC (Figure 9) by using the operating hours and rated power of the box fan, indicated on its nameplate, existing in the market.  Figure 11a illustrates the combined effect of control measures on the reduction of virus-laden aerosol exposure to uninfected occupants for each type of HVAC system. These control measures include mechanical ventilation using OA based on the RITE standard (i.e., an OA air supply of 8 l/s/person), enhanced recirculated air filtration using a MERV 13 filter in the CAHS and DOAS (90% efficiency for considered aerosol diameter of 1.0 μm), and PACs of several capacities (i.e., 150-900 m 3 /h of CADR). The annual electricity consumption of these control measures (Figure 11a) is depicted in Figure 11b.

Effect of Combined Control Measures
The NIE values and normalized electricity consumption reported in Figure 11a,b, respectively, were obtained by normalizing the corresponding IE and annual electricity consumption of each combined control measure with the base-case values.
As well as consuming electricity, CAHS and DOAS also consume 619.44 kWh/year and 234.67 kWh/year of natural gas, respectively, to supply the heating requirement of the conditioned indoor space.
The operational costs of the control measures can be estimated from the electricity and natural gas consumption reported in the present study and the current electricity and natural gas prices. The cost of electricity to run PACs is cheaper than other control measures, such as ventilation using outdoor air if the required aerosol mitigation is reached by using the PAC. For instance, the average yearly cost of electricity to run a portable HEPA air cleaner all the time is between USD 100 and 250 [61]. Filter replacement costs should also be included in estimating the annual overall cost of the PACs. Finally, it should be noted that other technologies such as the ultraviolet germicidal irradiation (UVGI) system and photocatalytic oxidation (PCO) for reducing aerosol transmission have not been included in this study, but they also have the potential to mitigate airborne virus propagation in indoor spaces. Nor does the modeling and simulation methodology include non-airborne routes of exposure (e.g., direct/indirect contact with fomites), so the control measures discussed should be implemented together with other strategies (e.g., hand washing).

Concluding Remarks
Three common types of HVAC system were analyzed for controlling the spread of virus-laden bioaerosols in a highly occupied indoor space (floor area of 50 m 2 and an occupancy density of 100 #/100 m 2 ) in Tarragona (Spain) weather conditions. The space was analyzed for 3 h/day of continuous exposure of occupants sharing the space with an infected subject. The HVAC system types were a central air handling system (CAHS), a dedicated outdoor air system (DOAS), and a wall unit system (WUS). The base-case occupant bioaerosol exposure for CAHS and DOAS was studied by considering the effect of mechanical ventilation with an outdoor airflow rate in compliance with the Spanish regulation (i.e., RITE standard, 8 l/s/person), whereas the WUS had no mechanical ventilation. The combined effect of mechanical ventilation (i.e., 400 l/s of outdoor air), recirculated air filtration using a MERV 13 that is 90% efficient (for an aerosol diameter of 1 μm), and a portable air cleaner (PAC) with a capacity in the range 150-900 m 3 /h of CADR mitigated the (normalized) integrated exposure of the occupant by ~0.66 to 0.51 and 0.43 to 0.36 for the CAHS and DOAS, respectively. In the case of WUS, the normalized integrated exposure was reduced by up to 0.2 when the PAC with a capacity of up to 900 m 3 /h was used. The electricity consumed by these combined control measures increased by approximately 297.37 kWh/year and 482.66 kWh/year for the CAHS and DOAS, respectively, in comparison to the base-case, while for the WUS, it increased by about 197.10 kWh/year. Even though the simulated occupant exposure (in terms of IE and NIE) is not a direct measure of the infection risk of occupants sharing a space with an infected subject, it can be used to effectively identify the best strategies for reducing building occupant exposure to airborne pathogens. Therefore, estimating the energy consumed by the control measures is useful for selecting cost-effective control strategies for the different types of HVAC systems in buildings.
In conclusion, the methodology in this study helps to select the best control options for the type of HVAC system in a building and reduce energy consumption. Thereby, control solutions can be implemented economically, and the indoor air quality can be guaranteed to be sufficient for human health. Thereby, risks related to the airborne transmission of COVID-19 (and alike diseases) can be effectively mitigated in highly occupied indoor spaces in hotels including bars and pubs by applying cost-effective control measures. Thus, other imposed non-engineering control measures can be relaxed/removed (such as reduction in occupancy level) due to COVID-19 preventive measures so that indoor activities can be increased to their full capacity.

Funding:
The research work presented in this paper is the outcome of a project funded by Diputació de Tarragona under the collaboration framework agreement between the Diputació de Tarragona and the Universitat Rovira i Virgili for the period 2020-2023 with the reference number 2020/19: reference number 2020/19. Acknowledgments: John F. Bates, coordinator of the language service at the Universitat Rovira i Virgili, is acknowledged for his editing contribution to the final version of the manuscript.

Conflicts of Interest:
The authors declare no conflict of interest.