Skip to main content
  • Original article
  • Open access
  • Published:

Modeling, testing, and parametric analysis of a parabolic solar cooking system with heat storage for indoor cooking



In an ever-changing world where needs increase daily due to economic growth and demographic progression, where prices are unstable, where reserves are running out, where climate change is topical, the energy issues are increasingly marked by the question of sustainability. In many developing countries, wood and subsidized butane are the main sources of energy used for cooking in households. The use of solar energy in domestic cooking becomes unavoidable. Several models of solar cookers have been proposed, but most of them dealt with box and oven types of solar cookers without storage.


This paper presents a dynamic thermodynamic model of a parabolic solar cooking system (PSCS) with heat storage, along with a comparison of the model solution with experimental measurements. The model uses various thermal resistances to take into account heat transfer between the different parts of the system.


The first experimental setup consists of a parabolic concentrator (0.80-m diameter and 0.08-m depth) and a 1.57-l cylindrical receiver. The second experimental setup is composed of a parabolic concentrator (1.40-m diameter and 0.16-m depth), the same receiver, and a 6.64-l heat storage. Tests were carried out in Rabat, Morocco, between April 24 and July 10, 2014, and between May 15 and June 18, 2015. Synthetic oil is used as a transfer fluid and a sensible heat storage.


Comparison between predicted and measured temperatures shows a good agreement with a relative error of ± 4.4%. The effects of important system design and operating parameters were also analyzed. The results show that a 50 W m−2 increase of the daily maximum solar radiation increases the storage temperature by 4 °C and a 5% increase of the receiver reflectance or absorptance improves the maximum storage temperature by 3.6 and 3.9 °C, respectively. Optimizing the aspect ratio of the receiver to 2 gives a maximum storage temperature of 85 °C. Increasing the thermal fluid mass flow rate from 0 to 18 kg h−1, or the receiver thermal insulation from 0.01 to 0.08 m, increases the maximum storage temperature by 65 and 17 °C, respectively.


With increasing population, economic growth, and environmental concerns, the use of solar energy in domestic cooking is becoming a good alternative for sustainable development which will greatly decrease mortality, deforestation, and soil erosion. The World Health Organization (WHO) reports that each year, 1.6 million people die from respiratory diseases caused by indoor air pollution due to solid fuel use for cooking [1]. A domestic solar cooker saves 100 trees in 15 years of life, prevents annually the release of 1.5 ton of CO2, increases the household purchasing power (by reducing the budget allocated to cooking), and gives more time to women and children who spend 15 h per week to the chore of wood [1].

While most solar cookers in use today do not have heat storage, this feature will alleviate the mismatch between solar heat energy supply and energy demand for cooking. Heat storage is important for indoor solar cooking requirements and will ensure continuity of service, reduce the use of conventional energy, and give a reasonable cooking time compared with conventional cooking [2].

Modeling solar concentrating systems including parabolic solar cooking systems (PSCS) is a key tool in order to increase their effectiveness and optimize their operating conditions. Several models of solar cookers have been proposed in the last years, but most of them dealt with box and oven types of solar cookers. Very little modeling work considered detailed dynamic temperature distribution and heat transfer in PSCS with storage. Existing models of parabolic solar systems—other than cooking applications—emphasize optimization of power production and not maximizing fluid temperature. In cooking systems, the fluid temperature determines not only the types of food that can be cooked but also the cooking time.

There is therefore a need to develop a detailed dynamic model of PSCS with heat storage, which will determine the temperature variations in all system components. The present work is focused on developing such a model and on its experimental validation. The “Brief literature review” section presents a brief literature review on the topic of modeling PSCS. The objective is to compare the different modeling approaches, the numerical solutions, the formulas used to assess heat losses, and their validation methods. The “System description and heat transfer processes” section describes the experimental system components and the heat transfer processes involved in its operation. Particular attention was given to the receiver, the key element which absorbs incoming solar radiation, converts it to heat, and transmits it to the heat transfer fluid. The “Governing equations and numerical solution” section gives the governing equations derived from heat balance relationships and heat transfer coefficient formulas and describes their numerical solution. The “Convergence and validation of the numerical solution” section presents the model validation by comparison with other known models and with the experimental results obtained from prototype testing. Finally, the “Parametric analysis” section gives the results of a parametric analysis on the effect of most relevant design elements and operating conditions.

Brief literature review

Solar cooking can be classified into four categories depending on the required temperature range: cooking (85 to 90 °C), boiling (100 to 130 °C), frying (200 to 250 °C), and grilling (over 300 °C) [3]. However, the most frequent classification on solar cooking systems distinguishes between direct and indirect systems [3]. In direct cookers, the solar heat is transferred directly from the reflective surface to the food container (casserole, pot, dish, etc.) whereas in indirect cookers, the pot is physically separate from the collector and a heat-transferring medium is required to convey the heat to the cooking pot. Both types can integrate or not a heat storage system whose size depends on the desired autonomy. Modern technologies often include internal large-scale kitchens for collective applications using fluid transfer with or without heat storage depending on the desired autonomy [3].

Most modeling of solar parabolic systems is concerned with Stirling engine applications, and only few studies have dealt with modeling such systems used in solar cooking. Mawire et al. [4] introduced the energy balance equations to model a solar energy capture (SEC) system and a thermal energy storage (TES) system for an indirect parabolic solar cooker. An oil–pebble bed is used as the TES material. A Simulink block was used to solve the equations and to perform energy and exergy analyses. The results indicate a greater degree of thermal stratification and energy stored when using constant-temperature charging than when using constant-flowrate charging. There are greater initial energy and exergy rates for the constant-flowrate method when the solar radiation is low. Energy efficiencies using both methods are comparable whilst the constant-temperature method results in greater exergy efficiency at higher levels of solar radiation. The same authors [5] used the system to carry out discharging simulations for the thermal energy storage system (TES). It was observed from the results that the TES system at a constant flow rate allows for a higher rate of heat utilization. However, this is not beneficial to the cooking process since the maximum cooking temperature is not maintained for the duration of the discharging period. On the other hand, the controlled load power discharging method (variable flow rate) has a slower initial rate of heat utilization but the maximum cooking temperature is maintained during the whole discharging process, and this is desirable for the cooking process. Prasanna [2] modeled and designed a hybrid solar cooking system consisting of a parabolic collector, a thermal storage tank, and a heat exchanger. The energy source is a combination of solar energy and liquefied petroleum gas (LPG). A bond graph modeling approach was used to build a dynamic model. He found an optimal flow rate allowed a 6% increase in system efficiency as compared to the thermosyphon flow rate. He also found that when the pipe diameter is decreased, the efficiency curve moves up. Mussard [6] developed a low-cost small-scale parabolic trough coupled with a thermal storage unit for higher temperature cooking. The system is built with a self-circulation loop and uses thermal oil. The thermal behavior of the system was simulated using the finite-volume method. He compared different sensible and latent heat storage materials and concluded on the relevance of latent heat-based systems. He also showed that a glass cover with an air gap around the absorber would not improve the efficiency at low temperatures, but when reaching high temperatures (around 220 °C), thermal insulation would be necessary. A storage mainly based on thermal oil is much more efficient than on aluminum crossed by thermal oil channels. Comparison between the current heat storage and the direct cooker for boiling water even with a standard pot shows that heat storage increases cooking time from 27 to 38 min. He also showed that the selective coating does not drastically improve the efficiency of the system, but the use of an evacuated tube around the absorber reduces by a factor of 2 the charging time of the heat storage.


System description and heat transfer processes

System description and operation

Figure 1 shows the schematics of the experimental system used in this study and described in more details in a previous paper [3]. The system is composed of the following elements: a solar concentrator, a receiver, a heat storage tank, and a circulation pump placed in the primary fluid circuit. Synthetic oil SAE-40 is used as the heat transfer fluid, and the system has a two-axis tracking mechanism. Figure 2 shows pictures of the second experimental setup.

Fig. 1
figure 1

Schematics of a solar cooking system with a heat storage [3]

Fig. 2
figure 2

Detailed pictures of the second experimental setup. a Complete system. b Concentrator and receiver. c Heat storage tank and circulation pump

The key part of the PSCS is the receiver also called the absorber which is composed of two black iron cylinders as shown in Fig. 3. The inner cylinder, with a volume of 1.57 l, has a thickness of 1.5 mm, a diameter of 0.1 m, and a length of 0.2 m and is black acrylic painted to maximize absorption of solar radiation. The outer cylinder is larger with a thickness of 1 mm, a diameter of 0.2 m, and a length of 0.25 m. The absorber is maintained at the focal point by three square sliding iron tube arms expandable from 0.4 to 0.6 m in length. Glass wool is placed between the two cylinders, as insulation to reduce heat losses. The front of the receiver can optionally be equipped with a glass cover. Table 1 gives the system size parameters and optical properties of each material in the two devices.

Fig. 3
figure 3

Longitudinal section of the receiver

Table 1 System size parameters and optical properties

Heat transfer modes

Incident sunlight reaching the parabolic dish is concentrated on the glass cover at the front face of the receiver. A first part of the concentrated solar radiation is reflected to the ambient, the glass absorbs a second part, and a third part is transmitted through the glass cover. The latter part is absorbed by the receiver plate. A small part is reflected back to the glass cover. Heat is transmitted to the fluid via the black-painted metal absorber. The selective coating layer has a high absorptance and a low emittance in order to reduce thermal radiation losses. Figure 4 shows a cross section of the receiver and all heat transfer processes involved. The absorber is considered to be very rigid, and its properties are not affected by temperature.

Fig. 4
figure 4

Cross section of the receiver with heat transfer processes

The thermal resistances and the different modes of heat transfer between the external environment (ambient and sky), the receiver, and the fluid are depicted in Fig. 5. We take into account the energy stored in each node and consider a one-dimensional space temperature variation in the receiver and heat storage tank. The different heat losses, which are conduction through the receiver insulation, convection from the receiver to the ambient air, and radiation from the receiver to the sky, are also considered in a similar manner as Guendouz [7] and Rongrong et al. [8]. For evaluating the different heat loss coefficients, we used the equations presented by Duffie and Beckman [9] and Incropera et al. [10].

Fig. 5
figure 5

Equivalent thermal resistance model

Results and discussion

Governing equations and numerical solution

The present model of the PSCS considers all the abovementioned heat transfer processes, takes into account the presence of the glazing on the receiver, and assumes one-dimensional variations of the temperature along the receiver and the storage tank. The first law of thermodynamics is applied between times t and t + Δt to various system components to obtain the governing energy balance equations in a convenient explicit finite difference form ready for numerical solution.

- Solar concentrator:

$$ {m}_c{C}_c\Delta {T}_c=\left[\left(1-\rho \right){I}_c-{h}_1\left({T}_c-{T}_a\right)-{h}_2\left({T}_c-{T}_s\right)-{h}_3\left({T}_c-{T}_g\right)\right]{A}_c\Delta t $$

- Glass cover:

$$ {m}_g{C}_g\Delta {T}_g=\left[{\upalpha}_g\gamma \rho {A}_c{I}_c+{h}_3\left({T}_c-{T}_g\right){A}_c-\left({h}_4\left({T}_g-{T}_a\right)+{h}_5\left({T}_g-{T}_s\right)\right){A}_g+\left({h}_6+{h}_7\right)\left({T}_p-{T}_g\right){A}_g\right]\Delta t $$

The intercept factor γ accounts for various imperfections in the system operation: shadowing, tracking system inaccuracy, geometry, mirror clearness, dust on the glass cover, and miscellaneous factors [11].

$$ \gamma =\prod_{i=1}^6{\gamma}_i $$

- Receiver plate:

$$ {m}_p{C}_p\varDelta {T}_p=\left[{\alpha}_p{\tau}_g\gamma \rho {A}_c{I}_c-{h}_8\left({T}_p-{T}_{f,r}^1\right){A}_p-\left({h}_6+{h}_7\right)\left({T}_p-{T}_g\right){A}_g\right]\varDelta t $$

- Fluid portion in position 1 of the receiver:

$$ {m}_{r,k}{C}_f\varDelta {T}_{f,r}^1=\left[\dot{m}{C}_f\left({T}_{f,r}^2-{T}_{f,r}^1\right)+{h}_8\left({T}_p-{T}_{f,r}^1\right){A}_p-\frac{\lambda_f}{\varDelta {X}_r}\left({T}_{f,r}^1-{T}_{f,r}^2\right){A}_P-{h}_g\left({T}_{f,r}^1-{T}_a\right){S}_{ur}\right]\varDelta t $$

- Fluid portion in intermediate position k (1 < k < kmax) of the receiver:

$$ {m}_{r,k}{C}_f\varDelta {T}_{f,r}^k=\left[\dot{m}{C}_f\left({T}_{f,r}^{k+1}-{T}_{f,r}^k\right)+\frac{\lambda_f}{\varDelta {X}_r}\left({T}_{f,r}^{k-1}-{T}_{f,r}^k\right){A}_P-\frac{\lambda_f}{\varDelta {X}_r}\left({T}_{f,r}^k-{T}_{f,r}^{\begin{array}{l}\\ {}k+2\end{array}}\right){A}_P-{h}_g\left({T}_{f,r}^k-{T}_a\right){S}_{ur}\right]\varDelta t $$

- Fluid portion in position kmax of the receiver:

$$ {m}_{r,k}{C}_f\varDelta {T}_{f,r}^{k\max }=\left[\dot{m}{C}_f\left({T}_{f,r}^1-{T}_{f,r}^{k\max}\right)+\frac{\lambda_f}{\varDelta {X}_r}\left({T}_{f,r}^{k\max -1}-{T}_{f,r}^{k\max}\right){A}_P-{h}_g\left({T}_{f,r}^{k\max }-{T}_a\right)\left({S}_{ur}+{A}_P\right)\right]\varDelta t $$

Similarly, the storage tank is divided into imax fluid zones.

- Fluid portion in position 1 of the storage tank:

$$ {m}_{r,k}{C}_f\varDelta {T}_{f,r}^1=\left[\dot{m}{C}_f\left({T}_{f,r}^2-{T}_{f,r}^1\right)+\frac{\lambda_f}{\varDelta {X}_r}\left({T}_{f,r}^2-{T}_{f,r}^1\right){A}_p-{h}_g^{\prime}\left({T}_{f,r}^1-{T}_a\right)\left({S}_{ur}+{A}_s\right)\right]\varDelta t $$

- Fluid portion in position i (1 < i < imax) of the storage tank:

$$ {m}_{s,i}{C}_f\varDelta {T}_{f,s}^i=\left[{C}_f\left({T}_{f,s}^{i+1}-{T}_{f,s}^i\right)+\frac{\lambda_f}{\varDelta {X}_s}\left({T}_{f,s}^{i+1}-{T}_{f,s}^i\right){A}_s-\frac{\lambda_f}{\varDelta {X}_s}\left({T}_{f,s}^i-{T}_{f,s}^{i-1}\right){A}_S-{h}_g^{\hbox{'}}\left({T}_{f,s}^i-{T}_a\right){S}_{us}\right]\varDelta t $$

- Fluid portion in position imax of the storage tank:

$$ {m}_{s,i}{C}_f\varDelta {T}_{f,s}^{i\max }=\left[{C}_f\left({T}_{f,r}^1-{T}_{f,s}^{i\max}\right)-\frac{\lambda_f}{\varDelta {X}_s}\left({T}_{f,s}^{i\max }-{T}_{f,s}^{i\max -1}\right){A}_s-{h}_g^{\hbox{'}}\left({T}_{f,s}^{i\max }-{T}_a\right)\left({S}_{us}+{A}_s\right)\right]\varDelta t $$

The global heat loss coefficients h g and \( {h}_g^{\prime } \) are calculated as follows:

$$ {h}_g={h}_g^{\prime }=\frac{1}{\frac{1}{h_i}+\frac{1}{h_e}+\frac{e_{\mathrm{ins}}}{\lambda_{\mathrm{ins}}}} $$

The above governing equations are solved numerically to determine the time variations of the temperatures of the concentrator, the glass cover, the plate on the front face of the receiver, and the fluid at different positions both in the receiver and the storage tank. The simulation parameters for the model are summarized in Table 2.

Table 2 Simulation parameters

Convergence and validation of the numerical solution

The time step is lowered until stability and convergence of the numerical solution is obtained. Figure 6 shows the time variation of the receiver plate temperature at time steps (∆t) 100 and 10 s; the latter coincides with the curves obtained using smaller value time steps (1 and 0.1 s). When the other temperatures are considered, similar behavior is observed, and total stability was obtained with 0.1 s which is the adopted value throughout the present work.

Fig. 6
figure 6

Time variation of the receiver plate temperature at different time steps

Our solution was compared with the results of the simpler model described by Newton [12] which consider that the receiver metal plate and the fluid have the same temperature. Zeghib [13] used the same model, but he considers that there are temperature gradients along the receiver. Figure 7 shows that the temperature difference between the receiver plate and the fluid in position 1 tends to vanish when the heat transfer coefficient is infinite (actually larger than 3000 W m−2 °C−1). Similarly, the temperature difference between the plate and the fluid in position kmax decreases to 0 when the fluid thermal conductivity is infinite (actually larger than 2000 W m−1 °C−1).

Fig. 7
figure 7

Time variation of the temperature difference between the plate and the fluid when we increase a the heat transfer coefficient and b the fluid thermal conductivity

Figure 8 presents the fluid temperature at position kmax in the receiver obtained by solving our model using a high heat transfer coefficient (a heat coefficient of 10,500 W m−2 °C−1 is taken, to guarantee that the plate temperature and the fluid temperature at position k = 1 are equal) and varying the fluid thermal conductivity. At higher values (λ f larger than 6000 W m−1 °C−1), the numerical solution matches the results obtained using the simple model of [12, 13], under the same operating conditions and with the same time step (0.1 s).

Fig. 8
figure 8

Comparison of the present model with the simple model of [12, 13]

The first experimental system was tested in the region of Rabat (Morocco) during the period from April 24 to July 10, 2014. Tests were conducted from 9:00 am to 5:30 pm local time, under clear sky conditions (see Fig. 9).

Fig. 9
figure 9

Time variation of solar radiation and ambient temperature in the closed circuit of the first experimental device

Figure 10 compares the measured and the theoretical fluid temperatures in the upper part of the receiver in the closed circuit of the first experimental device using SAE-40 synthetic oil at 15-min intervals which reached a maximum temperature of 153 °C after 5 h.

Fig. 10
figure 10

Measured and predicted fluid temperatures in the receiver in the closed circuit of the first experimental device

The second experimental system was tested at the same location during the period from May 15 to June 18, 2015, under clear sky conditions (see meteorological data in Fig. 11).

Fig. 11
figure 11

Time variation of solar radiation and ambient temperature in the closed circuit of the second experimental device

The measured and the theoretical fluid temperatures in the upper part of the receiver in the closed circuit of the second experimental device using SAE-40 synthetic oil at 5-min intervals are given in Fig. 12. The maximum fluid temperature was 150 °C reached after 1 h of heating.

Fig. 12
figure 12

Measured and predicted fluid temperatures in the receiver in the closed circuit of the second experimental device

A good agreement is noticed between the theoretical values and the experimental results. The relative error (RE) is between ± 2.4 and ± 4.3%, and the root mean square error (RMSE), which represents the arithmetic mean of the squares of the differences between the forecasts and the observations, is between 1.2 and 3.0 °C. These values are summarized in Table 3.

Table 3 Relative and root mean square errors in the closed circuit

The DNI and temperature data used in an open circuit are shown in Fig. 13.

Fig. 13
figure 13

Time variation of solar radiation and ambient temperature in the open circuit of the second experimental device

Figure 14 shows the oil temperature measurements in the upper part of the receiver and in the lower part of the storage tank using SAE-40 synthetic oil at 15-min intervals. The maximum fluid temperature in the storage was 75 °C.

Fig. 14
figure 14

Measured and predicted fluid temperatures in the receiver and in the storage tank in the open circuit of the second experimental device

Table 4 summarizes the relative error and the root mean square error which were, respectively, between ± 4.0 and 5.9% and between 1.3 and 1.5 °C in the receiver and between ± 4.4 and ± 7.5% and between 1.3 and 1.9 °C in the storage tank. The increase in RE noted on the day of May 27, 2015, is due to the temperature difference that was greater than 5 °C between the receiver and the storage tank due to the blocking of the pump shaft at the end of the afternoon.

Table 4 Relative and root mean square errors in the open circuit

Parametric analysis

The performance of a PSCS can be significantly affected by numerous parameters such as weather conditions (solar radiation, wind) that vary according to the site, material optical properties (reflectance, absorptance, emissivity), system design parameters (aspect ratio, rim angle, intercept factor, exposure ratio), and operating parameters (mass flow rate, glazing, air between the glass and the plate on the front face of the receiver, tracking mechanism, fluid nature, heat losses). The following paragraphs present the results of a parametric analysis of the effect of these parameters on the system performance measured by the maximum fluid temperature in the heat storage tank. This choice is dictated by the fact that this one will ensure the cooking (types of food to be cooked and cooking time).

Effect of maximal solar radiation

Figure 15 shows that the maximum heat storage temperature is greatly affected by maximal solar radiation, which may change depending on the season and geographic location. The daily maximal radiation varies between 500 and 950 W m−2, and a change of daily maximal solar radiation by 50 W m−2 increases the maximum storage temperature by about 4 °C. It should be noted that thermal losses also increase due to the increased receiver temperature. However, this increase is smaller than the enhanced absorbed solar energy. Rongrong et al. [8] reported that when solar radiation changes by 20 W m−2, the collector output thermal fluid decreases or increases by about 2 °C. Luo et al. [14] found when solar irradiation increases from 900 to 1100 W m−2 and decreases from 900 to 700 W m−2, the collector outlet temperature will increase or decrease by about 10 °C.

Fig. 15
figure 15

Maximum storage temperature as a function of solar radiation

Effects of the material reflectance and absorptance

Aluminum (0.7), silver mirrors (0.7 to 0.9), and Mylar (0.94) are commonly used as reflectors, whereas copper coated with black paint (0.75) and stainless steel coated with cermet (mixture of ceramics and metal) or black chrome (0.9 to 0.94) are used as absorbers. Figure 16 illustrates the effect of the reflectance and absorptance of the materials used on the system performance. For a reflectance of 0.75, an increase of the absorptance from 0.6 to 0.75 will cause an 18 °C increase of the heat storage temperature. Generally, an increase of 5% of the reflectance or the absorptance will increase the storage temperature by 3.6 and 3.9 °C, respectively.

Fig. 16
figure 16

Maximum storage temperature as a function of the receiver reflectance and absorptance

Effect of the aspect ratio

The aspect ratio defined as cavity length to diameter of the receiver is a design parameter which influences system performance. A large diameter increases the “dead area” created by the receiver shadow on the concentrator and the heat losses to ambient air, while a longer cylindrical creates a temperature gradient with slower heat transfer. For a given volume, it is crucial to optimize the aspect ratio as indicated in Fig. 17. The maximum heat storage temperature (85 °C) is reached when the aspect ratio of the receiver is near to 2. Beyond this value, the temperature decreases because the effect of heat losses particularly radiation in the receiver is greater than the amount of solar radiation intercepted.

Fig. 17
figure 17

Maximum storage temperature as a function of the aspect ratio L/d

This behavior is close to the results found by other authors [15,16,17] who also showed increased radiation loss with increasing receiver aperture diameter. Beltran et al. [16] found that the aperture diameter of the cavity reaches a point of maximum efficiency when the aperture diameter is equal to 0.13 m which corresponds to a receiver aspect of ratio of 1.5. Gil et al. [17] found with a heat loss coefficient of 12 W m−2 °C−1, the receiver best aspect ratio is 1.9 which is close to our case. Prakash et al. [18] found that the convective loss values increase with opening ratio or exposure ratio while the aspect ratio is equal to 1. The higher the opening ratio, the greater the convective zone leading to higher convective losses. An increase of about 30–50% in the convective loss values is observed when the opening ratio increases from 0.5 to 1.

Paitoonsurikarn et al. [19] used in their work an aspect ratio of 2.2 for the base case and show that the heat flux increases with decreasing cavity aspect ratio for different inclinations. Madadi et al. [20] used many receivers which have an aspect ratio of 2, 2.8, and 3.6 and show that the effect of receiver temperature on radiation loss for greater apertures is higher than smaller apertures. The radiation heat losses from the receiver in comparison with convection heat losses are so low. By decreasing the heat transfer fluid mass flow rate from 0.1 to 0.0083 kg s−1, with a receiver with an aspect ratio of 2, the average receiver temperature increases from 197 to 310 °C, the convective heat loss increases only 22% while the radiation heat loss increases up to 165%.

Abbasi-Shavazia et al. [21] used a receiver which has an aspect ratio of 2, and the comparison of radiation loss values based on bottom surface and area-weighted cavity average temperature shows no more than a 10–20% difference.

Effect of the thermal fluid mass flow rate

In high-temperature solar thermal systems with energy storage and automatic control, it is known that fluid mass flow rate is a very important operating parameter. A low mass flow rate causes higher temperatures in the receiver and thus bigger heat losses and lower storage tank temperatures. Conversely, with higher mass flow rate, thermal energy is carried away much faster from the receiver, but also increases the pump size and power consumption. It is therefore important to determine the “optimal” value that yields the highest storage temperature without oversizing the pump. In our case, Fig. 18 shows that increasing the fluid flow rate from 0 to 18 kg h−1 increases the storage temperature by 65 °C, and beyond this value, there is no significant improvement of the system performance. These results are very similar to those found by Rongrong et al. [8], Luo et al. [14], and Madadi et al. [20]. The outlet oil temperature in the receiver will be reduced when increasing flow rate in the same working condition.

Fig. 18
figure 18

Maximum storage temperature as a function of the thermal fluid mass flow rate

They found a curve which decreases when the flow rate increases. In this case, the curve of the receiver thermal efficiency presents a maximum at an optimal mass flow rate. The difference between fluid inlet and outlet temperatures is controlled and maintained at 5 °C in our case, which is commonly used in solar thermal systems.

Effect of insulation thickness

Heat losses must be carefully analyzed as their assessment can significantly affect the overall system performance accuracy. Thermal insulation of the receiver deserves a careful analysis as a thicker receiver wall reduces thermal losses but it increases the shadow area on the solar concentrator, as well as the receiver weight and system cost.

The receiver wall thermal resistance is the inverse of the global heat transfer coefficient resulting from conduction through the receiver and convection heat transfer on the interior and the exterior of the receiver.

An insulation thickness in the receiver between 10 and 80 mm, which corresponds to a global heat loss coefficient between 0.5 and 4.0 W m−2 °C−1 (or a receiver wall thermal resistance between 0.25 and 2 m2 °C W−1), is studied. Increasing the insulation thickness from 0.01 to .08 m improves the maximum storage temperature by 17 °C (see Fig. 19). To minimize conduction in the receiver, an insulation thickness of 75 mm has been suggested as an effective width by Fraser [15].

Fig. 19
figure 19

Maximum storage temperature as a function of receiver wall thermal resistance


The present model of parabolic solar cooking systems introduced with heat storage for continuous use allowed a valuable analysis of the performance of such systems. Improvement over previous simpler models included a non-uniform receiver temperature and temperature difference between the glass, receiver cover, and thermal fluid. The model-governing equations were solved using an explicit finite difference method, and the method was mathematically validated. The results of the simulation were compared with experimental results, which proved that the model predicts adequately the thermal behavior of the described system with a relative error ± 4.4% and a root mean square error of 3 °C. The model was used to carry out a detailed parametric study of the main design and operating parameters which affect the system energy performance. The results show that a change of daily maximum solar radiation of 50 W m−2 increases the heat storage temperature by about 4 °C and an increase of 5% of the reflectance or the absorptance improves the heat storage temperature by 3.6 and 3.9 °C, respectively. It was also shown that the best aspect ratio of the receiver is 2. Increasing the fluid mass flow rate from 0 to 18 kg h−1 leads to a maximum storage temperature improvement of 65 °C. Going from an insulation thickness of 0.01 to 0.08 m increases the maximum storage temperature by 17 °C.


A c

Concentrator surface area


A p

Receiver plate surface area


A g

Glass surface area


A s

Storage cross-section area



Concentrator specific heat

kJ kg−1 °C−1


Glass specific heat

kJ kg−1 °C−1

C f

Fluid specific heat

kJ kg−1 °C−1

C p

Receiver plate specific heat

kJ kg−1 °C−1


Receiver diameter



Insulation thickness


h 1

CHTC concentrator to ambient

W m−2 °C−1

h 2

RHTC concentrator to sky

W m−2 °C−1

h 3

RHTC concentrator to glass

W m−2 °C−1

h 4

CHTC glass to ambient

W m−2 °C−1

h 5

RHTC glass to sky

W m−2 °C−1

h 6

CHTC plate to glass

W m−2 °C−1

h 7

RHTC plate to glass

W m−2 °C−1

h 8

CHTC plate to fluid

W m−2 °C−1


GHTC receiver to ambient

W m−2 °C−1

\( {\mathrm{h}}_g^{\prime } \)

GHTC storage to ambient

W m−2 °C−1

I c

Direct normal irradiance

W m−2


Receiver length perpendicular to the aperture



Thermal fluid mass flow rate

kg s−1

m c

Concentrator mass


m g

Glass mass


m p

Receiver plate mass


m r , k

Receiver partial mass in position k


m s , i

Storage partial mass in position i


S ur

Receiver partial lateral surface area


S us

Storage partial lateral surface area


T a

Ambient temperature


T c

Concentrator temperature


T f

Fluid temperature


T g

Glass temperature


T p

Receiver plate temperature


T s

Sky temperature


Greek symbols


Glass absorptance



Receiver plate absorptance



Intercept factor



Time step


X r

Receiver length step


X s

Storage length step


λ f

Fluid thermal conductivity

W m−1 °C−1

λ ins

Insulation thermal conductivity

W m−1 °C−1


Mirror reflectance


τ g

Glass transmittance










Time step








  1. ASDER (2012) Solar cookers—technical presentation. Savoyard Association for the Development of Renewable Energies, Chambery.

  2. Prasanna UR (2010) Modeling, optimization and design of a solar thermal energy transport system for hybrid cooking application. Ph.D. thesis, Indian Institute of Science.

  3. Mbodji N, Hajji A (2016) Performance testing of a parabolic solar concentrator for solar cooking. ASME J Sol Energy Eng 138(4). doi:10.1115/1.4033501

  4. Mawire A, McPherson M, Van den Heetkamp RRJ (2008) Simulated energy and exergy analyses of the charging of an oil-pebble bed thermal energy storage system for a solar cooker. Sol Energy Mater Sol Cells 92:1668–1676. doi:10.1016/j.solmat.2008.07.019

    Article  Google Scholar 

  5. Mawire A, McPherson M, Van den Heetkamp RRJ (2010) Discharging simulations of a thermal energy storage (TES) system for an indirect solar cooker. Sol Energy Mater Sol Cells 94(6):1100–1106. doi:10.1016/j.solmat.2010.02.032

    Article  Google Scholar 

  6. Mussard M (2013) A solar concentrator with heat storage and self-circulating liquid. Ph.D. thesis, Norwegian University of Science and Technology, Trondheim, Norway.

  7. Guendouz B (2012) Use of solar energy for air conditioning needs. Chapter 3: state of the art of solar collectors and modeling. Magister’s thesis, Abou Bakr Belkaid University, Tlemcen, Algeria.

  8. Rongrong Z, Yongping Y, Qin Y, Yong Z (2013) Modeling and characteristic analysis of a solar parabolic trough system: thermal oil as the heat transfer fluid. Journal of Renewable Energy. doi:10.1155/2013/389514

  9. Duffie JA, Beckman WA (1980) Solar engineering of thermal processes, 2nd edn. John Wiley & Sons, Inc, New York

    Google Scholar 

  10. Incropera F, Dewitt DP, Bergman TL, Lavine AS (2007) Fundamentals of heat and mass transfer, 6th edn. Wiley, New York.

  11. Padilla RV (2011) Simplified methodology for designing parabolic trough solar power plants. Ph.D. thesis, South Florida University.

  12. Newton CC (2006) A concentrated solar thermal energy system. Master’s thesis, Florida State University, Tallahassee, FL.

  13. Zeghib I (2005) Design and construction of a parabolic solar concentrator. Master’s thesis. Mentouri – Constantine University, Algeria.

  14. Luo N, Yu G, Hou HJ, Yang YP (2015) Dynamic modeling and simulation of parabolic trough solar system. Energy Procedia 69:1344–1348. doi:10.1016/j.egypro.2015.03.137

    Article  Google Scholar 

  15. Fraser PR (2008) Stirling dish system performance prediction model. Master of Science, Wisconsin.

  16. Beltran R, Velazquez N, Espericueta AC, Sauceda D, Perez G (2012) Mathematical model for the study and design of a solar dish collector with cavity receiver for its application in Stirling engines. J Mech Sci Technol 26(10):3311–3321. doi:10.1007/s12206-012-0801-0

    Article  Google Scholar 

  17. Gil R, Monné C, Bernal N, Muñoz M, Moreno F (2015) Thermal model of a dish Stirling cavity-receiver. Energies 8:1042–1057. doi:10.3390/en8021042

    Article  Google Scholar 

  18. Prakash M, Kedare SB, Nayak JK (2009) Investigations on heat losses from a solar cavity receiver. Sol Energy 83(2):157–170. doi:10.1016/j.solener.2008.07.011

    Article  Google Scholar 

  19. Paitoonsurikarn S, Lovegrove K, Hughes G, Pye J (2011) Numerical investigation of natural convection loss from cavity receivers in solar dish applications. ASME J Sol Energy Eng 133(2):1–10. doi:10.1115/1.4003582

    Article  Google Scholar 

  20. Madadi V, Tavakoli T, Rahimi A (2015) Estimation of heat loss from a cylindrical cavity receiver based on simultaneous energy and exergy analyses. J Non-Equilib Thermodyn 40(1):49–61. doi:10.1016/j.solener.2008.07.011

    Article  Google Scholar 

  21. Abbasi-Shavazia E, Hughes GO, Pye JD (2015) Investigation of heat loss from a solar cavity receiver. Energy Procedia 69:269–278. doi:10.1016/j.egypro.2015.03.031

    Article  Google Scholar 

Download references

Author information

Authors and Affiliations



NM worked as a PhD student under the supervision of AH, wrote the paper, and did the field tests. AH supervised the research work, gave guidance on modeling, and corrected the English text. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Ndiaga Mbodji.

Ethics declarations

Competing interests

This work has received no external funding. I declare that there are no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Mbodji, N., Hajji, A. Modeling, testing, and parametric analysis of a parabolic solar cooking system with heat storage for indoor cooking. Energ Sustain Soc 7, 32 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: