International Journal of Engineering and Technical Research (IJETR)

ISSN: 2321-0869, Volume-2, Issue-12, December 2014

Mathematical modeling of heat transfer processes in

energy-saving solar energy systems
Morteza Piradl, Farzad Choobdar Rahim

mobility. As assumptions considered further mathematical

Abstract Interest in renewable energy sources (RES) has model will be assumed that the inner surface of the collector is
been steadily increasing throughout the world each year. It is a "black body" and does not reflect the solar radiation. The
not so much due to possible depletion of fossil fuels (coal, oil, gas),
upper wall of the manifold is transparent, made of glass. The
but hopes for the environmentally sound and sustainable
development of mankind in the future. Increasingly there is heating medium is passed through the collector due to the fan
widespread use of solar radiation. running on the photoelectric converter.

Index Terms solar energy, collector, temperature

distribution, heat transfer

I. INTRODUCTION
Solar energy system is the most promising and common
energy source for heat [1,2]. According to [3] it is currently
used in the world more than 180 million solar panels that
provide heat to various customers. The most common of these Fig. 1. Schematic diagram of the solar collector with active
technologies in China (59%), in second place - Europe (14%). circulation of coolant: 1 - blowing fan; 2 -clear glass; 3 - the
186 large firms produce solar collectors in 41 countries. air flow within the reservoir; 4 - - black coating absorbing
Characterized by high values of the average daily radiation in power metal coherency; 5 - insulating layer
the winter the amount of incoming solar energy is reduced by
the latitudinal location of the installation several times [3,4]. One way to reduce convective heat loss through the
The efficiency of solar collectors depends on the design transparent glass is the use of a cellular structure, located
features and climatic factors [5]. The most common solar air above the absorbing surface [7]. In this recirculation
collectors are used for heating buildings and rooms for installation is not performed, the coolant (air) is forced
various purposes, as in agriculture. [2,3,6] Despite the large through the manifold once, which simplifies the design of the
number of scientific publications on the use of solar energy, apparatus.
research local characteristics of reservoirs, taking into For the formulation of the computational problem and heat
account the actual characteristics of their operation, such as a transfer equation follows we assume that the air - it's viscous
higher degree of turbulence of the coolant at the inlet, and incompressible medium. In order to construct a
different wall temperature and so neglected.. mathematical model primarily formulated boundary value
problem and the boundary conditions. In addition discusses
the main heat balance equation.
II. MATHEMATICAL MODELLING
Complete the equation describing heat transfer in the
The role of mathematical modeling in research and reservoir with the free convection of the form:
development of this kind is very high. This is due to the need v 0 (1)
for deeper penetration into the essence of the objects and aims
v 1
of reducing the cost and time of development. To create a (v ) v p
simulation of the mathematical model the authors have t (2)
considered a flat solar collector with active circulation of g (T T 0) v
coolant (Figure 1). The main structural elements are the air
T 1
duct without absorber, whose side facing the sun blackens and (v ) T T q v (3)
the other - Insulate. Collector frame is made of durable and t c
lightweight metal, which reduce its weight and improve

(4)

Manuscript received December 03, 2014.
Morteza Piradl, Department of Mechanical Engineering, Payame Noor w r 1 -t ont nu t qu t on v 0 2 -
University, Islamic Republic of Iran, +98-9144475944. qu t on o mot on 3 - t n r qu t on u,v,w) - the
Farzad Choobdar Rahim, Nuclear Science and Technology Research air velocity vector in the channel; T - temperature in the duct;
Institute, Atomic Energy Organization of Iran (AEOI), Tehran, Iran,
+98-9141495975. p - pressure caused by the movement of air; g - acceleration of

39 www.erpublication.org
 Mathematical modeling of heat transfer processes in energy-saving solar energy systems

r v t v tor - expansion coefficient of air; - the calculated area take the inner surface of the collector for the
k n m t v s os t - air density; qv - heat flux density, c - approximation of the area and construct an orthogonal grid,
sp to r 1.005 J/k , - o nt o t rm l instead of the area further consider the set of nodes formed by
on u t v t o r 0.241 /m2 )- operator of the the intersection of the lines parallel to the coordinate axes.
form: Fig. 3 shows an orthogonal grid and a plurality of nodes for
calculating the temperature distribution in the reservoir.
( ) v v v ;
xx yy zz

Taking assumptions, consider a two-dimensional problem,

and we believe the process to be stationary, hence, the time
derivative is zero, the air velocity in the vertical v = 0. The
initial portion of the channel is small and the velocity profile
in the channel is known from the classical theory [8] (the
initial distribution - parabola) (see Fig. 2).

Fig. 3. Grid elements with dummy nodes

To set the boundary conditions G1 and G2 with the

accuracy with which the original equation is approximated,
we supplement the grid with dummy nodes, placing them at
the top and bottom borders of the G1 and G2 at a distance hy.
The desired function of the temperature distribution T (x,
Fig. 2. The velocity profile in the air manifold y) is approximated grid function that is a set of values of the
function at the grid points. The magnitude of the size hx and hy
In this case, the initial system of equations is reduced to a respectively can be calculated from the following relations:
boundary value problem: yj y2 j
T T 2 H
Hy ;
L
hx ; ui , j 6u 0 ( );
u 2 m n hy hy
x y
T (0, y) T0 Approximation of derivatives on the four-point look at the
template, the original equation can be rewritten in a Cartesian
T
(T T ) coordinate system. The original equation:
y When y = 0 T 2T
T ( x,0)
u 2 ;
q x y
y When y = h, where Approximated by the first and second derivatives with
2 difference quotients can be carried out as follows:
y y
u 6u 0 ( 2 ) Ti Ti , j Ti 1, j
h h q1=q2 u ;
x hx
2Ti 1 Ti 1 Ti Ti Ti 1 Ti , j 1 2Ti , j Ti 1
Where q is thermal flow directed from the top wall to the .
bottom. y 2 hy hy hy hy
2

This system is solved by one of the numerical methods. To Replace the partial derivatives of and at
simulate this process was chosen finite difference method [9, node 1 through the difference relations, as a result we obtain
10]. Differential equation as a result of transformations the equation in general form for i = 1. In the future, we form
replaced by an equivalent value in the finite difference these equations for all nodes in which unknown values of the
solution of which is to implement simple algebric operations. unknown function.
The final outcome of the decision given by the expression in
Ti , j Ti 1, j Ti , j 1 2Ti , j Ti , j 1
2
yj yj
which the value of the "future" of the potential (temperature) 6u0 ( 2
) a 2
;
at a given point, (node) is determined by the "real" potential hy hy hx hy
and the "real" potential adjacent nodal points. Repeatability
i=1n; j=1m;
of the same operations in the calculation of temperature fields
a great convenience for the application of modern computer
The boundary conditions can be represented as:
technology, so the efficiency of the many times increases. For

40 www.erpublication.org
 International Journal of Engineering and Technical Research (IJETR)
ISSN: 2321-0869, Volume-2, Issue-12, December 2014

Ti ,1 Ti , 1 - For the node (i, m).

(Ti ,0 T ); The resulting system is solved by the Gauss method.To
2h y solve this system, a program that allows us to calculate the
- The node (i,0)
temperature distribution in the collector channel. To do this,
Ti ,m Ti ,m1
q1 ; the program stores the value of T1,j (j = 0,1,2 ..., m), then move
2h y to the next value of i, and so on until i = m. The result is an
- The node (i,m).
array of temperatures, which is the last column - is the
temperature at the outlet of the reservoir. When i = 1, Ti-1,j=
As a result, we obtain the following system of equations for
T0,j (j = 0,1,2 ... m), T0,j (j = 0,1,2 ... m), the values set in other
i = 1:
i as Ti-1,j (j = 0,1,2 ... m), value is the temperature obtained in
the previous step.
Ti ,m Ti ,m1
q1 ; In addition the program takes into account the influence of
2h y parameters such as the dimensions and layout of the reservoir
capacity of the solar insolation, the slope of the collector,
ambient air temperature, and air velocity in the channel.
y0 y0 T1,0 T0,0 T 2T1,0 T1, 1
2

6u0 ( 2 ) a 1,1 2
; Table 1. The results of calculation of temperatures for steps
hy hy hx hy
0-3
i = 0, j = 0;

y1 y1 T1,1 T0,1 T 2T1,1 T1,0

2
6u0 ( 2 ) a 1, 2 2
;
hy hy hx hy
i =1, j = 1;

y1 y 2 T1, 2 T0, 2 T 2T1, 2 T1,1

2
6u0 ( 2 ) a 1,3 2
;
hy hy hx hy
i = 1, j = 2;

..
y m y m T1,m T0,m 2T1,m T1,m1
2
T
6u0 ( 2 ) a 1,m1 2
;
hy hy hx hy
i = 1, j = m;

Ti ,1 Ti , 1
(Ti ,0 T);
2h y

The totality of the above equations written in the general

form for the node (i, j) form a discrete mathematical model of As seen from Table 1, the temperature is increased
the original problem. throughout the reservoir. With the help of the developed
A further step is to compile the simulation algorithm and program can perform calculations for different initial
writing software to solve the resulting lower system of linear conditions and obtain the values of the temperature at
algebraic equations On the basis of the above discrete different points of the manifold. Varying the number of steps
mathematical model of the original problem to make the can increase the accuracy of calculations. For the most
algorithm of formation of the coefficient matrix and realistic data about the temperature in the collector must be
right-hand side column. set as the initial value of the measured quantity of solar
For each i, starting with i=1, a system of linear algebraic radiation on the lower surface of the collector, after passing
qu t ons ons st n o m + 1 qu t on j 0 1 2 m o t the radiation through a transparent surface of the upper wall.
form: The results of calculation of the temperature distribution
Ti , j Ti 1, j Ti , j 1 2Ti , j Ti , j 1 along the length of the reservoir can also be presented
ui , j a 2
; graphically (see Fig. 4). Red line indicates the temperature in
hy hy the first step of the program (on the inlet manifold), blue
- For the node (i, 0); dotted line - in the middle collector orange dashed line - to the
Ti ,m Ti ,m1 final step of the program, i.e. output from the manifold.
q1 ;
2h y

41 www.erpublication.org
 Mathematical modeling of heat transfer processes in energy-saving solar energy systems

Fig. 4. The temperature distribution along the length of the

reservoir under normal conditions

III. CONCLUSION
As a result of a simulation of a mathematical model
obtained data on the temperature distribution within the
reservoir, as well as to determine the temperature
characteristics of the coolant at the outlet of the reservoir. In
addition, the model allows us to obtain data on the influence
of external conditions, such as outdoor temperature, the
power of solar insolation, as well as the angle of the collector
to the horizon and the percentage of heat loss through the
walls of the reservoir.

REFERENCES
[1] 1. J. A. Duffy, and W.A. Beckman. "Fundamentals of solar thermal
power" Recently: Intellect. - 2013 - 886 c.
[2] 2. F. Kreith, and J.F. Kreider "Principle of solar engineering"
Washington: Hemisphere Pub. Corp. 1978 778 p.
[3] 3. V.A. Butuzov "Solar heating in Russia: status and regional
characteristics" Energosovet, 2011, No. 5 (18).
[4] 4. V.N. Petrenko, and N.V. Mokrova. "Develop a system of hot water
supply from renewable energy sources", "Engineering Gazette Don",
2013, No. 2.
[5] 5. A. Hematian, Ya. Ajabshirchi, and A.A. Bakhtiari "Experimental
analysis of flat plate solar air collector efficiency," Indian Journal of
Science and Technology. 2012, V. 5. P. 3183-3187.
[6] 6. F. Chabane, N. Moummia, and S. "Benramachec Experimental
analysis on thermal performance of a solar air collector with
longitudinal fins" Algeria J. of Power Technologies, 2013, V. 93(1).
[7] 7. M.I. Romanov, and V.V. Sherstyukov "Energy-efficient method of
using the excess heat of the solar collector" Engineering Gazette Don.
2012, No 4-2.
[8] 8. A.V. Klimenko and V.M. "Theoretical, basics of heat, Thermal
experiment": Handbook. Publisher MEI, 2001 - 564 p.
[9] 9. Y.A. Bahvalov. "Mathematical modeling: Proc. Manual for schools"
Novocherkassk: SRSTU (NPI), 2010. - 142c.
[10] 10. A.N. Tikhonov, A.A. "Samara Equations of mathematical physics"
M.: Publisher State University, 2004 - 799 p.

42 www.erpublication.org