Experimental Validation of Water Flow Glazing: Transient Response in Real Test Rooms

The extensive use of glass in modern architecture has increased the heating and cooling loads in buildings. Recent studies have presented water flow glazing (WFG) envelopes as an alternative building energy management system to reduce energy consumption and improve thermal comfort in buildings. Currently, commercial software for thermal simulation does not include WFG as a façade material. This article aims to validate a new building simulation tool developed by the authors. Simulation results were compared with real data from a scale prototype composed of two twin cabins with different glazing envelopes: a Reference double glazing with solar-control coating and a triple water flow glazing. The results showed a good agreement between the simulation and the real data from the prototype. The mean percentage error of the indoor temperature cabin was lower than 5.5% and 3.2% in the WFG cabin and in the Reference glazing one, respectively. The indoor air temperature of the WFG cabin was 5 ◦C lower than the Reference one in a free-floating temperature regime when the outdoor air temperature was 35 ◦C and the maximum value of solar radiation was above 700 W/m2. WFG has energy-saving potential and is worthy of further research into the standardization of its manufacturing process and its ability to increase building occupants’ comfort.


Introduction
Energy consumption in buildings shows a growing trend worldwide and is of primary concern for the world population [1]. Over the last decade, nearly 60% of total net electricity consumption in Organization for Economic Co-operation and Development (OECD) economies, was in the building sector, both residential and commercial [2]. The residential building sector is responsible for more than half of the electricity consumption in developing countries [3].
In the frame of the Paris agreement in 2015, 195 countries adopted 17 sustainable development goals (SDGs) as the outcome of the UN's inclusive and comprehensive negotiations in the frame of the 2030 agenda [4]. The seventh goal states that using clean and sustainable energy sources is an opportunity to transform economies and lives, especially in developing countries.
Annual power consumption depends on the use of the building, construction year, number of floors, building structure, and building location [5]. When it comes to heating and cooling consumption, the heating, ventilation, and air conditioning (HVAC) system, exterior walls, and glazing are the

Materials and Methods
This section aims to provide a simplified model that helps designers to analyze the energy strategy of a dynamic envelope. Commercial BES tools do not include WFG as an option, so it is necessary to validate the hypothesis with data from real prototypes. The first subsection set the criteria to select the spectral and thermal parameters of the WFG. The second subsection described geometry, energy management, and the materials used in the prototype.

Simplified Model of Triple WFG
Water flow glazing (WFG) allows the flow of water between two glass panes. Water captures the solar NIR and increases its temperature through the window. The flow of the water enables the homogenization of the building envelope temperature so that designers can apply energy-saving strategies, such as energy storage or solar energy rejection. Table 1 shows the list of symbols that have been used in equations. Heat fluxes through the different layers of the glazing i 0 Normal incident solar irradiance (W/m 2 ) g OFF Solar heat gain coefficient without flow rate. g ON Solar heat gain coefficient with high flow rate. θ i Interior temperature (K) θ e Exterior temperature (K) θ j Temperature of the glass layer (K) θ IN Inlet temperature of the water chamber (K) θ OUT Outlet temperature of the water chamber (K) θ w Temperature of the water (K) θ STAGNATION Temperature of the water whenṁ = 0 (K) U Thermal transmittance (W/m 2 K) U i Interior thermal transmittance (W/m 2 K) U e Exterior thermal transmittance (W/m 2 K) U w Thermal transmittance (water chamber-interior) (K) T Transmittance of the glazing R Reflectance of the glazinġ m Mass flow rate (kg/s m 2 ) c Specific heat of the water (J/Kg K) P Heat gain in the water chamber (W) Figure 1 shows the heat flux and the temperature distribution when the outdoor temperature (θ e ) is higher than the indoor (θ i ) through a triple WFG. The thermal transmittance, U, is the parameter that explains the heat transfer. The European Standard determines its value and a calculation method [58,59]. Previous studies explained the thermal and spectral properties of WFG and its behavior [60,61]. This research used a simplified set of equations from those studies, along with data from commercial software, to assess the performance of the prototype defined in Section 2.2. Equations (1) to (10) show the heat fluxes through the different layers of the glazing. They consider the energy balance at each layer and the Newton's definition of heat flux.
= ℎ ( − ), = , = ℎ ( − ), = ℎ ( − ), = ℎ ( − ), The U values depend on the interior heat transfer coefficient, hi, the exterior heat transfer coefficient, he, the air chamber heat transfer coefficient, hg, and the water heat transfer coefficient. Previous studies explained the thermal and spectral properties of WFG and its behavior [60,61]. This research used a simplified set of equations from those studies, along with data from commercial software, to assess the performance of the prototype defined in Section 2.2. Equations (1) to (10) show the heat fluxes through the different layers of the glazing. They consider the energy balance at each layer and the Newton's definition of heat flux.
The U values depend on the interior heat transfer coefficient, h i , the exterior heat transfer coefficient, h e , the air chamber heat transfer coefficient, h g , and the water heat transfer coefficient.
Sustainability 2020, 12, 5734 6 of 23 A 1 is the absorptance of the exterior glass pane, A 2 , is the absorptance of the middle glass pane, and A 3 is the absorptance of the interior one. A w is the absorptance of the water chamber. The absorptance, A v , depends on the energy absorbed by the glass panes and by the water: Solving the Equations (1) to (10) and using the values of Equations (11) to (13): Water heat gain is a power magnitude, and it is measured in watts (W). Equation (2) shows the analytical expression.
where P is the power absorbed by the water, θ OUT and θ IN the temperature of water leaving and entering the glazing, respectively,ṁ is the mass flow rate, and c is the specific heat of the water. The mass flow rate is a measurement of the amount of mass passing by a point over time. The goal of absorbing the same power, P, can be achieved with a highṁ, which results in a low-temperature increase or a lowṁ, which results in a high-temperature difference between the inlet and outlet. Equation (16) results by combining Equations (14) and (15).
where θ e and θ i are, respectively, the temperature outdoors and indoors of the room. Analyzing the Equation (16), the power absorbed by the water varies with θ IN . Considering the rest of the equation as a constant, P 0 , P linearly decreases with slope (U i + U e ). Figure 2 shows that at a specific value ofṁ, there is a maximum absorbed power, P, depending on θ IN .
If boundary conditions do not change with time, the solution becomes constant when the system reaches a steady state. In this set of test cases, sunrays are perpendicular to the glazing. This assumption eliminates uncertainties associated with the dependence of each layer's absorptance with the angle of incidence. Two case studies are considered: (a) and (b). Table 2 defines the outdoor and indoor boundary conditions. Case (a) studies the influence of θ IN with a fixed absorptance, and case (b) considers the impact of the absorptance when θ IN is fixed.   Figure 2 illustrates the power variations with water flow rate,ṁ. There is a maximum water heat gain whenṁc >> Ue + Ui. The maximum power absorption occurs when A v is high and θ IN is low. If the goal is to reject the solar energy, A v must be as low as possible, with solar control coatings in the outermost glass panes. In this case (ṁc >> Ue + Ui), the absorbed power (P) is the maximum power (P max ) Sustainability 2020, 12, x 7 of 24 the goal is to reject the solar energy, Av must be as low as possible, with solar control coatings in the outermost glass panes. In this case (ṁc >> Ue + Ui), the absorbed power (P) is the maximum power (Pmax) Combining Equations (14) and (19), the value of θOUT is: θSTAGNATION is the temperature of the water when the mass flow rate is stopped. Figure 3 illustrates the relationship between θOUT and ṁ in two cases: (a) and (b). Case (a) shows that θSTAGNATION is the same if the water absorption Av does not change; case (b) shows the variations of θSTAGNATION with different conditions of Av. When ṁ is close to zero, then θOUT gets the maximum value, which is θSTAGNATION. When ṁ is very high, then θOUT = θIN. Equation (21) shows the total heat flux, q: Combining Equations (14) and (19), the value of θ OUT is: θ STAGNATION is the temperature of the water when the mass flow rate is stopped. Figure 3 illustrates the relationship between θ OUT andṁ in two cases: (a) and (b). Case (a) shows that θ STAGNATION is the same if the water absorption A v does not change; case (b) shows the variations of θ STAGNATION with different conditions of A v . Whenṁ is close to zero, then θ OUT gets the maximum value, which is θ STAGNATION . Whenṁ is very high, then θ OUT = θ IN .
Sustainability 2020, 12, x 7 of 24 the goal is to reject the solar energy, Av must be as low as possible, with solar control coatings in the outermost glass panes. In this case (ṁc >> Ue + Ui), the absorbed power (P) is the maximum power (Pmax) Combining Equations (14) and (19), the value of θOUT is: θSTAGNATION is the temperature of the water when the mass flow rate is stopped. Figure 3 illustrates the relationship between θOUT and ṁ in two cases: (a) and (b). Case (a) shows that θSTAGNATION is the same if the water absorption Av does not change; case (b) shows the variations of θSTAGNATION with different conditions of Av. When ṁ is close to zero, then θOUT gets the maximum value, which is θSTAGNATION. When ṁ is very high, then θOUT = θIN. Equation (21) shows the total heat flux, q: Equation (21) shows the total heat flux, q: where g describes the proportion of solar energy transmitted indoors, U(θ e −θ i ) expresses the heat exchange between the room and outdoors, and U w (θ i −θ IN ) represents the heat exchange between the water chamber and indoors. As per Equations (34) and (35) in [50]: Since the WFG is a dynamic envelope, the g factor depends on the flow rate. When the water flows, the g factor decreases, and when the water flow stops, the solar energy enters the building, and the g factor increases. The thermal transmittance, U, is almost zero when the flow rate is the design flow rate because the water chamber isolates the building from outdoor conditions. When the water flow stops, the thermal transmittance depends mainly on the components that make up an insulated glass unit: the glass panes, coatings, spacer, sealing, and the gas filling the sealed space.

Description of the Prototype
A prototype has been built to compare the thermal performance of WFG with a double-glazing solution. It is located near Madrid, Spain (latitude 41 • 22 6" N, longitude 3 • 29 57" W, altitude 1037 MAMSL). Conceptually, this mock-up is a mobile and autonomous prototype consisting of two cabins named WFG and Reference. The prototype design allows both the primary and secondary circuits to be housed within the demonstrator structure, although in independent sectors. Besides, the mock-up has four wheels at the bottom of the structure for easy transport and orientation.
The prototype has two levels: the lower level of 500 mm high and the upper level of 1000 mm high. The lower level holds the primary circuit, which is composed of a "Peltier" unit, and a lithium battery that feeds it. The upper level includes the two cabins. Finally, a photovoltaic panel has been installed on the top roof allowing to store electric energy in the battery for the operation of the cooling device. A circulating device is made up of a 10 W solar pump. The primary circuit connects the circulating system with a "Peltier" device. The secondary circuit goes from the circulating system to the WFG, with a design flow rate of 0.5 L/min. Figure 4 illustrates a schematic of the prototype with its main components, and Figure 5 shows the layout of the prototype. It consists of two semi-detached insulated cabins, one with WFG facing South and the other with a Reference glazing in the same orientation. The dimensions of the windows are 1 m × 0.5 m. where g describes the proportion of solar energy transmitted indoors, U(θe−θi) expresses the heat exchange between the room and outdoors, and Uw (θi−θIN) represents the heat exchange between the water chamber and indoors. As per Equations (34) and (35) in [50]: Since the WFG is a dynamic envelope, the g factor depends on the flow rate. When the water flows, the g factor decreases, and when the water flow stops, the solar energy enters the building, and the g factor increases. The thermal transmittance, U, is almost zero when the flow rate is the design flow rate because the water chamber isolates the building from outdoor conditions. When the water flow stops, the thermal transmittance depends mainly on the components that make up an insulated glass unit: the glass panes, coatings, spacer, sealing, and the gas filling the sealed space.

Description of the Prototype
A prototype has been built to compare the thermal performance of WFG with a double-glazing solution. It is located near Madrid, Spain (latitude 41°22′6″ N, longitude 3°29′57″ W, altitude 1037 MAMSL). Conceptually, this mock-up is a mobile and autonomous prototype consisting of two cabins named WFG and Reference. The prototype design allows both the primary and secondary circuits to be housed within the demonstrator structure, although in independent sectors. Besides, the mock-up has four wheels at the bottom of the structure for easy transport and orientation.
The prototype has two levels: the lower level of 500 mm high and the upper level of 1000 mm high. The lower level holds the primary circuit, which is composed of a "Peltier" unit, and a lithium battery that feeds it. The upper level includes the two cabins. Finally, a photovoltaic panel has been installed on the top roof allowing to store electric energy in the battery for the operation of the cooling device. A circulating device is made up of a 10 W solar pump. The primary circuit connects the circulating system with a "Peltier" device. The secondary circuit goes from the circulating system to the WFG, with a design flow rate of 0.5 L/min. Figure 4 illustrates a schematic of the prototype with its main components, and Figure 5 shows the layout of the prototype. It consists of two semi-detached insulated cabins, one with WFG facing South and the other with a Reference glazing in the same orientation. The dimensions of the windows are 1 m × 0.5 m.  Regarding the construction materials, the prototype is made up of a steel structure formed by welded tubular profiles. Furthermore, for the cladding of the opaque parts, a white aluminum sandwich panel with 100 mm of XPS has been selected. Hence, for the glazed facade, a high selective double glazing has been chosen for the Reference cabin and a highly selective triple glazing with a water chamber towards the inside, for the WFG cabin. Figure 6 shows the configuration of both glazing facades. Regarding the construction materials, the prototype is made up of a steel structure formed by welded tubular profiles. Furthermore, for the cladding of the opaque parts, a white aluminum sandwich panel with 100 mm of XPS has been selected. Hence, for the glazed facade, a high selective double glazing has been chosen for the Reference cabin and a highly selective triple glazing with a water chamber towards the inside, for the WFG cabin. Figure 6 shows the configuration of both glazing facades.
The prototype has undergone an exhaustive commissioning process from the design stage to the manufacturing, factory assembly, on-site assembly, and commissioning of all the systems involved. Hence, high reflective WFG is chosen to reject sun energy and use the water chamber of the glazing facing indoors to eliminate heat loads when required. Since the energy absorption should be minimized, a high reflective coating (Xtreme 60.28) is positioned close to the outer glass pane (face 2). This coating allows for the reduction of the U-value because it can be considered a low emissivity coating at the same time (Planitherm XN). The maximum g value (g OFF ) and the minimum g value (g ON ) are almost the same and around 0.2. Table 3 shows the spectral and thermal properties of the glazing defined in Figure 6. WFG presents different values of U and g, depending on the mass flow rate. U ON and U OFF have been calculated with the equation (23)  The prototype has undergone an exhaustive commissioning process from the design stage to the manufacturing, factory assembly, on-site assembly, and commissioning of all the systems involved. Hence, high reflective WFG is chosen to reject sun energy and use the water chamber of the glazing facing indoors to eliminate heat loads when required. Since the energy absorption should be minimized, a high reflective coating (Xtreme 60.28) is positioned close to the outer glass pane (face 2). This coating allows for the reduction of the U-value because it can be considered a low emissivity coating at the same time (Planitherm XN). The maximum g value (g OFF ) and the minimum g value (g ON ) are almost the same and around 0.2. Table 3 shows the spectral and thermal properties of the glazing defined in Figure 6. WFG presents different values of U and g, depending on the mass flow rate. U ON and U OFF have been calculated with the equation (23) using ṁ = 1 L/min m 2 and ṁ = 0 L/min m 2 , respectively. The steady values of the reference glazing have been placed on the columns U OFF and g OFF .  (23) with ṁ = 1 L/min m 2 ; 2 equation (23) with ṁ = 0 L/min m 2 .
The WFG was selected to eliminate internal heat loads by circulating cool water through the water chamber facing indoors. This cool isothermal envelope allows insulating the building from external boundary conditions. Furthermore, the "Peltier" device connected to a buffer tank produces cool water for the glazing. To minimize the electrical consumption of the "Peltier" device, the temperature of the buffer tank should be close to the cool water of the glazing. When the outdoor conditions made it possible, evaporative cooling and cooling by night were considered to cool down the buffer tank. Besides, the prototype allows the understanding of the real issues of glass facades, analyzing deep concepts such as overheating and thermal mass. Figure 7 shows the prototype and the position of sensors. A short description of the monitoring devices is listed below: • Pyranometer: The Delta Ohm LP PYRA 03 Pyranometer measures the irradiance on a horizontal surface (W/m 2 ). The measured irradiance is the sum of the direct irradiance of the sun and the diffuse irradiance.  The WFG was selected to eliminate internal heat loads by circulating cool water through the water chamber facing indoors. This cool isothermal envelope allows insulating the building from external boundary conditions. Furthermore, the "Peltier" device connected to a buffer tank produces cool water for the glazing. To minimize the electrical consumption of the "Peltier" device, the temperature of the buffer tank should be close to the cool water of the glazing. When the outdoor conditions made it possible, evaporative cooling and cooling by night were considered to cool down the buffer tank. Besides, the prototype allows the understanding of the real issues of glass facades, analyzing deep concepts such as overheating and thermal mass. Figure 7 shows the prototype and the position of sensors. A short description of the monitoring devices is listed below: • Pyranometer: The Delta Ohm LP PYRA 03 Pyranometer measures the irradiance on a horizontal surface (W/m 2 ). The measured irradiance is the sum of the direct irradiance of the sun and the diffuse irradiance.

•
One wire temperature probe: The DS18B20-PAR uses Dallas' exclusive 1-Wire bus protocol that implements bus communication using one control signal. • One wire temperature probe: The DS18B20-PAR uses Dallas' exclusive 1-Wire bus protocol that implements bus communication using one control signal.

Results
The empirical tests, conducted by the experimental setup, were run over the years 2018 and 2019 to collect data in all the possible weather conditions. Some logged parameters were used as input data to the simulation tool, and other measurements compared the performance of the selected glazing and its importance on the indoor temperature. Continuous experimental data were taken to reflect the influence of variable weather conditions. The experimental data logging time step was set at 1 min.

Analysis of the Prototype. Free-Floating Temperature Regime
In 2019, the "Peltier" device was switched off, in a free-floating temperature regime. Figure 8 presents the interior temperature of both cabins and the exterior temperature over the last week of April, May, June, July, September, and October of 2019. In all the cases, the interior temperature of the WFG cabin (Tw_air) was 5 °C below the internal temperature of the Reference cabin (Tr_air) and reached similar or slightly lower values compared to the outside temperature (T_ext1). When analyzing 27/04/2019, Tw_air almost reached 20 °C, while Tr_air was 25 °C. Likewise, during the night on that day, the minimum Tw_air remained at 7 °C, while Tr_air dropped to 3 °C when the outdoor temperature reached 0 °C. Furthermore, in July, the maximum values for Tw_air were below Tr_air and T_ext1. The minimum temperature inside the WFG cabin did not drop as much as that inside the Reference cabin.

Results
The empirical tests, conducted by the experimental setup, were run over the years 2018 and 2019 to collect data in all the possible weather conditions. Some logged parameters were used as input data to the simulation tool, and other measurements compared the performance of the selected glazing and its importance on the indoor temperature. Continuous experimental data were taken to reflect the influence of variable weather conditions. The experimental data logging time step was set at 1 min.

Analysis of the Prototype. Free-Floating Temperature Regime
In 2019, the "Peltier" device was switched off, in a free-floating temperature regime. Figure 8 presents the interior temperature of both cabins and the exterior temperature over the last week of April, May, June, July, September, and October of 2019. In all the cases, the interior temperature of the WFG cabin (T w_air ) was 5 • C below the internal temperature of the Reference cabin (T r_air ) and reached similar or slightly lower values compared to the outside temperature (T _ext1 ). When analyzing 27/04/2019, T w_air almost reached 20 • C, while T r_air was 25 • C. Likewise, during the night on that day, the minimum T w_air remained at 7 • C, while T r_air dropped to 3 • C when the outdoor temperature reached 0 • C. Furthermore, in July, the maximum values for T w_air were below T r_air and T _ext1 . The minimum temperature inside the WFG cabin did not drop as much as that inside the Reference cabin. On 24/07/2019 the graph shows that the maximum T w_air was 34 • C, while the T r_air reached 40 • C, and T _ext1 remained at 38 • C. During the night, T w_air reached 25 • C, while T r_air was 21 • C, when T _ext1 dropped to 18 • C.
Besides, Figure 8 shows the curve of the indoor temperature of the WFG cabin as a damped wave shape compared to the interior temperature curve of the Reference cabin or the outside temperature. The temperature difference between day and night inside the WFG cabin did not exceed 10 • C in most cases. However, the temperature difference between day and night inside the Reference cabin was around 20 • C. Analyzing the same parameters on 24 July 2019, it can be observed that the difference between the maximum (34 • C) and the minimum (25 • C) temperature inside the WFG cabin was 9 • C, while the difference between the maximum (40 • C) and the minimum (21 • C) temperature inside the Reference cabin was 19 • C. However, the difference between the maximum (38 • C) and the minimum (18 • C) outside temperature was 20 • C. The interior temperature curve of the WFG cabin is flatter than the reference glazing one. Sustainability 2020, 12, x 12 of 24 Besides, Figure 8 shows the curve of the indoor temperature of the WFG cabin as a damped wave shape compared to the interior temperature curve of the Reference cabin or the outside temperature. The temperature difference between day and night inside the WFG cabin did not exceed 10 °C in most cases. However, the temperature difference between day and night inside the Reference cabin was around 20 °C. Analyzing the same parameters on 24 July 2019, it can be observed that the difference between the maximum (34 °C) and the minimum (25 °C) temperature inside the WFG cabin was 9 °C, while the difference between the maximum (40 °C) and the minimum (21 °C) temperature inside the Reference cabin was 19 °C. However, the difference between the maximum (38 °C) and the minimum (18 °C) outside temperature was 20 °C. The interior temperature curve of the WFG cabin is flatter than the reference glazing one. Figure 9 presents the indoor temperature of both cabins (Tr_air, Tw_air), the outdoor temperature (T_ext1), and solar radiation from 24 July 2019 to 30 July 2019. The maximum value of solar radiation Besides, Figure 8 shows the curve of the indoor temperature of the WFG cabin as a damped wave shape compared to the interior temperature curve of the Reference cabin or the outside temperature. The temperature difference between day and night inside the WFG cabin did not exceed 10 °C in most cases. However, the temperature difference between day and night inside the Reference cabin was around 20 °C. Analyzing the same parameters on 24 July 2019, it can be observed that the difference between the maximum (34 °C) and the minimum (25 °C) temperature inside the WFG cabin was 9 °C, while the difference between the maximum (40 °C) and the minimum (21 °C) temperature inside the Reference cabin was 19 °C. However, the difference between the maximum (38 °C) and the minimum (18 °C) outside temperature was 20 °C. The interior temperature curve of the WFG cabin is flatter than the reference glazing one. Figure 9 presents the indoor temperature of both cabins (Tr_air, Tw_air), the outdoor temperature (T_ext1), and solar radiation from 24 July 2019 to 30 July 2019. The maximum value of solar radiation  Figure 9 presents the indoor temperature of both cabins (T r_air , T w_air ), the outdoor temperature (T _ext1 ), and solar radiation from 24 July 2019 to 30 July 2019. The maximum value of solar radiation was above 720 W/m 2 , almost every day. Despite these high values of solar radiation, the interior temperature of the WFG cabin remained under 35 • C on 24 July 2019, while the Reference cabin temperature was above 40 • C.
Sustainability 2020, 12, x 13 of 24 was above 720 W/m 2 , almost every day. Despite these high values of solar radiation, the interior temperature of the WFG cabin remained under 35 °C on 24 July 2019, while the Reference cabin temperature was above 40 °C. The next step was to analyze a typical summer day in both cabins. Figures 10 and 11 show the glazing temperatures, measured in different points. There is a little difference between the surface temperatures of the outer glass pane in the WFG cabin (Tw_extU and Tw_extD) with the Reference cabin (Tr_extD). The effect of the water flowing through the glazing affected the indoor surface temperatures in both cabins. Figure 10 illustrates that the measured indoor glass surface temperatures in the Reference cabin (Tr_up, Tr_down) are remarkably variable during the daily hours, according to the relevant variations of the outdoor thermal conditions and the low thermal inertia of the glazing. However, there is no difference between Tr_up and Tr_down. The indoor surface temperatures in the WFG was measured at the bottom and the top of the indoor glass pane. The water heats up as it moves through the glazing, and the figure explains the influence of the water flow in the temperature distribution.  The next step was to analyze a typical summer day in both cabins. Figures 10 and 11 show the glazing temperatures, measured in different points. There is a little difference between the surface temperatures of the outer glass pane in the WFG cabin (T w_extU and T w_extD ) with the Reference cabin (T r_extD ). The effect of the water flowing through the glazing affected the indoor surface temperatures in both cabins. Figure 10 illustrates that the measured indoor glass surface temperatures in the Reference cabin (T r_up , T r_down ) are remarkably variable during the daily hours, according to the relevant variations of the outdoor thermal conditions and the low thermal inertia of the glazing. However, there is no difference between T r_up and T r_down . The indoor surface temperatures in the WFG was measured at the bottom and the top of the indoor glass pane. The water heats up as it moves through the glazing, and the figure explains the influence of the water flow in the temperature distribution.
Sustainability 2020, 12, x 13 of 24 was above 720 W/m 2 , almost every day. Despite these high values of solar radiation, the interior temperature of the WFG cabin remained under 35 °C on 24 July 2019, while the Reference cabin temperature was above 40 °C. The next step was to analyze a typical summer day in both cabins. Figures 10 and 11 show the glazing temperatures, measured in different points. There is a little difference between the surface temperatures of the outer glass pane in the WFG cabin (Tw_extU and Tw_extD) with the Reference cabin (Tr_extD). The effect of the water flowing through the glazing affected the indoor surface temperatures in both cabins. Figure 10 illustrates that the measured indoor glass surface temperatures in the Reference cabin (Tr_up, Tr_down) are remarkably variable during the daily hours, according to the relevant variations of the outdoor thermal conditions and the low thermal inertia of the glazing. However, there is no difference between Tr_up and Tr_down. The indoor surface temperatures in the WFG was measured at the bottom and the top of the indoor glass pane. The water heats up as it moves through the glazing, and the figure explains the influence of the water flow in the temperature distribution.   Figure 11 shows the surface temperatures of the WFG facade, both on the outside and the inside face of the glazing. The external pane reached temperatures above 41 • C, while the maximum surface temperature on the inner pane was 34 • C. There was no significant difference between T w_ext_U and T w_ext_D . When it comes to the interior pane, the water absorbed part of the solar radiation as it flowed. There were two sensors in the upper part (T w_up_L , T w_up_R ), and two sensors in the lower part (T w_down_L , T w_down_R ). A 2 • C difference was observed between the lower and upper probe of the inside surface of the glazing. The lower part of the inner glass reaches a maximum temperature of 32 • C, while the upper part reaches 34 • C.
Sustainability 2020, 12, x 14 of 24 Figure 11 shows the surface temperatures of the WFG facade, both on the outside and the inside face of the glazing. The external pane reached temperatures above 41 °C, while the maximum surface temperature on the inner pane was 34 °C. There was no significant difference between Tw_ext_U and Tw_ext_D. When it comes to the interior pane, the water absorbed part of the solar radiation as it flowed. There were two sensors in the upper part (Tw_up_L, Tw_up_R), and two sensors in the lower part (Tw_down_L, Tw_down_R). A 2 °C difference was observed between the lower and upper probe of the inside surface of the glazing. The lower part of the inner glass reaches a maximum temperature of 32 °C, while the upper part reaches 34 °C.  Figure 12 shows the interior temperature of both cabins and the exterior temperature from 22 July 2018 to 27 July 2018. The "Peltier" cell was in operation according to a simple control logic based on the interior temperature programmed in the monitoring control unit. It was set to operate every summer day from 12:30 to 20:00. The goal was to keep the inlet temperature between 15 and 17 °C to test the indoor air conditions inside the WFG cabin. The mean maximum temperature reached inside the WFG cabin (Tw_air) was 26.5 °C, with a mean maximum solar radiation of 720 W/m 2 . However, the mean maximum temperature of the Reference cabin (Tr_air) exceeds 37 °C, when the mean maximum outdoor temperature (T_ext1) is 34.5 °C. Therefore, there were more than 10 °C of difference inside both cabins. The Reference glazing replicated the thermal oscillations of the outside temperature, generating discomfort inside the building and contributing to overheating. The same behavior was observed when analyzing the minimum temperatures. The minimum value of T_ext1 was 10 °C, while the Reference cabin remained at 12 °C, and the WFG cabin reached 15 °C. Therefore, the thermal inertia that characterizes the WFG facade managed to dampen the oscillation of the interior temperature. The temperature oscillations of the Reference cabin were similar to the outside temperature. Hence, the WFG cabin allowed for the bringing of the maximum and minimum values close to the comfort temperature. Specifically, analyzing 23 July 2018, the mean value between the maximum Tw_air daytime temperature (26.5 °C) and the mean minimum value of night time temperature (16.5 °C) was 21.5 °C, which was very close to the comfort temperature. 3.2. Analysis of the Prototype. "Peltier" Device ON Figure 12 shows the interior temperature of both cabins and the exterior temperature from 22 July 2018 to 27 July 2018. The "Peltier" cell was in operation according to a simple control logic based on the interior temperature programmed in the monitoring control unit. It was set to operate every summer day from 12:30 to 20:00. The goal was to keep the inlet temperature between 15 and 17 • C to test the indoor air conditions inside the WFG cabin. The mean maximum temperature reached inside the WFG cabin (T w_air ) was 26.5 • C, with a mean maximum solar radiation of 720 W/m 2 . However, the mean maximum temperature of the Reference cabin (T r_air ) exceeds 37 • C, when the mean maximum outdoor temperature (T _ext1 ) is 34.5 • C. Therefore, there were more than 10 • C of difference inside both cabins. The Reference glazing replicated the thermal oscillations of the outside temperature, generating discomfort inside the building and contributing to overheating. The same behavior was observed when analyzing the minimum temperatures. The minimum value of T _ext1 was 10 • C, while the Reference cabin remained at 12 • C, and the WFG cabin reached 15 • C. Therefore, the thermal inertia that characterizes the WFG facade managed to dampen the oscillation of the interior temperature. The temperature oscillations of the Reference cabin were similar to the outside temperature. Hence, the WFG cabin allowed for the bringing of the maximum and minimum values close to the comfort temperature. Specifically, analyzing 23 July 2018, the mean value between the maximum T w_air daytime temperature (26.5 • C) and the mean minimum value of night time temperature (16.5 • C) was 21.5 • C, which was very close to the comfort temperature.

Discussion
When the glazing is part of an insulated room, the thermal problem of the glazing is coupled with the thermal problem of the room, and the indoor temperature should be determined. The indoor boundary condition disappears to be part of the solution to the thermal problem. The prototype was a rectangular room with glazing facing south. The dimensions of the room, the near and far-infrared absorption (αNIR, αFIR), and the thermal transmittance of the opaque envelope are defined in Table 4.

Discussion
When the glazing is part of an insulated room, the thermal problem of the glazing is coupled with the thermal problem of the room, and the indoor temperature should be determined. The indoor boundary condition disappears to be part of the solution to the thermal problem. The prototype was a rectangular room with glazing facing south. The dimensions of the room, the near and far-infrared absorption (α NIR , α FIR ), and the thermal transmittance of the opaque envelope are defined in Table 4. In these following test cases, outdoor temperature and solar irradiance varied during the day, and thermal performances depended on time. The indoor temperature was unknown, and it should have been obtained at the same time as the glazing temperature profile. Regarding the water flow glazing, the flow rate and the inlet temperature were constant values given by Table 5. Transient behavior occurred when the outdoor temperature and solar irradiance varied during the day. Besides, each wall had a different temperature due to the luminance of the direct beam solar radiation. Since the WFG was facing south, the north indoor wall absorbed solar radiation. The rest of this energy was diffusely reflected and created the indoor diffuse irradiance. Later, this irradiance was absorbed in each indoor surface. Hence, in this test case, the water flow glazing absorbed extra energy from the indoor irradiance. Figure 13 shows the validation of the Software Tool using real data from both cabins. The figure illustrates six days, from 22 July 2018 to 27 July 2018, in which we can see how the measured curves replicate the simulation curves in all the cases. radiation. Since the WFG was facing south, the north indoor wall absorbed solar radiation. The rest of this energy was diffusely reflected and created the indoor diffuse irradiance. Later, this irradiance was absorbed in each indoor surface. Hence, in this test case, the water flow glazing absorbed extra energy from the indoor irradiance. Figure 13 shows the validation of the Software Tool using real data from both cabins. The figure illustrates six days, from 22 July 2018 to 27 July 2018, in which we can see how the measured curves replicate the simulation curves in all the cases.

Analysis of the Reference Cabin
The reference glazing is considered a high-performance transparent envelope. It is made up of glass panes with coatings and gas cavities. Figure 14 illustrates real data and simulated indoor temperatures of the Reference cabin. In both measured and simulated curves, the interior peak temperature is slightly above 40 • C between 03:00 and 10:00 on 23 July 2018. Equation (24) represents the mean error (ME), which is the difference between the measured value and simulation results. Equation (25) represents the mean percentage error (MPE), with a total number of measurements of n = 1440.
where T Si is the simulated value, and T Ri is the measured value. By computing ME and MPE, the sample summer day, 23 July 2018, MEs and MPEs of the indoor temperature of the Reference cabin, T r_air , were lower than 0.6 • C and 3.2%, respectively. Predictions of T r_air were more accurate because the boundary conditions were more suitable to predict.
where TSi is the simulated value, and TRi is the measured value. By computing ME and MPE, the sample summer day, 23 July 2018, MEs and MPEs of the indoor temperature of the Reference cabin, Tr_air, were lower than 0.6 °C and 3.2%, respectively. Predictions of Tr_air were more accurate because the boundary conditions were more suitable to predict.

Analysis of the WFG Cabin
Two different mechanisms appeared to modify the thermal performances of the WFG. The first one was the value of the indoor temperature, which can be very high, depending on the wall insulation. The second one was the absorbed back irradiance. Depending on the insulation of the walls and the near-infrared absorptances, the water heat gain can differ when compared to the water heat gain of isolated glazing. Figure 15 shows the comparison between real data and simulation over the same day, 23 July 2018. MEs and MPEs of the indoor temperature of the WFG cabin were lower than 1.2 °C and 5.5%, respectively.

Analysis of the WFG Cabin
Two different mechanisms appeared to modify the thermal performances of the WFG. The first one was the value of the indoor temperature, which can be very high, depending on the wall insulation. The second one was the absorbed back irradiance. Depending on the insulation of the walls and the near-infrared absorptances, the water heat gain can differ when compared to the water heat gain of isolated glazing. Figure 15 shows the comparison between real data and simulation over the same day, 23 July 2018. MEs and MPEs of the indoor temperature of the WFG cabin were lower than 1.2 • C and 5.5%, respectively.  Figure 16 illustrates the time histories of indoor and outdoor air temperatures and the solar irradiance in W/m 2 on the outdoor horizontal roof surface on a sample summer day (23 July 2018). The indoor air temperature of the Reference cabin (Tr_air) varied according to the solar irradiance and the outdoor temperature (T_ext1). The measured inlet temperature (T_inletI) was input in the simulation tool as boundary conditions. By considering the activation of a solar fed "Peltier" device, the inlet  Figure 16 illustrates the time histories of indoor and outdoor air temperatures and the solar irradiance in W/m 2 on the outdoor horizontal roof surface on a sample summer day (23 July 2018). The indoor air temperature of the Reference cabin (T r_air ) varied according to the solar irradiance and the outdoor temperature (T _ext1 ). The measured inlet temperature (T _inletI ) was input in the simulation tool as boundary conditions. By considering the activation of a solar fed "Peltier" device, the inlet temperature was set between 18 and 22 • C. The operation schedule for the cooling system activation was from 12:30 to 20:00. Figure 16 reports the results of the daily analysis. The thermal effect of the "Peltier" device kept the indoor temperature below 27 • C over the hottest hours of the day.  Figure 16 illustrates the time histories of indoor and outdoor air temperatures and the solar irradiance in W/m 2 on the outdoor horizontal roof surface on a sample summer day (23 July 2018). The indoor air temperature of the Reference cabin (Tr_air) varied according to the solar irradiance and the outdoor temperature (T_ext1). The measured inlet temperature (T_inletI) was input in the simulation tool as boundary conditions. By considering the activation of a solar fed "Peltier" device, the inlet temperature was set between 18 and 22 °C. The operation schedule for the cooling system activation was from 12:30 to 20:00. Figure 16 reports the results of the daily analysis. The thermal effect of the "Peltier" device kept the indoor temperature below 27 °C over the hottest hours of the day. The daily cooling demands are shown in Figure 17. The daily energy that the flow of water can absorb was calculated using Equation (15). WFG prevented this energy from entering the cabin. The "Peltier" device provided the cooling power that explains the difference in temperatures between the The daily cooling demands are shown in Figure 17. The daily energy that the flow of water can absorb was calculated using Equation (15). WFG prevented this energy from entering the cabin. The "Peltier" device provided the cooling power that explains the difference in temperatures between the WFG cabin and the Reference cabin. The energy absorbed by the WFG was 117.5 Wh in 0.5 m 2 , so the ratio of energy per area was 0.235 kWh/m 2 .

Final Energy and Cost Considerations
This section has considered best available technologies (BAT), which are innovative and economically viable [48]. Hydronic technologies are compatible with radiant systems, such as WFG and radiant floors and walls. Air-to-air heat pumps are used to compare the final energy consumption

Final Energy and Cost Considerations
This section has considered best available technologies (BAT), which are innovative and economically viable [48]. Hydronic technologies are compatible with radiant systems, such as WFG and radiant floors and walls. Air-to-air heat pumps are used to compare the final energy consumption and cost of different cooling systems. WFG can be a part of hydronic HVAC systems, and it is compatible with water-to-water heat pumps. Tables 6 and 7 illustrate the final energy consumption with different energy generators. It takes into account the effect of the operative temperature of the system [63]. The performance of water-to-water heat pumps (WWHP) depends on the inlet temperature of the WFG (θ IN ), and the source inlet temperature (T s,i ) in the heat pump. A typical value of T s,i ranges from 20 • C in ground source heat pumps (GSHP) to 35 • C in other WWHPs. The source and load sides are relevant when it comes to calculating the performance of the cooling device. Air-to-air heat pumps (AAHP) for heat recovery on ventilation are also analyzed. The parameters that influence air-to-air heat pumps' performance are the dry bulb exterior air temperature (T e_db ) and the dry bulb interior return air temperature (T ri_db ).  PEF stands for the primary energy factor from final energy (FE) to non-renewable final energy (NRFE). The Spanish code for thermal systems in buildings (RITE) recommends a value of 1.954. The factor of CO 2 emissions for electricity was 0.331 [64]. The RITE aims to establish primary energy factors and CO 2 emission factors, for each final energy consumed by buildings in Spain and for each geographic area with a different electricity generation source. Table 6 shows that the best performance belongs to WWHP when the inlet temperature (θ IN ) is close to 18 • C, and the source inlet temperature (T s,i ) in the heat pump is 20 • C. Figure 16 shows that it is possible to keep the cabin's indoor temperature within the comfort range by setting θ IN between 18 • C and 22 • C. Other systems, such as fan-coils, need lower operating temperatures. The lower the difference between T s,i and θ IN the higher the EER coefficient. The NRFE consumption of GSHP is 0.065 kWh per m 2 of WFG, and the average NRFE use of air-to-air heat pumps is 0.155 kWh per m 2 of Sustainability 2020, 12, 5734 20 of 23 WFG. When it comes to CO 2 emissions, air-to-air heat pumps generate an average of 0.027 KgCO 2 , twice as much as the CO 2 emitted by GSHP.

Conclusions
This paper has studied a model to assess innovative building envelopes' energy performance, such as water flow glazing (WFG), the system equations for load calculation, and the relationships with steady and transient parameters. Some of the magnitudes can be measured accurately in the prototype. The presented tool has been developed and tested by the authors. Details of the prototype and the on-site measures have been used to validate the tool. The analysis included a free-floating temperature regime and a cooling system with simple logic to keep the prototype within a comfort range. Results included simulated indoor air and glazing temperatures along with the potential final energy savings.
1. When the WFG cabin's interior temperature was below the exterior temperature, the WFG facade cooled down the room. As the "Peltier" device was not in operation over 2019, it can be concluded that WFG working on a free-floating temperature regime, without auxiliary energy systems, results in smaller indoor temperature fluctuations.
2. Circulating water increased the temperature gap between external and internal glass panes. The outer pane reached temperatures above 41 • C, while the maximum surface temperature on the inner pane was 34 • C. The reference glazing showed a smaller gap between outer pane (41 • C) and inner pane (38 • C). The reduction of radiant temperature of indoor envelopes can improve the occupant's comfort.
3. The damping effect on the WFG cabin's temperature is shown in Figure 12. The WFG system provided the facade with thermal inertia, and the cabin did not suffer large thermal oscillations between day and night. This effect can increase the thermal comfort inside the building and reduce energy consumption.
4. The WFG increased the thermal inertia of the facade. Once the maximum temperature was reached, the interior of the WFG cabin cooled down more slowly than the Reference cabin did.
There was a good agreement between the simulation and the real data from the prototype. MEs and MPEs of the indoor temperature in the WFG cabin, were lower than 1.2 • C and 5.5%, respectively. The simulation results of the Reference cabin were more accurate because the boundary conditions were more suitable to predict.
The weather and indoor conditions impact the efficiency coefficients of heat pumps. The EER values of the cooling systems were evaluated for different combinations of indoor and outdoor conditions. Ground source heat pumps (GCHP) coupled with borehole heat exchangers in a closed loop are very effective and make the most of the near-constant ground temperature over the year.
Finally, if the electricity is supplied from solar cells, it is possible to use a renewable and CO 2 free energy source to provide thermal comfort.