KnE Materials Science | Sino-Russian ASRTU Conference Alternative Energy: Materials, Technologies, and Devices | pages: 146-158

, , , and

1. Introduction

Volumetric receiver with high-porosity material appears to be a promising technology for converting solar flux into thermal energy [1]. As one of the most important heat transfer processes in the receiver, the solar radiation transmission and absorption process is important to the overall system performance [2]. The thermal behavior and efficiency of the solar receiver have been extensively investigated. The incoming concentrated solar radiation is modeled using two major approaches. In the first approach, the heat flux at the boundary is considered as a surface phenomenon based on the assumption of very large optical thickness for the porous media. The solar flux distribution is treated as the secondary type thermal barrier coating (TBC) in the simulation [3]. This assumption neglects the gradual absorption characteristics of the solar radiation in the receiver and cannot reflect the volumetric effect.

The second approach incorporates the solar radiation transport within the volumetric receiver, and the concentrated solar radiation is considered as a radiative heat source in the whole volume of receiver [4]. Nevertheless, instead of appropriate modeling, the concentrated solar flux distribution, the impinging solar radiation on the front surface of receiver has been treated as a Cox-Ingersoll-Ross (CIR) beam with a uniform or Gaussian distribution in several studies. In fact, the concentrated solar radiation is extremely nonuniform in spatiality and direction. This distribution has a significant impact on the temperature field in the solar receiver [5]. Therefore, the performance evaluation of the solar receiver should be based on the coupling between solar flux modeling and heat transfer modeling [6]. Monte Carlo Ray Tracing (MCRT) simulations have been carried out to obtain the realistic concentrated heat flux distribution boundary conditions by some researchers.

An integrated model of solar radiation propagation in a concentrator-window-absorber system is established. The local solar radiative source distribution within the porous absorber and the optical efficiency of the system are predicted with the MCRT method. The flow and heat transfer simulation of a solar receiver with double-layer ceramic foam is performed in this study. The effects of geometric parameters of each porous layer on the thermal performance are mainly discussed.

2. Methods

Model description

A key component in the solar thermal system is the volumetric receiver, which is located at the focal plane of the concentrator. As shown in Figure 1, the solar radiation is collected and redirected to the receiver by the dish concentrator. The surface equation for the dish is:


(1) where f is the focal length of the dish.

Figure 1

Schematic of the volumetric receiver with solar dish concentrator system.


The optical errors may enlarge the solar image and reduce the optical performance of system [7]. For a solar concentrator, the sources of optical errors are typically a slope error σslope , a tracking error σtrack , a nonspecular reflection σscatt , and some alignment error σalign [8]:



The heat transfer model for numerical analysis is shown in Figure 2. A local cylindrical coordinate system o-rφz is established for the porous media. Besides, the reflection and absorption also occurs at the receiver wall. Geometrical and physical parameters of the concentrator and volumetric solar receiver are presented in Table 1.

Figure 2

Schematic diagram of heat transfer in the volumetric solar receiver across a gray slab.

Table 1

Geometrical and physical parameters of the concentrator and volumetric solar receiver.

Concentrator and Volumetric Receiver Value
Focal length of the dish concentrator (f ) 3.25 m
Aperture radius of the dish concentrator (R C ) 1.3 m
Length of the porous media (L) 0.05 m
Diameter of the receiver (2R) 0.05 m
Thickness of the plane window (h g ) 8 mm
Distance between the window and the porous media (h 1 ) 5 mm
Refractive index of the window (n g ) 1.42
Absorption coefficient of the window (k g ) 1.4 m -1
Emissivity of the porous media (ε) 0.92
Emissivity of the receiver wall (ε w ) 0.3
Reflectivity of the dish surface (ρ C ) 0.9
Slope error of the dish concentrator (σ slope ) 2.0 mrad

Solar radiation transport simulation

MCRT method is a random simulation method based on the probability statistics and is widely used in concentrated solar research [9]. The pillbox distribution is chosen to describe the sunshape effect, which means that the solar radiation is constant within a cone with a half angle of θs=4.65 mrad [10]. Besides, the value of the solar radiation heat flux incident on the concentrator used in this study is qs=1000 W/m 2 [11]. The slope error σslope is applied to indicate the deviation of a real dish surface from a perfect one. In MCRT process, the slope deviation can be defined by a azimuth angle φse and a zenith angle θse , which can be determined by the following expressions [12]:




(4) where Rsθ and Rsφ are random numbers, which are uniformly distributed between 0 and 1.

Once the ray penetrates into the porous media, based on media radiation transfer theory, a possible transfer distance may be obtained from the following probability expression [13]:


(5) where β is the extinction coefficient and Rβ is a random number, which ranges from 0 to 1 uniformly. When the location lβ is reached, the second random number Rω is required to decide whether the ray is absorbed or scattered, with ω being the scattering albedo ( ω=σs/β , where σs is the scattering coefficient.) [14]:

whiteRωω, scattering whiteRω>ω, absorption

(6) If the ray is scattered, it will travel into a new direction. Otherwise, the ray is absorbed by the local element. The solar radiative source of the element i caused by the concentrated solar radiation can be computed as:


(7) where Ni denotes the number of rays absorbed within the volume element, ΔVi is the volume of the element, which can be determined based on the element division, and e represents the energy carried by each ray sampling, e=qs/nray , where nray is the ray sampling number per unit area.

Heat transfer simulation within the receiver

The governing equations are described as follows.

Continuity equation


where ρf is the fluid density and V denotes the superficial velocity.

Momentum equation

Flow in the porous media is modeled by adding a momentum source term to the standard fluid flow equation, which can be expressed as [4]:


(6) where p is the pressure of fluid, μf is the dynamic viscosity, ϕ is the porosity, and F is the momentum source term to calculate the pressure drop resulting from the porous media, which is calculated by the following equation [15]:


(7) where dc is the mean cell size. This equation is valid for 0.66<ϕ<0.93 , 10<Re<400 ( Re=ρfudc/μf ) and a cross-section of the pore channel that approaches to a circle.

Energy equation

The energy equations with LTNE condition for the fluid phase and solid phase can be expressed as follows [16].

For the fluid phase:


For the solid phase:


where Tf and Ts are the fluid and solid temperatures, respectively, cp is the specific heat of fluid, Sr is the volumetric heat source term due to radiative heat transfer, and the effective thermal conductivity for the fluid phase and solid phase could be determined using the Schuetz–Glicksman empirical formulas [9]:


The symbol hv is the volumetric convection heat transfer coefficient between the fluid phase and the solid phase. The empirical correlation proposed by Wu et al. [17] is used:


This correlation is valid for 0.66<ϕ<0.93 and 70<Re<800 .

Boundary conditions

The front surface of the porous media is subjected to the concentrated solar radiation. The boundary conditions are as follows:









Thermophysical parameters

The porous media is assumed to be SiC foam with isotropic properties. The thermal properties are assumed to be constant. The thermal conductivity is 80 W/(m·K). The thermal capacity and density of bulk SiC are 750 J/(kg·K) and 3200 kg/m 3 , respectively [4]. Air is regarded as ideal gas and the viscosity is calculated through Sutherland Law. The thermal capacity and conductivity can be described as functions of temperature by the following polynomial functions [15]:





3. Results

MCRT model validation

To examine the validity of the numerical simulation, two comparative studies have been carried out. First, for the solar flux modeling, the solar flux distributions at the focal plane for two dish concentrators with the case of f=5.0 m at σslope=2.0 mrad are simulated. The results are compared with the data obtained by Lee et al. [18], where ψrim is the rim angle. To further validate the radiative transfer simulation using the P1 approximation, the results are compared with those in [14].

Heat transfer results of receiver

Effects of porous structure parameters

The porosity and pore diameter are the two main parameters to characterize the porous structure. The impacts of porosity and pore diameter on the receiver performance are studied, respectively. A local coordinate system O-XYZ is aligned with the porous absorber, where the XY plane is set on the bottom surface of the absorber.

According to the simulation, it can be seen that the radiative source decreases steeply along the center line and the maximum value decreases with the porosity increasing. The optical performance can be evaluated for different porosities, as illustrated in Table 2.

Table 2

Optical performance of the solar system with different porosities.

ϕ η opt (%) η sw (%) η bs (%) η con (%)
0.7 72.35 7.12 10.53 10.0
0.8 72.35 7.12 10.53 10.0
0.9 72.36 7.12 10.52 10.0
0.95 72.34 7.12 10.54 10.0

On the other hand, the extinction coefficient decreases with an increase in the pore diameter. The optical performance with different pore diameters is given in Table 3. The optical efficiency and optical losses have no noticeable change for different pore diameters. Besides, the structure parameters affect the distribution of absorbed solar energy within the absorber greatly and have no strong effect on the optical efficiency and optical losses.

Table 3

Optical performance of the solar system with different pore diameters.

d p (mm) η opt (%) η sw (%) η bs (%) η con (%)
1.5 72.35 7.12 10.53 10.0
2.0 72.36 7.12 10.52 10.0
2.5 72.37 7.12 10.51 10.0
3.0 72.37 7.12 10.51 10.0

Effects of slope error of the solar concentrator

Slope error of the concentrator is a crucial factor in determining the concentrated solar flux and direction distribution affecting the receiver aperture. Figure 3 exhibits the effects of slope error on the distribution of solar radiative source and the optical performance. The optical efficiency decreases as the slope error increases, while the optical loss of dish concentrator increases gradually from 10.0% to 58.55%. The optical efficiency is 72.36% for σslope=0 mrad, while it is 31.62% for σslope=3.0 mrad.

Figure 3

Distribution of the solar radiative source along the center line of porous absorber with different slope errors.


Simulation methods

Comparison with the approach using TBC

To solve the conservation equations, FLUENT software with user defined functions is used implementing a two-dimensional axisymmetric model. The inlet velocity is set to 1.5 m/s. In this section, the TBC approach stating that the concentrated solar radiation is considered as thermal boundary condition is discussed.

The maximum temperature occurs at the front surface and the solid temperature decreases along the z-axis direction with the TBC approach. Since the extinction coefficient decreases reversely with the porosity, solar radiation can penetrate into the receiver much deeper with a high porosity. Besides, the heat absorbed in the front region is removed by convection effectively [19]. These factors result in the temperature distribution. The simplified TBC approach neglects the penetration of solar radiation into the receiver. This approximation may be suitable for a very low porosity with small mean cell size, since the extinction coefficient is extremely large. Thus, this TBC approach is not appropriate for predicting the temperature distribution in the receiver, especially in the front region.

Comparison with the approach using CIR

The corresponding percent distribution of incident solar radiation on the front surface along the zenith angle θ is presented in Figure 4. It can be seen that the maximum incident angle is about 23 , which contrasts with the approximation of collimated incidence. The absorbed solar radiation becomes concentrative in the center area near the front surface, and then the direction is redistributed due to the scattering of the porous media. Since isotropic scattering is considered, equal amounts of energy are scattered into all directions. As a result, the temperature away from the front surface shows no discernible difference between the two approaches. Therefore, the CIR approach can get acceptable results compared to the coupling simulation.

Figure 4

Directional characteristics of incident solar radiation on the front surface of the porous receiver.


Effects of slope error on the thermal performance

For a solar concentrator, the concentrating performance is significantly affected by the optical errors. The sources of optical errors are typically a slope error, a tracking error, nonspecular reflection, and some alignment error [8]. The predominant one is the slope error ranging from 1 to 6 mrad [8]. The previous studies mostly focus on the effects of characteristic parameters of the receiver on the thermal performance, whereas the effect of concentrating performance is excluded. Slope error mainly causes the reflected ray to deviate from the specular reflection direction and eventually from the receiver. Figure 5 demonstrates that the slope error has a great influence on the temperature distribution.

Figure 5

Effect of the slope error on the temperature distribution: (a) ϕ=0.8 , dc=2.0 mm (b) ϕ=0.9 , dc=2.0 mm.


4. Conclusion

In this study, a numerical model of volumetric solar receiver is presented by coupling the solar flux modeling and heat transfer modeling. As the first step, the solar radiation transport from concentrator to the interior of receiver is simulated with the MCRT method. The following conclusions have been drawn:

  • With the porosity increasing, the maximum radiative source within the porous absorber decreases and the distribution becomes uniform gradually. Moreover, the effect of porosity on the optical efficiency and losses is negligible. The effects of pore diameter on the solar radiative source and optical performance are similar to those of porosity.

  • The TBC approach obviously overestimates the solid temperature in the front region of receiver, while the outlet air temperature is underestimated. The temperature deviation increases as the porosity or mean cell size increases. The deviations of inlet solid and outlet air temperatures are up to 76.4% and 13.2%, respectively.

  • Using the CIR approach, the temperature distribution is almost the same as the result of coupling simulation except that the solid temperature near the front surface of receiver is slightly underestimated. The maximum deviation in this study is an acceptable value of 3.4%.

  • The slope error of solar concentrator has a remarkable influence on the temperature field in receiver. The fluid and solid temperatures both greatly decrease as the slope error increases. For example, the outlet air temperature decreases from 1053.9 K to 651.2 K as the slope error increases from 2.0 to 4.0 mrad for ϕ=0.8 and dc=2.0 mm.


These investigations are supported by the National Natural Science Foundation of China (No. 51536001) and the Fundamental Research Project of China Industrial and Information Ministry (No. B2320132001).



Roldán M. I., Zarza E., Casas J. L., Modelling and testing of a solar-receiver system applied to high-temperature processes, Journal of Renewable Energy, Year: 2015, Volume: 76, Page: 608-618. DOI: 10.1016/j.renene.2014.11.075


Cui F. Q., He Y. L., Cheng Z. D., Li D., Tao Y. B., Numerical simulations of the solar transmission process for a pressurized volumetric receiver, Energy, Year: 2012, Volume: 46, Issue: 1, Page: 618-628. DOI: 10.1016/


Wang F., Shuai Y., Tan H., Zhang X., Mao Q., Heat transfer analyses of porous media receiver with multi-dish collector by coupling MCRT and FVM method, Solar Energy, Year: 2013, Volume: 93, Page: 158-168. DOI: 10.1016/j.solener.2013.04.004


Villafán-Vidales H. I., Abanades S., Caliot C., Romero-Paredes H., Heat transfer simulation in a thermochemical solar reactor based on a volumetric porous receiver, Applied Thermal Engineering, Year: 2011, Volume: 31, Issue: 16, Page: 3377-3386. DOI: 10.1016/j.applthermaleng.2011.06.022


Wang F., Shuai Y., Tan H., Yu C., Thermal performance analysis of porous media receiver with concentrated solar irradiation, International Journal of Heat and Mass Transfer, Year: 2013, Volume: 62, Issue: 1, Page: 247-254. DOI: 10.1016/j.ijheatmasstransfer.2013.03.003


Lee H.-J., Kim J.-K., Lee S.-N., Kang Y.-H., Consistent heat transfer analysis for performance evaluation of multichannel solar absorbers, Solar Energy, Year: 2012, Volume: 86, Issue: 5, Page: 1576-1585. DOI: 10.1016/j.solener.2012.02.020


Shuai Y., Xia X., Tan H., Radiation performance of dish solar concentrator/cavity receiver systems, Solar Energy, Year: 2008, Volume: 82, Issue: 1, Page: 13-21. DOI: 10.1016/j.solener.2007.06.005


Buie D., Dey C. J., Bosi S., The effective size of the solar cone for solar concentrating systems, Solar Energy, Year: 2003, Volume: 74, Issue: 5, Page: 417-427. DOI: 10.1016/S0038-092X(03)00156-7


Delatorre J., Baud G., Bézian J. J., Blanco S., Caliot C., Cornet J. F., Coustet C., Dauchet J., El Hafi M., Eymet V., Fournier R., Gautrais J., Gourmel O., Joseph D., Meilhac N., Pajot A., Paulin M., Perez P., Piaud B., Roger M., Rolland J., Veynandt F., Weitz S., Monte Carlo advances and concentrated solar applications, Solar Energy, Year: 2014, Volume: 103, Page: 653-681. DOI: 10.1016/j.solener.2013.02.035


Hasuike H., Yoshizawa Y., Suzuki A., Tamaura Y., Study on design of molten salt solar receivers for beam-down solar concentrator, Solar Energy, Year: 2006, Volume: 80, Issue: 10, Page: 1255-1262. DOI: 10.1016/j.solener.2006.03.002


Johnston G., Focal region measurements of the 20 m2 tiled dish at the Australian National University, Solar Energy, Year: 1998, Volume: 63, Issue: 2, Page: 117-124. DOI: 10.1016/S0038-092X(98)00041-3


Dai G.-L., Xia X.-L., Sun C., Zhang H.-C., Numerical investigation of the solar concentrating characteristics of 3D CPC and CPC-DC, Solar Energy, Year: 2011, Volume: 85, Issue: 11, Page: 2833-2842. DOI: 10.1016/j.solener.2011.08.022


Modest M. F., Radiative Heat Transfer (third edition, Year: 2013, San DiegoAcademic Press


Vafai K., Handbook of Porous Media (second edition)., Porland, Taylor and Francis, Year: 2005,


Wu Z., Caliot C., Flamant G., Wang Z., Coupled radiation and flow modeling in ceramic foam volumetric solar air receivers, Solar Energy, Year: 2011, Volume: 85, Issue: 9, Page: 2374-2385. DOI: 10.1016/j.solener.2011.06.030


Fuqiang W., Jianyu T., Lanxin M., Yong S., Heping T., Yu L., Thermal performance analysis of porous medium solar receiver with quartz window to minimize heat flux gradient, Solar Energy, Year: 2014, Volume: 108, Page: 348-359. DOI: 10.1016/j.solener.2014.07.016


Wu Z., Caliot C., Flamant G., Wang Z., Numerical simulation of convective heat transfer between air flow and ceramic foams to optimise volumetric solar air receiver performances, International Journal of Heat and Mass Transfer, Year: 2011, Volume: 54, Issue: 7-8, Page: 1527-1537. DOI: 10.1016/j.ijheatmasstransfer.2010.11.037


Lee H., The geometric-optics relation between surface slope error and reflected ray error in solar concentrators, Solar Energy, Year: 2014, Volume: 101, Page: 299-307. DOI: 10.1016/j.solener.2013.12.035


Kribus A., Grijnevich M., Gray Y., Caliot C., Parametric study of volumetric absorber performance, Proceedings of the International Conference on Solar Power and Chemical Energy Systems, SolarPACES 2013September 2013USA Page: 408-417.



  • Downloads 13
  • Views 62



ISSN: 2519-1438