# Effects of Nonuniform Incident Illumination on the Thermal Performance of a Concentrating Triple Junction Solar Cell

Authors: Fahad Al-Amri; Tapas Kumar Mallick
Journal: International Journal of Photoenergy
 (2014)

---

## Abstract

A numerical heat transfer model was developed to investigate the temperature of a triple junction solar cell and the thermal characteristics of the airflow in a channel behind the solar cell assembly using nonuniform incident illumination. The effects of nonuniformity parameters, emissivity of the two channel walls, and Reynolds number were studied. The maximum solar cell temperature sharply increased in the presence of nonuniform light profiles, causing a drastic reduction in overall efficiency. This resulted in two possible solutions for solar cells to operate in optimum efficiency level: (i) adding new receiver plate with higher surface area or (ii) using forced cooling techniques to reduce the solar cell temperature. Thus, surface radiation exchanges inside the duct and Re significantly reduced the maximum solar cell temperature, but a conventional plain channel cooling system was inefficient for cooling the solar cell at medium concentrations when the system was subjected to a nonuniform light distribution. Nonuniformity of the incident light and surface radiation in the duct had negligible effects on the collected thermal energy.


---

## Body

## 1. Introduction

Nonuniformity of the incident illumination in concentrated photovoltaic (CPV) systems occurs as a result of using concentrating optics that reflect or refract light flux on the solar cell surface. The optimal optical system design would concentrate radiation on the solar cell surface in a uniform distribution. However, existing concentrators condense the majority of the incident flux onto a limited area of the cell surface, and most areas of the cell receive the remaining small amount of radiation. This situation produces a highly nonuniform irradiance distribution on the photovoltaic cells. This nonuniformity has two main effects on the solar cell: (i) electrical impacts, such as high ohmic drops and the flow of an internal current, and (ii) thermal impacts, such as increased temperatures in some regions of the cell that cause hot spots and drastically reduce the overall efficiency and power output. This reduction of power limits the viability of photovoltaic cells as an alternative energy source.Previous studies reported thermal behavior of the CPV solar cells, such as Min et al. [1] and Cotal and Frost [2] who developed a theoretical thermal model for predicting the temperature of a multijunction solar cell based on passive cooling systems. Vincenzi et al. [3] used a silicon wafer with microchannels circulating water directly beneath the cells as an active cooling system. Moshfegh and Sandberg [4] performed both a numerical study and an experimental study of the flow and heat transfer characteristics of natural air convection behind solar cells. The results demonstrated that surface radiation significantly affects the temperature and conversion efficiency of the solar panel. Bhargava et al. [5], Garg and Adhikari [6], and Hegazy [7] analyzed the performance of hybrid photovoltaic/thermal (PV/T) air heating collectors and reported the effects of different design and controlling parameters on system performance. Al-Amri and Mallick [8, 9] recently developed a numerical heat transfer model for predicting the maximum cell temperatures of multijunction concentrating solar cell systems that were actively cooled by water-forced and air-forced convection. The maximum cell temperature was strongly dependent on the inlet velocity and the channel width. Teo et al. [10] experimentally investigated an active cooling system for photovoltaic modules. A parallel array of ducts with an inlet/outlet manifold was attached to the back of the PV panel. When an air-based active cooling mechanism was employed, the temperature dropped effectively, leading to a 12% to 14% increase in the efficiency of the solar cells. Kumar et al. [11] presented a numerical heat transfer model for a novel concentrating photovoltaic design for the active solar panel initiative system (ASPIS). The incorporation of microchannels is an innovative way to remove a large amount of heat from a small cell surface that is associated with a high concentration [12, 13]. Tuckerman and Pease [14] pioneered the use of the microchannel heat sink for cooling planar integrated circuits. Several additional studies [15–17] were performed to analyze the heat transfer performance and the pressure drop of flows in microchannels. Qu and Mudawar [18] investigated the pressure drop and heat transfer performance of a laminar flow in a microchannel both experimentally and numerically. The microchannel was 213 μm wide and 713 μm deep. The results of this study indicated that the flow in this system is still governed by conventional equations for the conservation of mass, momentum, and energy. Micheli et al. [19] recently presented an overview of micro- and nanotechnologies that are applicable to passive CPV cooling and associated manufacturing technologies (e.g., monolithic applications). The authors reported that carbon nanotubes and high-conductive coating are the technologies that offer the best CPV cooling performance. This technological review was also critically assessed.The problem of nonuniform illumination patterns has been studied extensively. However, most of these studies focused on the effects of nonuniformity on electrical performance [12]. Few studies investigated the thermal impacts of nonuniformity in radiation profiles. Franklin and Conventry [20] demonstrated both numerically and experimentally that a silicon solar cell under nonuniform focused illumination experiences a drop in both open-circuit voltage and efficiency compared to a cell under uniform illumination. In addition, the results indicated that the temperature profile closely resembles a Gaussian distribution. Coventry et al. [21] measured the thermal and electrical output from a concentrating PV/thermal collector at 25–30 suns. While the overall concentration ratio was approximately 30x, the illumination flux intensity across the width of the solar cell was far from even. Domenech-Garret [22] investigated the effects of a combination of nonuniform illumination and nonuniform temperature on a silicon solar cell. He used Gaussian and inverse Gaussian profiles to describe the temperature in the cell, which simulated a general cell cooling device. Chemisana and Rosell [23] studied the effects of Gaussian and anti-Gaussian temperature profiles on the cell’s electrical parameters under different types of concentrated illumination both numerically and experimentally. The electric conversion efficiency increases under a Gaussian temperature profile and decreases when the cell is subjected to an inverse Gaussian temperature curve. The authors concluded that the temperature profile could be tailored to maximize efficiency for determined radiation conditions.The studies presented above failed to adequately address the effects of nonuniformity of the incident flux on the temperature profiles of CPV systems. Thus, further studies of heat transfer characteristics under nonuniform flux are necessary to understand the effect of this nonuniformity on the overall efficiency of the solar cell. The objective of this investigation was to extend the work of Al-Amri and Mallick [8] by assessing the effect of nonuniform incident light on the temperature distribution of a concentrated triple junction solar cell.A heat transfer model was developed to simulate the temperature distribution on the concentrated solar cells under a nonuniform incident light. The results of this study are important for determining the overall efficiency of a CPV system and for promoting the development of a cooling system and new concentrators.


## 2. Mathematical Formulations

The geometry of the physical problem and the coordinate system considered in this study are presented in Figure1. The components consist of a triple junction GalnP/GaAs/Ge solar cell, a Cu-Ag-Hg front contact to the solar cell, a 2 mm glass cover, a 1.5 mm thick aluminum back plate, and a 0.5 mm adhesive material used to attach the solar cell to the cover glass and the rear aluminum plate. The cooling air is forced to flow within the ducts behind the back plate, whose external wall is assumed to be adiabatic. The two walls of the duct are assumed to be gray, opaque, and diffuse surfaces. The inlet and exit channel areas are assumed to be black at the ambient temperature (
T
∞
=
27
°
C
) and the exit bulk temperature (
T
e
), respectively. The efficiency of the triple junction solar cell (
η
) is 0.38 and the reflectivity is assumed to be zero for all solar cell components except for the front contact metal which is assumed to be 0.85. Table 1 shows the dimensions and the thermophysical properties of the components used in the model.Table 1
Thermophysical properties assumed in the simulation.
Component
Material
Length (mm)
Thickness (mm)
Thermal conductivity W/m·K
Absorptivity
Solar cell
GuInp
10
0.1
73
0.04
GuAs
10
0.2
65
0.40
Ge
10
0.2
60
0.50
Front contact
Cu + Ag + Hg
2
0.025
300
0.15
Glass cover
Low-iron glass
10
2
12
0.04
Adhesive
10
0.5
50
0.02
Back plate
Aluminium
10
1.5
238Figure 1
Schematic of the geometry of the employed system.The illumination distribution incident on the cell surface in a concentrator system consisting of reflective troughs is highly non-uniform and takes the shape of a Gaussian profile [12]. Thus, the illumination intensity formula adopted in the current model is
(1)
I
(
z
)
=
2
C
1
0
4
SD
2
π
e
-
(
2
(
z
-
z
0
)
2
)
/
S
D
2
,
where z
0 is the position of the maximum illumination on the cell surface and SD is the standard deviation, which measures how far incident illumination is spread out.The intensity of incident lightI
(
z
) can be written in the nondimensional form as follows:
(2)
I
-
(
Z
)
=
I
(
z
)
I
avg
,
where I
avg is the average incident illumination across the entire cell surface.The physical problem depicted in Figure1 is governed by equations of energy conservation in both the multiwalls (solar cell assembly) and the fluid inside the duct and by equations of continuity, momentum in the fluid region, and radiation constraint along each of the two duct surfaces as a result of surface radiation exchange between the two walls of the channel. Using typical boundary layer assumptions and neglecting the axial conduction of heat in both the solid and the fluid regions, the corresponding equations governing the conjugate laminar mixed flow and the heat transfer in the entry region of the two parallel plates can be expressed in the following dimensionless forms.Continuity Equation.Consider
(3)
∂
V
∂
Y
+
∂
U
∂
Z
=
0
.Momentum Equation.Consider
(4)
V
∂
U
∂
Y
+
U
∂
U
∂
Z
=
-
d
P
d
Z
+
G
r
*
Re
(
θ
-
θ
∞
)
+
∂
2
U
∂
Y
2
.Energy Equation for the Fluid.Consider
(5)
V
∂
θ
f
∂
Y
+
U
∂
θ
f
∂
Z
=
1
Pr
∂
2
θ
f
∂
Y
2
.Energy Equation for Each Solid Material in Solar Cell Assembly. Consider
(6)
∂
2
θ
s
∂
Y
2
=
0
.
The integral form of the continuity equation (1) can be written in the following form [24]:
(7)
F
=
∫
0
1
U
d
Y
=
1
.
The radiation constraint equations can be derived by applying the first law of thermodynamic per unit area of the surface on the channel plate element [25] which yields the following two equations.Surface 1  (
Y
=
0
). Consider
(8)
N
rad
θ
∞
4
[
1
2
-
(
Z
Re
2
)
(
1
+
Z
2
R
e
2
)
-
1
/
2
]
+
N
rad
θ
e
4
[
1
2
-
(
(
L
-
Z
)
Re
2
)
(
1
+
(
L
-
Z
)
2
R
e
2
)
-
1
/
2
]
+
∫
0
L
{
1
-
ε
2
ε
2
(
∂
θ
f
∂
Y
|
Y
=
1
-
-
∂
θ
s
∂
Y
|
Y
=
1
+
)
+
N
rad
θ
w
2
4
(
Z
)
}
*
1
2
[
1
+
R
e
2
(
Z
-
Z
′
)
2
]
(
-
3
/
2
)
Re
d
Z
′
=
1
-
ε
1
ε
1
[
-
∂
θ
f
∂
Y
|
Y
=
0
+
]
+
N
rad
θ
w
1
4
(
Z
)
-
∂
θ
f
∂
Y
|
Y
=
0
+
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
Where:
N
rad
=
σ
q
avg
3
b
4
k
f
4
.Surface 2
(
Y
=
1
). Consider
(9)
K
R
s
1
-
f
∂
θ
s
∂
Y
|
Y
=
1
+
+
N
rad
θ
∞
4
[
1
2
-
(
Z
Re
2
)
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
×
(
1
+
Z
2
R
e
2
)
-
1
/
2
1
2
-
(
Z
Re
2
)
]
+
N
rad
θ
e
4
[
1
2
-
(
(
L
-
Z
)
2
Re
)
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
×
(
1
+
(
L
-
Z
)
2
R
e
2
)
-
1
/
2
1
2
-
(
(
L
-
Z
)
2
Re
)
]
+
∫
0
L
{
1
-
ε
1
ε
1
(
-
∂
θ
∂
Y
|
Y
=
0
+
)
+
N
rad
θ
w
1
4
(
Z
)
}
*
1
2
[
1
+
R
e
2
(
Z
-
Z
′
)
2
]
(
-
3
/
2
)
Re
d
Z
′
=
1
-
ε
2
ε
2
[
∂
θ
f
∂
Y
|
Y
=
1
-
-
∂
θ
s
∂
Y
|
Y
=
1
+
]
+
N
rad
θ
w
2
4
(
Z
)
-
∂
θ
f
∂
Y
|
Y
=
1
-
.When surface radiation is absent (i.e.,ε
=
0), (7) and (8) are reduced to
(10)
∂
θ
f
∂
Y
|
Y
=
0
=
0
,
K
R
s
1
-
f
∂
θ
s
∂
Y
|
Y
=
1
+
=
∂
θ
f
∂
Y
|
Y
=
1
-
.
The conjugate convection field equations are subject to the following dimensionless boundary conditions.In the Fluid Region. For Z
=
0 and 0
<
Y
<
1,(11a)
U
=
1
,
V
=
P
=
0
,
θ
=
θ
∞
.ForZ
>
0 and Y
=
0,
(11b)
U
=
V
=
0
.ForZ
>
0 and Y
=
1,
(11c)
U
=
V
=
0
.In the Solid Region (solar cell assembly). In the present model, we assumed that some of the incident light is absorbed by each surface of solar cell assembly layers based on the value of absorptivity and converted into thermal and/or electrical energy. Then, heat is diffused to the next layer by conduction. Each layer of the solar cell is in the order of 100 to 200 microns; therefore, the skin depth of each layer is assumed to be thin and the radiation effects are confined to the surface layer. The boundary conditions in this region can be deduced by applying the continuity of temperature and heat flux at the solid-solid interface. Taking a surface energy balance at each solid-solid interface (see Figure 2) leads to the following boundary condition equations.Figure 2
Energy fluxes in solar cell assembly.ForZ
>
0 and Y
=
Y
s
(
i
)
-
(
i
+
1
), i
=
1 and 2,
(11d)
∂
θ
s
∂
Y
|
Y
=
Y
s
(
i
)
-
(
i
+
1
)
-
=
KR
s
(
i
+
1
)
-
(
i
)
∂
θ
s
∂
Y
|
Y
=
Y
s
(
i
)
-
(
i
+
1
)
+
.ForY
=
Y
s
(
i
)
-
(
i
+
1
) and 0
<
Z
<
Z
c or (
Z
c
+
L
c
)
<
Z
<
L, i
=
3 and 4,
(11e)
∂
θ
s
∂
Y
|
Y
=
Y
s
(
i
)
-
(
i
+
1
)
-
=
k
R
s
(
i
+
1
)
-
(
i
)
∂
θ
s
∂
Y
|
Y
=
Y
s
(
i
)
-
(
i
+
1
)
+
+
(
1
-
η
)
α
s
i
(
1
-
α
s
8
-
α
s
7
)
I
-
(
Z
)
K
R
s
i
-
f
.ForY
=
Y
s
5
-
7 and 0
<
Z
<
Z
c or (
Z
c
+
L
c
)
<
Z
<
L,
(11f)
∂
θ
s
∂
Y
|
Y
=
Y
s
5
-
7
-
=
k
R
s
7
-
5
∂
θ
s
∂
Y
|
Y
=
Y
s
5
-
7
+
+
(
1
-
η
)
α
s
5
(
1
-
α
s
8
-
α
s
7
)
I
-
(
Z
)
K
R
s
5
-
f
.ForY
=
Y
s
(
i
)
-
(
i
+
1
) and Z
c
≤
Z
≤
(
Z
c
+
L
c
), i
=
3
,
4, and 5,
(11g)
∂
θ
∂
Y
|
Y
=
Y
s
(
i
)
-
(
i
+
1
)
-
=
k
R
s
(
i
+
1
)
-
(
i
)
∂
θ
∂
Y
|
Y
=
Y
s
(
i
)
-
(
i
+
1
)
+
.ForY
=
Y
s
6
-
7 and Z
c
≤
Z
≤
(
Z
c
+
L
c
),
(11h)
∂
θ
∂
Y
|
Y
=
Y
s
6
-
7
-
=
k
R
s
7
-
6
∂
θ
∂
Y
|
Y
=
Y
s
6
-
7
+
+
α
s
6
(
1
-
α
s
8
-
α
s
7
)
I
-
(
Z
)
k
R
s
6
-
f
.ForZ
>
0 and Y
=
Y
s
7
-
8,
(11i)
∂
θ
∂
Y
|
Y
=
Y
s
7
-
8
-
=
k
R
s
8
-
7
∂
θ
∂
Y
|
Y
=
Y
s
7
-
8
+
+
α
s
7
(
1
-
α
s
8
)
I
-
(
Z
)
k
R
s
7
-
f
.ForZ
>
0 and Y
=
Y
s
8,
(11j)
∂
θ
∂
Y
|
Y
=
Y
s
8
=
α
s
8
I
-
(
Z
)
k
R
s
8
-
f
,
where k
R
s
(
i
+
1
)
-
i
=
k
s
i
+
1
/
k
s
i and k
R
s
i
-
f
=
k
s
i
/
k
f.


## 3. Method of Solution

The conjugate convection field equations (i.e., continuity, momentum, and energy), the two radiation constraint equations, and the corresponding boundary conditions are transformed into finite difference equations. The conjugate convection field equations are initially solved using a numerical marching technique. The two radiation constraint equations are then solved iteratively to update the wall temperatures using the Gauss-Seidel technique. These updated wall temperatures are used to resolve the conjugate convection field equations. The process is then repeated until convergence is achieved. The details of the solution procedure and the validation of the numerical code can be found in Al-Amri and Mallick [9].


## 4. Results and Discussion

Figures3(a) and 3(b) show the incident illumination profiles obtained using (1) for the concentration ratio C
=
100 at three selected standard deviation values (SD = 0.0005, 0.001, and 0.002) and three selected values for the position of the maximum illumination on the cell surface (
z
0
=
2
,
5
,
8
mm
), respectively. These profiles were used as input data in the code to study the effects of these parameters on the temperature distribution of the solar cell and the fluid bulk temperature. In addition, the effect of the emissivity of the two duct walls and the Reynolds number in the presence of nonuniform incident light was investigated.Incident illumination profiles forC
=
100 at (a) different values of SD and (b) different values of z
0.
(a)
(b)


### 4.1. Effect of a Nonuniform Light Profile on the Solar Cell Temperature

Figure4 shows the variations of solar cell temperature with z
0
=
5
mm and with the axial distance from the channel entrance for both uniform and nonuniform incident light at SD = 0.0005, 0.001, and 0.002. For a uniform distribution, the results demonstrate that the solar cell temperature first increases sharply near the entrance and then gradually increases until reaching its maximum value at the channel exit. In addition, there is a sudden increase in cell temperature at z
=
4
mm and a sudden decrease at z
=
6
mm. These changes occur due to the fact that no current is generated in this portion of the solar cell and all of the light is converted to heat due to the presence of the Cu-Ag-Hg front contact. However, this increment is small (≈7°C) and has minor effects on the overall efficiency. For a nonuniform distribution, the temperature of the first quarter of the solar cell remains constant at ambient temperature, then rapidly increases, and reaches its maximum value at the center of the cell. The temperature then drops sharply and gradually in the fourth quarter of the solar cell where the temperatures are higher than those observed in the first quarter. The profile of the temperature becomes thinner and the maximum temperature increases as the SD decreases, as shown in Figure 4. In addition, the maximum temperature of the solar cell is much higher when using a nonuniform light distribution than when using a uniform light distribution. In this example, the maximum temperature obtained using a uniform light distribution was 240°C, which indicates that the solar cell can be operated with a reasonable efficiency as the modern multijunction solar cells can adequately be operated at this level of temperature [25, 26]. In contrast, the results obtained using a nonuniform light distribution indicate that the cell efficiency will be drastically reduced due to the high temperature of some portions of the cell surface.Figure 4
Variations in solar cell temperature obtained using different values of SD, whenz
0
=
0.5, Re
=
500, ε
=
0.95, and C
=
100.The effects of the position of maximum incident light on the solar cell temperature distribution are shown in Figure5. A hot spot formed on the portion of the cell surface that was exposed to a large amount of light. This hot spot is slightly hotter when z
0
=
5
mm (i.e., at the center of the solar cell) due to the presence of the front contact metal in this region. Figures 6(a)–6(d) show the effect of the emissivity of the two duct walls on the temperature profile of the solar cell for both uniform and nonuniform incident light when z
0
=
2, 5, and 8 mm and when ε
=
0 (i.e., pure forced convection), 0.5, and 0.95. The results presented in these figures demonstrate that the maximum solar cell temperature decreases significantly as surface emissivity increases. However, the presence of surface radiation inside the channel does not change the shape of the temperature distribution for the nonuniform light profile. In this example, as emissivity increased from ε
=
0 to 0.95, the maximum solar cell temperature was reduced by approximately 35% when using nonuniform illumination and 20% when using uniform illumination. The large contribution of surface radiation to the nonuniform light profile is due to the increase in the effective temperature differential for surface radiation between the two duct walls. Nevertheless, the maximum temperature of the solar cell when using nonuniform light was still higher than that obtained using a uniform light distribution. For the uniform light profile, surface radiation exchange between the two duct walls and between the heated duct wall and the inlet and exit areas of the channel leads to a reduction of maximum cell temperature from 305°C to 240°C, which enables the solar cell to operate with a reasonable efficiency.Figure 5
Variations in solar cell temperature obtained using different values ofz
0, when SD = 0.001, Re = 500, ε
=
0.95, and C
=
100.Variations in solar cell temperature obtained using different values ofε, when SD = 0.001, Re = 500, and C
=
100.
(a)
z
0
=
2
mm
(b)
z
0
=
5
mm
(c)
z
0
=
8
mm
(d)
Uniform fluxThe results presented in Figures7(a)–7(c) illustrate the variations in solar cell temperature that occur with changes in axial distance from the channel entrance when using uniform and nonuniform light distributions, three selected values of Reynolds number (500, 1000, and 2000), and a given value of ε
=
0.95. The cell temperature decreases with increasing Re, but the profile maintains its general shape. The percent reductions of maximum solar cell temperature that occurred as Reynolds number increased by 300% from 500 to 2000 were 20.7%, 13.4%, and 13.7% for z
0
=
2, 5, and 8 mm, respectively. This value was approximately 22.5% when a uniform light distribution was employed. These results indicate that, for nonuniform light profiles, increasing flow rate is more useful when the maximum incident light is located near the channel entrance because the difference between the wall temperature and the fluid temperature is high in this region, which enhances convective heat transfer.Variations in solar cell temperature obtained using different values of Re, when SD = 0.001,ε
=
0.95, and C
=
100.
(a)
z
0
=
2
mm
(b)
z
0
=
5
mm
(c)
z
0
=
8
mm
(d)
Uniform flux


### 4.2. Effect of a Nonuniform Light Profile on the Fluid Bulk Temperature

Figures8 and 9 show the variation of fluid bulk temperature with increasing axial distance from the channel entrance under different incident light profiles. The results presented in these two figures demonstrate that the incident light distribution strongly influences the fluid bulk temperature profile inside the channel. However, the effect of the incident light profile on the fluid bulk temperature at the channel exit is insignificant. Figures 10(a), 10(b), and 10(c) show the effect of the emissivity of the two duct walls on the fluid bulk temperature when z
0
=
2, 5, and 8 mm, respectively. Surface radiation increased the temperature of the fluid closer to the entrance than pure convection (i.e., ε
=
0) due to radiative heat exchanges between the heated and the adiabatic walls, which redistribute the temperature profiles of the walls and hence the convective heat transferred to the flowing fluid. The effect of emissivity on the fluid bulk temperature at the channel exit is insignificant, as shown in Figures 10(a)–10(c).Figure 8
Variations in fluid bulk temperature obtained using different values of SD, whenz
m
=
0.5, Re = 500, ε
=
0.95, and C
=
100.Figure 9
Variations in fluid bulk temperature obtained using different values ofz
m, when SD = 0.001, Re = 500, ε
=
0.95, and C
=
100.Variations in fluid bulk temperature obtained using different values ofε, when SD = 0.001, Re = 500, and C
=
100.
(a)
z
0
=
2
mm
(b)
z
0
=
5
mm
(c)
z
0
=
8
mmThe effects of Re on the fluid bulk temperature are shown in Figures11(a)–11(d) for uniform and nonuniform light flux. The results demonstrate that the fluid bulk temperature decreases with increasing Re. In this example, the fluid bulk temperature decreases by approximately 43% as Re increases by 300% for given values of controlling parameters. However, the collected thermal energy increases with increasing Re, as a result of increasing flow rate.Variations in fluid bulk temperature obtained using different values of Re, when SD = 0.001,ε
=
0.95, and C
=
100.
(a)
z
0
=
2
mm
(b)
z
0
=
5
mm
(c)
z
0
=
8
mm
(d)
Uniform flux


## 4.1. Effect of a Nonuniform Light Profile on the Solar Cell Temperature

Figure4 shows the variations of solar cell temperature with z
0
=
5
mm and with the axial distance from the channel entrance for both uniform and nonuniform incident light at SD = 0.0005, 0.001, and 0.002. For a uniform distribution, the results demonstrate that the solar cell temperature first increases sharply near the entrance and then gradually increases until reaching its maximum value at the channel exit. In addition, there is a sudden increase in cell temperature at z
=
4
mm and a sudden decrease at z
=
6
mm. These changes occur due to the fact that no current is generated in this portion of the solar cell and all of the light is converted to heat due to the presence of the Cu-Ag-Hg front contact. However, this increment is small (≈7°C) and has minor effects on the overall efficiency. For a nonuniform distribution, the temperature of the first quarter of the solar cell remains constant at ambient temperature, then rapidly increases, and reaches its maximum value at the center of the cell. The temperature then drops sharply and gradually in the fourth quarter of the solar cell where the temperatures are higher than those observed in the first quarter. The profile of the temperature becomes thinner and the maximum temperature increases as the SD decreases, as shown in Figure 4. In addition, the maximum temperature of the solar cell is much higher when using a nonuniform light distribution than when using a uniform light distribution. In this example, the maximum temperature obtained using a uniform light distribution was 240°C, which indicates that the solar cell can be operated with a reasonable efficiency as the modern multijunction solar cells can adequately be operated at this level of temperature [25, 26]. In contrast, the results obtained using a nonuniform light distribution indicate that the cell efficiency will be drastically reduced due to the high temperature of some portions of the cell surface.Figure 4
Variations in solar cell temperature obtained using different values of SD, whenz
0
=
0.5, Re
=
500, ε
=
0.95, and C
=
100.The effects of the position of maximum incident light on the solar cell temperature distribution are shown in Figure5. A hot spot formed on the portion of the cell surface that was exposed to a large amount of light. This hot spot is slightly hotter when z
0
=
5
mm (i.e., at the center of the solar cell) due to the presence of the front contact metal in this region. Figures 6(a)–6(d) show the effect of the emissivity of the two duct walls on the temperature profile of the solar cell for both uniform and nonuniform incident light when z
0
=
2, 5, and 8 mm and when ε
=
0 (i.e., pure forced convection), 0.5, and 0.95. The results presented in these figures demonstrate that the maximum solar cell temperature decreases significantly as surface emissivity increases. However, the presence of surface radiation inside the channel does not change the shape of the temperature distribution for the nonuniform light profile. In this example, as emissivity increased from ε
=
0 to 0.95, the maximum solar cell temperature was reduced by approximately 35% when using nonuniform illumination and 20% when using uniform illumination. The large contribution of surface radiation to the nonuniform light profile is due to the increase in the effective temperature differential for surface radiation between the two duct walls. Nevertheless, the maximum temperature of the solar cell when using nonuniform light was still higher than that obtained using a uniform light distribution. For the uniform light profile, surface radiation exchange between the two duct walls and between the heated duct wall and the inlet and exit areas of the channel leads to a reduction of maximum cell temperature from 305°C to 240°C, which enables the solar cell to operate with a reasonable efficiency.Figure 5
Variations in solar cell temperature obtained using different values ofz
0, when SD = 0.001, Re = 500, ε
=
0.95, and C
=
100.Variations in solar cell temperature obtained using different values ofε, when SD = 0.001, Re = 500, and C
=
100.
(a)
z
0
=
2
mm
(b)
z
0
=
5
mm
(c)
z
0
=
8
mm
(d)
Uniform fluxThe results presented in Figures7(a)–7(c) illustrate the variations in solar cell temperature that occur with changes in axial distance from the channel entrance when using uniform and nonuniform light distributions, three selected values of Reynolds number (500, 1000, and 2000), and a given value of ε
=
0.95. The cell temperature decreases with increasing Re, but the profile maintains its general shape. The percent reductions of maximum solar cell temperature that occurred as Reynolds number increased by 300% from 500 to 2000 were 20.7%, 13.4%, and 13.7% for z
0
=
2, 5, and 8 mm, respectively. This value was approximately 22.5% when a uniform light distribution was employed. These results indicate that, for nonuniform light profiles, increasing flow rate is more useful when the maximum incident light is located near the channel entrance because the difference between the wall temperature and the fluid temperature is high in this region, which enhances convective heat transfer.Variations in solar cell temperature obtained using different values of Re, when SD = 0.001,ε
=
0.95, and C
=
100.
(a)
z
0
=
2
mm
(b)
z
0
=
5
mm
(c)
z
0
=
8
mm
(d)
Uniform flux


## 4.2. Effect of a Nonuniform Light Profile on the Fluid Bulk Temperature

Figures8 and 9 show the variation of fluid bulk temperature with increasing axial distance from the channel entrance under different incident light profiles. The results presented in these two figures demonstrate that the incident light distribution strongly influences the fluid bulk temperature profile inside the channel. However, the effect of the incident light profile on the fluid bulk temperature at the channel exit is insignificant. Figures 10(a), 10(b), and 10(c) show the effect of the emissivity of the two duct walls on the fluid bulk temperature when z
0
=
2, 5, and 8 mm, respectively. Surface radiation increased the temperature of the fluid closer to the entrance than pure convection (i.e., ε
=
0) due to radiative heat exchanges between the heated and the adiabatic walls, which redistribute the temperature profiles of the walls and hence the convective heat transferred to the flowing fluid. The effect of emissivity on the fluid bulk temperature at the channel exit is insignificant, as shown in Figures 10(a)–10(c).Figure 8
Variations in fluid bulk temperature obtained using different values of SD, whenz
m
=
0.5, Re = 500, ε
=
0.95, and C
=
100.Figure 9
Variations in fluid bulk temperature obtained using different values ofz
m, when SD = 0.001, Re = 500, ε
=
0.95, and C
=
100.Variations in fluid bulk temperature obtained using different values ofε, when SD = 0.001, Re = 500, and C
=
100.
(a)
z
0
=
2
mm
(b)
z
0
=
5
mm
(c)
z
0
=
8
mmThe effects of Re on the fluid bulk temperature are shown in Figures11(a)–11(d) for uniform and nonuniform light flux. The results demonstrate that the fluid bulk temperature decreases with increasing Re. In this example, the fluid bulk temperature decreases by approximately 43% as Re increases by 300% for given values of controlling parameters. However, the collected thermal energy increases with increasing Re, as a result of increasing flow rate.Variations in fluid bulk temperature obtained using different values of Re, when SD = 0.001,ε
=
0.95, and C
=
100.
(a)
z
0
=
2
mm
(b)
z
0
=
5
mm
(c)
z
0
=
8
mm
(d)
Uniform flux


## 5. Conclusions

The effects of nonuniformity of the incident light on the temperature and fluid bulk temperature of a triple junction solar cell were studied. Nonuniformity was simulated using a Gaussian distribution. The results indicated that nonuniform illumination significantly alters the temperature distribution of the solar cell, causing hot spots in some regions of the solar cell. Although increasing the emissivity of the two duct walls and increasing Re can reduce the maximum solar cell temperature by as much as 35% and 20%, respectively, solar cell temperatures in the region exposed to the most incident light are high enough to cause degradation in solar cell efficiency at medium concentrations. Thus, secondary optical concentrators are required to redistribute the incident light on the cell surface into a uniform shape. Alternatively, an innovative cooling system, such as the introduction of a secondary jet flow into a channel, should be used to redistribute solar cell temperatures and reduce the maximum temperature to an acceptable value. Surface emissivity and nonuniform illumination affect the fluid bulk temperature distribution inside the channel. However, the effect of these parameters on fluid bulk temperature at the channel exit is insignificant.



---

*Source: 642819-2014-06-16.xml*