Journal of Solar Energy
Volume 2015 (2015), Article ID 326536, 10 pages
http://dx.doi.org/10.1155/2015/326536

Research Article

Ray Tracing Study of Optical Characteristics of the Solar Image in the Receiver for a Thermal Solar Parabolic Dish Collector

Saša R. Pavlovic and Velimir P. Stefanovic

Faculty of Mechanical Engineering, Department of Energetics and Process Technique, University of Niš, Niš, Serbia

Received 15 June 2015; Revised 2 September 2015; Accepted 27 September 2015

Academic Editor: Santanu Bandyopadhyay

Copyright © 2015 Saša R. Pavlovic and Velimir P. Stefanovic. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This study presents the geometric aspects of the focal image for a solar parabolic concentrator (SPC) using the ray tracing technique to establish parameters that allow the designation of the most suitable geometry for coupling the SPC to absorber-receiver. The efficient conversion of solar radiation into heat at these temperature levels requires a use of concentrating solar collectors. In this paper detailed optical design of the solar parabolic dish concentrator is presented. The system has diameter  mm and focal distance  mm. The parabolic dish of the solar system consists of 11 curvilinear trapezoidal reflective petals. For the construction of the solar collectors, mild steel-sheet and square pipe were used as the shell support for the reflecting surfaces. This paper presents optical simulations of the parabolic solar concentrator unit using the ray tracing software TracePro. The total flux on the receiver and the distribution of irradiance for absorbing flux on center and periphery receiver are given. The goal of this paper is to present the optical design of a low-tech solar concentrator that can be used as a potentially low cost tool for laboratory scale research on the medium-temperature thermal processes, cooling, industrial processes, polygeneration systems, and so forth.

1. Introduction and Survey of Literature

This paper presents the numerical results of the optimization of the solar image in a receiver for a fixed absorber in a solar parabolic concentrator, which was a project supported by the Ministry of Education, Science and Technological Development of Republic of Serbia. The device which is used to transform solar energy to heat refers to a solar collector. Solar thermal collectors have been widely used to concentrate solar radiation and convert it into medium-high-temperature thermal processes. In addition, the list of possible alternative applications of this technology is growing, due to the problems of oil dependency and global warming. They can be designed as various devices including solar cooker [1], solar hydrogen production [2, 3], and Dish Stirling system of harvest electricity [4, 5]. The main types of concentrating collectors are parabolic dish, parabolic trough, power tower, a Fresnel collector with mirror or lens, and stationary concentrating collectors. The ideal optical configuration for the solar parabolic thermal concentrator is a parabolic mirror. The parabolic mirror is very expensive to fabricate and its cost escalates rapidly with increase of aperture area. The parabolic mirror can be designed with large number of elementary components known as reflecting petals or facets. Usually reflecting petals are made from glass and their thickness is from 0.7 to 1.0 mm. Traditionally, the optical analysis of radiation concentrators has been carried out by means of computer ray-trace programs. Recently, an interesting analytical solution for the optical performance of parabolic dish reflectors with flat receivers was presented by O’Neill and Hudson [6]. Their method for calculating the optical performance is fast and accurate but assumes that the radiation source is a uniform disk. Saleh Ali et al. [7] have presented a study that aims to develop a 3D static solar concentrator that can be used as a low cost and low energy substitute. Their goal was to design solar concentrator for the production of portable hot water in rural India. They used ray tracing software for evaluation of the optical performance of a static 3D Elliptical Hyperboloid Concentrator (EHC). Optimization of the concentrator profile and geometry is carried out to improve the overall performance of the system. Kaushika and Reddy [8] used satellite dish of 2.405 m in diameter with aluminium frame as a reflector to reduce the weight of the structure and the cost of the solar system. In their solar system the average temperature of water vapor was 300°C, when the absorber was placed at the focal point. The cost of their system was US$ 950. El Ouederni et al. [9] were testing parabolic concentrator of 2.2 m in diameter with reflecting coefficient 0.85. Average temperature in their system was 380°C. Rafeeu and Ab Kadir [10] have presented simple exercise in designing, building, and testing small laboratory scale parabolic concentrators. They made two dishes from acrylonitrile butadiene styrene and one from stainless steel. Three experimental models with various geometrical sizes and diameters were used to analyze the effect of geometry on a solar irradiation. Liu et al. [11] presented a procedure to design a facet concentrator for a laboratory scale research on medium-temperature thermal processes. The facet concentrator approximates a parabolic surface with a number of flat square facets supported by a parabolic frame and having two edges perpendicular to the concentrator axis. Authors [12] presented a physical and mathematical model of the new offset type parabolic concentrator and a numerical procedure for predicting its optical performances. Also the process of design and optical ray tracing analysis of a low cost solar concentrator for medium-temperature applications is presented. This study develops and applies a new mathematical model for estimating the intercept factor of the solar concentrator based on its geometry and optical behaviour. Pavlovic et al. [13] developed mathematical model of solar parabolic dish concentrator based on square flat facets applied to polygeneration system. The authors developed optimization algorithm for search of optimal geometric, optical, and cost parameters. They have applied Monte Carlo ray tracing methodology which is used for analysis of the optical performance of the concentrator and to identify the set of geometric concentrator parameters that allow for flux characteristics suitable for medium- and high-temperature applications in trigeneration and polygeneration systems. A 164-facet concentrator will deliver up to 8.15 kW of radiative power over 15 cm radius disk located in the focal plane. Their system had an average concentration ratio exceeding 100. Hasnat et al. [14] presented two prototype parabolic dishes. Dunn et al. [15] investigated experimental evaluation of ammonia receiver geometries with a 9 m² dish concentrator. The 20 m² dish is mirrored with around 2,000 flat mirror tile facets arranged in concentric rings on a parabolic fiber glass support structure. Size of mirror facets is from 5 cm to 10 cm. Johnston et al. [16] analyzed optical performance of spherical reflecting elements for use with parabolic dish concentrators. This concentrator consists of 54 triangular mirrors. The effective rim angle of the dish is 46°. The 54 units are composed of nine separate panel shapes, each of the shapes is duplicated six times. The focal length of system is 13.1 m.

The decision to make solar parabolic concentrator with 11 petals is based on large number of design concepts that are realized in the world. This concept already proved to be useful in solar techniques, especially in production of heat and electrical energy as well as in trigeneration and polygeneration systems.

The basic idea behind this research is to start with primary concept of solar parabolic concentrator which will generate from 10 to 25 kW in polygeneration systems. Only with employment of parabolic concentrating systems it is possible to obtain high temperatures in range from 200°C to 800°C and high optical and thermal efficiency of concentrating solar collectors.

2. Geometrical Model of the Solar Parabolic Concentrator and Receiver

The design of the solar parabolic thermal concentrator and operation are presented. Optical design is based on parabolic dish with 11 curvilinear trapezoidal petals. Solar dish concentrators are generally concentrators that concentrate solar energy in a small area known as focal point. Dimensions of reflecting surfaces in solar dish concentrator are determined by desired power at maximum levels of insolation and efficiency of collector conversion.

The ray tracing technique is implemented in a software tool that allows the modelling of the propagation of light in objects of different media. This modelling requires the creation of solid models, either by the same software or by any computer aided design (CAD) software. Once in the optical modelling software, portions of the rays of the light source propagate in the flow of light, in accordance with the properties assigned to the relevant objects, which may be absorption, reflection, transmission, fluorescence, and diffusion. The sources and components of the light rays, adhering to various performance criteria involving the system parameters, result in simulation of the spatial and angular distribution, uniformity, intensity, and spectral characteristics of the system. Mathematical representation of parabolic concentrator is paraboloid that can be represented as a surface obtained by rotating parabola around axis. Mathematical equations for the parabolic dish solar concentrator in Cartesian and cylindrical coordinate systems are defined as

where are and are coordinates in aperture plane and is distance from vertex measured along the line parallel with the paraboloid axis of symmetry; is focal length of paraboloid, that is, distance from the vertex to the focus along the paraboloid axis of symmetry. The relationship between the focal length and the diameter of parabolic dish is known as the relative aperture and it defines shape of the paraboloid and position of focal point. The shape of paraboloid can be also defined by rim angle . Usually paraboloids that are used in solar collectors have rim angles from 10 degrees up to 90 degrees. The relationship between the relative aperture and the rim angle is given by

The paraboloid with small rim angles has the focal point and receiver at large distance from the surface of concentrator. The paraboloid with rim angle smaller than 50° is used for cavity receivers while paraboloids with large rim angles are most appropriate for the external volumetric receivers (central receiver solar systems). = 0.59. = 0.59 is focus/diameter ratio—ratios of focal distance to diameter of aperture reflector. Increasing the ratio reduces the rim angle. Paraboloid, with marginal angle, is a little curved, and its focal point and the receiver must be far from the surface of the concentrator. Paraboloid with rim angle of less than 50° is used when reflective radiation passes into the cavity of the receiver, while the paraboloids with large edge angles are most suitable for external receivers.

The geometric concentration ratio can be defined as the area of the collector aperture divided by the surface area of the receiver and can be calculated by

The designed solar parabolic concentrator has geometric concentration ratio = 100.

2.1. Parameters Design of Solar Parabolic Concentrator

Mechanical design of the solar parabolic concentrator is done in 3D design software CATIA, Dassault Systemes, USA. Parabolic shape of solar concentrator is obtained by entering and coordinates for selected points. For calculation of necessary points that define parabola public domain software Parabola Calculator 2.0 [17] is used. The calculated coordinates ( and ) for designed parabola are shown in Table 1. The calculated values are performed for 22 points in the parabola curve. The equation for parabola is

Table 1: Coordinates of designed parabola.

(cm)
−190.0
−158.3
−126.6
−95.00
−63.33
−31.67
0.0

(cm)
39.9
27.7
17.78
9.97
4.43
1.12
0.0

(cm)
31.67
63.33
95.00
126.6
158.3
190.0
—

(cm)
1.12
4.43
9.97
17.78
27.7
39.9
—

Geometrical model of solar parabolic concentrator is parametrically designed from calculated coordinates and it is shown in Figure 1(a). Selected model of solar dish concentrator with 11 petals requires very precise definition of parameters during geometrical modelling of system. Results obtained by optical analysis of solar concentration system are very much dependent on the selected method of the CAD model generation. Figure 2 represents the rear side of solar parabolic dish collector and geometrical model of construction of paraboloid curve with 11 petals/liefs. The geometrical parameters from Mathematica is exported in CATIA software and have built 3D CAD model.

(a)

(b)

Figure 1: (a) Parabola Calculator 2.0 [17]. (b) Paraboloid generated from Mathematica Wolfram software.

Figure 2: Mirrored segmented parabolic dish concentrator.

Figures 3(a) and 3(b) show the profile of spiral absorber tube and whole assembly of solar parabolic dish collector with receiver in their focus.

(a)

(b)

Figure 3: 3D CAD model of thermal solar concentration system and sinusoidal profile of corrugated pipe of spiral absorber.

A truncated paraboloid of revolution (circular paraboloid) is obtained by rotating the parabola segment about its axis (Figure 4). Consider a concentrator consisting of 11 trapezoidal reflective petals of identical nonoverlapping trapezoidal segments. 3D model of trapezoidal reflective petal of solar parabolic concentrator is presented in Figure 4. The outer diameter of corrugated pipe is = 12.2 mm, inner diameter is = 9.3 mm, and thickness of wall pipes is = 0.25 mm. Figure 5 shows the paraboloidal lief/petal of solar paraboloidal concentrator system. The front side of the petals is reflective surface with coefficient of reflectance of 90% (mirrored specular reflectance surface).

Figure 4: Schematic of truncated parabola.

Figure 5: Trapezoidal reflective petal of solar parabolic concentrator.

is aperture radius of absorber, while is aperture radius of solar dish concentrator/reflector. The small hole is constructive approach.

Detailed design parameters of solar parabolic concentrator are given in Table 2.

Table 2: Design parameters of solar parabolic concentrator.

Parameters
Numerical value
Unit

Aperture radius
0.2
[m]

Radius of smaller hole
1.9
[m]

Gross collector area
11.82
[m²]

The cross section of the opening parabola
9.89
[m²]

A sheltered area of the concentrator
0.126
[m²]

The effective area of the concentrator
9.764
[m²]

Receiver diameter
0.40
[m]

Shape of receiver
Directly irradiated corrugated spiral thermal absorber pipe
—

Depth of the concentrator
0.399
[m]

Focal length
2.26
[m]

6
[°]

45.6
[°]

Receiver-absorber is placed in focal area where reflected radiation from solar concentrator is collected. In the process of designing parabolic solar concentrators one always seeks for the minimum size of the receiver. With small receiver size one can reduce heat losses as well as cost of whole system. Also small receiver size provides increase of absorbed flux on the surface of receiver. This is the way of obtaining greater efficiency in conversion of solar radiation to heat. In our system receiver-absorber is Archimedean spiral corrugated pipe type with diameter of 400 mm. It is shown in Figure 6.

Figure 6: Solar cavity receiver-spiral corrugated pipe type of solar absorber.

In this paper only optical properties of receiver are analyzed. In our further research we plan to model all necessary details of receiver’s geometry which are important for conversion of solar energy into heat of fluid that is used for transfer energy.

3. Numerical Simulations and Ray Tracing Studies to Determine the Optical Characteristics of the Solar Image

For optical ray tracing analysis of solar parabolic thermal concentrator software TracePro, Lamda Research Corporation, USA, is used. Simulations have been performed using the technique of ray tracing to describe the behaviour of parabolic dish solar concentrators and such simulations have been used to study the behaviour of such concentrators in nonimaging optics; however, there is no such study or simulation of ray tracing type to describe parabolic solar concentrators. In TracePro all material properties are assigned. 11 trapezoidal reflective petals are defined as standard mirrors with reflective coating. Reflection coefficient was 95%. Receiver was cylinder with diameter 400 mm placed on 2075 mm from vertex of parabolic dish (optimal focal distance from vertex of parabolic solar concentrator). Absorbing surface was defined as perfect absorber. After definition of geometry of solar parabolic concentrator radiation source was defined. Radiation source was circular with diameter the same as diameter of parabolic dish (3800 mm). Radiation source was placed 2000 mm from vertex of parabolic dish and had circular grid pattern for generating 119401 rays for Monte Carlo ray tracing. Spatial profile of generated rays was uniform and angular profile was solar radiation. Input parameter for optical analysis is solar irradiance 800 W/m². Experiential value for solar irradiation for town of Niš in Serbia is between 750 W/m² and 900 W/m². Optical system for solar parabolic concentrator with traced rays is given in Figure 7.

(a)

(b)

Figure 7: (a) Ray tracing model of solar parabolic concentrator with traced rays. (b) Ray tracing model of solar spiral heat pipe with housing, SolTrace optical software tool, NREL, USA.

Temperature of the spiral absorber pipe-stagnation surface temperature is about 680°C. Detailed heat transfer analysis of this solar spiral heat pipe absorber will be presented in my other future papers.

Optical analysis is done by generating and calculating Monte Carlo ray trace for 119401 ray. From all emitted rays only 103029 rays reached absorber surface of which 82% of emitted rays are absorbed on receiver. Calculated irradiance for absorbed rays on receiver is from  W/m² to  W/m². Total calculated flux on receiver was 6001.7 W. In Figure 8 total irradiance map for absorbed flux on receiver is shown. When the sunlight shines on the solar collector including the direct and scattered radiation, there are three conditions affecting the absorption properties of the solar collector: (1) direct solar radiation absorbed directly by the solar collector and the light energy after the specular reflection; () the light energy of the scatter solar radiation after the diffuse and specular reflection; () the light energy of the direct solar radiation after the diffuse reflection. The slope of incident light is different in different latitude and time.

Figure 8: Irradiance map for absorbed flux on spiral receiver (solar focal image).

Figure 9 represents path sun rays from the source to the concentration of absorbing surface corrugated profile of absorber tube receiver.

(a)

(b)

Figure 9: (a) Schematic diagram of optical path shinning on the sinusoidal corrugated pipe absorber. (b) Dual axis solar tracking system with slew drives and gear motor boxes.

We take city of Nis (Laboratory for Solar Thermal Engineering) as an example. Nis is located in southeast Serbia longitude 43.30° north latitude and 21.90° east longitude. Solar Parabolic Dish concentrator (Figure 9(b)) is tracking the sun with motorized slewing drive solar tracking systems in two axes: azimuth and elevation (dual axis solar tracking precision system).

The influences of absorption rate for the sinusoidal corrugated absorber plate by aspect ratio and the slope of incident light are plotted in Figure 10. From the optical point of view, the time of light reflection absorption increases with the aspect ratio, which causes the absorption rate of the absorber plate to increase. Particularly, the absorption rate is close to 1, when the aspect ratio tends to infinity.

(a)

(b)

Figure 10: (a) Ray tracing and radiation control volume grid (1000 elements) for part of corrugated pipe. (b) Spiral corrugated one coil heat absorber from total 13 coils.

Figure 10 right shows ray tracing model, radiation control volume grid (1000 elements) for part of corrugated pipe.

Figure 10 shows control volume grid (1000 elements) for part of corrugated pipe, small part from all spiral corrugated absorber; this figure presented incident solar flux in control volume—small control volume—3 wrinkles-folds.

According to the boundary conditions, the optical properties can be obtained after many simulation calculations, which is called the ray tracing method [11].

Volume radiation flux distribution (incident and absorbed flux at the wall surface of profiles), software TracePro generated this bmp. file.

Figure 11 shows distribution of incident and absorbed solar hear flux at the corrugated surfaces of spiral thermal absorber.

Figure 11: Volume radiation flux distribution (incident and absorbed flux).

Figure 12 shows scattered model for reflection of incidence angle for solar parabolic dish concentrator with 11 curvilinear trapezoidal reflective petals. In general, the reflectance of a surface increases with angle of incidence. But for a mirror, the reflectance is already very high when AOI = 0°. It cannot increase much more. This is why the peak BRDF does not change with AOI. However, the reflectance of a black surface increases with AOI, no matter how black the surface is.

Figure 12: ABg giving the BRDF as a function of incidence angle.

4. Conclusions

This paper presents optical analysis of the solar parabolic concentrator using the ray tracing software TracePro. One can see that results obtained from optical design of solar parabolic concentrator are satisfactory. Total flux in focal area is good. Irradiance distribution for absorbed flux is relatively uniform for small area for absorber. A detailed simulation and analysis were conducted to evaluate the absorption rate of sinusoidal corrugated absorber pipe. As a next step, various analysis and simulations of the model are considered. Among others are variation of number of petals, size of petals, and shape of petals. Optimization methodology is not the topic of this work. Here it is only shown that the application of the Monte Carlo method can be performed to determine optimal parameters for obtaining the maximum value of total solar flux density and average irradiance-solar radiant energy at the spiral corrugated absorber.

The significance of this work is primarily in finding the optimal parameters of the geometric and optical parameters in the concentrator and receiver/absorber of concentration radiation. The main objective of this work is to find the optimal parameters in order to develop a prototype of a solar parabolic concentrator aspects of polygeneration systems, production of heating, cooling, and electricity using solar energy in Southeastern Europe.

Nomenclature

Latin Symbols

: | Area, in [m]
: | Facets height, in [m]
: | Coordinate at the direction
: | Coordinate at the direction
: | Coordinate at the direction
: | Collector aperture surface [m²]
: | Receiver aperture surface [m²]
: | Diameter of receiver surface [m]
: | Average geometric concentration ratio [−]
: | Focal length [m]
: | Incident direct solar flux, in [Wm⁻²]
: | Total radiative intensity, in [Wm⁻² sr⁻¹]
: | Radiative power, in [kW]
: | Radiative flux, in [Wm⁻²]
: | Average radiative flux, in [Wm⁻²]
: | Radius, in [m]
: | Position vector [−]
: | Collection efficiency [−]
: | Incident angle, in
: | Rim angle, in
: | Azimuthal angle, in
: | Constant.

Greek Symbols

: | Facet tilt angle, rad; limb darkening parameter
: | Small positive number
: | Density of the fluid, in [kg m⁻³].

Subscripts

: | Circumsolar
: | Flat-facet concentrator
: | Corresponding to sun disk.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This paper is done within the research framework of research projects: III42006: Research and Development of Energy and Environmentally Highly Effective Polygeneration Systems Based on Renewable Energy Resources. Project is financed by Ministry of Education, Science and Technological Development of Republic of Serbia. The authors would like to acknowledge the help provided by the Lambda Research Corporation by allowing the use of the TracePro software for the Ph.D. research of Saša R. Pavlovic.

References

A. A. Badran, I. A. Yousef, N. K. Joudeh, R. A. Hamad, H. Halawa, and H. K. Hassouneh, “Portable solar cooker and water heater,” Energy Conversion and Management, vol. 51, no. 8, pp. 1605–1609, 2010.

A. S. Joshi, I. Dincer, and B. V. Reddy, “Solar hydrogen production: a comparative performance assessment,” International Journal of Hydrogen Energy, vol. 36, no. 17, pp. 11246–11257, 2011.

P. Furler, J. R. Scheffe, and A. Steinfeld, “Syngas production by simultaneous splitting of H₂O and CO₂ via ceria redox reactions in a high-temperature solar reactor,” Energy & Environmental Science, vol. 5, no. 3, pp. 6098–6103, 2012.

T. Mancini, P. Heller, B. Butler et al., “Dish-stirling systems: an overview of development and status,” Journal of Solar Energy Engineering, vol. 125, no. 2, pp. 135–151, 2003.

D. Mills, “Advances in solar thermal electricity technology,” Solar Energy, vol. 76, no. 1–3, pp. 19–31, 2004.

M. J. O'Neill and S. L. Hudson, “Optical analysis of paraboloidal solar concentrators,” in Proceedings of the Annual Meeting, U.S. Section of International Solar Energy Society, vol. 2.1, pp. 855–862, Denver, Colo, USA, August 1978.

I. M. Saleh Ali, T. S. O'Donovan, K. S. Reddy, and T. K. Mallick, “An optical analysis of a static 3-D solar concentrator,” Solar Energy, vol. 88, pp. 57–70, 2013.

N. D. Kaushika and K. S. Reddy, “Performance of a low cost solar paraboloidal dish steam generating system,” Energy Conversion & Management, vol. 41, no. 7, pp. 713–726, 2000.

A. R. El Ouederni, M. Ben Salah, F. Askri, and F. Aloui, “Experimental study of a parabolic solar concentrator,” Revue des Energies Renouvelables, vol. 12, no. 3, pp. 395–404, 2009.

Y. Rafeeu and M. Z. A. Ab Kadir, “Thermal performance of parabolic concentrators under Malaysian environment: a case study,” Renewable and Sustainable Energy Reviews, vol. 16, no. 6, pp. 3826–3835, 2012.

Z. Liu, J. Lapp, and W. Lipiński, “Optical design of a flat-facet solar concentrator,” Solar Energy, vol. 86, no. 6, pp. 1962–1966, 2012.

S. Pavlovic, D. Vasiljevic, V. Stefanovic, Z. Stamenkovic, and S. Ayed, “Optical model and numerical simulation of the new offset type parabolic concentrator with two types of solar receivers,” Facta Universitatis, Series: Mechanical Engineering, vol. 13, no. 2, pp. 169–180, UDC 535.2, 2015.

S. R. Pavlovic, V. P. Stefanović, and S. H. Suljković, “Optical modeling of a solar dish thermal concentrator based on square flat facets,” Thermal Science, vol. 18, no. 3, pp. 989–998, 2014.

A. Hasnat, P. Ahmed, M. Rahman, and K. A. Khan, “Numerical analysis for thermal design of a paraboloidal solar concentrating collector,” International Journal of Natural Sciences, vol. 1, no. 3, pp. 68–74, 2011.

R. Dunn, K. Lovegrove, G. Burgess, and J. Pye, “An experimental study of ammonia receiver geometries for dish concentrators,” Journal of Solar Energy Engineering—Transactions of the ASME, vol. 134, no. 4, Article ID 041007, 2012.

G. Johnston, K. Lovegrove, and A. Luzzi, “Optical performance of spherical reflecting elements for use with paraboloidal dish concentrators,” Solar Energy, vol. 74, no. 2, pp. 133–140, 2003.

http://mscir.tripod.com/parabola/