Energy Performances Assessment of Extruded and 3D Printed Polymers Integrated into Building Envelopes for a South Italian Case Study

: Plastic materials are increasingly becoming used in the building envelope, despite a lack of investigation on their effects. In this work, an extruded Acrylonitrile-Butadiene-Styrene panel has been tested as a second-skin layer in a ventilated facade system using a full-scale facility. The experimental results show that it is possible to achieve performances very similar to conventional materials. A numerical model has then been developed and used to investigate the performances of plastic and composite polymer panels as second-skin layers. The experimental data has been used to verify the behavior of the numerical model, from a thermal point of view, showing good reliability, with a root mean square error lower than 0.40 ◦ C. This model has then been applied in different refurbishment cases upon varying: the polymer and the manufacturing technology (extruded or 3D-printed panels). Eight refurbishment case studies have been carried out on a typical ofﬁce building located in Napoli (Italy), by means of a dynamic simulation software. The simulation results show that the proposed actions allow the reduction of the thermal and cooling energy demand (up to 6.9% and 3.1%, respectively), as well as the non-renewable primary energy consumption (up to 2.6%), in comparison to the reference case study. the T cavity is equal to 0.3 °C; (ii) the values of RMSE are equal to 0.4 °C and 0.2 °C for T indoor and T cavity , respectively. Therefore, the results show the ability of the Type 1230 to accurately predict the be- havior of the extruded ABS panels in an OVF system. Thus, the same methodology is used to carry out a complete numerical campaign on a set of case studies upon varying the polymer (selecting the ones more used in architecture, as highlighted in the literature re- view) and technology (extruded and


Introduction
Approximately 40% of the EU energy consumption can be directly attributed to the building sector, which is also responsible for about 36% of the greenhouse gas emissions [1,2]. In addition, in the EU-28, only 3% of the edifices have an efficient building envelope [3], mainly due to the fact that about 35% of the EU's buildings are over 50 years old and only around 1% of them are renovated each year [1]. Certainly, the constraints associated with the new buildings are fewer with respect to those associated with the refurbishment of existing constructions, so the new one allows for better-optimized design in terms of energy efficiency of the envelope [4]. However, in Italy, many buildings (about 4 million) were built in the early 1900s and about half of these have been classified as historical architectures and nowadays have been reused [5]. Therefore, in the Italian scenario, the improvement of the energetic performances of the existing building envelope represents a crucial aspect in the increasing of the building's energy efficiency and the indoor environmental quality on a large-scale [6]. In this context, different products and systems have been proposed to improve the buildings energy efficiency, visual and thermal comfort, as well as their sustainability [7][8][9][10][11][12] and, in recent years, the interest of the scientific community has seen an increase in the facade domain to improve the overall building energy efficiency [13]. In particular, the use of passive systems is raising more

Plastic and Composite Polymers in Building Facades
In recent years, the use of polymers in building and engineering has increased substantially, thanks to their: (i) ease of production, (ii) ease of installation, (iii) durability, (iv) low maintenance requirements, (v) lightweight nature and (vi) ability to be formed into complex shapes [34,35]. Moreover, polymers form good thermal and electrical insulators that are not affected by chemical and biological risks [35]. They have not only been used to replace the traditional construction materials (cement, brick, concrete, wood, metal, and glass), but these materials have also been used in a complementary way to improve the building envelope performance to satisfy the modern demands of both new projects and refurbishment ones [32][33][34][35]. From an aesthetic point of view, these materials are available in several colors and texture alternatives. Several applications of plastic and composite polymer walls in buildings were reported in the literature [33]. In [33], 23 examples of architecture were reported, demonstrating that plastic and composite polymers can be used in dwelling domes, large-span volumes or envelop large facade surfaces, and transparent sky-lights on the roofs of industrial buildings. Another example can be found in [36,37], where the designers have realized a temporary pavilion with an envelope made of double-walled transparent corrugated sheets of PolyEthylene Terephthalate (PET) recycled plastic. Similarly, in [38], a massive pavilion, designed as an exhibition hall for the 2010 Taipei International Flora Expo, has been built using recycled PET for the building envelope, also proving strong resistance to fires and earthquakes. As reported in [39], polycarbonate multi-sheet systems are increasing their share of the glazing market since they provide good performance while weighing and costing significantly less than glass. For these reasons, several studies have been conducted to assess their characteristics upon varying geometries and installation typology. In [39], the authors assessed the thermal and optical characteristics of different polycarbonate sheets, highlighting a strong angular dependence in polycarbonate panels' optical properties, significantly different to the conventional multiple-layered glass, and a good overall thermal behavior. For these reasons, the polycarbonate could be used as a valid alternative, fulfilling the energy requirements and improving the visual comfort, reducing the glaring problem by diffusing the light, while providing more flexibility in the design and appearance of the buildings. In [40], a detailed experimental analysis of the thermal behavior of different polycarbonate multi-sheet systems has been carried out by varying the installation angle. The results highlighted a very low incidence of the angle of inclination on the equivalent value of thermal conductivity, thus allowing the material to be equally used in every part of the building envelope. Thus, polycarbonate has already found large usage, as in the Bavaria Brewery Tocancipá headquarters, by Construcciones Planificadas [41], where the plastic has been used to get an industrial appearance while providing for thermal and visual comfort, or the property registration offices, by Irisarri + Piñera, where polycarbonate has been used to complement and balance the appearance of the existing structure.
Across these materials already implemented in the traditional architectural paradigms, there is also a strong boost in the usage of plastics from the additive manufacturing (AM) sector, as several plastic products can also be utilized in AM processes, providing great freedom of form and enhancing designers', architects', and engineers' freedom in creating complex designs [42][43][44][45][46]. In addition, the global 3D printing filament material market volume was 1.8 billion US$ in 2019, growing at a compound annual growth rate (CAGR) of 27% [47]. The most popular filament materials are PolyLactic Acid (PLA) and Acrylonitrile-Butadiene-Styrene (ABS), holding about 47% and 29% of the market, respectively. In contrast, materials as PolyEthylene Terephthalate Glycol-modified (PETG) can be useful alternatives, despite not being as popular, providing similar mechanical properties while also offering excellent recyclability and scalability [47]. The AM in the facade industry presents new significant potential and requires relevant research to be conducted [46]. Indeed, more and more 3D printing materials are concurrently becoming utilized in contemporary architecture design [42][43][44][45][46], thanks to the lightness and effortless installation procedure, which results in a design solution useful for both new projects and refurbishment ones [44][45][46]. The 3D printing technology has been often used to create everything, from prototypes [42,43,46] and simple parts of facades, to give a distinctive signature to the constructions [44,45]. The 3D printing materials prove themselves to offer quite unique characteristics from an architectural and economic point of view [43][44][45][46]. Several scientific papers have been conducted into the loadbearing capacities and/or other essential qualities of AM products for the building industry, such as durability, water vapor diffusion resistance, thermal conductivity, or fire-resistance [48][49][50][51]. The authors emphasized the potential applications of additive manufacturing to build honeycomb panels that optimize mechanical properties and heat transfer [48][49][50][51]. However, scientific research related to 3D printing materials in building energy efficiency applications is limited due to its relatively new nature as a technology. Only Sarakinioti et al. [46] aimed their research at developing an integrated 3D printed for thermal insulation and building physics. In particular, they presented a 3D printed facade panel design for thermal insulation and movable liquid heat storage [46], providing an overview of the development process. The authors tested the prototype and, at the same time, simulated to verify the thermal effects of the proposed facade system on indoor spaces in different climates. The simulation results reported in [46] showed the potential of the proposed 3D printed facade panel for reducing heating and cooling energy demand. Therefore, the effects of adopting 3D printing materials as a second-skin layer on the indoor environment have been scarcely investigated.
Moreover, there is a lack of experimental testing and numerical model development of these materials in building simulation, even more, if considered in a second-skin in front of the building envelope, in a 3D printed composite facade arrangement. Indeed, in a facade arrangement realized with these innovative materials (ABS, PLA, PETG, etc.), the difficulties lie in predicting the behavior of the various facade sections, as the second-skin layer, the resulting air cavity, and finally, the effects on the indoor environment. Therefore, from an experimental point of view, standardization bodies, experts, and researchers are continually developing new methodologies or new procedures to correctly calculate the performances of these envelope components in simple and economical ways [6].

Research Aims
In this work, extruded ABS panels have been tested as second-skin materials in order to verify their performances in an OVF system. This novel material for building envelopes has been investigated through in-situ measurements by using two outdoor comparative test cells. The experimental data have been used to calibrate and validate a numerical model in TRNSYS 18 [31], also verifying the ability of the simulation software to effectively reproduce the behavior of a light material in an OVF system, which usually is made of materials as porcelain gres. Then, the validated numerical methodology has been used to implement different plastic materials in a set of refurbishment case studies, compared to a reference office building, in order to assess the potential benefits. The comparison has been performed in terms of (i) heating and (ii) cooling energy demands, as well as (iii) non-renewable primary energy consumption, upon varying the plastic material. Finally, additional refurbishment case studies have also been implemented considering 3D printed panels as a second-skin layer.
The aims of this research can be summarized as reported below: • investigate the performances of extruded ABS panels as a second-skin layer for innovative building envelopes with experimental tests in-situ; • calibrate and validate a simulation model to predict the energy performance of the plastic and composite polymer panels used as a second-skin layer in an OVF system; • assess the potential energy saving achievable in office building refurbishment using the proposed materials (extruded and 3D printed polymers) through numerical simulation.

Structure of the Research
The research is structured as follows. Section 2 describes the methodology used to carry out the research, showing in detail how (i) the experimental data have been acquired, (ii) the numerical model has been implemented in TRNSYS 18, and how (iii) the experimental data have been used to calibrate and validate this numerical model. Section 3 reports the numerical results, in terms of the reduction of non-renewable primary energy consumption, achieved in an office building refurbishment through the installation of an OVF system, upon varying the material used as a second-skin layer, considering both extruded and 3D printed ones.
Finally, Section 4 discusses the integration capacity of the plastic and composite polymer panels in a second-skin layer of an OVF system, highlighting the advantages and limitations of such materials.

Methodology
This section describes in detail the measurement methodologies and the experimental results obtained during the in-situ test as well as the methods and results related to the validation of the implemented numerical model.

Description of the Gemini Facilities, Experimental Results, and Discussion
In this sub-section, a couple of experimental test cells and the experimental results are reported. Gemini facilities [11]  with the long side along the north-south axis. These dimensions correspond to the gross dimension of the main steel frame structure, on which the shell has been fixed externally and seamlessly in order to avoid thermal bridges. The shell has been realized in a single layer of 10 mm thick sandwich panels consisting of two galvanized steel sheets and a polyurethane rigid foam filling, with a thermal transmittance (U wall ) value of 0.23 W/m 2 K [11]. Then, for the floor, a 0.10 m air gap and a wood flooring have been added above the structure, while, for the ceiling, a sheet metal roof has been placed 0.10 m above the outer panels, with a 2% slope, to allow a natural rainwater outflow.
The Gemini's facilities are designed to allow the in-situ characterization of innovative layers to be applied in double-skin facades with different geometries, layout, materials, and technologies. The acquired data can be used to evaluate the in-situ performances of the system under investigation and to realize, calibrate and validate simulation models.
The Gemini is well-instrumented to acquire different indoor and outdoor physical quantities. Table 1 shows the measurement range, the accuracy, and the response time of the sensors used for outdoor and indoor climate characterization. In particular, with the aim to evaluate the real weather conditions, sensors for wind direction, wind speed, air temperature, air relative humidity, air pressure, global horizontal radiation, and diffuse horizontal radiation were placed at about 6.50 m from the ground, in the best position to minimize the influence of external obstructions (i.e., the obstructions angles are less than 10 • ). In order to acquire diffuse horizontal radiation, one of the pyranometers is equipped with a shadow ring (diameter of 0.574 m and thickness equal to 0.052 m), and the data were corrected following the methodology proposed in [53], to take into account both the isotropic and anisotropic conditions. Figure 1 shows the weather station, with all the aforementioned sensors. The air temperature inside the Gemini is monitored by a combined temperature-relative humidity sensor placed in the middle of the room. Also, when the test cells are configured to test a second-skin system, this is monitored through a set of ten thermocouples, placed on the significant interfaces and in the cavity. The sensor layout follows the layout shown in Figure 2a. In particular, it can be noted that (i) six thermocouples are placed in the middle of both the inlet (T1, T2, and T3) and the outlet (T8, T9, and T10) sections of the air cavity, (ii) four thermocouples are placed in line at the center of the second-skin system (T4 on the back surface of the second-skin, T5 in the middle of the cavity, T6 on the external surface of the south wall of the test cell, T7 on the internal surface of the south wall of the test cell, respectively); this last set of thermocouples falls in line with the same thermo-hygrometer which monitors the temperature of the air inside the test cell, in order to have all the sensors aligned at the center of the system.
In addition, all the thermocouples (Tx) have been shielded with high-reflective domes in order to avoid any direct solar radiation (Figure 2b).
The pyranometer (Ivert) has been installed to acquire the vertical solar radiation incident on the south surface and, finally, two hot-wire anemometers (Win, and Wout) are placed in the air cavity, one in the inlet section and the other one in the outlet section, in order to monitor the airflow in the second-skin cavity. The air temperature inside the Gemini is monitored by a combined temperaturerelative humidity sensor placed in the middle of the room. Also, when the test cells are configured to test a second-skin system, this is monitored through a set of ten thermocouples, placed on the significant interfaces and in the cavity.
The sensor layout follows the layout shown in Figure 2a. In particular, it can be noted that (i) six thermocouples are placed in the middle of both the inlet (T1, T2, and T3) and the outlet (T8, T9, and T10) sections of the air cavity, (ii) four thermocouples are placed in line at the center of the second-skin system (T4 on the back surface of the second-skin, T5 in the middle of the cavity, T6 on the external surface of the south wall of the test cell, T7 on the internal surface of the south wall of the test cell, respectively); this last set of thermocouples falls in line with the same thermo-hygrometer which monitors the temperature of the air inside the test cell, in order to have all the sensors aligned at the center of the system. In order to verify the measurement methodologies and characterize each test cell from the thermal point of view, preliminary data has been acquired in a standard configuration (both Gemini s without a second-skin system). These data are recorded with the aim of (i) verifying the operation of the different instruments and their correct positioning, (ii) comparing the thermal behavior of the two test cells, and (iii) defining a reference point for the following evaluation of the real performances of double-skin facades or smart win- In addition, all the thermocouples (Tx) have been shielded with high-reflective domes in order to avoid any direct solar radiation (Figure 2b).
The pyranometer (I vert ) has been installed to acquire the vertical solar radiation incident on the south surface and, finally, two hot-wire anemometers (W in , and W out ) are placed in the air cavity, one in the inlet section and the other one in the outlet section, in order to monitor the airflow in the second-skin cavity.
In order to verify the measurement methodologies and characterize each test cell from the thermal point of view, preliminary data has been acquired in a standard configuration (both Gemini s without a second-skin system). These data are recorded with the aim of (i) verifying the operation of the different instruments and their correct positioning, (ii) comparing the thermal behavior of the two test cells, and (iii) defining a reference point for the following evaluation of the real performances of double-skin facades or smart windows. The preliminary experimental data were acquired and stored every 1 min on a period of 1 month (from 1 June to 30 June) and, later, averaged on an interval of 15 min. Figure 3 reports the indoor air temperature of Gemini 1 and Gemini 2, the external air temperature, and the global horizontal radiation for three typical days in June. This figure highlights that the difference between the indoor air temperature of Gemini 1 and the indoor air temperature of Gemini 2 (T indoor, Gemini 1 -T indoor, Gemini 2 ) is negligible, varying within an interval with a maximum of 0.2 • C and a minimum of −0.2 • C.  After the preliminary experimental campaign, a second-skin system has been mounted and tested on a test cell (Gemini 1) with an air cavity gap equal to 0.10 m, while the other cell has been left unequipped and used as a reference (Gemini 2). In particular, the investigated second-skin system has been realized with extruded ABS panels [54]. The experimental data were acquired and stored every 1 min over a period of 1 month (from 8 December to 31 December) and, later, averaged on an interval of 15 min. This experimental campaign aims to verify the performances of plastic panels as a second-skin layer for innovative envelopes.
The extruded ABS panels have been selected with dimensions equal 600 mm × 1200 mm; such dimensions have been selected on the basis of the dimension of conventional panels in OVF systems, in order to guarantee an easy installation in a commercial OVF structure, as well as an easy substitution of conventional OVF system materials in a realistic scenario. In Figure 4a, a detailed view of the analyzed extruded ABS panel before the installation is displayed, while Figure 4b shows the Gemini 1 equipped with the second-skin realized with six extruded ABS panels.  The second-skin system equipped during the tests has been realized by mounting the six extruded ABS panels on a steel frame, then hanging the whole system to the brackets on the south facade of the test cell Gemini 1. Finally, the sides of the second-skin system were covered and sealed by means of panels, similar to those used for the test cells' enve- The second-skin system equipped during the tests has been realized by mounting the six extruded ABS panels on a steel frame, then hanging the whole system to the brackets on the south facade of the test cell Gemini 1. Finally, the sides of the second-skin system were covered and sealed by means of panels, similar to those used for the test cells' envelope, in order to allow only for a vertical airflow in the cavity (Figure 4b). In this configuration, the acquisition period lasted for almost a month, in the wintertime. During the acquisition period, the temperatures were monitored following the layout reported in Figure 2a. Figure 5 shows an overview of the whole acquisition period, reporting the external air temperature, the global horizontal radiation, and the total vertical radiation on the south facade.      Figure 5 highlights that the external air temperature was quite warm, despite being the winter season. Also, the radiation values, both global horizontal and total vertical, show mostly high values, thus confirming good weather and clear sky across the whole acquisition period. Figure 6a,b report a focus for four typical days on the weather conditions during the measurements with the Gemini 1 equipped with the second-skin system, and the Gemini 2 left uncovered as a reference. In particular, Figure 6a shows the temperature and solar radiation data (global horizontal solar radiation, diffuse horizontal solar radiation, and total vertical solar radiation on the south facade), while Figure 6b reports the wind characteristics acquired during the analyzed days; on the left axis the wind speed is reported, while on the right axis the wind direction is displayed, considering 0 • as north direction, 90 • as east direction, 180 • as south direction, and 270 • as west direction.  Figure 5 highlights that the external air temperature was quite warm, despite being the winter season. Also, the radiation values, both global horizontal and total vertical, show mostly high values, thus confirming good weather and clear sky across the whole acquisition period. Figure 6a,b report a focus for four typical days on the weather conditions during the measurements with the Gemini 1 equipped with the second-skin system, and the Gemini 2 left uncovered as a reference. In particular, Figure 6a shows the temperature and solar radiation data (global horizontal solar radiation, diffuse horizontal solar radiation, and total vertical solar radiation on the south facade), while Figure 6b reports the wind characteristics acquired during the analyzed days; on the left axis the wind speed is reported, while on the right axis the wind direction is displayed, considering 0° as north direction, 90° as east direction, 180° as south direction, and 270° as west direction.  Figure 6a better highlights that, during the measurement period, sunshine days were acquired with atypical temperatures for the period, ranging between a minimum of about 5.9 • C and a maximum of about 19.4 • C. Figure 6b shows a low wind speed in general during the measurement period and a slight wind predominance in the west/north-west direction. Also, the wind speed values acquired during the nighttime are equal to zero, due to a threshold value for the start/stop of the sensor equal to 0.15 m/s. Figures 7 and 8 report the experimental data associated with the cavity with a 1-h timestep for a single acquisition day. In more detail, Figure 7 shows the trends of the daily values of the air temperature inside the cavity upon varying the height from the ground, while Figure 8 reports the airspeed at the inlet and the outlet of the air cavity, for the same day.
The data reported in Figure 7 corresponds to the measures of the three thermocouples positioned in the middle point of the cavity of the second-skin system, more specifically T2, in the middle of the air cavity inlet, T5, in the middle of the air cavity geometrical center, and T9, in the middle of the air cavity outlet, as shown in Figure 2a. As a first observation, the overall temperature distribution is directly related to solar radiation throughout the day, where the temperatures rise during the morning and drop during the afternoon. Also, during the day, the temperature trend seems to be substantially constant from the air cavity inlet to the middle of the facade, and then to increase to the air cavity outlet; this behavior is due to the chimney effect that is created thanks to the OVF system. Then, in the evening (starting from 16.00), the temperatures at the outlet of the cavity are only slightly higher than those at the center and the inlet; this is due to the reduction of solar radiation happening in the evening. of the sensor equal to 0.15 m/s. Figures 7 and 8 report the experimental data associated with the cavity with a 1-h timestep for a single acquisition day. In more detail, Figure 7 shows the trends of the daily values of the air temperature inside the cavity upon varying the height from the ground, while Figure 8 reports the airspeed at the inlet and the outlet of the air cavity, for the same day.  The data reported in Figure 7 corresponds to the measures of the three thermocouples positioned in the middle point of the cavity of the second-skin system, more specifically T2, in the middle of the air cavity inlet, T5, in the middle of the air cavity geometrical center, and T9, in the middle of the air cavity outlet, as shown in Figure 2a. As a first observation, the overall temperature distribution is directly related to solar radiation throughout the day, where the temperatures rise during the morning and drop during the  of the sensor equal to 0.15 m/s. Figures 7 and 8 report the experimental data associated with the cavity with a 1-h timestep for a single acquisition day. In more detail, Figure 7 shows the trends of the daily values of the air temperature inside the cavity upon varying the height from the ground, while Figure 8 reports the airspeed at the inlet and the outlet of the air cavity, for the same day.  The data reported in Figure 7 corresponds to the measures of the three thermocouples positioned in the middle point of the cavity of the second-skin system, more specifically T2, in the middle of the air cavity inlet, T5, in the middle of the air cavity geometrical center, and T9, in the middle of the air cavity outlet, as shown in Figure 2a. As a first observation, the overall temperature distribution is directly related to solar radiation throughout the day, where the temperatures rise during the morning and drop during the    Figure 2a, where W in measured the airspeed value at the inlet of the cavity and W out measured the value at the outlet of the cavity. This figure shows how, during the day, the values acquired by W out is higher than those measured by W in , referable to a direct effect of the solar radiation on the rising of the air temperatures along the air cavity, thus causing an increase in airspeed. After 15:00, the chimney effect in the air cavity is reduced because the temperatures are gradually decreasing over time due to the reduction of the solar radiation on the south facade; this causes a significant drop in the airspeed values measured in the cavity.
Therefore, the analysis of the experimental results shows that a plastic material (i.e., ABS) can be used as a second-skin layer in OVF systems.

Description of the Numerical Model
The software TRNSYS 18 [31] has been used to model the Gemini test cells and to develop the second-skin model. TRNSYS software adopts a modular approach by using Fortran subroutines. Each Fortran subroutine is called a 'Type' and contains the model for a single system component. Several studies have been carried out in order to validate the numerical models developed in TRNSYS from the Colorado State University experimental houses and other researchers around the world [55][56][57][58]. In this study, the following main TRNSYS Types [31,[59][60][61]  Type 33e to determine the moist air properties [59,60].
At first, the Type 56 subroutine has been used to model the thermal behavior of the Gemini test cells; in particular, two thermal zones have been modeled, one for each test cell. Also, Type 56 contains information about the test cells' surroundings (buildings, trees, bushes), which are described as 'shading objects'. The geometrical modeling occurred in the SketchUp software [62], where it was possible to model the shapes and the position of each element accurately. Then, by means of the Trnsys3D plug-in, the geometries were imported into the Type 56 subroutine. The physical properties of each test cell's external surface have been defined, on the basis of the data provided by the manufacturers. Lastly, the internal thermal gains have been set for each test cell, which was determined on the basis of the equipment installed inside the facilities (notebook, data acquisition systems, and uninterruptible power supply units).
The OVF system has been modeled using the Type 1230 subroutine, which effectively reproduces the behavior of an external second-skin layer with an air cavity behind it. Using this TRNSYS Type, the behavior of the OVF system has been correlated to that of the Gemini test cell modeled through the Type 56 subroutine. In particular, the last external layer of the Type 56 wall acts as an interface layer between the Type 1230 and the Type 56, by coupling its temperature and thermal resistance to model the wall heat transfer. Figure 9 shows a schematic of the boundaries of the two coupled Types (56 and 1230), highlighting the resistive interface layer. The Type 1230 parameters have been set following the data provided by the manufacturer of the extruded ABS panels [54], taking into account thickness, density, and thermal conductivity, specifically.
During the simulations, Type 1230 takes into account: • the solar radiation, the longwave radiation, and the air convection on the external surface of the outside layer; • the energy storage and the conduction in the outside layer; • radiation exchange between the outside layer and the air cavity; • the convective exchanges from all the surfaces facing in the air cavity; • the conduction through the interface layer.
Type 16c has been implemented to model the solar radiation on all the external surfaces. This Type accepts global radiation, ambient temperature, and ambient relative humidity data as input, in order to output several quantities related to the position of the sun, as the diffuse radiation fraction on the horizontal, by estimating the cloudiness of the sky on the basis of the dry bulb temperature and the dew point temperature. Finally, the radiation on every external surface is computed, on the basis of their own orientation. Buildings 2021, 11, x FOR PEER REVIEW 12 of 27 Figure 9. Coupling between Type 56 and Type 1230 [61].
During the simulations, Type 1230 takes into account:  the solar radiation, the longwave radiation, and the air convection on the external surface of the outside layer;  the energy storage and the conduction in the outside layer;  radiation exchange between the outside layer and the air cavity;  the convective exchanges from all the surfaces facing in the air cavity;  the conduction through the interface layer.
Type 16c has been implemented to model the solar radiation on all the external surfaces. This Type accepts global radiation, ambient temperature, and ambient relative humidity data as input, in order to output several quantities related to the position of the sun, as the diffuse radiation fraction on the horizontal, by estimating the cloudiness of the sky on the basis of the dry bulb temperature and the dew point temperature. Finally, the radiation on every external surface is computed, on the basis of their own orientation.
The effective sky temperature is determined by means of Type 69b, which calculates the long-wave radiation exchanges between the external surfaces and the atmosphere. Type 69b calculates the cloudiness factor as well, on the basis of the dry bulb and the dew point temperatures.
Finally, Type 33e has been implemented in order to calculate the properties of the moist air, in particular, by taking the air temperature, the relative humidity, and the air pressure as input; it returns the density of the air mixture for every timestep, which is then used to calculate the airflow at the OVF inlet. In this way, the inlet airflow is not a fixed value, but it corresponds to the experimental data acquired through the hot wire anemometers, placed as reported in Figure 2a.
In this study, the experimental weather data acquired from 8 December to 31 December have been used as input data for the Type 16c, Type 69b, and Type 33e.
During the simulation, both the time base used to solve the differential equations and the simulation timestep has been set equal to 15 min in order to have a full correlation to the timestep of the experimental input data. The effective sky temperature is determined by means of Type 69b, which calculates the long-wave radiation exchanges between the external surfaces and the atmosphere. Type 69b calculates the cloudiness factor as well, on the basis of the dry bulb and the dew point temperatures.
Finally, Type 33e has been implemented in order to calculate the properties of the moist air, in particular, by taking the air temperature, the relative humidity, and the air pressure as input; it returns the density of the air mixture for every timestep, which is then used to calculate the airflow at the OVF inlet. In this way, the inlet airflow is not a fixed value, but it corresponds to the experimental data acquired through the hot wire anemometers, placed as reported in Figure 2a.
In this study, the experimental weather data acquired from 8 December to 31 December have been used as input data for the Type 16c, Type 69b, and Type 33e.
During the simulation, both the time base used to solve the differential equations and the simulation timestep has been set equal to 15 min in order to have a full correlation to the timestep of the experimental input data.

Validation of the Numerical Model
This sub-section reports the methods and results related to the validation of the numerical model. The model reliability has been verified in terms of indoor air temperature (T indoor ) and the average temperature of the air cavity (T cavity ) by comparing the experimental values with those obtained as an output of the simulation model above described, defining the following percentage differences ∆T indoor and ∆T cavity : ∆T cavity = T cavity,exp − T cavity,sim /T cavity,exp Figure 10a,b report the comparison between the simulation results and the experimental data acquired during the whole test period (from 8 December to 31 December) in terms of T indoor and T cavity , respectively. Figure 10a,b report the comparison between the simulation results and the experimental data acquired during the whole test period (from 8 December to 31 December) in terms of Tindoor and Tcavity, respectively.
These figures highlight how the developed model is quite reliable, with values of ΔTindoor ranging between a minimum of − . 2%, and a maximum of .22%, as well as the values of ΔTcavity between a minimum of − . % and a maximum of . %. The accuracy of the model has also been validated by calculating the mean error (ME), the mean absolute error (MAE), and the root mean square error (RMSE) as reported below: where Texp,i is the experimental value at time step i, Tsim,i is the simulated value at time step i, and N is the number of measurements. Table 2 reports the values of the ME, MAE, and RMSE for both Tindoor and Tcavity. These figures highlight how the developed model is quite reliable, with values of ∆T indoor ranging between a minimum of −9.62%, and a maximum of 6.22%, as well as the values of ∆T cavity between a minimum of −8.51% and a maximum of 7.85%.
The accuracy of the model has also been validated by calculating the mean error (ME), the mean absolute error (MAE), and the root mean square error (RMSE) as reported below: where T exp,i is the experimental value at time step i, T sim,i is the simulated value at time step i, and N is the number of measurements. Table 2 reports the values of the ME, MAE, and RMSE for both T indoor and T cavity . The values reported in Table 2 highlight that there is a slight difference between the measured and predicted results, in particular: (i) the ME associated to the T indoor is equal to −0.3 • C, which means that the simulation model slightly overestimates the indoor air temperature, while that associated to the T cavity is equal to 0.3 • C; (ii) the values of RMSE are equal to 0.4 • C and 0.2 • C for T indoor and T cavity , respectively. Therefore, the results show the ability of the Type 1230 to accurately predict the behavior of the extruded ABS panels in an OVF system. Thus, the same methodology is used to carry out a complete numerical campaign on a set of case studies upon varying the polymer (selecting the ones more used in architecture, as highlighted in the literature review) and manufacturing technology (extruded and 3D printed).

Materials and Numerical Modeling Implementation
The software TRNSYS 18 [24] is used to assess the potential energy saving achievable in an office building refurbishment using plastic and composite polymers as the second-skin layer material.
The office building modeled in this work is the same for all configurations and it consists of three identical floors. Each floor has a surface of 451 m 2 and a volume equal to 1503 m 3 , with a total window area (A w, total ) of 112.3 m 2 (A w, North = 24.5 m 2 , A w, South = 87.8 m 2 ). It is located in Napoli (latitude = 40 • 51 N; longitude = 14 • 16 E), and as such, in order to simulate the weather condition, the corresponding EnergyPlus weather data has been used [30]. The office is firstly modeled in the SketchUp 3D modeling software ( Figure 11). Then, the 3D model geometries were exported by means of the Trnsys3D plugin and successively imported into TRNSYS 18 in order to model the building envelope (stratigraphy of the opaque wall and window typology), to define the infiltration, the internal gains, the operation period of the heating and cooling systems as well as the operation of the electric equipment and lighting system. In particular, the same TRNSYS Types described in Section 4 have been used to simulate the case study. view) and manufacturing technology (extruded and 3D printed).

Materials and Numerical Modeling Implementation
The software TRNSYS 18 [24] is used to assess the potential energy saving achievable in an office building refurbishment using plastic and composite polymers as the secondskin layer material.
The office building modeled in this work is the same for all configurations and it consists of three identical floors. Each floor has a surface of 451 m 2 and a volume equal to 1503 m 3 , with a total window area (Aw, total) of 112.3 m 2 (Aw, North = 24.5 m 2 , Aw, South = 87.8 m 2 ). It is located in Napoli (latitude = 40°51′ N; longitude = 14°16′ E), and as such, in order to simulate the weather condition, the corresponding EnergyPlus weather data has been used [30]. The office is firstly modeled in the SketchUp 3D modeling software ( Figure 11). Then, the 3D model geometries were exported by means of the Trnsys3D plug-in and successively imported into TRNSYS 18 in order to model the building envelope (stratigraphy of the opaque wall and window typology), to define the infiltration, the internal gains, the operation period of the heating and cooling systems as well as the operation of the electric equipment and lighting system. In particular, the same TRNSYS Types described in Section 4 have been used to simulate the case study. The same typical three-story office building is investigated upon, varying the insulation layer thickness on the south facade and the typology of second-skin material, for a total of eight case studies. Table 3 summarizes the eight simulation cases investigated in this work. In particular, this table reports the reference case (Case 0) without the second- The same typical three-story office building is investigated upon, varying the insulation layer thickness on the south facade and the typology of second-skin material, for a total of eight case studies. Table 3 summarizes the eight simulation cases investigated in this work. In particular, this table reports the reference case (Case 0) without the second-skin ( Figure 12a) and seven refurbishment case studies with the OVF system (Figure 12b) upon varying the second-skin material. These cases are: • Case 1, with an OVF system made of a conventional second-skin material (Porcelain gres); • Cases 2-5, where the OVF systems have been implemented by using the extrude plastic and polymer materials more used in architecture (polycarbonate multi-wall sheets, ABS, PETG, and PLA); • Cases 3_3D-5_3D, where the second-skin materials used in the OVF are the most popular 3D printed polymers (ABS, PETG, and PLA); Table 3. Summary of case studies investigated.

Case Study Second-Skin Material Insulation Thickness (m) Air Gap (m)
Case 0 ---  (Figure 12a) and seven refurbishment case studies with the OVF system (Figure 12b) upon varying the second-skin material. These cases are:  Case 1, with an OVF system made of a conventional second-skin material (Porcelain gres);  Cases 2-5, where the OVF systems have been implemented by using the extrude plastic and polymer materials more used in architecture (polycarbonate multi-wall sheets, ABS, PETG, and PLA);  Cases 3_3D-5_3D, where the second-skin materials used in the OVF are the most popular 3D printed polymers (ABS, PETG, and PLA); In addition, an insulation layer has been added in each case study in order to reach the threshold values specified by the Italian Law [63] and equal to 0.36 W/m 2 K for the climatic zone considered in this work. The different insulation thicknesses are also reported in Table 3, upon varying the simulation case. In all the retrofit cases, the cavity inlet airspeed is directly related to the wind speed and direction, as only the wind coming from a similar orientation as the second-skin system (wind direction = 180° ± 45°) has been considered as input for the Type 1230. In addition, the second-skin system has a control logic for the air cavity shutters, which are considered open during the cooling period and closed during the heating period.  In addition, an insulation layer has been added in each case study in order to reach the threshold values specified by the Italian Law [63] and equal to 0.36 W/m 2 K for the climatic zone considered in this work. The different insulation thicknesses are also reported in Table 3, upon varying the simulation case.
In all the retrofit cases, the cavity inlet airspeed is directly related to the wind speed and direction, as only the wind coming from a similar orientation as the second-skin system (wind direction = 180 • ± 45 • ) has been considered as input for the Type 1230. In addition, the second-skin system has a control logic for the air cavity shutters, which are considered open during the cooling period and closed during the heating period. Table 4 shows the thermal-physical properties of the opaque walls of the envelope implemented in the case studies.  Table 5 reports the general simulation parameters adopted in the eight different case studies. In particular, this table highlights: (i) the values of the thermal transmittance for both opaque walls and windows, (ii) the air infiltration rate, (iii) the target of the indoor air temperature, the operation period, and the characteristics of the heating and cooling system, (iv) the occupancy schedule and (v) the internal gains. As can be noticed from Table 5, in all the retrofit cases the values of thermal transmittance for the opaque surfaces are equal to those of the reference case (Case 0), with the exception of the south wall, where the OVFs have been implemented and the thermal transmittance has been set equal to the threshold values specified by Italian law [63], considering the second-skin materials, the air cavity and the insulation thicknesses reported in Table 3.
Two parallel-connected electric heat pump (EHP) devices, model CRA/K 91 [64], coupled with a multi-split type air conditioning system, have been used to cover both the heating and cooling demands. The national grid has been used to cover all the electrical energy demand.
The simulation timestep has been set to 30 min. Finally, Type 1230 [61] has been used to model a second-skin layer in plastic and composite polymers. The parameters required by Type 1230 for each material (i.e., density, thermal capacity, and thermal conductivity), for Cases 1-5, have been derived on the basis of the manufacturers or literature data [40,54,[65][66][67]. With respect to the 3D printed panels (Cases 3_3D-5_3D), several specimens have been printed (Figure 13), in order to measure the final dimensions and density. (considered equal to 12 W/m 2 K [48]) and Vd is the volume fraction of the filler in the polymer matrix (calculated as 0.69 of the specimens' total volume). The internal geometries have been modeled as hexagons as suggested by [50], where hexagons specimens resulted as the most resilient to physical stress, thus more suitable for a building envelope integration.  Figure 13. The specimens made through the 3D printing process: (a) ABS [68]; (b) PETG [69]; (c) PLA [70]. Table 6 reports the parameters used to simulate the second-skin systems in the retrofit cases. The thermal conductivity of the 3D printed materials (k 3D ) has been calculated by means of the equation expressed by [48,51] and reported below: where k m is the thermal conductivity of the polymer as declared by their manufacturers [68][69][70], k d is the thermal conductivity of the filler consisting of the air in the hexagonal cavities of the printed panels (equal to 0.026 W/mK [71]), a is the filler radius measured from the specimens (measured as 0.0045 m), h c is the interfacial boundary conductance (considered equal to 12 W/m 2 K [48]) and V d is the volume fraction of the filler in the polymer matrix (calculated as 0.69 of the specimens' total volume). The internal geometries have been modeled as hexagons as suggested by [50], where hexagons specimens resulted as the most resilient to physical stress, thus more suitable for a building envelope integration.  Table 6 reports the parameters used to simulate the second-skin systems in the retrofit cases. In conclusion, two additional parameters have been set for each case study, required by the Type 1230 as highlighted in Section 4: the resistance of interface layer, to be set in Type 1230 itself, and the convective heat transfer coefficient of the interface layer, to be set in the construction south wall in Type 56 instead. The convective heat transfer coefficients of the back of the south wall value have been set equal to the thermal transmittance value of the insulation layers, which act as interface layers between the two TRNSYS Types, while the thermal resistance values of the interface layer have been simply calculated as the inverse of the convective heat transfer coefficients.
These plastic and polymer materials do not differ only in terms of thermo-physical properties but also in terms of cost. In this work, the capital cost for each retrofit action has been neglected; however, Table 7 provides an overview of the costs per square meter associated with each plastic and polymer material implemented as a second-skin layer in the OVF system [54,[77][78][79]. In general, the polycarbonate multi-wall sheets prove to be the cheaper material (10-25 €/m 2 ), being also the only one which is already used in the building sector; instead, the 3D printed panels are the more expensive ones (188-225 €/m 2 ). This cost limitation is typical for the 3D printing technology, especially compared to more traditional building materials. Despite this significant difference in cost, it must be noted that 3D printing is an emerging technology which usage is still not so widespread. However, the 3D printing technology is the only one that would allow for the obtaining of complex panels' shapes easily. Also, the price of the 3D printed panels reported in Table 7 are related to the brand-new spool, while it is possible to integrate also recycled filament spools in the production process [47]; indeed, the 3D printing manufacturing process is the only one where it's easy to fully integrate eco-compatible materials, like PLA. Table 7. Costs per square meter of plastic and polymer materials used as a second-skin layer in this research [54,[77][78][79].

Energy Analyses: Methods
According to [12,75], the energy comparison between the proposed case (PC) and the reference case (RC) has been carried out considering the non-renewable primary energy consumption through the index PES (non-renewable primary energy saving): where E RC p is the non-renewable primary energy associated with the reference case (Case 0, see Table 3), while E PC p is the non-renewable primary energy associated with the eight proposed cases (Cases 1-5 and Cases 3_3D-5_3D, see Table 3).
The values of the E RC p and E PC p are calculated as reported below: where η PP is the Italian power plants' average efficiency, including the transmission losses, and it is assumed equal to 0.42 [74]. A positive value of the index PES means that the proposed refurbishment allows reducing the non-renewable primary energy consumption with respect to the reference case.

Energy Analyses: Results
In this section, the simulation results of the refurbishment case study are reported and commented on. • all the proposed OVF systems return positive PES values in comparison to the reference case, which means a reduction of the non-renewable primary energy consumption ranging from 2.58% (Case 1) and 2.64% (Cases 2 and 5_3D); this is due to an average reduction of the thermal and cooling energy demands of about 6.9% and 3.0%, respectively; • the retrofit actions where the plastic and composite polymers materials are used as a second-skin layer (Cases 2-5 and Cases 3_3D-5_5D, see Table 3) allow for a slight performance improvement with respect to those realized with a conventional secondskin material (Case 1), thanks to a reduction in the space cooling energy demand ranging from 31 kWh (Case 4) and 120 kWh (Case 5_3D); • the results associated with the polycarbonate multi-wall sheets (Case 2) show a behavior similar to the Cases 3_3D-5_5D, mostly due to the fact that the polycarbonate panels have a structure assimilable to the 3D printing logic; • considering the cases with the same polymers (Case 3 vs. Case 3_3D, Case 4 vs. Case 4_3D, and Case 5 vs. Case 5_3D), the 3D printed panels allow for a slight improvement in performances.
uildings 2021, 11, x FOR PEER REVIEW Figure 14. Values of PES upon varying the case studies.
These figures highlight that:  all the proposed OVF systems return positive PES values in compariso ence case, which means a reduction of the non-renewable primary ene tion ranging from 2.58% (Case 1) and 2.64% (Cases 2 and 5_3D); thi average reduction of the thermal and cooling energy demands of ab 3.0%, respectively;  the retrofit actions where the plastic and composite polymers material second-skin layer (Cases 2-5 and Cases 3_3D-5_5D, see Table 3) allo performance improvement with respect to those realized with a con ond-skin material (Case 1), thanks to a reduction in the space cooling en ranging from 31 kWh (Case 4) and 120 kWh (Case 5_3D);  the results associated with the polycarbonate multi-wall sheets (Case havior similar to the Cases 3_3D-5_5D, mostly due to the fact that the panels have a structure assimilable to the 3D printing logic;  considering the cases with the same polymers (Case 3 vs. Case 3_3D, C 4_3D, and Case 5 vs. Case 5_3D), the 3D printed panels allow for a sl ment in performances.   all the proposed OVF systems return positive PES values in comparison to the refer ence case, which means a reduction of the non-renewable primary energy consump tion ranging from 2.58% (Case 1) and 2.64% (Cases 2 and 5_3D); this is due to an average reduction of the thermal and cooling energy demands of about 6.9% and 3.0%, respectively;  the retrofit actions where the plastic and composite polymers materials are used as a second-skin layer (Cases 2-5 and Cases 3_3D-5_5D, see Table 3) allow for a sligh performance improvement with respect to those realized with a conventional sec ond-skin material (Case 1), thanks to a reduction in the space cooling energy demand ranging from 31 kWh (Case 4) and 120 kWh (Case 5_3D);  the results associated with the polycarbonate multi-wall sheets (Case 2) show a be havior similar to the Cases 3_3D-5_5D, mostly due to the fact that the polycarbonate panels have a structure assimilable to the 3D printing logic;  considering the cases with the same polymers (Case 3 vs. Case 3_3D, Case 4 vs. Case 4_3D, and Case 5 vs. Case 5_3D), the 3D printed panels allow for a slight improve ment in performances.   In order to better investigate the performance of the plastic and polymer materials used as a second-skin layer in the proposed OVF system, the trends of the values of the air temperature inside the cavity and the airspeed at the inlet and the outlet of the air cavity for the Case 5_3D in a typical summer day (2 August) is also reported in Figures 17 and 18 with a 1 h timestep. In general, similar trends in the values of the temperatures in the cavity as well as the airspeed at the cavity inlet and the outlet have been predicted for all the other cases.
In more detail, Figure 17 shows the trends of the daily values of the air temperature inside the cavity upon varying the height from the ground, and, in particular reports, the temperature of the air at the cavity inlet (about 0.15 m from the ground), the average temperature of the air cavity (highlighted by the blue region in the building section) and the temperature of the air at the cavity outlet (about 10.00 m from the ground); Figure 18 reports the values of the airspeed at the cavity inlet (about 0.15 m from the ground) and the airspeed at the cavity outlet (about 10.00 m from the ground), as well as the total vertical radiation on the external surface of the second-skin layer.   In order to better investigate the performance of the plastic and polymer materials used as a second-skin layer in the proposed OVF system, the trends of the values of the air temperature inside the cavity and the airspeed at the inlet and the outlet of the air cavity for the Case 5_3D in a typical summer day (2 August) is also reported in Figures 17 and 18 with a 1 h timestep. In general, similar trends in the values of the temperatures in the cavity as well as the airspeed at the cavity inlet and the outlet have been predicted for all the other cases. In order to better investigate the performance of the plastic and polymer materials used as a second-skin layer in the proposed OVF system, the trends of the values of the air temperature inside the cavity and the airspeed at the inlet and the outlet of the air cavity for the Case 5_3D in a typical summer day (2 August) is also reported in Figures 17 and 18 with a 1 h timestep. In general, similar trends in the values of the temperatures in the cavity as well as the airspeed at the cavity inlet and the outlet have been predicted for all the other cases.
In more detail, Figure 17 shows the trends of the daily values of the air temperature inside the cavity upon varying the height from the ground, and, in particular reports, the temperature of the air at the cavity inlet (about 0.15 m from the ground), the average temperature of the air cavity (highlighted by the blue region in the building section) and the temperature of the air at the cavity outlet (about 10.00 m from the ground); Figure 18 reports the values of the airspeed at the cavity inlet (about 0.15 m from the ground) and the airspeed at the cavity outlet (about 10.00 m from the ground), as well as the total vertical radiation on the external surface of the second-skin layer.   The data reported in Figure 17 corresponds to the input and outputs of Type 1230. In particular, the inlet air temperature is an input for Type 1230, while the average air temperature in the whole cavity (Tcavity, the blue-edged markers in Figure 17) and the outlet air temperature are returned as output results by Type 1230 itself. These values represent the only two temperature values associated with the air cavity returned by Type 1230 [61]. As a first observation, the overall temperature distribution is directly related to solar radiation throughout the day, where the temperatures rise during the morning and drop during the afternoon. Also, during the day, the temperature trend seems to be constantly rising from the air cavity inlet to the air cavity outlet; this behavior is due to the chimney effect that is created thanks to the OVF system. Figure 18 reports the simulation data corresponding to the cavity inlet airspeed (input of the Type 1230) and the cavity outlet airspeed (output of the Type 1230), as well as the total vertical solar radiation. This figure shows how, during the day, the values predicted at the cavity air outlet are always higher than those at the cavity air inlet (with a difference between outlet and inlet ranging from 0.01 m/s and 1.46 m/s), gradually rising to a maximum peak at around 14:00.
Finally, in order to verify the potential benefits coming from the best case, an additional simulation case, not reported in Table 3, has been carried out. In this last simulation, the OVF system has been implemented on the whole building, following the same installation methodology of the previous cases. The material selected as a second-skin layer is the 3D printed PLA, which proved to be one of the most effective in improving the nonrenewable primary energy saving. The proposed OVF system returned a PES value equal to 8.10% if compared to the reference case.

Conclusions
The OVFs have been, more and more frequently, chosen for different building typologies (offices, schools, residential) and in different climates. Nowadays, there are always more innovative materials used in architecture and as a second-skin layer, even if the evaluation of their impact on the envelope's energy performance is a complex task. In particular, the use of polymers in building and engineering has increased substantially, thanks to their: (i) ease of production, (ii) ease of installation, (iii) durability, (iv) low maintenance 10 In more detail, Figure 17 shows the trends of the daily values of the air temperature inside the cavity upon varying the height from the ground, and, in particular reports, the temperature of the air at the cavity inlet (about 0.15 m from the ground), the average temperature of the air cavity (highlighted by the blue region in the building section) and the temperature of the air at the cavity outlet (about 10.00 m from the ground); Figure 18 reports the values of the airspeed at the cavity inlet (about 0.15 m from the ground) and the airspeed at the cavity outlet (about 10.00 m from the ground), as well as the total vertical radiation on the external surface of the second-skin layer.
The data reported in Figure 17 corresponds to the input and outputs of Type 1230. In particular, the inlet air temperature is an input for Type 1230, while the average air temperature in the whole cavity (T cavity , the blue-edged markers in Figure 17) and the outlet air temperature are returned as output results by Type 1230 itself. These values represent the only two temperature values associated with the air cavity returned by Type 1230 [61]. As a first observation, the overall temperature distribution is directly related to solar radiation throughout the day, where the temperatures rise during the morning and drop during the afternoon. Also, during the day, the temperature trend seems to be constantly rising from the air cavity inlet to the air cavity outlet; this behavior is due to the chimney effect that is created thanks to the OVF system. Figure 18 reports the simulation data corresponding to the cavity inlet airspeed (input of the Type 1230) and the cavity outlet airspeed (output of the Type 1230), as well as the total vertical solar radiation. This figure shows how, during the day, the values predicted at the cavity air outlet are always higher than those at the cavity air inlet (with a difference between outlet and inlet ranging from 0.01 m/s and 1.46 m/s), gradually rising to a maximum peak at around 14:00.
Finally, in order to verify the potential benefits coming from the best case, an additional simulation case, not reported in Table 3, has been carried out. In this last simulation, the OVF system has been implemented on the whole building, following the same installation methodology of the previous cases. The material selected as a second-skin layer is the 3D printed PLA, which proved to be one of the most effective in improving the non-renewable primary energy saving. The proposed OVF system returned a PES value equal to 8.10% if compared to the reference case.

Conclusions
The OVFs have been, more and more frequently, chosen for different building typologies (offices, schools, residential) and in different climates. Nowadays, there are always more innovative materials used in architecture and as a second-skin layer, even if the evaluation of their impact on the envelope's energy performance is a complex task. In particular, the use of polymers in building and engineering has increased substantially, thanks to their: (i) ease of production, (ii) ease of installation, (iii) durability, (iv) low maintenance requirements, (v) lightweight nature and (vi) ability to be formed into complex shapes. Several of these plastic products can also be utilized in additive manufacturing processes, providing excellent freedom of form, enhancing designers, architects, and engineers' freedom in creating complex designs.
In this work, the numerical model of extruded ABS panels in an OVF system has been developed and validated. Then, the simulation methods, suggested by the authors, have been applied in different refurbishment cases upon varying: the polymer and the manufacturing technology, extrusion (polycarbonate multi-wall sheets, ABS, PETG, and PLA), and 3D printing (ABS, PETG, and PLA).
The simulations have been carried out in order to assess the potential benefits achievable in terms of non-renewable primary energy saving, as well as thermal and cooling energy demand reduction. The simulation results highlight that: (i) all the proposed retrofit cases allow to achieve a benefit in terms of PES; (ii) the plastic and composite polymers materials allow for a slight performance improvement with respect to conventional second-skin material, such as porcelain gres; (iii) the best performances among the extruded polymers are returned by the polycarbonate multi-wall sheets (PES value equal to 2.64%); (iv) the best performances among the 3D printed polymers are achieved when using the PLA (PES value equal to 2.64%). However, the polymers' results show very similar performances, thus allowing building contractors, designers, and architects to select the material based on other project requirements, as mechanical strength, weather resistance, environmental impact, etc.
In this work, the thermal conductivity of the 3D printed materials has been calculated by means of an equation expressed in literature; therefore, in order to improve the accuracy of the model for the 3D printed material, in future works, the authors will carry out experimental investigations on full-scale 3D printed panels in OVF systems through the Gemini test cells. In addition, the capital cost for each retrofit action has been neglected. However, they represent an important parameter in the refurbishment typology choice; therefore, in future work, the authors will focus on a detailed economic analysis considering both the operating cost reduction and the simple payback period.

Data Availability Statement:
The data presented in this study are available on request from the corresponding authors.