Ray Tracing Study to Determine the Characteristics

of the Solar Image in the Receiver for a Thermal

Solar Concentration System
* * **
Saša R. Pavlovic PhD student. , Velimir P. Stefanović Dr Sci, Predrag Rajković Dr Sci, Emina P. Petrovic
PhD student,*** Sadoon Ayed MSc****
*Nis. Serbia, Faculty of Mechanical Engineering, Department for Energetics and Process technique, University in Niš, Serbia
** Nis. Serbia, Faculty of Mechanical Engineering, Department for Mathematics and Natural Sciences University in Niš, Serbia
*** Nis. Serbia, Faculty of Mechanical Engineering, Mechatronic and Control Systems, University in Niš, Serbia
****
University of Technology, Department of Mechanical Engineering, Baghdad ,Iraq,

e-mail : [email protected], [email protected], [email protected], [email protected],

[email protected]

Abstract-- This study presents the geometric aspects of the focal electricity [4,5].The main types of concentrating collectors are:
image for a solar parabolic concentrator (SPC) using the ray parabolic dish, parabolic trough, power tower, Fresnel
tracing technique to establish parameters that allow the designationcollector with mirror or lens and stationary concentrating
of the most suitable geometry for coupling the SPC to absorber – collectors. The ideal optical configuration for the solar
receiver. The efficient conversion of solar radiation in heat at these
temperature levels requires a use of concentrating solar collectors.
parabolic thermal concentrator is a parabolic mirror. The
In this paper detailed optical design of the solar parabolic dish parabolic mirror is very expensive to fabricate and it's cost
concentrator is presented. The system has diameter D = 3800 mm escalating rapidly with increase of aperture area. The
and focal distance f = 2260mm. The parabolic dish of the solar parabolic mirror can be designed with large number of
system consists from 11 curvilinear trapezoidal reflective petals. elementary components known as reflecting petals or facets.
For the construction of the solar collectors, mild steel- sheet and Usually reflecting petals are made from glass and their
square pipe were used as the shell support for the reflecting thickness is from 0.7 to 1.0 mm. Traditionally, the optical
surfaces. This paper presents optical simulations of the parabolic analysis of radiation concentrators has been carried out by
solar concentrator unit using the ray-tracing software TracePro. means of computer ray-trace programs. Recently, an
The total flux on receiver and the distribution of irradiance for
absorbed flux on center and periphery receiver are given. The goal
interesting analytical solution for the optical performance of
of this paper is to present optical design of a low-tech solar parabolic dish reflectors with flat receivers was presented by
concentrator that can be used as a potentially low-cost tool for O'Neill and Hudson [6]. Their method for calculating the
laboratory-scale research on the medium-temperature thermal optical performance is fast and accurate but assumes that the
radiation source is a uniform disk. Imhamed M. Saleh Ali et al.
processes, cooling, industrial processes, polygeneration systems etc.
[7] have presented study that aims to develop a 3-D static
Keywords: ray tracing analysis, optical modelling, thermal solar solar concentrator that can be used as low cost and low energy
concentration technology, corrugated Archimedean spiral pipe solar substitute.Their goal were to design solar concentrator for
absorber production of portable hot water in rural India. Their used ray
tracing software for evaluation of the optical performance of a
I. INTRODUCTION AND SURVEY OF static 3-D Elliptical Hyperboloid Concentrator (EHC).
LITERATURE Optimization of the concentrator profile and geometry is
carried out to improve the overall performance of system.
This paper presents the numerical results of the optimisation Kaushika and Reddy [8] used satellite dish of 2.405 m in
of the solar image in a receiver for a fixed absorber in a solar diameter with aluminium frame as a reflector to reduce the
parabolic concentrator, which was a project supported by the weight of the structure and cost of the solar system. In their
Ministry of Education, Science and Technological solar system the average temperature of water vapor was
Development of Republic of Serbia.The device which is used 3000C, when the absorber was placed at the focal point. Cost
to transform solar energy to heat is reffered to a solar collector. of their system was US$ 950. El Ouederni et al. [9] was
Solar thermal collectors have been widely used to concentrate testing parabolic concentrator of 2.2 m in diameter with
solar radiation and convert it into medium-high temperature reflecting coefficient 0.85. Average temperature in their
thermal processes. In addition, the list of possible alternative system was 3800C.Y. Rafeeu and M.Z.Z.AbKadir [10] have
applications of this technology is growing, due to the presented simple exercise in designing, building and testing
problems of oil dependency and global warming. They can be small laboratory scale parabolic concentrators. They made two
designed as various devices including solar cooker [1], solar dishes from acrylonitrile butadiene styrene and one from
hydrogen production [2,3] and Dish Stirling system of harvest stainless steel. Three experimental models with various
geometrical sizes and diameters were used to analyze the of symmetry. The relationship between the focal length and
effect of geometry on a solar irradiation. Zhiqiang Liu et al. the diameter of parabolic dish is know as the relative aperture
[11] presented a procedure to design a facet concentrator for a and it defines shape of the paraboloid and position of focal
laboratory-scale research on medium – temperature thermal point. The shape of paraboloid can be also defined by rim
processes. The facet concentrator approximates a parabolic angle ψrim. Usually paraboloids that are used in solar
surface with a number of flat square facets supported by a collectors have rim angles from 10 degrees up to 90 degrees.
parabolic frame and having two edges perpendicular to the The relationship between the relative aperture and the rim
concentrator axis. angle is given by:
The decision to make solar parabolic concentrator with 11 1
petals is based on large number of design concepts that are (2)
4 tan /2
realized in the world. This concept already proved useful in The paraboloid with small rim angles have the focal point and
solar techniques, especially in production of heat and receiver at large distance from the surface of concentrator.
electrical energy as well as in trigeneration and polygeneration The paraboloid with rim angle smaller than 500 (
systems. 45,6 is used for cavity receivers while paraboloids with large
The basic idea behind this research is to start with primary rim angles are most appropriate for the external volumetric
concept of solar parabolic concentrator which will generate receivers (central receiver solar systems). f/D = 0.59
from 10 to 25 kW in polygeneration systems. Only with The geometric concentration ratio can be defined as the
employment of parabolic concentrating systems it is possible area of the collector aperture Aapp divided by the surface area
to obtain high temperatures in range from 200C to 800C and of the receiver Arecand can be calculated by eq.3.
high optical and thermal efficiency of concentrating solar
 
1
collectors. CRg  sin 2  a  Ac Ar 1  Aapp / Arec (3)
II. GEOMETRICAL MODEL OF THE SOLAR The designed solar parabolic concentrator has geometric
PARABOLIC CONCENTRATOR AND RECEIVER concentration ratio CRg = 100.

The design of the solar parabolic thermal concentrator, and A. Parameters Design of Solar Parabolic Concentrator
operation are presented. Optical design is based on parabolic Mechanical design of the solar parabolic concentrator is
dish with 11 curvilinear trapezoidal petals. Solar dish done in 3D design software CATIA, Dassault Systems, USA.
concentrators are generally concentrators that concentrate Parabolic shape of solar concentrator is obtained by entering x
solar energy in a small area known as focal point. Dimensions and y coordinates for selected points. For calculation of
of reflecting surfaces in solar dish concentrator are determined necessary points that define parabola public domain software
by desired power at maximum levels of insolation and Parabola Calculator 2.0 [11] is used. The calculated
efficiency of collector conversion. coordinates (x and y) for designed parabola are shown in
The ray tracing technique is implemented in a software Table I. The calculated values is performed for 22 point in the
tool that allows themodelling of the propagation of light in parabola curve. The equation for parabola is:
objects of different media. This modelling requires the
creation of solid models, either by the same software or by (4)
any computer aided design (CAD) software. Once in the The coefficients describing this parabola are:
optical modelling software, portions of the rays of the light a = 1.10526 x 10 -4
source propagate in the flow of light, in accordance with the b=0
properties assigned to the relevant objects, which may be c = -2.428275 x 10-14
absorption, reflection, transmission, fluorescence, and
diffusion. The sources and components of the light rays, 1.10526 10 2.428275 10 (5)
adhering to various performance criteria involving the system
parameters, result in simulation of the spatial and angular
TABLE I COORDINATES OF DESIGNED PARABOLA
distribution, uniformity, intensity, and spectral characteristics
of the system. Mathematical representation of parabolic X(cm) -190.0 -158.3 -126.6 -95.00 -63.33 -31.67 0.0
concentrator is paraboloid that can be represented as a surface Y(cm) 39.9 27.7 17.78 9.97 4.43 1.12 0.0
obtained by rotating parabola around axis. Mathematical X(cm) 31.67 63.33 95.00 126.6 158.3 190.0 -
equations for the parabolic dish solar concentrator in Cartesian Y(cm) 1.12 4.43 9.97 17.78 27.7 39.9 -
and cylindrical coordinate systems is defined as:
4 /4 (1)
where are x and y are coordinates in aperture plane and z is
distance from vertex measured along the line parallel with the
paraboloid axis of symetry; f is focal length of paraboloid i.e.
distance from the vertex to the focus along the aparboloid axis
 Fig.4 Schematic of truncated Fig.5 Trapezoidal reflective petal
parabola of solar parabolic concentrator

Detailed design parameters of solar parabolic concentrator

Fig.1. Parabola Calculator 2.0 [11]
is given in Table II.
Geometrical model of solar parabolic concentrator is TABLE II DESIGN PARAMETERS OF SOLAR PARABOLIC CONCENTRATOR
parametrically designed from calculated coordinates and it is
shown on Fig. 1. Selected model of solar dish concentrator Parameters Numerical value Unit
Aperture radius R1 0.2 [m]
with 11 petals requires very precise definition of parameters
Radius of smaller hole R2 1.9 [m]
during geometrical modelling of system. Results obtained by Gross collector area Agross 11.82 [m2]
optical analysis of solar concentration system are very much The cross section of the opening
9.89 [m2]
dependent on the selected method of the CAD model parabola Aproj
generation. A sheltered area of the concentrator
0.126 [m2]
Ashadow
The effective area of the
9,764 [m2]
concentrator Aef=AprojAshadow
Receiver diameter 0.40 [m]
Corrugated spiral
Shape of receiver
pipe – circular disc -
Depth of the concentrator 0.399 [m]
Focal lenght 2.26 [m]
ψ1 6 [ 0]
ψ2 45,6 [ 0]

Receiver - absorber is placed in focal area where reflected

radiation from solar concentrator is collected. In the process of
Fig.2 CAD model of solar parabolic concentrator with dimensions and rear designing parabolic solar concentrators one always seek for
view of solar concentration system 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 provide 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 on Fig. 6.

Fig.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
[Fig.4]. Consider a concentrator consisting of 11 trapezoidal Fig.6. Solar cavity receiver – spiral corrugated pipe type of solar absorber
reflective petals of identical non-overlapping trapezoidal
segments. 3D model of trapezoidal reflective petal of solar In this paper only optical properties of receiver are analysed.
parabolic concentrator is presented on Fig. 4. The outer In our further research we plan to model all necessary details
diameter of corrugated pipe is De = 12.2 mm, inner diameter of receiver's geometry which are important for conversion of
Di = 9.3 mm and thickness of wall pipes is s = 0.25 mm. solar energy into heat of fluid that is used for transfer energy.
 III. NUMERICAL SIMULATIONS AND RAY (1) direct solar radiation absorbed directly by the solar
TRACING STUDIES TO DETERMINE THE collector and the light energy after the specular reflection; (2)
OPTICAL CHARACTERISTICS OF THE the light energy of the scatter solar radiation after the diffuse
SOLAR IMAGE and specular reflection; (3) the light energy of the direct solar
radiation after the diffuse reflection. The slope of incident
For optical ray tracing analysis of solar parabolic thermal light is different in different latitude and time, so we deal it
concentrator software TracePro, Lamda Research Corporation, with integral processing in the range of incident angle.
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 non – imaging
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 Fig.8. Irradiance map for absorbed flux on spiral receiver (solar focal image)
source was circular with diameter 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/m2. Experiential value for
solar irradiation for town of Niš in Serbia is between 750
W/m2and 900 W/m2. Optical system for solar parabolic
concentrator with traced rays is given in Fig. 7. Fig.9. Schematic diagram of optical path shinning on the sinusoidal
corrugated pipe absorber

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. According to the declination angle δ (–23.45°on
winter solstice), sunset hour angle ω, altitude of the sun h, the
installation angle of 45° and the North − south placement of
solar parabolic collector. In order to calculate the angle of
incidence, the sun angle is defined as follows:
δ = 23.45sinε, ɛ = 2πd/365, ϕ = 21.900, ω = 15t
sinh = sinφsinδ + cosφcosδcosω
sin , cos cos sin cosh ,
where d is the number day of 365, t is the time in hours.
TABLE III CALCULATION RESULTS IN NIS
Fig.7. Optical system of solar parabolic concentrator with traced rays
Sunset hour Altitude of the Angle of
Optical analysis is done by generating and calculating Monte Time
angle [ω] sun h incident i
Carlo ray trace for 119401 ray. From all emitted rays only 9:00 -450 180 550
103029 rays reached absorber surface which is 82% rays of 10:00 -300 250 450
11:00 -150 300 370
emitted rays are absorbed on receiver. Calculated irradiance
12:00 00 320 340
for absorbed rays on receiver is from 1.74·10-8 W/m2 to 62100 13:00 150 300 370
W/m2. Total calculated flux on receiver was 7800W. On Fig. 14:00 300 250 450
8 is shown total irradiance map for absorbed flux on receiver. 15:00 450 180 550
When the sunlight shines on the solar collector including the
direct and scattered radiation, there are three conditions The influences of absorption rate for the sinusoidal corrugates
affecting on the absorption properties of the solar collector: absorber plate by aspect ratio and the slope of incident light
are plotted in Fig. 10. From the optical point of view, the is good. Irradiance distribution for absorbed flux is relatively
times of light reflection absorption increases with the aspect uniform for small area for absorber. A detailed simulation
ratio, which causes the absorption rate of the absprber plate and analysis was conducted to evaluate the absorption rate of
increasing. Specially, the absorption rate is close to 1, when sinusoidal_corrugated absorber pipe. Lights with arbitrary
the aspect ratio tends to infinity. incident angle shinning on the sinusoidal_corrugates absorber
plate can occur second times of reflection, part of the light
have third times and more reflection. In future development of
optimization method is planned. This optimization method
will make it possible to find optimal geometrical and optical
parameters of the various types of solar parabolic dish
concentrators as well as geometrical, optical and thermal
parameters of different types of solar cavity receivers -
absorbers.

ACKNOWLEDGMENT
Fig.10. Ray tracing and radiation control volume grid (1000 elements)
for part of corrugated pipe This paper is done within the research framework of
research projects: III42006 – Research and development of
According to the boundary conditions, the optical properties energy and environmentally highly effective polygeneration
can be obtained after many simulation calculations, which is systems based on renewable energy resources. Project is
called the ray tracing method [11]. financed by Ministry of Education, Science and Technological
Development of Republic of Serbia. Authors acknowledge to
Lambda Research Corporation for allowing to use software
TracePro for PhD Thesis research of SašaPavlović.

REFERENCES
[1] Badran, A.A.Yousef,I.A.,et al.,“Portable solar cooker and water
heater”, Energy Conversion and Management 51 (8), 1605–1609, 2010.
[2] Joshi,A.S.,Dincer,I.,Reddy, B.V., “Solar hydrogen production: a
comparative performance assessment”, International Journal of
Fig. 11. Volume radiation flux distribution (incident and absorbed flux) Hydrogen Energy 36(17),11246–11257, 2011.
[3] Furler, P.,Scheffe,J.R.,Steinfeld,A., “Syngas production by
simultaneous splitting of H2O and CO2 via ceria redox reactions in
Figure 12. shows scattered model for reflection of ahigh-temperature solar reactor”, Energy & Environmental Science
incidence angle for solar parabolic dish concentrator with 11 5(3),6098–6103,2012.
curvilinear trapezoidal reflective petals. In general, the [4] Mancini,T.,Heller,P.,Butler,B.,Osborn,B.,Schiel,W.,Goldberg,V.,Buck,
R.,Diver,R.,Andraka,C.,Moreno,J..“Dish–Stirling systems:an over
reflectance of a surface increases with angle of incidence. But view of development and status“,Journal of Solar Energy Engineering-
for a mirror, the reflectance is already very high when the AOI Transactions of the ASME 125(2),135–151, 2003.
= 0 deg. It can’t increase much more. This is why the peak [5] Mills,D. “Advances in solar thermal electricity technology”, Solar
Energy 76 (1–3),19–31, 2004.
BRDF doesn’t change with AOI. However, the reflectance of
[6] O'Neill, M. J. and Hudson, S. L. “Optical Analysis of Paraboloidal
a black surface increases with AOI, no matter how black the Solar Concentrators." Proceedings 1978 Annual Meeting, U.S. Section
surface. of Int. Solar Energy Society, August 1978; Denver, CO. Vol. 2.1: p.855
Scatter Distribution (BRDF)
[7] Imhamed M. Saleh Ali, Tadhg S. O’Donovan , K.S. Reddy, Tapas K.
for various Angles of Incidence (sample mirror) Mallick, “An optical analysis of a static 3-D solar concentrator, ” Solar
1000 Energy 88, 57-70, 2013.
100
[8] N.D.Kaushika, K.S.Reddy,“Performance of low cost solar paraboloidal
AOI =10 deg
dish steam generating,” Energy Conversion & Management 41, 713-
10 726, 2000.
BRDF (1/sr)

AOI =20 deg

1
AOI =30 deg [9] A.R.El. Ouederni, M. Ben Salah, F.Askri and F. Aloui, “Experimental
AOI =40 deg study of a parabolic solar concentrator, Revue des Renouvelables ,Vol.
0.1 AOI =50 deg
12, 395-404, 2009.
AOI =60 deg
0.01 AOI =70 deg
[10] Rafeeu, M.Z.A.AbKadir, “Thermal performance of parabolic
concentrators under Malysian environment,”: A case study, Renewable
0.001 and Sustainable Energy Reviews, 16, pp. 3826-3835, 2012.
0 10 20 30 40 50 60 70 80 90
[11] Zhiqiang Liu, Justin Lapp, Wojciech Lipinski, “Optical design of a flat
θ (from surface normal)
solar concentrator, Solar Energy 86, pp. 1962-1966, 2012.
Fig. 12. ABg giving the BRDF as a function of incidence angle [12] http://mscir.tripod.com/parabola/
[13] Lubos Mitas, Quantum Monte Carlo, J. Current Opin ion Solid State
CONCLUSIONS Mater. Sci., 1997, vol. 6, no. 2, pp. 696–700.

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