Solar Concentrator Consisting of Multiple Aspheric Reﬂectors

: This paper presents an o ﬀ -axis-focused solar concentrator system consisting of 190 aspheric reﬂectors, where the aperture radius of each reﬂector is 10 cm, and vertices of all reﬂectors are orderly arranged in the same plane. The aspheric parameters controlling the curvature of the reﬂectors are determined using coordinate transformations and the particle swarm optimization (PSO) algorithm. Based on these aspheric parameters, the light distribution of focal plane was calculated by the ray tracing method. The analyses show that the designed concentrator system has a spot radius of less than 1 cm and the concentration ratio over 3300:1 is achieved using only one reﬂection. The design results have been veriﬁed with the optical design software Zemax.


Introduction
Solar energy is often considered as one of the most viable renewable energy sources because it is unpolluted, abundant, and holds great promises for utilization and sustainable development in the near future [1,2]. Solar thermal and photovoltaic (PV) power generations are the two major applications of solar energy currently in use. In the case of solar thermal power generation, steam is generated and then current is produced using the steam turbine, whereas in the case of solar PV power generation, solar radiation is directly converted into the electricity [3][4][5]. Use of solar concentrator is an effective approach in solar PV power generation in order to increase the output power and to reduce the investment cost [6].
In the last five decades, various forms of devices making use of solar concentrators to heat a heat transfer medium or working fluid were developed providing power to mechanical equipment. The models of the concentrators can be cataloged into parabolic dishes [7,8], Fresnel lenses [9,10], parabolic troughs [11], special bowls [12], etc. For example, a parabolic dish concentrator is mounted using sun-tracking mechanism and its motion is well controlled. The solar radiation is focused on the thermal receiver placed at the focal point. A trough concentrator is a solar linear focusing device, and the solar radiation is focused on tubular heat absorbers placed along its focal line. A trough concentrator could produce a temperature of about 100 • C, and dish systems could reach a concentration temperature of 1500 • C at the receiver [13]. The sunlight concentration ratio must be improved in order to increase the focusing temperature. The theoretical concentration ratio for trough solar concentrators and dish concentrators were 285:1 [14] and 2240:1 [15,16], respectively. At present, the concentration ratio for a large trough solar concentrator made in domestic enterprises is generally up to 80:1 and for a large dish solar concentrator, it is usually about 1000:1. Figure 1 depicts a new idea about a concentrator system. This system has multiple reflectors which reflectively focus sunlight on the same small areas. If the number of such reflectors is large enough, the concentration ratio will also be high enough, and the focused spot can reach very high temperatures.
Reflecting surfaces employ even-aspheric surfaces whose expression is shown in Equation (1). (1). The first term in the above equation corresponds to a quadric surface, C is the apex curvature of a quadric surface, 1 + a 2 is the coefficient of a quadric surface, a 2 relates with the eccentricity of a quadric surface, a 4 , a 6 , a 8 are higher order polynomial coefficients of an aspheric equation, respectively.
Energies 2019, 12, x FOR PEER REVIEW 3 of 16 (PSO) method. PSO algorithm was successfully implemented for designing the laser beam shaping lenses and LED collimator systems [39][40][41] in our recent works. Figure 1 depicts a new idea about a concentrator system. This system has multiple reflectors which reflectively focus sunlight on the same small areas. If the number of such reflectors is large enough, the concentration ratio will also be high enough, and the focused spot can reach very high temperatures.

Ideas of Focusing
Reflecting surfaces employ even-aspheric surfaces whose expression is shown in Equation (1). 8 in Equation (1). The first term in the above equation corresponds to a quadric surface, C is the apex curvature of a quadric surface, 1 + a2 is the coefficient of a quadric surface, a2 relates with the eccentricity of a quadric surface, a4, a6, a8 are higher order polynomial coefficients of an aspheric equation, respectively.

Analysis for Focusing Action of a Single Aspheric Reflector
A Cartesian coordinate system is established as shown in Figure 2; the coordinate origin O is located at the apex of an aspheric reflector. In Figure 2, we draw three incident-reflecting rays passing through the reflector, one of them travels through the vertex O of the aspheric surface, and the other two do not travel through O. Let the angle between an incident ray passing through the vertex O and the surface normal at O be θ (incident angle), then the angle of reflection must also be θ. As long as this incident ray is given, the reflected ray is also determined, which has nothing to do with the shape of the aspheric surface. Therefore, the focal point must be placed on the reflected light ray passing through the apex point O. Let the incident solar beam be parallel to the xoy plane, then the Z coordinate of any point on the reflected light ray through O must be zero. Take any point Pf (x′, y′,0) on this ray as the center of a focal spot, accordingly, the relationship between x′ and y′ is x′ = −y′tanθ. If a 4 = a 6 = a 8 = 0, then Equation (1) is simplified to y = , this is a rotational quadric surface around the symmetry Y-axis. Equation (1) respectively defines a hyperboloid, paraboloid, ellipsoid, and spherical with different a 2 values.

Analysis for Focusing Action of a Single Aspheric Reflector
A Cartesian coordinate system is established as shown in Figure 2; the coordinate origin O is located at the apex of an aspheric reflector. In Figure 2, we draw three incident-reflecting rays passing through the reflector, one of them travels through the vertex O of the aspheric surface, and the other two do not travel through O. Let the angle between an incident ray passing through the vertex O and the surface normal at O be θ (incident angle), then the angle of reflection must also be θ. As long as this incident ray is given, the reflected ray is also determined, which has nothing to do with the shape of the aspheric surface. Therefore, the focal point must be placed on the reflected light ray passing through the apex point O. Let the incident solar beam be parallel to the xoy plane, then the Z coordinate of any point on the reflected light ray through O must be zero. Take any point Pf (x , y , 0) on this ray as the center of a focal spot, accordingly, the relationship between x and y is x = −y tanθ.
is approximately obtained by the following method described in brief [14]. Given that the unit vector ) , , ( We can see from the above expression that 2 . As long as these parameters are properly selected, we can find a suitable 2 Q  directing the neighborhood of the focal point, thereby the objective of focusing the sunbeams on a small area can be achieved. Since it has been assumed that the incident solar beam is parallel to the xoy plane, thus γ1 = 0 in ) , , (

Analyses for Collecting Solar Light with a Single Reflector
As shown in Figure 2, take a point Pf (x′, y′, 0) on a reflected light ray passing through the vertex O as a focus point. P′′ (x′′, y′′, z′′) is a point on a reflected light ray not passing through O, it is known from Figure 2 where, x, y, and z are coordinates of intersection P between the light ray and the aspheric reflector. Let y′′ = y′, then the distance between point P'' and the focus point P' is given by  The reflected light vector Q 2 (α 2 , β 2 , γ 2 ) is approximately obtained by the following method described in brief [14]. Given that the unit vector Q 1 (α 1 , β 1 , γ 1 ) of incidence ray, a point on an incident ray and an aspheric equation, a ray intersection point P (x, y, z) at the aspheric reflecting surface can be gotten using approximate methods, thus the normal vector N(α N , β N , γ N ) at the intersection point can be determined by the aspheric equation. According to the vector law of reflection, the unit direction vector Q 2 (α 2 , β 2 , γ 2 ) of reflected ray can be found from Q 1 and N.
We can see from the above expression that Q 2 is only a function of Q 1 and N, N is only a function of aspheric parameters C, a 2 , a 4 , a 6 , a 8 and intersection point P, and intersection point P is determined by Q 1 and the aspheric parameters. When Q 1 is given, intersection point P is just determined by the aspheric parameters, therefore, N is also obtained only from aspheric parameters. In this case, given Q 1 , Q 2 is just a function of the aspheric parameters C, a 2 , a 4 , a 6 , a 8 . As long as these parameters are properly selected, we can find a suitable Q 2 directing the neighborhood of the focal point, thereby the objective of focusing the sunbeams on a small area can be achieved. Since it has been assumed that the incident solar beam is parallel to the xoy plane, thus γ 1 = 0 in Q 1 (α 1 , β 1 , γ 1 ).

Analyses for Collecting Solar Light with a Single Reflector
As shown in Figure 2, take a point P f (x , y , 0) on a reflected light ray passing through the vertex O as a focus point. P" (x", y", z") is a point on a reflected light ray not passing through O, it is known from Figure 2 that a relation of x", z", and y" is determined by the following formula: where, x, y, and z are coordinates of intersection P between the light ray and the aspheric reflector. Let y" = y , then the distance between point P" and the focus point P is given by R = (x − x ) 2 + (z − 0) 2 , the smaller the R value, the better the focusing effect. It can be seen from the above analysis that R is also only a function of aspheric parameters. The appropriate R-value can be gotten as long as we properly select these parameters. In this paper, let R be less than a certain value, this certain value is used as the maximum allowable value of the fitness function of the PSO algorithm, and aspheric parameters C, a 2 , a 4 , a 6 , a 8 can be found by PSO techniques. Here we would like to point out that in order to better understand our specific implementation of the algorithm, readers should have knowledge about PSO.
For example, let the incident angle θ = 25 • , and the focus point be P (400*tan25 • , 400, 0), then Evenly taking multiple incident rays on a reflecting surface, their corresponding reflected rays have the respective R-values, and all R-values are summed up, then we have: Here, i denotes the ith light ray.
i R i serves as the fitness function of PSO, aspheric parameters C, a 2 , a 4 , a 6 , a 8 can be obtained by seeking the particle with the smallest i R i using the PSO method [41].

The Construction of a New Coordinate System and Coordinate Transformations between the Old and the New Coordinate System
As shown in Figure 3, the vertex of another reflector of a solar concentrator system is not at the origin O of the coordinate system, but at point O (a, b, c). The point P f (−y f tanθ, y f , 0) is still the focus point, then a vector along the direction of O P f is gotten as long as we properly select these parameters. In this paper, let R be less than a certain value, this certain value is used as the maximum allowable value of the fitness function of the PSO algorithm, and aspheric parameters can be found by PSO techniques. Here we would like to point out that in order to better understand our specific implementation of the algorithm, readers should have knowledge about PSO.
For example, let the incident angle θ = 25°, and the focus point be P′ (400*tan25°, 400, 0), then Evenly taking multiple incident rays on a reflecting surface, their corresponding reflected rays have the respective R-values, and all R-values are summed up, then we have:

186.5231
Here, i denotes the i th light ray.  i i R serves as the fitness function of PSO, aspheric parameters can be obtained by seeking the particle with the smallest  i i R using the PSO method [41].

The Construction of a New Coordinate System and Coordinate Transformations between the Old and the New Coordinate System
As shown in Figure 3, the vertex of another reflector of a solar concentrator system is not at the origin O of the coordinate system, but at point O′ (a, b, c). The point Pf (−yf tanθ, yf, 0) is still the focus point, then a vector along the direction of O′Pf is The vector along the bisector of the angle is: The direction cosines of the vector of angle bisector are: Let the normal vector to the plane constituted by the vector Q 1 and the vector The direction cosines of Take the vector along the angle bisector as the positive direction of the Y axis of a new coordinate system, and − −−− → O Z as the positive direction of the Z axis of a new coordinate system, then the positive direction of the X axis is also determined. In the XYZ coordinate system, direction cosines of the X axis are: Let: cos α 1 = cos β 2 cos γ 2 cos β 3 cos γ 3 , cos β 1 = − cos α 2 cos γ 2 cos α 3 cos γ 3 , cos γ 1 = cos α 2 cos β 2 cos α 3 cos β 3 . Then: The coordinate transformation between the two coordinate systems satisfies the following transforming relationships: x = x cos α1 + y cos α2 + z cos α3 + a y = x cos β1 + y cos β2 + z cos β3 + b z = x cos γ1 + y cos γ2 + z cos γ3 + c Writing the above equations in matrix form we get:

Parameter Optimization of Aspheric Surfaces Whose Vertices Locate at the Origin of the New Coordinate System
As shown in Figure 4, the vertex of an aspheric reflector locates at the origin O of the new coordinate system. The incident angle θ 1 of light ray passing through the origin O can be calculated with Equation (4). In the new coordinate system, unit vectors Q 1 (α 1 , β 1 , γ 1 ) along the incident light rays can be gotten from θ 1 , α 1 = −sinθ 1 , β 1 = −cosθ 1 , γ 1 = 0. Focus point coordinates in the new coordinate system can be solved using Equations (11) or (12). Therefore, in the X Y Z coordinate system, given Q 1 and focus point coordinates, according to the analysis in Sections 2.2 and 2.3, just like in the XYZ coordinate system, we can also find aspheric parameters C, a 2 , a 4 , a 6 , a 8 with the smallest fitness value in the X Y Z reference system. In this way, we will obtain aspheric parameters of another series of optimal aspheric reflectors constituting a concentrator system.

Seeking Parameters of the Optimal Aspheric Surface Whose Vertex Is the XYZ Coordinate Origin
As shown in Figure 4, let θ = 25°, and the coordinates of focus point in XYZ space be Pf (−186.5231, 400, 0). In the XOZ plane, we uniformly pick multiple incident points on the circle with center O and radius 10 cm, as shown in Figure 5a,b. Figure 5b shows that a square with a side length of 20 cm is divided into 2N × 2N uniform grids, and these grid nodes are used as incident points at which light rays (or their virtual extension) strike the XZ plane. These nodes will be reserved while their distance from the origin O is less than 10 cm, otherwise, they will be deleted. Thus, multiple incident light rays evenly hitting the aspheric reflector will be gotten. Taking N = 5, 2 2 186.5231 of all the reflected light rays can be obtained from Equation

Seeking Parameters of the Optimal Aspheric Surface Whose Vertex Is the XYZ Coordinate Origin
As shown in Figure 4, let θ = 25 • , and the coordinates of focus point in XYZ space be P f (−186.5231, 400, 0). In the XOZ plane, we uniformly pick multiple incident points on the circle with center O and radius 10 cm, as shown in Figure 5a,b. Figure 5b shows that a square with a side length of 20 cm is divided into 2N × 2N uniform grids, and these grid nodes are used as incident points at which light rays (or their virtual extension) strike the XZ plane. These nodes will be reserved while their distance from the origin O is less than 10 cm, otherwise, they will be deleted. Thus, multiple incident light rays evenly hitting the aspheric reflector will be gotten. Taking N = 5, all the reflected light rays can be obtained from Equation (3).
i R i is used as the fitness function of the PSO algorithm. In this design example, the number of parameter variables C, a 2 , a 4 , a 6 , a 8 controlling the shape of a reflector is the dimension of each particle's position vector, the searching range of each variable is based on a sensitivity analysis of effects of function parameters, and respectively set as

Seeking Parameters of the Optimal Aspheric Surface Whose Vertex Is at the Origin of the X′Y′Z′ Coordinate System
For the O′-X′Y′Z′ coordinate system, coordinates of point O' in the XYZ coordinate system are (a, b, c), and the positive direction cosines of X′, Y′, and Z′-axis in the XYZ coordinate system can be calculated with Equations (6), (8), and (9), respectively. The coordinates of focus point Pf in the X′Y′Z′ coordinate system can be worked out through Equation (10), and the angle θ1 between an incident light ray and Y′ -axis can be calculated using Equation (4). An incident ray vector ) , , ( in the X′Y′Z′ coordinate system is obtained from θ1, α1 = −sinθ1，β1 = −cosθ1，γ1 = 0. We can also use the PSO algorithm to seek aspheric parameters to suit our needs in the X′Y′Z′ coordinate system, just as in the XYZ coordinate system.  Figure 6 shows the focused results of an designed aspheric reflector. We can imagine from this figure that if there are many such aspheric reflectors whose vertex locates at different locations, and each aspheric reflector focuses the sunlight beam into the same small area, then the purpose of improving the concentration ratio will be achieved.

Seeking Parameters of the Optimal Aspheric Surface Whose Vertex Is at the Origin of the X Y Z Coordinate System
For the O -X Y Z coordinate system, coordinates of point O' in the XYZ coordinate system are (a, b, c), and the positive direction cosines of X , Y , and Z -axis in the XYZ coordinate system can be calculated with Equations (6), (8), and (9), respectively. The coordinates of focus point P f in the X Y Z coordinate system can be worked out through Equation (10), and the angle θ 1 between an incident light ray and Y -axis can be calculated using Equation (4). An incident ray vector Q 1 (α 1 , β 1 , γ 1 ) in the X Y Z coordinate system is obtained from θ 1 , α 1 = −sinθ 1 , β 1 = −cosθ 1 , γ 1 = 0. We can also use the PSO algorithm to seek aspheric parameters to suit our needs in the X Y Z coordinate system, just as in the XYZ coordinate system.
For example, the coordinates of point O in the XYZ coordinate system are (−100, 0, 40), namely a = −100, b = 0, c = 40, then we establish the X Y Z coordinate system according to Equations (6), (8), and (9). Thus, we can calculate from Equation (10) that P f (−186.5231, 400, 0) in the XYZ coordinate system is transformed into P f (−132.3867, 389.3071, 0) in the X Y Z coordinate system. From Equation (4), we have θ 1 = 18.7810 • . Therefore, in the X Y Z coordinate system, Q 1 (α 1 , β 1 , γ 1 ) is gotten from α 1 = −sin18.7810 • , β 1 = −cos18.7810 • , γ 1 = 0. The known Q 1 and P f , we can use the same method as in the XYZ coordinate system to find the aspheric parameters in the X Y Z coordinate system: C = −1/828.2705, a 2 = 34.6162, a 4 = −2.7765 × 10 −9 , a 6 = −1.8671 × 10 −11 , a 8 = −5.0820 × 10 −15 . The focused spot radius max (R i ) = 0.5580 cm. Figure 6 shows the focused results of an designed aspheric reflector. We can imagine from this figure that if there are many such aspheric reflectors whose vertex locates at different locations, and each aspheric reflector focuses the sunlight beam into the same small area, then the purpose of improving the concentration ratio will be achieved. Figure 7 shows the relative position of each aspheric reflector vertex in a concentrator system. First, all reflector vertices are located in the XOZ plane. Second, they all are at the grid nodes as shown in Figure 7. Watching along the Z direction, all reflectors are placed in 10 rows, each row has 21 reflectors from first to eighth row, the ninth row has 15 reflectors, and the tenth row has 7 reflectors; thus this concentrating system consists of a total of 190 reflectors. The distance between two vertices of adjacent aspheric reflectors in the same row is 20 cm.  Figure 7 shows the relative position of each aspheric reflector vertex in a concentrator system. First, all reflector vertices are located in the XOZ plane. Second, they all are at the grid nodes as shown in Figure 7. Watching along the Z direction, all reflectors are placed in 10 rows, each row has 21 reflectors from first to eighth row, the ninth row has 15 reflectors, and the tenth row has 7 reflectors; thus this concentrating system consists of a total of 190 reflectors. The distance between two vertices of adjacent aspheric reflectors in the same row is 20 cm. One by one, we use the above method to find aspheric parameters  Table 1 lists the optimization results of aspheric parameters of 21 reflectors in the eighth row, and the last column of Table 1 is the maximum spot radius being brought about by each aspheric reflector. Because of the limitation of length, the optimization results for 169 other groups of aspheric parameters are not listed.
This concentrator system has a total of 190 aspheric reflectors, the sunlight receiving area of each reflector in the plane vertical to the light beam is   Figure 7 shows the relative position of each aspheric reflector vertex in a concentrator system. First, all reflector vertices are located in the XOZ plane. Second, they all are at the grid nodes as shown in Figure 7. Watching along the Z direction, all reflectors are placed in 10 rows, each row has 21 reflectors from first to eighth row, the ninth row has 15 reflectors, and the tenth row has 7 reflectors; thus this concentrating system consists of a total of 190 reflectors. The distance between two vertices of adjacent aspheric reflectors in the same row is 20 cm. One by one, we use the above method to find aspheric parameters  Table 1 lists the optimization results of aspheric parameters of 21 reflectors in the eighth row, and the last column of Table 1 is the maximum spot radius being brought about by each aspheric reflector. Because of the limitation of length, the optimization results for 169 other groups of aspheric parameters are not listed.
This concentrator system has a total of 190 aspheric reflectors, the sunlight receiving area of each reflector in the plane vertical to the light beam is One by one, we use the above method to find aspheric parameters C, a 2 , a 4 , a 6 , a 8 of each reflector. Table 1 lists the optimization results of aspheric parameters of 21 reflectors in the eighth row, and the last column of Table 1 is the maximum spot radius being brought about by each aspheric reflector. Because of the limitation of length, the optimization results for 169 other groups of aspheric parameters are not listed.
This concentrator system has a total of 190 aspheric reflectors, the sunlight receiving area of each reflector in the plane vertical to the light beam is π[10 cos (25 0 )] 2 cm 2 as shown in Figure 5a, the focused spot radius is no more than 1 cm. So the concentration ratio of the entire concentrator systems is calculated to be 15606:1. Of course, the number of aspheric reflectors is still appropriately added to thereby increase the concentration ratio of the concentrator systems.
In the above analysis, we take the sun as the point sun. In fact, solar rays are not strictly parallel to each other, but with a 2σ included angle as shown in Figure 8, and σ = 16 . In the XY plane, the theoretical diameter of the sun's image formed by reflection of light from any point on the reflector is calculated by the following equation, D = d cos θ = 2R tan 16 cos θ , here R is the distance between point P and the intersection P f of a reflected ray and the focal plane, and its value changes slightly with the change of reflection point P. R can be obtained by using the two-point distance formula, and the incident angle θ can be calculated from vectors α 2 and β 2 of a reflected ray at point P, cos θ = α 2 α 2 2 +β 2 2 . In the present embodiment, R ≈ 4.5 m, θ ≈ 24 • , then D ≈ 4.6 cm, therefore, there is 4.6 times the error between the actual spot and the spot obtained by using the point sun. Hence, the corrected concentration ratio is 3393. focused spot radius is no more than 1 cm. So the concentration ratio of the entire concentrator systems is calculated to be 15606:1. Of course, the number of aspheric reflectors is still appropriately added to thereby increase the concentration ratio of the concentrator systems. Figure 8. While the solar beam is actually a cone beam with a solid angle of 32′, the focus spot will be extended.
In the above analysis, we take the sun as the point sun. In fact, solar rays are not strictly parallel to each other, but with a 2σ included angle as shown in Figure 8, and σ = 16′. In the XY plane, the theoretical diameter of the sun's image formed by reflection of light from any point on the reflector is calculated by the following equation,   cos 6 1 tan 2 cos , here R is the distance between point P and the intersection Pf of a reflected ray and the focal plane, and its value changes slightly with the change of reflection point P. R can be obtained by using the two-point distance formula, and the incident angle θ can be calculated from vectors and of a reflected ray at point P, In the present embodiment, R ≈ 4.5 m, θ ≈ 24°, then D ≈ 4.6 cm, therefore, there is 4.6 times the error between the actual spot and the spot obtained by using the point sun. Hence, the corrected concentration ratio is 3393. Table 1. The parameters of aspheric reflectors whose vertices are located at different positions obtained with the particle swarm optimization algorithm. Pf (x′, y′) Figure 8. While the solar beam is actually a cone beam with a solid angle of 32 , the focus spot will be extended.  Figure 9a is the focusing simulation of three aspheric reflectors, and the three local coordinate systems for three reflectors and the global coordinate system XYZ are also shown in the figure. In the simulation process, we first calculate the coordinate values of each desired point in the local coordinate system, and then convert these local coordinate values to the global coordinate system using Equations (11) or (12).  Figure 9a is the focusing simulation of three aspheric reflectors, and the three local coordinate systems for three reflectors and the global coordinate system XYZ are also shown in the figure. In the simulation process, we first calculate the coordinate values of each desired point in the local coordinate system, and then convert these local coordinate values to the global coordinate system using Equations (11) or (12).  Figure 9b shows the optical path simulation of the concentrator system consisting of 190 aspheric reflectors. It is shown that sunlight beams irradiating on 190 aspheric reflectors are reflected and focused in the same small area, so the idea of focusing sun beams in section 2 of this paper can be realized.

Computer Simulations of Focusing Effect
It can also be seen from Figure 9a,b that such a concentrator system can have better off-axis focus than a rotation paraboloid. If the light receiving aperture of each reflector is reduced, then in the meantime the number of reflectors is increased; the smaller focused spot will be gotten. If necessary, the smaller focused spot can be turned into a parallel beam with a powerful energy, so that it can be further exploited.  Figure 9b shows the optical path simulation of the concentrator system consisting of 190 aspheric reflectors. It is shown that sunlight beams irradiating on 190 aspheric reflectors are reflected and focused in the same small area, so the idea of focusing sun beams in Section 2 of this paper can be realized.

Verifying the Correctness of Designed Datum with Zemax Software
It can also be seen from Figure 9a,b that such a concentrator system can have better off-axis focus than a rotation paraboloid. If the light receiving aperture of each reflector is reduced, then in the meantime the number of reflectors is increased; the smaller focused spot will be gotten. If necessary, the smaller focused spot can be turned into a parallel beam with a powerful energy, so that it can be further exploited.

Verifying the Correctness of Designed Datum with Zemax Software
The datum in the first row of Table 1 is selected. It can be known from this row that the vertex coordinate of the reflector in the old coordinate system (global coordinate system) is (−60, 0, 200), establishing a new coordinate system (local coordinate system) whose origin is the vertex of the reflector. In this local coordinate system, the coordinate of focus point is obtained from Equation (10) and is (−190.2248, 424.0550, 0). Designed aspheric parameters are C = −1/960.0, a 2 = 112.3796, a 4 = −3.8535 × 10 −8 , a 6 = 5.00 × 10 −11 , a 8 = −1.2403 × 10 −14 , respectively. The angle θ 1 between the incident solar rays and Y' axis of the local coordinate system can be calculated from Equation (4) and is 24.1603 • . Let the beam diameter be 10 cm. The aspheric parameters, the viewing angle θ 1 , and the aperture value (beam diameter) are inputted to Zemax software as shown in Figure 10a, and the optical path through the reflector is also shown in Figure 10a. Figure 10b is the optical path simulated with a self-compiled Matlab program. By comparing Figure 10a,b, the spot radius obtained in Matlab and Zemax are about 1 mm and 990 µm, respectively, which means that both design results and self-compiled programs are correct. In this way, we can verify the correctness of the design of each reflector, and the correctness of the whole system design is also verified.
Let the beam diameter be 10 cm. The aspheric parameters, the viewing angle θ1, and the aperture value (beam diameter) are inputted to Zemax software as shown in Figure 10a, and the optical path through the reflector is also shown in Figure 10a. Figure 10b is the optical path simulated with a selfcompiled Matlab program. By comparing Figure 10a,b, the spot radius obtained in Matlab and Zemax are about 1 mm and 990 μm, respectively, which means that both design results and self-compiled programs are correct. In this way, we can verify the correctness of the design of each reflector, and the correctness of the whole system design is also verified.

Conclusions
In this study, we chose 25° special incident angle (between the incident ray and Y-axis) and designed an off-axis focused solar concentrator system consisting of 190 aspheric reflectors with 10 cm in radius and 3 mm in thickness. These aspheric reflectors were arranged in 10 rows. The detailed arrangement is listed as follows: row 1 to 8 contained 21 reflectors, row 9 contained 15 reflectors, and row 10 contained 7 reflectors. The vertex of each reflector lied in the same plane and located in the area defined by the four coordinates (

Conclusions
In this study, we chose 25 • special incident angle (between the incident ray and Y-axis) and designed an off-axis focused solar concentrator system consisting of 190 aspheric reflectors with 10 cm in radius and 3 mm in thickness. These aspheric reflectors were arranged in 10 rows. The detailed arrangement is listed as follows: row 1 to 8 contained 21 reflectors, row 9 contained 15 reflectors, and row 10 contained 7 reflectors. The vertex of each reflector lied in the same plane and located in the area defined by the four coordinates ( The distance between two vertices of adjacent aspheric reflectors was 20 cm and the focusing point of each reflector was located on coordinate (−186.5231, 400, 0). The system had theoretical concentration ratio of 3393:1 and the focused spot radius was less than 1 cm. This type of concentrating solar power system can be used in the solar thermal utilization under high temperature, concentrating solar power (CSP) plants, and solar-pumped lasers. This system has an important scientific significance for the development of solar energy utilization technology and has an important practical significance for environmental protection. In our future work, we will further optimize our current systems with reflectors of different designs, parameters, and even different optimization algorithms [42][43][44].