Performance Analysis of a Hybrid Thin Film Photovoltaic (PV) Vacuum Glazing

Hasila Jarimi; Ke Qu; Shihao Zhang; Qinghua Lv; Jun Liao; Benyuan Chen; Hui Lv; Chunfu Cheng; Jin Li; Yuehong Su; Shijie Dong; Saffa Riffat

1. Introduction

In recent years, researchers have shown interest in the study of vacuum glazing technology. Innovative designs were developed and researched. Among them are; Bao, Liu () who have introduced a novel hybrid vacuum/triple glazing system in which a vacuum glazing unit is enclosed by two glass panels to form a triple glazing unit system. To eliminate additional loads due to pressure difference which can cause unwanted breakage, they designed an equalized air pressure arrangement on both sides of the vacuum glazing unit. Fang, Hyde () have performed simulation model using a finite volume model on thermal performance of triple vacuum glazing (TVG) sealed using indium alloy, and 4 mm thick glass panes. In their validated computer simulation work, the emissivity values of the glass sheets were varied from 0.18 to 0.03. The simulated centre U-values are reduced from 0.41 Wm^–2K^–1 to 0.22 Wm^–2K^–1 for a 0.4 m by 0.4 m size triple vacuum glazing. Memon () introduced a triple vacuum glazing using Cerasolzer CS186 alloy as the edge sealant. From the 3D finite element model, the centre of glass U-value of 0.33 Wm^–2K^–1 is found achievable. Memon and Eames () investigated the building heating performance when retrofitted with the composite edge sealed triple vacuum glazing. An annual space-heating energy savings of 14.58% and 15.31% are expected to be achieved. Meanwhile, if single glazing windows are to be replaced with the designed triple vacuum glazing, a heat loss reduction of 12.92% and 2.69% are achievable. Research on vacuum glazing technology has also led to integrating the technology with other advanced glazing technology. Fang, Hyde () have introduced an integration between electrochromic (EC) with vacuum glazing (VG) technology. The EC component was arranged facing the outdoor environment. For an incident solar irradiance of 300 W m^–2, simulations show that the EC layer is opaque and the inside temperature of the glass pane is higher than the indoor temperature. Meanwhile, for solar irradiance of 1000 W m^–2, the outdoor glass pane temperature exceeds the indoor glass pane temperature, of which as a consequence heat is transferred from outdoors to indoors. Ghosh, Norton () integrated suspended particle device (SPD) with vacuum glazing. In their study, a strong correlation between the measures of the clearness of the atmosphere or clearness index was evaluated for south facing vertical glazing. For clearness index below 0.5, the SPD-vacuum glazing transmission was 17% and 1.1% for transparent and opaque states respectively. A group of researchers from China Zhang, Lu and Chen () proposed an integration between a thin film PV glazing with a vacuum glazing. The glazing prototype was designed and manufactured in such a way that a lamination layer was sandwiched between the thin film PV glazing unit and a vacuum glazing unit, both are independently manufactured. According to the researchers, although semi-transparent photovoltaic windows can generate electricity in situ, they also increase the cooling load of buildings significantly due to the waste heat as a by-product. Experiment has been conducted and the results have indicated that the prototype can not only generate electricity, but also help reduce the cooling load as well as improve the indoor thermal comfort. Recently, researchers Ghosh, Sundaram () investigated a combined semi-transparent multi-crystalline PV vacuum glazing. In their research, both thermal and electrical performance have been investigated. They have found that the combination of multi-crystalline PV cells with vacuum glazing provide low overall heat transfer coefficient, reduces solar heat gain, generates clean electricity and allows comfortable daylight. This paper introduces a design of thin film PV glazing combined with a vacuum glazing. The novelty of the glazing lies in the designs, manufacturing and application of the PV VG-4L. The heat transfer analysis based on the analytical analysis is discussed and solved numerically in MATLAB. Experimental works are presented and the model was then validated against indoor and outdoor experimental results.

2. Design concept and heat transfer analysis of thin film PV vacuum insulated glazing (PV VG-4L)

The design of PV VG-4L is illustrated in Figure 1a. PV VG-4L involves an integration between a thin film PV glazing with a double vacuum glazing (both manufactured independently) and an additional layer of self-cleaning coated glass which totalling 4 layers of glass units which gives the total thickness of 14 mm. The layers were combined together via Ethylene Vinyl Acetate (EVA) film using an autoclave facility set at its optimum temperature and pressure.

Figure 1

a) Schematics of PV VG-4L and b) the thermal resistance network.

The thermal resistance network for the PV VG-4L is shown in Figure 1b. Both the width and the length of the PV VG-4L have the same length and denoted as W. The heat transfer was analysed based on the following temperature nodes which are; the average glass temperature facing indoors denoted by Tg₁, the average internal glass temperature denoted by Tg₂, the average temperature of the thin film PV glazing denoted by T_TF and the average temperature of the self-cleaning coated glass denoted by Tg₄. The lateral heat transfer ends at the point right at the edge sealing and the temperature was denoted as T_edge₁. Due to high thermal conductivity of the edge sealing, it acts as a short circuit for the heat transfer at the edges of the glazing. Meanwhile, due to high thermal conductivity between g2 and g4, the lateral heat conduction in the glass slab g2 and thin film PV glazing (TF) are assumed negligible. The temperature of the position in g4, of which the edge heat transfer takes place is denoted as T_edge₂. For the uninsulated edges, there will be heat transferred from the internal ambient to the area around the edges.

The thermal resistance network implies that, at certain level of solar irradiance and ambient temperature, the internal surface temperature is in general lower than the external surface temperature. This can be explained as follows. When the PV panel absorbed solar irradiance, a fraction of the absorbed solar irradiance will be converted into electricity meanwhile a fraction will be wasted in the form of heat. This heat released by the PV panel is insulated from the building via the vacuum layer of which, the main heat transfer occurs in the gap are mainly due to radiation and conduction through the support pillars. In the summer, this can be considered as an advantage since the installation of PV glazing would normally cause additional heating during summer of which in return, increase the cooling load of the building. Based on the thermal nodal networks, the energy balance equations were developed. In order to simplify the analysis, the following assumptions have been made:

The heat transfer involved is assumed symmetrical. Therefore, we only consider a quarter of the PV VG-4L area in the heat transfer analysis.
In this study we have considered the lateral heat transfer through the glass slab in g1. This is due to the fact that, the vacuum glazing is highly thermally insulated which makes the lateral heat transfer to become prominent.
The edge boundaries of the PV VG-4L are well insulated and hence the edge losses are assumed negligible.
Due to high thermal conductivity of the edge sealing, the edge seal is assumed as a thermal short circuit.
Meanwhile, due to high thermal conductivity between g2 and g4, the lateral heat conduction in the glass slab g2 and thin film PV glazing (TF) are assumed negligible.

To simulate the performance of the PV VG-4L, the following energy balance equations were developed for each of the temperature nodes.

For g1:

(1)

S g 1 ︸ 1 + h a_g 1 T a_i − T g 1 ︸ 2 = h g 1_g 2 T g 1 − T g 2 ︸ 3 + h k_edge T g 1 − T edge 1 ︸ 4

M1 \documentclass[10pt]{article} \usepackage{wasysym} \usepackage[substack]{amsmath} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage[mathscr]{eucal} \usepackage{mathrsfs} \usepackage{pmc} \usepackage[Euler]{upgreek} \pagestyle{empty} \oddsidemargin -1.0in \begin{document} \[ \begin{array}{l} \underbrace {{S_{g1}}}_1\,\, + \,\,\underbrace {{h_{a\_g1}}\left({{T_{a\_i}} - {T_{g1}}} \right)}_2\,\, = \,\,\underbrace {{h_{g1\_g2}}\left({{T_{g1}} - {T_{g2}}} \right)}_3\;\\ \qquad\qquad\qquad\qquad\qquad\qquad + \,\,\underbrace {{h_{k\_edge}}\left({{T_{g1}} - {T_{edge1}}} \right)}_4 \end{array} \] \end{document}

The heat transfer terms are defined as follows:

1: The rate of the solar energy received by the glass cover or surface facing indoors after transmission through different glass layers per unit area; 2: The rate of heat transfer from the indoor ambient to g1 per unit area; 3: The rate of heat transfer from g1 to g2 per unit area. The heat transfer includes heat conduction through the glass slabs of g1, which then followed by heat transfer in the vacuum gap which are due to the heat conduction of the gas particles, heat transfer via radiation and heat conduction through the support pillars, and then followed by heat conduction through the glass slab g4; 4: The rate of lateral heat conducted along the x- direction and y-direction of the PV VG-4L from the centre of the PV VG-4L.

For edge1:

(2)

h k_edge T g 1 − T edge 1 ︸ 4 + h iedge γ edge T a_i − T edge 1 ︸ 5 = h edge γ edge T edge 1 − T edge 2 ︸ 6

M2 \documentclass[10pt]{article} \usepackage{wasysym} \usepackage[substack]{amsmath} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage[mathscr]{eucal} \usepackage{mathrsfs} \usepackage{pmc} \usepackage[Euler]{upgreek} \pagestyle{empty} \oddsidemargin -1.0in \begin{document} \[\begin{array}{l} \underbrace {{h_{k\_edge}}\left({{T_{g1}} - {T_{edge1}}} \right)}_4\,\, + \,\,\underbrace {{h_{iedge}}{\gamma _{edge}}\left({{T_{a\_i}} - {T_{edge1}}} \right)}_5\\ \,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\underbrace { = \,\,{h_{edge}}{\gamma _{edge}}\left({{T_{edge1}} - {T_{edge2}}} \right)}_6 \end{array}\] \end{document}

The heat transfer terms are defined as follows:

5: The rate of heat transfer from the indoor ambient to the edge area per unit glazing area (if not insulated); 6: The rate of heat transfer from edge 1 to edge 2 through g2, TF and EVA layer per unit area.

For g2:

(3)

h g 1_g 2 T g 1 − T g 2 ︸ 3 + S g 2 ︸ 7 = h g 2_TF T g 2 − T TF ︸ 8

M3 \documentclass[10pt]{article} \usepackage{wasysym} \usepackage[substack]{amsmath} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage[mathscr]{eucal} \usepackage{mathrsfs} \usepackage{pmc} \usepackage[Euler]{upgreek} \pagestyle{empty} \oddsidemargin -1.0in \begin{document} \[\underbrace {{h_{g1\_g2}}\left({{T_{g1}} - {T_{g2}}} \right)}_3\,\, + \,\,\underbrace {{S_{g2}}}_7 = \,\,\underbrace {{h_{g2\_TF}}\left({{T_{g2}} - {T_{TF}}} \right)}_8\] \end{document}

The heat transfer terms are defined as follows:

7: The rate of the solar energy received by g2 after transmission through different glass layers per unit area; 8: The rate of heat transfer from g2 to TF through the EVA layer per unit area.

For TF:

(4)

h g 2_TF T g 2 − T TF ︸ 8 + S TF ︸ 9 = h TF_g 4 T TF − T g 4 ︸ 10

M4 \documentclass[10pt]{article} \usepackage{wasysym} \usepackage[substack]{amsmath} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage[mathscr]{eucal} \usepackage{mathrsfs} \usepackage{pmc} \usepackage[Euler]{upgreek} \pagestyle{empty} \oddsidemargin -1.0in \begin{document} \[\underbrace {{h_{g2\_TF}}\left({{T_{g2}} - {T_{TF}}} \right)}_8\;\, + \underbrace {{S_{TF}}}_9 = \underbrace {{h_{TF\_g4}}\left({{T_{TF}} - {T_{g4}}} \right)}_{10}\] \end{document}

The heat transfer terms are defined as follows:

9: The rate of the solar energy received by TF after transmission through different glass layers per unit area; 10: The rate of heat transfer from TF to g4 through the EVA layer per unit area.

For Tedge2:

(5)

h edge γ edge T edge 1 − T edge 2 ︸ 6 = h k_edge T edge 2 − T g 4 ︸ 11

M5 \documentclass[10pt]{article} \usepackage{wasysym} \usepackage[substack]{amsmath} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage[mathscr]{eucal} \usepackage{mathrsfs} \usepackage{pmc} \usepackage[Euler]{upgreek} \pagestyle{empty} \oddsidemargin -1.0in \begin{document} \[\underbrace {{h_{edge}}{\gamma _{edge}}\left({{T_{edge1}} - {T_{edge2}}} \right)}_6\,\,\, = \,\,\,\underbrace {{h_{k\_edge}}\left({{T_{edge2}} - {T_{g4}}} \right)}_{11}\] \end{document}

The heat transfer terms are defined as follows:

11: The rate of lateral heat conducted along the x-direction and y-direction of the PV VG-4L from the centre of the PV VG-4L.

For g4:

(6)

h TF_g 4 T TF − T g 4 ︸ 10 + h k_edge T edge 2 − T edge 4 ︸ 11 + S g 4 ︸ 12 = h o T g 4 − T a_o ︸ 14

M6 \documentclass[10pt]{article} \usepackage{wasysym} \usepackage[substack]{amsmath} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage[mathscr]{eucal} \usepackage{mathrsfs} \usepackage{pmc} \usepackage[Euler]{upgreek} \pagestyle{empty} \oddsidemargin -1.0in \begin{document} \[ \begin{array}{l} \underbrace {{h_{TF\_g4}}\left({{T_{TF}} - {T_{g4}}} \right)}_{10} + \,\,\underbrace {{h_{k\_edge}}({{T_{edge2}} - {T_{edge4}}} )}_{11} \underbrace {\, + {S_{g4}}}_{1}\\ \qquad\qquad\qquad\qquad\ = \,\,\underbrace {{h_o}\left({{T_{g4}} - {T_{a\_o}}} \right)}_{14} \end{array} \] \end{document}

The heat transfer terms are defined as follows:

12: The rate of the solar energy received by g4 after transmission through different glass layers per unit area; 14: The rate of heat transfer from g4 to T_a_{_}_o per unit area

The thin film PV VG-4L introduced in this study is new and has never been discussed in existing research. Therefore, in order to theoretically predict the glazing’s electrical performance, the equation that has been widely used in the study of PV/Thermal solar collector is utilized. The following correlation developed by Schott () and Evans () is the most common correlation implemented by researches who study the PV/T –type solar collector such as Sarhaddi, Farahat () and Tonui and Tripanagnostopoulos ():

(7)

η eff = η Tref 1 − β ref T TF − T ref

M7 \documentclass[10pt]{article} \usepackage{wasysym} \usepackage[substack]{amsmath} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage[mathscr]{eucal} \usepackage{mathrsfs} \usepackage{pmc} \usepackage[Euler]{upgreek} \pagestyle{empty} \oddsidemargin -1.0in \begin{document} \[{\eta _{eff{\rm{\;}}}} = \,\,\,{\eta _{Tref}}\left({1 - {\beta _{ref}}\left({{T_{TF}} - \,\,{T_{ref}}} \right)} \right)\] \end{document}

Where β_ref is the temperature coefficient of the thin film PV cells and η_Tref is the electrical efficiency at reference temperature (25°C). In this study, by referring to (), for the a-Si thin film PV layer, the value of β_ref is –0.36%, meanwhile, the η_Tref is 4% by referring to the preliminary experimental results conducted by the current authors. It is important to note that, under real solar condition, the amount of the irradiance reaching the thin-film PV glass component will be influenced by the optical properties of the glazing which depends on the solar angle incidence that varies with the time of the day, and also the day of the year.

From the simulation, the values of the average temperature of the glass sheets (layers) were used to compute the heat transfer coefficients summarised in Table 1. By referring to Figure 1b the centre thermal resistance R_centre and the edge thermal resistance R_edge can be computed as in equation (8) and (9) respectively.

(8)

R centre = 1 h a_g 1 + 1 h g 1 g 2 + 1 h g 2 TF + 1 h T F g 4 + 1 h o

M8 \documentclass[10pt]{article} \usepackage{wasysym} \usepackage[substack]{amsmath} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage[mathscr]{eucal} \usepackage{mathrsfs} \usepackage{pmc} \usepackage[Euler]{upgreek} \pagestyle{empty} \oddsidemargin -1.0in \begin{document} \[{R_{centre}} = \,\,\frac{1}{{{h_{a\_g1}}}} + \frac{1}{{{h_{g{1_{g2}}}}}} + \frac{1}{{{h_{g{2_{TF}}}}}} + \frac{1}{{{h_{T{F_{g4}}}}}} + \frac{1}{{{h_o}}}\] \end{document}

(9)

R edge_tot = 1 h edge + 1 h k_edge + 1 h iedge

M9 \documentclass[10pt]{article} \usepackage{wasysym} \usepackage[substack]{amsmath} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage[mathscr]{eucal} \usepackage{mathrsfs} \usepackage{pmc} \usepackage[Euler]{upgreek} \pagestyle{empty} \oddsidemargin -1.0in \begin{document} \[{R_{edge\_tot}} = \,\,\,\frac{1}{{{h_{edge}}}} + \frac{1}{{{h_{k\_edge}}}} + \frac{1}{{{h_{iedge}}}}\] \end{document}

Table 1

The definition of the symbol used.

Symbol	Definition	Symbol	Definition

α	Solar absorptance	τ	Solar transmittance
γ_edge	Correction factor due to the total area of edge sealing to the centre area	G_rad	Total global solar irradiance incident upon the glazing structure.
δ_edge	The edge sealing thickness	h_edge	Conductive heat transfer coefficient due to the edge sealing
h_P	Heat transfer coefficient through the support pillars	h_iedge	Total heat transfer coefficient from the indoor ambient air to the uninsulated edge sealing area
h_{k_edge}	The lateral heat transfer to the glass edge	h_{g1_g2}	The total heat transfer coefficient from g1 to g2 of vacuum glazing
h_{a_g1}	Total heat transfer coefficient from the internal ambient to g1	h_{g2_TF}	The total coefficient of heat transfer by conduction from g2 to the thin film PV glazing layer TF.
h_{TF_g4}	The total heat transfer coefficient by conduction from the thin film PV glazing layer TF to g4.	k_edge	Thermal conductivity through the edge sealing
h_o	Total heat transfer coefficient from the external glass g4 to the ambient	T_{a_i}	Indoor ambient temperature
T_{a_o}	Outdoor ambient temperature

Hence, the centre U-value (U_centre) and the total U-value (U_total) may be computed using equation 10 and 11 respectively:

(10)

U centre = 1 R centre

M10 \documentclass[10pt]{article} \usepackage{wasysym} \usepackage[substack]{amsmath} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage[mathscr]{eucal} \usepackage{mathrsfs} \usepackage{pmc} \usepackage[Euler]{upgreek} \pagestyle{empty} \oddsidemargin -1.0in \begin{document} \[{U_{centre}} = \,\,\,\frac{1}{{{R_{centre}}}}\] \end{document}

(11)

U total = γ c U centre + γ edge U edge

M11 \documentclass[10pt]{article} \usepackage{wasysym} \usepackage[substack]{amsmath} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage[mathscr]{eucal} \usepackage{mathrsfs} \usepackage{pmc} \usepackage[Euler]{upgreek} \pagestyle{empty} \oddsidemargin -1.0in \begin{document} \[{U_{total}} = \,\,{\gamma _c}{U_{centre\, + \,}}{\gamma _{edge}}{U_{edge}}\] \end{document}

In this study, to solve the energy balance equations, we have used the inverse matrix method. MATLAB is used to carry out the iteration process. Newton-Raphson iteration technique is used to estimate the temperature and hence the temperature-dependant heat transfer coefficients of the variables.

3. Experimental validation of the theoretical model

3.1. U-value measurement using a TEC driven calibrated hot box

A PV VG-4L prototype using an amorphous silicon (α-Si) solar cell as shown in Figure 2 was manufactured. The U-value of the prototype was evaluated using the TEC-driven calibrated hot box built at the University of Nottingham. Interested readers may refer to Jarimi () for further details. As can be seen in Figure 3, by following closely ISO 12567 standards, the sample was installed at the specimen area of the calibrated hot box. It was tested under three different air temperature conditions summarised in Table 2. However, the air speed in the hot and cold side were fixed at 0.3 m/s and 1.5 m/s respectively. Using the calibrated hot box, we could estimate the total heat transfer coefficient from the hot and cold surface of the PV VG-4L prototype. The values were then used as the input parameters for the computer simulation. To derive the absolute error, the Kline–McClintock second power law as given in NCEES (National Council of Examiners for Engineering and Surveying) (2001) is used. These errors were represented by the error bars of the associated curves. Additionally, the guideline in ISO 12567 was also being referred to evaluate the error from indoor testing.

Figure 2

The PV VG-4L prototype.

Figure 3

The calibrated hot box with the installed sample.

Table 2

Comparison between the Theoretical and Experimental Results for the glazing surface temperature and U-value using the calibrated hot box under different sets of controlled ambient temperature conditions.

Cond.	Text (°C)	Tint (°C)	Tint (Exp) (°C)	Text (Exp) (°C)	Tint (Theo) (°C)	Text (Theo) (°C)	U-value total (Exp)	U-value total (Theo)

1	12.7	32.70	30.34	13.50	29.90	13.37	0.56	0.57
2	17.5	27.5	27.40	18.65	26.42	17.84	0.57	0.57
3	7.6	27.83	25.19	9.13	25.67	8.29	0.56	0.56

The mathematical model validation method is performed by comparing the results obtained experimentally and theoretically based on the trends shown on the related graphs. In this study, the mathematical model has been validated against the indoor experimental data with the input parameters recorded in the experiment were used in the computer simulation for all the three different conditions. In addition to the direct comparison between the simulation and theoretical curves, the validation of the mathematical model is further justified using root mean square percentage deviation (RMSPD). As shown in Figure 4, and summarised in Table 3, the evaluated glazing surface temperatures and U-value are found to be in good agreement such that the trend of the theoretical curves are consistent with the experimental curves and the computed RMSPD for the temperatures and U-value are 4.02% and 0.92% respectively.

Figure 4

Comparison between the theoretical and experimental results for the glazing surface temperature and U-value using the calibrated hot box.

Table 3

Measured parameters and uncertainties.

Parameters	Sensor	Uncertainties

Temperatures	K-type Thermocouples	±1.5°C
Heat flux	Hukesflux Thermal Sensors	±1.9 × 10^–6 V/(W/m²)
Solar irradiance	Pyranometer	±5%
Maximum power produced by the PV VG-4L	RO4 with Keysight 34972A	Electric current (I) (±1.5 μA) Voltage (V) (±190 μV)

3.2. Performance analysis under real conditions

The developed mathematical model has been validated against indoor experimental analysis. Nonetheless, the true performance of the PV VG-4L under real sky conditions still needs to be investigated in order to further justify the validity of the mathematical model especially that the electrical performance of the thin film PV glazing could not be evaluated indoors. That said, this section first discusses the performance of the PV VG-4L under real conditions. To carry out the testing, the prototype was installed at E.ON 2016 research house at the University of Nottingham, United Kingdom with latitude of 52.9438°N, and longitude of 1.1934°W. It is worth emphasising that the outdoor monitoring of the PV VG-4L has been conducted under two conditions; during the day and during the night (at zero solar irradiance). The thermal and electrical characteristics of the PV VG-4L under real conditions were monitored using the sensors as summarised in Table 3. During the day, the reading given by the heat flux sensor for the thermal transmittance of the PV VG-4L is not reliable due to the error from the solar irradiance to the heat flux sensor and also due to the absorbed solar irradiance by the PV layer. Therefore, during the day, the only parameters that are being considered are the surface temperature difference of the PV VG-4L and its electrical power produced. The typical monitoring period is between 9:30 a.m. to 6:00 p.m. in a typical day of May and June. Meanwhile, the focus of the experiment during the night is the thermal transmittance measurement or U-value. A 1kW radiator was used as the heat source for indoors.

Figure 5 shows the solar irradiance, surface temperatures of the PV VG-4L, and both internal and external; glazing surface temperature during the day and night. It should be noted that for the data taken on the 23^rd to 24^th of May 2019, the internal heat source was switched off during the day meanwhile, the average heating temperature was set at 30°C during the night (i.e. from 8 p.m. to 4 a.m. the next day) to heat the ambient room at 24 to 25°C. On the day of testing, the sky was in a clear sky condition and hence clear pattern of solar irradiance curve with the time of the day was obtained. Another set of experiment was conducted to evaluate the performance of the PV VG-4L on the 20^th to 21^st of June 2019. However, as can be seen on Figure 6, during the day of testing, the sky turned cloudy and hence periods of fluctuating solar irradiance were recorded. During the aforementioned period, steady state condition of the PV VG cannot be well justified. This is attributed to the response time of the glass (time constant) which is caused by the thermal mass of the PV VG components. Furthermore, the influence of the thermal time response of the PV VG to the sudden drop in solar irradiance is clearly shown by the change in the recorded temperature profiles from 1:45 p.m. to 2:15 p.m. During the aforementioned period, the average intensity of the incident solar irradiance dropped for a step change from 400 W/m² to 60 W/m² and the temperature profiles clearly follows an exponential decay lagging behind the step change. Additionally, by carefully examining Figure 6, at points in which there was a sudden drop in the value of solar irradiance due to the movement of clouds, the temperature profiles are also observed to lag behind the values of solar irradiance.

Figure 5

Variation in the external glazing surface temperature (Texternal), internal glazing surface temperature (Tinternal), Outdoor ambient temperature (To_ambient), indoor ambient temperature (Ti_ambient) and solar irradiance with time of the day (23^rd of May (7:00 a.m) to 24^th of May 2019 (4:00 a.m.)).

Figure 6

Variation in the external glazing surface temperature (Texternal), internal glazing surface temperature (Tinternal), Outdoor ambient temperature (To_ambient), indoor ambient temperature (Ti_ambient) and solar irradiance with time of the day (20^th of June from 9:15 a.m. to 21^st of June 2019 9:15 a.m.).

It is worth noting that, in both graphs, at low–zero solar irradiance, the indoor surface temperature of the PV VG-4L is in general higher than its external surface temperature. However, as the incident solar irradiance increases, the external surface temperature increases. The trend of the graph is explained as follows; when the PV component of the PV VG-4L absorbed the incident solar irradiance, a fraction was converted into electrical energy meanwhile the rest was wasted in the form of heat. However, the vacuum layer behind the PV component of the PV VG-4L act as the insulation layer or barrier to the wasted heat from being transferred indoors. As a result, the external surface temperature of the PV VG-4L became higher compared to its internal surface temperature. In the summer, this will be an advantage in comparison to the typical installation of BIPV in double glazed configuration. For the data analysis, the temperatures, solar irradiance and heat flux were recorded for every 1 s meanwhile, the electrical parameters were recorded for every 10 s.

The outdoor experimental results obtained from the 23^rd to 24^th of May, as discussed previously, were compared with the theoretical results using the developed mathematical model. Due to the varying condition of the ambient climate, the electrical performance and thermal characteristics of PV VG-4L in a steady state condition are analysed as per time constant of the PV VG-4L. From our analysis, it is concluded that the computed time constant for the PV VG-4L is approximately 30 minutes. The experimental results were compared based on the PV VG-4L surface temperature difference and electrical performance at quasi-steady state during the day and thermal transmittance or U-value during the night. The comparisons between the outdoor experimental and theoretical results are represented in Figures 7 and 8 for the reading during the day and night respectively. The trend given by both outdoor experimental and theoretical curves are in good agreement. Figure 7 shows that the maximum power produced at average solar irradiance of 600 W/m² is approximately 14W/m² for 0.4 m × 0.4 m PV VG-4L made of amorphous silicon solar cells. It is worth emphasizing that, a different power produced is expected when a different type of thin film PV is used as the prototype. For example, using the validated mathematical model, it is predicted that, if amorphous/microcrystalline silicon solar cells at 20% transparency is used as the thin film PV layer, the power output at the same average solar irradiance can achieve as high as 32 W/m². In order to further justify the validity of the mathematical model, error analysis using RMSPD analysis was performed. The average RMSPD for the glazing surface temperatures, the power produced and U-value are 2.90%, 8.7% and 2.89% respectively. It is worth noting that the derived absolute errors for the power produced are too small to be included in the plotted curves.

Figure 7

Top: Outdoor Experimental and Theoretical Curves for the Glazing surface temperatures (Tint for the internal surface and Text for the external surface) and bottom: Power produced per m² for 0.4 m × 0.4 m prototype as per time constant of the PV VG-4L with time of the day as (data was taken on the 23^rd of May 2019 during the day).

Figure 8

Top: Outdoor Experimental and Theoretical Curves for the Glazing surface temperature and bottom: U-value for 0.4 m × 0.4 m prototype as per time constant of the PV VG-4L with time of the day as (data was taken on the 23^rd of May 2019 to 24^th of May 2019 at low-zero solar irradiance).

4. Parametric studies of PV VG-4L

Using the validated mathematical model, the U-value of the PV VG-4L was evaluated at different pairs of width and length as summarised in Table 4. The U-value improves (decreases with the increase in the glazing size) due to the lateral heat conduction across the glass due to the edge seal. The finding is in agreement with the findings in Fang ().

Table 4

The parameters and the computed U-values.

Lw (Glazing width) (m)	Lw (Glazing length) (m)	U-value (W/m²K)

0.3	0.3	0.66
0.5	0.5	0.55
1.0	1.0	0.51
1.5	1.5	0.50

At 1.0 m × 1.0 m in size, the influence of solar irradiance G_rad to the increase in PV temperature and hence the PV performance of the thin film PV glass was evaluated. The thin film PV performance was simulated with the variation in solar irradiance at fixed outdoor temperature of 15°C. The size of the window was fixed at 1 m by 1 m. The aim is to keep the indoor temperature at 23°C. Figure 9 shows that, due to the electricity generation, the by-product of waste heat leads to the increase in the average temperature of the glazing which influences the electrical efficiency of the PV VG. Figures 9 and 10 shows that, the influence of the absorbed solar irradiance to the increase on PV temperature and hence the electrical efficiency of the PV VG was only obvious at solar irradiance above 300 W/m². On average, the drop in the electrical efficiency due to the presence of vacuum level is only 0.1%.

Figure 9

The variation in, the electrical efficiency ηele with the increase in solar irradiance at fixed indoor and ambient temperature for PV VG-4L at high and low vacuum level.

Figure 10

The variation in, temperature of the PV cells, with the increase in solar irradiance at fixed indoor and ambient temperature for PV VG-4L at high and low vacuum level.

Additionally, the temperature of the PV VG was also predicted with the change in the outdoor ambient temperature ranging from –10°C to the extreme of 40°C at fixed solar irradiance and indoor ambient condition of 700 W/m² and 23°C respectively Figure 11 shows that the electrical efficiency of the PV VG decreases by 0.032% with the increase in each degree of the outdoor ambient temperature. From Figures 9, 10, 11, we may conclude that influence of the excellent heat insulating properties of the thin film PV VG to the electrical efficiency of the thin film PV glass is significantly small and can be neglected.

Figure 11

The variation in the electrical efficiency with the change in outdoor ambient temperature at fixed solar irradiance and indoor temperature.

5. Discussions and Conclusions

An innovative Thin Film Photovoltaic Glazing with Vacuum Insulated Layer (PV VG-4L) is presented in this paper. A mathematical model was developed by taking into account all the parameters related to the individual component of a typical vacuum glazing unit which is the dominant in the design and also the fraction the solar irradiance absorbed by the different layers of PV VG-4L. To validate the mathematical model, a lab-scale prototype was manufactured and tested indoors using a calibrated hot box, and outdoors by installing the sample at a research house. Under controlled conditions, the overall U-value of the PV VG-4L was measured to be as low as 0.6 W/m²K. When investigated under real conditions, an obvious trend in glazing surface temperature variation with solar irradiance was obtained. During low to zero solar irradiance, the internal glazing surface temperature is on average higher in comparison to the external glazing surface temperature. However, as the solar irradiance increases, the by-product of the absorbed heat by the thin film PV glazing layer, has led to an increase in the external glazing surface temperature as the heat is hindered from being transferred indoors by the vacuum layer. At average solar irradiance of 600 W/m² the PV VG-4L can produced in total of 14W of power per m² of panel. Meanwhile, at low to zero solar irradiance (i.e. during the night), at outdoor ambient temperature of 14°C, the average U-value of the typical PV VG-4L was found to be as low as 0.6 W/m²K while maintaining the indoor ambient temperature at 30°C. Please note that high internal ambient temperature was obtained due to the use of heater at its maximum setting. For a conventional thin film PV glazing, to improve the thermal performance of the thin film PV glazing, a combination with a double-glazing unit is possible with the estimated U-value of 2.5–2.8 W/m²K depending on the type of gas used to provide the insulation in the air gap. Clearly, the vacuum layer introduced in the PV VG-4L design presented in this paper is better in performance with the slim configuration of the glazing unit as an additional benefit. The results also show that the PV VG-4L not only can produce power but also has high insulation properties when compared to a single thin film PV glazing with a typical U-value of 5 W/m²K; the U-value is higher by almost 90%. The promising U-value implies its range of potential applications can be improved depending on the energy needs and applications, such as for BIPV solar façade (PV curtain walling) in commercial buildings, greenhouses, skylight and conservatory.

Future Cities and Environment

Technical Articles

Performance Analysis of a Hybrid Thin Film Photovoltaic (PV) Vacuum Glazing

Abstract

1. Introduction

2. Design concept and heat transfer analysis of thin film PV vacuum insulated glazing (PV VG-4L)

3. Experimental validation of the theoretical model

3.1. U-value measurement using a TEC driven calibrated hot box

3.2. Performance analysis under real conditions

4. Parametric studies of PV VG-4L

5. Discussions and Conclusions

Acknowledgements

Competing Interests

References