Towards Energy-Positive Buildings through a Quality-Matched Energy Flow Strategy

: Current strategies for net-zero buildings favor envelopes with minimized aperture ratios and limiting of solar gains through reduced glazing transmittance and emissivity. This load-reduction approach precludes strategies that maximize on-site collection of solar energy, which could increase opportunities for net-zero electricity projects. To better leverage solar resources, a whole-building strategy is proposed, referred to as “Quality-Matched Energy Flows” (or Q-MEF): capturing, transforming, buffering, and transferring irradiance on a building’s envelope—and energy derived from it—into distributed end-uses. A mid-scale commercial building was modeled in three climates with a novel Building-Integrated, Transparent, Concentrating Photovoltaic and Thermal fenestration technology (BITCoPT), thermal storage and circulation at three temperature ranges, adsorption chillers, and auxiliary heat pumps. BITCoPT generated electricity and collected thermal energy at high efﬁciencies while transmitting diffuse light and mitigating excess gains and illuminance. The balance of systems satisﬁed cooling and heating demands. Relative to baselines with similar glazing ratios, net electricity use decreased 71% in a continental climate and 100% or more in hot-arid and subtropical-moderate climates. Total EUI decreased 35%, 83%, and 52%, and peak purchased electrical demands decreased up to 6%, 32%, and 20%, respectively (with no provisions for on-site electrical storage). Decreases in utility services costs were also noted. These results suggest that with further development of electriﬁcation the Q-MEF strategy could contribute to energy-positive behavior for projects with similar typology and climate proﬁles.


Introduction
Measures employed in the development of net-zero and net-generating building projects are determined through a design framework, whether implicit or explicit, and any such framework must address multiple criteria to develop traction and ultimately enjoy uptake. These criteria are occupant-related, physics-based, engineering and manufacturing, logistical, and economic in nature. One such framework, termed Quality-Matched Energy Flows (Q-MEF), is explored in this study, and focuses on maintaining and applying the inherent value of climatic energy resources as they pertain to the service demands of the built environment-namely, lighting, electricity supply, and heating and cooling. Specifically considering on-site solar resources, Q-MEF is explored through simulation with a specific generalized building type, in multiple climates. Measures proposed through the Q-MEF rubric comprise a novel Building-Integrated, Transparent, Concentrating Photovoltaic and through the Q-MEF rubric comprise a novel Building-Integrated, Transparent, Concentrating Photovoltaic and Thermal fenestration system (BITCoPT) coupled to an interconnected network of thermal storage and heat pumps (both thermally and electrically powered) which are distributed throughout the building. BITCoPT (investigated in earlier works through modeling, simulation, and experimentation [1][2][3]) uniquely enables the Q-MEF strategy by providing daylighting (illumination and glare control) in addition to generating electricity and collecting thermal energy from solar energy incidents on the building envelope. As Q-MEF and BITCoPT are architecturally grounded, they draw on multiple realms of expertise, and the context of these realms-and the state of their technology-is explored.

Addressing Limits of Load-Reduction Strategies with Active Envelope and Thermal Distribution
Traditionally, net-zero energy building (NZEB) strategies follow a three-tier approach to energy efficiency: load reduction is foremost, followed by the application of passive systems such as shading, and lastly the use of active systems such as energy harvesting [4]. Although load reduction is important, in contemporary architecture this approach begets façade designs that supply good daylighting and natural ventilation, but do not attempt energy harvesting [5,6]. The proliferation of on-site NZEB has been stalled by these legacy design frameworks, which remain grounded in efficient use of grid-supplied electricity and fuel, and employ envelope designs that offload solar energy gains through insulating, limiting apertures, shading, and glazing emissivity control. Moreover, prevalent building-integrated energy harvesting systems (such as roof-top and façadeintegrated PV) exhibit sub-optimal exergetic efficiency and limited installable area, limiting the addressable fraction of a building's service demands, even though in many scenarios solar exergy available to a building is far greater than its own consumption.
Here, we propose an alternative approach, manipulating incident solar energy into daylighting (with minimal glare production) by intercepting that energy with an optically-concentrating photovoltaic fenestration system and hydronic mechanisms that capture and transfer "waste" heat. In this approach, solar gain is re-cast from a low-grade "load" to be remediated through energy intensive heating, ventilation, and air-conditioning (HVAC) systems into a high-grade (thermodynamically useful) energy source for cooling (and other) systems. This approach (Q-MEF) suggests re-distributing a building's systems out across its envelope to more-efficiently interface with the distributed solar resource. This distribution represents a revision of the traditional system topology that has evolved around grid-sourced energy supplies, of centralized plants and trunk-and-branch systems ( Figure 1).

Figure 1.
Transitioning from centralized (grid-source optimized) systems to distributed, integrated systems (climate-sourced optimized), matching resource (sunlight) to service demands (daylighting, power, cooling, heating). If active control is thus shifted outward, the building envelope is empowered to redistribute the abundant light, heat, and power resources available, transferring solar energy towards the building's particular demands and, when necessary, transforming incident solar energy into consumable forms. In Q-MEF, this shift is implemented through integrating distributed coolth storage, (thermally driven) adsorption chillers, and a transparent, optically concentrating photovoltaic and thermal collecting fenestration system (BITCoPT). BITCoPT transmits diffuse irradiance for daylighting and views, and strips out direct normal irradiance, either transforming insolation into electricity or collecting it as thermal energy. Rather than exacerbating cooling demands, insolation can therefore be used directly, as electricity, or as thermal energy that supplies small, localized chillers. Local, distributed use of heat is important for avoiding losses from piping that energy to a central plant. The exergetic value of collected thermal energy is sensitive to such temperature losses, and the surface area-mass ratio of the system (potential for losses through insulation) is inherently high compared to centralized distribution networks with larger pipes. Generated heating and cooling capacities are stored as warm and cold water and circulated through the building via thermally massive hydronic circuits, redistributing these resources towards demands for cooling as well as heating and other process energy, maintaining small but useful temperature differences rather than immediately dispersing generated potential into the built environment. In concept, this integrated system of the multifunctional envelope collector and activated thermal mass distribution can uniquely maintain the quality of the solar resource (through daylighting), and its exergetic content, by collecting a high fraction of direct irradiance either as electricity or as thermal energy at elevated temperatures.

Precedents: Active Facades
The Q-MEF strategy builds on precedent and contemporary building-scale integrated strategies and is dependent on multifunctional envelope technologies. As such, a review of precedent active integrated envelopes and facades (AIFs) is relevant. Buildingintegrated technologies have been reviewed in such focus areas as photovoltaics [7], thermal collection and life cycle analysis [8], concentrating technologies [9,10], and architectural integration [11]. To the authors' knowledge, an investigation coupling a highconcentration/daylighting system such as BITCoPT to a full-building scale application and distribution strategy-designed specifically to leverage envelope-wide collection-has not yet been undertaken.

Benefits of Distributed Systems
Beyond deficits in energy harvesting capacity, the drawback of load-reduction NZEB strategies is in the utilization of harvested energy: Although high quality flows of electricity, light, and thermal energy can be developed with envelope-integrated systems, the building's mechanical infrastructure must be designed to apply these flows effectively. Specifically, generated electricity and thermal resources must be managed. DC microgrids can be employed to efficiently transfer envelope-generated electrical power to equipment and lighting, which often requires DC power natively, using available DC-DC conversion [12]. Harvested thermal resources can be distributed with hydronic systems if the network topology is designed to maintain the useful temperatures (and therefore exergy) that can be developed by a concentrating-type envelope collection system.
Traditionally, building mechanical systems are centrally organized as a logical result of the availability of grid electricity and fuels. Current NZEB strategies therefore involve centralized systems as well (save for active shading/daylighting and natural ventilation). Even on-site harvesting of electricity and heat is typically routed through the same wiring and plumbing that are optimized for the centralized resource-generating systems. To take advantage of envelope-generated energy resources, the mechanical processes that use these resources can be scaled down, parallelized, and located at the envelope, distributed in the same fashion as energy resources are harvested. This is especially relevant for thermal energy collected at high temperatures (beyond roughly 25 • C above the comfort range), which endures larger transport losses than electricity. Recent research into distributed ventilation and humidity control [13] and DC microgrids [14] indicates the inherent efficiency and demand-response benefits of co-locating processes with both their sources of power and the demands they service within a building. Analysis of coupling these end uses to a polygenerating envelope such as BITCoPT has not yet been done.
Many technological elements of load-reduction NZEB strategies can be successfully reconfigured to the distributed Q-MEF strategy. Radiant and mixed-mode heating and cooling, for example, is efficient at managing sensible thermal loads when coupled with effective ventilation and latent-controlled mechanisms such as dedicated outdoor air systems (DOAS) [15,16]. In standard practice, the heat and coolth these systems distribute are generated by centralized chillers and boilers. However, they integrate just as well with distributed thermal control and are employed in the Q-MEF framework.

Thermal Energy Storage for Resource-Matching and Dispatchability
In contrast to grid-sourced energy, which is available on-demand, climatic energy resources must be collected when they are available. To make these resources useful, they must be made dispatchable. Therefore, storing energy in a recoverable way is useful. The benefits of energy storage additionally extend to load factor reduction, grid efficiency and stability, primary fuel costs and pollution reduction, equipment longevity, and resiliency [17]. A range of technologies and materials are actively investigated to these ends [18].
On-site fluidic thermal energy storage (TES) at temperatures suitable for heating (over 30 • C) and cooling (under 18 • C) is attractive because of the potential for good round-trip efficiency: Losses can be minimized through proper insulation, and fluid pumps can be used to dispatch thermal power with little power consumption [19]. Although broadly applicable [20], on-site thermal storage can be limited because of space requirements as well as the relative unfamiliarity of Architecture, Engineering, Construction, and Operation (AECO) sector industries. The current commercial solution to this spatial limitation is to deploy ice storage [21]. Research into alternative phase-changing materials and systems that melt at higher temperatures than water [22][23][24] could improve system efficiencies (by reducing required temperature lifts and maintaining optimal temperatures for process efficiencies) and design desirability (by reducing system footprint). For the sake of conceptual clarity, in the present study, water is chosen as the heat transfer fluid and storage media. The media is well-characterized, but the system geometric configuration-with thermal mass distributed throughout the building's floorplate, along the fluid loop-is uniquely configured to distributed generation and, as such, adds a novel dimension to this investigation.
TES at three temperature ranges is important in the Q-MEF strategy, corresponding to the supply, re-cooling, and output circuits of thermally driven chillers. Because commercial buildings have significant cooling loads and electricity is typically used to create cooling, on-site coolth TES can reduce site energy use, as well as peak demands on the grid, by phase-shifting cooling [25].
In addition to the cool temperature range, TES at two higher ranges is important. In support of the envelope-integrated thermal collection, a buffering mechanism is required to match high-temperature collection to the operation of a thermally driven chiller. Hightemperature storage has been deployed for utility-scale thermal power cycles such as solar electric plants [26], but the lower temperatures useful for driving adsorption chillers (above 60 • C) are easier to maintain in the building context.
A mid-range TES (at temperatures above the range of thermal comfort but below 50 • C) is also required and considered. This storage can act as both a sink for heat rejected from the chillers' operation and a source for heating processes, either as pre-heating for service hot water (SHW) demands or as space heating. Additionally, it has been shown that by storing rejected heat until night-time, when the temperature of the atmospheric sink is lower, the rejection process (and overall energy use) can be more efficient [27].

Entropic Efficiency of Thermal Cascade
By defining multiple temperature (or quality) ranges for thermal storage, a cascade develops through which harvested thermal resources can be matched to a building's service demands ( Figure 2). Besides the demands discussed here, other uses such as humidity control through desiccant regeneration have been investigated. This cascading framework has been recognized in other research, such as the Building as Power Plant concept [28], and in studies and optimizations of latent thermal storage in the building context [23,29]. However, to the authors' knowledge, including daylighting as a primary product and analysis through the mechanism of a multifunctional envelope has not been undertaken.
from the chillers' operation and a source for heating processes, either as pre-heating for service hot water (SHW) demands or as space heating. Additionally, it has been shown that by storing rejected heat until night-time, when the temperature of the atmospheric sink is lower, the rejection process (and overall energy use) can be more efficient [27].

Entropic Efficiency of Thermal Cascade
By defining multiple temperature (or quality) ranges for thermal storage, a cascade develops through which harvested thermal resources can be matched to a building's service demands ( Figure 2). Besides the demands discussed here, other uses such as humidity control through desiccant regeneration have been investigated. This cascading framework has been recognized in other research, such as the Building as Power Plant concept [28], and in studies and optimizations of latent thermal storage in the building context [23,29]. However, to the authors' knowledge, including daylighting as a primary product and analysis through the mechanism of a multifunctional envelope has not been undertaken.

Figure 2.
Cascade of thermal potentials for servicing built environment demands. Optimal exergy application performed by cascading heat through and into processes that require different grades (temperatures) of energy. In the current Q-MEF investigation, thermal storage in three discrete ranges (driving, rejection/heating, chilled) governs energy flow between processes.

Valuing Different Forms of Energy
To evaluate Q-MEF as a strategy, a rubric is devised for the value of the different energy flows-lighting, electricity, heat, and coolth-that can be gleaned from solar energy by BITCoPT. The relative values of these different flows must be grounded in reasonable goals for a given context [30]. Cascade of thermal potentials for servicing built environment demands. Optimal exergy application performed by cascading heat through and into processes that require different grades (temperatures) of energy. In the current Q-MEF investigation, thermal storage in three discrete ranges (driving, rejection/heating, chilled) governs energy flow between processes.

Valuing Different Forms of Energy
To evaluate Q-MEF as a strategy, a rubric is devised for the value of the different energy flows-lighting, electricity, heat, and coolth-that can be gleaned from solar energy by BITCoPT. The relative values of these different flows must be grounded in reasonable goals for a given context [30].
Because well-distributed lighting with good color rendering is difficult to reproduce in the built environment, effective daylighting is a priority in a high-performance building. Diffuse sunlight is multi-directional, exhibits good color rendering, and has good luminous efficacy (relative to common tube fluorescents, although LED lamps are targeted to surpass [31,32]). Diffuse sunlight is easier to apply as daylighting than direct sunlight, which is typically strong and narrowly directional, causing overheating and glare through high contrasts and distracting reflections. Of the resources into which direct irradiance can be transformed, electricity is more valuable than thermal energy, as it is easy to transduce into other forms of work, and its exergy content is not directly subject to the Carnot efficiency. DC electricity is objectively more efficient to use than AC power in many commercial building contexts, as most modern electrical equipment is DC-native, and DC-DC conversion efficiency has reached parity with AC-DC rectification. Although DC distribution is not yet common, standards are mature and gaining acceptance for integration of multiple sources, storage mechanisms, and use devices into sub-building scale microgrids [12].
Thermal energy stored as a heat transfer fluid (HTF) at high temperatures (here defined as over 60 • C, in the context of building systems) can be applied to drive processes such as adsorption chilling, so it can both do work and be applied as heat, in sequence. The higher the temperature developed in the working fluid (relative to an available thermal sink such as the atmosphere or earth), the greater the exergy content, if losses incurred in storage and transportation are managed.
In this study, cooling power is more valuable than heating power for two reasons: More cooling is typically required in medium-sized office buildings, and the coefficient of performance (COP) of a reversible heat pump acting as a chiller is lower than that of the same pump in heating mode.

Exergy-Efficient Solar Collection Aligned with Façade Criteria
Successful on-site energy strategies will internalize a large fraction of the solar exergy available to a building by integrating with the multiple objectives of a successful envelopeobjectives such as thermal regulation, daylighting, privacy, and design identity. The method proposed here for gleaning maximum utility from solar energy is represented by the thermal/energy cascade (Figure 2), and, to functionalize it within the Q-MEF framework, a model is incorporated of BITCoPT, which has previously undergone development [33], simulation [34], and prototype characterization [2,3] (Figure 3).
Because well-distributed lighting with good color rendering is difficult to reproduce in the built environment, effective daylighting is a priority in a high-performance building. Diffuse sunlight is multi-directional, exhibits good color rendering, and has good luminous efficacy (relative to common tube fluorescents, although LED lamps are targeted to surpass [31,32]). Diffuse sunlight is easier to apply as daylighting than direct sunlight, which is typically strong and narrowly directional, causing overheating and glare through high contrasts and distracting reflections.
Of the resources into which direct irradiance can be transformed, electricity is more valuable than thermal energy, as it is easy to transduce into other forms of work, and its exergy content is not directly subject to the Carnot efficiency. DC electricity is objectively more efficient to use than AC power in many commercial building contexts, as most modern electrical equipment is DC-native, and DC-DC conversion efficiency has reached parity with AC-DC rectification. Although DC distribution is not yet common, standards are mature and gaining acceptance for integration of multiple sources, storage mechanisms, and use devices into sub-building scale microgrids [12].
Thermal energy stored as a heat transfer fluid (HTF) at high temperatures (here defined as over 60 °C , in the context of building systems) can be applied to drive processes such as adsorption chilling, so it can both do work and be applied as heat, in sequence. The higher the temperature developed in the working fluid (relative to an available thermal sink such as the atmosphere or earth), the greater the exergy content, if losses incurred in storage and transportation are managed.
In this study, cooling power is more valuable than heating power for two reasons: More cooling is typically required in medium-sized office buildings, and the coefficient of performance (COP) of a reversible heat pump acting as a chiller is lower than that of the same pump in heating mode.

Exergy-Efficient Solar Collection Aligned with Faç ade Criteria
Successful on-site energy strategies will internalize a large fraction of the solar exergy available to a building by integrating with the multiple objectives of a successful envelope-objectives such as thermal regulation, daylighting, privacy, and design identity. The method proposed here for gleaning maximum utility from solar energy is represented by the thermal/energy cascade ( Figure 2), and, to functionalize it within the Q-MEF framework, a model is incorporated of BITCoPT, which has previously undergone development [33], simulation [34], and prototype characterization [2,3] (Figure 3).

Comparison to Alternative Building Energy Schema
Existing schema cross over in purpose and methods with Q-MEF. Although wholebuilding polygeneration is typically considered component by component, as evidenced by review work on net-zero strategies [35], the impact of intersectional modeling is recog- nized [36]. Notably, high-performing built works achieve performance through combined strategies such as occupant-driven energy savings and multi-functional envelope structures. The Bullitt Center's canopy exemplifies this, reducing gains on the vertical envelope while increasing functional area for traditional photovoltaics [37]. Buildings such as Power Plants, mentioned earlier, emphasize the efficiency gains and building demand-matched resources from on-site cogeneration [28]. Exergy analysis provides a rubric for identifying appropriate quality matches between resources and demands [38] and for valuing both warming and cooling storage mechanisms according to temperature differences against an ambient state [39]. The Low Exergy (or LowEx) design emphasizes the application of thermal resources to the small-potential thermal demands of the built environment (heating, cooling, humidity control), the efficiency of heat pumps at maintaining small, useful temperature differences, and sequencing heat pumps to obtain larger temperature lifts [40,41]. LowEx has been extended to district scales [42]. Thermally Active Building Systems (TABS) emphasize the utility of hydronic thermal redistribution, activating thermal mass within the occupied regions of a building, and the efficiency benefits of transferring heat across the large surface areas that comprise indoor built spaces [43,44]. A Passive House is specific regarding the reduction of thermal loads, which can limit solar harvesting opportunities, although work is ongoing to integrate the strategy with on-site generation [45]. These frameworks can be employed to design systems that maximize wise use of available energy resources, but they do not explicitly tackle the provision of desirable lighting conditions in the built environment.
More holistically, Integrated Design emphasizes the consideration of an architectural design as a dynamic totality that shifts between states in a knowable fashion. Attention is paid to the thermal properties (particularly diffusivity) of a building design at the material, construction, and space scales [46,47]. Bioenergetic modeling emphasizes the pathways of potential, flow, and transformation for a broad category of resources, distributed sensing and response, and dynamic systems at the boundary between built and natural environments [48]. Furthermore, thermo-economic analysis is ongoing to ascertain the value of exergetic efficiency and on-site coupled generation and energy storage [49][50][51][52][53], and novel interpretive tools are developed towards these ends [54]. An intention of developing the Q-MEF framework is to leverage the insights of these other strategies to provide for a wide reach of building service demands, including quality lighting.

Research Objective: Modeling Q-MEF Building Energy Behavior
Overlapping with these complementary schemas, Q-MEF incorporates similar concepts. Q-MEF's combination of daylighting-primary operation, distributed high temperature buffering and active chilling, and distributed active thermal mass is potentially novel. Hence, this present work investigates Q-MEF's capacity to more fully leverage building-incident solar energy-and contribute to net-zero and energy-positive designthrough integrating BITCoPT with thermal management at multiple temperature ranges and controlled redistribution.
As a platform for this investigation, a representative medium-scale office building with Q-MEF-aligned systems was modeled and tested in diverse climates, by first generating a pair of baseline (BL) building energy models with both moderate and high window-wall ratios (representing a desirable condition in commercial architecture), and then revising the high-ratio configuration to represent a full Q-MEF implementation with experimentally validated transparent active solar façade energy collectors, distributed thermally driven chillers, sequential thermal redistribution, and thermal storage at three temperature ranges. A fourth building energy model was also generated and tested, wherein the geometry of the building was modified to harvest additional solar energy. The Energy Use Intensity (EUI), peak demands, thermal load profiles, and economic benefits (based on energy consumption) that resulted from this matrix of full-building-scale simulations were compared, and salient behaviors were highlighted.

Methods
To evaluate the effects of Q-MEF on a building's energy flows and consumption, a model was constructed by combining (through post-processing) precursor building energy models of a medium-scale commercial office building (generated in the EnergyPlus environment) with an analytical model of an active façade system (BITCoPT). In addition to these two primary precursors, the Q-MEF model includes separate functions representing: hydronic thermal storage and redistribution, thermally driven chillers, electrically driven auxiliary heat pumps, a deep-mullion curtain wall cavity, fan-powered air volume exchanges in that cavity, and reactive (not predictive) controls. The models and functions were integrated through post-processing in a time-step fashion. The block diagram of the Q-MEF model is presented in Figure 4, and is referred to subsequently.
perimentally validated transparent active solar façade energy collectors, distributed thermally driven chillers, sequential thermal redistribution, and thermal storage at three temperature ranges. A fourth building energy model was also generated and tested, wherein the geometry of the building was modified to harvest additional solar energy. The Energy Use Intensity (EUI), peak demands, thermal load profiles, and economic benefits (based on energy consumption) that resulted from this matrix of full-building-scale simulations were compared, and salient behaviors were highlighted.

Methods
To evaluate the effects of Q-MEF on a building's energy flows and consumption, a model was constructed by combining (through post-processing) precursor building energy models of a medium-scale commercial office building (generated in the EnergyPlus environment) with an analytical model of an active façade system (BITCoPT). In addition to these two primary precursors, the Q-MEF model includes separate functions representing: hydronic thermal storage and redistribution, thermally driven chillers, electrically driven auxiliary heat pumps, a deep-mullion curtain wall cavity, fan-powered air volume exchanges in that cavity, and reactive (not predictive) controls. The models and functions were integrated through post-processing in a time-step fashion. The block diagram of the Q-MEF model is presented in Figure 4, and is referred to subsequently. The model was configured in four ways, with two baseline options (low and high glazing ratios) and two options including Q-MEF components, which were both highly glazed, but with either normal (vertical) or outward-tilted fenestration. Three climates were analyzed, for twelve total model configurations ( Figure 5). The model was configured in four ways, with two baseline options (low and high glazing ratios) and two options including Q-MEF components, which were both highly glazed, but with either normal (vertical) or outward-tilted fenestration. Three climates were analyzed, for twelve total model configurations ( Figure 5).

Systems for Thermal Energy Collection, Redistribution and Use
In the Q-MEF rubric, solar energy on a building's glazed envelope areas is either transmitted as daylight, or concentrated and converted to electricity, or collected as thermal energy, where it is then passed through subsystems, applied to do work (such as chilling), and released as space or water heating, or directly out to climatic sinks. The modeled HVAC strategy (shared between all four configurations) is based on a four-pipe distribution that services all the building's conditioned zones and envelope collection areas ( Figure 6). A chilled distribution loop is maintained at sub-ambient temperatures. A rejection-stage (heating) loop was maintained (through system control algorithms) at super-ambient or near-ambient temperatures, depending on the average demands for heating or cooling which varied seasonally.
Two parallel heat pump systems operate between the loops. In the Q-MEF configurations, small-capacity (<10 kW) adsorption chillers are distributed in the region of the building's envelope and driven by collected solar thermal energy, which is stored adjacently in small buffer tanks. Additionally, in all four configurations, electric water-water heat pumps act in parallel between the two loops.

Systems for Thermal Energy Collection, Redistribution and Use
In the Q-MEF rubric, solar energy on a building's glazed envelope areas is either transmitted as daylight, or concentrated and converted to electricity, or collected as thermal energy, where it is then passed through subsystems, applied to do work (such as chilling), and released as space or water heating, or directly out to climatic sinks. The modeled HVAC strategy (shared between all four configurations) is based on a four-pipe distribution that services all the building's conditioned zones and envelope collection areas ( Figure 6). A chilled distribution loop is maintained at sub-ambient temperatures. A rejection-stage (heating) loop was maintained (through system control algorithms) at super-ambient or near-ambient temperatures, depending on the average demands for heating or cooling which varied seasonally.
Two parallel heat pump systems operate between the loops. In the Q-MEF configurations, small-capacity (<10 kW) adsorption chillers are distributed in the region of the building's envelope and driven by collected solar thermal energy, which is stored adjacently in small buffer tanks. Additionally, in all four configurations, electric water-water heat pumps act in parallel between the two loops.
In both the baseline and Q-MEF configurations, heating and cooling in a zone was modeled as transferring energy through baseboard units or chilled beams, from the rejectionrange loop or to the chilled loop, respectively. Sustainability 2022, 14, x FOR PEER REVIEW 10 of 30 In both the baseline and Q-MEF configurations, heating and cooling in a zone was modeled as transferring energy through baseboard units or chilled beams, from the rejection-range loop or to the chilled loop, respectively. In the Q-MEF model configurations, thermal energy collected by BITCoPT was both used to drive adsorption chilling, and as heat. The glazed areas of the precursor BEMs were defined as BITCoPT collectors [34]. In post-processing, each collector fed a thermal buffer tank, which in turn was plumbed to a chiller ("Buffer/Chiller Operation" block in Figure 4). When a buffer's temperature (T buffer ) rose high enough to drive the connected chiller at a reasonable COP (over 0.5, requiring T buffer ≥ 60 • C, depending on rejection loop temperature), the chiller engaged, pumping heat from the chilled loop to the rejection loop. If the chilled loop became too warm, or the rejection loop cooled off too much, auxiliary heat pumps engaged, pumping heat in parallel to the thermally driven chillers.
Pre-heating for service hot water demands and the DOAS air intake were modeled to occur through water-water heat exchangers, subtracting energy from the rejection-range loop balance according to the difference between the supply temperature reported in the BEM and the loop temperature, less a 2.0 • C approach (Figure 6). The quantity of pre-heating was subtracted from boiler usage of the precursor BEM on a time step basis.
By modeling envelope-integrated solar thermal collection with distributed thermal systems and circulating thermal mass, it was possible to represent the application of thermal energy both as a driving force (for adsorption chillers, at high temperatures) and as useful heat (for zone demands and pre-heating), taking advantage of both the exergy and energy value of the collected energy.

Inputs to Q-MEF Model
Inputs to the overall Q-MEF analysis are described in this section, including precursor models (building energy model baselines, the BITCoPT model, and daylighting method) and boundary conditions (climates).

Climates Considered
The Q-MEF model was analyzed in three climates, exploring a range of system response patterns (Table 1, where T db,outd is outdoor dry bulb temperature). New York City (NYC, LaGuardia Airport TMY3 data) is a fluctuating climate: continental and seasonally humid, with varying weather and strong seasonal swings. Although the city is characteristically dense, and sites are typically shaded by adjacent structures, because this study parametrized climate context, not site, no external shading was defined in NYC or elsewhere. Phoenix (PHX) is an arid subtropical desert with a strong solar resource, high mean temperatures, and large diurnal swings. Mountain View (MTV) is in a semi-arid Mediterranean climate with mild temperatures and a significant solar resource.  [56]), and used in two ways: directly, as two baseline configurations; and as precursors for two Q-MEF configurations. Each energy model, whether it was used directly as a baseline or as a precursor to the Q-MEF model, was configured for three distinct climates (Table 1) Table 2. Stepped dimming with centered zone sensor The first baseline configuration (BL40) was a building energy model with no postprocessing, based on the United States Department of Energy's (DOE) Medium Office Commercial Reference, with ASHRAE 90.1 (2013) specifications [57]. Default operation and occupancy schedules were used, including 5-day work weeks throughout the year, with reduced occupancy and thermostat setbacks outside working hours. The window-wall ratio was 40%. BL40 s HVAC systems were assigned per ASHRAE 189.1: hydronic thermal distribution, with zoned chilled beams, baseboard heating, and DOAS. These updates slightly reduced energy use, relative to the DOE reference.
Representing architectural trends towards highly glazed facades, the second baseline configuration (BL95) had a 95% window-wall ratio. The glazing specifications and systems behaviors were held constant, though systems were re-sized (through the auto-sizing algorithms in EnergyPlus). BL95 was the precursor to the Q-MEF configuration.
For the precursor BEM to the Q-MEF + Tilt configuration, the morphology of the Q-MEF BEM was revised to increase the solar gain on the building. At each floor, facades were swung upward at 20 • (see Figure 5). The chosen angle maximized acquirable irradiance by maximizing solar flux through the exterior glazing on a per-area basis, while limiting the shading of a floor by the floor above.

BITCoPT Envelope Cavity Model
A model of the BITCoPT façade collector [34] was integrated into the Q-MEF model, contributing electrical generation, thermal collection, and direct solar gain reduction based on inputs of climate data and envelope orientation. A constant efficiency (η conv,Egen = 0.96) was applied to the model's electrical output to represent the transformation of the collector's variable-voltage DC output to a useful form (constant-voltage DC for tying into to zonelevel microgrids, or AC for tying into building distribution).
The thermal collection output was likewise modified to represent the heat transfer fluid inlet and envelope cavity temperatures in the Q-MEF model, which were distinct from those values in the BITCoPT simulation. Because the cavity temperatures (T cav ) in the Q-MEF model were held lower (T outd,db < T cav < 38 • C, less than 45 • C in the precursor model, to be realistic for mechanical control components), the Q-MEF-modeled array lost more heat, and thermal collection efficiency was lower than in the precursor model.

Daylighting Modeling Method
All baseline and precursor BEMs employed daylighting controls. Daylighting algorithms in EnergyPlus modified interior lighting schedules according to illuminance measured by virtual sensors. In a prior study [58], a threshold was noted in the room depth for which sufficient work-plane illuminance could be achieved when BITCoPT was introduced, but the perimeter zone depth here (4.57 m) was shallower than that threshold.
Precursor BEMs employed in this study did not include skylights, in keeping with the referenced DOE models. Daylighting from skylights would reduce the baseline lighting electrical loads in the core of the third floor (which comprises 20% of the building's floor area) but would likely contribute to net cooling loads due to solar gains. Although it is surmised that the zones' daylighting and thermal circumstances would improve with glazed roof expanses and installed BITCoPT, in this study the roof collection was treated as stand-alone equipment which did not interact (through lighting, or thermally) with the adjacent zones (in contrast to the vertical glazing expanses, for which light and thermal interactions were modeled).

Q-MEF Model Components and Functions
To represent the impacts of the Q-MEF strategy, components were incorporated with the collection of inputs through additional modeled functions.

Thermal Energy Storage Elements Model
Two types of thermal energy storage elements (TES) were modeled: high-range buffer tanks, and distribution loops. One buffer was defined for each solar collector. One cooling loop and one rejection-range loop were defined to service the building. Both TES types were defined to be well-mixed, rectangular water tanks with a fixed height and width. (Although higher-capacity thermal storage media such as PCMs might increase the overall performance of the Q-MEF system, water was modeled in this study for the sake of conceptual clarity.) As indicated in Figure 4 ("Buffer/Chiller . . . " and "Distribution Loops" modules) the TES exchanged heat with the building's systems (controlled) and with the interior environment (uncontrolled, due to losses through their insulated surface areas). Losses were modeled as the heat transferred across the TES boundary (tank walls) assuming insulation of 100 mm of polyisocyanurate foam, and negligible film coefficients on both interior and exterior surfaces, for an effective thermal resistance of R TES,wall = 0.4 W/m 2 -K. The heat lost from the buffers transferred to the curtain wall cavity in which BITCoPT was installed. The heat lost from the distribution loops transferred to the occupied zones. Pumping power required to equalize the temperature between the storage elements in each loop was defined as equivalent to the power used by the equivalent circulation pumps in the contributing building energy model. Buffer TESs were each sized (through iterative testing) to maximize the yearly output of their paired chillers. The distribution loop capacities were sized to minimize whole-building electrical use (Table 3). It can be noted that in an earlier (unreported) configuration of the Q-MEF model, a high-temperature-range distribution loop was considered as well, to simplify the work extraction from that source. However, the energy losses through the system insulation to indoor ambient temperatures were too high to justify, so a distributed configuration and 4-pipe distribution (reported here) were adopted, with paired buffers and adsorption chillers.

Deep-Mullion Cavity Energy Balance
The building envelopes for the baseline and precursor models were modeled natively in EnergyPlus. For the Q-MEF configurations, the deep-mullion curtain wall cavity (into which BITCoPT integrated) was represented by an energy balance in post-processing. Energy transfers across the glazing determined in the precursor BEMs were replaced by equivalent transfers as determined by this energy balance. The function did not account for variations of temperature or fluid movement within (or external to) the cavity control volume, or non-homogeneous masses. Overall conductivity values between the cavity and adjacent environments were constant, whereas these relationships were variable in the precursor BEM. More accurate heat transfers would be expected from a cavity model that included convection, surface emissivity, material diffusivity, and thermal bridging effects, but one-dimensional, steady state assumptions were deemed sufficient to contrast behavior between baseline and Q-MEF model configurations.
The energy balance comprised: the transmittance of direct and diffuse solar energy into and out from the cavity; the thermal and electrical energy generated by BITCoPT; the heat transferred (via conduction and convection) across glazing surfaces with the two adjoining environments; the heat extracted from the cavity by a flushing function; and the heat lost from high-range buffers (Figure 4: "Cavity balance" block). The balance was ∑ m cav c p,cav dT dt where m cav c p refers to the cavity's thermal mass (in kJ/K); G to transmitted solar power (in W); Q to thermal flow (in W); E gen to electrical generation (in W); DN to direct irradiance; Diff to diffuse irradiance; ind and outd to the indoor and outdoor environments; cond to non-irradiation thermal transfer; cav to the cavity; and flush to the cavity flushing function (see Section 2.3.4).

BITCoPT Area-Based Gap Transmittance
The direct irradiance transmitted through the cavity to a building zone was modeled by multiplying the transmitted irradiance reported by the BEM by T gap , a function that represents the collector's area-based transmittance of direct irradiance (the "Gap function" block in Figure 4). When the solar vector is near to normal with the surface of an envelope region that incorporates a collector, insolation passes through the gaps between BITCoPT modules. As the solar vector moves away from envelope-normal and the collector modules track around pitch and yaw axes, the gaps decrease, falling to zero width at an excursion angle determined by the collector's geometry. T gap was defined as where c vert and c hz are vertical and horizontal components of the fractions of envelope area filled by BITCoPT lens modules; φ and λ are the rotations of the BITCoPT modules around their pitch and yaw axes (in radians); and φ full and λ full (radians) are the respective threshold angles where the gaps between lenses decreased to nil when observed parallel to the solar vector. The floor for T gap was 0, when modules are rotated beyond the pitch and yaw thresholds, and no direct irradiance was transmitted. T gap was applied to only direct insolation since, as modeled, BITCoPT does not intercept diffuse insolation.

Cavity Flushing Function
A cavity flushing function was implemented to simulate the removal of heat from the BITCoPT cavity, adjusting cavity temperature to complement the heating or cooling demands of the adjacent zone. The function simulated a fan moving air between the cavity and the environment (Figure 7). Flushing with outdoor air only was modeled, and the function remained decoupled from the building DOAS.
where mcavcp refers to the cavity's thermal mass (in kJ/K); G to transmitted solar power (in W); Q to thermal flow (in W); Egen to electrical generation (in W); DN to direct irradiance; Diff to diffuse irradiance; ind and outd to the indoor and outdoor environments; cond to non-irradiation thermal transfer; cav to the cavity; and flush to the cavity flushing function (see Section 2.3.4).

BITCoPT Area-Based Gap Transmittance
The direct irradiance transmitted through the cavity to a building zone was modeled by multiplying the transmitted irradiance reported by the BEM by Tgap, a function that represents the collector's area-based transmittance of direct irradiance (the "Gap function" block in Figure 4). When the solar vector is near to normal with the surface of an envelope region that incorporates a collector, insolation passes through the gaps between BITCoPT modules. As the solar vector moves away from envelope-normal and the collector modules track around pitch and yaw axes, the gaps decrease, falling to zero width at an excursion angle determined by the collector's geometry. Tgap was defined as where cvert and chz are vertical and horizontal components of the fractions of envelope area filled by BITCoPT lens modules; ϕ and λ are the rotations of the BITCoPT modules around their pitch and yaw axes (in radians); and ϕfull and λfull (radians) are the respective threshold angles where the gaps between lenses decreased to nil when observed parallel to the solar vector. The floor for Tgap was 0, when modules are rotated beyond the pitch and yaw thresholds, and no direct irradiance was transmitted. Tgap was applied to only direct insolation since, as modeled, BITCoPT does not intercept diffuse insolation.

Cavity Flushing Function
A cavity flushing function was implemented to simulate the removal of heat from the BITCoPT cavity, adjusting cavity temperature to complement the heating or cooling demands of the adjacent zone. The function simulated a fan moving air between the cavity and the environment (Figure 7). Flushing with outdoor air only was modeled, and the function remained decoupled from the building DOAS.  Flushing was implemented both in the precursor BEM (as a window:Airflow property in EnergyPlus, with a constant volumetric rate of 0.6 m 3 /s-m and assumption of no fan power) and in post-processing. In post-processing, the maximum allowable cavity temperature was T cav,lim,high = 40 • C, but if zone cooling demand was significant, the target temperature was set to the outdoor (dry bulb) temperature). The post-processing flushing function was modeled as: where Q cav,flush was the resulting thermal flow (in W); K flush was the (non-dimensional) proportional gain tuned for each zone; C cav was the cavity thermal mass (in kJ/K); T cav was the cavity temperature; T cav,target was set to either T outd,db[n −1] or T ind depending on heating demand; and T cav,offset = 2 • C was used to establish a dead band, preventing the system from operating if heat removal would be inconsequential. K flush,IDN was a separate gain used to modify the flush rate according to zone's direct insolation at the current-time step.
Fan power required to flush the cavity was determined assuming a constant pressure head (and therefore a constant power draw) multiplied by an hourly duty cycle. 100% of fan (electrical) power was designated to be taken up by the airflow (as increased temperature). These flushing controls maintained the cavity temperature close to the chosen target, minimizing unwanted non-insolation thermal transfer from the envelope cavity to the building's interior, and contributing energy to under-heated zone conditions. Gains were tuned to optimize net generation at the building meter.

Thermally Driven Chillers Model
Adsorption chillers were modeled by fitting a solution surface to a manufacturer's COP data [59] relative to temperatures of buffers, chilled loop (T chilled ), and rejection loop (T rej ). Upper and lower limits on COP chiller were implemented to keep it between 0 and 0.56. The chillers were controlled to activate when the temperature of their paired buffer exceeded a threshold (T buffer > 60 • C). Chillers drew energy according to the excess in the buffer over the threshold, with no capacity limit and no dependence of COP on the fractional capacity, representing multiple chillers ganged in parallel.

Auxiliary Heat Pumps and Heat Rejection Model
Water-water heat pumps were incorporated in the thermal control strategy to maintain the temperatures of the chilled and rejection-range distribution loops (see Figure 4, "Distribution Loops" block). If the chilled loop became too warm, energy was pumped from it to the rejection-range loop. Pumps are modeled steady state, at a COP varying with temperature lift, with functions sourced from the literature [41], capped at 16.0, with an 0.5 exergetic efficiency.
To prevent the rejection-range loop from getting too warm, heat removal from the loop to the outdoor environment was simulated by modeling a dry cooler, or fan-assisted water-to-air heat exchanger. The heat removal function, Q rej,env (in W) was: where K rej was the (non-dimensional) proportional gain; C loop,rej was the thermal mass of the rejection loop (in kJ/K); T rej is the rejection distribution loop temperature; and T rej,target was the temperature set point. Like with the cavity flush fans, the fan power required to remove heat was determined with a constant power draw and an hourly duty cycle.

Loop Temperature Controls
The temperature of the chilled loop (T chilled ) was maintained by the activity of the adsorption chillers and auxiliary heat pump (when necessary), removing the energy gained from space cooling. If T chilled was ever driven too low (T chilled > 11.4 • C), both the chillers and heat pumps were deactivated in the next time step. If the chilled loop temperature rose too high (T chilled > 17.0 • C, representing over-loaded adsorption capacity) the auxiliary heat pump activated in the next time step.
The rejection-range loop temperature (T rej ) was controlled by the activity of the auxiliary boiler, the auxiliary heat pump, and the heat rejection system. A target temperature T rej,target was set, either to T out,db (if heating demands were expected to be low) or to T rej = 27 • C. If T rej dropped below T rej,target , the auxiliary heat pump activated at the next time step. If the heat pump was unavailable (due to a low chilled loop temperature), the boiler activated. If T rej was too high, the heat rejection system was called, proportional to the difference between T rej and T rej,target . An additional check of the current heat rejection COP (COP rej > 3.0) prevented operation when the temperature difference was small, to favor night flushing of excess energy in the rejection loop. These controls together maintained loop temperatures within useful bounds, providing cooling and heating to the zones, and heat rejection capacity for the adsorption chillers.

Utility Cost Metrics
Costs for electrical use, capacity, and demand and gas use charges were initially computed within the BEM. Simulation outputs for the baseline configurations were used and scaled on an hourly basis for the Q-MEF configurations. Electricity costs were determined with monthly peak demand charges (USD 17.00/kW for the winter and USD 38.15/kW for the summer) and hourly energy charges (USD 0.125/kWh on-peak and USD 0.105/kWh off-peak), with surplus electricity generation net-metered at 100% of the current rate. The cost for natural gas was USD 1.30 per therm.

Results
Simulating a medium-scale commercial office building according to the Q-MEF framework demonstrated implications for lighting demands, cooling and heating loads, peak demands, net energy use, operational costs, and design considerations.

Daylighting Impact on Lighting Energy Use
Modeled lighting energy use in the high-glazing configurations improved between 13% and 17% over BL40 (Table 4). Daylighting controls were active in all configurations. Because daylighting controls were modeled in the precursor BEMs, the controls did not respond to the reduction of direct irradiance due to BITCoPT. Prior daylighting and glare analyses of BITCoPT [58] determined sufficient illuminance through the same depth as the perimeter zones modeled in this study (4.57 m), and as other equipment (such as blinds) was not modeled to control for over-lit moments, it was determined that the daylighting behavior would translate from the precursor BEMs to the Q-MEF configurations.
Daylighting potentially has a great effect on the energy use profile of a building, but because it is highly contingent on occupant behavior [60] representation in the simulation method employed here is difficult. It is possible that the daylighting energy benefits of BITCoPT are over-predicted, as cloudy moments would cause under-lighting. However, experimentation has suggested that lighting through the system increases with partial clouds, as there is more diffuse light to transmit [3]. It is also possible that benefits are underpredicted due to glare, as, during brighter moments in the baseline buildings, occupants who experience excessive brightness or glare might deploy blinds and electric lighting. More differentiation between the configurations might be noted if active technologies (blinds, BITCoPT) and occupant behavior were modeled dynamically with the sensors and dimmers.

Heating and Cooling Loads
The incorporation of Q-MEF systems resulted in various responses in heating and cooling loads across the climate types, demonstrating complex interdependencies between envelope loads and building demands ( Table 5). Heating loads comprise the sum of modeled baseboard heating and DOAS preheating, while the cooling loads consist only of the modeled chilled beam responses. Due to the preponderance of cooling demands in the moderate MTV climate, the heating increases were low (though not negligible) relative to total demands. The NYC climate showed sensitivity to glazing area and type, as the BL95 and QMEF cases showed marked changes in the total loads. Installation of BITCoPT (the Q-MEF case vs. the BL95 case) resulted in reductions in cooling loads in all climates, although the Q-MEF + Tilt case did not decrease loads further.
The difference in loads occurred mainly in the south, east, and west perimeter zones, where decreased heat gain during direct solar conditions due to BITCoPT caused more frequent net-heating loads. In effect, the Q-MEF configurations resulted in more "skindominated" behavior of the building, where heating and cooling demands correlate to the difference between indoor and outdoor temperatures, as opposed to "core-dominated" behavior, where internal energy gains cause persistent cooling demands throughout the year (shown in the evolving Heating-Cooling ratio, in Table 5).

On-Site Thermal Collection and Application
In addition to driving chilling processes, thermal energy collected by BITCoPT was applied to building demands for zone heating and SHW preheating ( Table 6). The final column in Table 6 is a sum of the zone heating loads, SHW preheating, collected thermal energy (negative sign), and DOAS heating (not described in a separate column). Via the rejection loop, collected thermal energy was applied to zone heating loads and SHW preheating (although with low demand in office-dominated buildings, the latter factor was small). Boiler (and fuel) usage was observed to increase in response to heating demands, despite net thermal collection in some cases exceeding net heating demands. This disparity indicates non-optimal behaviors in the thermal storage mechanisms, including mismatches between the times of collection and demand. This mismatch occurred largely over the annual cycle (Figure 8), indicating the usefulness of ground-source thermal storage, a function which was not implemented in the current Q-MEF model.

Solar Cooling with Adsorption Chilling
Solar thermal energy, collected by BITCoPT, was applied to drive adsorption chillers (Table 7). This cooling was additive to any passive reductions in direct solar gains from the system, with the added benefit of dispatchability since the capacity was stored in the thermally massive, chilled distribution loop.

Solar Cooling with Adsorption Chilling
Solar thermal energy, collected by BITCoPT, was applied to drive adsorption chillers (Table 7). This cooling was additive to any passive reductions in direct solar gains from the system, with the added benefit of dispatchability since the capacity was stored in the thermally massive, chilled distribution loop.
Cooling loads in the Q-MEF configurations were lower than in the BL95 configuration. The work done by adsorption chillers further reduced the cooling required from the auxiliary heat pump. The magnitude of collected thermal energy (Table 6) relative to loads suggests that more systems (such as night-flush controls and ground-source heat exchange) would be useful to perform more controlled storage and release of heat over both diurnal and annual cycles.
Cooling power produced from solar energy was not sufficient to provide 100% of modeled demands. Chilling processes capable of higher COPs might close that gap. This requires higher operating temperatures, and therefore higher solar concentration ratios, to boost exergy collection and offset losses from the thermal collection stage. A higher rejection temperature would also then be allowable, reducing the gap between that target and the target for zone heating (the two services being provided by the common loop). This tradeoff was not explored here, but it is noted that the required insulation of the hydronic system is complex, and at higher cell operating temperatures (roughly 100 • C and above) radiation from the cell becomes a significant thermal loss factor.

Energy Use Profile Comparison
Broadly, Q-MEF simulations improved over the baselines for all observed metricsnet electrical demand, net energy demand, and peak electrical demand (defined as the maximum observed electrical demands between noon and 5 pm during the summer season). There were conditional exceptions: BL40 demonstrated the lowest peak demand of the NYC models, and the demand reductions of QMEF over BL40 were trivial. In both these cases, reductions in electrical EUI were still significant. Results are summarized in Table 8. The Q-MEF + Tilt configuration demonstrated the highest generation and lowest total site EUI across climates. It is notable that electricity EUI reaches net zero in the Phoenix and Mountain View climates, although there remains significant consumption of gas for heating-higher than that of the baselines.
Although it is significant that Q-MEF + Tilt demonstrated net zero electrical use, benefits were incremental (or negative) over the straight Q-MEF configuration, suggesting further analysis comparing the marginal utility of the energy benefits with the marginal costs of increasing the complexity of the building's design. It's notable that in the Phoenix condition, the Q-MEF + Tilt configuration showed higher peak demand than the straight Q-MEF configuration, differing from the other two climates. This is due to an attenuation of the cavity flush function from the combination of elevated outdoor temperatures and the increased available insolation.

Utility Cost Analysis
Annual energy costs were calculated for all configurations according to rates for electricity demand (kW), supply (kWh) and natural gas supply (therms or kWh) (Figure 9, Figure 10, Table 10). Q-MEF configurations reduced total use costs, while demand costs were generally lower than costs in BL95, and similar to costs in BL40.  Table 10). Q-MEF configurations reduced total use costs, while demand costs were generally lower than costs in BL95, and similar to costs in BL40.    Observed in Table 8, Q-MEF configurations showed increased natural gas consumption, but due to lower electricity use, the total EUI was reduced in all Q-MEF cases, relative to BL40 and BL95. The analysis of costs demonstrates that (due to the low cost of natural gas compared to electricity) Q-MEF configurations showed annual reductions (Table 10).  HVAC, which incurs additional effects on a building project's value, such as reduced inhabitable room height due to the depth of duct work. The overall cost analysis of Q-MEF is highly contextual, due to the interaction of these factors and additional localized factors, such as capitalization rates expected on monies obtained to finance a project, which is an in-part function of a perceived risk. A detailed cost analysis is therefore usefully done at the scale of individual projects, or for broader applicability, by parameterizing these factors and undertaking the resulting matrix of sensitivity analyses. That breadth of analysis is outside the bounds of this study, which, for the sake of its own broader applicability, considered generalized circumstances. The present results, however, are a precursor to such a technoeconomic analysis, which merits further investigation.

Discussion
The modeled application of the specific systems described in this study according to the Q-MEF strategy resulted, in simulation, in significant reductions in energy use in the three modeled climates. On-site net-zero electricity was demonstrated in two of the climates-a significant result for buildings in the modeled size class. The results further suggest that additional modifications might show additional energy-use benefits: the baseline configurations were designed according to current minimum efficiency codes but did not incorporate the full gamut of currently available high-performance building strategies, such as DOAS enthalpy recovery (important in high-humidity climates such as New York) or ground-source heat exchange (particularly useful in climates with steady cycles such as Phoenix). Pursuing the current industrial interest in full electrification might also reveal further reductions, but as the modeling process dictates systems similarity between proposed and baseline models (which employ gas heating), these impacts were not represented. The application of these and other strategies would decrease the baselines' EUI, and therefore the EUI of the Q-MEF implementations as well.
Not all high-performance design strategies are synergistic with Q-MEF. External shading devices, for example, are a passive-design strategy to reduce fenestration gains. Q-MEF attempts to internalize these gains, which increases cooling loads, particularly over the BL40 baseline, but overall reduces both site energy use (all climates) and peak demands (in Phoenix, Table 8). Although counter-intuitive from the passive design perspective, these results suggest the benefits of engaging the solar resource.
A primary benefit to coolth storage systems is peak demand reduction, which commercial systems bank on for their value proposition [16]. The demand reductions demonstrated in this study, though significant, may not fully realize the benefits of this storage. The gap may be due to the simplified controls in the Q-MEF model, which were chosen in part for compatibility with the post-processing modeling method. Optimizing controls has been shown to benefit thermal storage applications [61] and might improve the utilization of the modeled storage.
A general shift from cooling loads to heating loads was observed with Q-MEF vs. baseline configurations (Table 5), which is in keeping with the Q-MEF concept of utilizing available solar resources, rather than mitigating and compensating for them. In the Mountain View climate, the cooling-dominated baselines were flipped to more heating operation. This flexibility around the balanced point suggests that introducing ground-source heat exchange, which benefits from that neutrality, might complement the Q-MEF strategy. Ground-sourcing would be worth investigating in the New York climate as well, which shifted from balanced to heating-dominated, as such heat pump-enabled strategies are intrinsically more effective at heating than cooling. Ground-source systems likely could not "keep up" with the required quantity of heat rejection in the more extreme Phoenix climate, suggesting other strategies might be useful, such as radiative night-sky sinking, which has recently advanced through material investigations [62].
In this study, the thermal collection efficiency of BITCoPT averaged lower than in precedent studies. This stems from conflicting demands on the envelope cavity: sensitivity analysis performed through multiple simulation runs showed that allowing elevated cavity temperatures reduced transport losses in BITCoPT, but increased perimeter zone cooling demands and overall net energy use. This trade-off reinforces how a building with Q-MEF is a coupled system of components that experience unique forcing functions. Optimizing overall objectives in such a system requires subordinating the peak performance of specific sub-systems.

Conclusions
To test the integration of technologies through an architectural, whole-system approach to design, the production of benefits from on-site solar resources was explored through simulations with a parametric group of building energy models. The models were assembled from precursors: a set of building energy models, a model of an active envelope technology, and representations of adsorption chillers, hydronic thermal distribution and storage elements, water-to-water heat pumps, and ancillary systems. The overall strategy was described as a quality-matched energy flow (Q-MEF) network. Simulations of a 5000 m 2 , three-floor office building demonstrated reductions in electrical use over 70% from the baseline in a humid-continental climate (New York City), and on-site net zero electricity in arid subtropical (Phoenix) and semi-arid Mediterranean (Mountain View) climates. Peak purchased electrical demands decreased up to 6%, 32%, and 20% respectively. The magnitude of these results suggests the usefulness of deeply integrating the multifunctional envelope technology with the balance of a building's systems that process and distribute collected thermal energy. Demand for (purchased) electricity remained significant, indicating the potential for future investigation of on-site electrical storage, in addition to the modeled thermal storage. Cost outcomes were reported, with reductions (in summed energy and demand charges) of 35% in New York City, 64% in Phoenix, and 66% in Mountain View, relative to the highly glazed baseline model configurations. A shift was noted, from demand for cooling in the baseline configurations to demand for heating in the Q-MEF configurations. This shift indicates potential benefits from additional thermal storage technologies beyond those modeled in this study (such as ground-source heat exchange), and, overall, the possibility of energy-positive performance for this common class of buildings, in a range of climates.