Flexible thin-film solar cell is an efficient energy system on the surface of stratospheric airship for utilizing the solar energy. In order to ensure the normal operation of airship platform, the thermal control problem between the flexible thin-film solar cell and the airship envelope should be properly resolved. In this paper, a multilayer insulation material (MLI) is developed first, and low temperature environment test is carried out to verify the insulation effect of MLI. Then, a thermal heat transfer model of flexible thin-film solar cell and MLI is proposed, and the equivalent thermal conductivity coefficients of flexible thin-film solar cell and Nomex honeycomb are calculated based on the environment test and the temperature profile of flexible thin-film solar cell versus each layer of MLI. Finally, FLUENT is used for modeling and simulation analysis on the flexible thin-film solar cell and MLI, and the simulation results agree well with the experimental data, which validate the correctness of the proposed heat transfer model of MLI. In some way, our study can provide helpful support for further engineering applications of flexible thin-film solar cell.
1. Introduction
As a high-altitude platform, stratospheric airships are widely used in many important fields, especially in communication, broadcasting, remote sensing, scientific research, and so forth. Currently, United States, Japan, and South Korea are the major countries to develop stratospheric airships [1–10].
Solar energy is an ideal choice to provide power for high-altitude and long-endurance airships. This type of power system is actually a photovoltaic (PV) array coupled to an energy storage system [2]. The efficiency, mass, surface shape, and position of the PV array significantly influence the airship performance, whereas the PV array performance depends on its operation conditions such as latitude, time, temperature, and geometry. Therefore, the PV array and its operation conditions are critical to the airship design and capabilities. Besides, the PV array also has significant effects on the temperature profile of the airship hull, which can enhance the “superheat” or “supercool” of airship [9].
In the past decade, many investigations have been carried out on the energy system of the stratospheric airship. Harada et al. [5] conducted experiments on the thermal characteristics of a 35-meter-long low-altitude airship equipped with PV array, which was developed as a flying test bed for the development of stratospheric airship. Naito et al. [6] proposed the design and analysis of solar energy power system for stratospheric airship operations. Shi et al. [7] studied the solar energy collected by stratosphere airship’s surface. Zheng et al. [8] studied the solar power and power output characteristics of the solar cell array of a stratospheric airship. Li et al. [9] investigated the thermal characteristics of photovoltaic array of a stratospheric airship. In addition, Eskind et al. [10] mentioned a thermal control measure by using 7 mm thick insulation panels between the solar cell array and the envelope. Overall, current investigations mostly focus on the power system but seldom pay attention to the insulation of stratospheric airship solar array.
According to the characteristics of stratospheric airship envelope and flexible thin-film solar cell, this paper develops a multilayer insulation material (MLI) and then proposes a heat transfer model for the MLI. Further, in order to verify the insulation effect of MLI and get the temperature profile of flexible thin-film solar cell and each MLI layer, an experimental test is carried out to calculate the equivalent thermal conductivity coefficient of flexible thin-film solar cell and Nomex honeycomb. Finally, FLUENT is used for modeling and simulation analysis on the flexible thin-film solar cell and MLI.
2. Physical Model
MLI is a commonly used material in multilayer insulation for space vehicles. It can be composed of multilayer reflector screens, or laminated together with spacer. Since the flexible thin-film solar cell and the hull of stratospheric airship are both of flexible structures, considering the environment of stratosphere, MLI can be applied to insulation of the flexible thin-film solar cell. Yinbing et al. [11] developed an Al/PI/SiO2 multilayer insulation film with Al film as the heat reflective screen, silicone nano-composite coating as the low thermal conductivity layer, and polyimide (PI) as the flexible substrate. Weidong et al. [12] obtained the patent of a thin flexible insulation composite fabric. Ahluwalia et al. [13] developed a flexible multilayer composite insulation material which was made of metal foil (film) and flexible foam layer. These MLIs have good flexibility, light weight, and good heat insulation performance, but their mechanical strength is very weak. For example, the tensile strength of PI is 200 MP at 20°C and temperature rise can further decrease its tensile strength.
Based on existing flexible insulation structures and enhancement structures, we develop an MLI by virtue of PET, fiberglass mesh, epoxy resin film, and Nomex honeycomb. The PET film can be used as the substrate to which the fiberglass and Nomex honeycomb are attached. The fiberglass has high mechanical strength and can be used as the enhancement structure. The epoxy resin film can be used as adhesive to combine each layer together. The Nomex honeycomb is a hollow flexible insulation material and has a good capability of thermal insulation. The structure of MLI is shown as Figure 1.
Structure of MLI.
3. Analysis of Heat Transfer Process3.1. Heat Transfer Model
The structure is composed of the thin film solar cell, PET, fiberglass mesh, epoxy resin film, Nomex honeycomb, and airship hull. Thermal property of each layer is different, and solar radiation absorbed on different parts of the thin film solar cell is nonuniform. Thus, the thin film solar cell can be treated as a heat conductor with nonuniform inner heat source. In view of this, a separate heat-balance equation should be written for each layer to obtain a better estimation of the temperature. The heat transfer mechanism of the structure is shown as Figure 2, which can be expressed by the following five differential time-dependent equations:
(1)m1c1dT1dτ=Qsun-Qref-Qrad-Qconv-Qpower-Q1,2m2c2dT2dτ=Q1,2-Q2,3m3c3dT3dτ=Q2,3-Q3,4m4c4dT4dτ=Q3,4-Q4,5m5c5dT5dτ=Q4,5-Q5,base,
where
(2)Qi,i+1=A(Ti-Ti+1)δi/λi+δi+1/λi+1Q5,base=A(T5-Tbase)δ5/λ5+δbase/λbaseQsun=A*I0Qref=A*r*I0Qpower=A*a*I0,
where A is the area of flexible thin-film solar cell, m2; I0=1367w/m2 is the solar constant; r is the reflectivity of flexible thin-film solar cell; a is the conversion efficiency of solar cell; mi is the mass of Layer I, kg; and ci is the specific heat capacity of the material of Layer i, J/(kg·K).
MLI heat transfer model.
The radiation heat loss on the surface of thin film solar cell is given by
(3)Qrad=Aεσ(T14-T∞4),
where ε is the emissivity of flexible thin-film solar cell; σ is the constant of Stefan-Boltzmann, σ=5.67×10-8 W/(m2·K4); T1 is the temperature of the upside of flexible thin-film solar cell, K; and T∞ is the sky equivalent temperature, K.
The conductive heat loss is given by
(4)Qconv=Ah(T1-T∞),
where h is the convection heat transfer coefficient. h can be given by
(5)h=NuλairL,
where L is the characteristic length, m; λair is the thermal conductivity of air,
(6)λair=0.0241Tm273.150.9;Tm=0.5×(T1+T∞) is the characteristic temperature; and Nu is the Nusselt number:
(7)Nu=0.54(Gr·Pr)1/42×104≤Gr·Pr≤8×1060.15(Gr·Pr)1/38×106≤Gr·Pr≤1×1011,
where Pr is the Prandtl number, Pr=0.804-3.25×10-4×Tm; and Gr is the Grashof number, Gr=g·α·(T1-T∞)·L3/ϑ2, in which g is the acceleration of gravity; α is the coefficient of cubic expansion, 1/K; and ϑ is the kinematic viscosity of air, m2/s.
To calculate the thermal conductivity of flexible thin-film solar cell and Nomex honeycomb accurately, an experiment is designed to obtain the equivalent thermal conductivity coefficient.
3.2. Thermal Conductivity of Flexible Thin-Film Solar Cell
The surface of flexible thin-film solar cell is exposed to ambient air. The exposed area is A and the temperature of ambient air is T∞. When the solar radiation intensity is I0, the upside and downside temperatures of flexible thin-film solar cell are Tup and Tdown when the solar radiation intensity is, respectively.
The structure of flexible thin-film solar cell is shown in Figure 3.
Structure of flexible thin-film solar cell.
The flexible thin-film solar cell is composed of PET, EVA, and single-junction amorphous silicon cells, which can be treated as a multilayer model. Parameters of the flexible thin-film solar cell are listed in Table 1.
Parameters of each layer of the thin-film solar cell.
Material
Thickness (10−6 m)
Thermal conductivity(W/(m·K))
Emissivity
Transmittance
Photoelectric conversionefficiency
PET
18
0.22
0.92
90%
—
EVA
7
0.25
—
91%
—
A-Si
0.2
—
—
—
2.22%
EVA
10
0.25
—
91%
—
PET
65
0.22
0.92
90%
—
To simplify the energy model of flexible thin-film solar cell, it is assumed that a portion of solar flux is absorbed by the solar cell while the other part is reflected back to the atmosphere. The surface has convective heat transfer with the atmosphere. The flexible thin-film solar cell has conductive heat transfer with the MLI. The mechanism of heat transfer of the flexible thin-film solar cell is shown in Figure 4.
Mechanism of heat transfer of the flexible thin-film solar cell.
According to the mechanism of heat transfer of the flexible thin-film solar cell, steady-state energy equation of flexible thin-film solar cell can be established. Thus, when the Tup, Tdown, T∞, and I0 are obtained, we can calculate the thermal conductivity of flexible thin-film solar cell:
(8)AI0=Qref+Qrad+Qconv+Qpower+Qcond,(9)Qcond=AλTup-Tdownδ,
where λ is the equivalent thermal conductivity of flexible thin-film solar cell; and δ is the thickness of flexible thin-film solar cell.
Since each layer of the flexible thin-film solar cell is tightly laminated together, the thermal contact resistance can be ignored. Thus, the thermal conductivity can be given by
(10)δλ=δ1λ1+δ2λ2+δ3λ3+δ4λ4+δ5λ5,
where δ1, δ2, δ3, δ4, and δ5 are the thickness and λ1, λ2, λ3, λ4, and λ5 are the thermal conductivity of each layer of the flexible thin-film solar cell. Then, the thickness of the flexible thin-film solar cell can be given by
(11)δ=δ1+δ2+δ3+δ4+δ5.
The conductive heat can be calculated from the steady-state energy equation, and then we can get the thermal conductivity of the flexible thin-film solar cell according to (5)–(8).
3.3. Equivalent Thermal Conductivity of Nomex Honeycomb
Nomex honeycomb is used at the bottom of the MLI. The Nomex honeycomb layer has conductive heat transfer with the adjacent epoxy resin film layer. To establish the calculation model of equivalent thermal conductivity of Nomex honeycomb, it is necessary to study the cell unit of honeycomb, as shown in Figures 5 and 6.
Nomex honeycomb.
Cell unit of honeycomb.
The side length of Nomex honeycomb is l=1.83×10-3 m, the height is H=5.5×10-3 m, and the wall thickness of single layer is t=0.2×10-3 m. According to the results of Liang et al. [14], convection heat transfer of air inside honeycomb can be ignored. Hence, heat transfer model of the honeycomb mainly consists of solid heat conduction of the sandwich layer of cell wall, air heat conduction inside the honeycomb, and radiation heat transfer between honeycomb cell wall and the upper and lower panels. The equivalent thermal conductivity can be given by
(12)λn=λfΔAA+λa1-ΔAA+λr,
where λf is the thermal conductivity of the thin foil parent material of honeycomb; λa is the thermal conductivity of air inside the honeycomb; λr is the radiation thermal conductivity of honeycomb cell wall and the upper and lower panels; ΔA is the cross-sectional area of honeycomb core material, ΔA=16lt; and A is the cross-sectional area of cell unit, A=63l2.
Thermal conductivity of air inside the honeycomb is given by
(13)λa=0.0241Ta273.150.9,
where Ta=T1+T2/2 is the characteristic temperature, in which T1 and T2 are the temperatures of the upside and down side of honeycomb.
Radiation thermal conductivity of honeycomb is given by Swarm and Pittman [15]:
(14)λr=0.664η+0.3-0.69Hε1.63η+1-0.89×σ(T1+T2)(T12+T22),
where η=H/l is the ratio of the height and side length of the honeycomb; ε is the emissivity inside the honeycomb; σ=5.67×10-8 W/m2·K4 is Stefan-Boltzmann constant; and T1 and T2 are unknown. In order to determine T1 and T2, iterative calculations are required. In this paper, according to the experimental test data, T1=T5 and T2=Tbase.
4. Experimental Test
In this section, an experimental test is carried out to verify the insulation effect of the MLI, as well as to calculate the thermal conductivity of flexible thin-film solar cell and Nomex honeycomb accurately. Besides, the heat transfer model is also validated by the experimental test.
The dimensions of environmental cabin are D×L=1.8m×3.5m. Two pieces of flexible thin-film solar cell are put in the environmental cabin, and their dimensions are 0.24m×0.22m. The temperature of the environmental cabin is −56°C, and the pressure is 5 kPa. Radiant intensity of solar simulator is 1260W/m2. Test specimen 1 is just flexible thin-film solar cell, which is used to test the temperatures of the upside and downside of flexible thin-film solar cell, as shown in Figure 7. Test specimen 2 is composed of flexible thin-film solar cell and MLI, which is used to test the temperatures of the upside and downside of flexible thin-film solar cell as well as the temperature of the downside of MLI, as shown in Figure 8.
Test specimen 1.
Test specimen 2.
Taking test specimen 2 as an example, temperature sensors 1 and 2 are at the upside of the thin film solar cell, temperature sensors 3 and 4 are at the downside of the thin film solar cell, and temperature sensors 5 and 6 are at the downside of the MLI. Here, sensor 1 corresponds to sensors 3 and 5, and sensor 2 corresponds to sensors 4 and 6. Figure 9 illustrates the test circuit diagram.
Test circuit diagram.
Temperature data are recorded by temperature inspection instrument until temperature change <3°C/10 min. The test data of test specimen 1 and test specimen 2 are shown in Figures 10 and 11, respectively.
Test temperature data of test specimen 1.
Test temperature data of test specimen 2.
T1up1 and T1up2 are the temperatures of the upside of test piece 1. T1down1 and T1down2 are the temperatures of the downside of test piece 1. T1up(T1up=(T1up1+T1up2)/2) and T1down(T1down=(T1down1+T1down2)/2) are the average temperatures of the upside and downside of thin film solar cell, respectively.
T2up1 and T2up2 are the temperatures of the upside of test piece 2. T2middle1 and T2middle2 are the temperatures of the downside of thin film solar cell of test piece 2. T2down1 and T2down2 are the temperatures of the downside of the MLI of test piece 2. T2up(T2up=(T2up1+T2up2)/2) and T2middle(T2middle=(T2middle1+T2middle2)/2) are the average temperatures of the upside and downside of thin film solar cell, respectively. T2down(T2down=(T2down1+T2down2)/2) is the average temperature of the downside of the MLI.
5. Computational Simulation and Analysis
According to the test data of test specimen 1, the steady temperatures of the upside and downside of the flexible thin-film solar cell are 65.5°C and 63°C, respectively. According to (5)–(8), the thermal conductivity of flexible thin-film solar cell is λ=0.2246 W/(m·k). Figure 11 shows that the steady temperature of the downside of the Nomex honeycomb is 25.5°C. According to the test data of test specimen 2 and (1)–(7) and (12)–(14), through iterative calculation, the final thermal conductivity of Nomex honeycomb is λ=0.021 W/(m·k), as shown in Figure 12.
Thermal conductivity of Nomex honeycomb.
Figure 13 shows the simulation temperature profile of the multilayer insulation. T1 = 61°C is the average steady temperature of flexible thin-film solar cell, which is measured from test specimen 2. T2, T3, T4, and T5 are the simulation temperatures of PET, fiberglass mesh, epoxy resin film, and Nomex honeycomb, respectively, and their final steady temperatures are 59.5°C, 58°C, 39.3°C, and 24.4°C. It can be seen that the temperatures of flexible thin-film solar cell, PET, and fiberglass mesh are more or less the same, whereas the temperatures of the epoxy resin film and Nomex honeycomb decrease dramatically due to their large thermal resistance caused by the hollow structure.
Simulation temperature of each layer.
Figure 14 shows the comparison of simulation results with relevant test data of Tbase. The absolute error of simulation results relative to relevant test data is less than 2°C; namely, the relative error of the temperature of Tbase is within 5%. Therefore, the simulation results are accurate and reliable.
Simulation result and relevant test data of Tbase.
Finally, to investigate the applicability of the numerical model of MLI, FLUENT is used for modeling and simulation analysis on the flexible thin-film solar cell and MLI. To simplify the model, a 24cm×22cm×0.6cm three-dimensional model of flexible thin-film solar cell and MLI is put in a fixed computational domain.
Solar load model is used to calculate the effect of solar irradiation on temperature distribution of the thin film battery. This model includes direct solar irradiation, diffuse solar irradiation, and ground-reflected solar radiation. Two options, Fair Weather Conditions method and Theoretical Maximum method, are, respectively, put forward for computing the solar load in the software, and Theoretical Maximum method is proved to be more appropriate for stratosphere and outer surface where little atmospheric extinction happens. The equation for this method is
(15)Ed=cI,
where c is the correction factor used to modify the reduction of solar irradiation due to aerosphere.
The material of the wall around the computational domain is supposed to be copper (absorptance and emittance are all 0.98). Temperature of the wall is −56°C. Solar radiation enters the computational domain from the top, and the bottom is set as outflow. Air in the computational domain is set as incompressible-ideal-gas. DO radiation model and DO irradiation model are chosen. Direct solar radiation and diffuse solar radiation are calculated by solar calculator. After a 50 min calculation, the accurate temperatures of the upside and downside of the flexible thin-film solar cell and MLI can be obtained, as shown in Figure 15.
Temperature distribution on upside of the thin film cell and downside of the MLI.
Figure 15 shows that the highest temperature of the surface of flexible thin-film solar cell is 336 K (63°C), and the highest temperature of the downside of MLI is 295 K (22°C), with the temperature decreasing from center to around. For the upside of flexible thin-film solar cell, the absolute error of simulation result relative to test data is about 2.4°C, and the relative error is within 10%. For the downside of MLI, the absolute error of simulation result relative to test data is about 2°C, and the relative error is within 3%. Therefore, the simulation results are accurate and reliable.
6. Conclusions
This paper developed a multilayer insulation material (MLI) according to the characteristics of stratospheric airship hull and flexible thin-film solar cell. The insulation effect of MLI was verified by environmental test, and its heat transfer model and thermal characteristics were analyzed. The main conclusions are as follows.
The equivalent thermal conductivity coefficients of flexible thin-film solar cell and Nomex honeycomb are calculated as 0.2246 W/(m·k) and 0.021 W/(m·k), respectively.
MLI has a good heat insulation performance of 36.5°C. The temperature difference between the layers of flexible thin-film solar cell, PET, and fiberglass mesh is small due to the small thickness and heat resistance, whereas the temperatures of the epoxy resin film and Nomex honeycomb decrease dramatically due to their large thermal resistance caused by the hollow structure.
The heat transfer model established in the paper and simulation results of MLI agree well with the environmental data. It can be concluded that the simulation results are accurate and that the heat transfer model is reliable.
FLUENT is used for modeling and simulation analysis on the flexible thin-film solar cell and MLI, and the results agree well with the environmental test data.
In some way, our study can provide support for further engineering application of flexible thin-film solar cell. Our future research will focus on the effects of MLI on the thermal characteristics and output efficiency of the photovoltaic array, and the effects of MLI on thermal characteristics of the airship.
Conflict of Interests
The authors declare that they do not have any commercial or associative interest that represents a conflict of interests in connection with the work submitted.
Acknowledgments
This work was supported by the National Natural Science Foundation of China under Grant no. 51307004 and Shanghai Aerospace Science and Technology Innovation Fund no. SAST201268.
LiaoL.PasternakI.A review of airship structural research and development2009454-5839610.1016/j.paerosci.2009.03.0012-s2.0-68949206411ColozzaA.Initial feasibility assessment of a high altitude long endurance airship2003NASA CR-21272EguchiK.YokomakuY.MoriM.Feasibility study program on stratospheric platform airship technology in Japan199999-3912KimD.-M.LeeY.-G.KangW.-G.LeeJ.-W.YeomC.-H.Korea stratospheric airship program and current resultsProceedings of the 3rd Annual Aviation Technology, Integration, and Operations (ATIO) Forum (AIAA '03)November 2003200367822-s2.0-84894487291HaradaK.EguchiK.SanoM.SasaS.Experimental study of thermal modeling for stratospheric platform airships20032003-6833NaitoH.EguchiK.HoshinoT.Design and analysis of solar power system for SPF airship operationsInternational Balloon Technology Conference1999AIAA Paper 1999-3351ShiJ.ZhuX.SuiS.Modeling and analysis of solar energy collected by stratosphere airship's surface20072486467ZhengW.SongQ.LiY.WuZ. H.HuL. Y.Computation and analysis of power generated by the solar cell array of a stratospheric airship20103141224123010.3873/j.issn.1000-1328.2010.04.0462-s2.0-77953282733LiX.FangX.DaiQ.Research on thermal characteristics of photovoltaic array of stratospheric airship20114841380138610.2514/1.c0312952-s2.0-79961061752EskindLarryG.AanestadMultilayered flexible fire-retardant isolative sheets comprising laminates of insulation layers and fire-retardant glass fiber layers and aluminum foil filmsWO Pat, 2001043972, 2001YinbingC.LanfangY.GuangmingW.Study on Al/PI/SiO2 multilayer insulation2003212, article 221WeidongY.QinhuaY.HeH.ZhipingM.NanliangC.QuanZ.Preparation and applications of an adiabatic thin soft compound fabricCN Patent ZL 200410066650.2.2004-09-24, 2004 (Chinese)AhluwaliaY.KiikM.KarolT.Composite material2004US20040229053 A1LiangW.LiuZ.MaiH.YangJ.Theoretical calculation for thermal properties of honeycomb structure based on the radiation and conduction coupling model20083543136SwarmR. T.PittmanC. M.Analysis of effective thermal conductivities of honeycomb-core and corrugated-core sandwich panels1961D-714