The present invention relates to heating, ventilation, and air conditioning (HVAC) systems and, more specifically, to an approach for controlling an HVAC system to minimize indoor infection risk and the concentrations of air contaminants of interest while maximizing energy savings through optimal ventilation, tunable air cleaning and temperature control and their integration.
An indoor environment is traditionally maintained by the heating, ventilation, and air conditioning (HVAC) system of using appropriate temperature and ventilation controls. Conventional control approaches focus on how to improve indoor environment through ventilation controls or save energy through model predictive controls. These control approaches do not, however, take into account the need to minimize infection risk indoors and in order to reduce infection risk most buildings will simply increase outdoor air intake and/or install higher performance air filters, both of which increases the energy consumption of HVAC system. Accordingly, there is a need for an HVAC control system that can minimize indoor infection risk and the concentrations of air contaminants of interest while maximizing energy savings through optimal ventilation, tunable air cleaning and temperature control.
The present invention is a system and method for minimizing indoor infection risk and improving indoor air quality by reducing concentrations of pollutants of interest while maximizing energy savings that integrates real-time occupancy detection and forecasting that determines the current and future occupancy in a space (presence and people counting), outdoor weather condition forecasting that provides current and future outdoor weather information, indoor infection risk and air quality models that model the infection risk assuming an imperfectly mixed realistic indoor air environment as well as the concentrations of pollutants of interest for indoor air quality (IAQ) control, tunable air filtration/purification/disinfection technologies with different efficiency that provides different filtration/purification/disinfection efficiency and clean air delivery rate in the HVAC system for removing and diluting the virus-containing particles, and any portable air purifier/cleaner devices that provide the infectious particle removal rate (or fresh air supply rate) by air cleaners. The present invention outputs the total amount of outdoor air intake, the air temperature of the supply air into the space, the supply air flow rate into the space, the operation mode of tunable air filtration/purification/disinfection, the operation mode of the in-room air cleaner, and space/room temperature set-points, and thus can serve as the central controller for an entire HVAC system. The present invention can thus transform existing model predictive control capability to meet indoor infection risk control and indoor air quality (IAQ) requirements while maximize the energy savings.
The present invention will be more fully understood and appreciated by reading the following Detailed Description in conjunction with the accompanying drawings, in which:
Referring to the figures, wherein like numeral refer to like parts throughout, there is seen in
Control system 12 is based on a predictive control strategy to integrate occupancy prediction, weather forecasting, indoor infection risk modeling and building automation system. The disclosed overall architecture leverages building automation system information through BACnet™, and develops MPC platform to minimize indoor infection risk while maximizing energy savings. The system architecture is designed to compute and implement the optimized control strategy by forecasting future states of occupancy (presence and number of occupants), forecasting future states of local weather (temperature and solar radiation), forecasting future states of ambient air quality, forecasting future room heating and cooling loads, forecasting future states of room temperature and indoor air quality, communicating with building automation system to obtain current and historical room temperature and indoor air quality, and communicating with building automation system to obtain current operation mode of tunable air filtration/purification/disinfection. All computations may be performed by cloud computing or by local computers or controllers or any types of processors. System implement may be performed Internet of Things (IOTs) enabled components. Control system 12 is designed to fulfill the following functions: optimize total amount of outdoor air intake amount and its schedule; optimize supply air temperature into the space and its schedule; optimize supply air flow rate into the space and its schedule; optimize room/space temperature setpoint and its schedule; optimize operation mode of the tunable air filtration/purification/disinfection and its schedule; and optimize operation mode of the in-room air cleaner on/off and its schedule.
Control system 12 includes a multi-objective model predictive control design, energy and IAQ management and optimization using modeling components in overall model predictive control design. The overall model predictive control modeling components include an occupancy prediction model, a local weather (temperature and solar radiation forecasting) model, a local ambient air quality forecasting model, a room/space heating and cooling load forecasting model, an indoor infection risk model, an indoor air quality model, a building physics-based model, and an HVAC physics-based model. These components are seen in
Referring to
The integrated control, energy and IAQ management, and optimization system makes decisions based on: infection risk and IAQ models that predict the infection risk and the concentration of the indoor air contaminants of interest based on outdoor air intake flow rate and air quality, total supply airflow rate, operation mode and efficiency of the tunable air filtration/purification/disinfection technologies, in room air cleaner and type of mask; room/space load forecasting; occupants' thermal comfort preferences; occupants' schedule forecasting; room/space control state estimation; room/space air quality states; room/space operational constraints; and weather forecasting. The integrated control, energy and IAQ management, and optimization system communicates these decisions and constraints by dispatching heating, cooling and ventilation set-points (e.g., thermostat set-points, air handling unit airflow rate, and supply air temperature set-points) to the local device controllers. The system has decision-making capabilities to control heating, cooling and ventilation system over a time horizon.
The system has an optimization engine that computes schedules of outdoor air intake airflow rate, supply air temperature, supply airflow rate, room set points for temperature, RH and target pollutant concentrations, operation mode of the tunable air filtration and operation mode of the in-room air cleaner based on the information it collects. The main features of the optimization engine include: a multi-objective model predictive control architecture, and a mixed integer programming formulation to solve both continuous and discrete equipment controls such as in-room air cleaners. The disclosed system has three different approaches to solve the MPC problem and implement control outputs in real-time. The user can choose any of those three approaches based on the availability of onsite computing resources.
First, is an online centralized MPC at each AHU level as shown in
Second, an online distributed MPC at zone level as shown in
Third, offline learning as shown in
Component Level Modeling according to the present invention includes a Room/Space Thermal Model, where A room l can be represented using a 3R2C model, as seen in
Where {dot over (Q)}room is the total cooling or heating load injected into the room by the HVAC system, calculated as:
{dot over (Q)}
room
=c·{dot over (m)}
room(Tsupply,room−Tset) (2)
Further, the mass flow rate and air supply temperature to the room are functions of the damper position δdamp, and reheating coil valve position δvav_rh:
{dot over (m)}
room
=f
damp(δdamp) (3)
T
supply,room
=f
T(δvav_rh,δdamp) (4)
Substituting (3), (4) into (2) results in the room cooling or heating load being expressed in terms of three control variables δvav_rh, δdamp, Tset:
{dot over (Q)}
room
=f
Q(δvav_rh,δdamp,Tset)=c·fdamp(δdamp)·(fT(δvav_rh,δdamp)−Tset) (5)
Thus, the dynamics of the room (1) can be expressed in a semi-linear form:
These room dynamics can be written in state space form as follows:
{dot over (x)}
l
t
=A
l
x
l
t
+B
l
f(ult)+Elwlt∀l∈,∀t∈ (7)
where xlt=[]lT is the state of the room (i.e, wall and zone temperatures); f(ult)=fQ(δvav_rh,δdamp,Tset)lt is a nonlinear function of the three room-specific control input variables (VAV reheating coil position, damper position, setpoint temperature); wlt=[]lT the uncontrollable input at time t, comprised of ambient temperature, heat gains due to solar radiation, and internal heat gains due to occupants, lights, and equipment; and the system matrices Al, Bl, El are as follows:
where R1, R2, Rwin, C, Czone are the building thermal resistance and capacitance values and μ is the coefficient of performance of the AHU.
The HVAC model is based on the AHU which serves the entire building provides a total cooling or heating load according to its mass flow rate and temperature differential:
{dot over (Q)}
AHU
=c·{dot over (m)}
AHU(Tsupply,AHU=Tmix) (9)
The mass flow rate is a mixture of outdoor air intake and return air, with α being the fraction made up by outdoor air or outdoor air intake percentage:
{dot over (m)}
AHU
=α·{dot over (m)}
out+(1−α)·{dot over (m)}return (10)
Therefore by substituting (10) into (9), the air side total cooling or heating load can be expressed in terms of AHU-specific control variables uAHU=[Tsupply,AHU,{dot over (m)}AHU,α]:
{dot over (Q)}
AHU
=f
AHU(Tsupply,AHU,α)=c·(α·{dot over (m)}out+(1−α)·{dot over (m)}return)·(Tsupply,AHU−Tmix) (11)
The AHU model and room models are coupled through a mass flow balance (i.e., the total AHU flow rate is the sum of all n room flow rates):
{dot over (m)}
AHU
={dot over (m)}
1
+{dot over (m)}
2
+ . . . +{dot over (m)}
l
+ . . . +{dot over (m)}
n (12)
Finally, the total HVAC fan power use is calculated as a function of the total flow rate:
p
fan
=f
fan({dot over (m)}AHU) (13)
The Infection Risk Model is used to determine infection risk, which can be quantified by the well-known Wells-Riley equation [2]:
where NC=the number of new cases in the space; NS=the number of susceptible people; I=the number of virus carrier at the start of the exposure period; q=the infectious quantum generation rate per virus carrier (quanta/h per person); p=pulmonary ventilation rate (m3/h); t=exposure time (h); V=space air volume (m3); Λ=the fresh air change rate in the room (1/h).
Equation (5) establishes the relationship between the probability of infection (i.e., infection risk) in a perfectly mixed air space and the average viral dose exposure as:
P=1−e−D (15)
Where D=viral dose exposure (quanta). Note that 1 quantum is the amount of viral dose exposure needed to result in a probability of infection of 63% per Equation (5). The dose exposure for an individual in the space is calculated as follows:
Where RS and RL=fraction of infectious particles passing through the mask worn by the virus carriers and the succeptibles, respectively; and Λ=the “total equivalent” clean air change rate for the space in terms of infection virus dilution due to the air cleaning devices in central HVAC system and/or within the room space as well as that from the outdoor air intake flow rate (outdoor ventilation rate). Equations (6) and (7) are used to predict the risk of infection for an individual in an indoor space accounting for the effects of ventilation, HVAC system supply and/or recirculated air cleaning, in-room air cleaners, and mask wearing [3].
The fraction of infectious particle penetrated through the mask or respirator for susceptible (RS) and infected (RI) population can be calculated by Eqn. 8 and 9, respectively. And both depend on the mask filtration efficiency (ηs or ηI). The penetration fraction (R) equals 1 when no mask or respirator is used during the exposure period. An additional fractional factor (fR) is multiplied by the original filtration efficiency of the mask to represent the fraction of time using a mask/respirator over the entire exposure period. It equals 1 when the mask is worn during the entire exposure period.
R
S=1−fR,SηS (17)
R
I=1−fR,IηI (18)
The “total effective” air change rate depends on the “effective” air change rates due to ventilation rate (λvent), pathogen inactivation rate (kUV) by ultraviolet germicidal irradiation (UVGI) systems, infectious particle deposition rate (kdeposition) and pathogen natural inactivation rate in the air (kinactivation), as shown in Eqn. 10. The “effective” air change rate includes the air change rate due to fresh air supply rate by the HVAC system (λHVAC or {dot over (m)}AHU), natural ventilation rate (λNV) and infectious particle removal rate by air purifiers (kpurifier). The fraction (f) of operation time over the entire exposure period is applied to each term in Eqn. 11 to determine the net overall ventilation rate. The fresh/clean air change rate supplied by the HVAC system (λHVAC or in {dot over (m)}AHU) includes the outdoor part and the recirculated part. The recirculated fresh/clean air change rate (Eqn. 12) depends on the recirculated air change rate (λrecirculated) and the filtration efficiency of the filters in the HVAC system for the virus-containing particles (ηfilter). The natural ventilation air change rate (λNV) results from the airflows through openings and cracks on the building envelope.
Λ=λvent+fUVkUV+kdeposition+kinactivation (19)
λvent=fHVACλHVAC+fNVλNV+fpurifierkpurifier (20)
λHVAC=λoutdoor+λrecirculatedηfilter (21)
An in-room air purifier can supply additional fresh/clean air to the space. The infectious particle removal rate (or fresh air supply rate) by air purifiers (kpurifier) can be estimated by its airflow rate (λpurifier) and filter efficiency (ηpurifier), or based on its clean air delivery rate (CADR) and room volume (V):
The actual ventilation rate (λvent) is a spatial-variable that depends on the particular location in the space. Thus, an additional factor is applied to the original equation in Eqn. 11 to adapt it to imperfect mixing scenarios. The infectious particle removal rates due to UVGI systems, deposition, and natural inactivation are assumed to be uniform in the whole space. Therefore, these terms do not have to be modified for imperfect mixing.
The pathogen removal rate by the UVGI system depends on the fraction of UVGI operation time (fUV) and the pathogen inactivation rate due to UV irradiation (kUV). The infectious particle deposition rate (kdeposition) relies on an approximate estimate of gravitational settling (Eqn. 14) from Nicas et al. [4], which depends on the particle diameter (dp) and room height (H). It is assumed that the deposited particles will not be resuspended into the air space again.
The pathogen natural inactivation rate is not considered, in part because of the lack of existing data on the size-resolved natural inactivation rate of SARS-CoV-2 and in part because quanta generation rates (q), when back-calculated using Eqn. 1, will inherently account for any inactivation that occurred during the case study period [5].
The IAQ model predicts the concentrations of pollutants of interest in the occupied zone and estimate the health risks associated with the occupant exposure to the various pollutants. The concentration of an individual pollutant (i) of interest is governed by the following zone mass balance equation:
Where V=zone volume (m3); C=concentration (ug/m3); Ae=emission surface area (m2); E(t) emission factor (i.e., emission rate per unit surface area, ug/h/m2); Q=supply air volumetric flow rate to the zone (m3/h) which is
Cs=concentration of the supply air which include the effects of ventilation and air cleaning in the central HVAC system (ug/m3); As=Area of sink surfaces (m2); S(t) sorption/filtration/cleaning rate by sinks including in-room air cleaners (ug/h/m2).
The model is further extended to multizone buildings to account for the inter-zone pollutant transport [6]. The model is used to predict the concentrations over a time of horizon in the MPC model, and the predicted concentrations are compared to pre-established threshold limits for each pollutant of interest such as formaldehyde, PM2.5, CO2, etc. for maintaining satisfactory IAQ in the buildings, or for estimating the health risk associated with the cumulative exposure to the pollutants.
The multi-objective model predictive control can implement a centralized MPC approach where a joint optimization objective can be designed to minimize the weighted sum of total AHU fan energy, total air side heating/cooling energy in an AHU, total zone reheating energy (if applicable), zone infection risk, and the concentrations of pollutants of interest (if applicable) for a prediction horizon h:
J
t=ρf+ρcoil+ρphz+ρrisk +ρcpct (25)
Each term is weighted by a cost scalar ρ(⋅), and the overall objective is a time average of the sum of the five weighted terms. These weights and those in subsequent objective functions are tuned to achieve optimal control performance based on training data. Once tuned, the weights may remain static unless the underlying system models (i.e., room or AHU models) are changed. In addition, any cost function from utility companies could be added into all energy terms in Eq. (25). Thus, a centralized MPC controller will solve the following multi-objective problem for all room at an AHU level, over some prediction horizon h:
where =[] are new control variables at the AHU level, and where =[,,,, Rs,l,Rl,l] are new control variables at each room 1.
In the distributed MPC approach, a joint optimization objective can be designed to minimize the weighted sum of total air side heating/cooling energy in each room, total zone reheating energy (if applicable), zone infection risk, and the concentrations of pollutants of interest (if applicable) for a prediction horizon h, for each room l:
Jlt=ρroom+ρrhz+ρrisk+ρcpct (27)
Each term is weighted by a cost scalar ρ(⋅), and the overall objective is a time average of the sum of the four weighted terms. Thus, for each room, a MPC controller will solve the following multi-objective problem over some prediction horizon h:
where =[,,,,,Rs,l,Rl,l] are new control variables at each room l. is the required outdoor air portion of the total fresh air into each zone/room l.
As described above, the present invention may be a system, a method, and/or a computer program associated therewith and is described herein with reference to flowcharts and block diagrams of methods and systems. The flowchart and block diagrams illustrate the architecture, functionality, and operation of possible implementations of systems, methods, and computer programs of the present invention. It should be understood that each block of the flowcharts and block diagrams can be implemented by computer readable program instructions in software, firmware, or dedicated analog or digital circuits. These computer readable program instructions may be implemented on the processor of a general purpose computer, a special purpose computer, or other programmable data processing apparatus to produce a machine that implements a part or all of any of the blocks in the flowcharts and block diagrams. Each block in the flowchart or block diagrams may represent a module, segment, or portion of instructions, which comprises one or more executable instructions for implementing the specified logical functions. It should also be noted that each block of the block diagrams and flowchart illustrations, or combinations of blocks in the block diagrams and flowcharts, can be implemented by special purpose hardware-based systems that perform the specified functions or acts or carry out combinations of special purpose hardware and computer instructions.
Number | Date | Country | |
---|---|---|---|
63122460 | Dec 2020 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 17544383 | Dec 2021 | US |
Child | 18209167 | US |