The geometrical optics model about the offset effect of solar rays by the thickness of concentrating mirror and the diametric solar model were established. The radiant flux density on the surface of the receiver for parabolic trough solar concentrator was obtained by numerical calculation with the established models. Charge-coupled device (CCD) was used for testing gray image on the surface of the receiver for parabolic trough solar concentrator. The image was analyzed by Matlab and the radiant flux density on the surface of the receiver for parabolic trough solar concentrator was achieved. It was found that the result of the theory is consistent with that of the experiment, and the relative deviation on the focal length width was 8.7%. The geometrical structure of receiver based on parabolic trough solar concentrator was optimized, a new parabolic receiver has been proposed, and it has been shown that the optimized geometrical structure of receiver was beneficial to improve the working performance of the entire system.
1. Introduction
The environmental issue is a significant challenge for humanity; therefore many countries are vigorously developing renewable energy technology in the world at present [1]. Solar energy has attracted much attention over the past decades with its renewable, clean, and abundant characteristics, so the concentrating solar power system (CSP) is getting more attention. Concentrating solar power system mainly consists of parabolic trough solar concentrator, linear Fresnel solar concentrator, compound parabolic solar concentrator, and solar tower, among which, the parabolic trough solar concentrator, as the most mature technology, is most widely researched and its cost is gradually reduced day by day [2–6]. At present, the parabolic trough solar concentrator, in which operating temperature is more than 400°C, has been put into commercial operation in many countries [7] such as America, Germany, Spain, and China.
Parabolic trough solar concentrating system is mainly composed of concentrating device and thermal converting equipment. Concentrating device is the key to the entire system, and its characteristic, which was investigated by many researches, is very important for the performance of system [8]. The optical performance of half parabolic trough solar concentrator has been studied [9], and the results showed that the width f focal line for half parabolic trough solar concentrator was associated with the rim angle and aperture of system. Paper [10] has studied the thermal and optical performance in the parabolic trough solar concentrator, and the results showed that the stress of concentrating mirror was large at the edge of system and the deformation of concentrating mirror was bigger than that of anywhere in the whole system, so the optical performance at the edge was not well and the system efficiency was reduced [10]. Study [11] demonstrated that the optical loss of concentrating mirror and thermal converting equipment accounted for about 20% of the total energy losses in parabolic trough solar concentrator. It was found that optical concentrating ratio was small influence by focal length in [12]. The method of measuring the error of surface shape with concentrating mirror was used to study the concentrating performance of parabolic trough solar concentrator [13]; it has been carried out that the performance of concentrated solar ray was reduced with the mirror surface shape as tiny change.
In the real situation, the solar rays arriving at the parabolic trough solar concentrator are partly reflected by the concentrating mirror, and the left get into concentrating mirror. However, there is rare investigation on this topic. In this study, the geometrical optics model about the offset effect of solar rays by the thickness of concentrating mirror was established. The diametric solar model was established for the convenience of calculation. The radiant flux density on the focal plane for parabolic trough solar concentrator was calculated by the established models, and the numerical result was verified by the experiment. Based on the photovoltaic and thermal system with the parabolic trough solar concentrator, the geometrical structure of receiver was optimized, and it was found that the optimized geometrical structure of receiver was beneficial to improve the working performance for whole system.
2. The Optical Model of Concentrator Mirror2.1. Optical Characteristic of Concentrating Mirror
The solar rays arrive at the concentrating mirror, which are reflected, refracted, and transmitted by the concentrating mirror under the regulation of tracking device. In the actual parabolic trough solar concentrator, in order to improve the performance of the system, the material of concentrating mirror is ultraclear glass, and the thickness of concentrating mirror is not more than 5 mm.
Next, we calculate the optic performance of the ultraclear glass. Assuming that the thickness L of concentrating mirror is 1 mm, 2 mm, 3 mm, 4 mm, and 5 mm, respectively, the refractive index of the concentrating mirror is 1.6, and the extinction coefficient K of the concentrating mirror is 4.0 m−1 [14]. When the solar rays reach the 5 types of ultraclear glass vertically, the computational formula of absorbance, reflectance, and transmittance is as follows [15]:
(1)r=n-1n+12,α=1-e-KL1-r1-re-KL,τ=e-KL1-r21-re-KL2,ρ=r1+e-KLτ,
where r is surface reflectance of glass, α is absorbance, τ is transmittance, and ρ is reflectance. The calculated results are shown in Table 1.
The absorbance, reflectance, and transmittance of ultraclear glass at different thicknesses.
Thickness of glass
α
τ
ρ
α+τ+ρ
1 mm
0.003991
0.895186
0.100823
1.00000
2 mm
0.007965
0.891592
0.100443
1.00000
3 mm
0.011920
0.888013
0.100067
1.00000
4 mm
0.015859
0.884449
0.099693
1.00000
5 mm
0.019779
0.880899
0.099322
1.00000
It can be seen in Table 1 that the absorption is very small for any thickness of glass. In general, the reflecting layer of ultraclear glass is silver, and its reflectivity is more than 95% [16]. So the main path of sunlight for parabolic trough solar concentrator is shown in Figure 1 when the solar rays arrive at the concentrating mirror. In Figure 1, WN and VM are first reflecting surface and second silver reflecting surface, respectively. SA is the incoming solar ray, AR is reflected solar ray, EG is the refracted solar ray, AC and CE are inner solar rays of ultraclear glass, and VW is the thickness of ultraclear glass.
The main path of sunlight arriving at the concentrating mirror.
2.2. The Equation of Sunlight Path
The concentrating mirror of parabolic trough solar concentrator is a curved cylindrical mirror, which is curved by plane mirror usually. As shown in Figure 1, the plane mirror is pure bending deformation, so the top and bottom surfaces of curved plane mirror are equidistant surface by mechanics of materials. During the manufacturing process of concentrating mirror, the bottom surface of concentrating mirror is parabolic cylinder generally. The top surface is equidistant surface of bottom surface, and the distance between top surface and bottom surface of concentrating mirror is the thickness of the plane mirror.
Assuming that the parametric equation of surface M in Figure 1 is
(2)x=4ft,y=4ft2,
where t is parameter and f is focal length, based on the hypothesis of the equidistance surface, the parametric equation of surface N is
(3)x=4ft-2δt1+4t2,y=4ft2+δ1+4t2;δ is the thickness of concentrating mirror. The refracted ray AC can be gained by Snell’s law:
(4)y-4ft2+δ1+4t2W=4t2+u2t-2tux-4ft+2δt/1+4t2,(5)u=4n2t2+n2-4t2.
Combining formulas (2) and (4), we have
(6)xCt=-2f4t2+u2tu+A1+A2A3,A1=2f4t2+u2tu-2t2,A2=4f4ft2+δ1+4t2+4t2+u2t-2tu,A3=2δt1+4t2-4ft,kCD=-2fxCt,
where xCt is the parametric abscissa of point C. k is slope. The ray CE can be expressed as follows:
(7)y-yCt=kCEx-xCt,(8)kCE=kAC-2kCD-k2CDkACk2CD-2kCDkAC-1,(9)kAC=4t2+u2t-2tu.
From formula (3), the general equation of surface N is
(10)x=4f-2δA4A42-12,A4=231+y/f3cosα3,α=arccos-3q-3p2p2,p=-1+yf,q=δf.
So yE can be written as follows:
(11)yE-yC=kCE4f-2δA4A42-12-xC.
Formula (11) is a transcendental equation which can be programmed for solving. So the refracted ray EG can be obtained as follows:
(12)y-yE=kEGx-xE,(13)kEG=tanθ9+kEF1-kEFtanθ9,(14)kEF=-12tE,(15)tanθ9=nsinθ81-n2sin2θ8,(16)sinθ8=tanθ81+tanθ82,(17)tanθ8=kCE-kEF1+kCEkEF.
From Figure 1, the refracted ray AR can be expressed as follows:
(18)y-4ft2+δ1+4t2=4t2-14tx-4ft+2δt1+4t2.y is replaced by f in formulas (12) and (18), the horizontal ordinate x can be calculated when solar rays arrive on focal plane of system. The solar subtending angle is not taken into consideration and solar ray is assumed to come from the centre of the sun in formulas (12) and (18). When the solar ray comes from the other position, the θ1 will be replaced by θz, and the value of θz will be discussed below.
3. Diametric Solar Model3.1. The Numerical Calculation of Buie Model
There are many optical models of the sun, among which, Buie model is widely used. The model was established based on the observer standing on the earth and the face of the observer towards the centre of the sun. In the model, the luminance value of the sun’s centre is set to be 1, and the luminance value of other positions on the sun is the relative luminance value of center position. The advantage of the model is not affected by geographical location of the observer; the model is axial symmetry and center symmetry, and it can be written as [17–20]
(19)Φθs=cos0.326θscos0.308θs,0≤θs≤4.65mrad,ekθsγ,θS>4.65mrad,k=0.9ln13.5χχ-0.3,γ=2.2ln0.52χχ0.43-0.1,
where Φ(θs) is factor of luminance value. θs is the polar angle of the sun; the value of θs is 16′ at the edge of the sun. χ is circumsolar ratio (CSR) [18]; CSR is defined as the radiant flux contained within the circumsolar region of the sky divided by the incident radiant flux from the direct beam. Under the sunny weather condition, χ is 0.05, and the result of numerical calculation for Buie model is shown in Figure 2.
The calculated result of Buie model.
It can be seen from Figure 2 that luminance value of the sun is discontinuous in Buie model, and the luminance value of inside the sun is greater than that of outside the sun, so the outside energy of the sun is negligible in the concentrating solar system, the polar angle of the sun just only 16′ at the edge of the sun [21]; therefore, the range of polar angle θz can be expressed as follows:
(20)θ1-16′≤θz≤θ1+16′.
3.2. Optical Performance between Receiver and Concentrating Mirror
There is a straight line on the sun’s surface, which is parallel to the generatrix of parabolic trough solar concentrator, and the rays of the straight line arrive to the concentrating mirror, the main path of sunlight arriving at the concentrating mirror as shown in Figure 3. Figure 3 shows the path of two rays, in green and pink, which are reflected and refracted; CB is parallel to the generatrix of parabolic trough solar concentrator. AJ is the normal of concentrating mirror at point D. Based on the law of reflection, we have that AB is equal to AS and AC is equal to AR. It can be obtained that RS is parallel to CB. So the solar rays of every point, in the arbitrary straight line which is parallel to the generatrix of parabolic trough solar concentrator on the sun’s surface, arrive to the surface of concentrating mirror. Some of them are reflected directly by concentrating mirror to the receiver; these solar rays form a straight line at surface of receiver, and the straight line is parallel to the arbitrary straight line on the sun’s surface.
The main sunlight path for parabolic trough solar concentrator.
From Figure 3, KH is normal of reflection ray, Plane ABD is perpendicular to plane BCD, AB is perpendicular to BC and AD, respectively, BD and CD are incident rays, DG and DH are refracted rays, G,H are points of intersection between DG, DH and bottom surface of concentrating mirror. Assuming that GH is parallel to BC, so it can be obtained that
(21)BDcos∠ADB=DCcos∠ADC,(22)DFcos∠FDI=DEcos∠EDI,(23)BDDC=ADcos∠ADBcos∠ADCAD=cos∠ADCcos∠ADB,(24)DFEC=DIcos∠FDIcos∠EDIDI=cos∠EDIcos∠FDI,(25)DItan∠FDIDItan∠HDI=DItan∠EDIDItan∠GDI,(26)tan∠ADBtan∠HDI=tan∠ADCtan∠GDI.
Obviously, formula (26) is contrary to Snell’s law. So the assumption is not right, GH is not parallel to BC. But the value of tangent function is proximately equal to sine function when the angle value is less than 5°. The sun subtends an angle of 32′, so it is high order approximate that GH is parallel to BC. We have the conclusion that the solar rays are refracted by concentrating mirror to the receiver; these solar rays form a straight line at surface of receiver, and the straight line is highly approximatively to parallel to the arbitrary straight line on the sun’s surface.
3.3. The Diametric Solar Model
The geometric relationship of Buie model is shown in Figure 4, so
(27)θS=arctanx2+y2D,0≤y≤R2-x2,
where D is average distance between the sun and the earth, and it is 1.495 × 1011 m. R is radius of the sun, and it is 6.950 × 108 m [21].
The geometric relationship of Buie model.
Based on the optical performance of solar rays reflected and refracted by concentrating mirror, Buie model can be simplified. Buie model could be compressed with a diameter which is parallel to the generatrix of parabolic trough solar concentrator as shown in Figure 5.
The diametral solar model.
So the luminance factor I(x) of diametric solar model is
(28)Ix=∫0y0xΦx,ydyU,U=∫0Rcos0.326arctanr/Dcos0.308arctanr/Ddr.
Based on the geometric relationship from Figure 5, we have
(29)y=xtanα,dy=xsec2αdα,α0=arccosxR.
So
(30)Ix=∫0arccosx/RL1xsec2α/L2dαU,(31)L1=cos0.326arctanxsecαD,(32)L2=cos0.308arctanxsecαD.
The primitive function of formula (30) is hard to find. We can use high order algebraic accuracy of the Gauss formula to calculate as follows [22]:
(33)∫0bfxdx≈b2fb2+b23+fb2-b23.
4. The Numerical Method of Diametric Solar Model
Based on the two-dimensional parabolic trough solar concentrator, Fortran language is used to calculate the radiant flux density on the focal plane as shown in Figure 6. In the programming, the arc length of parabolic trough solar concentrator and the solar subtending angle are scattered, and the luminance factor could be obtained from the diametric solar model. The position and energy of the solar rays arrived on the surface of receiver could be calculated with formulas (12) and (18) and the radiant flux density could be gained by frequency statistics method.
Schematic of the numerical method for energy density distribution.
The flow chart of the program is shown in Figure 7, step size of concentrating mirror and solar subtending angle are 0.5 mm and 0.2′, respectively. It can be seen in Figure 1 that yE is less thanyA with geometric relationship, therefore the initial value of yE could take the place of yA, and then the value of yE decreased time by time with step size until calculation error to meet the requirement.
The flow chart of the numerical process.
5. Experimental Verification
The calculated results were as verified by experiment. CCD was used for testing concentrating characteristic of parabolic trough solar concentrator, and the schematic is shown in Figure 8. The CCD is put on the top surface of parabolic trough solar concentrator. CCD and computer are connected by data line. Camera lens of CCD is toward the Lambert target to make the grayscale image on the surface of Lambert target clearly visible on the computer. In order to prevent saturation of the gray images taken by CCD, a neutral attenuation was put in front of CCD lens. The specifications of the CCD are shown in Table 2.
The specifications of the CCD.
Parameters
Specification
Unit
Array size
1280×1024
Pixel
Pixel size
5.2
μm
Signal to noise ratio
45
dB
Sensitivity
2.1
V/Lux-s
Time of exposure
50
μs
The testing schematic of CCD.
Figure 9 is experimental system of parabolic trough solar concentrator, the aperture width of the system is 161 cm, and focal length is 106 cm. Concentrating mirror is ultraclear glass with a thickness of 2 mm. The reflectivity of the concentrating mirror is 93.7%, which was measured by spectrophotometer. The refractive index of concentrating mirror is 1.57.
The parabolic trough solar concentrator of measuring energy density distribution.
In the experiment, a two-dimensional tracking device was used in parabolic trough solar concentrator. There was a pyrheliometer near the parabolic trough solar concentrator to measure solar irradiance. The solar irradiance is measured by pyrheliometer which can track the sun automatically. The relative error of the pyrheliometer is no more than 3%.
The camera was used to take the width of focal line on the surface of receiver (14 September, 2013, 12:51, Kunming, China) as shown in Figure 10, and the scale division of dividing ruler is millimeter. It can be seen in Figure 10 that the focal line width of experiment system is about 2 cm.
The width of focal line for parabolic trough solar concentrator.
Based on the experimental system, the concentrating equation of concentrating mirror and the diametric solar model were used to calculate the luminance factor of focal plane as shown in Figure 11; in order to prevent the influence of the direct irradiance, the centre luminance factor of focal plane is 1, and the other position of the focal plane is relative luminance factor.
Calculated result of the luminance factor of focal plane.
In order to get the radiant flux density of focal plane at different direct irradiance, we can use formula as follows:
(34)DpId=2∫0x0Icfxdx,Ix=fxIc,
where Dp is the aperture width of system, Id is direct irradiance, Ic is irradiance constant with Dp and Id, I(x) is focal plane irradiance at the position of x, f(x) is scale factor of relative luminance, and x0 is the integral upper limit as shown in Figure 11.
The grayscale images of the radiation with focal plane of system and the sky were shot by CCD as shown in Figure 12; the direct irradiance of the sky was 829 W/m2 measured by pyrheliometer (14 September, 2013, 12:52, Kunming, China), and Matlab was used for analyzing the gray value of the images. The average gray value of the image of direct irradiance on the sky is 3.08, and the value of irradiance is proportional to the gray level, so we have
(35)Ix=829Hx3.08,
where H(x) is the gray level of image.
The grayscale image taken by CCD.
The grayscale image of the focal plane of system
The grayscale image of the direct irradiance
The results of radiant flux density on the focal plane of parabolic trough solar concentrator by theoretical calculation and experimental measurement are shown in Figure 13. It can be seen in Figure 13 that the widths of focal line by calculation and experiment are 18.04 mm and 19.76 mm, respectively, and the relative deviation is 8.7%. The centre irradiance on focal plane of experimental result is smaller than that of numerical result. This is due to the fact that the experimental system has some machining errors, so the concentrating mirror is not standard parabolic mirror. The tracking device has certain hysteresis, and it is difficult to track the sun continuously. From Figure 11, it can be seen that the endpoint in the irradiance is zero by theoretical calculation, but the endpoint in the irradiance is not zero in the experiments. This is due to the reason that the irradiance on the surface of receiver was calculated by theory that just only considered the solar rays coming from the sun, but in the practical situation of operating parabolic trough solar concentrator, a few part of the irradiance on the surface of receiver comes from the diffuse radiation.
Radiant flux density of calculation and experiment.
6. Optimization of the Receiver
In the solar concentrating system, the radiant flux density on the surface of concentrating solar cell is more uniform, and the operating status of PN junction with solar cell is more consistent, so the loss of circulation current in solar cell is lower, the operating temperature of solar cell is lower, the temperature rise of cooling medium is much less, and the fill factor and output efficiency are much higher. It is necessary to optimize radiant flux density characteristics on the surface of the receiver. During the operating process of parabolic trough solar concentrator, the solar rays are concentrated by concentrating mirror; if the position of focal plane has receiver, the radiant energy will be converted into electricity or thermal energy; otherwise the solar rays will diverge. By geometric optical properties of parabola and the basic principle of reversible optical path, we can adopt appropriate ways to improve the uniformity of the radiant flux density on the surface of receiver; the receiver surface can be designed as parabolic and the focal point of parabolic receiver and parabolic trough solar concentrator is set to the same point as shown in Figure 14.
The parabola absorber is adopted on parabolic trough solar concentrator.
The radiant flux density on the surface of parabolic receiver and flat receiver is shown in Figure 15; the direct radiation is 1000 W/m2, the rim angle of parabolic trough solar concentrator is 45°, the focal length is 1 m, the thickness of concentrating mirror is 2 mm, and reflectivity and refractivity of ultraclear glass are 95% and 1.6, respectively.
The radiant flux density on the surface of parabolic receiver and flat receiver.
It can be seen in Figure 15 that the centre energy density on the focal plane of parabolic receiver is smaller than that of plane receiver, and the variances of the radiant flux density on the focal plane of parabolic receiver and plane receiver are 7.2 × 108 and 5.3 × 108, respectively. The widths of focal line for parabolic receiver and plane receiver are 21.42 mm and 17.92 mm, respectively. It is beneficial for operating performance of solar concentrating system.
7. Conclusion
The optical model of concentrator mirror and the diametric solar model were established. Based on the parabolic trough solar concentrator, the geometric construction of receiver was optimized. The main conclusion can be obtained as follows.
The optical model of concentrator mirror and the diametric solar model were verified by experiment and it was found that results of theory and experiment were consistent with each other. The models could be supplied to the design for parabolic trough solar concentrator.
Based on the geometric optical properties of parabola and the principle of reversible optical path, a new parabolic receiver has been proposed. It was found that the optimized and designed parabolic receiver was beneficial to improve the characteristic of concentrating solar cells.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The present study was supported by National Natural Science Foundation, China (Grant no. U1137605), the Program of Changjiang Scholars and Innovative Research Team in Ministry of Education, China (Grant no. IRT0979), and the National and international scientific and technological Cooperation projects, China (Grant no. 2011DFA60460).
WeiJ.HuangB.WangP.Photocatalytic properties of nitrogen-doped Bi12TiO20 synthesized by urea addition sol-gel method20122012813513210.1155/2012/1351322-s2.0-84863073986WangF.-Q.ShuaiY.TanH.-P.Thermal analysis of solar cavity type receiver20113258438462-s2.0-79959715841SolankiC. S.SanganiC. S.GunashekarD.AntonyG.Enhanced heat dissipation of V-trough PV modules for better performance200892121634163810.1016/j.solmat.2008.07.0222-s2.0-55349141055TaoT.HongfeiZ.KaiyanH.MayereA.A new trough solar concentrator and its performance analysis201185119820710.1016/j.solener.2010.08.0172-s2.0-78650307248ZhaoH.LiuH.JingL.LiuY.LiC.LuZ.SunQ.Discrete spectral local measurement method for testing solar concentrators20122012829047410.1155/2012/2904742-s2.0-84863037313Fernández-GarcíaA.ZarzaE.ValenzuelaL.PérezM.Parabolic-trough solar collectors and their applications20101471695172110.1016/j.rser.2010.03.0122-s2.0-77956093649HeY.-L.XiaoJ.ChengZ.-D.TaoY.-B.A MCRT and FVM coupled simulation method for energy conversion process in parabolic trough solar collector201136397698510.1016/j.renene.2010.07.0172-s2.0-77958181556GongG.HuangX.WangJ.HaoM.An optimized model and test of the China's first high temperature parabolic trough solar receiver201084122230224510.1016/j.solener.2010.08.0032-s2.0-78349312145JiangS. L.2009Heifei, ChinaUniversity of Science and Technoligy of ChinaCoventryJ.2004Canberra, AustraliaAustralian National UniversityCuiY. H.BeiZ. H.ZhaoX.Research on the heat loss of parabolic trough solar collector field201028559ShuaiY.ZhangX.-F.TanH.-P.Simulation for concentrating characteristic of parabolic solar collector system20062734844862-s2.0-33744952712MaccariA.MontecchiM.An optical profilometer for the characterisation of parabolic trough solar concentrators200781218519410.1016/j.solener.2006.04.0042-s2.0-33846183600ZhangH. F.2007Xi'an, ChinaNorthwestern Polytechnical University PressDuffieA.BeckmanA.1980Madison, Wis, USWiley-InterscienceDuJ.-L.TangD.-W.LiT.SuG.-P.Design and experiment research of a solar simulator in Dish/Stirling Solar Power Generation System20103111188318852-s2.0-78649721774XuC. M.LiM.JiX.LuoX.Frequency statistics analysis for energy-flux-density distribution on focal plane of parabolic trough solar concentrators2013334040800110.3788/AOS201333.04080012-s2.0-84877886943BuieD.MongerA. G.DeyC. J.Sunshape distributions for terrestrial solar simulations200374211312210.1016/S0038-092X(03)00125-72-s2.0-0037742417LiuY.DaiJ.LangZ.XinC.Finite-element analysis for flux distribution on focal plane of rotating parabolic concentrators20072710177517782-s2.0-35948944903DaiJ. M.LiuY.The study of flux distribution on focal plane in parabolic trough concentrators2008299109611002-s2.0-55149117985DuffieJ. A.BeckmanW. A.20134thHoboken, NJ, USAJohn Wiley & SonsJiangJ. F.HuL. J.TangJ.2004Beijing, ChinaScience Press