Next Article in Journal
Experimental Analysis and Optimization of an R744 Transcritical Cycle Working with a Mechanical Subcooling System
Next Article in Special Issue
On the Performance of Night Ventilation in a Historic Office Building in Nordic Climate
Previous Article in Journal
Optimal Coordination Strategies for Load Service Entity and Community Energy Systems Based on Centralized and Decentralized Approaches
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application and Validation of a Dynamic Energy Simulation Tool: A Case Study with Water Flow Glazing Envelope

by
Belen Moreno Santamaria
1,
Fernando del Ama Gonzalo
2,*,
Danielle Pinette
2,
Roberto-Alonso Gonzalez-Lezcano
3,
Benito Lauret Aguirregabiria
1 and
Juan A. Hernandez Ramos
4
1
Department of Construction and Architectural Technology, Technical School of Architecture of Madrid, Technical University of Madrid (UPM), Av. Juan de Herrera, 4, 28040 Madrid, Spain
2
Keene State College, 229 Main St, Keene, NH 03435, USA
3
Escuela Politécnica Superior, Universidad San Pablo-CEU, Montepríncipe Campus, Boadilla del Monte, 28668 Madrid, Spain
4
Department of Applied Mathematics, School of Aeronautical and Space Engineering, Technical University of Madrid (UPM), Plaza Cardenal Cisneros 3, 28040 Madrid, Spain
*
Author to whom correspondence should be addressed.
Energies 2020, 13(12), 3203; https://doi.org/10.3390/en13123203
Submission received: 30 April 2020 / Revised: 11 June 2020 / Accepted: 13 June 2020 / Published: 19 June 2020
(This article belongs to the Special Issue Thermal Energy Storage in Building Integrated Thermal Systems)

Abstract

:
The transparent materials used in building envelopes significantly contribute to heating and cooling loads of a building. The use of transparent materials requires to solve issues regarding heat gain, heat loss, and daylight. Water flow glazing (WFG), a disruptive technology, includes glazing as part of the Heating, Ventilation and Air Conditioning (HVAC) system. Water is transparent to visible wavelengths, but it captures most of the infrared solar radiation. As an alternative to fossil fuel-based HVAC systems, the absorbed energy can be transferred to the ground through borehole heat exchangers and dissipated as a means of free-cooling. Researchers of the Polytechnic University of Madrid have developed a software tool to calculate the energy balance while incorporating the dynamic properties of WFG. This article has studied the mathematical model of that tool and validated its ability to predict energy savings in buildings, taking spectral and thermal parameters of glazing catalogs, commercial software, and inputs from the measurements of the prototypes. The results found in this article showed that it is possible to predict the thermal behavior of WFG and the energy savings by comparing the thermal parameters of two prototypes. The energy absorbed by the water depends on the mass flow rate and the inlet and outlet temperatures.

Graphical Abstract

1. Introduction

The percentage of glazing in buildings affects significantly heat gains and losses [1,2]. The visual and thermal properties of glass, the fluid used to fill the space between glass panes, and the use of coatings alter the energy absorption in each layer of the glazing [3]. Article 6 of the Energy Performance of Buildings Directive (EPDB 2010) states that building designers have to consider renewable sources of energy for Heating, Ventilation and Air Conditioning (HVAC) and high-performance envelopes [4]. The first step to accomplish the goal of zero energy buildings is to minimize the heating loads. The second step is to select and size innovative and more efficient devices for heating and cooling, and the third step includes the use of renewable energies [5].

1.1. Closed-Loop Ground Source Systems

Borehole heat exchangers have been introduced as an alternative to fossil fuel-based HVAC systems in Europe and North America over the last decades [6]. Open-loop systems use groundwater as a geothermal heat source and may affect local temperature evolution [7]. Closed-loop systems use a fluid that circulates through borehole heat exchangers to transfer the thermal energy from the subsurface to a hydronic HVAC system. Some articles report tests on closed-loop borehole heat exchangers and their capability to extract energy from the ground in winter. A refrigerant fluid transports the heat. This solution is optimum in cold countries when the systems allow heat injection into the ground in summer [8,9]. At a certain depth, the temperature of the ground remains steady over the year, so borehole heat exchangers are considered an optimum solution when the available area is limited [10]. Plastic pipes made of polyethylene or polypropylene are deployed in the boreholes and the irregular space that remains between the pipe and the surrounding soil has to be filled with bentonite, a fluid that absorbs water and expands. This property enhances heat exchange and performance as a thermal conductor [11,12]. Ground Source Heat Pumps (GSHP) exploit the nearly constant ground temperature over the year to increase their performance coefficients. GSHP fit well in the low-CO2 energy supply system [13]. The soil mass is considered a great means of thermal energy storage, and the temperature oscillations are reduced at 10-20 m depth over the seasons. [14].

1.2. The Water Flow Glazing (WFG)

The Architectural Engineering and Construction (AEC) industry has focused on the building envelope performance to reduce the heating and cooling loads. The elements used in the building envelope contribute decisively to the energy performance of buildings. Using transparent materials requires first, to understand their spectral properties and second, to develop systems that can solve some of the issues regarding heat gain, heat loss, and daylight [15]. Solar heat gain from infrared radiation is the factor that contributes the most to the cooling loads of a building [16]. Water flow glazing (WFG) is a disruptive technology that includes glazing as part of the HVAC system [17]. WFG includes a circuit that allows fluid to flow through space between two glass panes. In a double glazing, the transmissivity coefficient depends on the radiation’s wavelength. Clear glass is transparent to visible and near-infrared (NIR) wavelengths. However, water is opaque to NIR wavelengths, while its visible transmittance is very high [18]. Therefore, water captures solar NIR and increases its temperature through the window. WFG reduces solar gains through transparent envelopes without hindering light transmission. The flow of the water enables the building envelope to apply energy-saving measures, such as energy storage in buffer tanks, and facade homogenization [19].
The lack of thermal inertia affects the thermal behavior of glazing [20] because the daily oscillations in the outdoor temperature are replicated indoors. The thermal inertia buffers the daily fluctuations in the outside temperature. By decreasing these oscillations, the HVAC systems need less power. In warm dominant climates, variable thermal inertia enables better management of the energy flow to and from the building [21,22]. Using WFG, the strong sunlight of the summer is kept out of the building. Releasing the excess of heat in the ground during the summer allows to increase the thermal inertia of a light and transparent envelope. Therefore, the regulation of the glazing’s thermal inertia can be achieved by changing the flow rate that circulates through the panels. The lower the flow rate, the lesser the damping effect.
Active water flow glazing raises its thermal inertia compared to conventional double or triple glazing. The water inside is in continuous movement as it manages to release the absorbed heat. Therefore, it reaches to thermal inertia comparable to other traditional construction elements (brick and concrete walls) by simply connecting the system to energy sinks, such as borehole heat exchangers or buffer tanks.

1.3. Dynamic Building Energy Modeling

Accurate prediction of the behavior of glazing facades has to include the quality of construction and permeability to the passage of air and vapor. Those aspects depend on the design skills, as much as on the quality of the construction. The dynamic thermal performance of an envelope produces variations in thermal properties over the day, such as variable thermal inertia, which the Fourier model does not take into account [23]. Thus, developing and validating new mathematical models that can simulate dynamic properties is an important goal. Energy Plus is the most popular engine for calculating energy load and energy system performance of buildings [24,25], but it does not include a module to simulate WFG.

1.4. Innovation and Objectives

Engineers and researchers of the Polytechnic University of Madrid have developed a tool to calculate preliminary energy balances of a building. Under the two main tabs, thermal behavior and sun energy, the interface allows the introduction of parameters, such as glass types, coatings, fluids filling the space between glass panes, thermal mass of the building, and internal heat gains. The software library is an open-source project and can be used by developers to validate results. According to some studies, architects and building designers are prone to use simple energy simulation tools than others that might include more options but are more complicated [26,27].
The transient state of the building envelopes is affected by changes in temperature and solar radiation. The steady-state model is not accurate when it comes to dynamic forms of heat transfer, additionally affected by indoor HVAC systems. Remote sensing systems can be used in future studies to compare the actual energy performance with the results of simulation models [28,29]. Other technologies, such as switchable glazing, can adapt to different weather conditions. Polymer dispersed liquid crystal (PDLC) devices can adapt to different weather conditions, but their expected range of U value and g-factor is not as broad as expected for WFG [30]. Suspended Particle Devices (SPD) switches from transparent to colored due to the suspension of freely arranged particles between two glass panes. The electrochromic (EC) glass varies transmission and reflection parameters by electrical stimuli to mobile ions in the electrochromic layer. The solar heat gain coefficient of SPD and EC ranges from 0.5 in the transparent state to 0.06 in the non-transparent state [31]. The addition of double glazing to an SPD can vary the U value from 2.98 W/m2 K to 1.99 W/m2 K [32]. The heat transfer coefficient of SPD ranges between 5.02 W/m2 K and 5.2 W/m2 K for opaque and transparent states [33]. The critical aspects remain in the high cost of the products and a lack of standardization in the manufacturing process [34].
This paper aims to study the dynamic thermal parameters of WFG. The prototype with WFG in facades and the roof, connected to borehole heat exchangers, was built. The only additional energy source was the electricity needed to run the circulation water pumps. The flow of water requires a circulating water pump that consumes electric energy, though the consumption of that pump is low compared with the heating and cooling savings over the year. Some articles show that photovoltaic panels integrated into the WFG facade can provide the water pump with energy [35]. A second prototype was made with the same dimensions and double glazing with an air cavity.
There were three main objectives in the analysis of the WFG prototype. First, it allowed comparing the indoor temperatures of the WFG cabin and the Reference cabin. Second, the simulation tool based on the mathematical model to predict the performance of WFG was validated using real data. Finally, the effect of the borehole heat exchangers connected to the WFG cabin and its influence on the heat gains in the water chamber were taken into account. By comparing the indoor temperatures of two prototypes, the influence of WFG as a means of energy manager was tested. Not only has the indoor temperature in both prototypes been analyzed, but also the capability of WFG to absorb energy with a specific mass flow rate. The results of this study show that the solution of double glazing with a circulating water chamber is a less polluting and more efficient option than the systems currently used, which are, most of them, fed by fossil fuels. This study aims to validate a new building energy simulation tool by comparing numerical results and real data obtained in test rooms. The empirical tests were carried out over a year in a free-floating indoor temperature regime with variable outdoor conditions.
The challenges of WFG are related to the ability to produce large panels and their behavior in the long term. Although considerable efforts have been made in the industrialization process, there is still low market penetration because of the limited information between professionals and consumers. It is not the aim of this article to study the efficiency of closed-loop shallow geothermal systems. Data from the borehole heat exchangers have been used as input to validate the results of the mathematical model.

2. Materials and Methods

Since commercial building energy simulation tools do not include WFG as an option for energy management, it is necessary to validate the outputs of simulation with data from real prototypes. The first subsection described the geometry, the energy management, and the materials used in two cabins. The second one showed the equations and the mathematical models of the simulation tool. Finally, the last subsection set the criteria to select the spectral and thermal parameters of the WFG.

2.1. Description of the Prototype

For this prototype, active water flow glazing is combined with borehole heat exchangers. On one hand, the active envelope circulates the water at a temperature close to comfort conditions with little energy cost, making the most of the ground’s thermal inertia. Another beneficial effect is that the ground temperature increases during the summer as heat absorbed by the water flowing, through the facades and roof, is transferred to the ground by the borehole heat exchangers.
The prototype is located near Madrid, Spain (latitude 40°36′42” N, longitude 2°26′57” O, altitude 1111 MAMSL). The dimensions of the prototypes are 2 m × 2 m × 2 m, with a total volume and floor area of 8 m3 and 2 m2, respectively. The prototypes are made of aluminum frames without thermal breakage. The North facades consist of sandwich panels with a thermal transmittance of U = 0.6 W/m2 C. The makeup of the panels is: 1 mm of aluminum, 5 cm extruded polystyrene, 1 mm aluminum. The floor of the prototype is constructed with the same panels as the North facade. The first prototype cabin incorporates WFG. This prototype will be referred to as WFG cabin. While the second prototype cabin, of the same dimensions, orientation, and glazing, has an air chamber in lieu of WFG. This prototype will be referred to as Reference cabin (RC). For the WFG prototype, four borehole heat exchangers are buried 50 m below the surface. There are two different closed-circuit circulating systems. For the primary circulating system, pipes with refrigerant fluid connect from the borehole heat exchangers to the circulating pump. There are two secondary systems, one for the WFG on the vertical facade and another for the WFG on the horizontal roof. For this circulating system, water connects the circulating pump to the WFG. Figure 1 represents the schematics of the energy management system of the WFG cabin. The primary circulating system goes from the plate heat exchangers (1) to the borehole heat exchangers (4). The two secondary circulating systems go from the circulating pump (2) to the WFG facade (3). Both secondary circuits contain a precision flow meter and a meter for monitoring the inlet and outlet temperature.
Calculating the required length of the borehole heat exchangers depends upon the geology and the thermal conductivity of the soil [16]. Typical values for specific heat extraction range between 40 and 70 watts per meter of borehole length.
Weather data over the year were recorded and compared with the building energy simulation (BES) model. The indoor temperature sensors were placed below the roof frame in both prototypes. The outdoor and indoor air temperatures were monitored by several thermocouples (error ±0.5 °C). A pyranometer with a spectral response from 280 to 2800 nm, a linear response up to 3000 W/m2, and an absolute error of 3% measured the solar radiation on the horizontal roof.
Figure 2 shows the real prototypes on location. The climate of the location of the prototypes experiences severe seasonal changes in temperature. The cabin on the left is the WFG cabin and the right is the Reference cabin (RC).

2.2. Mathematical Model

Previous studies have shown a mathematical model that explains the behavior of WFG [36,37,38]. This study used the equations and the mathematical models from those articles, along with data from commercial software, to assess the performance of the prototypes described in Section 2.1. The internal temperature of the prototype obtained experimentally was compared with the results from the mathematical model. Commercial catalogs provided thermal and spectral properties of glazing.

2.2.1. Double Glazing with a Gas Cavity

Heat flux inside double glazing with a gas cavity is caused by the direct and diffuse solar radiation, i0, and the difference between the indoor and outdoor temperatures, e − θi). The thermal transmittance, U, is the parameter that explains the heat transfer. The units of measurement are W/m2 K, and it is determined by the European Standard [39,40]. Two components explain the value of the g factor. The first one is the direct solar transmittance, T. The second one is the fraction of the solar radiation absorbed by the glazing, which is later transferred inside by natural convection. Equation (1) represents the heat flux that goes through the glazing. The heat flux q depends on U and g [41].
q = U ( θ e θ i ) + g i 0 ,
the thermal transmittance, U, and the g factor are defined by:
1 U = 1 h e + 1 h i + 1 h g ,
g = T + A I ,
A I = U [ A 1 ( 1 h e ) + A 2 ( 1 h e + 1 h g ) ] ,
where hi is the interior heat transfer coefficient, he is the exterior heat transfer coefficient, hg is the heat transfer coefficient of the air chamber, and A1 and A2 are the thermal absorptances of the exterior and the interior glass panes, respectively.

2.2.2. The Influence of Water Flow Glazing

With a fluid circulating between two glass panes, there are more parameters to take into account to evaluate the thermal performance. The temperature and the mass flow rate of the fluid are variable and this affects the heat flow through the glazing. The result of adding the transmitted solar energy radiation (T i0) and the secondary internal heat flux qI is the total indoor heat flux, q
q = T   i 0 + q I .
A general solution of heat flux for WFG is given in Reference [36], adequate for all glazing configuration, is presented below. The algebraic equation is valid, assuming that the mass flow rate is uniform, the values for convective and radiative coefficients do not change with time, and the thermal resistance and thermal mass of the water and glass panes are negligible.
q I = U ( θ e θ i ) + U w ( θ I N θ i ) + A I i 0 ,
where Uw is the thermal transmittance between the interior and the water chamber, U is the thermal transmittance of the glazing, and θIN is the inlet temperature of the fluid that enters the WFG. Taking into account Equations (3), (5), and (6), the expression for the total indoor heat flux, q, is:
q = U ( θ e θ i ) + U w ( θ I N θ i ) + g i 0 ,
where:
U = U i   U e m ˙ c + U e + U i ,
U w = U i   m ˙ c m ˙ c + U e + U i .
Water flow glazing has the potential to modify its thermal properties by changing the mass flow rate per unit of surface, . U and Uw depend on the mass flow rate , which is considered uniform within the glass pane. The product ṁc is the capacity of the water flow to absorb energy. Ui and Ue are thermal transmittances that can be obtained using the convective heat coefficients, he, hi, hg, and hw. Ue is the thermal transmittance from the water chamber to the outdoors. Ui is the thermal transmittance from the water chamber to indoors. The definition of the g factor for WFG is the same as stated in Equation (3) for double glazing. However, the expression for determining AI depends on the flow rate.
A I = ( U i   m ˙ c + U e + U i ) A v + A i .
In a WFG with 2 glass panes and a water chamber:
1 U e = 1 h e + 1 h w ,
1 U i = 1 h i + 1 h w .
The absorptance, Av, depends on the energy absorbed by the glass panes and by the water:
A v = A 1 ( U e h e ) + A 2 ( U i h i ) + A w ,
where A1 and A2 are the absorptances of the outdoor and the indoor glass pane, respectively, and Aw is the absorptance of the water chamber. The Equations (14) and (15) determine absorptance Ai and Ae.
A i = A 2 ( 1 U i h i ) ,
A e = A 1 ( 1 U e h e ) .
Figure 3 illustrates two panes of glass that make up a water chamber with an ascending flow of water, where θIN is the inlet temperature of the water chamber and θOUT is the temperature at the outlet of the water chamber. Out of the incoming solar radiation (i0), the g factor represents the percentage of the heat that ends up going through the glazing. T is the transmittance of the glass. The absorbed energy is divided into three fractions. A constant part of the secondary internal heat transfer factor, Ai, is transferred indoors. Another constant part of the secondary external heat transfer factor, Ae, is transferred outdoors. Both, Ai and Ai, do not depend on the flow rate. The remaining net absorptance of the water chamber, Av, is transferred indoors, (Av int), outdoors, (Av ext), and transported by the water chamber (Av water). The three factors depend on the flow rate. Using the hypothesis stated at the beginning of this section, the thermal profile inside each glass pane does not change with time. T denotes the transmittance of the glazing. The spectral properties depend on the incoming solar radiation, i0.
The g factor and transmittance of double glazing are constants. When it comes to WFG, these values can be actively controlled by modifying the flow rate, .
g = T + A I .
Introducing the value of AI from Equation (10):
g = T + A i + A v   i n t .
Av int impacts the g factor. Equation (8) shows that Av int is affected by the flow rate.
A v   i n t = ( U i   m ˙ c + U e + U i ) A v .
The g factor for WFG is:
g = ( U i   m ˙ c + U e + U i ) A v + A i + T .
If the flow rate is much higher than the sum of indoor and outdoor thermal transmittances (ṁc >> Ue+ Ui), the values of g and U are denoted with the superscript ON. The g factor of Equation (20) becomes:
g O N = A i + T .
The condition ṁc >> Ue+ Ui establishes the range of validity for the ON values. By using the thermal parameters detailed above, this condition becomes > 0.0073 kg/m2s ( > 0.4 l/m2min) for double glazing with a water chamber. Figure 4 shows the U value from Equation (8) and the g factor from Equation (19) for different spectral parameters. U ranges from 4.8 W/m2 K when there is no flow (UOFF) to 0.66 W/m2 K when = 0.039 Kg/m2s. The g-factor depends on , but also on the spectral properties of the glass.
The superscripts OFF and ON are used when the flow rate is stopped and when it is very high, respectively. When the flow rate is high (ṁc >> Ue + Ui), the Equation (7) becomes:
q = U w O N ( θ I N θ i ) + g O N i 0 ,
and Equation (9) becomes:
U w O N = U i .
Considering the water chamber as an open system, the power per unit of surface that a flow of water can transport is:
P = m ˙ c ( θ O U T θ I N ) .
The heat flux that goes through the two glass panes and the water is shown in Figure 3 and can be split into q1, q2, q3, q4. The unknown θOUT is obtained by solving the Equations (25)–(31).
q 2   = q 1   + A 1   i 0   ,
q 3   = q 2   + A w   i 0   + m c ˙ ( θ I N θ O U T ) ,
q 4   = q 3   + A 2   i 0   ,
q 1   = h e   ( θ e θ 1 ) ,
q 2   = h w   ( θ 1 θ w ) ,
q 3   = h w   ( θ w θ 2 ) ,
q 4   = h i   ( θ 2 θ i ) .
Thus, the outlet temperature of the water chamber is:
θ O U T = ( A v   i 0   + U i   θ i + U e   θ e + m ˙ c θ I N m ˙ c + U e + U i )
where is the mass flow rate per unit of area, c is the specific heat of water, Av is the net absorptance of the water chamber, i0 is the direct sun radiation, θIN is the temperature of the inlet water, Ui is the interior thermal transmittance, and Ue is the exterior thermal transmittance. These last two thermal transmittances Ui and Ue measure the heat transfer between the water chamber and indoors and outdoors, respectively. (Avi0) is the sum of the energy absorbed by the water chamber and the energy transferred by convection because of the glass panes absorption. Av comes from Equation (13). It measures the fraction of energy that is absorbed by the water. Regardless of the flow rate, the maximum amount of energy that the water chamber can absorb is Avi0. When reaches a specific value, gw becomes negligible and the absorbed energy is not transferred indoors. Therefore, the absorbed energy, Avi0, is transported by the water chamber.

2.3. Glass Selection

The simulation tool, used to select the thermal and spectral properties, is based on previous papers that have shown a mathematical model to predict the performance of WFG, along with other radiant systems [17,42]. The selection of the glazing impacts the design of the HVAC system of buildings. This subsection aims to select the WFG parameters by using the simulation tool. The four tabs include:
1.
Energy balance considerations based on the location, including potential sun energy.
2.
Spectral properties of glass panes and coatings.
3.
Thermal simulator of the WFG modules.
4.
Thermal simulator of simplified rooms combining WFG and non-transparent walls, roofs, and floors.
Designing a WFG facade starts with the analysis of the potential solar energy. A few factors include location, usage, and the orientation of different facades. This information is shown in graphs that represent the energy balance considerations for a particular day, week, month, or year. To accomplish the energy balance between the demand and the energy absorption, the dynamic g factor and thermal transmittance are considered. The wavelength of solar radiation, the orientation, and the angle are the factors that affect spectral properties of glass, such as transmittance, reflectance, and absorptance. The second phase is to define the simulation inputs, such as exterior temperature, direct and diffuse solar radiation, the selected flow rate, and inlet temperature in the WFG. The outputs of the simulation tool are the indoor temperature, the outlet temperature in the WFG, and the water heat gains calculated over some time. Finally, the simulation tool allows the user to introduce the dimensions of rooms, the thermal properties of opaque facades, and the interior absorption properties of walls, roofs, and floors. The spectral properties of the glass have been taken from the software Optics and Window [43], which are based on the International Glazing Database (IGDB). Manufacturers provide optical and spectral data of their products with standards of quality. The composition of the WFG comprises a laminated glass pane of 6 + 6 mm, a water chamber of 16 mm, and another laminated glass pane of 8 + 8 mm. The glass panes selected were Planiclear 6 and 8 mm thick (manufactured by Saint Gobain), and the Poly-Vinyl Butyral (PVB) layers were standard (1 × 0.38 mm). Table 1 presents the thermal and spectral parameters of the WFG, along with the Reference glazing.

3. Results

This section aims to validate the model presented in Section 2.2 with real data from the prototype described in 2.1. Regarding the thermal parameters, the coefficients hi and he are taken from the European Standard [39,40], being hi = 10 W/m2 K and he = 25 W/m2 K for horizontal glazing with heat flux going upwards. The value of hw is the result of using Equation (11) of Reference [17], taking into account the Nusselt number of 7.541. The thermal resistance of a 16 mm water layer is 0.026 m2K/W. Finally, hw = 452 W/m2 K.
The values of Ue and Ui are obtained with Equations (11) and (12).
(1/Ue) = (1/he) + (1/hw) = 1/10 + 1/452 = 0.102; Ue = 9.80.
(1/Ui) = (1/hi) + (1/hw) = 1/25 + 1/452 = 0.042; Ui = 23.69.
Typical values for Av are given in previous articles [36] depending on the sum of the absorptances of each layer, which makes up the glass pane. A range of values from 0.331 to 0.429 has been considered in this case study. The flow rate should be high enough to transport the absorbed energy. Then, if ṁc is much higher than Ui + Ue, the proportion of the absorbed energy that is transferred indoors, AI, and outdoors, AE, is minimized.
In this section, real data of monitoring temperatures and power efficiency of the model described in Section 2.2 are presented. By analyzing the performance of the WFG cabin, the energy strategy may be refined or improved, achieving energy savings. Hence, the power performance of the WFG envelope is obtained by measuring the inlet and outlet temperature and the flow rate of each facade.

3.1. Data from the Prototype

Real data, obtained from the prototype in three consecutive days of Summer, allows the comparison between the WFG cabin and the Reference cabin. It is essential to highlight that the prototype is not designed to assure comfort temperature inside but to evaluate the potential solar energy that WFG envelopes can absorb for Energy Management strategies. When the sun radiation impinges in the glazing units, an amount of energy is absorbed in the water chamber. This harvested energy can be distributed through the rest of the glazing facades of the cabin by homogenizing temperature gradients. Besides, energy surplus can be stored in the geothermal boreholes heat exchangers.
Figure 5 shows the curves of the outdoor temperature (printed in grey), the indoor temperature of the WFG cabin (printed in cyan), and the indoor temperature of the Reference cabin (printed in red). As can be seen, the peak outdoor temperature reaches 31.5 °C, and the peak indoor temperature in the WFG cabin rises to 38.5 °C. In contrast, the peak indoor temperature of the Reference cabin amounts to 53.5 °C. Hence, the temperature inside the WFG cabin is reduced in 15 °C in relation to the Reference cabin as a direct consequence of the fact that the solar energy absorbed by the water chamber is transferred to the ground before being transmitted indoors. Therefore, although the composition of the glass facade is the same for both prototypes, the solar factor of the WFG cabin when the system is operating at the design flow rate ṁ = 0.9 l/min m2, is reduced to 0.54, compared to the 0.67 provided by the manufacturer, as it is shown in Table 1.
Furthermore, the effect of thermal inertia derived from geothermal boreholes contributes to cushioning the temperature curve of the WFG prototype compared with the Reference one, as it is shown in the graphs of Figure 4. Even, the benefits of the effect of thermal inertia are also reflected when the peak temperature of the WFG prototype is one hour behind from the Reference cabin.
Figure 6 shows the inlet and outlet temperatures of the fluid taken from the horizontal glazing. As can be seen in the graph, when the inlet temperate is lower than the outlet temperature, the solar energy is absorbed in the water chamber over the daylight hours and transferred to the ground through the borehole heat exchangers. Likewise, over the night, the temperature of the glazing is cooled down, and the inlet temperature is higher than the outlet temperature.

3.2. Validation of the Mathematical Model

The design mass flow rate of the water flow glazing was taken from real data from the prototype. The design mass flow rate was 3.2 L/min through 3.5 m2 of horizontal WFG, so the steady value = 0.9 L/min m2 (0.015 kg/s m2). This value has been introduced in the tool, and it is high enough to validate the hypothesis in Equation (21). The heat capacity of the fluid, a mixture of water and glycol, is c = 3600 J/kg K. The boundary conditions were the same in both the mathematical model and the real prototype. The inlet temperature of the fluid was taken from real data, and the mass flow rate per unit of surface used in the simulation was = 0.015 kg/m2s.
A pyranometer installed on the roof of the prototype measures global solar irradiance. Data acquisition is made through the monitoring system. Figure 7 shows the irradiance curve corresponding to the selected days (23/08/2010, 24/08/2010, and 25/08/2010). As it is shown, the selected days were clear. When the day is clear, direct beam radiation predominates over diffuse radiation. However, when the day is cloudy, almost all the irradiance is diffuse. Likewise, the simulation curve (printed in cyan) shows the typical irradiance curve reaching maximum levels of around 800 W/m2. Over the first hours in the morning (from 09:00 to 13:00 h), and in the late afternoon (after 18:00 h), the roof is shaded and the irradiance is mainly diffuse reaching values around 200 W/m2. However, throughout the central hours of the day, the roof is exposed to direct solar radiation.
Figure 8, Figure 9 and Figure 10 show the validation of the prototype and the Software Tool using real data. These figures show the main results of three days in August (23/08/2010, 24/08/2010, and 25/08/2010) where the measured curves (printed in a continuous line) replicate the simulation curves (printed in a continuous line with marked dots) in all the cases.

4. Discussion

Double glazing has constant g and U values. In WFG, these values can be actively controlled or modified by the flow rate . A WFG allows controlling the incoming solar energy into a building by changing its flow rate. A primary system should be designed for the building to manage this absorbed energy. The studied dynamic simulation tool assessed the energy performance of new building technologies that are not implemented yet in commercial building energy simulation software. Real data obtained in the summer period over the selected days (23/08/2010, 24/08/2010, and 25/08/2010). These days allow validating the simulation tool, considering the highest potential of the WFG over the summer, when the sun irradiance reaches maximum levels. The outdoor temperature and solar irradiance vary during the day and thermal performances are functions of time. Hence, transient behavior is expected.
Figure 7 shows how the measured curves of the solar irradiance (printed in a continuous line) replicate the simulation curve (printed in a continuous line with marked dots), following the same shape and remarkably similar values slightly over 800 W/m2, which would be expected for the summer period in this location. This prototype is simulated in the actual location and the weather file is the standard EPW file (EnergyPlus Weather).

4.1. Comparison between the WFG Cabin Data and the Simulation

Concerning the temperatures of the WFG cabin, Figure 9 shows that real data curves do not precisely replicate the behavior according to the simulation both in terms of the interior and exterior temperatures. The interior temperature of the WFG cabin is undoubtedly the most critical variable to validate with the built prototype since it can be affected by infiltrations or thermal bridges coming from the aluminum profiles or the connection between the slab and the glazed walls. These would explain the temperature oscillations in the curves, especially on 23 August. However, in both measured and simulated curves, the interior peak temperature remains below 40 °C. In contrast, at minimum temperatures, it is the actual data curve that registers a peak value of less than 14 °C, compared to the 15.5 °C provided by the simulation tool. Likewise, the exterior temperature presents very slight differences between the real data curve and the simulated curve. Some oscillations in the former curve between 09:00 and 13:00 on 23 August and between 17:00 and 19:00 on 25 August would indicate that although the days were bright and sunny, some shadow coming from the prototype itself could be slightly influencing the exterior temperature probe. The difference between the measured value and simulation results is the Mean Error (ME). The Mean Percentage Error (MPE) is calculated to validate the simulation results with the real data. The total number of measurements was n = 10,167. Both values are shown in Equations (32) and (33).
M E = 1 n i = 1 n | T S i T R i | ,
M P E = 1 n i = 1 n | T S i T R i T R i | 100.
where TSi is the simulated value, and TRi is the measured value. By computing ME and MPE, the sample summer days: 23/08/2010, 24/08/2010, and 25/08/2010, MEs and MPEs of the outlet temperature, θOUT, were lower than 0.3 C and 2.1%, respectively. When it comes to the indoor temperature, θi, the ME was 0.7, and the MPE was 2.8%. The reason for this might be the uncertainties about infiltration, which might affect the indoor temperature. Predictions of θOUT were more accurate because the boundary conditions were more suitable to predict.

4.2. Energy Management of the WFG Cabin

When analyzing the WFG cabin, the glazing is part of a room, and the thermal problem of the glazing is coupled with the thermal problem of the room, consequently, the indoor temperature is unknown, and it should be obtained at the same time as the glazing temperature profile. Furthermore, when the system is circulating, the flow rate is set to the design flow rate 0.9 l/min m2 and the inlet temperature, θIN, varies between 16 °C and 20 °C. The outlet temperature, θOUT, is also an unknown and should be obtained, as it is shown in Equation (31). Likewise, when the system is stopped, the flow rate is set to zero, and the outlet temperature of the water chamber, θOUT, is called the stagnation temperature. These outputs are the water heat gain or thermal power of the WFG. Hence, as it is shown in Figure 8, the outlet temperature ranges between 24.5 °C during the day and 15 °C during the night, matching with the simulation curve. Both measured and simulated curves reveal that when the inlet temperature is above the outlet temperature, the WFG is cooled down. This effect takes place during the night and allows to fully grasp the energy strategy that the WFG envelope follows, absorbing solar energy in the water chamber and transferring that heat into the ground through the borehole heat exchangers.
Using Equation (23) and the values taken from Figure 8, the daily energy per unit of surface that the flow of water can transport was calculated. WFG stops this energy from heating up the prototype. It explains the difference in temperatures between the WFG cabin and the Reference cabin. This energy is released through the borehole heat exchangers. Figure 10 shows the energy absorbed by the horizontal WFG. The daily average absorbed energy was 5.8 kWh in 3.5 m2, so the ratio of energy per area is 1.66 kWh/m2. Table 2 illustrates the energy balance between energy absorption and the electric consumption of the circulation water pumps. The circulation pump was working 24 h per day, and its steady power was 25 W. The average daily use of the circulation pump to provide the WFG with a flow rate of 0.9 L per minute was 0.58 kWh, which represents a 10% of the energy absorbed by the water. The net energy savings were 1.49 kWh/m2.
Table 3 compares the final energy consumption of a heat pump with a seasonal coefficient of performance (SCOP) of 3 with the energy consumed by the water pump. The primary energy factor (PEF) from final energy (FE) to non-renewable final energy (NRFE) was 1.954. The factor of CO2 emissions for electricity was 0.331 [44].
WFG allows considering the envelope as a part of the HVAC system. The building envelope is the first shield that receives solar radiation and protects the building against the cold and hot weather. Traditional curtain walls use passive measures, such as coatings and gas cavities. WFG can be used as a component of an active energy management strategy. This article has not focused on analyzing the behavior of the borehole heat exchangers. Further development of this research may analyze the design of an appropriate shallow geothermal system, taking into account the predicted energy absorbed by WFG envelopes.

4.3. Cost Considerations

Borehole heat exchangers are considered low valued and environmentally sustainable energy sources for heating and cooling of buildings that contribute to reducing the primary energy demand. In addition, borehole heat exchangers are considered low exergy (LowEx) systems that allow designing high-performance buildings using the freely available dispersed energy in the environment [45].
The costs of the system made up of WFG and borehole heat exchangers are high compared to passive glazing systems. The cost of boreholes ranges between € 30 and € 50 per meter of drilling. The expected energy yield is 50 W per meter, depending on the characteristics of the soil [11,12]. Figure 10 shows that the average daily energy absorbed by the water is 5.8 kWh over 10 h. The borehole length to dissipate 5.8 kWh is 11.6 m, so the extra cost of this system is € 580 for a 3.5 m2 of WFG. Table 3 shows that the difference in NRFE compared with an air-to-air heat pump is 2.66 Kwh/day. The average price of electricity in Spain is 0.2165 €/Kwh [46], and the daily savings are 0.58 €/day, so the expected return of investment time is less than 3 years.
Furthermore, the commissioning process of the studied system includes preventive and predictive maintenance protocols. The commissioning process is an effective method to ensure that buildings reach their operating potential and that advanced components and systems arrive at their technical specifications. The relationship between investment in commissioning and the cost reduction of the entire project is a topic for further research.

5. Conclusions

Dynamic envelopes, such as WFG, supplied by renewable energy sources like boreholes heat exchangers, can adapt to the building environment instead of creating an insulated barrier. This article has shown that it is possible to simulate dynamic parameters for WFG, such as g factor and U value. Besides, it appears a second thermal transmittance, Uw, associated with the flow of water. Section 2.2.1 provided a mathematical model to calculate U and g-factor as dynamic properties of double WFG. U ranges from 4.8 W/m2 K when the mass flow rate = 0 (UOFF) to 0.66 W/m2 K when = 0.039 Kg/m2s. The g-factor depends on but also the spectral properties of the glass. The g-factor can be modified up to 60% by changing the mass flow rate from = 0 (gOFF) to = 0.039 Kg/m2s.
The potential of WFG as a means to reduce the cooling loads in summer has been tested. A prototype with WFG in facades and the roof, connected to borehole heat exchangers, was built along with another prototype with the same dimensions and made of traditional double glazing. The effect of thermal inertia derived from geothermal boreholes contributes to flatten the temperature curve of the WFG prototype compared with the Reference one. Besides, the non-renewable final energy (NRFE) consumption that is needed to compensate for the absorbed solar radiation is 11.39 Kwh in 3.5 m2 of horizontal WFG. In contrast, the NRFE required to run the circulating pump over the same period is 3.39 Kwh. It represents 70% of NRFE savings. The difference in CO2 emissions between the two options is 1.36 KgCO2. It represents a 70% reduction in CO2 emissions.
The reliability of the mathematical model to simulate the behavior of the prototypes was tested by developing real prototypes. In this paper, the results of the experimental data were compared with the simulation. Two parameters were simulated in the WFG cabin: the indoor air temperature and the outlet temperature in the WFG roof panels. The mean error (ME) and the mean percentage error (MPE) were measured.
  • The ME of the outlet temperature, θOUT, was 0.29 C, and the MPE was 2.1%.
  • The ME of the indoor temperature, θi, was 0.7, and the MPE was 2.8%.
Finally, a number of new topics can be explored in future studies. The results of this article may lead to the following future research suggestions:
  • The performance of WFG in the wintertime.
  • The effect of WFG in comfort by analyzing the mean radiant temperature (MRT).
  • The study of different energy generation systems in buildings (fuel-based or heat pumps) and its potential for energy savings using WFG envelopes.

Author Contributions

Conceptualization, B.M.S., F.d.A.G., and J.A.H.R.; methodology, B.M.S. and F.d.A.G.; software, J.A.H.R.; formal analysis, R.-A.G.-L and F.d.A.G.; data curation, F.d.A.G. and J.A.H.R.; writing—original draft preparation, B.M.S., F.d.A.G., and J.A.H.R.; writing—review and editing, D.P. and B.M.S; visualization, D.P., B.L.A., and R.-A.G.-L.; supervision, J.A.H.R. and B.L.A.; project administration, F.d.A.G. and R.-A.G.-L.; funding acquisition, F.d.A.G. and R.-A.G.-L. All authors have read and agreed to the published version of the manuscript.

Funding

This article has been funded by the KSC Faculty Development Grant (Keene State College, New Hamshire, USA). The authors wish to thank CEU San Pablo University Foundation for the funds dedicated to the Project CEU-Banco Santander (Ref: MVP19V14) provided by CEU San Pablo University and financed by Banco Santander.

Acknowledgments

This work was supported by program Horizon 2020-EU.3.3.1: Reducing energy consumption and carbon footprint by smart and sustainable use, project Ref. 680441 InDeWaG: Industrialized Development of Water Flow Glazing Systems. The authors wish to thank the municipality of Peralveche, Spain, for its generous support.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Sun, Y.; Wilson, R.; Wu, Y. A Review of Transparent Insulation Material (TIM) for building energy saving and daylight comfort. Appl. Energy 2018, 226, 713–729. [Google Scholar] [CrossRef]
  2. Cao, X.; Dai, X.; Liu, J. Building energy-consumption status worldwide and the state-of-the-art technologies for zero-energy buildings during the past decade. Energy Build. 2016, 128, 198–213. [Google Scholar] [CrossRef]
  3. Hermanns, M.; del Ama, F.; Hernández, J.A. Analytical solution to the one-dimensional non-uniform absorption of solar radiation in uncoated and coated single glass panes. Energy Build. 2012, 47, 561–571. [Google Scholar] [CrossRef] [Green Version]
  4. EU. Directive 2010/31/EU of the European Parliament and of the Council of 19 May 2010 on the energy performance of buildings (recast). Off. J. Eur. Union 2010, 153, 13–35. [Google Scholar]
  5. Priarone, A.; Silenzi, F.; Fossa, M. Modelling Heat Pumps with Variable EER and COP in EnergyPlus: A Case Study Applied to Ground Source and Heat Recovery Heat Pump Systems. Energies 2020, 13, 794. [Google Scholar] [CrossRef] [Green Version]
  6. Sanner, B.; Karytsas, C.; Mendrinos, D.; Rybach, L. Current status of ground source heat pumps and underground thermal energy storage in Europe. Geothermics 2003, 32, 579–588. [Google Scholar] [CrossRef]
  7. Vienken, T.; Kreck, M. Dietrich, P. Monitoring the impact of intensive shallow geothermal energy use on groundwater temperatures in a residential neighborhood. Geotherm. Energy 2019, 7, 8. [Google Scholar] [CrossRef]
  8. Boockmeyer, A.; Bauer, S. Efficient simulation of multiple borehole heat exchanger storage sites. Environ. Earth Sci. 2016, 75, 1021. [Google Scholar] [CrossRef]
  9. Bauer, S.; Pfeiffer, T.; Boockmeyer, A.; Dahmke, A.; Beyer, C. Quantifying induced effects of subsurface renewable energy storage. Energy Procedia 2015, 76, 633–641. [Google Scholar] [CrossRef] [Green Version]
  10. Shao, H.; Hein, P.; Sachse, A.; Kolditz, O. Geoenergy Modeling II: Shallow Geothermal Systems; Springer: Berlin, Germany, 2016. [Google Scholar]
  11. Pan, S.; Kong, Y.; Chen, C.; Pang, Z.; Wang, J. Optimization of the utilization of deep borehole heat exchangers. Geotherm. Energy 2020, 8, 6. [Google Scholar] [CrossRef]
  12. Nalla, G.; Shook, G.; Mines, G.; Bloomfield, K. Parametric sensitivity study of operating and design variables in wellbore heat exchangers. Geothermics 2005, 34, 30–46. [Google Scholar] [CrossRef] [Green Version]
  13. Erbs Poulsen, S.; Alberdi-Pagola, M.; Cerra, D.; Magrini, A. An Experimental and Numerical Case Study of Passive Building Cooling with Foundation Pile Heat Exchangers in Denmark. Energies 2019, 12, 2697. [Google Scholar] [CrossRef] [Green Version]
  14. Li, Z.; Zhu, W.; Bai, T.; Zheng, M. Experimental study of a ground sink direct cooling system in cold areas. Energy Build. 2009, 41, 1233–1237. [Google Scholar] [CrossRef]
  15. Li, D.; Wong, S.L.; Tsang, C.L.; Cheung, H.W. A study of the daylighting performance and energy use in heavily obstructed residential buildings via computer simulation techniques. Energy Build. 2006, 38, 1343–1348. [Google Scholar] [CrossRef]
  16. Gueymard, C.; DuPont, W. Spectral effects on the transmittance, solar heat gain, and performance rating of glazing systems. Sol. Energy 2009, 83, 940–953. [Google Scholar] [CrossRef]
  17. Chow, T.T.; Li, C.; Lin, Z. Thermal characteristics of water-flow double-pane window. Int. J. Therm. Sci. 2010, 50, 140–148. [Google Scholar] [CrossRef]
  18. Zevenhoven, R.; Fält, M.; Gomes, L.P. Thermal radiation heat transfer: Including wavelength dependence into modelling. Int. J. Therm. Sci. 2014, 86, 189–197. [Google Scholar] [CrossRef]
  19. Del Ama Gonzalo, F.; Moreno, B.; Hernandez, J.A. Dynamic Solar Energy Transmittance for Water Flow Glazing in Residential Buildings. Int. J. Appl. Eng. Res. 2018, 13, 9188–9193. [Google Scholar]
  20. Fokaides, P.A.; Kylili, A.; Kalogirou, S.A. Phase change materials (PCMs) integrated into transparent building elements: A review. Mater. Renew. Sustain. Energy 2015, 4, 6. [Google Scholar] [CrossRef] [Green Version]
  21. Feijo Munoz, J.; Poza Casado, I.; González-Lezcano, R.A.; Pardal, C.; Echarri, V.; Assiego De Larriva, R.; Fernández-Agüera, J.; Dios-Viéitez, M.J.; Del Campo-Díaz, V.J.; Montesdeoca Calderín, M.; et al. Methodology for the study of the envelope airtightness of residential buildings in Spain: A case study. Energies 2018, 11, 704. [Google Scholar] [CrossRef] [Green Version]
  22. Evangelisti, L.; Battista, G.; Guattari, C.; Basilicata, C.; de Lieto Vollaro, R. Influence of the Thermal Inertia in the European Simplified Procedures for the Assessment of Buildings’ Energy Performance. Sustainability 2014, 6, 4514–4524. [Google Scholar] [CrossRef] [Green Version]
  23. Mateus, N.M.; Pinto, A.; da Graça, G.C. Validation of EnergyPlus thermal simulation of a double skin naturally and mechanically ventilated test cell. Energy Build. 2014, 75, 511–522. [Google Scholar] [CrossRef]
  24. Torregrosa-Jaime, B.; Martínez, P.J.; González, B.; Payá-Ballester, G. Modelling of a Variable Refrigerant Flow System in EnergyPlus for Building Energy Simulation in an Open Building Information Modelling Environment. Energies 2019, 12, 22. [Google Scholar] [CrossRef] [Green Version]
  25. Nardi, I.; Lucchi, E.; de Rubeis, T.; Ambrosini, D. Quantification of heat energy losses through the building envelope: A state-of-the-art analysis with critical and comprehensive review on infrared thermography. Build. Environ. 2018, 146, 190–205. [Google Scholar] [CrossRef] [Green Version]
  26. Bambardekar, S.; Poerschke, U. The architect as performer of energy simulation in the early design stage. In Proceedings of the IBPSA 2009—International Building Performance Simulation Association, Glasgow, Scotland, 27–30 July 2009; pp. 1306–1313. [Google Scholar]
  27. Fernandez-Antolin, M.; del-Rio, J.M.; del Ama Gonzalo, F.; Gonzalez-Lezcano, R. The Relationship between the Use of Building Performance Simulation Tools by Recent Graduate Architects and the Deficiencies in Architectural Education. Energies 2020, 13, 1134. [Google Scholar] [CrossRef] [Green Version]
  28. Echarri, V.; Espinosa, A.; Rizo, C. Thermal Transmission through Existing Building Enclosures: Destructive Monitoring in Intermediate Layers versus Non-Destructive Monitoring with Sensors on Surfaces. Sensors 2017, 17, 2848. [Google Scholar] [CrossRef] [Green Version]
  29. Pires, I.M.; Garcia, N.M.; Pombo, N.; Flórez-Revuelta, F.; Rodríguez, N.D. Validation Techniques for Sensor Data in Mobile Health Applications. J. Sens. 2016, 2016, 2839372. [Google Scholar] [CrossRef] [Green Version]
  30. Abdulmohsin, H.; Aritra, G.; Senthilarasu, S.; Tapas, K.M. Evaluation of thermal performance for a smart switchable adaptive polymer dispersed liquid crystal (PDLC) glazing. Sol. Energy 2020, 195, 185–193. [Google Scholar] [CrossRef]
  31. Allen, K.; Connelly, K.; Rutherford, P.; Wu, Y. Smart windows—Dynamic control of building energy performance. Energy Build. 2017, 139, 535–546. [Google Scholar] [CrossRef]
  32. Ghosh, A.; Norton, B.; Duffy, A. Measured overall heat transfer coefficient of a suspended particle device switchable glazing. Appl. Energy 2015, 159, 362–369. [Google Scholar] [CrossRef]
  33. Ghosh, A.; Norton, B.; Duffy, A. Behaviour of a SPD switchable glazing in an outdoor test cell with heat removal under varying weather conditions. Appl. Energy 2016, 180, 695–706. [Google Scholar] [CrossRef] [Green Version]
  34. Casini, M. Smart windows for energy efficiency of buildings. Int. J. Civ. Struct. Eng. IJCSE 2015, 2, 230–238. [Google Scholar] [CrossRef]
  35. Moreno, B.; Del Ama Gonzalo, F.; Ferrandiz, J.A.; Hernandez, J.A. A Building energy simulation methodology to validate energy balance and comfort in zero energy buildings. J. Energy Syst. 2019, 3, 168–182. [Google Scholar] [CrossRef]
  36. Sierra, P.; Hernandez, J.A. Solar heat gain coefficient of water flow glazing. Energy Build. 2017, 139, 133–145. [Google Scholar] [CrossRef] [Green Version]
  37. Chow, T.T.; Li, C. Liquid-filled solar glazing design for buoyant water-flow. Build. Environ. 2013, 60, 45–55. [Google Scholar] [CrossRef]
  38. Chow, T.T.; Li, C.; Lin, Z. Innovative solar windows for cooling-demand climate. Sol. Energy Mater. Sol. Cells 2010, 94, 212–220. [Google Scholar] [CrossRef]
  39. EN 673, Glass in Building. In Determination of Thermal Transmittance (U Value). Calculation Method; British Standards Institution: London, UK, 2011.
  40. EN 410, Glass in Building. In Determination of Luminous and SolarCharacteristics of Glazing; British Standards Institution: London, UK, 2011.
  41. ISO 18292. Energy Performance of Fenestration Systems for Residential Buildings—Calculation Procedure; ISO: Geneva, Switzerland, 2011. [Google Scholar]
  42. Li, C.; Chow, T.T. Water-filled double reflective window and its year-roundperformance. Proc. Environ. Sci. 2011, 11, 1039–1047. [Google Scholar] [CrossRef] [Green Version]
  43. Finlayson, E.; Arasteh, D.; Huizenga, C.; Rubin, M.; Reilly, M. WINDOW 4.0:Documentation of Calculation Procedures; University of California, LawrenceBerkeley Laboratory: Berkeley, CA, USA, 1993. [Google Scholar]
  44. Spanish Regulation of Thermal Installations in Buildings (RITE). Factores de Emisión de CO2 y Coeficientes de paso a Energía Primaria de Diferentes Fuentes de Energía Final Consumidas en el Sector de Edificios en España; Instituto para la Diversificación y Ahorro de la Energía (IDAE): Madrid, Spain, 2016. [Google Scholar]
  45. Meggers, F.; Ritter, V.; Goffin, P.; Baetschmann, M.; Leibundgut, H. Low exergy building systems implementation. Energy 2012, 41, 48–55. [Google Scholar] [CrossRef]
  46. Statistical Requirements Compendium; Eurostat: Luxembourg, 2019.
Figure 1. Schematics of the energy system. 1. Plate heat exchangers, 2. Water circulation pumps and precision flow meters, 3. Water flow glazing (WFG), 4. Borehole heat exchangers.
Figure 1. Schematics of the energy system. 1. Plate heat exchangers, 2. Water circulation pumps and precision flow meters, 3. Water flow glazing (WFG), 4. Borehole heat exchangers.
Energies 13 03203 g001
Figure 2. Real prototypes on their location.
Figure 2. Real prototypes on their location.
Energies 13 03203 g002
Figure 3. (A) Temperature distribution in a WFG and description of heat fluxes and temperatures of layers at a specific height of the envelope. (B) Distribution of solar radiation on the glazing. It is divided into the reflectance (R), the energy transmittance (T), the constant part of the secondary external and internal heat transfer factor (Ae and Ai respectively), and the net absorptance of the water (Av).
Figure 3. (A) Temperature distribution in a WFG and description of heat fluxes and temperatures of layers at a specific height of the envelope. (B) Distribution of solar radiation on the glazing. It is divided into the reflectance (R), the energy transmittance (T), the constant part of the secondary external and internal heat transfer factor (Ae and Ai respectively), and the net absorptance of the water (Av).
Energies 13 03203 g003
Figure 4. (a) U value of WFG depending on the mass flow rate, . (b) g-factor of WFG depending on the mass flow rate, , with three different glazing types.
Figure 4. (a) U value of WFG depending on the mass flow rate, . (b) g-factor of WFG depending on the mass flow rate, , with three different glazing types.
Energies 13 03203 g004
Figure 5. Measured data of the outdoor temperature and the indoor temperature of the WFG cabin and Reference cabin. Sample summer days: 23/08/2010, 24/08/2010, and 25/08/2010.
Figure 5. Measured data of the outdoor temperature and the indoor temperature of the WFG cabin and Reference cabin. Sample summer days: 23/08/2010, 24/08/2010, and 25/08/2010.
Energies 13 03203 g005
Figure 6. Measured data of the inlet and outlet temperature of the WFG cabin. Sample summer days: 23/08/2010, 24/08/2010, and 25/08/2010.
Figure 6. Measured data of the inlet and outlet temperature of the WFG cabin. Sample summer days: 23/08/2010, 24/08/2010, and 25/08/2010.
Energies 13 03203 g006
Figure 7. Solar irradiance. Comparison results of simulation data and measured data of the solar irradiance of the prototype. Sample summer days: 23/08/2010, 24/08/2010, and 25/08/2010.
Figure 7. Solar irradiance. Comparison results of simulation data and measured data of the solar irradiance of the prototype. Sample summer days: 23/08/2010, 24/08/2010, and 25/08/2010.
Energies 13 03203 g007
Figure 8. Comparison between measured data of the inlet and outlet temperature of the WFG cabin and results from simulation data. Sample summer days: 23/08/2010, 24/08/2010, and 25/08/2010.
Figure 8. Comparison between measured data of the inlet and outlet temperature of the WFG cabin and results from simulation data. Sample summer days: 23/08/2010, 24/08/2010, and 25/08/2010.
Energies 13 03203 g008
Figure 9. Performance of WFG. Room indoor temperature in the WFG cabin and outdoor temperature. Comparison between simulation results calculated with the tool and real data from the prototype. Sample summer days: 23/08/2010, 24/08/2010, and 25/08/2010.
Figure 9. Performance of WFG. Room indoor temperature in the WFG cabin and outdoor temperature. Comparison between simulation results calculated with the tool and real data from the prototype. Sample summer days: 23/08/2010, 24/08/2010, and 25/08/2010.
Energies 13 03203 g009
Figure 10. Energy absorbed in the horizontal glazing. Area of glazing = 3.5 m2; mass flow rate = 3.2 L per minute. Sample summer days: 23–26 August.
Figure 10. Energy absorbed in the horizontal glazing. Area of glazing = 3.5 m2; mass flow rate = 3.2 L per minute. Sample summer days: 23–26 August.
Energies 13 03203 g010
Table 1. Thermal and spectral properties of glazing.
Table 1. Thermal and spectral properties of glazing.
A1A2AwU U w O N     1 ggON 1
Reference Glazing0.5850.037 2.6-0.67-
WFG0.5850.0370.429–0.5240.7625.802-0.544
1 These values have been obtained with the hypothesis: ṁc >> Ue + Ui.
Table 2. Energy balance of the horizontal WFG.
Table 2. Energy balance of the horizontal WFG.
Absorbed Energy, kWh (Day Time)Dissipated Energy, kWh (Night Time)Water Pump, kWh (24 h)
Day 15.421.470.58
Day 25.841.550.58
Day 36.231.510.58
Table 3. Final energy analysis. Sample summer days: 23/08/2010, 24/08/2010, and 25/08/2010.
Table 3. Final energy analysis. Sample summer days: 23/08/2010, 24/08/2010, and 25/08/2010.
WFG Absorbed Energy (3-Day Time)Water Pump
(3-Day Time)
Energy consumption17.491.74
SCOP3
FE consumption (Kwh)5.831.74
NRFE consumption (Kwh)5.83 × 1.954 = 11.391.74 × 1.954 = 3.39
CO2 emissions (KgCO2)5.83 × 0.331 = 1.931.74 × 0.331 = 0.57

Share and Cite

MDPI and ACS Style

Moreno Santamaria, B.; del Ama Gonzalo, F.; Pinette, D.; Gonzalez-Lezcano, R.-A.; Lauret Aguirregabiria, B.; Hernandez Ramos, J.A. Application and Validation of a Dynamic Energy Simulation Tool: A Case Study with Water Flow Glazing Envelope. Energies 2020, 13, 3203. https://doi.org/10.3390/en13123203

AMA Style

Moreno Santamaria B, del Ama Gonzalo F, Pinette D, Gonzalez-Lezcano R-A, Lauret Aguirregabiria B, Hernandez Ramos JA. Application and Validation of a Dynamic Energy Simulation Tool: A Case Study with Water Flow Glazing Envelope. Energies. 2020; 13(12):3203. https://doi.org/10.3390/en13123203

Chicago/Turabian Style

Moreno Santamaria, Belen, Fernando del Ama Gonzalo, Danielle Pinette, Roberto-Alonso Gonzalez-Lezcano, Benito Lauret Aguirregabiria, and Juan A. Hernandez Ramos. 2020. "Application and Validation of a Dynamic Energy Simulation Tool: A Case Study with Water Flow Glazing Envelope" Energies 13, no. 12: 3203. https://doi.org/10.3390/en13123203

APA Style

Moreno Santamaria, B., del Ama Gonzalo, F., Pinette, D., Gonzalez-Lezcano, R. -A., Lauret Aguirregabiria, B., & Hernandez Ramos, J. A. (2020). Application and Validation of a Dynamic Energy Simulation Tool: A Case Study with Water Flow Glazing Envelope. Energies, 13(12), 3203. https://doi.org/10.3390/en13123203

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop