1. Introduction
The paper puts forward an in-house mathematical model of a parabolic trough collector (PTC). The collector is fitted with a sun-tracking system with two rotation axes to increase the energy gain. The solar collector comprises evacuated tubes with a triple-wall design of the absorber and two U-tubes inside. It should be highlighted that commonly used PTC collectors are equipped with a single U-tube. Thus, the case considered seems to be an innovative construction.
Solar technologies are becoming increasingly popular these days. Sunlight as a source of energy is used to generate direct current electricity in photovoltaic (PV) panels and to produce heat in solar collectors. There are also hybrid photovoltaic panels with flat solar collectors coupled with the PV panel to increase solar energy harvesting. Two basic types of solar collectors may be distinguished as non-concentrating and concentrating collectors. Among concentrated collectors, compound parabolic collectors (CPC) and parabolic trough collectors (PTC) draw the attention of researchers and industrial developers. The greatest assets of the concentrated collectors constitute its unique features of capturing the solar rays, including diffuse rays, no tracking mechanism at low to moderate concentrations, minimal heat loss, and higher collector efficiency. A vast review of concentrated collectors in terms of their development, experimental and theoretical study, and recent applications was presented by Pranesh et al. [
1].
PTC belongs to the group of solar collectors, of which development began with their first practical application in 1898 [
2]. So far, the most significant use of this type of technology is the Solar Electric Generating System (SEGS) power plant (345 MWe) [
3] and Plataforma Solar de Almeria in Spain, equipped with PTCs of total capacity equal to 1.2 MW. PTC collectors can effectively generate heat between 50 °C and 400 °C [
1]. Additionally, there is reliable information showing that PTC technology is adequate for generating electricity using medium temperature (~400 °C) of HTF (heat transfer fluid) [
4].
The work presented in this publication refers to PTC for low-temperature residential applications. PTC has numerous applications, and research work was conducted in this area. Thus, Sayed and Sivaraman presented research on a hot water generation system using a PTC collector [
5]. The twisted tape inserts were used to enhance the thermal performance of the collector. Jamadi et al. [
6] compared flat-plate solar collectors (FPC) and PTC systems to heat a building space. Authors showed experimentally that the PTC system with lower occupied space could produce higher thermal energy quality compared to the FPC heating system. Zou et al. performed experimental work concerning small-sized PTC for water heating in cold areas [
7]. The authors concluded that the investigated PTC could collect solar radiation efficiently, demonstrating an event with an efficiency of 67% under solar radiation of less than 310 W/m
2. A solid review of prototypes of PTCs used for low-temperature demands was presented in Reference [
8].
PTC is commonly used in practice. Thus, it is given constant and extensive attention in the literature. Many studies were devoted to PTC efficiency improvement [
9], both experimentally and numerically [
10]. The literature also offers a vast number of results for experimental and theoretical computational testing carried out in the range of optical [
11] and thermal parameters of the collectors [
12]. The effect of the research is that PTCs are now a common and advanced structure intended for solar energy utilization. Most of the studies described in the literature related to the modeling and simulation of solar collectors were based on models making use of computational fluid dynamics (CFD) techniques. This is due to the full availability of CFD software and the relatively simple geometry of the collectors (constant along their length). Two major groups of studies on CFD models of solar collectors can be distinguished. The first comprises analyses conducted with a view to minimize heat losses to the environment. The second aims to enhance the solar collector performance via changes in design or solar fluid modifications [
13], using nanofluids, for instance.
As an example of the studies from the first group, one may highlight Reference [
14], which presented an analysis of compound parabolic concentrator (CPC) solar collectors mounted on a tracker with a single rotation axis. The CPC collectors, tested by the authors, used a novel solution, implementing new absorber tubes and multi-parabolic receivers with a total aperture surface area of 15.4 m
2. One of the tested collectors was lined with ethylene tetrafluoroethylene (ETFE) foil. The absorber tube was covered with a foil to minimize heat losses to the environment. The aim of the research was to establish the impact of the tested foil on the collector heat loss and efficiency. The authors not only carried out an experimental investigation but also built a CFD model of the new-design solar collector. The presented model reflected the full three-dimensional (3D) geometry of one row of the collector. The authors used the surface radiation model. The relative error obtained between the calculated and measured values of the fluid temperature was minimal. Even though the authors referred to their computational model as a simplified one, it contained 1.8 million computational cells mapping the collector geometry in great detail. The developed model confirmed the justification for the use of the proposed foil. However, due to the specificity of CFD calculations, the model is more suitable for steady-state testing (as performed by the authors) rather than for analyzing the collector operation in transient states.
In Reference [
15], Kaloudis et al. presented the results of their research on a PTC. The work was computational in nature. Nanofluids were used as the solar fluid to investigate the impact of the content of nanoparticles on the PTC efficiency. An extended model of the collector was proposed, comprising an outer glass tube, the space in between the tubes, the tube with an absorber, and the solar fluid. The model took account of the conjugate heat transfer, and radiation was modeled using the surface-to-surface model. Heat losses to the surroundings were taken into account by setting the convective boundary condition, where the heat transfer coefficient was a function of the wind speed (a similar approach was used in many other studies, e.g., References [
14,
16]). The model proposed in Reference [
15] was validated through a comparison with the experimental testing results available from Reference [
17] (the results are widely quoted in issues related to research on PTCs). The authors’ numerical model showed that a solar fluid containing up to 4% of aluminum (III) oxide nanoparticles improved the collector efficiency by about 10%.
A computational analysis using nanofluids in a direct absorption solar collector (DASC) intended for household/non-industrial applications was conducted in References [
18,
19]. In Reference [
19], the authors took account of the phenomenon of an increase in the absorptivity of the heat transfer fluid (HTF) depending on the amount of suspension containing carbon nanohorn particles. The variation in collector efficiency was obtained as a function of the nanoparticle concentration for two design variants of the collector (double-tube and triple-tube options). Like in other studies, the outer diameter of the solar tube or of the system built of a mirror and a solar tube was small compared to the total length of the collector under analysis. Due to the use of nanoparticles and the conjugate thermal and flow analysis solution to account for the energy source terms with due accuracy, testing of flows through the annular cross-section required a higher density of the numerical mesh in the solar fluid domain [
18,
19]. The considerable length of the analyzed collectors and the need to increase the numerical mesh density in the fluid domain resulted in a high number of computational elements, which makes the model less useful when it comes to transient-state analyses.
Evacuated tube solar collectors (ETSC) in a system with a water tank were analyzed in Reference [
20]. The authors performed a CFD optimization study of various system parameters and their impact on the system’s thermal performance. Among the other parameters, nanoparticle ratios in water were studied. Results proved that nanoparticle copper and an increase in their volume ratio gave the highest thermal performance among the nanoparticles used in the study. Farhana et al. [
21] presented a ground review of works showing the influence of nanofluid use on the thermal yield of direct absorbing solar collectors, where the latest results of experimental and computational works with CFD tools were presented.
Missirlis et al. [
22] analyzed the influence of different configurations of solar collectors on heat transfer through the DASC polymeric solar collector. The author took into account a few cases, allowing redesign of the collector without increasing production costs. The application of CFD methods enabled choosing one of the possible solutions for improving the heat transfer properties of the polymeric DASC.
Reference [
23] presented the results of CFD modeling of a PTC, in which annular receivers were partially filled with air. According to the authors, the solution is dedicated to PTC low-temperature applications. Non-vacuum PTCs are much cheaper than vacuum ones, but they are also affected by greater heat losses to the surroundings. Using numerical tools, Chandra et al. analyzed intermediate solutions, where the upper part of the annular space was filled with insulation, and the lower part on the side of the incident reflected concentrated radiation was filled with air. They validated their model against the testing data presented in Reference [
17]. According to the authors of Reference [
23], the most important conclusion drawn from their research was the reduction in convective losses by about 12.5% compared to the standard structure filled with air.
Apart from 3D models that take into account geometrical details and complex flows, one-dimensional (1D) models are also popular in the literature. These models were created a few decades ago and were developed ever since. Some of them evolved into commercial software intended for collector operation analysis, also in transient-state conditions. Among others, Hongbo et al. [
24] performed calculations using four 1D models available from the literature based on the energy and mass conservation equations. The models concerned tube-in-tube PTCs. The calculations results were compared to each other and experimental data. Nevertheless, the study was related to the steady-state problem only. The authors concluded that the tested 1D model is an excellent trade-off between the simplicity of the models and their accuracy. The difference between the calculated and measured values of the fluid temperature obtained in Reference [
24] was below 1 °C.
Yılmaz and Soylemez [
25] presented a mathematical model of a single-pass collector mounted on a tracker with a single-rotation axis. The model contained a module for thermal and optical calculations. After the model validation (against experimental data available from the literature), the model was used to analyze the collector performance characteristics under different operating conditions. The model enabled analyzing the impact of many factors (e.g., the wind speed and optical parameters) on the collector efficiency. However, the model is intended for steady-state analyses only.
Liang et al. [
26] presented a 1D mathematical model of a filled-type evacuated tube solar collector with a U-tube. The filled-type evacuated tube consisted of a glass evacuated tube, a U-tube, and a heat transfer component. To describe physical phenomena, the model proposed made use of thermal resistance equations instead of the energy equation. The model concerned steady states and was used to find the optimal value of thermal conductivity for the filling material, which ensured the highest efficiency of the collector.
A model of a classical collector with a concentrating (converging) mirror with a tube absorber in a glass envelope, enabling a transient-state analysis, was presented by Lamrani et al. [
27]. The collector was fitted with a sun-tracking system. The model input data wwre meteorological graphs. Like in most other studies, the model was validated against experimental data [
17]. Using the results, the authors analyzed the collector performance for a location in Morocco. They obtained conclusions concerning, among others, the maximum collector efficiency and the minimum and maximum amount of the useful heat gain per day.
Another 1D model intended for PTC transient-state performance analysis was presented by Xu et al. [
28]. Like in the papers referred to above, the PTC described had an evacuated-type receiver that consisted of an absorber, a glass envelope, and a vacuum annulus. It operated in an industrial plant. The model was validated using data from a 1-MW solar plant. The authors of the study proposed a fast and straightforward model to simulate the transient-state performance of PTCs. It was observed that the model could be helpful in identifying a safe operation mode of the PTC to achieve higher useful power outputs.
The presented literature review indicates that no analysis was made so far of the operation of a sun-tracked PTC in which each repetitive segment has two U-tubes containing the solar fluid. This paper proposes a 1D distributed parameter PTC model, including the solar fluid two-pass flow through a double U-tube. The model takes account of the collector operation in unsteady-state conditions. One arm of the U-tube is in the direct radiation area and the other is in the area of reflected radiation (concentrated on the absorber). The U-tubes are installed collision-free inside a solar glass tube, the inner surface of which is in direct contact with the absorber. The collector tube has the shape of a Dewar vessel, where space in between the tubes is filled with vacuum, to reduce heat losses from the absorber. The collector design is presented together with a mathematical model and a discretized CFD model constructed according to best practice resulting from the authors’ experience and the literature review. The CFD model results are used for preliminary computational verification of the developed in-house mathematical model. The verification consists of a comparison of the values of the solar fluid temperature at the U-tube outlet obtained for different boundary conditions. The in-house model is less demanding computationally, and less time is needed to perform the calculations compared to the CFD model. The model is suitable for both steady-state analyses and simulations of the collector unsteady-state performance (real solar collectors always operate in transient conditions). The developed mathematical model is similar to the model of the flat-plate solar collector (FPC) presented in Reference [
29].
The next stage of the mathematical model validation, presented herein, is its extensive experimental validation, which is planned to be performed further. A test stand is now being built for this purpose. It comprises a three-segment PTC equipped with a sun-tracking system with two axes of rotation and a data acquisition system. For the purpose of this work, the installation was put into operation for one day, and preliminary measurements, as well as comparisons, were performed.
3. Model Development
This section presents a proposal for a mathematical model of the PTC described in
Section 2. It is a one-dimensional (1D) distributed parameter model that is used to simulate the solar collector unsteady-state operation under any variable boundary conditions. Five nodes are taken into consideration in each cross-section of every segment, both in the direct and in the concentrated radiation area (
Figure 4). The nodes cover the inner and the outer solar tube, two layers of the absorber (in the form of thin aluminum sheets), and the solar fluid flowing through the U-tubes. The inner solar tube is divided into two control volumes. For the air inside the collector, the impact of direct and concentrated radiation is taken into account. It is assumed that the U-tubes are a part of absorber 2. Thus, in each analyzed cross-section of the collector, 11 nodes are obtained, for which energy balance equations are formulated. In this way, an 11 ×
N-node model is obtained, where
N denotes the number of nodes (the number of control volumes) in the direction of the fluid flow inside the U-tubes.
A basic simplification of the developed model is that the heat conduction in the circumferential direction by the aluminum lamellae, constituting the collector absorber, is omitted. The simplification is adopted to substantially reduce the number of differential equations solved at the same time, which makes it possible to obtain results much faster. The effect of this approach is uniform “heating” of absorbers to different temperatures in the area of direct and concentrated radiation. However, these temperatures will have a constant value on the respective circumference parts (falling on a given radiation area). Considering the high thermal conductivity of aluminum, in real conditions, the absorber in the direct radiation area will heat up quickly due to the heat conduction from the area of concentrated radiation. This will cause a temperature
ta1 and
ta2 increase in the circumference of absorbers (
Figure 4). In the proposed model, each of the absorber
N control volumes has a constant temperature on the part of the circumference exposed to the direct and concentrated radiation. As a result, the absorbers and the fluid temperature distributions along the U-tube length, obtained from the developed model, will be different compared to real conditions. The approach described in
Section 5 is proposed in the developed model in order to partially take into account the heat conduction in the circumferential direction.
Nonetheless, the total heat flux absorbed by the fluid will be the same in both cases. Therefore, the solar fluid temperature obtained from the proposed model at the collector outlet will be the same as in real conditions. The 3D model of the solar collector was also developed (
Section 4) to confirm this assumption. The 3D model takes entirely into account the heat conduction in the circumferential, radial, and axial direction (in all computational domains). Selected results are compared in
Section 5.
In the case of non-uniform heating on the circumference of thick-walled tubes made of a material with a much lower thermal conductivity coefficient (e.g., steel tubes), high thermal stresses may arise and affect operational safety. The stresses are the effect of significant temperature differences on the tube circumference and the tube wall thickness. It is then necessary for the mathematical model to take account of heat conduction in the circumferential and radial direction [
30,
31]. However, this lengthens the computation time considerably with no effect on the accuracy of determination of the fluid temperature at the tube outlet.
The other assumptions adopted in the collector mathematical model are as follows:
- -
Each segment of the collector operates in the same thermal and flow conditions;
- -
The fluid flow is uniform in all U-tubes of the collector composed of any number of segments;
- -
The heating of the walls of the U-tubes (being a part of absorber 2) is uniform on the circumference;
- -
The solar radiation reflected from the reflector (and concentrated on the absorber) is related to the outer surface of absorber 1 with diameter
dab (
Figure 4);
- -
The fluid and the absorber material physical properties are calculated on-line as a function of temperature;
- -
Physical properties of solar tubes are assumed as constant;
- -
On-line calculations are carried out of the coefficient of the heat transfer from the U-tube inner surface to the fluid;
- -
The collector loses heat to the environment via convection and radiation (conditions set on the outer surface of the outer glass tube);
- -
The radiative heat transfer occurs between the inner and the outer solar tube (due to vacuum, there is no convection);
- -
The air inside the collector is heated due to convection from absorber 2 (the air temperature is the effect of both direct and concentrated radiation);
- -
The radiative heat transfer occurs between the inner surfaces of absorber 2 (between the areas with temperatures ta4 and ta2, in the area of concentrated and direct radiation, respectively).
The assumptions formulated above do not simplify the model and do not change the physics phenomena occurring in the solar collector. They rather indicate the model complexity, e.g., taking into account radiative heat transfer between the inner surfaces of the absorber 2, as well as the on-line calculation of physical properties of the heat transfer fluid and absorbers. It should be noted that it is assumed that the temperature-independent physical properties of solar tubes do not affect the calculation results. The physical properties of the solar tubes change slightly, and the thermal conductivity remains constant and reaches 1.0 W/(m·K) within the working temperature range. Considering the small number of segments installed (
Figure 3), inequalities of the medium flow distribution through the separate U-tubes are negligible. Moreover, such a high thermal conductivity value of aluminum causes the U-tube to heat up quickly and equalize its temperature around the perimeter. Therefore, a one-dimensional heat transfer model is adopted in this case.
When considering the symmetry of the system, the energy balance equations for all elements of the collector (all the collector nodes) are written below in relation to one U-tube of one segment. The equations are derived for the area of direct radiation and the area of radiation reflected from the reflector (concentrated on the receiver).
The energy balance equations are written first for the control volume of one U-tube (
Figure 4) in the direct radiation area.
– Outer solar tube:
where
Cg1—thermal capacity of the outer glass tube, J/K;
tg1—temperature of the outer glass tube, °C;
τ—time, s;
Φam—heat flux lost to the environment, W;
Φr1—heat flux exchanged between solar tubes due to radiation, W;
Φs1—solar heat flux absorbed by the outer glass tube, W.
– Inner solar tube (control volume with node
tg2):
where
Cg2—thermal capacity of the analyzed control volume of the inner glass tube, J/K;
tg2—temperature of the inner glass tube, °C;
Φc1—heat flux conducted in the radial direction by the inner glass tube, W;
Φs2—solar heat flux absorbed by the analyzed control volume of the inner glass tube, W.
– Inner solar tube/absorber 1 (control volume with node
ta1):
where
Ca1—thermal capacity of the analyzed control volume, J/K;
ta1—temperature of absorber 1, °C;
Φc2—heat flux conducted in the radial direction by both layers of the absorber, W;
Φs3—solar heat flux absorbed by absorber 1, W.
– Absorber 2:
where
Ca2—thermal capacity of absorber 2, J/K;
ta2—temperature of absorber 2, °C;
Φconv—heat flux transferred by convection to the air inside the collector, W;
Φr2—heat flux exchanged by radiation between the surfaces of absorber 2 (in the area of concentrated and direct radiation), W;
Φ conv,f—heat flux transferred by convection to the fluid inside the U-tube, W.
– Solar fluid:
where
Cf—thermal capacity of the fluid, J/K;
tf—fluid temperature, °C;
Φf1—heat flux flowing into the control volume with the fluid, W;
Φf2—heat flux flowing out of the control volume with the fluid, W.
The energy balance equations below relate to the control volume of one U-tube (
Figure 4) in the area of radiation reflected from the reflector (concentrated on the absorber—subscript:
cd).
– Outer solar tube:
where
Cg3—thermal capacity of the outer glass tube, J/K;
tg3—temperature of the outer glass tube, °C;
Φam,cd—heat flux lost to the environment, W;
Φr1,cd—heat flux exchanged between solar tubes due to radiation, W;
Φs1,cd—solar heat flux absorbed by the outer glass tube, W.
– Inner solar tube (control volume with node
tg4):
where
Cg4—thermal capacity of the analyzed control volume of the inner glass tube, J/K;
tg4—temperature of the inner glass tube, °C;
Φc1,cd—heat flux conducted in the radial direction by the inner glass tube, W;
Φs2,cd—solar heat flux absorbed by the analyzed control volume of the inner glass tube, W.
– Inner solar tube/absorber 1 (control volume with node
ta3):
where
Ca3—thermal capacity of the analyzed control volume, J/K;
ta3—temperature of absorber 1, °C;
Φc2,cd—heat flux conducted in the radial direction by both layers of the absorber, W;
Φs3,cd—solar heat flux absorbed by absorber 1, W.
– Absorber 2:
where
Ca4—thermal capacity of absorber 2, J/K;
ta4—temperature of absorber 2, °C;
Φconv,cd—heat flux transferred by convection to the air inside the collector, W;
Φconv,f,cd—heat flux transferred by convection to the fluid inside the U-tube, W.
– Solar fluid:
where
Cf,cd—thermal capacity of the fluid, J/K;
tf,cd—fluid temperature, °C;
Φf1,cd—heat flux flowing into the control volume with the fluid, W;
Φf2,cd—heat flux flowing out of the control volume with the fluid, W.
The last energy balance equation relates to the air inside the collector. The air temperature is the effect of heat transferred from absorber 2 in the direct and concentrated radiation area. The energy balance equation has the following form:
where
Cair—thermal capacity of air, J/K;
tair—air temperature, °C.
In each analyzed cross-section of the collector, 11 differential equations are obtained. The equations are solved by approximating derivatives using appropriate difference schemes. The relations derived in the process make it possible to calculate temperature histories in 11 nodes of each control volume of the collector, and they are characterized by simple notation. Therefore, results can be obtained very quickly. The in-house computational program, written based on the derived relations, additionally makes it possible to conduct a detailed analysis of the impact of each term of the equation on the results. The method is iterative. All thermophysical properties are calculated on-line as a function of temperature using the relations given in
Section 4. The heat transfer coefficient values are also calculated on-line.
In order to obtain a stable numerical solution, the time and spatial steps are selected to meet the Courant–Friedrichs–Lewy condition for one-dimensional problems as follows [
32]:
where Δ
τ—time step for the calculations, s; Δ
z—spatial step (length of the collector control volume), m;
w—fluid flow velocity in U-tubes, m/s.
The values of concentrated solar radiation are calculated using the following balance equation (
Figure 4):
which gives
where
dgl—diameter of the outer solar tube, m;
dab—outer diameter of the absorber, m;
Gd—intensity of direct solar radiation, W/m
2;
Gcd—intensity of concentrated solar radiation, W/m
2;
L—absorber length, m;
W—reflector width, m (
W = 0.418 m—
Figure 2).
In the model developed, the influence of mutual shading of individual tubes was omitted, since the solar collector is assembled with a two-axis sun-tracking system. It becomes a significant improvement over industrial applications where one-axis sun-tracking is commonly applied. In the two-axis tracking, the collector plane is always set perpendicularly to the incident solar radiation, which avoids shading. On the other hand, in a one-axis tracking system, shading may still occur.
The derived relations are valid for the analyzed vacuum solar PTC built of any number of repetitive segments. The proposed model is also suitable (after necessary corrections to the equations) for a collector with a single solar tube (no vacuum) or with a single U-tube inside the solar tube. After some structure-related simplifications, the model is also valid for the most popular concentrating solar collectors with the absorber in the form of a single straight tube.
Firstly, in order to verify the proposed 1D model, an additional 3D model of the analyzed collector was built in the ANSYS Fluent environment (
Section 4). The 3D model includes the following domains: absorbers, U-tube, heat transfer fluid, and air inside the collector. Therefore, in the 3D model, solar tubes are omitted, since the model was developed to compare, with the 1D model results, the transient temperature courses of the heat transfer. Accuracy of the 1D model depends significantly on the adopted time and spatial steps, Δ
τ and Δ
z, respectively. Thus, a comparison of the 1D and 3D model results is aimed at finding such time–space step values, for which a satisfactory correspondence of the calculation results is obtained.
In addition to the computational verification, a pre-experimental verification using the Δτ and Δz steps identified was also performed. The test stand, currently still under construction, was started for one day, to enable the necessary measurements to be made.
Extensive experimental verification of the 1D model is planned to be carried out in the future, using the Δτ and Δz values verified in this article.
4. 3D CFD Model
This section presents the simulation model and methodology used to simulate the operation of the U-tube parabolic trough collector (PTC). The collector model was built using the ANSYS Fluent 19.3 code [
33], where thermal and fluid dynamic analyses can be performed simultaneously. This software is the leading tool for flow and heat simulations. It was also used in previous studies to simulate various solar thermal collectors [
12,
18,
23,
34]. The developed model is three-dimensional. It is assumed that the model is symmetrical, and the symmetry plane is the
YZ surface passing through the center of the solar pipe. Water flows through the upper tube of the absorber, then turns on the U-tube and flows in the second pass through the lower tube, as shown in
Figure 1 and
Figure 4. A modification is introduced here, whereby the tube is shifted to keep the symmetry plane. In the real object, the first and the second pass are in different halves of the solar pipe, as described in
Section 2 (see
Figure 1). In the collector model geometry, the second and the first pass are in the same half. This approach is correct because the upper pipe is in the direct radiation area, and the lower pipe is in the concentrated radiation area. This simplification makes it possible to achieve a more economical model without compromising the model reliability.
The model omits the glass tubes and the vacuum space between them. The thermal load on the absorber surface in the direct and in the concentrated radiation area is reduced by losses to the environment and losses due to thermal resistance of the glass tubes and vacuum. The direct and concentrated radiation areas are marked in
Figure 5 as
Wall_d and
Wall_c, respectively. The geometric model is 1600.0 mm long, and the diameter of the outer domain corresponding to the outer diameter of the tube formed by the sheet metal of absorber 1 is 47 mm (
Figure 1). The other dimensions of the U-tube and metal sheets are as shown in
Section 2 in
Figure 1. The model takes into account that the space between the U-tube and the metal sheets is filled with thermal paste.
The governing equations of the computational fluid dynamics (CFD) analysis include the mass conservation equation, the momentum conservation equation, and the energy conservation equation. The CFD model calculates the conjugate heat transfer between the fluid and the solid material, the convective heat transfer in the fluid regions, and heat radiation. Radiation
qrad in the model occurs on the inner surface of the inner metal sheet due to temperature non-uniformity (
Figure 5). The non-uniformity is the effect of the heat flux values affecting the upper and the lower half of the outer metal sheet and of the HTF interaction. The radiative heat transfer inside the collector was calculated using the surface-to-surface (S2S) radiation model.
The HTF physical properties are assumed to be temperature-dependent. The following equations for density, specific heat, thermal conductivity, and viscosity calculations with respect to temperature were introduced into the CFD model:
where
T is the fluid temperature in K.
The above equations are valid in the temperature range from 0 °C to 100 °C and are the same as those used in the in-house program. Similarly, the properties of aluminum (specific heat and thermal conductivity) were introduced into the model according to the following equations [
35]:
where
t is the solid temperature in °C.
The velocity of the HTF flowing through the collector tubes is low, from about a few mm/s to about a dozen cm/s [
36]. The Reynolds number values range only from a dozen to several hundred when considering such a small diameter of the tubes. Thus, the Reynolds number value is much lower than the critical value of
Recr = 2300 [
29,
36]. Due to the above, it is assumed that the HTF flow is laminar.
The boundary conditions assigned to the CFD model are as follows:
- ▪
Inlet—a uniform velocity inlet boundary condition is adopted at the U-tube inlet where vz = v, vy = vy = 0; T = Tin;
- ▪
Outlet—an outflow boundary condition is used at the U-tube outlet;
- ▪
Wall—no-slip conditions on the inner absorber surfaces;
- ▪
The top outer wall (Wall_d) of the absorber is subjected to heat flux Gd, which corresponds to the heat load from the direct radiation area;
- ▪
The bottom outer wall (Wall_c) of the absorber is subjected to heat flux Gcd, which corresponds to the heat load from the concentrated radiation area.
The calculations start with the initial condition of zero velocity and uniform temperature throughout the calculation domains. The calculations are performed using the coupled algorithm and PRESTO! discretization for pressure, while the second-order upwind approach is followed for other quantities. The numerical solutions are convergent when the residuals of continuity and of the momentum equations are less than 10−4, and the residuals of the energy equation are less than 10−7. Furthermore, the convergence criterion must ensure that the mass-weighted average temperature at the outlet is constant.
A hexahedral-dominant mesh is created to discretize the computational domain to obtain results with high levels of accuracy. There is a substantial difference between the dimension of the tube length (1600 mm) and the absorber tube outer diameter (47 mm). Moreover, the HTF flow in the U-tube has a significant impact on collector thermal efficiency. Due to that, a fine mesh is required in the cross-section of the computational domain with a refined mesh along the U-tube perimeter. The refined boundary-layer mesh is also applied to the inner surface of absorber 2 on the air-washed sides (to capture the air movement caused by natural convection as a result of air heating and the difference in density).
Figure 6 illustrates the front cross-section divided into control volumes. The model has 1.1 × 10
6 elements. It is built by sharing the frontal surface of the model (
Figure 6) in the first place. The grid was compacted in the areas of expected large gradients in the wall-fluid contact areas. Mesh tools are adopted to create the hexa dominant mesh. Then, the mesh is created on the frontal surface, and the periodic mesh operation is used to generate the entire mesh of the PTC collector model. This way, the dominant hexa mesh is obtained, and the division of the model into subdomains connected by interfaces is avoided.
The CFD model is a model where coupled heat transfer is considered, and the results are affected by the division of the domain into control volumes. Therefore, the grid independence study was performed to make sure that the numerical solution is independent of the mesh resolution. The parameters for the grid-independent verification test include the overall heat transfer rate from the absorber tube wall to heat transfer fluid, and the HTF outlet temperature. The grid independence study was carried out for four grids with the following element numbers: 0.5 × 106, 0.71 × 106, 1.1 × 106, and 1.7 × 106. As a result, the outlet temperature of the HTF remained constant when the number of elements was greater than 1.0 × 106. Therefore, the model with the number of elements equal to 1.1 × 106 was selected for CFD analysis based on the comprehensive consideration of computational time and accuracy of the results obtained.
5. Results and Discussion
Firstly, a number of calculations were performed for different boundary conditions to verify the accuracy of the results obtained from the in-house mathematical model of the collector. The same conditions were also used to carry out CFD calculations. Selected results of the calculations and comparisons are presented in
Table 2 and in
Figure 7,
Figure 8,
Figure 9 and
Figure 10.
The 3D CFD model (
Section 4) takes no account of the two solar tubes, which is due to the fact that the model was developed to verify the in-house model comprising only the most essential elements of the solar collector, i.e., absorbers, U-tubes, and the solar fluid. It was not the task of the 3D model to precisely map the entirety of the occurring phenomena, which was the case in other works, e.g., References [
23,
34]. For this reason, the heat flux set on the absorber outer walls (
Wall_d and
Wall_c) was the same as in the in-house model, as were the mass flux and the HTF temperature. The calculations aim to verify the correctness of the fluid temperature at the collector outlet determined using the in-house model (the 1D model should satisfy the mass and energy balance). Considering the above, the aspects of solar radiation transmissivity through the glass division walls are omitted in the developed 3D model.
The assumed volume fluxes of the fluid shown in
Table 2 relate to the entire collector, built of three segments (
Figure 3). Considering that, in every segment, there are two U-tubes, it is assumed that the fluid flows uniformly through six U-tubes (
Table 2 shows the assumed flow velocity in one U-tube). The calculations were performed for three selected heat flux values of solar radiation. i.e., 600, 800, and 1000 W/m
2. In each analyzed case, it is assumed that the initial temperature distributions in the two absorbers, in the fluid and in the air inside the collector, are the same and total 40 °C. Moreover, in the developed in-house model, it is assumed that the initial temperature distribution of the outer solar tube is equal to the ambient temperature of 20 °C. The initial temperature distribution of the inner solar tube is assumed as uniform and equal to 40 °C. The assumed wind speed totals 3 m/s.
Satisfactory agreement between the developed 1D model and the 3D CFD model results was obtained for calculations performed using the in-house model with the time step of 0.001 s and the spatial step of 0.04 m. In total, 80 control volumes were obtained in the direct and concentrated radiation area. Considering that there are 11 nodes in every control volume (
Figure 4), 880 differential equations were solved in every time instant. The computational program is written in the Fortran code [
37].
When analyzing
Table 2, it can be seen that the fluid temperature at the collector outlet obtained from the in-house model differs from the CFD model value by 0.1 °C maximum, which means that the error totals about 0.21%. The obtained agreement between the results should be considered as fully satisfactory. The time of the in-house model calculations (until the steady state was achieved) totaled about one minute. Such a short computation time, at the fully satisfactory accuracy of the result, is a fundamental advantage of the developed in-house model.
The equivalent value of the heat flux absorbed partially by the absorbers in the direct radiation area (
Figure 4) is determined to take into account partial heat conduction of both absorbers in the circumferential direction. In real conditions, these parts of the absorbers will heat up quickly due to heat conduction in the circumferential direction. The heating is simulated using the equivalent heat flux value found from the balance Equation (21). In the equation, it is taken into account that, in real conditions in the direct radiation area, solar radiation is absorbed by the surface
, whereas, in the concentrated radiation area, it is absorbed by the surface
. Assuming that the heat flux with equivalent intensity
Gequiv is absorbed in the concentrated radiation area, and that half of this value is absorbed in the direct radiation area, the following balance equation can be written:
After relevant transformations, the following relation is obtained (the relation was used to calculate the equivalent value of the solar radiation intensity):
Selected results of calculations and comparisons for the assumed solar radiation intensity
G = 600 W/m
2 are shown in
Figure 7 and
Figure 8, and those for
G = 1000 W/m
2 are shown in
Figure 9 and
Figure 10. In both cases, the volume flux of the fluid flowing through the entire solar collector (built of three segments) totals 1 L/min. The computational mass flux of the fluid flowing through one U-tube is 0.00285 kg/s.
Figure 7 and
Figure 9 present a comparison of the fluid temperature distribution along its flow path through the U-tube for the final steady state. The U-tube total length in contact with the absorber is 3.2 m (1.6 m in the direct and in the concentrated radiation area, respectively). As expected, the distributions obtained from the developed in-house model differ from those produced by the 3D CFD model.
Since the heat conduction is omitted in the 1D model, the fluid in the direct radiation area gains temperature more slowly than when considering the heat conduction (3D model). Consequently, the high temperature of the absorbers in the concentrated radiation area causes a rapid increase in the fluid temperature for the 1D model.
The situation is entirely different for the 3D model, where heat conduction in the circumferential, axial, and longitudinal directions is fully considered. In this case, the working fluid heats up faster in the direct radiation area (compared to model 1D) because the absorbers reach a higher temperature in this section (due to heat conduction from the concentrated radiation area). At the same time, the heat drain in the absorbers from the concentrated radiation area causes the fluid to heat slowly in this area. Thus, the temperature distributions of the fluid, obtained by means of the 3D model (
Figure 7 and
Figure 9), are characterized by a smooth course. The distributions obtained using the 1D model indicate specific refraction of these distributions over the U-tube length of 1.6 m, i.e., where the direct radiation ends and the concentrated radiation starts to act. The analysis of
Figure 7 and
Figure 9 shows that, in spite of different distributions along the U-tube length, the temperature of the fluid at the outlet, obtained by means of both models, is the same (at 3.2 m—outlet cross-section).
To confirm the adopted assumption that, regardless of different temperature distributions on the circumference of absorbers, the total heat flux absorbed by the fluid will be the same in both models (especially when absorbers are made of a high thermal conductivity material),
Figure 8 and
Figure 10 present comparisons of the fluid temperature histories at the U-tube outlet. A slightly bigger difference can be noticed in
Figure 8 for the time interval of about 200 s to about 300 s. The maximum difference is about 0.25 °C, which is negligibly small. It may, therefore, be stated that the agreement observed in
Figure 8 and
Figure 10 between the histories of the fluid temperature at the collector outlet obtained from the in-house 1D model and the 3D CFD model is fully satisfactory. This confirms the correctness of the adopted assumption that the total heat flux absorbed by the fluid does not depend on the temperature distribution on the absorber circumference (omitting heat conduction in the circumferential direction of absorbers in the in-house model does not generate errors in the results of the calculations of the fluid temperature at the U-tube outlet). It also confirms the effectiveness of the developed mathematical model and its usefulness for simulations of the operation of vacuum solar parabolic trough collectors. In calculations of thermal systems equipped with solar collectors, the user’s/designer’s attention is focused on the achieved value of the fluid temperature at the collector outlet and not on the temperature distribution on the collector length.
As mentioned earlier, in both cases (1D and 3D models), the total heat flux absorbed by the fluid along the entire length of the U-tube is the same. Therefore, the temperature of the fluid at the collector outlet (only at the outlet) is the same. To further confirm this phenomenon, the computed values of the heat flux absorbed by the fluid along the length of the U-tube in both models are given below in
Table 3 and
Table 4. These values refer to steady state. Since, in the case of the 1D model, the assumed spatial step of the calculation is
z = 0.04 m, 40 cross-sections are obtained in both radiation areas. Computed values of heat fluxes were summed up and compared with the values obtained from the 3D model (
Table 4). The tables refer to the data based on which
Figure 7 and
Figure 8 were developed (
G = 600 W/m
2,
).
For the 1D model, in each analyzed section, the heat flux absorbed by the fluid was calculated as follows:
where
h—heat transfer coefficient from the inner surface of the U-tube to the fluid, W/(m
2K);
din—internal diameter of the U-tube, m;
ta—absorber temperature, °C;
tf—fluid temperature, °C.
The temperature of the fluid at the outlet of the U-tube is, therefore, equal to
In order to illustrate the differences in the heat fluxes absorbed by the fluid in the individual sections of both models, the 3D CFD U-tube model was also divided into 40 control volumes in both radiation areas. The calculation results are presented in
Table 4.
For the 3D CFD model, the fluid temperature at the U-tube outlet is equal to
In the equations above, it is assumed that the specific heat of 50% propylene glycol is 3600 J/(kg·K). The obtained temperature values are in accordance with the values given in
Table 2 and shown in
Figure 7 and
Figure 8 (for the steady-state operation in the outlet section of the collector tube).
An analysis of
Table 3 and
Table 4 indicates that different values of the heat flux absorbed by the working fluid along the solar tube length, respectively, in 1D and 3D models, occurred. However, the total (summary) values of the heat flux absorbed along the whole U-tube length are almost the same in both cases. The temperatures of the fluid at the outlet of the U-tube are, therefore, the same. In the case of the 1D model (
Table 3), in the individual cross-sections of the direct radiation area, the heat flux absorbed by the fluid changes slightly, similarly, in the area of concentrated radiation, where the heat flux is larger. This is due to the fact that the heat conduction in the absorbers from the concentrated radiation area to the direct radiation area is not included in the 1D model. The inclusion of heat conduction in the 3D model causes a visible decrease in the heat flux absorbed by fluid along both arms of the U-tube (
Table 4).
In order to confirm the conclusions of the computational verification, a preliminary experimental verification was also carried out. The research stand is not yet fully completed, but it was possible to run it for one day and carry out the measurements necessary for the purpose of the article. Selected results of the measurements and comparisons are shown in
Figure 11,
Figure 12,
Figure 13 and
Figure 14.
Figure 11 shows the intensity of solar radiation (heat flux of solar radiation) measured with a pyranometer. The total volume flow of the working fluid flowing through the collector (
Figure 12) and the ambient temperature were also measured.
Figure 13 shows the measured inlet and outlet temperatures of the working fluid along with the outlet temperature, computed with the proposed 1D model. The measured solar radiation intensity, fluid volume flow, and its temperature at the collector inlet, as well as ambient temperature, were the input data for the computations. The analysis of
Figure 13 shows that the measured and computed outlet collector temperature is fully satisfactory. To illustrate this comparison in more detail, a section of
Figure 13, covering a shorter period, is shown in
Figure 14.
These comparisons, being the preliminary experimental verification of the developed 1D model, fully confirm the conclusions of the computation verification. The proposed model reflects with high accuracy the measured course of the fluid temperature at the collector outlet. It should be emphasized that the measurements were made within a single day when cloudiness occurred, which caused high variability of the measured solar radiation intensity (
Figure 11) and the working fluid temperature (especially at the outlet—
Figure 13 and
Figure 14). Therefore, it was a demanding verification of the mathematical model, which confirmed its usefulness for simulating the operation of the analyzed collector in transient conditions.