Elliptical-Shaped Fresnel Lens Design through Gaussian Source Distribution

: A novel three-dimensional elliptical-shaped Fresnel lens (ESFL) analytical model is presented to evaluate and maximize the solar energy concentration of Fresnel-lens-based solar concentrators. AutoCAD, Zemax and Ansys software were used for the ESFL design, performance evaluation and temperature calculation, respectively. Contrary to the previous modeling processes, based on the edge-ray principle with an acceptance half-angle of ± 0.27 ◦ as the key deﬁning parameter, the present model uses, instead, a Gaussian distribution to deﬁne the solar source in Zemax. The results were validated through the numerical analysis of published experimental data from a ﬂat Fresnel lens. An in-depth study of the inﬂuence of several ESFL factors, such as focal length, arch height and aspect ratio, on its output performance is carried out. Moreover, the evaluation of the ESFL output performance as a function of the number/size of the grooves is also analyzed. Compared to the typical 1–16 grooves per millimeter reported in the previous literature, this mathematical parametric modeling allowed a substantial reduction in grooves/mm to 0.3–0.4, which may enable an easy mass production of ESFL. The concentrated solar distribution of the optimal ESFL conﬁguration was then compared to that of the best ﬂat Fresnel lens conﬁguration, under the same focusing conditions. Due to the elliptical shape of the lens, the chromatic aberration effect was largely reduced, resulting in higher concentrated solar ﬂux and temperature. Over 2360 K and 1360 K maximum temperatures were found for ESFL and ﬂat Fresnel lenses, respectively, demonstrating the great potential of the three-dimensional curved-shaped Fresnel lens on renewable solar energy applications that require high concentrations of solar ﬂuxes and temperatures.


Introduction
Optical concentration provides strong cost leverage for photovoltaic cells [1]. High concentrated photovoltaic technology uses relatively inexpensive optics, such as mirrors and lenses, to concentrate sunlight from a broad area into a much smaller area of active semiconductor cell and converts sunlight directly to electricity [2,3]. The Fresnel lens has been widely used in the concentrated photovoltaic field, with the advantages of simple structure, light weight, low cost, easy processing, etc. [4]. However, it has a limited concentration ratio due to its strong chromatic aberration.
To overcome this issue, non-flat Fresnel lenses have been widely studied by many researchers, proposing new modeling methods and configurations with the aim of increasing the solar energy concentration. The shaped Fresnel lens conception was initiated in 1977 by Cosby [5]. One year later, a patent was filed by O'Neill [6], while Kritchman published his finding [7]. Since then, many other researchers throughout the world have been proposing their own non-flat Fresnel lens models and modeling processes [8][9][10][11][12][13]. Currently, the most Fresnel lens with published experimental data [36], whose results were in accordance with those obtained through the proposed model. An in-depth study of the influence of the ESFL parameters, such as focal length, arch height and aspect ratio on its output performance, in terms of concentrated solar flux, optical efficiency and full width at half maximum (FWHM), is also carried out. Furthermore, to the best of our knowledge, this is the first time that the ESFL concentration efficacy is evaluated as a function of the number of grooves, with optimal 0.4-0.3 groove numbers per millimeter being found. This study reveals that the number of grooves necessary to maximize the concentrated solar flux could be significantly reduced in relation to that reported by the literature and market, which may enhance the cost efficiency of the manufacturing process of Fresnel lens solar concentrators. From the abovementioned studies, the best ESFL design is found. The optimal concentrated solar flux value within its focal cone is then analyzed and compared to that of a flat Fresnel lens, which proved the effectiveness of the ESFL in substantially reducing the chromatic aberration and, consequently, on maximizing the solar flux. The temperature analysis of both concentrators is also performed by Ansys. The maximum temperatures of 2360 K and 1360 K are attained from the ESFL and the flat Fresnel lens, respectively, demonstrating once more the potential of the ESFL in many solar energy research and applications that require high solar flux and temperature.

Modeling of an Elliptical-Shaped Fresnel Lens
The ESFL model has an elliptical-shaped arch with a smooth surface facing the Sun and grooves on the opposite side, as shown in Figure 1. This guarantees the easier cleaning of the lens in case of dirt, without damaging the grooves.

Analytical Method
The design of the ESFL followed a generic ellipse Equation (1), represented in Figure 2.
As shown in Figure 2, the conception of the elliptical arch depends on the radius or aperture of the concentrator (r), its focal height (h f ) and the height of the arch (h l ). The angle ω is the aperture angle or rim angle of the concentrator, with h f and r defining its size. The ellipse's major axis "a" and minor axis "b" can be both calculated with those variables through Equation (2). The minor axis "b" is the combination of h f and h l and, once it is found, the major axis "a" can be acquired.
The ESFL can be considered as a set of well-defined and positioned prisms, as shown in Figure 3.  The prisms were modeled through Equation (3), based on Snell's law, starting with the outmost one, as indicated in Figure 3a. Each prism is defined by three facets: the input, output and back facets, represented in Figure 3b. The input facet is defined by the segment between point A and point B. α represents the angle between the vertical line (solar ray path), with the normal line of the input facet AB. The solar ray is refracted as it hits AB. Depending on the refractive index of the medium (n λ ), the ray is deviated by an angle γ. It is then shifted further as its leaves the output facet of the prism, making an angle of ω in relation to the vertical axis. The output facet is either defined by the segment AC or AD depending on the pitch angle. In this paper, the refractive index n λ of poly (methyl methacrylate) (PMMA) material (n λ = 1.492) was considered.
The number of prisms/grooves (N) was calculated through Equation (4), where δω is the groove division angle, which represents a small angular segment of the elliptical Fresnel lens.
Equations (5) and (6) give the distances of point A and B to the focal point (O, as the origin point), respectively, as shown in Figure 3a. The size of the input facet AB depends on δω. A larger δω means a larger input facet, thus a larger prism. n is the current groove number.
Knowing the distances AO n and BO n and the angle ω, then all coordinates can be found with basic trigonometry. The Cartesian coordinates of points A and B are then obtained through Equations (7) and (8), respectively.
x A, n = AO n sin(ω − (n − 1) δω) y A, n = AO n cos(ω − (n − 1) δω) , x B, n = BO n sin(ω − (n) δω) The length of the facet AB at any given prism number is calculated by Equation (9): The input angle α at any given prism number is given by Equation (10): The angle γ is obtained by isolating it from the remaining parameters of Snell's Equation (3), as shown in Equation (11): The calculation of angle β n is acquired by directly applying Snell's law into Equation (10), resulting in Equation (12): The exiting angle ω n on the nth prism can be calculated through Equation (13): The inclusion of the pitch angle (θ p ) changes the length of the output facet from AC to AD. Equation (14) represents the calculation of the coordinates of point D considering an isolated prism, where point A is defined as the origin. This facilitates the calculation of the real coordinates of point D. The output facet length (AD) is then found by Equation (15).
x D,n = x A,n − AD n cos(β n ) y D,n = y A,n − AD n sin(β n ) , Finally, a single prism can be drafted with Equations (7), (8) and (16), forming a triangle of points ABD. The grooves of the ESFL's output surface were the first components to be modeled by drawing all the facets from coordinates A n and D n , starting from the outmost prism (n = 1) to the innermost prism N. This chaining process followed a sequence of A 1 D 1 A 2 D 2 . . . A N D N . The input surface of the ESFL was then drawn from the last prism (Nth prism) to the first one, i.e., N to n = 1, with the Y coordinate offset by an increment of d t , the thickness of the concentrator from point A N . Hence, the coordinates of point A N take the form of (x A,N , y A,N + d t ). The output drawing is shown in Figure 4, and the procedural chaining process is summarized in a flowchart, as presented in Figure 5.   Figure 6 shows the sequence of the numerical simulation process. The sequence was the same for each ESFL configuration, but each output was unique in terms of concentrated solar flux, optical efficiency and focal size at FWHM.

ESFL Modeling
Each individual coordinate of the ESFL was calculated with the Mathematica programming language. Then, the coordinates of all the relevant Cartesian points were exported into AutoCAD and linked as straight lines by the "polyline" command. The enclosed object was then revolved into a solid with the highest polygon count possible by adjusting the "facetres" command to its maximum. The object was then exported to Zemax as an IGES file. Both AutoCAD and Zemax shared the same coordinates, as shown in Figure 7.

Solar Source Modeling
An "ellipse source" with the same size of the concentrator was used to simulate the solar rays. The total power attributed is given by the product of the concentrator's collection area and the irradiance, i.e., a 1 m diameter concentrator would have a source power of about 785 W at 1000 W/m 2 irradiance. Figure 8 shows both the global and the direct reference spectra ASTM (American Society for Testing and Materials) G173 at Air Mass 1.5 (AM1.5) [37]. In Zemax, 21 wavelengths were selected as the solar spectrum data, each normalized as a function of its weight. The weight determines the intensity/power of each solar ray, which is dependent on the power attributed to the emitting source. The Gaussian distribution of the solar source was defined by Zemax through the G x and G y parameters, i.e., the Gaussian distribution parameters in X and Y axes, respectively. These parameters are based on the solar irradiance measured by Vittitoe and Biggs [35].
The root-mean-square width (δ RMS ) adopted by Vittitoe and Biggs, as shown in Equation (17), describes the Sun shape as a function of the irradiance (I).
Equation (18) represents the Gaussian distribution of the solar rays from the object "elliptical source" in Zemax. At 1000 W/m 2 solar irradiance, the Gaussian source would have a value of 35,276.6 at both G x and G y axes.
Alternatively, the Gaussian distribution of G x and G y can be manually adjusted in Zemax based on experimental measurements. For example, in a numerical study of a threedimensional ring-array concentrator [28], where its output performance was compared to that of the medium sized solar furnace (MSSF) concentrator of PROMES-CNRS (Procédés Matériaux et Energie Solaire-Centre National de la Recherche Scientifique) [38], G x = G y = 36,000 was considered [28]. In this case, the solar distribution at the focal zone of the MSSF concentrator had the same characteristics as those described by [38][39][40].
The terrestrial solar half-angle (θ E ) can be calculated by determining the effective size of the solar-terrestrial image (d E ) at a distance (L), as represented in Equation (19).
θ E can also be found in Zemax by using a detector at a certain distance from a small solar light source, as shown in Figure 9a. The length L is set as 10 m from the source of 0.002 mm diameter, and the size of the detector (d) as 150 × 150 mm 2 with a precision of 1001 × 1001 pixels. The calculation of θ E is independent of the distance L since d E is adjusted by the inverse-square law.  Figure 9b shows the solar distribution of the solar source obtained from both Equations (17) and (18), considering an irradiance of 1000 W/m 2 with G x = G y = 35,276. For an accurate calculation of the terrestrial solar angle, 95% of the total d E distribution was considered and θ E = 0.52 • was calculated. It is important to note that the acceptance half-angle (θ a ) of 0.27 • can be found at the effective size (d a ) by considering 56% of the focal Gaussian distribution.

Output Solar Distribution at the Focal Zone of the ESFL
To analyze the solar distribution at the focal zone of each ESFL configuration, over 60 million rays were employed per simulation. The output data were obtained from an absorbing detector with dimensions 20 × 20 mm 2 and resolution 150 × 150 pixels, positioned at the mathematical origin. The obtained data contained the solar distribution characteristics in terms of concentrated solar flux, optical efficiency (the total number of rays that strike the detector over the total emitted rays from the source) and FWHM. Figure 10 shows the comparison of the three-dimensional focal distributions from a Gaussian source with G x = G y = 35,276 and a collimated source with G x = G y = 0 (or infinite). In this case, an arbitrary ESFL model with 1 m diameter (D), h f = 400 mm, h l = 400 mm, d t = 3 mm, and δω = 0.20 • (256 grooves), with θ p = 12 • , was used. The focal shape formed by the Gaussian source ( Figure 9a) has a wide normal distribution with 10.8 mm FWHM and a concentrated solar flux of 5.0 W/mm 2 , while the collimated source ( Figure 9b) has a needle-shaped distribution with 0.3 mm FWHM, resulting in a peak concentrated solar flux of more than 100 W/mm 2 . The collimated source offers the highest concentration intensity and its FWHM is nearly as tight as a laser beam. However, this is not possible to obtain with incoherent light provided from the Sun.

Comparative Study of the ESFL Output Performance with the Measured Output
Performance of a Fresnel Lens Figure 11 shows the simulated layout and the focal output of the ESFL and a flat Fresnel lens [36] at the same focusing conditions. Ferriere et al. used a flat Fresnel lens of 900 mm diameter with 757 mm focal length, 20 grooves/cm (1800 grooves), and 31.7 mm thickness [36]. The measurement was conducted in PROMES-CNRS, where the irradiance can reach 1000 W/m 2 [36,38,41]. The remaining Fresnel lens parameters were adjusted to achieve the same output performance of the publication, by assuming a conic constant of −1.95 and a modest pitch angle of 2 • . By using the solar source Gaussian distribution from [35], the simulated focal output of the flat Fresnel lens in Zemax matched well with the measured data [36]. An ESFL with same collection size and focal length, and h l = 300 mm, δω = 0.0017 • (1807 grooves), d t = 3 mm and θ p = 12 • , was numerically simulated and compared. Figure 11a,b present the cross-sectional view of the light rays from five concentric annulus solar sources (with same area and power) passing through the ESFL and the Ferriere Fresnel lens, respectively. Area 1 represents the solar rays from the outmost concentric annulus source, while Area 5 corresponds to the solar rays of the innermost circular source. Figure 11c,d show the contribution of each source on the concentrated solar flux of the ESFL and Ferriere Fresnel lens, respectively, as well as the combined concentrated solar fluxes distribution. As shown in Figure 11b, it is noticeable that the solar rays at the external annulus area of the flat Fresnel lens (Area 1 and 2) are barely detected, leading to solar fluxes close to zero (Figure 11d). The curved shaped of the ESFL overcomes this problem, allowing the solar rays from the external annulus area to be more efficiently focused, as demonstrated in Figure 11a,c. Since the external annulus areas collect a majority of the incoming solar power, the advantage of the ESFL concentrator in attaining higher solar flux becomes evident. Chromatic aberration is also significantly reduced, leading finally to a higher concentrated solar flux of 4.5 W/mm 2 (Figure 11c) compared to that of the flat Fresnel lens with 2.6 W/mm 2 (Figure 11d) [36]. In addition to the elliptical shape of the Fresnel lens, there are also important factors that can contribute to the better performance of this concentrator, such as the aspect ratio and the number/size of grooves. The influence of these factors on the ESFL output performance is addressed in Section 3.  Table 1 shows a summary of some predominant and documented analytical and measured data of other types of concentrators. It is important to note that the experimental data from the existing concentrators should be used as reference in the design and simulation of other concentrators, in order to obtain an accurate and fair evaluation of their performances.

Comparative Study of the ESFL Output Performance with Other Concentrators
The parabolic mirror of the MSSF had a measured peak concentration of 16 W/mm 2 with 2000 mm diameter, 850 mm focal length and a focal size of 10 mm FWHM [39]. This parabolic concentrator has currently the highest known concentrated solar flux. However, based on analytical predictions, some authors have claimed higher concentrations with focal width similar to that presented in Figure 10b [36]. The present ESFL, at the same conditions of the flat Fresnel lens from [36], had a better focal solar flux of 4.5 W/mm 2 .

ESFL Configurations with Fixed Total Height of 700 mm
As any lens, the shape of the ESFL concentrator has direct influence on its output performance. Figure 12 shows some of the possible configurations of the ESFL with a combined h f and h l of 700 mm total height. The ESFL has a similar configuration of a flat Fresnel lens at h f = 700 mm, h l = 0 mm. As h l increases, the lens becomes more curved, and the aspect ratio becomes larger. The aspect ratio (AR) of the ESFL is given by the quotient of h l and D (diameter of the lens): Figure 13 shows the variation of the concentrated solar flux as a function of the aspect ratio of ESFL and h l , with 700 mm total height. All the ESFL configurations had D = 1 m, δ ω = 0.28 • , θ p = 12 • and d t = 3 mm. The concentrated solar flux starts with 3.9 W/m 2 at h l = 0 mm, which closely resembles that of a flat Fresnel lens. It gradually increases with h l , at which point the Fresnel lens is shaped into a parabola. The concentrated solar flux is maintained over 5.2 W/mm 2 from h l of 500 mm to 700 mm, i.e., with aspect ratio from 0.603 to 0.702.   Figure 14 shows the concentrated solar flux, optical efficiency and FWHM at various combinations of h f and h l with ESFLs of D = 1 m, δ ω = 0.28 • , θ p = 12 • and d t = 3 mm. h f ranges from 200 mm to 600 mm and h l ranges from 50 mm to at least 500 mm.

Concentrated Solar Flux, Optical Efficiency and FWHM as a Function of the Aspect Ratio and Several Combinations of h f + h l
Each h f has its own peak concentrated solar flux at different h l (Figure 14a). For example, in the ESFL configuration with h f = 200 mm, the concentrated solar flux varies from 4.2 W/mm 2 at h l = 50 mm to 5.2 W/mm 2 at h l = 500 mm, above which it starts to decrease. The highest concentrated solar flux was attained by the lowest h f = 200 mm at h l = 500 mm (aspect ratio of 0.603). The increase in h f resulted in a lower concentrated solar flux, whose peak value shifted to a lower h l , decreasing the aspect ratio of the concentrator.  The optical efficiency tends to stagnate at larger h l , as shown in Figure 14b. The highest optical efficiency of about 60% was attained at the higher h f of 600 mm. The FWHM enlarges in a near linear fashion with the increase in h l, and it also expands with the increase in h f , as shown in Figure 14c.

Influence of Size/Number of Grooves on the ESFL Output Performance
The influence of different sizes of δω on the ESFL output performance was analyzed for the most favorable combinations of h f and h l that achieved the highest concentrated solar fluxes (Figure 14a). The variation of δω changes the size of each groove/prism within the ESFL concentrator, mainly changing the size of the input (A n B n ) and the output (A n D n ) facets, as calculated in Equations (9) and (15). A range of δω, from 0.05 • to 1 • , was used for the conception of the ESFL model. Figure 15a shows the number of grooves per δω at different h f . The number of grooves drops exponentially independently of the h f with the increase in δω. A gap in the number of grooves between h f is observed at a smaller δω, while it converges closer into the same number of grooves at a larger δω. Figure 15b shows the concentrated solar flux as a function of δω. Regardless of h f , the maximum concentrated solar flux was found at δω = 0.3 • , which is equivalent to 210 to 132 grooves for h f 200 mm to 600 mm, respectively. As shown in Figure 15c, the increase in δω lowers the optical efficiency from 66% to 56% and 53% to 50% at h f of 600 mm and 200 mm, respectively. The FWHM shows a minimum within a certain δω range, as shown in Figure 15d. For example, at h f = 600 mm, more than 10.7 mm FWHM was found at δω = 0.35 • , while at h f = 200 mm, the minimum FWHM of about 8 mm was found at δω = 0.2 • . In all cases, the largest FWHM is located at δω = 1 • . The number of grooves has a direct influence on the focal characteristics of the concentrator, as demonstrated in Figure 16. On the one hand, with the increasing number of grooves (lower δω), each prism becomes smaller, which reduces the refraction space within the prism. Consequently, above a certain number of grooves, the number of effective rays that could be refracted onto the targeted focal position is diminished, causing the decrease in the concentrated solar flux, as illustrated in Figure 16b. On the other hand, the reduction in the number of grooves (higher δω) leads to an increment in the size of the prism, hence broadening the final output focal shape, as shown in Figure 16d. Therefore, the concentration solar flux also diminishes with the decrease in the number of grooves (Figure 16b). Maximum concentrated solar flux was numerically found at δω = 0.30 • for all the ESFLs, regardless the h f , as demonstrated in both Figures 15b and 16b. However, the optimum number of grooves varies with h f . Consequently, for δω = 0.30 • , the optimum number of grooves varied from 227 grooves with h f = 200 mm to 132 grooves with h f = 600 mm, corresponding to 0.42 grooves/mm and 0.26 grooves/mm, respectively. This variation is even more pronounced with smaller δω (higher number of grooves). The optimum groove number per size is by far lower than that of other Fresnel lenses (Table 1), such as the shaped Fresnel lens with 4.3 groove/mm at D = 460 mm [22], flat Fresnel lens with 2 grooves/mm at D = 900 mm [36], 2 grooves/mm at D = 889 mm [45], and a range of 15.7 grooves/mm to 1 grooves/mm at various sizes of Fresnel lenses found at the market, such as the Fresneltech [25]. The higher performance of the ESFL with larger prisms and, consequently, lower number of grooves is attributed to the fact that the design process was not under the influence of the common modeling processes based on the edge-ray principle, but based on the solar Gaussian distribution. The reduction in the number of grooves and subsequent increase in the prism size could be favorable for the manufacturer in neglecting the limit losses due to manufacturing inaccuracies, such as blunt tips and deformed grooves forming.   Figures 17a and 18a, were represented by a detector volume with 100 × 100 × 100 voxels in vacuum, where each voxel accumulates and stores the energy data from each ray that passes through it, with the associated wavelength and power. Figure 17b shows the top view focal distribution of the ESFL in the detector positioned at the origin (Z = 0 mm), with a concentrated solar flux of 5.08 W/mm 2 , while Figure 17c shows the focal distribution at Z = −5 mm in relation to the origin, with a maximum concentrated solar flux of 5.48 W/mm 2 . Figure 18b shows the focal distribution of the flat Fresnel lens at Z = 0 mm, with a concentrated solar flux of 0.75 W/mm 2 . At Z = +30 mm in relation to the origin, the maximum concentrated solar flux of only 1.86 W/mm 2 was reached, as shown in Figure 18c.

Optimal Focal Position Analysis of Both the ESFL and a Flat Fresnel Lens
For both the ESFL and the flat Fresnel lens, the focal distributions at the origin (Z = 0 mm) do not correspond to the positions where the concentrated solar flux is maximum, as shown in Figures 17a and 18a. This phenomenon occurs due to the chromatic aberration, which is more abundant in the flat Fresnel lens. It is important to note that a parabolic concentrator (a concentrator with no chromatic aberration effect) has its maximum concentration exactly at the origin point [28].

Temperature Analysis of Both the ESFL and the Flat Fresnel Lens
Zemax non-sequential ray-tracing and Ansys finite element analyses were both used to evaluate the temperature of both the ESFL and the flat Fresnel lens at the Z position with the highest concentration solar flux found in Section 4.1. Ansys allows thermo-optical calculations that deal with complex geometric shape and boundary conditions, enabling the approximation of variables in a volume or surface element that changes across the matrix [46][47][48].
In Zemax, a square absorbing black body (emissivity ε = 1 [49]) detector of 20 × 20 mm 2 and 150 × 150 pixels was used to calculate the concentrated solar flux at the optimal focal positions of the ESFL and the flat Fresnel lens. The matrixial data were then exported to the Ansys workbench through the "External data" component and loaded as a heat flux source.
In Ansys 2021 finite element analysis, a graphite disk receiver of 20 mm diameter and 5 mm thickness was used to obtain the temperature of both concentrators. The graphite of 2250 kg/m 3 constant density, 24 W/m K thermal conductivity and 709 J/Kg K specific heat were chosen from the Ansys internal material library. The disk receiver was divided by the tetrahedrons meshing method with a sizing element of 0.4 mm. It contained approximately 1600 elements, which were enough for FEA calculations, with good approximation. The boundary conditions set for the convection applied onto the disk were the same as the natural stagnant air convection, representing a heat transfer coefficient of 5 × 10 −6 W/mm 2 /K. The radiation exchange between surfaces was restricted by a gray-diffused surface and the emissivity for the graphite disk surface was confined to ε = 0.85. A room temperature of 295.15 K was considered.
The respective temperatures of the ESFL and the flat Fresnel lens are shown in Figure 19. Both temperatures were generated from the focal distributions of Figures 17c and 18b, respectively. The ESFL attained maximum and minimum temperatures of 2362 K and 1945 K, respectively, which were 1.73 and 1.60 times more than that of the Fresnel lens with the maximum temperature of 1363 K and minimum temperature of 1217 K, respectively.

Conclusions
To overcome the aberration chromatic issue of Fresnel lens concentrators, a novel parametric model of a three-dimensional elliptical-shaped Fresnel lens was provided and analyzed. The modeling process took into account the solar Gaussian distribution based on the measured parameters by Vittitoe and Biggs [35], instead of the classical edge-ray principle method and the solar acceptance half-angle of 0.27 • . The design was performed by CAD software and then imported into Zemax for numerical calculations. The accuracy of this model was confirmed by the numerical analysis of the output performance of a flat Fresnel, whose results matched well with the experimental data [36]. The ESFL output performance was compared to these results, under the same focusing conditions, revealing the advantage of the ESFL in focusing more efficiently the solar rays from the external annulus area, at the external annulus collection area, where a significant part of the incoming solar power is collected, resulting in a significant increase in the concentrated solar flux.
The study of different combinations of ESFL focal length and arch height, and their influence on the output performance was carried out. The highest peak solar flux of about 5.2 W/m 2 was attained for the ESFL with a shorter focal length of h f = 200 mm and high arch height of h l = 500 mm, resulting in a large aspect ratio of 0.603. As h f increased, the solar concentration decreased, but maximum solar fluxes were attained with a smaller h l , hence lowering the aspect ratio.
The optimal concentrated solar flux value within the focal cone of the best ESFL was then studied and compared to that of a flat Fresnel lens at its best output performance configuration. This study demonstrated the effectiveness of the ESFL in reducing the chromatic aberration, leading to a significant enhancement of the concentrated solar flux and temperature, as compared to that of the flat Fresnel lens.
In addition, the present work also provided a comprehensive study of the influence of the number of grooves and size on the ESFL performance. It was found that the optimal number of groves per millimeter could be substantially reduced in relation to that reported by the previous literature and market. This could greatly facilitate the manufacturing process of Fresnel lens concentrators while increasing its solar concentration capacity, revealing the promising potential of ESFL in many solar energy research and applications.