# Effect of Tracking Error of Double-Axis Tracking Device on the Optical Performance of Solar Dish Concentrator.

1. IntroductionSolar energy is a clean and environment-friendly renewable energy source, which is plentiful and can be widely distributed. Developing and utilizing solar energy are important ways to solve energy shortages and environmental pollution problems, for example, the H2 production by solar thermochemical reactions [1], solar cooker [2], and solar energy thermal electric power system [3]. Among all those typical solar energy applications in medium-high temperature fields, the solar concentrator was commonly used in the abovementioned fields to collect the sunlight as heat source and provide heat source to the heat receiver because the solar radiation density received by the earth surface was low. There are four main types of solar concentrators: Fresnel concentrator [4,5], parabolic trough concentrator [6,7], tower heliostat [8, 9], and parabolic dish concentrator [10-12]. The single-dish concentrator could realize the high concentration ratio and high temperature so it is superior to other concentrators. Typical application, such as the solar dish concentrator photovoltaic system and the dish-Stirling solar power system, consists of a dish concentrator and a Stirling engine [10] such as the 38 kW XEM-Dish system as shown in Figure 1.

The dish concentrator via an elevation-azimuth double-axis tracking device (EATD) was used to accurately track the sun's position so it could focus the sunlight to the precalculated position of the receiver [10-12]. The EATD consists of an elevation tracking device and an azimuth tracking device (Figures 1(b) and 1(c)), which is widely used in dish-concentrating system. The azimuth tracking axis in EATD is vertical to the ground level, and the elevation tracking axis in EATD is parallel to the ground level and vertical to the focal axis of a dish concentrator. The program tracking process of a dish concentrator system is as follows: calculate the position information of the sun (the elevation angle and azimuth angle of the sun) by astronomical algorithm [13]; according to the sun position information, the output rotation angle of an azimuth motor and elevation motor in the EATD is controlled and then the elevation tracking axis [n.sub.1] and azimuth tracking axis [n.sub.2] in the EATD can rotate a corresponding angle (Figure 2(b)), which will change the azimuth and elevation angle of the dish concentrator for realizing the focal axis aiming at the sun's center. Therefore, the tracking accuracy of a dish concentrator could be influenced by astronomical algorithm, driving error of the motor and EATD, installation error of the EATD, and so on, where driving error of the EATD is the main factor (i.e., output rotation angle error of the azimuth and elevation tracking axis in the EATD, called EATD tracking error). The EATD tracking error could lead the dish concentrator to produce tracking error so it would cause some adverse phenomenon, such as the optical losses and nonuniform flux in the receiver. Thus, improving the tracking accuracy of a dish concentrator and reducing the manufacturing precision requirement of the EATD (i.e., reducing the manufacturing cost) have always been the pursuit of commercial applications. Therefore, researching the influence of the EATD tracking error on the tracking performance and optical performance of dish concentrator system is very important.

Many works focus on the research for the influence of tracking error on the optical performance or flux distribution of dish concentrator system [14-21], but they do not consider the influence of EATD tracking error. Hughes [14] used statistical methods to establish an expression for the expected value of the intercept factor for various configurations and control law implementations. The analysis assumes that a radially symmetric flux distribution is generated at the focal plane. Badescu [15] researched the different tracking error distributions and their effects on the long-term performances of parabolic dish solar power systems; he adopted four different tracking error distributions to evaluate the intercept factor and the thermal and engine efficiencies, three of them are bounded distributions and the last one is the widely used classical Gaussian distribution. Further, he also studied the theoretical derivation of heliostat tracking error distribution using methods of integral geometry and geometrical probabilities [16]. Hu and Yao used MCRT method to research the flux distribution on the solar cell plane in a two-stage concentration photovoltaic system with concentrator tracking error (dish concentrator is the first level) [17]. Shanks et al. [18] researched the flux distribution and optical efficiency in a high-concentrating photovoltaic system with concentrator tracking error (using the Cassegrain concentrator). Xia et al. [19] researched the influence of the concentrator tracking error on the intercept factor of a sixteen-dish concentrator. Li et al. [20] and Castellanos et al. [21] through the theoretical model research the influence of the total optical error on the intercept factor in the dish concentrator system. In the experimental measurement, the measured results of the flux distribution on the focal plane have been widely used to evaluate the surface slope error, tracking error, and optical efficiency of the dish concentrator system [19, 22, 23], which could also be used for adjusting the tracking error of dish concentrators. For example, Johnston [22] developed a flux mapping equipment and measured the flux distribution on the focal plane of a 20 [m.sup.2] dish concentrator. This concentrator mirror slope error has been determined is 2.0 mrad. Xia et al. [19] developed a type of concentrated solar flux distribution measurement formed by a thermal infrared imager and water-cooled Lambert target and measured the flux distribution on the focal plane of a sixteen-dish concentrator. Through the comparison of the measured results and ray tracking calculation results of the flux distribution, they determined the concentrator mirror slope error is 2.2 mrad and concentrator total tracking error is 0.8 mrad. But the single tracking error of the concentrator has not been analyzed in their paper. The measured results of the flux distribution on the focal plane of the dish concentrator could be used to evaluate or calibrate the tracking error of the dish concentrator. But this process would be influenced by the concentrator mirror slope error and solar direct normal irradiance (DNI), which would change the geometry shape or flux distribution of the focal spot on the focal plane [24, 25]. Therefore, it is necessary to find a new focal spot characteristic to establish the functional relationship between the focal spot characteristic and tracking error of dish concentrator to realize the evaluation and calibration of concentrator tracking error.

The above-mentioned literatures all focus on the influence of the total tracking error of dish concentrator systems on the flux distribution and optical performance, while the researches of the influence of the tracking error of the EATD on the tracking performance and optical performance of dish concentrator system are limited. In this paper, we focused on the above two aspects; a flux distribution model of the focal plane in dish concentrator systems has been established based on ray tracking method. The concentrator mirror and solar cone in nonuniform energy distribution are discreted by square grids in this model (in Section 2). This model was adopted for researching the influence of the mirror slope error and DNI and EATD tracking error on the focal spot characteristics (i.e., flux distribution, geometrical shape, centroid position, and intercept factor), and the tracking error transmission law of the EATD transferred to dish concentrator was studied (in Section 4). Moreover, the transmission law of the EATD azimuth tracking error in solar heliostats is analyzed, a dish concentrator system using a spin-elevation double-axis tracking device is proposed, and the transmission law of the tracking error from the spin axis to the dish concentrator was studied (in Section 4.4). This work could provide fundamental for calibrating the tracking error of dish concentrator system, manufacturing precision allocation of a double-axis tracking device (appropriately decrease the manufacturing precision of the azimuth or spin tracking device), and developing a new type of tracking device.

2. Mathematical Model

The flux distribution model of the plane receiver in dish concentrator systems is developed based on the ray tracking method. The ray tracking method needs to discrete the mirror surface of dish concentrators, plane receiver, and sun cone. There are two methods to discrete the mirror surface and sun cone, a random MCRT method [4, 12] and a grid discrete method [26]. In this work, we adopted the grid discrete method to establish the sun cone discretization model considering a nonuniform energy distribution and flux distribution model of the focal plane in dish concentrator systems (i.e., optical model). The mirror slope error of the dish concentrator and tracking error of the EATD are included in this optical model.

2.1. Optical Model. As shown in Figure 2, a coordinate system O-xyz was built on the vertex point O of a parabolic mirror surface and the z-axis points to the focal point F. The mirror surface function of a dish concentrator can be defined as [x.sup.2] + [y.sup.2] = 4fz, where f is the focal length of the dish concentrator. In Figure 2(a), [R.sub.1] is aperture radius of the dish concentrator, and angle [[theta].sub.m] and radius [R.sub.2] are no mirror area in the dish concentrator. The steps for optical discretization of the dish concentrator are as follows: (1) use a square with side length of 2[R.sub.1] (blue wire frame in Figure 2(a)) to surround the dish concentrator mirror and then the discrete into a lot of square grid with side length of 2[R.sub.1]/K (M = K, i.e., the daylighting area [A.sub.km] of each grid is equal, [A.sub.km] = [(2[R.sub.1]/K).sup.2]); (2) according to the coordinate and distance [R.sub.km] of the point [p.sub.km] (i.e., grid center point in k line m column) to determine whether the grid km is in the reflection mirror. If the point [p.sub.km] is in the reflection mirror, the identifier [f.sub.km] = 1 (i.e., the grid km is effective); otherwise, [f.sub.km] = 0 (i.e., the grid km is a noneffective); (3) the coordinates and unit normal vector [N.sub.km] of the center point [p.sub.km] of all effective grid are calculated and stored for subsequent ray tracing simulation.

In practical application, the reflection mirror of a dish concentrator would have microslope error, which would impact the normal vector of the reflection mirror (Figure 2(b)). The mirror slope error could be described as a radial angle [[epsilon].sub.3] and tangential angle [[epsilon].sub.4] [17,24]:

[[epsilon].sub.3] = [square root of (-2[[sigma].sup.2])ln [r.sub.1] [[epsilon].sub.4] = 2[pi][r.sub.2], (1)

where [r.sub.1] and [r.sub.2] are uniformly distributed random numbers within the range of 0~1, [sigma] is standard deviation of the mirror slope error.

Considering the mirror slope error of the dish concentrator, the vector [N.sup.1.sub.km] on point [p.sub.km] is calculated by (2), where Rot([n.sub.3], [[epsilon].sub.3]) and Rot([n.sub.km], [[epsilon].sub.4]) are mirror slope error matrix which could all be determined by function Rot(e, [beta]) as shown in (3). Rot(e, [beta]) could be used to realize the vector P [member of] [R.sup.1x3] rotate [beta] angle around the unit vector e = [[e.sub.x], [e.sub.y], [e.sub.z]]; [n.sub.3] = ([N.sub.km] x [e.sub.x])/[absolute value of [N.sub.km] x [e.sub.x]]; = [1,0,0].

[N.sup.1.sub.km] = [N.sub.km] x Rot([n.sub.3], [[epsilon].sub.3]) x Rot([N.sub.km], [[epsilon].sub.4]), (2)

[mathematical expression not reproducible]. (3)

In ideal tracking condition, the focal axis of the dish concentrator is parallel to the center line S (S = [0,0, - 1.0]) in the sun cone (sun half-apex angle [delta] = 4.65 mrad). However, when the EATD exists tracking error (elevation and azimuth tracking error of the EATD is [[epsilon].sub.1] and [[epsilon].sub.2], respectively, Figure 2(b)); the point [p.sub.km] moving to the position [p.sup.1.sub.km] and the corresponding normal vector is [N.sup.2.sub.km]. As well as the included angle [[omega].sub.err] between the focal axis of the dish concentrator and center line S in sun cone, this angle [[omega].sub.err] is the total tracking error of the dish concentrator. This could be calculated by (4), (5), (6), (7), where n is elevation tracking axis vector of the EATD, [n.sub.1] = [1, 0, 0]; [n.sub.2] is the azimuth tracking axis vector of the EATD, [n.sub.2] = [0, -cos[beta], sin[beta]]; the angle [beta] is solar elevation angle (i.e., working elevation angle of the dish concentrator); g is a position vector of the crosspoint g of the elevation tracking axis and azimuth tracking axis in the EATD (Figures 1 and 2(b)), g = [0, gy, gz]. Noteworthiness, the tracking error of the dish concentrator originated in the EATD tracking error. This tracking error transmission related to the working elevation angle [beta] of the dish concentrator, which could be clearly seen in (6).

[p.sup.1.sub.km] = ([p.sub.km] - g) x Rot([n.sub.1], [[epsilon].sub.1]) x Rot([n.sub.2], [[epsilon].sub.2]) + g, (4)

[N.sup.2.sub.km] = ([N.sup.1.sub.km] x Rot([n.sub.1], [[epsilon].sub.1]) x Rot([n.sub.2], [[epsilon].sub.2]), (5)

[N.sup.1.sub.F] = ([N.sub.F] x Rot([n.sub.1], [[epsilon].sub.1]) x Rot([n.sub.2], [[epsilon].sub.2]), (6)

[[omega].sub.err]([[epsilon].sub.1], [[epsilon].sub.2], [beta]) = arcos[sin([[epsilon].sub.1] + [beta]) x sin [beta] x (1 - cos [[epsilon].sub.e]) + cos [[epsilon].sub.1] x cos [[epsilon].sub.2]]. (7)

When the dish concentrator system exists tracking error, the plane receiver (located at the focal plane of dish concentrators) would rotate with the dish concentrator. This time, the focal point F moves to position [F.sub.1] and the corresponding plane receiver normal vector is [N.sup.1.sub.F], which can be calculated by (8) and (6), respectively, where F = [0, 0, f]. These two vectors can determine the space equation of the plane receiver.

[F.sub.1] = (F - g) x Rot([n.sub.1], [[epsilon].sub.1]) x Rot([n.sub.2], [[epsilon].sub.2]) + g. (8)

In this work, the plane receiver is a square with side length is L, which can be discreted into a lot of square grid (discrete parameter H = U). The calculation process of the flux distribution on the plane receiver as is as follows: (1) based on the law of specular reflection, the reflected ray equation of the incident sunlight [I.sub.wn] on the mirror point [p.sup.1.sub.km] is calculated by (4) and (5); (2) the intersection point q of the reflected ray and plane receiver is calculated, q = [[x.sub.1], [y.sub.1], [z.sub.1]; (3) the grid number hu containing point q is found and then increasing the sunlight [I.sub.wn] carrying energy (Section 2.2) to the grid hu; (4) following these steps, calculating every effective grid in the dish concentrator whole sunlight transmission process could get the flux distribution results on the plane receiver.

2.2. Discretization Model of the Sun Cone. The solar radiation intensity in a sun cone gradually decreases from the center to the border. The nonuniform energy distribution model is given as [27]:

[mathematical expression not reproducible], (9)

where [gamma] is the intersection angle of the discrete ray and center ray S (Figure 3); [I.sub.1] is the radiation intensity of the center ray.

The sun cone received by different positions in the dish concentrator mirror was assumed to be exactly identical. Thus, the sun cone was established on the point O and discretization as shown in Figure 3(b). In details, the solar disk on the focal plane was established (the radius is [R.sub.s], [R.sub.s] = f x tan([delta])) and it adopts the square envelope and then discrete to square grid (W = N). Thus, the position vector of the point [Q.sub.wn] is

[Q.sub.wn] = (w - 0.5) 2[R.sub.s]/W - [R.sub.s], [R.sub.s] - (n - 0.5) 2/[R.sub.s]N, f]. (10)

Based on the principle of energy conservation equation, the dish concentrator and sun cone before and after the discretization should satisfy the following relationship:

[mathematical expression not reproducible], (11)

where [e.sub.0] is the radiation intensity of a single sun cone; [[gamma].sub.wn] = arctan([d.sub.wn]/f) is the intersection angle of the sunlight ray wn and central sunlight; [d.sub.wn] = [parallel][Q.sub.wn] - F[parallel]; [W.sub.0] is the DNI value (W/[m.sup.2]); [A.sub.dish] = [[summation].sup.K.sub.k=1] [[summation].sup.M.sub.m=1] ([f.sub.km] x [A.sub.km]) is aperture total area of the dish concentrator. If [d.sub.wn] [less than or equal to] [R.sub.s], [f.sub.wn] = 1; if [d.sub.wn] > [R.sub.s], [f.sub.wn] = 0. Thus, [I.sub.1] can be written as follows:

I = [W.sub.0]/[[summation].sup.W.sub.w=1] [[summation].sup.N.sub.n=1] [f.sub.km] [1 - 0.5138 [([[gamma].sub.wn]/[delta]).sup.4]]. (12)

Finally, the sunlight vector [I.sub.wn] = -[Q.sub.wn]/[parallel] [Q.sub.wn] [parallel] and the energy carried by the sunlight can be determined by (9) and (12).

2.3. Model Validation. Based on the optical model and sun cone model in this paper, the ray tracing codes were compiled in Visual C++ 6.0 software for calculating the flux distribution on the focal plane of the dish concentrator. The Jeter theoretical results [28] were adopted to check the validity of the optical model and ray tracing codes and the comparison results as shown in Figure 4(a). We can see that the ray tracing calculation results in good agreement to Jeter's theoretical results, which verifies the validity of the flux distribution model in this paper. The flux distribution contour map calculated by ray tracing method in this paper is shown in Figure 4(b).

3. Focal Spot Characterization

Due to the reason that the tracking error of the EATD would cause the tracking error of the dish concentrator, the position and geometry of the focal spot on the plane receiver would all change. The focal spot moves and extends to the upper right direction, and the flux contour is approximately ellipse shaped but the center position of the ellipse is different, such as in Figure 5. In this paper, using a flux distribution map, flux distribution curve along the x-axis and y-axis, intercept factor, centroid position, and ellipse geometry (a flux contour curve) are used to describe the focal spot. Where centroid position and ellipse geometry could describe the movement direction and geometrical extent of the focal spot, respectively.

The intercept factor is defined as the ratio of the energy on the plane receiver (radius is [R.sub.in]) to the total received energy [14, 19]. The centroid coordinates of the focal spot could be calculated by (13), where f(h, u) is the flux density of the grid hu on the plane receive. Ellipse equation as shown in (14), including five parameters: center point coordinate s([x.sub.s], [y.sub.s]), minor semiaxis b, major semiaxis a, and intersection angle [[phi].sub.2] with the major axis and axis [x.sub.1]. The grid of the flux density equal or approaching to [C.sub.t] x [W.sub.0] was found on the receiver and then the center point coordinates of the grid were extracted (as the ellipse boundary point) for fitting (16) to get the ellipse parameters using the least square method ([C.sub.t] is a flux density extraction threshold value). The searching process of the ellipse boundary point is that according to the plane receiver discrete grid, from top to bottom (y-axis direction), to find the two adjacent grid for the first time, satisfy f (h, u) [less than or equal to] [C.sub.t] x [W.sub.0] [less than or equal to] f (h, u + 1. Then we selected one grid center as the boundary point (i.e., ellipse up boundary point) which is the flux density closest to the [C.sub.t] x [W.sub.0]. In the same way, from bottom to the top to find the ellipse bottom boundary point.

[mathematical expression not reproducible], (13)

[mathematical expression not reproducible]. (14)

4. Result and Discussion

4.1. Concentrator Tracking Error. The tracking error of the EATD would cause the tracking error of the dish concentrator, but the tracking error of the dish concentrator may not be equal to the tracking error of the EATD. This tracking error transmission (i.e., tracking error of the EATD transmitted to dish concentrator) is related to the working elevation angle [beta] of the dish concentrator. In Figure 6, when the EATD only has azimuth tracking error, the total tracking error [[omega].sub.err] of the dish concentrator decreases with the increase of elevation angle [beta] and the error [[omega].sub.err] decreases to 0 rad when the angle [beta] is 90[degrees](curve 1~5 in Figure 6(a)). When the EATD has both elevation and azimuth tracking error, the error [[omega].sub.err] decreases with the increase of the elevation angle [beta] and the error [[omega].sub.err] is equal to the elevation tracking error [[epsilon].sub.1] of the EATD when the angle [beta] is 90[degrees], which could be clearly seen in Figure 6(a) (curve 7~11) and Figure 6(b). That is to say, the increase of the working elevation angle [beta] of the dish concentrator would vanish or reduce the influence from the azimuth tracking error of the EATD to the dish concentrator. This phenomenon is called vanish-reduction effect in tracking error transmission by us. However, the elevation tracking error of the EATD transmit to the dish concentrator has no vanish-reduction effect, that is, the elevation tracking error of the dish concentrator equal to the elevation tracking error of the EATD. Thus, the concentrator total tracking error in (6) can be rewritten as follows:

[mathematical expression not reproducible], (15)

where [[epsilon].sub.2r] is azimuth tracking error of the dish concentrator from the azimuth tracking error [[epsilon].sub.2] of the EATD. Obviously, only when the angle [beta] = 0, then [[omega].sub.err]([[epsilon].sub.1], [[epsilon].sub.2], 0) = [square root of [([[epsilon].sub.2]).sup.2] + [([[epsilon].sub.1]).sup.2], which is worth to be noticed.

The vanish-reduction effect produce essence is because the focal axis of the dish concentrator and the tracking axis of the EATD are not vertical. The intersection angle between the focal axis of the dish concentrator and azimuth tracking axis of the EATD is (90[degrees] - [beta]) (Figure 2(b)). The vanish-reduction effect is more obvious with the decrease of the angle (90[degrees] - [beta]). Meanwhile, the influence of the elevation tracking error of the EATD on the dish concentrator was reduced. Due to the elevation tracking axis of the EATD always vertical to focal axis of the dish concentrator, there is no vanish-reduction effect in the elevation tracking error of the EATD transmission to the dish concentrator. Therefore, we could appropriately decrease the manufacture precision of the azimuth tracking device in the EATD to reduce production costs. Besides, it also could develop a new type of double-axis tracking device to increase the vanish-reduction effect, such as a spin-elevation double-axis tracking device for dish concentrator system in Section 4.4. On the other hand, it is critical to evaluate or predict the long-term performance of the concentrator system [14, 15]. If tracking error of the dish concentrator all comes from output angle error of the EATD, then according to the probability distribution of EATD output angle error model (based on statistical methods) and tracking error transmission law by 15, we can evaluate the long-term performance of dish concentrator system, which will be our future research work.

4.2. Flux Distribution on the Focal Plane. In this paper, we take the typical 38 kW XEM-Dish system (Figure 1) as an example to study the effect of the tracking error of the EATD on the focal spot characteristics of the focal plane. The calculation parameters include concentrator radius [R.sub.1] = 8.85 m, focal length f = 9.49 m, [[theta].sub.m] = 30[degrees], [R.sub.2] = 0.9m, mirror reflectivity is 0.90, absorption rate of the plane receiver is 100%, DNI value [W.sup.0] = 1000 W/[m.sup.2], [C.sub.t] =50, mirror discrete parameter K = 275, sun cone discrete parameter W = 201, and total tracking rays are 1.46 x [10.sup.9].

The influence of the tracking error of EATD on the flux curve and flux contour map of the focal spot are shown in Figures 7 and 8, respectively. From Figure 7(a), the movement and extent of the focal spot along the y-axis negative direction intensify with increase of the elevation tracking error of the EATD, but the flux peak value changes small. Although the mirror slope error of the dish concentrator would increase the geometry size of the focal spot and decrease the flux peak value of the focal spot [24], it would not influence the centroid position of the focal spot. From the comparison of Figures 8(a) and 8(b) and Figures 8(e) and 8(f), respectively, we could clearly see that. The DNI value would only influence the flux density value of the focal spot [25] but would not influence the centroid position of the focal spot (comparison of Figures 8(d) and 8(e)). In conclusion, the centroid position of the focal spot on the focal plane depends on the tracking error of the dish concentrator but has nothing to do with the mirror slope error and DNI. Therefore, the centroid position of the focal spot on the focal plane could be used for estimate and calibration of the tracking error of the dish concentrator (Section 4.3).

The azimuth and elevation tracking error of the EATD only have influence on the centroid coordinates of the focal spot along x-axis and y-axis, respectively. For example, the coordinate of the focal spot centroid along y-axis is -23.45 mm in both Figure 8(a) ([beta] = 30[degrees], [[epsilon].sub.1] = 2 mrad, and [[epsilon].sub.2] = 8mrad) and Figure 8(g) ([beta] = 0[degrees], [[epsilon].sub.1] = 2mrad, and [[epsilon].sub.2] = 0mrad); the coordinate of the focal spot centroid along x-axis is -83.20 mm in Figure 8(a) (double tracking error) and Figure 8(h) (single tracking error); the azimuth and elevation tracking error of the EATD are all 6.9 mrad ([beta] = 0) in Figures 8(d)-8(f), so the coordinate of the focal spot centroid along x-axis and y-axis are the same -83.2 mm. These show that the influence of the azimuth and elevation tracking error of the EATD on the centroid position of the focal spot is independent. Moreover, the position of the point g (Figure 2(b)) has no influence for the focal spot characteristics, that is, to say the EATD could be installed on every position in a dish concentration system and does not need to modify the tracking program of the EATD.

4.3. Centroid Moving Function and Intercept Factor. The functional relationship between the centroid movement distance of the focal spot on focal plane and EATD tracking error is called centroid movement function; this function could be used for evaluating and quantitatively calibrating the tracking error of the EATD. The influence curve of the azimuth tracking error of the EATD on the centroid movement distance is given in Figure 9. We can see that the movement distance of the focal spot centroid has linear relationship with the azimuth tracking error of the EATD. The centroid moving function is x = -12.06[[epsilon].sub.2] and y = -12.06[[epsilon].sub.1] when the elevation angle of the dish concentrator is 0[degrees]. Considering the vanish-reduction effect, the centroid movement function of the focal spot in 38 kW XEM-Dish system is as follows:

x = -12.06arcos[cos2(90 - [beta]) x (1 - cos [[epsilon].sub.2]) + cos [[epsilon].sub.2]],

y = -12.06[[epsilon].sub.1]. (16)

Noteworthiness, the simplified formula y = -f x tan([[epsilon].sub.1]) is usually not equal to (16). For example, when the [[epsilon].sub.1] = 10 mrad, the centroid movement distance is -120.60 mm calculated by (16), while y = -f x tan([[epsilon].sub.1]) = -9490 x tan(0.010) = -94.90 mm and there is 25.70 mm difference between them. Therefore, the centroid movement function of the focal spot on focal plane needs to be accurately determined by ray tracing method.

The centroid movement function is used to guide the calibration of the tracking error of the dish concentrator system as the following: measuring the focal spot on the focal plane in the dish concentrator system and extracting the centroid coordinates of the focal spot; adopting the centroid coordinates to calculate the double-axis tracking error of the EATD by (16) and then calibrating the tracking error of the dish concentrator. Of note, the assessment and calibration of the tracking error of the dish concentrator system are not suitable for the large working elevation angle occasion (especially close to the elevation angle is 90[degrees]). Because the vanish-reduction effect is becoming intense during this time, it would lead the centroid movement distance of the focal spot along the x-axis to be too small to be discovered. For example, the centroid movement distance of the focal spot equals 0 mm when the working elevation [beta] equals 90[degrees] (Figure 9), so the effect of the azimuth tracking error of the EATD on the dish concentrator is completely covered, but this effect appears again after the working elevation angle ft decreases such as [beta] = 30[degrees].

The double-axis tracking error of the EATD and mirror slope error has significant influence on the intercept factor (Figure 10). When the influence of the EATD tracking error only exists, the intercept factor curve would be the same as long as the total tracking error [[omega].sub.err] of the dish concentrator is the same. Although the single tracking error of the EATD is different, such as curve 5, 6, and 7 in Figure 10, the intercept factor curves are the same. Another example is curve 2 and 3 which are the same.

4.4. Application of the Vanish-Reduction Effect. The solar heliostat in tower power plant also adopts the EATD [9] as shown in Figure 11. The characteristic axis of the heliostat is the mirror normal vector N, and azimuth tracking axis of the EATD is vertical to the ground level. Thus, the intersection angle between the characteristic axis N and azimuth tracking axis of the EATD attains the maximum and minimum when the solar elevation angle [beta] is 0[degrees] and 90[degrees], respectively, given by (17), where the intersection angle [[beta].sub.1] between the reflection ray (aimed at target point on the receiver tower) and horizontal plane satisfied tan([[beta].sub.1]) = H/[L.sub.1]; H is the height between the target point and heliostat center; [L.sub.1] is horizontal distance between the target point and heliostat center. The angle fi1 increases with the decrease of the distance between the heliostat and the receive tower, and the angle [[beta].sup.1.sub.min] and [[beta].sup.1.sub.max] will decrease at the same time, which could increase the vanish-reduction effect in the azimuth tracking error of the EATD transmit to heliostat for improving the tracking performance of the heliostat. According to this regulation, on the basis of guarantying the tracking accuracy of the heliostat, this could decrease the manufacturing cost of the heliostat according to gradual or gradient increase of the precision requirements of the azimuth tracking device in the EATD near or far from the receiver tower.

[[beta].sup.1.sub.min] = 0.5(90[degrees] - [[beta].sub.1]) [beta] = 90[degrees],

[[beta].sup.1.sub.max] = 90[degrees] - 0.5 [[beta].sub.1]) [beta] = 0[degrees]. (7)

In order to increase the vanish-reduction effect to reduce the effect of the tracking error of the double-axis device on the dish concentrator, a dish concentrator system using a spin-elevation double-axis tracking device (SETD) is proposed (Figure 12). The spin-elevation double-axis tracking was put forward by Chen et al. [8] and they defined the spin axis is a line between the heliostat center and fixed target and the elevation tracking axis in the SETD vertical to the spin axis. The SETD is widely used in the solar concentrator field where the solar concentrator can move but the heat absorber receiver is fixed [29]. However, the dish concentrator and receiver move together in the dish concentrator system and there is no fixed target in this system. Thus, we defined any direction as the spin axis. In Figure 12, angle [[beta].sub.2] is the inclination between the azimuth tracking axis of the SETD and horizontal plane. The spin tracking error [[epsilon].sub.2r] of the dish concentrator caused by spin tracking error of the SETD transmission to the dish concentrator could be calculated by (18). From Figure 13, we could clearly see that the error [[epsilon].sub.2r] reduced to 0mrad first and then increased with the increase of the working elevation angle [beta] of the dish concentrator, and the error [[epsilon].sub.2r] equals 0mrad when the angle p equals to the inclination angle [[beta].sub.2]. Therefore, selecting the appropriate angle [[beta].sub.2] could make the vanish-reduction effect adopt the SETD's obvious intensity than the EATD during the tracking error transmission, in order to increase the influence of the tracking error of the tracking device on the dish concentrator. Such as the angle [[beta].sub.2] = 60[degrees] in Figure 13, within the elevation angle [beta] range of 0[degrees]~75[degrees], the spin tracking error [[epsilon].sub.2r] of the dish concentrator adopted the SETD form obviously less than the EATD (i.e., [[beta].sub.2] = 90[degrees]). Thus, appropriately reducing the manufacture precision requirements of the spin tracking device in SETD could decrease the manufacturing cost.

[[epsilon].sub.2r]([[beta].sub.2], [beta]) = arcos[[cos.sup.2]([beta] - [[beta].sub.2]) x (1 - cos [[epsilon].sub.2]) + cos [[epsilon].sub.2]]. (18)

5. Conclusion

In this paper, a flux distribution model of the focal plane of a dish concentrator system has been developed based on ray tracking method for researching the effect of the mirror surface slope error, DNI and EATD tracking error on the focal spot characteristics, and the error transmission law of the tracking error of the double-axis tracking device transmission to dish concentrator system. The main conclusions are as follows:

(1) The azimuth tracking error [[epsilon].sub.2r] of the dish concentrator is transmitted by the azimuth tracking error [[epsilon].sub.2] of the EATD. When the error [[epsilon].sub.2] is constant, the error [[epsilon].sub.2r] decreases with the increase of the working elevation angle [beta] of the dish concentrator and the error [[epsilon].sub.2r] equals to 0 mrad when [beta] = 90[degrees]. We call this phenomenon the called vanish-reduction effect in tracking error transmission. The vanish-reduction effect-produced essence is caused by the nonvertical of the focal axis being nonvertical to the dish concentrator and tracking axis of the EATD, and the vanish-reduction effect is increased with the decrease of the nonvertical angle, which could improve the tracking accuracy of the dish concentrator. Due to the elevation tracking axis of the EATD always vertical to the focal axis, the elevation tracking error of the dish concentrator equals to the elevation tracking error of the EATD. In addition, the intercept factor only depends on the total tracking error of the dish concentrator when only considering the tracking error.

(2) The centroid position of the focal spot on the focal plane depends on the tracking error of the dish concentrator without the effects of the mirror slope error and DNI. Only the azimuth tracking error and elevation tracking error of the EATD have influence on the centroid's movement distance of the focal spot along x-axis and y-axis, respectively, and the centroid's moving distance has linear function relationship with the EATD tracking error. For 38 kW XEM-Dish system, the centroid's moving function is x = -12.06[[epsilon].sub.2r] andy = -12.06[[epsilon].sub.1]. This function could be used to evaluate and quantitative calibration of the tracking error of the dish concentrator system but is not suitable for the large working elevation angle occasion.

(3) The vanish-reduction effect also exists in the azimuth tracking error of the EATD transmission to the solar heliostat. Through decreasing the distance between the heliostat and receiver tower, increasing the height of the receiver tower and solar elevation angle could increase the vanish-reduction effect in order to improve the tracking performance of the heliostat. Besides, a dish concentrator system using a spin-elevation double-axis tracking device (SETD) is proposed, which could increase the influence of the spin tracking error of the SETD on the dish concentrator system and improve the tracking performance of the dish concentrator system.

Nomenclature [A.sub.km]: Daylighting area of the mirror grid km ([m.sup.2]) f: Focal length of the dish concentrator (m) K, M: Discretion numbers of the dish concentrator mirror N, W: Discretion numbers of the sun cone [R.sub.1]: Radius of the dish concentrator (m) [R.sub.2]: Void field radius of the dish concentrator (m) [W.sub.0]: DNI value (W[m.sup.-2]) [[omega].sub.err]: Total tracking error of the dish concentrator (mrad), [[omega].sub.err] = [square root of [([[epsilon].sub.2r]).sup.2] + [([[epsilon].sub.1]).sup.1]. Greek symbols [delta]: Solar half angle (mrad) [sigma]: Standard deviation of the mirror slope error (mrad) [[epsilon].sub.1]: Elevation tracking error of the EATD (mrad) [[epsilon].sub.2]: Azimuth tracking error of the EATD (mrad) [[epsilon].sub.2r]: Azimuth angle or spin angle tracking error of the dish concentrator (mrad) [beta]: Concentrator working elevation angle (degrees) [[beta].sub.1]: The intersection angle between the reflection ray (at heliostat center) and the horizontal plane (degrees) [[beta].sub.2]: Angle between spin axis and horizontal plane (degrees) [[theta].sub.m]: Lack angle of the dish concentrator (degrees). Abbreviations EATD: Elevation-azimuth double-axis tracking device DNI: Direct normal irradiation. Subscripts km: Concentrator mirror grid at position (k, m) wn: Sun cone grid at position (w, n).

https://doi.org/10.1155/2018/9046127

Conflicts of Interest

The authors declare that there is no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work is supported by the National Natural Science Foundation of People's Republic of China (nos. 51576061 and 51641504), Hunan Provincial Natural Science Foundation of China (no. 2016JJ2052), and Hunan Province Postgraduate Innovation Project Foundation of People's Republic of China (nos. CX2016B549 and CX2017B628).

References

[1] H. Sayyaadi and M. S. Boroujeni, "Conceptual design, process integration, and optimization of a solar Cu-Cl thermochemical hydrogen production plant," International Journal of Hydrogen Energy, vol. 42, no. 5, pp. 2771-2789, 2017.

[2] L. Nkhonjera, T. Bello-Ochende, G. John, and C. K. King'ondu, "A review of thermal energy storage designs, heat storage materials and cooking performance of solar cookers with heat storage," Renewable and Sustainable Energy Reviews, vol. 75, pp. 157-167, 2017.

[3] D. A. Baharoon, H. A. Rahman, W. Z. W. Omar, and S. O. Fadhl, "Historical development of concentrating solar power technologies to generate clean electricity efficiently--a review," Renewable and Sustainable Energy Reviews, vol. 41, pp. 996-1027, 2015.

[4] Y. Qiu, Y. L. He, M. Wu, and Z. J. Zheng, "A comprehensive model for optical and thermal characterization of a linear Fresnel solar reflector with a trapezoidal cavity receiver," Renewable Energy, vol. 97, pp. 129-144, 2016.

[5] G. Li and Y. Jin, "Optical simulation and experimental verification of a Fresnel solar concentrator with a new hybrid second optical element," International Journal of Photoenergy, vol. 2016, Article ID 4970256, 8 pages, 2016.

[6] F. Chen, M. Li, R. Hassanien Emam Hassanien et al., "Study on the optical properties of triangular cavity absorber for parabolic trough solar concentrator," International Journal of Photoenergy, vol. 2015, Article ID 895946, 9 pages, 2015.

[7] J. Wang, J. Wang, X. Bi, and X. Wang, "Performance simulation comparison for parabolic trough solar collectors in China," International Journal of Photoenergy, vol. 2016, Article ID 9260943, 16 pages, 2016.

[8] Y. T. Chen, K. K. Chong, T. P. Bligh et al., "Non-imaging, focusing heliostat," Solar Energy, vol. 71, no. 3, pp. 155-164, 2001.

[9] O. Behar, A. Khellaf, and K. Mohammedi, "A review of studies on central receiver solar thermal power plants," Renewable and Sustainable Energy Reviews, vol. 23, pp. 12-39, 2013.

[10] J. Coventry and C. Andraka, "Dish systems for CSP," Solar Energy, vol. 152, pp. 140-170, 2017.

[11] H. Ma, G. Jin, X. Zhong, K. Xu, and Y. Li, "Optical design of a solar dish concentrator based on triangular membrane facets," International Journal of Photoenergy, vol. 2012, Article ID 391921, 5 pages, 2012.

[12] Q. Mao, L. Zhang, and H. Wu, "Study on concentrating characteristics of a solar parabolic dish concentrator within high radiation flux," International Journal of Photoenergy, vol. 2015, Article ID 787280, 7 pages, 2015.

[13] I. Reda and A. Andreas, "Solar position algorithm for solar radiation applications," Solar Energy, vol. 76, no. 5, pp. 577-589, 2004.

[14] R. O. Hughes, "Effects of tracking errors on the performance of point focusing solar collectors," Solar Energy, vol. 24, no. 1, pp. 83-92, 1980.

[15] V. Badescu, "Different tracking error distributions and their effects on the long-term performances of parabolic dish solar power systems," International Journal of Solar Energy, vol. 14, no. 4, pp. 203-216, 1994.

[16] V. Badescu, "Theoretical derivation of heliostat tracking errors distribution," Solar Energy, vol. 82, no. 12, pp. 1192-1197, 2008.

[17] Y. Hu and Y. Yao, "Optical analysis and output evaluation for a two-stage concentration photovoltaic system by using Monte Carlo ray-tracing method," Optik-International Journal for Light and Electron Optics, vol. 131, pp. 713-723, 2017.

[18] K. Shanks, N. Sarmah, J. P. Ferrer-Rodriguez et al., "Theoretical investigation considering manufacturing errors of a high concentrating photovoltaic of cassegrain design and its experimental validation," Solar Energy, vol. 131, pp. 235-245, 2016.

[19] X. L. Xia, G. L. Dai, and Y. Shuai, "Experimental and numerical investigation on solar concentrating characteristics of a sixteen-dish concentrator," International Journal of Hydrogen Energy, vol. 37, no. 24, pp. 18694-18703, 2012.

[20] H. Li, W. Huang, F. Huang, P. Hu, and Z. Chen, "Optical analysis and optimization of parabolic dish solar concentrator with a cavity receiver," Solar Energy, vol. 92, pp. 288-297, 2013.

[21] L. S. M. Castellanos, G. E. C. Caballero, V. R. M. Cobas et al., "Mathematical modeling of the geometrical sizing and thermal performance of a dish/Stirling system for power generation," Renewable Energy, vol. 107, pp. 23-35, 2017.

[22] G. Johnston, "Focal region measurements of the 20m 2 tiled dish at the Australian National University," Solar Energy, vol. 63, no. 2, pp. 117-124, 1998.

[23] S. Ulmer, W. Reinalter, P. Heller et al., "Beam characterization and improvement with a flux mapping system for dish concentrators," in ASME Solar 2002: International Solar Energy Conference, pp. 285-292, Reno, NV, USA, 2002.

[24] H. Lee, "The geometric-optics relation between surface slope error and reflected ray error in solar concentrators," Solar Energy, vol. 101, pp. 299-307, 2014.

[25] Q. Mao, Y. Shuai, and Y. Yuan, "Study on radiation flux of the receiver with a parabolic solar concentrator system," Energy Conversion and Management, vol. 84, pp. 1-6, 2014.

[26] M. H. Tan and K. K. Chong, "Influence of self-weight on electrical power conversion of dense-array concentrator photovoltaic system," Renewable Energy, vol. 87, pp. 445-457, 2016.

[27] M. M. Elsayed, K. A. Fathalah, and O. M. Al-Rabghi, "Measurements of solar flux density distribution on a plane receiver due to a flat heliostat," Solar Energy, vol. 54, no. 6, pp. 403-411, 1995.

[28] S. M. Jeter, "The distribution of concentrated solar radiation in paraboloidal collectors," Journal of Solar Energy Engineering, vol. 108, no. 3, pp. 219-225, 1986.

[29] B. H. Lim, K. K. Chong, C. S. Lim, and A. C. Lai, "Latitude-orientated mode of non-imaging focusing heliostat using spinning-elevation tracking method," Solar Energy, vol. 135, pp. 253-264, 2016.

Jian Yan (iD), Zi-ran Cheng, and You-duo Peng (iD)

Hunan Provincial Key Laboratory of Health Maintenance for Mechanical Equipment, Hunan University of Science and Technology, Xiangtan, Hunan 411201, China

Correspondence should be addressed to You-duo Peng; ydpeng@hnust.edu.cn

Received 25 July 2017; Revised 9 November 2017; Accepted 29 November 2017; Published 14 March 2018

Academic Editor: Alberto Alvarez-Gallegos

Caption: Figure 1: A 38 kW dish-Stirling solar thermal power system (38 kW XEM-Dish system). The diameter is 17.70 m and the focal length is 9.49 m of the parabolic dish concentrator and was built in Xiangtan Electric Manufacturing Group, China.

Caption: Figure 2: (a) Geometry parameters and mirror discrete of the dish concentrator and (b) schematic diagram of the tracking error and mirror slope error of the dish concentrator.

Caption: Figure 3: Discretization of the (a) sun cone and (b) solar disk.

Caption: Figure 4: (a) The comparison of the flux distribution result on the focal plane and (b) flux distribution contour map using ray tracing method in this paper (W/[m.sup.2]).

Caption: Figure 5: Geometrical parameters of the focal spot on the plane receive.

Caption: Figure 6: Influence of the tracking error of the EATD on the dish concentrator. 1: [[epsilon].sub.1] = 0 mrad, [[epsilon].sub.2] = 2 mrad; 2: [[epsilon].sub.1] = 0 mrad, [[epsilon].sub.2] = 4 mrad; 3: [[epsilon].sub.1] = 0 mrad, [[epsilon].sub.2] = 6 mrad; 4: [[epsilon].sub.1] = 0 mrad, [[epsilon].sub.2] = 8 mrad; 5: [[epsilon].sub.1] = 0 mrad, [[epsilon].sub.2] = 10 mrad; 6: [[epsilon].sub.1] = 10 mrad, [[epsilon].sub.2] = 0 mrad; 7: [[epsilon].sub.1] = 10 mrad, [[epsilon].sub.2] = 10 mrad; 8: [[epsilon].sub.1] = 8 mrad, [[epsilon].sub.2] = 8 mrad; 9: [[epsilon].sub.1] = 6 mrad, [[epsilon].sub.2] = 6 mrad; 10: [[epsilon].sub.1] = 4 mrad, [[epsilon].sub.2] = 10 mrad; 11: [[epsilon].sub.1] = 2 mrad, [[epsilon].sub.2] = 2 mrad.

Caption: Figure 7: Influence of the EATD tracking error on the flux distribution of the focal plane.

Caption: Figure 8: The flux distribution of the focal plane under the EATD tracking error (W/[m.sup.2]). (a) [sigma] = 0, [beta] = 30[degrees], [[epsilon].sub.1] = 2 mrad, and [[epsilon].sub.2] = 8 mrad; (b) [sigma] = 2 mrad, [beta] = 30[degrees], [[epsilon].sub.1] = 2 mrad, and [[epsilon].sub.2] = 8 mrad; (c) [sigma] = 0, [beta] = 0[degrees], [[epsilon].sub.1] = 6.9 mrad, and [[epsilon].sub.2] = 2 mrad; (d) [sigma] = 0, [beta] = 0[degrees], [[epsilon].sub.1] = 6.9 mrad, and [[epsilon].sub.2] = 6.9 mrad; (e) [sigma] = 0, [beta] = 0[degrees], [[epsilon].sub.1] = 6.9 mrad, [[epsilon].sub.2] = 6.9 mrad, and [W.sub.0] = 800 W/[m.sup.2]; (f) [sigma] = 2 mrad, [beta] = 0[degrees], [[epsilon].sub.1] = 6.9 mrad, [[epsilon].sub.2] = 6.9 mrad, and [W.sub.0] = 800 W/[m.sup.2]; (g) [sigma] = 0, [beta] = 0[degrees], [[epsilon].sub.1] = 2 mrad, and [[epsilon].sub.2] = 0; (h) [sigma] = 0, [beta] = 0[degrees], [[epsilon].sub.1] = 0, [[epsilon].sub.2] = 6.9 mrad or [beta] = 30[degrees], [[epsilon].sub.1] = 0, and [[epsilon].sub.2] = 8 mrad. Note: only Figures 8(e) and 8(f) are with [W.sub.0] = 800 W/[m.sup.2] and the other figures are with [W.sub.0] = 1000 W/[m.sup.2].

Caption: Figure 9: The relationship between the centroid movement distance of the focal spot and azimuth tracking error of the EATD.

Caption: Figure 10: Influence of the EATD tracking error and mirror slope error on the intercept factor. 1: [[omega].sub.err] = 0; 2: [beta] = 0[degrees], [[epsilon].sub.1] = 3.6, [[epsilon].sub.2] = 0, [[omega].sub.err] = 3.6; 3: [beta] = 0[degrees], [[epsilon].sub.1] = 3, [[epsilon].sub.2] = 2, [[omega].sub.err] = 3.6; 4: [beta] = 0[degrees], [[epsilon].sub.1] = 6, [[epsilon].sub.2] = 0; 5: [beta] = 30[degrees], [[epsilon].sub.1] = 2, [[epsilon].sub.2] = 8, [[omega].sub.err] = 7.184; 6: [beta] = 0[degrees], [[epsilon].sub.1] = 3.95, [[epsilon].sub.2] = 6, [[omega].sub.err] = 7.184; 7: [beta] = 0[degrees], [[epsilon].sub.1] = 7.184, [[epsilon].sub.2] = 0; 8: [beta] = 0[degrees], [[epsilon].sub.1] = 6, [[epsilon].sub.2] = 0, [sigma] = 2; 9: [beta] = 0[degrees], [[epsilon].sub.1] = 10, [[epsilon].sub.2] = 0; 10: [beta] = 0[degrees], [[epsilon].sub.1] = 10, [[epsilon].sub.2] = 0, [sigma] = 2. Note: only curve 10 and 14 have the mirror slope error of [sigma] = 2 mrad, other curves have the mirror slope error of [sigma] = 0 mrad.

Caption: Figure 11: A schematic diagram of the elevation-azimuth double-axis tracking of a heliostat in a tower plant.

Caption: Figure 12: A schematic diagram of the spin-elevation double-axis tracking in dish concentrator system.

Caption: Figure 13: Error transmission law of the tracking error of the spin axis transmission to dish concentrator (spin axis tracking error is 4.0 mrad).

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Yan, Jian; Cheng, Zi-ran; Peng, You-duo |

Publication: | International Journal of Photoenergy |

Date: | Jan 1, 2018 |

Words: | 8308 |

Previous Article: | Renewable Energy Integration: Economic Assessment of Solar Energy to Produce Biodiesel at Supercritical Conditions. |

Next Article: | Review on Substrate and Molybdenum Back Contact in CIGS Thin Film Solar Cell. |

Topics: |