Journal of Renewable Energy

Volume 2013 (2013), Article ID 389514, 8 pages

http://dx.doi.org/10.1155/2013/389514

Research Article

Modeling and Characteristic Analysis of a Solar Parabolic Trough System: Thermal Oil as the Heat Transfer Fluid

Zhai Rongrong, Yang Yongping, Yan Qin, and Zhu Yong

School of Energy Power and Mechanical Engineering, North China Electric Power University, Beijing 102206, China

Received 31 January 2013; Revised 27 May 2013; Accepted 7 July 2013

Academic Editor: Shuhui Li

Copyright © 2013 Zhai Rongrong et al. 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.

The thermal oil is applied as the heat transfer fluid in a solar parabolic trough collector system. Firstly, the system dynamic model was established and validated by the real operating data in typical summer and spring days in references. Secondly, the alteration characteristics of different solar radiation, inlet water temperature and flow rate, and collectors’ area and length are analyzed and compared with the normal working condition. The model can be used for studying, system designing, and better understanding of the performance of parabolic trough systems.

1. Introduction

Parabolic trough solar collector is a kind of line-focus concentrating solar utilizations for industry, with medium temperature, maximum working temperature 400–500°C. According to the heat transfer fluid types in the absorber, thermal oil and water/steam are the basic sorts. Both the two systems are wildly used in generating or heating in modern industry in recent years. The solar energy generating system with thermal oil as the heat transfer fluid, called SEGS, has been developed in several countries, especially in the USA. Therefore, the generated green electricity has increased the impetus of future research and development.

Some of the published papers have studied the modeling and analysis of the heat transfer process of the thermal oil in parabolic trough collector systems. The one-dimension parabolic trough collector model had been developed by Odeh et al. [5], including synthetic oil and water/steam as the working fluid. The created model can be used to predict the performance of the collector. Rolim et al. [6] introduced an analytic model for a solar thermal electric generating system with parabolic trough collectors. The energy conversion of solar radiation into thermal power along the absorber tube of the parabolic collector is studied, taking into account the nonlinearity of heat losses and its dependence on the local temperature. The simulated results have been compared with experimental data. García-Valladares and Velázquez [7] developed detailed numerical simulations of thermal and fluid-dynamic behaviors of a single-pass and double-pass solar parabolic trough collector. The single-pass solar device numerical model has been carefully validated by experimental data. The parabolic trough collector and the main generating systems of the SEGS VI solar thermal oil power plant in California have been modeled by Stuetzle et al. [1, 2]; the control algorithms were established in MATLAB to control outlet temperature and flow rate. The one-dimension and two-dimension thermal oil solar parabolic collectors’ models were developed using EES software and validated by experimental data [8]. The detailed numerical simulation model of a single-pass and double-pass solar parabolic trough collector has been developed and validated by experimental data obtained from Sandia National Laboratories [7].

However, the heat transfer has been analyzed in the above papers, but it has not been studied in detail in variable working conditions in these papers. The novelty of this paper lies in that it is focused on the details of heat transfer process with times along the flow direction, especially in the variable working conditions.

In this paper, the Euler transformation method is introduced based on the forward modeling approach to simplify the dynamic model of parabolic trough with thermal oil as the heat transfer fluid. Meanwhile, the alteration characteristics in different input parameters are analyzed. The alterations include solar radiation, inlet oil temperature, and inlet oil mass flow rate and collector area. The model can be used for the analysis of parabolic trough system performances in variable working conditions.

2. Mathematic Modeling

As shown in Figure 1, the linear parabolic concentrators are used to focus sunlight to heat the pumped heat transfer fluid in the absorber. The good light transmittance and heat-durability glass envelope are used to maintain the vacuum space between the absorber and the envelope in order to reduce heat loss and prevent corroding [9]. On the one hand, there should be vacuum between the receiver and glass envelope to reduce the heat loss. On the other hand, the luminousness should be increased by enhancing the optical performance of glass envelope. It means that the performance of luminousness and anticorrosion of glass envelope should be very good; otherwise there will be no vacuum when the glass is corroded. The heat and mass balance of absorber is listed in Figure 2.

Figure 1: The schematic diagram of a parabolic trough solar collector system.

Figure 2: Schematic diagram of the heat transfer analysis [1–3, 8].

2.1. Energy and Mass Balance Equations

According to [1–3, 8], the following energy and mass balance equations are used to establish the model. Equations (1) and (2) are for thermal oil; (3)–(5) and (6)–(8) about absorber and glass envelop are the same as the heat balance in direct steam generating system.

Consider

Combining (2) and (3),

Combining (5) and (6),

The above mentioned inlet or outlet heat and for absorber and glass envelope all consisted of two parts: radiation and convection. The solar radiation absorbed by per unit area absorber [1, 3, 8] can be expressed as follows:

The heat absorbed by thermal oil is ; .

2.2. Equations Solving

The explicit Euler method is introduced to solve the above mentioned partial differential equations to avoid the complex solving, if the partial differential solution is used to all the individual collectors.

2.2.1. Introduction of the Explicit Euler Method

For normal ordinary differential equation (10),

where as the solution interval is partitioned equally, and as the step length. Then, the solution is marked as , , separately. The solution of at solved by forward difference method is as follows:

At the same time, ; then,

The approximation of is expressed as (13) or (14),

or

Combining (15) and (16), the approximate forward Euler solving solution of is as follows:

As shown in (17), the explicit Euler method is defined, and the value of is calculated by directly.

2.2.2. The Euler Transformation for Thermal Oil

The energy balance equation for thermal oil based on (1) and (2) can be transformed into the following:

Each collector is defined as a whole, and the predigesting is as follows:

Combining (19) and (18), then,

Equation (20) is sorted as follows:

For a whole collector, the outlet temperature of thermal oil is transformed by explicit Euler method as follows:

2.2.3. The Euler Transformation for Absorber

The heat transfer expression for a whole collector absorber is as follows:

The transformed equation is as follows:

2.2.4. The Euler Transformation for the Glass Envelope

The heat transfer expression of glass envelope is as follows:

The transformed equation is as follows:

The model and the program are established according to the aforesaid equations and solving approach. The schematic chart of the computer program is shown in Figures 3 and 4.

Figure 3: The schematic flow chart of the program to calculate the absorbed heat by absorber.

Figure 4: The schematic flow chart of the program to simulate the thermal oil system.

3. Model Validation

The operating data from real solar power plant in SEGS-VI [1–3] is used to validate the developed model. The solar collectors are LS-2, manufactured by LUZ Company. Other parameters of the power plant are listed in Table 1. The synthesized oil, THERMINOL VP-1, supplied by SOLUTIA Company, is selected as the heat transfer fluid. The operating temperatures of the thermal oil are from 12°C to 400°C in liquid and steam operating conditions. The main physical and chemical characteristics are shown in Table 2 [4].

Table 1: The system main parameters of the SEGS VI power plant [1–3].

Name
Parameters
Value
Unit

Power block
Name
SEGS VI in Kramer Junction, USA

Net capacity
30
MWe

Designed oil temperature
349
°C

Solar collector area
188000
m²

Solar collectors
Length of single collector
47.1
m

Width of single collector
5
m

Area of single collector
235
m²

Focal length of the collector
1.84
m

Focal ratio of the collector
71
—

Absorber external diameter
70
mm

Envelope external diameter
115
mm

Table 2: Main physical and chemical characteristics of THERMINOL VP-1 [4].

Composition
Diphenyl oxide/diphenyl

Appearance
Clear/sediment free liquid

Maximum temperature
400°C

Kinematic viscosity (40°C)
2.48 mmm²/s

Density (15°C)
1068 kg/m³

Fire point
127°C

Autoignition temperature
621°C

Pour point
12°C

Normal boiling point (1013 mbar)
257°C

Coefficient of thermal expansion (200°C)
0.00097/°C

Average molecular weight
166

In SEGS VI power plant, 50 groups of paratactic single circle parabolic trough collectors are the main heat resource. Each single circle is composed of 16 collectors connected in series. The outlet heated thermal oil enters into power block to heat water and generate steam, and then the cooled oil is back to each single circle by oil pump.

It is obvious that the thermal oil temperature will change along the collectors. The physical characters will change in different temperatures, according to the data from SOLUTIA company [4]. The fitted equations between parameters and temperature are gained. The relation of density (kg/m³) and temperature is as follows:

The relation of heat capacity (kJ/kg·°C) and temperature is as follows:

The relation of heat conductivity (W/m·°C) and temperature is as follows:

The relation of Kinematic viscosity (mm²/s) and temperature is as follows:

Two typical operating data in spring (May 20th) and summer (June 11) days are selected as the base working conditions to validate the model. The simulated results will be compared with measured values in integral hours from 8:00 am to 7:00 pm. All the input parameters are set the same as the real SEGS VI power plant. The comparisons are drawn in Figures 5 and 6.

Figure 5: Comparison of simulated and measured results in spring working condition.

Figure 6: Comparison of simulated and measured results in summer working condition.

Figures 5 and 6 show the comparison in spring and summer days, respectively. The solar radiation rate, oil input temperature, and flow rate are all input parameters in the model; the oil output temperature is the only output parameter. It can be seen in the two figures that the precisions are acceptable for engineering, with maximum temperature differences about 5°C. But it is negligible in steady solar radiation days for both spring and summer working conditions. The main sources of differences are caused by measuring hysteresis. For the real system, the real output oil is heated along the collectors, but the measured input solar radiation is only one point. The differences will be more obvious in multivariate solar radiation working days.

4. Characteristic Analysis

The working condition at 12:00 am in a summer day is selected as the base condition to analyze the relations with different alterations. The alterations include solar radiation, oil inlet temperature, oil inlet flow rate, and collector length (area). The details of the base conditions are listed in Table 3.

Table 3: Parameters of the base condition.

Parameters
Unit
Value

Time
H:m
12:00

Solar radiation
W/m²
925

Inlet oil temperature
°C
285

Output oil temperature
°C
360

Output oil mass flow rate
kg/s
500

Collectors total area
m²
188000 (235 × 16 × 50)

Collectors total length
m
37600 (47 × 16 × 50)

4.1. Solar Radiation Alteration

Solar radiation alteration simulates the possible working conditions of radiation sharp changes in cloudy days or when the parabolic trough solar collectors are barred by surrounding buildings. Then based on the base working condition in Table 3, the output oil temperature will change following the changed solar radiation but maintaining other inlet parameters. Figure 7 shows the relation between output oil temperature and solar radiation alterations from the range of 805 W/m² to 1045 W/m². In Figure 7, every changed solar radiation 20 W/m² will increase or reduce the output oil temperature about 2°C. Therefore, the figure can be used to analyze or forecast the operating performances in the range of solar radiation changing working conditions.

Figure 7: The relation between output oil temperature and solar radiation alterations.

4.2. Inlet Oil Temperature Alteration

In the thermal oil heat transfer system, the heated oil will be cooled in the oil/water heat exchanger to heat water and generate steam. Then, the cooled oil will return to the entry of solar collector as a new beginning point for a circle. Therefore, the performance of cooling system will affect the output parameters of heat transfer oil directly. Figure 8 shows the changes of output oil temperature with the inlet temperature from 258°C to 294°C at the interval of 3°C. Therefore, the steady heat transfer is an important operation evaluation point.

Figure 8: The relation between output oil temperature and oil inlet temperature alterations.

4.3. Inlet Oil Flow Rate Alteration

Flow rate influences the heat transfer effect directly, and the output oil temperature will be reduced when increasing flow rate in the same solar radiation working condition. Therefore, adjusting the inlet flow rate is the main way to control and set the output parameters in operation. Figure 9 shows the changed output oil temperatures when changing the inlet mass flow rate in the range of 625 kg/s–375 kg/s. It is clear that the rate speed from 350 kg/s to 500 kg/s is much faster than that in the range of 525 kg/s–625 kg/s. The maximum temperature difference between 350 kg/s and 625 kg/s is about 80°C.

Figure 9: The relation between output oil temperature and inlet oil mass flow rate alterations.

4.4. Collectors Area (Length) Alteration

As mentioned earlier, the collectors in SEGS VI power plant consisted of 50 groups of single circles. And each single circle is composed of 16 LS-2 collectors connected in series. The alteration sources of collectors’ area and length are one or more collectors fault. Figure 10 shows the relation between output oil temperature and the area of collectors. The range of area and length in Figure 10 means the removed one to nine series-wound collectors (each 47 m length and 235 m²) in a single circle. The average changed temperature with every removed collector is about 8°C. In addition, the changing speed rate from 1645 m² to 2350 m² is sharper than the other range. For example, the output temperature difference between 1645 m² and 1880 m² is 5°C, but it is 9.4°C between 3525 m² and 3760 m².

Figure 10: The relation between output oil temperature and collector area (length) alterations.

5. Conclusions

In the paper, the thermal oil heat transfer process of the solar parabolic trough collector from input solar radiation to heated oil output is analyzed. The process is sorted into the three following separate sections: glass envelop heat transfer with solar radiation, absorber heat transfer with solar radiation, and absorber heat transfer with thermal oil. Then the partial differential equations are simplified by, Euler method to develop the dynamic model. The model has been validated by real operating data in typical spring and summer days. In the end, the alteration characteristics of different solar radiation, inlet water temperature and flow rate, and collectors’ area and length are analyzed and compared with the selected normal working condition. The results can be used to research and forecast the performance of solar collectors in special working conditions.

Nomenclature

: | The absorbed heat by per unit length absorber from the solar radiation, kJ
: | The absorbed heat by per unit area absorber from the solar radiation, kJ
: | The transferred heat by per unit length absorber with the glass envelope, kJ
: | The transferred heat by per unit length glass envelope with atmosphere, kJ
: | The transferred heat by per unit length absorber with thermal oil, kJ
: | The absorbed heat by the absorber, kJ
: | The absorbed heat by the glass envelope, kJ
: | The transferred radiant heat by per unit length absorber with the glass envelope, kJ
: | The transferred convectional heat by per unit length absorber with the glass envelope, kJ
: | The transferred radiant heat by per unit length glass envelope with atmosphere, kJ
: | The transferred convectional heat by per unit length glass envelope with atmosphere, kJ
: | The temperature of the absorber, °C
: | The temperature of the glass envelope, °C
: | The cross-sectional area of the inner absorber, , m²
: | The cross-sectional area of the absorber, , m²
: | The cross-sectional area of the glass envelope, , m²
: | The coefficient of heat transfer between absorber and water/steam, W/(m²·°C)
: | The density of thermal oil in current condition, kg/m³
: | The density of the absorber in current condition, kg/m³
: | The density of the glass envelope in current condition, kg/m³
: | The specific heat of thermal oil in current condition, kJ/(g·°C)
: | The specific heat of the absorber in current condition, kJ/(g·°C)
: | The specific heat of the glass envelope in current condition, J/(g·°C)
: | The exterior emissivity of the absorber
: | The inner emissivity of the glass envelope
: | The thermal diffusivity of the space gas, m²/s
: | The coefficient of cubical expansion, 1/°C
: | The gas conductivity factor in the current condition
: | The coefficient of kinematics viscosity of the space gas, m²/s
: | The acceleration of gravity, m/s²
: | The inner diameter of the absorber, m
: | The exterior diameter of the absorber, m
: | The inner diameter of the glass envelope, m
: | The exterior diameter of the glass envelope, m
: | The length of the collector, m
: | The area of the whole collector, m²
| The Prandtl number of the space gas
| The Rayleigh number of the absorber based on the exterior surface diameter
: | The solar direct radiation, W/m²
: | The incidence angle of the sunlight, degree
: | The modifier of the incidence angle
: | The shadow barrier factor
: | The end loss factor
: | The optics efficiency of collectors’ mirrors
: | The optics efficiency of the glass envelope and absorber
: | The coefficient of kinematics viscosity of the fluid (water or steam) in average fluid temperature, m²/s
: | The coefficient of kinematics viscosity of the fluid (water or steam) in absorber temperature, m²/s
: | The factor of operating and tracking.

Acknowledgments

The authors are grateful for the Program of China National Natural Science Fund Project (no. 51106048) and National Basic Research Program of China, 973 Project (no. 2009CB219801).

References

T. A. Stuetzle, Automatic control of the 30MWe SEGS VI parabolic trough plant [M.S. thesis], The University of Wisconsin-Madison, Madison, Wis, USA, 2002.

T. Stuetzle, N. Blair, J. W. Mitchell, and W. A. Beckman, “Automatic control of a 30 MWe SEGS VI parabolic trough plant,” Solar Energy, vol. 76, no. 1–3, pp. 187–193, 2004.

A. M. Patnode, Simulation and performance evaluation of parabolic trough solar power plants [M.S. thesis], The University of Wisconsin-Madison, Madison, Wis,USA, 2006.

http://www.therminol.com/pages/products/vp-1.asp

S. D. Odeh, G. L. Morrison, and M. Behnia, “Modelling of parabolic trough direct steam generation solar collectors,” Solar Energy, vol. 62, no. 6, pp. 395–406, 1998.

M. M. Rolim, N. Fraidenraich, and C. Tiba, “Analytic modeling of a solar power plant with parabolic linear collectors,” Solar Energy, vol. 83, no. 1, pp. 126–133, 2009.

O. García-Valladares and N. Velázquez, “Numerical simulation of parabolic trough solar collector: improvement using counter flow concentric circular heat exchangers,” International Journal of Heat and Mass Transfer, vol. 52, no. 3-4, pp. 597–609, 2009.

R. Forristall, Heat Transfer Analysis and Modeling of a Parabolic Trough Solar Receiver Implemented in Engineering Equation Solver, National Renewable Energy Lab., Denver, Colo, USA, 2003.

H. Price, E. Lüpfert, D. Kearney et al., “Advances in parabolic trough solar power technology,” Journal of Solar Energy Engineering, Transactions of the ASME, vol. 124, no. 2, pp. 109–125, 2002.