Predicting spherical irradiance for volume disinfection

Abstract

Given the complexity of architectural spaces and the need to calculate spherical irradiances, it is difficult to determine how much ultraviolet radiation is necessary to adequately kill airborne pathogens. An interior environment with luminaires is modeled. Spherical irradiance meters are positioned in the model and the direct and indirect spherical irradiance is calculated for each sensor. From this, an irradiance field is interpolated for a volume of interest, and using known fluence response values for killing pathogens, a reduction in the pathogens is predicted. Based on the predicted reduction, spaces are built accordingly, and ultraviolet luminaires are installed and controlled.

Claims

1. A method for reducing pathogens in an interior environment, the method comprising the steps of: positioning, by a processor, spherical irradiance meters in a model of the interior environment; calculating, by the processor, for each spherical irradiance meter, a total spherical irradiance; interpolating, by the processor, the total spherical irradiances to result in a three-dimensional spherical irradiance field; calculating, by the processor, using the three-dimensional spherical irradiance field and a fluence response value of a pathogen, a predicted reduction in the pathogen in the interior environment; and controlling at least one ultraviolet luminaire in the interior environment to provide the predicted reduction.

2. The method of claim 1, comprising prior to the positioning step: preparing the model as a three-dimensional computer-aided drafting model; and calculating a radiosity solution for the interior environment.

3. The method of claim 1, wherein: the spherical irradiance meters are positioned in an array in a volume of interest for air flow calculations; the predicted reduction is calculated for the volume of interest; the interior environment comprises the at least one ultraviolet luminaire and zero or more further ultraviolet luminaires; and for each spherical irradiance meter, the total spherical irradiance is calculated by summing: direct spherical irradiance on the spherical irradiance meter due to all of the ultraviolet luminaires that are either fully or partially visible to the spherical irradiance meter; and indirect spherical irradiance on the spherical irradiance meter due to all surface patches of the model that are fully or partially visible to the spherical irradiance meter.

4. The method of claim 3, wherein calculating the predicted reduction comprises calculating, using the three-dimensional spherical irradiance field, a total radiant fluence for each of multiple finite volumes of air in the volume of interest.

5. The method of claim 1, wherein: the interpolation step is performed trilinearly or tricubically; and the predicted reduction is expressed as a log.sub.10 reduction.

6. The method of claim 1, comprising, after calculating each total spherical irradiance: marking the spherical irradiance meter as having been processed; and searching for a spherical irradiance meter that is unmarked.

7. The method of claim 1, wherein the total spherical irradiance is ultraviolet total spherical irradiance.

8. The method of claim 1, wherein the interior environment comprises air or water.

9. A method for predicting spherical irradiance for disinfecting an interior environment, the method comprising the steps of: positioning, by a processor, spherical irradiance meters in a model of the interior environment; calculating, by the processor, for each spherical irradiance meter, a total spherical irradiance; interpolating, by the processor, the total spherical irradiances to result in a three-dimensional spherical irradiance field; calculating, by the processor, using the three-dimensional spherical irradiance field and a fluence response value of a pathogen, a predicted reduction in the pathogen in the interior environment; mounting an irradiation sensor in an ultraviolet luminaire in the interior environment; periodically measuring, using the irradiation sensor, a radiant flux output of the ultraviolet luminaire; using the radiant flux output and the model to recalculate the three-dimensional spherical irradiance field; and comparing the recalculated three-dimensional spherical irradiance field or a recalculated predicted reduction in the pathogen to a threshold that results in an alarm being triggered if met.

10. A system for reducing pathogens in an interior environment, the system comprising: a processor; a non-transient computer-readable memory storing instructions, which, when executed by the processor, cause the processor to: position spherical irradiance meters in a model of the interior environment; calculate, for each spherical irradiance meter, a total spherical irradiance; interpolate the total spherical irradiances to result in a three-dimensional spherical irradiance field; and calculate, using the three-dimensional spherical irradiance field and a fluence response value of a pathogen, a predicted reduction in the pathogen in the interior environment; and one or more ultraviolet luminaires that provide the predicted reduction in the interior environment.

11. The system of claim 10, comprising an irradiation sensor in one or more of the one or more ultraviolet luminaires, wherein the processor is configured to recalculate the total spherical irradiances using a measurement from each irradiation sensor.

12. A system for predicting spherical irradiance for disinfecting an interior environment, the system comprising: a processor; a non-transient computer-readable memory storing instructions, which, when executed by the processor, cause the processor to: position spherical irradiance meters in a model of the interior environment; calculate, for each spherical irradiance meter, a total spherical irradiance; interpolate the total spherical irradiances to result in a three-dimensional spherical irradiance field; and calculate, using the three-dimensional spherical irradiance field and a fluence response value of a pathogen, a predicted reduction in the pathogen in the interior environment; and a controller which is connected to the processor and controls one or more ultraviolet luminaires in the interior environment to provide the predicted reduction.

13. The system of claim 12, connected to an occupancy sensor, an information feed, a timer, a user interface, a power meter, an HVAC (heating, ventilation and air conditioning) system, an energy storage system, an energy management system, an artificial intelligence engine, a utility company, a fenestration device or a visible luminaire, or any combination selected therefrom, and configured to provide minimum occupant requirements for occupants in the interior environment.

14. The method of claim 9, comprising prior to the positioning step: preparing the model as a three-dimensional computer-aided drafting model; and calculating a radiosity solution for the interior environment.

15. The method of claim 9, wherein: the spherical irradiance meters are positioned in an array in a volume of interest for air flow calculations; the predicted reduction is calculated for the volume of interest; the interior environment comprises the ultraviolet luminaire and zero or more further ultraviolet luminaires; and for each spherical irradiance meter, the total spherical irradiance is calculated by summing: direct spherical irradiance on the spherical irradiance meter due to all of the ultraviolet luminaires that are either fully or partially visible to the spherical irradiance meter; and indirect spherical irradiance on the spherical irradiance meter due to all surface patches of the model that are fully or partially visible to the spherical irradiance meter.

16. The method of claim 15, wherein calculating the predicted reduction comprises calculating, using the three-dimensional spherical irradiance field, a total radiant fluence for each of multiple finite volumes of air in the volume of interest.

17. The method of claim 9, wherein: the interpolation step is performed trilinearly or tricubically; and the predicted reduction is expressed as a log.sub.10 reduction.

18. The method of claim 9, comprising, after calculating each total spherical irradiance: marking the spherical irradiance meter as having been processed; and searching for a spherical irradiance meter that is unmarked.

19. The method of claim 9, wherein the total spherical irradiance is ultraviolet total spherical irradiance.

20. The method of claim 9, wherein the interior environment comprises air or water.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1, prior art, shows an empty rectangular room divided into an upper room with a UVGI luminaire and an occupied lower room.

(2) FIG. 2, prior art, shows the ultraviolet hazard function as a representation of the DNA spectral absorptance distribution from the German Institute for Standardization (DIN) and the Illuminating Engineering Society of North America (IESNA).

(3) FIG. 3, prior art, illustrates the difference between the irradiance of a point on the surface of a plane and the spherical irradiance (fluence rate) of an elementary sphere.

(4) FIG. 4 shows an empty rectangular room with one example surface subdivided into patches.

(5) FIG. 5 shows the spherical irradiance received at a point from a single patch.

(6) FIG. 6 shows a dual cubic tetrahedron as a component of a light buffer.

(7) FIG. 7 shows the subdivision of a cubic tetrahedron face into an array of elements.

(8) FIG. 8 shows the orientation of a light buffer with respect to the u-v-n Cartesian coordinate system.

(9) FIG. 9 shows a patch being projected onto the face of a cubic tetrahedron.

(10) FIG. 10 shows a flowchart of a method of predicting a three-dimensional spherical irradiance field in a virtual environment, according to an embodiment of the present invention.

(11) FIG. 11 is a flowchart for various implementations of the method, according to an embodiment of the present invention.

(12) FIG. 12 is a schematic block diagram of a system for predicting UV spherical irradiance, according to an embodiment of the present invention.

(13) FIG. 13 is a schematic block diagram of a system for disinfecting an interior environment, according to an embodiment of the present invention.

(14) FIG. 14 is a flowchart of a method for ongoing monitoring of the spherical irradiance field, according to an embodiment of the present invention.

DETAILED DESCRIPTION

(15) Referring to FIG. 4, an empty rectangular room 400 has its surfaces subdivided into an array of patches 410 that are irradiated by an ultraviolet luminaire 405. Known radiosity methods as described in, for example, Ashdown, I. 1994. Radiosity: A Programmer's Perspective. New York, N.Y.: Wiley, may be employed to calculate the radiant exitance M of each patch 410 due to direct irradiance E.sub.o from one or more luminaires 405 and indirect irradiance E.sub.i received from all other patches visible to the patch and multiplied by the patch diffuse reflectance ρ.

(16) A virtual spherical irradiance meter 415 is positioned at an arbitrary point P in the space. If the meter is at a distance r that is approximately two or more times the diagonal width of the emission area of luminaire 405, the luminaire can be represented as a point source of irradiation at the center of the emission area, and the direct irradiance E.sub.o of the spherical irradiance meter calculated using the inverse square law:
E.sub.o=I.sub.P/r.sup.2  (3)
where I.sub.P is the radiant intensity of the luminaire in the direction of the point P.

(17) If the distance r is less than the twice the diagonal width of the emission area of luminaire 405, the emission area may be subdivided into a one- or two-dimensional array of smaller subareas such that the diagonal width of each subarea satisfies the distance criterion, and the total radiant flux Φ emitted by the luminaire is evenly divided between them.

(18) The inverse square law calculation is repeated for each luminaire 405 that is directly visible from the point P, to determine the total direct spherical irradiance E.sub.o at point P.

(19) To calculate the indirect spherical irradiance E.sub.i from each patch 410 visible to point P, and referring to FIG. 5, it is first noted that the radiant flux E received at point 505 from a patch 510 with surface area A and radiant exitance M.sub.e that is at a distance r and at an angle φ to the surface normal Nis:

(20) $\begin{array}{cc}E=\frac{1}{\pi }{\int }_{\Omega }{M}_{e}d\omega & \left(4\right)\end{array}$
where the differential solid angle dω is defined as:
dω=cos(φ)dA/r.sup.2  (5)

(21) The integration can be performed by numerical quadrature using a “light buffer.” This consists of two back-to-back cubic tetrahedrons measuring two units on each side that are centered on the point P (FIG. 6). Each face of the cubic tetrahedrons is divided into an array of n=m×m/2 cells with area ΔA.sub.m, as shown in FIG. 7. (The cells along the adjoining edges of the two tetrahedrons have area ΔA.sub.m/2.) Thus, for a given patch j, and a given face of the convex polyhedron from which the patch j is visible:

(22) $\begin{array}{cc}{E}_i\approx \frac{1}{\pi }\underset{n}{.Math.}{M}_e\cos \left({\phi }_j\right)\Delta A_m/r_j^2& \left(6\right)\end{array}$

(23) FIG. 8 shows an orthographic projection of the upper cubic tetrahedron of a light buffer aligned on a horizontal plane with a u-v-n Cartesian coordinate system centered on point P and oriented such that each axis intersects a face of the cubic tetrahedron at a distance of unity from point P.

(24) Taking the cubic tetrahedron face intersected by the n-axis as an example, we have for each cell with variable coordinates u.sub.n and v.sub.n:
r.sub.n=√{square root over (u.sub.n.sup.2+v.sub.n.sup.2+1)}  (7)
and:
cos(φ.sub.n)=1/r.sub.n  (8)
Thus, for a given patch, its contribution of indirect spherical irradiance to a given cell in a face of the cubic tetrahedron from which the patch is visible is:
ΔE.sub.i≈M.sub.eΔA.sub.m/π(u.sub.n.sup.2+v.sub.n.sup.2+1).sup.3/2  (9)
and similarly for all six faces of the light buffer.

(25) This can be confirmed by assuming a uniform sphere with unit radiant exitance surrounding the light buffer and centered on point P. A diffusely-emitting surface with radiant exitance M has radiance M/π in any direction, which gives:
ΔE.sub.i≈ΔA.sub.m/(u.sub.n.sup.2+v.sub.n.sup.2+1).sup.3/2  (10)

(26) Summing ΔE.sub.i over all n cells (in the limit ΔA.sub.n.fwdarw.0) for each side yields 4π, which is the area of a unit sphere.

(27) Referring to FIG. 8, the cubic tetrahedron faces are oriented with respect to the u-v-n Cartesian coordinate system such that:
A={−√{square root over (⅙)},√{square root over (½)},0}
B={−√{square root over (⅙)},−√{square root over (½)},0}
C={√{square root over (⅔)},0,0}
D={0,0,0}
u.sub.CT={√{square root over (⅙)},√{square root over (½)},±√{square root over (⅓)}}
v.sub.CT={√{square root over (⅙)},−√{square root over (½)},±√{square root over (⅓)}}
n.sub.CT={−√{square root over (⅔)},0,±√{square root over (⅓)}}
apex={0,0,±√{square root over (⅓)}}
where the ± symbols are positive for the upper cubic tetrahedron, and negative for the lower one.

(28) FIG. 9 shows a triangular patch 905 being projected onto a face of a cubic tetrahedron 910, where the projection covers the set of face cells 915, wherein the projection process is disclosed in for example Chapter 5 of Ashdown, I. 1994. Radiosity: A Programmer's Perspective. New York, N.Y.: Wiley.

(29) Referring to FIG. 5, the differential solid angle dω is mathematically and conceptually equivalent to a ray that extends from the point 505 to a point on the surface of patch 510. Assuming that the patch p has constant radiant exitance M.sub.p and emits, transmits or reflects radiant flux in an ideal diffuse manner, the radiance L.sub.p of the patch will be M.sub.p/π. Further, the luminance of the ray anywhere along its length will also be L.sub.p. From this, it is evident that Equation 1 can be solved by summing the radiance of all rays with origins at point 505.

(30) More practically, an approximate solution to Equation 1 is:

(31) $\begin{array}{cc}{E}_{e,o}\approx \frac{1}{\pi }\underset{n}{.Math.}{M}_p\Delta \Omega & \left(11\right)\end{array}$
where n is the number of rays “cast” from point 505 into the environment in regular or random directions and M.sub.p is the radiant exitance of the patch intersected by the ray. Assuming a sufficient number of rays, ΔΩ is equal to the area of a unit sphere (4π) divided by n.

(32) A more accurate determination of the finite area ΔΩ can be determined by performing a Voronoi tessellation of the unit sphere using the ray-sphere intersections as seeds and calculating the area of each Voronoi cell. The Voronoi tessellation may be determined using for example the algorithm disclosed in Na, H.-S., C.-N. Lee, and O. Cheong. 2001. “Voronoi Diagrams on the Sphere,” Computational Geometry: Theory and Applications 23(2):183-194.

(33) A method of predicting the three-dimensional spherical irradiance field in an interior environment such as a complex architectural space is shown in FIG. 10. In step 1005, a three-dimensional computer-aided drafting (CAD) model of the interior environment is prepared, including ultraviolet reflectance values for each material and radiant intensity distributions for each ultraviolet luminaire positioned within the interior environment. The model exists as a virtual environment that represents the interior environment.

(34) In step 1010, a radiosity solution for the virtual environment is calculated, as described in detail in Ashdown (1994). This will produce an ultraviolet radiant exitance value for each patch of the meshed surfaces of the model. (As will be understood by one skilled in the art, other global illumination techniques such as photon mapping may also be used to determine the spatial distribution of radiant exitance over the set of surface patches.)

(35) In step 1015, a three-dimensional array of spherical irradiance meters is positioned within the virtual environment, with at least one within each volume of interest for air flow calculations using CFD techniques. The array may be regular or partially random for example, or the spherical irradiance meters may be positioned randomly. The volumes of interest may be, for example, those making up an upper area of a room, i.e. above about 2.1-2.5 m. In other cases, the volumes of interest may cover a portion of an upper area of a room. In yet other cases the volumes of interest may include the whole of the space in the interior environment.

(36) In step 1020, a spherical irradiance meter is selected from the set of one or more arrays.

(37) In step 1025, the direct spherical irradiance due to all of the ultraviolet luminaires either fully or partially visible to the spherical irradiance meter is calculated in accordance with the inverse square law (Equation 3). The result is a total direct spherical irradiance for the selected irradiance meter.

(38) In step 1030, the indirect spherical irradiance due to all surface patches fully or partially visible to the spherical irradiance meter is calculated using numerical integration with a light buffer as disclosed herein (Equation 9) and described in further detail in Ashdown (1994). Alternatively, the indirect spherical irradiance can be calculated using the ray casting technique (Equation 11). The result is a total indirect spherical irradiance for the selected irradiance meter.

(39) In step 1032, the total direct spherical irradiance and total indirect spherical irradiance for the selected irradiance meter are added together to result in a total spherical irradiance for the selected irradiance meter.

(40) In step 1035, the spherical irradiance meter is marked as having been processed.

(41) In step 1040, the set of spherical irradiance meters is searched for unmarked meters. If one is found, control returns to step 1020; otherwise control proceeds to step 1045. In other embodiments, there may be other ways of selecting the next spherical irradiance meter, for example they may be numbered.

(42) In step 1045, the three-dimensional spherical irradiance field is trilinearly or tricubically interpolated from the virtual spherical irradiance meter measurements in accordance with the spatial resolution requirements of a CFD program. In other embodiments, other methods of interpolation may be used, where interpolation includes any way of estimating new data points between or within a discrete set of known data points.

(43) In step 1050, the CFD program is used to calculate the total radiant fluence (that is, the sum of spherical irradiances multiplied by time step) for each finite volume of air in the volume of interest. The CFD calculation divides the volume of interest into small, elemental volumes in order to perform its calculation. The total radiant fluence is used to calculate the log.sub.10 reduction in the pathogen species of interest within each of the finite volumes of air. The log.sub.10 reductions may then be combined for all volumes of interest in the virtual environment to result in a global log.sub.10 reduction for the virtual environment. Different pathogens have different susceptibilities to the total radiant fluence (i.e. fluence response values), and so the log.sub.10 reduction will be dependent on both the total radiant fluence and the type of pathogen.

(44) In another embodiment, a different spectral range of optical radiation may be employed, including broadband ultraviolet radiation from pulsed xenon lamps (e.g., Nerandzic, M. M., et al. 2015. “Evaluation of a Pulsed Xenon Ultraviolet Disinfection System for Reduction of Healthcare-Associated Pathogens in Hospital Rooms,” Infection Control & Hospital Epidemiology 36(2):192-197), visible light (e.g., Gillespie, J. B., et al. 2017. “Efficacy of Pulsed 405-nm Light-Emitting Diodes for Antimicrobial Photodynamic Inactivation: Effects of Intensity, Frequency, and Duty Cycle,” Photomedicine and Laser Surgery 35(3):150-156), and infrared radiation from infrared heaters.

(45) In yet another embodiment, the medium may be a substantially transparent and optically non-diffusing fluid such as water, wherein ultraviolet light is employed for example to disinfect flowing water (e.g., Bolton, J. R., and C. A. Cotton. 2008. Ultraviolet Disinfection Handbook. Denver, Colo.: American Water Works Association), Water absorbs UV-C radiation, and so the Beer-Lambert law relating the attenuation of optical radiation to the distance it travels through an absorbing medium must be included in the calculation of both direct and indirect spherical irradiance. The distance from the spherical irradiance meter to each visible luminaire and each surface patch is known, and so this is easily incorporated into the calculations. For example, the inverse square law for direct spherical irradiance (Equation 3) become:
E.sub.o=I.sub.pe.sup.−μr/r.sup.2  (12)
where μ is the absorption coefficient of the fluid for the specified wavelength range. (This coefficient is close to zero for pure distilled water, but is typically 0.01 to 0.2 cm.sup.−1 for municipal drinking water due to dissolved iron, nitrate, manganese and organic matter.)

(46) Similarly, for indirect spherical irradiance, Equation 9 becomes:
ΔE.sub.i≈M.sub.eΔA.sub.me.sup.−μd/π(u.sub.n.sup.2+V.sub.n.sup.2+1).sup.3/2  (13)
where d is the distance from the spherical irradiance meter to the nearest patch in the direction of the cubic tetrahedron face cell. (The distance d is available from the patch projection technique described in Ashdown 1994.)

(47) Referring to FIG. 11, a flowchart is shown of various applications of the method for predicting spherical irradiance. In step 1105, the interior environment is designed, and then modeled in step 1110 using a CAD program. In step 1115, the UV fluence is predicted, according to the method in FIG. 10, for example. In step 1120, a determination is made as to whether there will be sufficient UV fluence in the interior environment to cause the required reduction in pathogen density. If not, the design of the interior environment is adjusted in step 1125, for example by adding in more UV luminaires, by adjusting the placement of one or more of the UV luminaires, or by modifying the UV reflectivity of one or more of the surfaces in the interior environment. After adjusting the design, the interior environment model is rerun in step 1110. If there is sufficient UV fluence to kill the required amount of pathogens, then the process moves on. An action that may be taken may be, for example, building the environment, in step 1130, or modifying the construction of an existing environment. Other actions that may be taken may include the installation of UV luminaires, in step 1135, in an existing building, or upgrading existing UV luminaires in the building. In other cases, the UV luminaires may be adjusted to a different setting, or newly installed luminaires may be set, in step 1140. In step 1145, another possible action is to service one or more of the luminaires, since the output of them may have deteriorated over time, or dust may have accumulated on them, or on the reflective surfaces within the environment that contribute to the UV fluence. In one example, the model may be rerun in step 1110 later in time with downgraded luminaire specifications that take into account the expected deterioration in their performance, which may then result in the resetting of the luminaires in step 1140.

(48) Referring to FIG. 12, a system 1210 is shown for predicting UV spherical irradiance. A processor 1224, which may include multiple constituent processors, is connected to one or more interfaces 1220, 1222 and also to one or more computer-readable memories 1226 storing a program 1230 in the form of computer-readable instructions, which, when executed by the processor, cause the system to follow the processes described in relation to FIGS. 10 and 11. The memory 1226 also stores computer-readable data 1232, which is used by the processor 1224 to run the program 1230. The data 1232 is created by the program 1230, by an external program, or both. Data 1232 includes a model of a virtual environment 1234, which corresponds to the interior environment for which the spherical irradiance is being calculated. Data 1232 includes ultraviolet reflectance values for each material and radiant intensity distributions for each ultraviolet luminaire positioned within the interior environment. User interface 1222 and/or display 1236, allow a user to input or create a model of an environment for which the spherical irradiance is to be calculated, and to receive an output from the system that indicates the spherical irradiance is sufficient to meet desired targets of pathogen reduction, or that indicates the expected log.sub.10 reduction in pathogens in the interior environment. Interface 1220 may be connected to a controller 1240 for one or more luminaires and/or fenestration devices 1244 in the interior environment, which are controlled by the system to provide a desired level of UV fluence at a given time, and/or a desired lighting level for the occupants of the interior environment. The interface 1220 may be connected to the controller 1240 temporarily or permanently, or the connection between the controller and the interface 1220 may be for the transfer of data as and when needed, via one or more other computing devices or non-transient data storage devices.

(49) The system 1210 includes a computer with sufficient calculation and graphics capability to prepare and display the virtual environment and calculate the total radiant fluence for each finite volume of air and the log.sub.10 reduction in the pathogen species of interest.

(50) Referring to FIG. 13, an interior environment 1302 is shown in which there are fenestration devices 1304, visible light luminaires 1306 and UV luminaires 1308. The controller 1240 may be part of a larger control system 1312 connected to multiple other components or systems.

(51) The controller 1240 may read input data from a variety of sensors 1318 and information feeds 1320, the sensors and feeds including at least a plurality of sensors and information feeds from among the class of sensors and information feeds that includes UV photosensors, temperature sensors, occupancy sensors, humidity sensors, timers 1322, personal lighting controls (e.g. interfaces 1324), utility smart power meters 1326, HVAC 1328 and energy storage controllers 1330. The controller 1240 may also calculate the effects of variable building design parameters on building environment characteristics, such as upper air flow, and may output building design parameter setting command signals, in order to maximize energy savings while maintaining selected minimal occupant, including plant or animal, requirements for the building environment characteristics respectively. Minimum occupant requirements may include, for example, a temperature within a preset range, lighting within a preset range, pathogen density below a preset level, etc. The controller reads input data from a variety of sensors 1318 and information feeds 1320, including but not limited to UV photosensors, temperature sensors, occupancy sensors, humidity sensors, timers 1322, personal lighting controls 1324, utility power meters 1326, and other energy management systems 1332, including HVAC 1328 and energy storage controllers 1330. The controller 1240 receives and processes information 1334 about UV light fixtures 1308 and light sources (luminaires) 1306 located in a building's interior environment 1302, including photometric and electrical properties of the luminaires.

(52) In one elementary form, the UVGI control system 1312 also includes: a) at least one controller 1240 that reads input data from a variety of sensors and information feeds, and includes an artificial intelligence engine 1336 or access to an artificial intelligence engine; b) at least one ambient condition sensor 1318 and at least one information feed 1320; and c) an interconnect system 1340 operatively coupling the controller to the sensor and the information feed, configured to provide output data suitable for dimming or switching ultraviolet luminaires 1308 and operating automated fenestration 1304 and other environmental devices.

(53) The controller may further maintain communication with other building automation subsystems, including but not limited to HVAC 1328 and energy storage systems 1330. It may also maintain communication with external systems such as electrical power utilities 1338 and smart power grids.

(54) In an exemplary mode of operation, the controller 1240 continually reads data from the external input and communication devices, calculates the optimal settings for the ultraviolet luminaires 1308, fenestration, and other environmental control devices, and controls those devices to achieve maximal annual energy savings while providing an interior luminous and thermal environment that is consistent with predefined goals and occupant preferences or plant or animal requirements. The “what-if” scenarios capability of the invention deriving from its simulating, using a virtual interior, a building interior environment on a regular basis (for example, hourly), enable a physical UVGI controller system to be designed and programmed accordingly. For example, the design may include an optimal layout of daylight photosensors. The controller may then further access the virtual representation of the interior environment during operation to refine its behavior in response to the building performance by means of “what-if” simulations.

(55) Irradiation of aerosolized particles are possible with well-designed upper air UVGI systems incorporating UV-C, even as standalone controllers that function independently of the HVAC systems and other lighting systems. As such, the controller 1240 may control only UV luminaires in some embodiments.

(56) In a specific embodiment, referring to step 1405 of FIG. 14, a UV irradiance sensor is mounted within each radiation source in the interior environment in order to monitor its radiant flux output. In step 1410, the radiant flux output of each UV luminaire is measured. The measurement may be an absolute measurement or a relative measurement from which an absolute value may be determined. The radiant flux output will decrease over time due to radiation source degradation and dirt accumulation. Periodically, the relative radiant flux outputs of the radiation sources are measured, and the resulting values are used in the model of the interior environment to calculate the spherical irradiance field within the environment, in step 1415. In step 1420, the spherical irradiance field is compared to one or more threshold values of the spherical irradiance field. A threshold may be a minimum threshold that is necessary to provide at least a minimum amount of pathogen reduction in the interior environment. In other cases, a threshold may be a maximum threshold above which the amount of UV in the interior environment is dangerous to the occupants. If no thresholds are met, the process moves to step 1425, after which the radiant flux outputs are remeasured in step 1410. An alarm is issued in step 1430 when the spherical irradiance distribution does not satisfy the minimum and maximum design requirements, i.e. when the spherical irradiance field meets or exceeds a threshold. While the radiation sources may be monitored without using a model of the interior environment, there is no direct relationship between the performance of individual radiation sources and the resultant spherical irradiance distribution. As it is important to know the spherical irradiance distribution, or field, in order to predict the reduction in pathogens, use of the model of the interior environment is essential to calculate the spherical irradiance field. In other embodiments, the predicted reduction in pathogen is recalculated, and the comparison is made with a threshold level of pathogen reduction.

(57) While the description has been given mainly in relation to UV spherical irradiance prediction, the invention also may be used for the prediction of infrared spherical irradiance.

(58) Predicted reductions in pathogen levels are approximate, and desired values of such are to be understood to include tolerances that are normal in the field.

(59) The embodiments of the invention may be varied in many ways. Such variations are not to be regarded as a departure from the scope of the invention, and all such modifications as would be obvious to one skilled in the art are intended to be included within the scope of the claims.