Abstract

The scope of this study is to present a novel methodology for the validation of a solar parabolic trough collector simulation model, developed in System Advisor Model software, using solar field thermal performance experimental model functions that duly account for the contributions of errors that are inherent in the collection of input experimental data. The study focuses on the uncertainty associated with the predicted thermal energy output and thermal efficiency, which were calculated by using standardized testing procedures that were developed by the National Renewable Energy Laboratory and the American Society of Mechanical Engineers. Errors due to imperfections of the experimental models used throughout the test are also counted in, as well as uncertainties attributed to the variability of meteorological conditions. The experimental tests were conducted under clear-sky and steady-state conditions in Kenya, using a prototype parabolic trough collector. The uncertainty analysis provided a realistic evaluation of the thermal performance of the prototype during testing, resulting in expanded uncertainties of 9.05% (0.104 kWhth) for thermal output and 3.66% (0.0258) for thermal efficiency. Notably, the predicted thermal output and thermal efficiency from the experimental models did not exceed observed levels, demonstrating a strong correlation between predicted and observed values, supported by R2 regression coefficients of 0.972 for thermal output and 0.989 for thermal efficiency. The comparison of the two experimental model results with the simulation outcomes validated the performance of the simulation model, as the simulation results fell within the experimental error margins. Additionally, the statistical analysis yielded significant results, with root mean square error, mean bias error, and t-statistics values for thermal energy output of 0.057, 0.033, and 2.135, respectively, and for thermal efficiency, 0.018, 0.010, and 1.993, respectively, indicating the accuracy and reliability of the simulation model. Thus, the simulation model has been successfully validated, proving its capability to accurately predict the thermal output and efficiency of the parabolic trough collector.

1. Introduction

In recent years, industrial activities have played a significant role in the economic growth of several countries, resulting in substantial increases in productivity and improved living standards. However, the expansion of these industrial activities has led to a notable surge in global energy consumption [1, 2], while contributing 40% of global carbon dioxide (CO2) emissions [3]. In Kenya, the industrial sector relies significantly on imported petroleum and wood fuel for heat energy. The tea factories, in particular, depend heavily on the burning of wood fuel, which has become increasingly expensive, to generate process heat for tea drying. This reliance on wood fuel also leads to the emission of greenhouse gases (GHG).

Renewable energy offers a significant opportunity to reduce pollution and address the environmental imbalances associated with biomass in the tea sector [1]. Solar energy, in particular, holds great promise for addressing multiple challenges, including the growing demand for energy [4], global warming [5], and environmental pollution [6]. In the industrial sector, a 52% share of heat is utilized within low and medium temperature range (<400°C) [7]. This is a temperature regime that is particularly accessible for commercially available solar–thermal technologies [8]. The parabolic trough collector (PTC) has demonstrated its potential as an excellent choice for industrial applications, establishing itself as the most promising and well-developed concentrated solar–thermal (CST) technology [9]. PTC technology currently dominates the market, with an 81% share [10], and can be a better technological alternative to the traditional wood-fired steam boilers in the tea industry in Kenya. Even though it has considerable potential for industrial heat, the development of the PTC field is still very limited, particularly in regions with high solar potential such as Africa.

Modeling is a powerful tool at the level of technical improvements that provides an efficient means for characterizing the performance of PTC systems by providing an approximation of their behavior [1, 11]. Models or dynamic simulation tools, such as the System Advisor Model (SAM), which was developed by the National Renewable Energy Laboratory (NREL) in the USA, can be used for assessing the performances of PTC systems under variable conditions. Other simulation software in solar–thermal applications include the Green Energy System Analysis Tool (Greenius) that was developed by the Institute of Solar Research (DLR) in Germany. This software has been used in several studies; for instance, one study [12] used it to simulate the performances of different heat-transfer fluids (HTFs) for a 1-MWe parabolic trough in Oujda (Eastern Morocco) and another study [13] used it to simulate the performances of a parabolic trough in Maan (Jordan). Transient System (TRNSYS) is another popular simulation tool that is used to simulate the behavior of transient systems, including CST technologies. For instance, one study [14] validated a TRNSYS model for a thermodynamic plant with parabolic trough solar–thermal power in Ain-Témouchent (Algeria), while another study [15] used is to analyse the performance of a 30-MWe SEGS VI parabolic trough plant in a Californian desert (USA). Nevertheless, for solar–thermal systems, SAM is the preferred option as it provides a wide range of performance metrics that are specific to solar–thermal systems, enabling users to gain insights into various aspects of system performance, including thermal output, efficiency, and more. SAM also allows a high degree of customization in defining system configurations, making it adaptable to different technologies and project requirements [16, 17].

There is no doubt that an accurate model for solar prediction allows improved designs of energy systems [18, 19]. In the context of the current rapid development of large-scale solar energy projects, modeling allows the saving of time and effort while providing a significant amount of information. Yet, the accuracy of the modeled radiation datasets is of utmost importance in providing a high degree of confidence to solar project developers while demonstrating the feasibility of industrial thermal PTC projects in countries such as Kenya, where CST activities are still very nascent. This calls for careful validation, normally against high-quality measurements [20]. Validation can be defined as the process of building an acceptable level of confidence that an inference about a simulated process is correct or valid inference for the actual process [21]. It is a critical step in establishing the credibility of a simulation model among its users [22]. It is therefore necessary to subject the models to experimental facilities or model prototypes to permit the validation of the results.

However, few studies have been published in the literature regarding the development of outdoor validation tests [23] or even prototypes for the purposes of PTC models validation. Some of the studies include, for instance, the work by Schiricke et al. [24], who made one of the first attempts at the experimental validation of CST systems by analysing the optical modeling of PTCs that was developed by using commercial Monte Carlo ray-tracing software. The experimental measurements confirmed the reliability of the model. Alfellag [25] presented an experimental investigation of a PTC that was developed in MATLAB® software by using a test rack containing a circulating pump, a storage tank, thermocouples, a heat exchanger, and a flow meter. The results showed acceptable agreement even though there were some variances. Mateo et al. [26] conducted the experimental validation of a mathematical model of a PTC for a solar adsorption refrigerator using experimental data gathered from similar collectors in two different locations. The model was satisfactorily validated by experimental data. Similarly, Soudani et al. [27] experimentally validated a mathematical model of a PTC and registered close results between theoretical and experimental data, with the difference in thermal efficiency recorded as being <6%. Mouaky et al. [2] studied a PTC analytical model and validated it by using experimental data from a PTC test loop, with the result that the model used for PTC simulation was accurate and the daily average deviation was ~4.8% under clear-sky conditions. In 2020, Boretti et al. [28] conducted the first validation of SAM for both PTC and solar tower models by using high-frequency data collected on existing facilities and concluded that computational results were relatively close to measured values for the PTC but very far away for solar tower systems.

It is evident from the review of the literature that most CST validation assessments that employ the use of experimental results rarely account for errors from the most influential factors in the PTC field measurement parameters in order to increase the accuracy of the measured data. This leads to less reliable field data output results. Extensive analysis of the uncertainties associated with experimental data—also referred to as acceptance tests—helps to address and quantify the contributions of errors that are inherent in the collection of input data. CST performance is subject to a number of uncontrollable effects, including solar position, level of direct normal irradiation (DNI), wind velocity, and ambient temperature. While the effect of some variables such as the DNI is unquestionably more prominent than those of others, each significant effect needs to be quantified in order to obtain a sufficiently precise prediction of solar field performance. Because no measurement technique can perfectly determine the quantity it measures (every measurement device has inherent uncertainty associated with the measurement values), the goal of uncertainty analysis is to determine whether the performance of the designed plant model in question will meet or exceed the target, given a particular uncertainty and required confidence level [29]. This is because a realistic assessment of model suitability is not possible without uncertainty analysis [30].

This paper seeks to demonstrate a novel approach to the validation of a solar PTC model developed in SAM for industrial process heat application in tea drying in Kenya by using experimental data from a PTC prototype. The novel validation approach of the SAM model begins with quantification of the uncertainties in the experimental data and further establishes thermal energy output and efficiency model equations that are then utilized for validation of the PTC simulation model. The validation equations are derived from a data-fitting technique through linear regression based on the least squares method. Experimental data are derived from performance measurements of the PTC model prototype that was designed, fabricated, and tested in Kajiado, Kenya. This is the first study in Kenya, to the best of the authors’ knowledge, to have used this approach to perform the validation of a PTC plant simulation model.

2. System description and experimental test set-up

A parabolic trough plant simulation model for industrial process heat application was initially developed in SAM by using meteorological and process heat demand data from the Toror tea factory in Kericho, Kenya [1]. A small-scale PTC prototype was built for the purpose of validating the predictions of this model—specifically the solar field component. The system was installed at a site in Kimuka Kajiado, Kenya (latitude 1.418426° S and longitude 36.59749° E) for experimental testing under outdoor weather conditions, using water as the system HTF.

A collection of governing equations is employed in determining the design parameters for the prototype. The initial parameter under consideration is the curve length S of the reflective surface, which is identical to the width of the aluminum reflector mirror, which is 1150 mm. By using this information, we can calculate the latus rectum of the parabola, which, in turn, will assist in determining the remaining parameters. The formula for calculating the curve length is derived from the work of Kalogirou and it is expressed as follows [31]:

(1)

where Hp represents the latus rectum of the parabola, which corresponds to the width of the parabolic shape at its focal point, and φr denotes the rim angle, which is the angle formed between the center line of the collector and the line connecting the outer rim of the reflector to the focal point [32].

It is possible to obtain different rim angles for the same aperture; for the different rim angles, the focus-to-aperture ratio, which defines the curvature of the parabola, changes. The optimum rim angle is between 70° and 110° [33]. Values of <70° reduce the aperture width and the collector working temperature while values of >110° will increase the reflective area without increasing the effective area [32]. Kalogirou also demonstrated that, for a 90° rim angle, the mean focus-to-reflector distance, and hence the reflected beam spread, is minimized so that the slope and tracking errors are also minimized, leading to higher optical efficiency [31]. For these reasons, this study adopts a rim angle of 90°. Therefore, Hp was established to be ~1000 mm by using Equation (1).

Another important parameter related to the rim angle φr is the aperture of the parabola Wa, which is given by the following equation [34]:

(2)

where f is the focal length of the parabola, derived from Equation (3):

(3)

The parabola latus rectum Hp is equal to Wa when φr = 90° [31]. Therefore, the value of Wa is equal to 1000 mm and, from Equation (3), our focal length f is 250 mm, as illustrated in Fig. 1.

The characteristic dimensions of the fabricated PTC.
Figure 1.

The characteristic dimensions of the fabricated PTC.

Another important parameter is the collector geometrical concentration ratio Cg, defined as the ratio of the effective area of the aperture Aa to the area of the receiver Ar, given by [35]:

(4)

where Dout is the external diameter of the glass envelope, dout is the external diameter of the absorber tube, and Lc and Lr are the lengths of the collector and the receiver, respectively.

In this research, the external diameter of the glass envelope is measured as 95 mm, while the external diameter of the absorber tube is 40 mm. The length of the absorber is consistent with that of the collector mirror, which is 2000 mm. Consequently, the concentration ratio, calculated as the ratio of the mirror width to the absorber diameter, was calculated to be 14.41. This value is considered desirable, as it exceeds 10, as recommended by [36].

The manufacturer provides several receiver parameters, including the absorber tube transmittance, absorptance, and glass envelope surface absorptance. For other parameters, such as the glass envelope surface emittance [37] and the absorber tube surface emittance [38], material composition data from the manufacturer were used to obtain values from existing literature. Table 1 lists both the selected and computed parameters for the collector mirror and absorber tube, while Fig. 2 shows the fabricated PTC prototype.

Table 1.

Design parameters and considerations.

ParameterValue
Collector
 Reflective aperture area, Aa1.81 m2
 Length of collector, Lc2.00 m
 Width of collector (curve length), S1.15 m
 Aperture of parabola, Wa1.00 m
 Focal length, f0.25 m
 Latus rectum, Hp1.00 m
 Rim angle, φr90°
 Mirror reflectance, 𝜌m0.95
 Geometric concentration ratio, Cg14.41
Receiver tube
 Length of receiver tube, Lr2.0 m
 Absorber tube internal diameter, din0.037 m
 Absorber tube external diameter, dout0.040 m
 Glass envelope internal diameter, Din0.091 m
 Glass envelope external diameter, Dout0.095 m
 Absorber tube surface emittance, abs0.32
 Absorber tube surface absorptance, αabs0.96
 Glass envelope surface emittance, env0.88
 Glass envelope surface absorptance, αenv0.02
 Glass envelope transmittance, 𝜏env0.92
ParameterValue
Collector
 Reflective aperture area, Aa1.81 m2
 Length of collector, Lc2.00 m
 Width of collector (curve length), S1.15 m
 Aperture of parabola, Wa1.00 m
 Focal length, f0.25 m
 Latus rectum, Hp1.00 m
 Rim angle, φr90°
 Mirror reflectance, 𝜌m0.95
 Geometric concentration ratio, Cg14.41
Receiver tube
 Length of receiver tube, Lr2.0 m
 Absorber tube internal diameter, din0.037 m
 Absorber tube external diameter, dout0.040 m
 Glass envelope internal diameter, Din0.091 m
 Glass envelope external diameter, Dout0.095 m
 Absorber tube surface emittance, abs0.32
 Absorber tube surface absorptance, αabs0.96
 Glass envelope surface emittance, env0.88
 Glass envelope surface absorptance, αenv0.02
 Glass envelope transmittance, 𝜏env0.92
Table 1.

Design parameters and considerations.

ParameterValue
Collector
 Reflective aperture area, Aa1.81 m2
 Length of collector, Lc2.00 m
 Width of collector (curve length), S1.15 m
 Aperture of parabola, Wa1.00 m
 Focal length, f0.25 m
 Latus rectum, Hp1.00 m
 Rim angle, φr90°
 Mirror reflectance, 𝜌m0.95
 Geometric concentration ratio, Cg14.41
Receiver tube
 Length of receiver tube, Lr2.0 m
 Absorber tube internal diameter, din0.037 m
 Absorber tube external diameter, dout0.040 m
 Glass envelope internal diameter, Din0.091 m
 Glass envelope external diameter, Dout0.095 m
 Absorber tube surface emittance, abs0.32
 Absorber tube surface absorptance, αabs0.96
 Glass envelope surface emittance, env0.88
 Glass envelope surface absorptance, αenv0.02
 Glass envelope transmittance, 𝜏env0.92
ParameterValue
Collector
 Reflective aperture area, Aa1.81 m2
 Length of collector, Lc2.00 m
 Width of collector (curve length), S1.15 m
 Aperture of parabola, Wa1.00 m
 Focal length, f0.25 m
 Latus rectum, Hp1.00 m
 Rim angle, φr90°
 Mirror reflectance, 𝜌m0.95
 Geometric concentration ratio, Cg14.41
Receiver tube
 Length of receiver tube, Lr2.0 m
 Absorber tube internal diameter, din0.037 m
 Absorber tube external diameter, dout0.040 m
 Glass envelope internal diameter, Din0.091 m
 Glass envelope external diameter, Dout0.095 m
 Absorber tube surface emittance, abs0.32
 Absorber tube surface absorptance, αabs0.96
 Glass envelope surface emittance, env0.88
 Glass envelope surface absorptance, αenv0.02
 Glass envelope transmittance, 𝜏env0.92
Solar PTC prototype.
Figure 2.

Solar PTC prototype.

3. Experimental methodology and validation approach

This study adopts four significant steps that are necessary for successful model validation by using experimental data, namely:

  • analysis and quantification of the uncertainties associated with experimental data;

  • development of experimental model functions with regression coefficients determined by using the least squares method;

  • evaluation of the quality of fit of the model functions;

  • validation of the PTC model by using thermal output and efficiency experimental model functions.

3.1 Approach to uncertainty analysis

The parabolic trough uncertainty analysis, or performance acceptance test, aims to determine whether the measured performance metrics meet the guaranteed results. To achieve this, accurate data measurement tools are crucial. However, any real measurement system is prone to both random and systematic errors, which can be attributed to various sources such as calibration errors, instrument inaccuracies, and methodological limitations [29].

Both NREL and the American Society of Mechanical Engineers (ASME) have developed guidelines for uncertainty analysis in parabolic trough systems. These guidelines provide recommendations for procedures that can yield accurate results, consistently with good engineering practices. The guidelines are specifically designed for PTC systems with high-temperature synthetic oil but the principles are applicable to other CST systems [39]. Additionally, the international testing standard (ISO 9806:2013) provides guidance on uncertainty analysis in low-temperature collectors.

This study adopts the approach presented by Kearney [39] and Wagner et al. [29] to carry out the uncertainty analysis. This approach adopts both NREL and ASME procedures. Generally, the following steps are observed:

  • (i) Define the uncertainty problem by identifying the specific uncertainty aspect through the parameters to be analysed. In this study, uncertainties in the solar field thermal output and thermal efficiency are selected for uncertainty analysis, as the same have been selected as the validation parameters.

  • (ii) Identify and quantify various sources of uncertainty, such as measurement errors and parameter variations.

  • (iii) Formulate mathematical models of validation parameters to propagate the uncertainties through the system, accounting for correlations between variables.

  • (iv) Measure and evaluate the performance metrics (e.g. direct normal irradiance, reduced temperature difference, mass flowrates) for each mathematical model.

  • (v) Analyse and interpret the results to understand the impact of each source of uncertainty on the overall system performance.

3.2 Measurement parameters and instruments

The measurements required for the uncertainty tests are the HTF inlet and outlet temperatures, the HTF flow rate, and real-time weather data, including ambient temperature, wind velocity, and direct normal irradiance. The real-time weather data are recorded on-site by using a digital weather station equipped with a color light-emitting diode display. This station includes a solar panel, transmission module, wind vane, high-speed anemometer, thermo-hygrometer, and a rain funnel. It measures the ambient temperature (–40°C to 60°C), wind direction and speed (0–45 m/s), atmospheric pressure (600–1100 hPa), humidity (1%–95%), and rainfall levels (0–9999 mm).

For SAM, the key input data include the ambient temperature, wind speed, wind direction, relative humidity, atmospheric pressure, and solar irradiance. For the PTC prototype, essential measurements are the inlet and outlet temperatures of the collector, as well as the mass flow rate. The mass flow rate and inlet temperature were monitored by using a digital liquid crystal display on a water flow sensor meter, which operates within a range of 1–25 liters/min (0.017–0.42 kg/s). The collector outlet temperature was measured by using a GM320 infrared digital thermometer that was capable of measuring temperatures from –50°C to 400°C. Global irradiance was recorded by using an EKO MS-602, while diffuse radiation was measured by using another MS-602 pyranometer with a shading ball to block direct sunlight. Consequently, the direct normal irradiance was calculated by using the relation [40]:

(5)

where Gg, Gd, and Gb are the global, diffuse, and direct normal radiation, respectively, while θz is the zenith angle of the Sun.

Table 2 lists the types of parameters measured and the accuracy of each instrument used in the solar field measurements.

Table 2.

Parameters and instrument uncertainties.

ParameterInstrument uncertainty
HTF temperature±1.5°C
Flow rate±0.1%
Direct normal irradiance±2.0%
Ambient temperature±0.1°C
Relative humidity±1.0%
Wind speed±0.5 m/s
Atmospheric pressure±3 hPa
ParameterInstrument uncertainty
HTF temperature±1.5°C
Flow rate±0.1%
Direct normal irradiance±2.0%
Ambient temperature±0.1°C
Relative humidity±1.0%
Wind speed±0.5 m/s
Atmospheric pressure±3 hPa
Table 2.

Parameters and instrument uncertainties.

ParameterInstrument uncertainty
HTF temperature±1.5°C
Flow rate±0.1%
Direct normal irradiance±2.0%
Ambient temperature±0.1°C
Relative humidity±1.0%
Wind speed±0.5 m/s
Atmospheric pressure±3 hPa
ParameterInstrument uncertainty
HTF temperature±1.5°C
Flow rate±0.1%
Direct normal irradiance±2.0%
Ambient temperature±0.1°C
Relative humidity±1.0%
Wind speed±0.5 m/s
Atmospheric pressure±3 hPa

3.3 Test conditions

A reliable uncertainty analysis requires testing under steady-state conditions and thermal equilibrium. However, the constantly changing energy input from the Sun and other factors in the solar field make it difficult to achieve this goal. Failing to achieve thermal equilibrium or steady-state test conditions can contribute to uncertainty in the test results. To minimize this uncertainty, it is essential to keep variations in key test parameters low so that they do not significantly impact the uncertainty band in the results. This can be achieved by ensuring that the solar system is under stable test conditions and in a thermal equilibrium state before testing begins [41].

Taking into account the impact of variations in test parameters on the overall uncertainty of the test results, a set of stabilization criteria and test conditions has been established to ensure a steady state for PTC energy tests, as presented below [29, 39].

Stabilization criteria are applicable to four parameters with allowable variability over the test period provided in brackets below:

  • HTF volumetric flow rate (≤0.5%)

  • Direct normal irradiance (≤0.5%)

  • HTF inlet temperature (≤0.2%)

  • HTF outlet temperature (≤0.2%)

The steady-state test, on the other hand, is defined by the following conditions:

  • Tests should be carried out between 9 a.m. and 4 p.m.

  • Duration of test run should be ~15–30 minutes with test-run data points collected every 10 seconds on average. Therefore, the number of test points in each test run should be between 90 and 180

  • Maximum wind speed should be ≤13 m/s

  • DNI on the solar field area has to be >500 W/m2 with DNI variations smaller than ±20 W/m2

In the experimental validation of PTCs, steady-state conditions are established to ensure that the system reaches stable and uniform operating conditions, allowing accurate measurement and analysis of the performance of the collector. The steady-state conditions during experimental validation ensures that results accurately represent the long-term performance of the collector under stable operating conditions. This is crucial for the design and optimization of PTCs for practical applications.

The given stabilization criteria and test conditions have been chosen in such a way that the variations in the key test parameters are low enough to contribute in only a minor way to the uncertainty band in the results. At the same time, the set of stabilization criteria parameters were settled upon because their variability has a strong influence on the total uncertainty of the test results. The tests were conducted by using a short duration of 10-second time steps to capture the impact of transient effects resulting from a sudden change in the solar field inlet temperature, otherwise the final stabilization criteria for a specific project will be strongly influenced by the design of the solar system and associated instrumentation.

3.4 Calculation of measurement uncertainty

Due to the natural variability of resources and the imperfections in control systems, variations in all measured parameters are unavoidable. Uncertainties are typically classified into two types: systematic errors and random errors. While random errors tend to average out when multiple measurements are taken over a period of time, systematic errors persist regardless of the number of measurements. Systematic errors can arise from various sources, including the calibration process, instrument malfunctions, transducer errors, and fixed errors in the measurement method. Similarly, random errors can be based on the manufacturer’s specifications. For repeated measurements, the random standard uncertainty can be defined by:

(6)

where 𝜎X is the standard deviation of a series of sampled data and N is the number of data points collected over the test interval.

Calculated results, such as the delivered thermal energy and the solar–thermal efficiency, are not typically measured directly but rather are based on parameters that are measured during the course of one or multiple acceptance tests. Equation (7) describes how the uncertainties in each of the measured variables X propagate into the value of a calculated resulting quantity R. In this case, the result R is a function of the individual or average values of these independent parameters:

(7)

where the subscript i describes the number of parameters used in the calculation of the result and X¯ is either the value of a single measurement of the parameter or the average value of the parameter based on a number of N repeated measurements.

The expression for the combined standard measurement uncertainty of a calculated result based on multiple error sources can in many cases be calculated from the root-sum-square of the total uncertainty of the individual random and systematic error sources, as shown in Equation (8):

(8)

where SR is the standard random uncertainty component of a result and bR is the systematic standard uncertainty component of a result, calculated as follows:

(9)
(10)

where SX¯i is the random standard uncertainty of the mean of N measurements (Equation (6)) and bX¯i is defined as the systematic standard uncertainty component of a parameter, calculated as follows:

(11)

where 𝛼 is the estimated accuracy of the measurement or the instrument uncertainty. When the uncertainty is presented as a percentage, the systematic standard uncertainty is the product of the measurement uncertainty (as a percentage) and the nominal measurement [28].

The combined standard uncertainty in Equation (8) implies that the calculated result will capture the true result within a 68% confidence level (or 1 standard deviation). Typically, a confidence level of 95% (or 2 standard deviations) is desirable. Therefore, the expanded uncertainty in the result is given by:

(12)

3.5 Thermal energy output and efficiency calculations

Once the thermal equilibrium and test condition stability have been established, the criteria for thermal energy output and efficiency measurements are primarily based on the level of uncertainty in the test results that are calculated by using the equations discussed under section 3.4 on ‘Calculation of measurement uncertainty’.

The solar field thermal output, or useful heat gain, of the prototype collector is determined from the measured inlet and outlet temperatures, mass flow rate, and specific heat of water (HTF) by using the following equation:

(13)

where Tin and Tout are the inlet and the outlet temperatures, respectively, m˙ is the mass flow rate, and Cp is the specific heat of the water.

The partial derivative terms symbolically represent the manner in which the result R varies with the changing independent variable. For example, the response of the thermal energy output to the mass flow rate given Equation (13) is:

(14)

By extending the principle shown in Equation (14) to the entire systematic uncertainty calculation for solar field thermal energy output, the absolute standard random uncertainty are given in Equations (15) and (16):

(15)
(16)

Similarly, the absolute standard systematic uncertainty can be derived as shown in Equation (16).

(17)

The solar field thermal efficiency is evaluated experimentally by using the equation below [33]:

(18)

where Aa is the effective aperture area of the collector and Gb is the direct solar irradiance component in the aperture plane of the collector.

The methodology for obtaining Equations (16) and (17) can be applied identically to the solar field efficiency in Equation (18) to arrive at the estimated uncertainties. The resulting equations for the absolute standard random systematic uncertainty and absolute standard systematic uncertainty associated with the solar field efficiency are presented in Equations (19) and (20), respectively:

(19)
(20)

3.6 Experimental model functions and data fitting

Data fitting is the process of utilizing the uncertainty analysis to develop experimental model functions with parameters that represent observed data with the greatest accuracy. A model function depicts the underlying functional relationship between the outcome of the measurement and the corresponding conditions, whether known or assumed. Therefore, the goal of data fitting is to find those parameters of the experimental model function that best describe the connection between observations and conditions [42].

The least squares method is a widely used approach in statistics and numerical analysis to find the best-fitting parameters of experimental model functions to a given set of data points. The method aims to minimize the sum of the squared differences between the observed data points and the corresponding values that are predicted by the experimental model [43]. In other words, any least squares curve- or line-fitting algorithm optimizes the constants of a fitting equation by minimizing the sum of the squares of the deviations of the actual (data) values from the values that are predicted by the equation [44].

In this study, we develop two experimental model functions that are later used for validation of the SAM PTC simulation model. These are the ‘solar–thermal output model function’ and the ‘thermal efficiency model function’.

3.7 Solar–thermal output model function

The model functions are derived by using the following procedure:

(i) The parameters for the steady-state model equation for the solar–thermal output of the collector prototype is established. The solar–thermal output model is a linear three-parameter function, adapted from [45–47]:

(21)

where Qgain is the solar–thermal output or dependent variable, and is a function of the direct normal irradiance Gb and the difference between the average collector temperature Tm and the ambient temperature Ta, denominated as ΔT (in this simplified model, the average collector temperature Tm is calculated by using (Tin + Tout)/2). The parameter 𝛼1 is the universal constant or y-intercept, otherwise known as the general mean [48], while 𝛼2 and 𝛼3 are the regression coefficients measuring the average effect of the respective meteorological (independent) variables, in this case Gb and ΔT(Tm – Ta) respectively.

Equation (21) is the fitted trend line and can be used to forecast the future solar energy output.

(ii) The next step is to find the linear regression coefficients 𝛼1, 𝛼2, and 𝛼3 together with their standard uncertainties by computing a least squares solution of Equation (21) as explained by [49]. That is, Equation (21) is treated as a system of equations in the unknowns 𝛼1, 𝛼2, and 𝛼3. In matrix form, this can be written as Ax = b, where:

As a least squares solution, x^ provides the best-fit parameters of Equation (21) and is determined by multiplying both sides of Ax = b by the transpose to form an augmented matrix (ATA. ATb), which is then row reduced to obtain x^:

The standard uncertainties of the parameters are obtained from the matrix (ATA)–1 where the diagonal elements are the squared uncertainties of the regression parameters and the off-diagonal elements are the co-variances between them.

(iii) Once the regression parameters are established, it is important to test the goodness-of-fit of the solar energy output model function. A satisfactory measure of the “goodness-of-fit” of the model is crucial to deem the combined standard uncertainties in model parameters trustworthy and to validate the suitability of employing the model [29]. A number of statistical indicators are available for evaluating the model functions that govern the experimental behavior of measured data. In this study, the coefficient of determination, or regression, R2 is used for evaluating the goodness-of-fit.

R2 represents the proportion of the variance in the dependent variable that is predictable from the independent variable(s). It is therefore a measure of how well the regression line approximates and fits the real data points. Hence, a coefficient of determination with a value of 1 indicates that the regression line perfectly fits the data [50]. It can be estimated by using Equation (22) [51]:

(22)

where n is the number of data points, and Xi and Yi are the model and measured values, respectively, at time step i.

3.8 Thermal efficiency model function

The thermal efficiency model is a three-parameter quadratic function adapted from [52]:

(23)

where the quantity T* is referred to as the reduced temperature difference and is expressed as (Tm – Ta)/Gb.

The parameters U0, U1 and U2 that best fit the collector efficiency model function are the solutions of a multiple linear regression problem [53]. The model function is derived by using a similar procedure to that for solving the solar–thermal output model function, discussed above.

4. Results and discussion

4.1 Test conditions

Tests were conducted on a clear-sky day (31 March 2024) and the period that was selected for this analysis falls between 12 p.m. and 12:15 p.m. (East African Time). This time period was selected because it corresponds to a relatively steady solar resource level and the incident conditions of almost solar noon (at 12:00 solar time). This also allows the tests to be performed under almost the same incidence angle values. The data readings were 10 seconds apart for a total of 90 readings over the 15-minute test duration. The test duration of 15 minutes was selected based on the data-reading time step of 10 seconds and the number of test readings that enables the acceptable uncertainty test conditions that are described under Section 3.1.2 to be met. The irradiance during the tests was in the range of between 890 and 908 W/m2, resulting in DNI variations of only 18 W/m2. The maximum wind speed and relative humidity recorded were 2.6 m/s and 22%, respectively. The mass flow rate was kept constant at 0.09 kg/s. Table 3 shows the comparison between the allowable steady-state variability of the four critical test parameters and the results from this study. The study has met the stabilization criteria for a steady-state performance test.

Table 3.

Stabilization criteria for short-duration steady-state tests.

ParameterAllowable variability over test period
SX¯/X¯(%)
Actual test variability over test period
SX¯/X¯(%)
HTF flow rate, kg/s0.5%0.00%
Direct normal irradiance, W/m20.5%0.05%
Solar field inlet temperature, °C0.2%0.14%
Solar field outlet temperature, °C0.2%0.18%
ParameterAllowable variability over test period
SX¯/X¯(%)
Actual test variability over test period
SX¯/X¯(%)
HTF flow rate, kg/s0.5%0.00%
Direct normal irradiance, W/m20.5%0.05%
Solar field inlet temperature, °C0.2%0.14%
Solar field outlet temperature, °C0.2%0.18%
Table 3.

Stabilization criteria for short-duration steady-state tests.

ParameterAllowable variability over test period
SX¯/X¯(%)
Actual test variability over test period
SX¯/X¯(%)
HTF flow rate, kg/s0.5%0.00%
Direct normal irradiance, W/m20.5%0.05%
Solar field inlet temperature, °C0.2%0.14%
Solar field outlet temperature, °C0.2%0.18%
ParameterAllowable variability over test period
SX¯/X¯(%)
Actual test variability over test period
SX¯/X¯(%)
HTF flow rate, kg/s0.5%0.00%
Direct normal irradiance, W/m20.5%0.05%
Solar field inlet temperature, °C0.2%0.14%
Solar field outlet temperature, °C0.2%0.18%

The prototype collector system was set up and tested approximately 1 hour before the performance test began at 12 p.m. This precaution was taken to ensure that heat loss by the system remained relatively consistent throughout the 15-minute experimental period. Otherwise, if the solar field is not in equilibrium—for instance, the temperatures of the solar field HTF and the piping are still changing over time—the delivered energy will be less than the absorbed power because some of the excess energy is used to heat the solar field HTF and piping. The largest HTF inlet and outlet temperature difference observed during the 15-minute test was 3.3°C, with the highest recorded inlet and outlet temperatures being 43.8°C and 47.1°C, respectively. The system was operated at lower temperatures to minimize heat-transfer fluid losses and achieve optimal efficiency.

4.2 Uncertainty evaluation

Tables 4–7 summarize the uncertainty data and results that were derived from the steady-state tests based on the methodology and equations discussed earlier as applied to the solar field energy output and efficiency calculations, respectively. These results provide the uncertainty contributions of the various measured data parameters over the course of the acceptance test time period. The tables are presented in the standard ASME PTC format for reporting uncertainty. The primary sources of uncertainty in the experimental results are attributed to the measurement of the solar irradiance and temperature difference. To improve the accuracy of the data, it is essential to focus on precise measurements of these variables, as they are critical factors in determining the thermal output and efficiency of the solar PTC.

Table 4.

Solar field energy output uncertainty for measured data.

Parameter information (in parameter units)Uncertainty contribution of parameters to the result (in units squared)
SymbolDescriptionUnitsNominal valueStandard deviation
σi
Number of measurements
Ni
Absolute systematic standard uncertainty
bX¯i
Absolute random standard uncertainty
SX¯i
Absolute sensitivity
RX¯i
Absolute systematic standard uncertainty contribution
(RX¯ibX¯i)2
Absolute random standard uncertainty contribution
(RX¯iSX¯i)2
m˙Mass flow ratekg/s0.090.0000900.00090.000012.601.29 × 10–4########
CpHTF specific heatkJ/kg °C4.200.0000900.00000.00000.270################
ToutHTF outlet temperature°C46.10.4189900.05800.04420.3784.81 × 10–42.79 × 10–4
TinHTF inlet temperature°C43.10.2280900.05800.02400.3784.81 × 10–48.25 × 10–5
Parameter information (in parameter units)Uncertainty contribution of parameters to the result (in units squared)
SymbolDescriptionUnitsNominal valueStandard deviation
σi
Number of measurements
Ni
Absolute systematic standard uncertainty
bX¯i
Absolute random standard uncertainty
SX¯i
Absolute sensitivity
RX¯i
Absolute systematic standard uncertainty contribution
(RX¯ibX¯i)2
Absolute random standard uncertainty contribution
(RX¯iSX¯i)2
m˙Mass flow ratekg/s0.090.0000900.00090.000012.601.29 × 10–4########
CpHTF specific heatkJ/kg °C4.200.0000900.00000.00000.270################
ToutHTF outlet temperature°C46.10.4189900.05800.04420.3784.81 × 10–42.79 × 10–4
TinHTF inlet temperature°C43.10.2280900.05800.02400.3784.81 × 10–48.25 × 10–5
Table 4.

Solar field energy output uncertainty for measured data.

Parameter information (in parameter units)Uncertainty contribution of parameters to the result (in units squared)
SymbolDescriptionUnitsNominal valueStandard deviation
σi
Number of measurements
Ni
Absolute systematic standard uncertainty
bX¯i
Absolute random standard uncertainty
SX¯i
Absolute sensitivity
RX¯i
Absolute systematic standard uncertainty contribution
(RX¯ibX¯i)2
Absolute random standard uncertainty contribution
(RX¯iSX¯i)2
m˙Mass flow ratekg/s0.090.0000900.00090.000012.601.29 × 10–4########
CpHTF specific heatkJ/kg °C4.200.0000900.00000.00000.270################
ToutHTF outlet temperature°C46.10.4189900.05800.04420.3784.81 × 10–42.79 × 10–4
TinHTF inlet temperature°C43.10.2280900.05800.02400.3784.81 × 10–48.25 × 10–5
Parameter information (in parameter units)Uncertainty contribution of parameters to the result (in units squared)
SymbolDescriptionUnitsNominal valueStandard deviation
σi
Number of measurements
Ni
Absolute systematic standard uncertainty
bX¯i
Absolute random standard uncertainty
SX¯i
Absolute sensitivity
RX¯i
Absolute systematic standard uncertainty contribution
(RX¯ibX¯i)2
Absolute random standard uncertainty contribution
(RX¯iSX¯i)2
m˙Mass flow ratekg/s0.090.0000900.00090.000012.601.29 × 10–4########
CpHTF specific heatkJ/kg °C4.200.0000900.00000.00000.270################
ToutHTF outlet temperature°C46.10.4189900.05800.04420.3784.81 × 10–42.79 × 10–4
TinHTF inlet temperature°C43.10.2280900.05800.02400.3784.81 × 10–48.25 × 10–5
Table 5.

Total solar field energy output uncertainty.

SymbolDescriptionUnitsCalculated mean value
R
Absolute systematic standard uncertainty of result
bR
Absolute random standard uncertainty of result
SR
Combined standard uncertainty
uR
Expanded uncertainty of the result
UR,95%
Expanded uncertainty of the result
UR,95% (%)
QgainSolar field energy outputkWthh1.1460.0000.0000.0000.0000.00%
SymbolDescriptionUnitsCalculated mean value
R
Absolute systematic standard uncertainty of result
bR
Absolute random standard uncertainty of result
SR
Combined standard uncertainty
uR
Expanded uncertainty of the result
UR,95%
Expanded uncertainty of the result
UR,95% (%)
QgainSolar field energy outputkWthh1.1460.0000.0000.0000.0000.00%
Table 5.

Total solar field energy output uncertainty.

SymbolDescriptionUnitsCalculated mean value
R
Absolute systematic standard uncertainty of result
bR
Absolute random standard uncertainty of result
SR
Combined standard uncertainty
uR
Expanded uncertainty of the result
UR,95%
Expanded uncertainty of the result
UR,95% (%)
QgainSolar field energy outputkWthh1.1460.0000.0000.0000.0000.00%
SymbolDescriptionUnitsCalculated mean value
R
Absolute systematic standard uncertainty of result
bR
Absolute random standard uncertainty of result
SR
Combined standard uncertainty
uR
Expanded uncertainty of the result
UR,95%
Expanded uncertainty of the result
UR,95% (%)
QgainSolar field energy outputkWthh1.1460.0000.0000.0000.0000.00%
Table 6.

Solar field thermal efficiency uncertainty for measured data.

Parameter information (in parameter units)Uncertainty contribution of parameters to the result (in units squared)
SymbolDescriptionUnitsNominal valueStandard deviation
σi
Number of measurements
Ni
Absolute systematic standard uncertainty
bX¯i
Absolute random standard uncertainty
SX¯i
Absolute sensitivity
RX¯i
Absolute systematic standard uncertainty contribution
(RX¯ibX¯i)2
Absolute random standard uncertainty contribution
(RX¯iSX¯i)2
m˙Mass flow ratekg/s0.090.0000900.000 900.00000.007 743 414.84 × 10–110.000 000 0
CpHTF specific heatkJ/kg °C4.200.0000900.00000.00000.000 165 930.000 000 00.000 000 0
ToutHTF outlet temperature°C46.10.4189900.05800.04420.000 232 301.81 × 10–101.05 × 10–10
TinHTF inlet temperature°C43.10.2280900.05800.02400.000 232 301.81 × 10–103.10 × 10–11
GbDirect normal irradianceW/m28994.81629018.0200.50770.000 696 911.58 × 10–41.25 × 10–7
AaEffective aperture aream21.810.0000900.00000.00000.000 696 910.000 000 00.000 000 0
Parameter information (in parameter units)Uncertainty contribution of parameters to the result (in units squared)
SymbolDescriptionUnitsNominal valueStandard deviation
σi
Number of measurements
Ni
Absolute systematic standard uncertainty
bX¯i
Absolute random standard uncertainty
SX¯i
Absolute sensitivity
RX¯i
Absolute systematic standard uncertainty contribution
(RX¯ibX¯i)2
Absolute random standard uncertainty contribution
(RX¯iSX¯i)2
m˙Mass flow ratekg/s0.090.0000900.000 900.00000.007 743 414.84 × 10–110.000 000 0
CpHTF specific heatkJ/kg °C4.200.0000900.00000.00000.000 165 930.000 000 00.000 000 0
ToutHTF outlet temperature°C46.10.4189900.05800.04420.000 232 301.81 × 10–101.05 × 10–10
TinHTF inlet temperature°C43.10.2280900.05800.02400.000 232 301.81 × 10–103.10 × 10–11
GbDirect normal irradianceW/m28994.81629018.0200.50770.000 696 911.58 × 10–41.25 × 10–7
AaEffective aperture aream21.810.0000900.00000.00000.000 696 910.000 000 00.000 000 0
Table 6.

Solar field thermal efficiency uncertainty for measured data.

Parameter information (in parameter units)Uncertainty contribution of parameters to the result (in units squared)
SymbolDescriptionUnitsNominal valueStandard deviation
σi
Number of measurements
Ni
Absolute systematic standard uncertainty
bX¯i
Absolute random standard uncertainty
SX¯i
Absolute sensitivity
RX¯i
Absolute systematic standard uncertainty contribution
(RX¯ibX¯i)2
Absolute random standard uncertainty contribution
(RX¯iSX¯i)2
m˙Mass flow ratekg/s0.090.0000900.000 900.00000.007 743 414.84 × 10–110.000 000 0
CpHTF specific heatkJ/kg °C4.200.0000900.00000.00000.000 165 930.000 000 00.000 000 0
ToutHTF outlet temperature°C46.10.4189900.05800.04420.000 232 301.81 × 10–101.05 × 10–10
TinHTF inlet temperature°C43.10.2280900.05800.02400.000 232 301.81 × 10–103.10 × 10–11
GbDirect normal irradianceW/m28994.81629018.0200.50770.000 696 911.58 × 10–41.25 × 10–7
AaEffective aperture aream21.810.0000900.00000.00000.000 696 910.000 000 00.000 000 0
Parameter information (in parameter units)Uncertainty contribution of parameters to the result (in units squared)
SymbolDescriptionUnitsNominal valueStandard deviation
σi
Number of measurements
Ni
Absolute systematic standard uncertainty
bX¯i
Absolute random standard uncertainty
SX¯i
Absolute sensitivity
RX¯i
Absolute systematic standard uncertainty contribution
(RX¯ibX¯i)2
Absolute random standard uncertainty contribution
(RX¯iSX¯i)2
m˙Mass flow ratekg/s0.090.0000900.000 900.00000.007 743 414.84 × 10–110.000 000 0
CpHTF specific heatkJ/kg °C4.200.0000900.00000.00000.000 165 930.000 000 00.000 000 0
ToutHTF outlet temperature°C46.10.4189900.05800.04420.000 232 301.81 × 10–101.05 × 10–10
TinHTF inlet temperature°C43.10.2280900.05800.02400.000 232 301.81 × 10–103.10 × 10–11
GbDirect normal irradianceW/m28994.81629018.0200.50770.000 696 911.58 × 10–41.25 × 10–7
AaEffective aperture aream21.810.0000900.00000.00000.000 696 910.000 000 00.000 000 0
Table 7.

Total solar field thermal efficiency uncertainty.

SymbolDescriptionUnitsCalculated mean value
R
Absolute systematic standard uncertainty of result
bR
Absolute random standard uncertainty of result
SR
Combined standard uncertainty
uR
Expanded uncertainty of the result
UR,95%
Expanded uncertainty of the result
UR,95% (%)
ηgSolar field efficiency0.7060.01260.00040.01290.02583.66%
SymbolDescriptionUnitsCalculated mean value
R
Absolute systematic standard uncertainty of result
bR
Absolute random standard uncertainty of result
SR
Combined standard uncertainty
uR
Expanded uncertainty of the result
UR,95%
Expanded uncertainty of the result
UR,95% (%)
ηgSolar field efficiency0.7060.01260.00040.01290.02583.66%
Table 7.

Total solar field thermal efficiency uncertainty.

SymbolDescriptionUnitsCalculated mean value
R
Absolute systematic standard uncertainty of result
bR
Absolute random standard uncertainty of result
SR
Combined standard uncertainty
uR
Expanded uncertainty of the result
UR,95%
Expanded uncertainty of the result
UR,95% (%)
ηgSolar field efficiency0.7060.01260.00040.01290.02583.66%
SymbolDescriptionUnitsCalculated mean value
R
Absolute systematic standard uncertainty of result
bR
Absolute random standard uncertainty of result
SR
Combined standard uncertainty
uR
Expanded uncertainty of the result
UR,95%
Expanded uncertainty of the result
UR,95% (%)
ηgSolar field efficiency0.7060.01260.00040.01290.02583.66%

The calculated total uncertainty, expressed as the expanded uncertainty, of the mean observed solar–thermal output during the performance test at the 95% confidence level is 9.05%, which is equivalent to 0.104 kWhth. The total expanded uncertainty of the solar field thermal efficiency is in the order of 3.66% (equivalent to 0.0258) at the same confidence level of 95%. The mean collector efficiency that was determined from experimental tests is 70.6%. The maximum efficiency of the solar collector is ~76%, which is attributed to the good-quality components with the properties listed in Table 1 and the operation of the system at lower temperatures. The uncertainties in the thermal output and efficiency are taken into account when analysing the data, allowing a more accurate assessment of the experimental models and validation of the SAM PTC model, in the following sections.

The experimental results indicate that the main source of uncertainty in calculating the solar field thermal output and efficiency is due to errors in measuring the solar irradiance and temperature difference. These errors can be minimized by employing high-quality sensors with very low uncertainty margins, thereby enhancing measurement reliability. Additionally, a more robust design of the solar PTC prototype can significantly improve the overall measurement accuracy, optimizing the experimental set-up to mitigate potential sources of error such as vibration and misalignment. Furthermore, maintaining a controlled experimental environment can help to reduce external influences on the performance measurements. In this study, factors such as wind and the presence of dust or soiling on the collector may have negatively impacted the solar irradiance measurements and hence the overall energy output and efficiency, particularly given the significant wind conditions at the measurement site.

4.3 Analysis of experimental models

Table 8 is a summary of the results of measurement data from the short-duration performance acceptance test that was used in developing the thermal output and efficiency experimental models by using the methodology and equations that were discussed in Section 3.2. The regression coefficients for the two models are presented in Table 9.

Table 8.

Measurement data for experimental analysis.

Count
no.
Gb (W/m2)Tm
(°C)
ΔT
(°C)
T*
(°C/W m2)
Count
no.
Gb
(W/m2)
Tm
(°C m2/W)
ΔT
(Wthh)
T*
(°C/W m2)
189343.7316.10.016 464689544.4515.50.015 75
289643.7816.60.016 854789744.6015.80.016 05
389443.7015.90.016 224889944.6516.30.016 46
489243.8316.20.016 594990244.7016.80.016 96
589243.8316.10.016 485090044.7016.70.016 89
689543.8016.20.016 545190244.7317.00.017 18
789343.8816.60.016 915290144.7016.70.016 87
889443.8015.90.016 225389744.7317.00.017 28
989443.8015.90.016 225489844.7016.80.017 04
1089443.8516.40.016 675589844.7016.90.017 15
1189143.7815.70.016 055689744.8317.00.017 28
1289443.9316.20.016 555790044.8316.90.017 11
1389544.0516.40.016 655890344.8817.50.017 61
1489044.0817.00.017 425990444.8517.20.017 31
1589144.1317.70.018 186090444.8517.30.017 37
1689344.0816.70.017 026190244.7816.50.016 63
1789744.1016.90.017 176290144.9317.00.017 20
1889644.1317.20.017 526389944.9016.70.016 91
1989644.2317.20.017 526489844.9317.10.017 32
2089744.2317.10.017 396589744.9317.00.017 28
2189844.1816.50.016 706690245.0517.20.017 35
2289644.1316.10.016 416790245.0817.50.017 63
2389044.2017.20.017 646890345.1017.60.017 72
2489244.2317.20.017 606990345.0517.20.017 28
2589744.3317.00.017 287090445.1017.60.017 70
2689744.3517.40.017 617190545.2017.60.017 68
2789844.3317.00.017 267290545.2317.90.018 01
2889844.4817.50.017 717390445.2317.90.018 03
2989944.4016.80.017 027490345.2017.60.017 72
3089644.3016.00.016 297590445.2017.70.017 81
3189844.4317.10.017 377690545.2317.80.017 90
3289944.4317.00.017 247790745.2017.70.017 75
3390044.4016.80.017 007890745.2017.60.017 64
3489944.4016.70.016 917990645.1817.40.017 44
3589544.4517.30.017 548090645.3317.90.017 99
3689644.4816.60.016 858190745.3017.50.017 53
3789944.4516.30.016 468290445.3017.50.017 59
3890044.5016.70.016 898390645.3317.90.017 99
3989744.5317.00.017 288490845.3017.60.017 57
4089644.5817.60.017 868590845.3317.80.017 84
4189444.5317.00.017 348690745.3017.60.017 64
4290044.6817.50.017 678790845.3317.70.017 68
4390044.6016.80.017 008890745.3518.10.018 08
4490144.6016.70.016 878990845.4317.90.017 95
4589944.5516.40.016 579090845.4317.70.017 73
Count
no.
Gb (W/m2)Tm
(°C)
ΔT
(°C)
T*
(°C/W m2)
Count
no.
Gb
(W/m2)
Tm
(°C m2/W)
ΔT
(Wthh)
T*
(°C/W m2)
189343.7316.10.016 464689544.4515.50.015 75
289643.7816.60.016 854789744.6015.80.016 05
389443.7015.90.016 224889944.6516.30.016 46
489243.8316.20.016 594990244.7016.80.016 96
589243.8316.10.016 485090044.7016.70.016 89
689543.8016.20.016 545190244.7317.00.017 18
789343.8816.60.016 915290144.7016.70.016 87
889443.8015.90.016 225389744.7317.00.017 28
989443.8015.90.016 225489844.7016.80.017 04
1089443.8516.40.016 675589844.7016.90.017 15
1189143.7815.70.016 055689744.8317.00.017 28
1289443.9316.20.016 555790044.8316.90.017 11
1389544.0516.40.016 655890344.8817.50.017 61
1489044.0817.00.017 425990444.8517.20.017 31
1589144.1317.70.018 186090444.8517.30.017 37
1689344.0816.70.017 026190244.7816.50.016 63
1789744.1016.90.017 176290144.9317.00.017 20
1889644.1317.20.017 526389944.9016.70.016 91
1989644.2317.20.017 526489844.9317.10.017 32
2089744.2317.10.017 396589744.9317.00.017 28
2189844.1816.50.016 706690245.0517.20.017 35
2289644.1316.10.016 416790245.0817.50.017 63
2389044.2017.20.017 646890345.1017.60.017 72
2489244.2317.20.017 606990345.0517.20.017 28
2589744.3317.00.017 287090445.1017.60.017 70
2689744.3517.40.017 617190545.2017.60.017 68
2789844.3317.00.017 267290545.2317.90.018 01
2889844.4817.50.017 717390445.2317.90.018 03
2989944.4016.80.017 027490345.2017.60.017 72
3089644.3016.00.016 297590445.2017.70.017 81
3189844.4317.10.017 377690545.2317.80.017 90
3289944.4317.00.017 247790745.2017.70.017 75
3390044.4016.80.017 007890745.2017.60.017 64
3489944.4016.70.016 917990645.1817.40.017 44
3589544.4517.30.017 548090645.3317.90.017 99
3689644.4816.60.016 858190745.3017.50.017 53
3789944.4516.30.016 468290445.3017.50.017 59
3890044.5016.70.016 898390645.3317.90.017 99
3989744.5317.00.017 288490845.3017.60.017 57
4089644.5817.60.017 868590845.3317.80.017 84
4189444.5317.00.017 348690745.3017.60.017 64
4290044.6817.50.017 678790845.3317.70.017 68
4390044.6016.80.017 008890745.3518.10.018 08
4490144.6016.70.016 878990845.4317.90.017 95
4589944.5516.40.016 579090845.4317.70.017 73
Table 8.

Measurement data for experimental analysis.

Count
no.
Gb (W/m2)Tm
(°C)
ΔT
(°C)
T*
(°C/W m2)
Count
no.
Gb
(W/m2)
Tm
(°C m2/W)
ΔT
(Wthh)
T*
(°C/W m2)
189343.7316.10.016 464689544.4515.50.015 75
289643.7816.60.016 854789744.6015.80.016 05
389443.7015.90.016 224889944.6516.30.016 46
489243.8316.20.016 594990244.7016.80.016 96
589243.8316.10.016 485090044.7016.70.016 89
689543.8016.20.016 545190244.7317.00.017 18
789343.8816.60.016 915290144.7016.70.016 87
889443.8015.90.016 225389744.7317.00.017 28
989443.8015.90.016 225489844.7016.80.017 04
1089443.8516.40.016 675589844.7016.90.017 15
1189143.7815.70.016 055689744.8317.00.017 28
1289443.9316.20.016 555790044.8316.90.017 11
1389544.0516.40.016 655890344.8817.50.017 61
1489044.0817.00.017 425990444.8517.20.017 31
1589144.1317.70.018 186090444.8517.30.017 37
1689344.0816.70.017 026190244.7816.50.016 63
1789744.1016.90.017 176290144.9317.00.017 20
1889644.1317.20.017 526389944.9016.70.016 91
1989644.2317.20.017 526489844.9317.10.017 32
2089744.2317.10.017 396589744.9317.00.017 28
2189844.1816.50.016 706690245.0517.20.017 35
2289644.1316.10.016 416790245.0817.50.017 63
2389044.2017.20.017 646890345.1017.60.017 72
2489244.2317.20.017 606990345.0517.20.017 28
2589744.3317.00.017 287090445.1017.60.017 70
2689744.3517.40.017 617190545.2017.60.017 68
2789844.3317.00.017 267290545.2317.90.018 01
2889844.4817.50.017 717390445.2317.90.018 03
2989944.4016.80.017 027490345.2017.60.017 72
3089644.3016.00.016 297590445.2017.70.017 81
3189844.4317.10.017 377690545.2317.80.017 90
3289944.4317.00.017 247790745.2017.70.017 75
3390044.4016.80.017 007890745.2017.60.017 64
3489944.4016.70.016 917990645.1817.40.017 44
3589544.4517.30.017 548090645.3317.90.017 99
3689644.4816.60.016 858190745.3017.50.017 53
3789944.4516.30.016 468290445.3017.50.017 59
3890044.5016.70.016 898390645.3317.90.017 99
3989744.5317.00.017 288490845.3017.60.017 57
4089644.5817.60.017 868590845.3317.80.017 84
4189444.5317.00.017 348690745.3017.60.017 64
4290044.6817.50.017 678790845.3317.70.017 68
4390044.6016.80.017 008890745.3518.10.018 08
4490144.6016.70.016 878990845.4317.90.017 95
4589944.5516.40.016 579090845.4317.70.017 73
Count
no.
Gb (W/m2)Tm
(°C)
ΔT
(°C)
T*
(°C/W m2)
Count
no.
Gb
(W/m2)
Tm
(°C m2/W)
ΔT
(Wthh)
T*
(°C/W m2)
189343.7316.10.016 464689544.4515.50.015 75
289643.7816.60.016 854789744.6015.80.016 05
389443.7015.90.016 224889944.6516.30.016 46
489243.8316.20.016 594990244.7016.80.016 96
589243.8316.10.016 485090044.7016.70.016 89
689543.8016.20.016 545190244.7317.00.017 18
789343.8816.60.016 915290144.7016.70.016 87
889443.8015.90.016 225389744.7317.00.017 28
989443.8015.90.016 225489844.7016.80.017 04
1089443.8516.40.016 675589844.7016.90.017 15
1189143.7815.70.016 055689744.8317.00.017 28
1289443.9316.20.016 555790044.8316.90.017 11
1389544.0516.40.016 655890344.8817.50.017 61
1489044.0817.00.017 425990444.8517.20.017 31
1589144.1317.70.018 186090444.8517.30.017 37
1689344.0816.70.017 026190244.7816.50.016 63
1789744.1016.90.017 176290144.9317.00.017 20
1889644.1317.20.017 526389944.9016.70.016 91
1989644.2317.20.017 526489844.9317.10.017 32
2089744.2317.10.017 396589744.9317.00.017 28
2189844.1816.50.016 706690245.0517.20.017 35
2289644.1316.10.016 416790245.0817.50.017 63
2389044.2017.20.017 646890345.1017.60.017 72
2489244.2317.20.017 606990345.0517.20.017 28
2589744.3317.00.017 287090445.1017.60.017 70
2689744.3517.40.017 617190545.2017.60.017 68
2789844.3317.00.017 267290545.2317.90.018 01
2889844.4817.50.017 717390445.2317.90.018 03
2989944.4016.80.017 027490345.2017.60.017 72
3089644.3016.00.016 297590445.2017.70.017 81
3189844.4317.10.017 377690545.2317.80.017 90
3289944.4317.00.017 247790745.2017.70.017 75
3390044.4016.80.017 007890745.2017.60.017 64
3489944.4016.70.016 917990645.1817.40.017 44
3589544.4517.30.017 548090645.3317.90.017 99
3689644.4816.60.016 858190745.3017.50.017 53
3789944.4516.30.016 468290445.3017.50.017 59
3890044.5016.70.016 898390645.3317.90.017 99
3989744.5317.00.017 288490845.3017.60.017 57
4089644.5817.60.017 868590845.3317.80.017 84
4189444.5317.00.017 348690745.3017.60.017 64
4290044.6817.50.017 678790845.3317.70.017 68
4390044.6016.80.017 008890745.3518.10.018 08
4490144.6016.70.016 878990845.4317.90.017 95
4589944.5516.40.016 579090845.4317.70.017 73
Table 9.

Regression coefficients.

Thermal energy outputThermal efficiency
a1a2a3U0U1
(W/m2/°C)
U2
(W/m2/°C)
Parameter value–6.760.00770.056–1.53131.80–38.34
Thermal energy outputThermal efficiency
a1a2a3U0U1
(W/m2/°C)
U2
(W/m2/°C)
Parameter value–6.760.00770.056–1.53131.80–38.34
Table 9.

Regression coefficients.

Thermal energy outputThermal efficiency
a1a2a3U0U1
(W/m2/°C)
U2
(W/m2/°C)
Parameter value–6.760.00770.056–1.53131.80–38.34
Thermal energy outputThermal efficiency
a1a2a3U0U1
(W/m2/°C)
U2
(W/m2/°C)
Parameter value–6.760.00770.056–1.53131.80–38.34

4.4 Thermal output experimental model

The thermal output experimental model equation is presented below:

(24)

To determine whether the system has met the acceptance criteria, we need to compare the actual thermal output of the system (including its associated uncertainty) with the predicted thermal output that was calculated by using Equation (24). In statistical terms, we are testing the “null hypothesis” that the observed thermal output does not exceed the predicted thermal output. The purpose of the acceptance test is to verify that the thermal output of the system meets or exceeds a predetermined level so, if our test shows that the null hypothesis is false, then it indicates that the system has successfully passed the acceptance test [27, 39].

From Table 5, the mean thermal output observed during the test is 1.150 kWhth, with an uncertainty of 9.05% (or 0.104 kWhth) at the 95% confidence level. In order to prove the null hypothesis, the predicted thermal output value must be greater than the observed value plus its uncertainty, which is 1.150 + 0.104 = 1.254 kWhth. This means that, if the modeled thermal output is >1.254 kWhth, then we can be 95% confident that the solar field is not performing at the guaranteed level that is predicted by the experimental model. Fig. 3 shows the observed and modeled thermal output values over time, as well as the average values (dotted lines) for the duration of the acceptance test.

Comparison of observed and predicted thermal output.
Figure 3.

Comparison of observed and predicted thermal output.

The mean experimental model thermal output value for the test period is 1.117 kWhth, which is significantly lower than the 1.254 kWhth required to prove the null hypothesis. As a result, the null hypothesis is disproved, as the predicted thermal output does not exceed the observed level. Therefore, the system has successfully passed the acceptance test based on its thermal output performance. The coefficient of determination, or regression, R2 is used to assess the ability of the experimental model to accurately predict the thermal output. From calculations, the R2 value is 0.972, which is statistically significant and acceptable. This suggests that the experimental model is a reliable tool for predicting the thermal output of the collector.

4.5 Thermal efficiency experimental model

Based on the regression coefficient parameters that are presented in Table 9, the thermal efficiency curve of the experimental model is presented below:

(25)

The predicted thermal efficiency of the model is 72.5%, which is lower than the observed value of 73.2% (0.706 + 0.0258) obtained from the data (see Table 7). This also disproves the null hypothesis, as it suggests that the predicted thermal efficiency does not exceed the observed level. As a result, the system has passed the acceptance test based on its thermal efficiency performance.

Additionally, the R2 value of 0.989 indicates that the model provides a good fit to the observed data, suggesting that it is reliable for predicting thermal efficiency by using measured data. The plot in Fig. 4 shows the comparison between the observed and predicted thermal efficiency values, along with the average values (dotted lines), for the duration of the acceptance test.

Comparison of observed and predicted thermal efficiency.
Figure 4.

Comparison of observed and predicted thermal efficiency.

4.6 Validation of the SAM PTC model

Following the verification of the ability of the experimental models to accurately predict thermal output and efficiency of the collector, the models are used to validate the performance of the parabolic trough collector as predicted by using the SAM software. The SAM model is tested by simulating various scenarios using the same boundary conditions as the experimental data (collector properties, solar irradiance, mass flow, inlet, and outlet temperatures, ambient temperature, wind speed, and relative humidity). The data for this analysis were collected on a separate day between 9 a.m. and 6 p.m., at an interval of 1 hour.

The simulations were carried out hourly, replicating operational conditions that were similar to those in the experimental set-up. Each simulation time step incorporated specific input conditions that were derived from the experimental set-up, as follows:

  • Meteorological data—the solar irradiance, ambient temperature, wind speed, wind direction, atmospheric pressure, and relative humidity.

  • Design point parameters—mainly inlet and outlet temperatures of the HTF and the HTF flow rate.

  • Collector and receiver parameters as detailed in Table 1.

Since SAM usually requires extensive annual data for robust simulations, the short-term experimental data are integrated within a broader annual dataset that includes the selected experimental day along with other relevant weather and performance data from the entire year at the location of the experiment. These annual data were sourced from a nearby weather station. However, for validation purposes, the specific SAM simulation results for the days of the experiment were required.

The experimental results are then compared to the simulated outcomes in order to validate the accuracy of the SAM model. The error bars of the experimental results were utilized to validate the SAM model by assessing how closely the simulation results aligned with these error margins. Additional validation was conducted by using three statistical parameters: the root mean square error (RMSE), the mean bias error (MBE), and the t-statistics method (tsta). These performance indices are calculated by using the following equations [54, 55]:

(26)
(27)
(28)

where Imod.,i is the value of the simulation model result at time step i, Imeas,i is the value of the experimental result at time step i and n is the number of values under investigation.

The comparison between the experimental and simulated data, as shown in Figs 5 and 6, reveals that the predictions of the SAM model fall within the error bars of the experimental results. This implies that the SAM model accurately predicts the thermal output and efficiency of the PTC throughout the day, as reflected in Table 10.

Table 10.

Experimental and simulation values for validation of SAM.

Thermal energy outputThermal efficiency
Ta
(°C)
Gb
(W/m2)
Qgain expt
(kWthh)
Qgain SAM
(kWthh)
ηg exptηg SAM
22.6328(0.331 ± 0.104)0.301(0.39 ± 0.026)0.37
23.8433(0.499 ± 0.104)0.438(0.51 ± 0.026)0.5
25.9564(0.646 ± 0.104)0.591(0.53 ± 0.026)0.53
26.3701(0.977 ± 0.104)0.901(0.61 ± 0.026)0.62
27.6746(1.054 ± 0.104)1.001(0.68 ± 0.026)0.69
28.1859(1.108 ± 0.104)1.0970.72 ± 0.026)0.7
27.8788(1.007 ± 0.104)1.093(0.70 ± 0.026)0.71
26.3744(0.933 ± 0.104)0.9640.68 ± 0.026)0.69
23.8484(0.688 ± 0.104)0.632(0.67 ± 0.026)0.69
21.1302(0.494 ± 0.104)0.446(0.54 ± 0.026)0.58
Thermal energy outputThermal efficiency
Ta
(°C)
Gb
(W/m2)
Qgain expt
(kWthh)
Qgain SAM
(kWthh)
ηg exptηg SAM
22.6328(0.331 ± 0.104)0.301(0.39 ± 0.026)0.37
23.8433(0.499 ± 0.104)0.438(0.51 ± 0.026)0.5
25.9564(0.646 ± 0.104)0.591(0.53 ± 0.026)0.53
26.3701(0.977 ± 0.104)0.901(0.61 ± 0.026)0.62
27.6746(1.054 ± 0.104)1.001(0.68 ± 0.026)0.69
28.1859(1.108 ± 0.104)1.0970.72 ± 0.026)0.7
27.8788(1.007 ± 0.104)1.093(0.70 ± 0.026)0.71
26.3744(0.933 ± 0.104)0.9640.68 ± 0.026)0.69
23.8484(0.688 ± 0.104)0.632(0.67 ± 0.026)0.69
21.1302(0.494 ± 0.104)0.446(0.54 ± 0.026)0.58
Table 10.

Experimental and simulation values for validation of SAM.

Thermal energy outputThermal efficiency
Ta
(°C)
Gb
(W/m2)
Qgain expt
(kWthh)
Qgain SAM
(kWthh)
ηg exptηg SAM
22.6328(0.331 ± 0.104)0.301(0.39 ± 0.026)0.37
23.8433(0.499 ± 0.104)0.438(0.51 ± 0.026)0.5
25.9564(0.646 ± 0.104)0.591(0.53 ± 0.026)0.53
26.3701(0.977 ± 0.104)0.901(0.61 ± 0.026)0.62
27.6746(1.054 ± 0.104)1.001(0.68 ± 0.026)0.69
28.1859(1.108 ± 0.104)1.0970.72 ± 0.026)0.7
27.8788(1.007 ± 0.104)1.093(0.70 ± 0.026)0.71
26.3744(0.933 ± 0.104)0.9640.68 ± 0.026)0.69
23.8484(0.688 ± 0.104)0.632(0.67 ± 0.026)0.69
21.1302(0.494 ± 0.104)0.446(0.54 ± 0.026)0.58
Thermal energy outputThermal efficiency
Ta
(°C)
Gb
(W/m2)
Qgain expt
(kWthh)
Qgain SAM
(kWthh)
ηg exptηg SAM
22.6328(0.331 ± 0.104)0.301(0.39 ± 0.026)0.37
23.8433(0.499 ± 0.104)0.438(0.51 ± 0.026)0.5
25.9564(0.646 ± 0.104)0.591(0.53 ± 0.026)0.53
26.3701(0.977 ± 0.104)0.901(0.61 ± 0.026)0.62
27.6746(1.054 ± 0.104)1.001(0.68 ± 0.026)0.69
28.1859(1.108 ± 0.104)1.0970.72 ± 0.026)0.7
27.8788(1.007 ± 0.104)1.093(0.70 ± 0.026)0.71
26.3744(0.933 ± 0.104)0.9640.68 ± 0.026)0.69
23.8484(0.688 ± 0.104)0.632(0.67 ± 0.026)0.69
21.1302(0.494 ± 0.104)0.446(0.54 ± 0.026)0.58
Comparison of predicted and simulated thermal output results.
Figure 5.

Comparison of predicted and simulated thermal output results.

Comparison of predicted and simulated efficiency results.
Figure 6.

Comparison of predicted and simulated efficiency results.

Moreover, the analysis that employs the three statistical metrics that were derived from the experimental and simulation values presented in Table 10 demonstrates statistically significant results. For the thermal energy output, the calculated values for the RMSE, MBE, and t-statistics were 0.057, 0.033, and 2.135, respectively. Similarly, for the thermal efficiency, the corresponding values for the RMSE, MBE, and t-statistics were 0.018, 0.010, and 1.993, respectively. As a result, the SAM PTC model has been successfully validated, demonstrating its reliability in simulating the performance of the collector.

It is essential to recognize that the experimental testing of the solar PTC has two significant limitations, and those are soiling and dust deposition, which can hinder the performance measurement of the parabolic trough by blocking the solar radiation from reaching the surface of the collector. This accumulation reduces the effective aperture area, leading to lower solar irradiance and thermal energy absorption, and thus diminishing the efficiency. Moreover, it complicates the interpretation of the performance data, as the output variations may be wrongly attributed to the system design rather than environmental conditions. Regular cleaning and maintenance are necessary for accurate assessments, but they can also introduce additional variability in the results. Despite this challenge, experimental testing is deemed essential for validating simulation models and ensuring trust in their applications.

5. Conclusion

This study outlines a detailed acceptance test process that integrates uncertainty analysis to validate experimental data in predicting the performance of a PTC prototype. The findings are used to create experimental performance models, which help to validate a simulation model that is developed in SAM for industrial heating applications. The acceptance test focuses on analysing the thermal energy output and efficiency of the solar field under clear-sky conditions during a steady state, reducing the risk of modeling errors. The results from the experimental models and their uncertainties are compared with the performance projections from the SAM model for these parameters.

The uncertainty analysis provided a realistic evaluation of the thermal performance of the prototype during testing, resulting in expanded uncertainties of 9.05% (0.104 kWhth) for thermal output and 3.66% (0.0258) for thermal efficiency. The predictions of the experimental models showed strong correlations with the observed values, indicated by R2 regression coefficients of 0.972 for thermal output and 0.989 for thermal efficiency. Comparison of the experimental results with the simulation outcomes confirmed the validity of the simulation model, as the results fell within the experimental error margins. Additionally, statistical analysis yielded significant values: RMSE of 0.057, MBE of 0.033, and t-statistics of 2.135 for thermal energy output, and RMSE of 0.018, MBE of 0.010, and t-statistics of 1.993 for thermal efficiency, underscoring the accuracy and reliability of the simulation model. Consequently, the simulation model has been validated and can effectively predict the thermal output and efficiency of parabolic trough collectors, and it can be utilized to forecast their thermal behavior under various operating conditions in Kenya, even without experimental data.

The four outlined simulation model validation steps—uncertainty analysis, regression model development, quality evaluation of experimental models, and simulation model validation—have produced several key research outputs that can enhance the design and implementation of solar PTC systems in industrial applications. First, the uncertainty analysis tables provide a reference for future studies, helping to refine experimental set-ups by identifying sources of variability. Second, the regression models offer a robust mathematical framework for predictive analytics, linking input variables such as solar irradiance and HTF temperatures to performance metrics such as thermal efficiency. Third, the quality-of-fit assessment using metrics such as R-squared values enhances the credibility of the models, allowing researchers to ascertain their representation of experimental data. Finally, validated performance predictions ensure that simulations accurately reflect real-world performance. Together, these outputs generate critical insights and can facilitate the development of more efficient, reliable, and scalable solar–thermal solutions for industrial applications.

Author contributions

Philip Akello (Conceptualization [Lead], Formal analysis [Supporting], Resources [Equal], Software [Equal], Writing—original draft [Lead]), Churchill Saoke (Formal analysis [Equal], Software [Supporting], Visualization [Lead], Writing—review & editing [Supporting]), Joseph Kamau (Formal analysis [Supporting], Methodology [Equal], Validation [Equal], Writing—review & editing [Equal]), and Jared Ndeda (Formal analysis [Supporting], Methodology [Equal], Validation [Equal], Writing—review & editing [Equal])

Conflict of interest statement

None declared.

Funding

None declared.

Data availability

The main data underlying this article are available in the article. Any remaining data will be shared on reasonable request to the corresponding author.

References

[1]

Akello
POO
,
Saoke
CO
,
Kamau
JN
et al.
Modeling and performance analysis of solar parabolic trough collectors for hybrid process heat application in Kenya’s tea industry using system advisor model
.
Sustain Energy Res
2023
;
10
:
7
. https://doi.org/

[2]

Mouaky
A
,
Merrouni
AA
,
Laadel
NE.
Simulation and experimental validation of a parabolic trough plant for solar thermal applications under the semi-arid climate conditions
.
Sol Energy
2019
;
194
:
969
85
. https://doi.org/

[3]

IEA
,
International Energy Agency
.
2019
. https://iea.blob.core.windows.net/assets/a846e5cf-ca7d-4a1f-a81b-ba1499f2cc07/Renewables_2019.pdf 
(30 September 2023, date last accessed)
.

[4]

Akbarzadeh
S
,
Valipour
MS.
Heat transfer enhancement in parabolic trough collectors: a comprehensive review
.
Renew Sustain Energy Rev
2018
;
92
:
198
218
. https://doi.org/

[5]

Conrado
LS
,
Rodriguez-Pulido
A
,
Calderón
G.
Thermal performance of parabolic trough solar collectors
.
Renew Sustain Energy Rev
2017
;
67
:
1345
59
. https://doi.org/

[6]

Wu
T
,
Zhang
H
,
Shang
L.
Optimal sizing of a grid-connected hybrid renewable energy systems considering hydroelectric storage
.
Energy Sources Part A
2020
;
46
:
5043
59
. https://doi.org/.

[7]

IEA, International Energy Agency
.
World Energy Statistics 2016, Online Tables
. www.iea.org/statistics/ 
(6 December 2018, date last accessed)
.

[8]

Häberle
A.
Concentrating solar technologies for industrial process heat and cooling
. In:
Lovegrove
K
,
Stein
W
(eds),
Concentrating Solar Power Technology
, 2nd edn.
Cambridge, UK
:
Woodhead Publishing
,
2021
,
602
19
. https://doi.org/

[9]

Hussain
CI
,
Norton
B
,
Duffy
A.
Technological assessment of different solar-biomass systems for hybrid power generation in Europe
.
Renew Sustain Energy Rev
2017
;
68
:
1115
29
. https://doi.org/

[10]

World Bank
.
Concentrating solar power: Clean Solar Power 24/7
.
Washington, DC
,
2021
. https://pubdocs.worldbank.org/en/849341611761898393/WorldBank-CSP-Report-Concentrating-Solar-Power-Clean-Power-on-Demand-24-7-FINAL.pdf
(21 June 2023, date last accessed)
.

[11]

Quezada–García
S
,
Sánchez–Mora
H
,
Polo–Labarrios
MA
et al.
Modeling and simulation to determine the thermal efficiency of a parabolic solar trough collector system
.
Case Stud Thermal Eng
2019
;
16
:
100523
. https://doi.org/

[12]

Ouali
HAL
,
Merrouni
AA
,
Moussaoui
MA
et al.
Electricity yield analysis of a 50 MW solar tower plant under Moroccan climate
. In: 2015 International Conference on Electrical and Information Technologies (ICEIT), Marrakech, Morocco,
25–27 March
.
New York
:
IEEE
,
2015
,
252
6
. https://doi.org/

[13]

Liqreina
A
,
Qoaider
L.
Dry cooling of concentrating solar power (CSP) plants, an economic competitive option for the desert regions of the MENA region
.
Sol Energy
2014
;
103
:
417
24
. https://doi.org/

[14]

Remlaoui
A
,
Benyoucef
M
,
Assi
D
et al.
A TRNSYS dynamic simulation model for a parabolic trough solar thermal power plant
.
Int J Energetica
2020
;
4
:
36
41
. https://doi.org/

[15]

Jones
SA
,
Pitz-Paal
R
,
Schwarzboezl
P
et al.
TRNSYS modeling of the SEGS VI parabolic trough solar electric generating system
. In: International Solar Energy Conference, Washington, DC, USA,
21–25 April 2001
, Vol.
16702
.
New York
:
American Society of Mechanical Engineers
,
2001
,
405
12
. https://doi.org/

[16]

Blair
N
,
DiOrio
N
,
Freeman
J
et al.
System Advisor Model (SAM) General Description (Version 2017.9. 5) (No. NREL/TP-6A20-70414)
.
Golden, CO
:
National Renewable Energy Laboratory (NREL)
,
2018
. https://doi.org/

[17]

Naaim
S
,
Ouhammou
B
,
Aggour
M
et al.
Multi-utility solar thermal systems: harnessing parabolic trough concentrator using SAM software for diverse industrial and residential applications
.
Energies
2024
;
17
:
3685
. https://doi.org/

[18]

de Simón-Martín
M
,
Díez-Mediavilla
M
,
Alonso-Tristán
C.
Modelling solar data: reasons, main methods and applications
. In: International Conference on Renewable Energies and Power Quality (ICREPQ’13) Proceedings, Bilbao, Spain,
20–22 March 2013
, Vol.
1
, pp.
1
5
,
2013
. https://doi.org/

[19]

Barhmi
K
,
Heynen
C
,
Golroodbari
S
van Sark
W
A review of solar forecasting techniques and the role of artificial intelligence
.
Solar
2024
;
4
:
99
135
.  https://doi.org/

[20]

Gueymard
CA.
A review of validation methodologies and statistical performance indicators for modeled solar radiation data: towards a better bankability of solar projects
.
Renew Sustain Energy Rev
2014
;
39
:
1024
34
. https://doi.org/

[21]

Jagdev
HS
,
Browne
J
,
Jordan
P.
Verification and validation issues in manufacturing models
.
Comput Ind
1995
;
25
:
331
53
. https://doi.org/

[22]

Collins
A
,
Koehler
M
,
Lynch
C.
Methods that support the validation of agent-based models: an overview and discussion
.
J Artif Soc Soc Simul
2024
;
27
:
Article 11
. https://doi.org/

[23]

Fernández-García
A
,
Zarza
E
,
Valenzuela
L
et al.
Parabolic-trough solar collectors and their applications
.
Renew Sustain Energy Rev
2010
;
14
:
1695
721
. https://doi.org/

[24]

Schiricke
B
,
Pitz-Paal
R
,
Lüpfert
E
et al.
(2009). Experimental verification of optical modeling of parabolic trough collectors by flux measurement
.
J Sol Energy Eng
0110
;
131
:
04
. https://doi.org/

[25]

Alfellag
MAA.
Modeling and experimental investigation of parabolic trough solar collector
.
Master’s thesis
,
2014
Florida
Embry-Riddle Aeronautical University
. https://commons.erau.edu/edt/11

[26]

Mateo
M
,
Echarri
R
,
Samsón
I.
Thermal analysis and experimental validation of parabolic trough collector for solar adsorption refrigerator
.
Energy Power Eng
2017
;
09
:
687
702
. https://doi.org/

[27]

Soudani
ME
,
Aiadi
KE
,
Bechki
D
et al.
Experimental and theoretical study of parabolic trough collector (PTC) with a flat glass cover in the region of Algerian Sahara (Ouargla)
.
J Mech Sci Technol
2017
;
31
:
4003
9
. https://doi.org/

[28]

Boretti
A
,
Nayfeh
J
,
Al-Kouz
W.
Validation of SAM modeling of concentrated solar power plants
.
Energies
2020
;
13
:
1949
. https://doi.org/

[29]

Wagner
MJ
,
Mehos
MS
,
Kearney
DW
et al.
Modeling of a parabolic trough solar field for acceptance testing: a case study
. In: ASME 2011 5th International Conference on Energy Sustainability, Washington, DC, USA,
7–10 August 2011
, 54686.
2011
,
595
603
. https://doi.org/

[30]

Mathioulakis
E
,
Voropoulos
K
,
Belessiotis
V.
Assessment of uncertainty in solar collector modeling and testing
.
Sol Energy
1999
;
66
:
337
47
. https://doi.org/

[31]

Kalogirou
SA.
Solar Energy Engineering: Processes and Systems
.
Amsterdam
:
Elsevier
,
2023
.

[32]

Akhbari
H
,
Bidi
M
,
Bakhtiari
A
et al.
Design and fabrication of a parabolic trough collector and experimental investigation of wind impact on direct steam production in Tehran
.
Int J Mech Ind Eng
2019
;
13
:
140
6
.

[33]

Guenther
M
,
Michael
J
,
Csambor
S
et al.
Parabolic trough technology
. In:
Advanced CSP Teaching Materials
.
Cologne
:
Deutsches Zentrums für Luft- und Raumfahrt
,
2011
,
1
106
. https://www.researchgate.net/publication/289941948_Parabolic_Trough_Technology

[34]

Jaramillo
OA
,
Venegas-Reyes
E
,
Aguilar
JO
et al.
Parabolic trough concentrators for low enthalpy processes
.
Renew Energy
2013
;
60
:
529
39
. https://doi.org/

[35]

El-Husseini
SM
,
Mostafa
M
,
El-Gindy
AM
et al.
Solar heating system using parabolic collector for thermal optimum conditions of biogas production in winter
.
Arab Univ J Agric Sci
2019
;
27
:
123
34
. https://doi.org/

[36]

Prasad
PI
,
Basha
SG
,
Janjanam
N.
Design, fabrication and evaluation of thermal performance of parabolic trough collector with elliptical absorber using Zn/H2O nanofluid
.
Int J Adv Sci Technol
2020
;
29
:
7805
18
. https://www.researchgate.net/publication/343988690

[37]

Shao
G
,
Wu
X
,
Cui
S
et al.
High emissivity MoSi2–TaSi2–borosilicate glass porous coating for fibrous ZrO2 ceramic by a rapid sintering method
.
J Alloys Compd
2016
;
690
:
63
71
. https://doi.org/

[38]

Nam
HY
,
Lee
KY
,
Kim
JM
et al.
Experimental study on the emissivity of stainless steel
. In: Proceedings of the Korean Nuclear Society Spring Meeting, Cheju, Republic of Korea,
24–25 May 2001
, Vol.
33
, No.
37
.
2001
.https://inis.iaea.org/search/search.aspx?orig_q=RN:33040336

[39]

Kearney
D.
Utility-Scale Parabolic Trough Solar Systems: Performance Acceptance Test Guidelines, April 2009-December 2010 (No. NREL/SR-5500-48895)
.
Golden, CO
:
National Renewable Energy Laboratory (NREL)
,
2011
. https://doi.org/

[40]

Pashiardis
S
,
Pelengaris
A
,
Kalogirou
SA.
Geographical distribution of global radiation and sunshine duration over the Island of Cyprus
.
Appl Sci
2023
;
13
:
5422
. https://doi.org/

[41]

Elmaihy
A
,
El Weteedy
A.
Short-duration performance acceptance test of parabolic trough solar field: a case study
.
J Sol Energy Eng
20170510
;
139
:
10
. https://doi.org/

[42]

Strutz
T.
Data Fitting and Uncertainty: A Practical Introduction to Weighted Least Squares and Beyond
, Vol.
1
.
Wiesbaden, Germany
:
Vieweg+ Teubner
,
2011
. https://link.springer.com/book/9783658114558

[43]

Pandey
AK.
Least Squares Curve Fitting: Original Data, Fitted Curve, and Residuals
.
2023
https://arunp77.medium.com/least-squares-curve-fitting-original-data-fitted-curve-and-residuals-5349f1b1c892 
(30 March 2024, date last accessed)
.

[44]

Shine
AD.
Fitting Experimental Data to Straight Lines (Including Error Analysis)
.
University of Delaware
.
2006
https://cbe.udel.edu/wp-content/uploads/2019/03/FittingData.pdf
(16 September 2023, date last accessed)
.

[45]

Mathioulakis
E
,
Panaras
G
,
Belessiotis
V.
Uncertainty in estimating the performance of solar thermal systems
.
Sol Energy
2012
;
86
:
3450
9
. https://doi.org/

[46]

Belessiotis
V
,
Mathioulakis
E
,
Papanicolaou
E.
Theoretical formulation and experimental validation of the input–output modeling approach for large solar thermal systems
.
Sol Energy
2010
;
84
:
245
55
. https://doi.org/

[47]

Kratzenberg
MG
,
Beyer
HG
,
Colle
S.
Analytic Determination of the Expanded Uncertainty for Steady State and Quasi Dynamic Collector Tests Under Outdoor Conditions.
2023
. https://lepten.ufsc.br/publicacoes/solar/eventos/2006/EUROSUN/Manfred.pdf. (
28 November 2024
, date last accessed).

[48]

Tasie
NN
,
Sigalo
FB
,
Omubo-Pepple
VB
et al.
Modelling of solar power production in dry and rainy seasons using some selected meteorological parameters
.
Energy Power Eng
2022
;
14
:
274
90
. https://doi.org/

[49]

Margalit
D
,
Rabinoff
J.
Interactive Linear Algebra
.
2019
. https://textbooks.math.gatech.edu/ila/ila.pdf 
(15 April 2024, date last accessed)
.

[50]

Zaversky
F
,
García-Barberena
J
,
Sanchez
M
et al.
Probabilistic modeling of a parabolic trough collector power plant: an uncertainty and sensitivity analysis
.
Sol Energy
2012
;
86
:
2128
39
. https://doi.org/

[51]

Kim
DH
,
Jenkins
BM
,
Rumsey
TR
et al.
Simulation and model validation of a horizontal shallow basin solar concentrator
.
Sol Energy
2007
;
81
:
463
75
. https://doi.org/

[52]

Sabatelli
V
,
Marano
D
,
Braccio
G
et al.
Efficiency test of solar collectors: uncertainty in the estimation of regression parameters and sensitivity analysis
.
Energy Convers Manage
2002
;
43
:
2287
95
. https://doi.org/

[53]

Press
WH
,
Teukolsky
SA
,
Vetterling
WT
et al.
Numerical Recipes
.
Cambridge, UK
:
Cambridge University Press,
1992
https://www.researchgate.net/publication/265107397_Numerical_Recipes
(30 May 2024, date last accessed)
.

[54]

Mecibah
MS
,
Boukelia
TE
,
Benyahia
NE.
Management and exploitation of direct normal irradiance resources for concentrating
.
Int J Energy Environ Eng
2015
;
6
:
65
73
. https://doi.org/

[55]

Kosmopoulos
P
,
Dhake
H
,
Melita
N
et al.
Multi-layer cloud motion vector forecasting for solar energy applications
.
Appl Energy
2024
;
353
:
122144
. https://doi.org/

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.