## Abstract

A volumetric solar receiver receives the concentrated radiation generated by a large number of heliostats. Turbulent heat transfer occurs from the solid matrix to the air as it passes through the porous receiver. Such combined heat transfer within the receiver, including radiation, convection and conduction, is studied using a local thermal non-equilibrium model. Both the Rosseland approximation and the P1 model are applied to consider the radiative heat transfer through the solar receiver. Furthermore, the low Mach approximation is exploited to investigate the compressible flow through the receiver. Analytic solutions are obtained for the developments of air and ceramic temperatures as well as the pressure along the flow direction. Since the corresponding fluid and solid temperature variations generated under the Rosseland approximation agree fairly well with those based on the P1 model, the Rosseland approximation is used for further analysis. The results indicate that the pore diameter must be larger than its critical value to obtain high receiver efficiency. Moreover, it has been found that optimal pore diameter exists for achieving the maximum receiver efficiency under the equal pumping power. The solutions provide effective guidance for a novel volumetric solar receiver design of silicon carbide ceramic foam.

### Keywords

- turbulent heat transfer
- thermal non-equilibrium
- Rosseland approximation
- P1 model
- volumetric solar receiver
- porous media
- ceramic foam

## 1. Introduction

A solar volumetric receiver is required to have the resistance to temperature as high as 1000 degree Celsius, high porosity for sufficiently large extinction volume such that the concentrated solar radiation penetrates through the receiver, high cell density to achieve large specific surface area and sufficiently high effective thermal conductivity to avoid possible thermal spots. Extruded monoliths with parallel channels (i.e. honeycomb structure) are being used in some solar power plants in Europe, including the solar power tower plant of 1.5 MW built in 2009, in Julich in Germany [1, 2]. However, in such conventional receivers, both thermal spots [3] and flow instabilities [4] have been often reported. In the monolith receiver, locally high solar flux leads to a low mass flow with high temperature, whereas locally low solar flux leads to a high mass flow with low temperature. This causes the absorber material to exceed the design temperature locally, which then leads to its destruction although the average temperature is comparatively low. These difficulties encountered in the receiver must be overcome to run the power plant safely.

In consideration of these requirements, ceramic foams have come to draw attention as a possible candidate to replace the conventional extruded monoliths with parallel channels. Many researchers including Becker et al. [4], Fend et al. [5] and Bai [6] focused on porous ceramic foams as a promising absorber material. Recently, Sano et al. [7] carried out a local non-thermal equilibrium analysis to investigate the receiver efficiency under the equal pumping power. For the first time, the complete set of analytical solutions based on the two-energy equation model of porous media was presented, so as to fully account for the combined effects of tortuosity; thermal dispersion and compressibility on the convective, conductive and radiative heat transfer within a ceramic foam receiver. In their analysis, however, the Rosseland approximation was applied to account for the radiative heat transfer through the solar receiver. It is well known that the Rosseland approximation ceases to be valid near boundaries. Although no wall boundaries exist for the case of the one-dimensional analysis of the solar volumetric receiver, the validity of applying the Rosseland approximation near the inlet boundary of the receiver has not been investigated yet. Furthermore, the effects of turbulence mixing on the heat transfer were not considered.

In this study, the validity of the Rosseland approximation [7] will be examined by comparing the results based on the Rosseland approximation and the results obtained from solving the irradiation transport equation based on the P1 model. The set of the equations will be reduced to a fifth-order ordinary differential equation for the air temperature. Once the air temperature distribution is determined, the pressure distribution along the flow direction can readily be estimated from the momentum equation with the low Mach approximation. Thus, the receiver efficiency, namely, the ratio of the air enthalpy flux increase to the concentrated solar heat flux, can be compared under the equal pumping power, so as to investigate the optimal operating conditions. Some analytical and numerical investigations [3, 4, 5, 6, 7, 8] have been reported elsewhere. However, none of them appeared to elucidate well the combined effects of turbulence, compressibility, radiation, convection and conduction within the volumetric receiver on the developments of air and ceramic temperatures as well as the pressure along the flow direction. This study appears to be the first to provide the complete set of analytical solutions based on the two-energy equation model of porous media [9], fully accounting for the combined effects of turbulence, tortuosity, thermal dispersion, compressibility and radiative heat transfer within a ceramic foam receiver.

## 2. Volume averaged governing equations

As illustrated in Figure 1, the structure of silicon carbide ceramic foam volumetric receiver may be treated as homogeneous porous medium. Since the dependence of the Darcian velocity on the transverse direction can only be observed in a small region very close to the walls of the passage, we may neglect the boundary effects (i.e. Brinkman term).

Based on a theoretical derivation of Darcy’s law, Neuman [10] pointed out that the application of Darcy’s law to compressible fluids is justified as long as Knudsen numbers are sufficiently small to ensure the no-slip conditions at the solid–gas interface. This is usually the case for the volumetric receivers. Thus, allowing the density to vary through the receiver, the following Forchheimer extended Darcy law should hold:

where *K* and *b* are the permeability and the inertial coefficients, respectively. Furthermore, by virtue of the volume averaging procedure [11, 12, 13], the microscopic energy equations of the compressible fluid flow phase and the solid phase may be integrated over an elemental control volume *V*, so as to derive the corresponding macroscopic energy equations. Since the porous medium is considered to be homogeneous, the integration of the two distinct energy equations gives:

For the air:

For the solid matrix:

where the intrinsic volume average of a certain local variable

Note that subscripts *f* and *m* refer to the fluid phase and solid matrix phase, respectively. The decomposition of the local variable

Moreover,

In order to simplify the foregoing set of the equations, the low Mach approximation is applied due to the relatively low Mach number when the air flows through a porous medium. Thus, the dynamic pressure change is sufficiently small as compared to the absolute pressure prevailing over the system, such that the stagnant enthalpy is approximated by

For the time being, let us assume

where

is the Darcian average of the variable

or

The term

such that

Using the effective porosity

for the solid matrix phase as:

Note that the assumption of equal temperature gradients,

while the interfacial heat transfer between the solid and fluid phases is modeled using Newton’s cooling law:

where

and

where the exponent n is 0.7 according to [4]. The specific heat capacity of the air

## 3. One-dimensional analysis for volumetric receiver

In this section, we perform one-dimensional analysis to obtain analytic solutions for convective-radiative heat transfer in volume receiver. Prior to that, the radiative heat flux

### 3.1 Analysis based on the Rosseland approximation

In the Rosseland approximation, the radiative heat flux is given by

where

As schematically shown in Figure 1, the air is flowing through a passage of length *L* at the rate of the mass flux

According to Calmidi and Mahajan [15, 16], Dukhan [17], Kuwahara et al. [18] and Yang et al. [19, 20], the permeability and inertial coefficient of foams are given by

and

respectively, where

With respect to the stagnant thermal conductivity and the volumetric heat transfer coefficient for foams, Calmidi and Mahajan [15, 16] empirically provided the following correlations:

Kamiuto et al. [21] experimentally affirmed that the Rosseland model is quite effective. Therefore, it can be deduced that the Rosseland model is also applicable for the present case of silicon carbide ceramic foam. Based on the measurements made on cordierite ceramic foams by Kamiuto et al., the mean extinction coefficient

For a given mass flux

such that

where

Likewise we shall define the solid phase average temperature as follows:

The two energy equations, that is, Eqs. (21) and (22) may be added together and integrated using the boundary conditions in Eqs. (29) and (31) to give

This equation is substituted into Eq. (21) to eliminate

This ordinary differential equation, with the boundary conditions in Eqs. (29), (30) and (31) and the auxiliary asymptotic condition

and

where

where

The solid phase temperature at the inlet

and

respectively. Usually, the receiver length *L* is sufficiently long to reach the local thermal equilibrium. Thus, the average air and solid temperatures are evaluated from

As one of the most important performance parameters, the receiver efficiency is defined by

Having established the temperature development, the momentum equation, that is, Eq. (20) along with the equation of state can easily be solved to find out the pressure distribution along the receiver as

Under the low Mach approximation, the required pumping power per unit frontal area may be evaluated from

Note that the dynamic pressure change is sufficiently small as compared to the absolute pressure such that

### 3.2 Analysis based on the P1 model

Since the Rosseland approximation used in the previous analysis ceases to be valid near boundaries, the validity of applying the Rosseland approximation near the inlet boundary of the receiver should be investigated. In order to examine the validity of the Rosseland approximation, the results based on the Rosseland approximation will be compared with the results obtained from solving the irradiation transport equation based on the P1 model. Since the silicon carbide ceramic foam is optically thick, the radiant energy emitted from other locations in the domain is quickly absorbed such that the radiative heat flux is given by

where the diffuse integrated intensity

where

Moreover, the effects of turbulence mixing on the heat transfer are also considered. Therefore, the energy equation for the air will be written as

where turbulent Prandtl number

Under the low Mach number approximation, namely, we may reduce the macroscopic governing equations namely Eqs. (1), (49), (14) and (48) to a one-dimensional set of the equations as follows:

The turbulence kinetic energy is dropped from the momentum equation since it stays nearly constant within the receiver.

Nakayama and Kuwahara [22] established the macroscopic two-equation turbulence model, which does not require any detailed morphological information for the structure. The model, for given permeability and Forchheimer coefficient, can be used for analyzing most complex turbulent flow situations in homogeneous porous media. For the case of fully developed turbulent flow in an isotropic porous structure, the eddy viscosity is given by

Note that

such that the dispersion thermal conductivity usually overwhelms the eddy thermal conductivity.

For absorption coefficient

The boundary conditions of

and

Furthermore, the streamwise gradients of the dependent variables *x* = *L*.

The two energy equations, namely, Eqs. (51) and (52) may be added together and integrated using the boundary conditions in Eqs. (29) and (57) to give

Eqs. (52) and (53) are combined to give

This equation, Eq. (59), and Eq. (51) are substituted into Eq. (58) to eliminate

This ordinary differential equation, with the boundary conditions in Eqs. (29), (56) and (57) and the zero derivative conditions far downstream (*L* is sufficiently large), yields Eqs. (36) and (37). Note that

where

The solid phase temperature at the inlet

where the boundary condition in Eq. (56) is utilized. Usually, the receiver length *L* is sufficiently long to reach the local thermal equilibrium. Thus, the average air and solid temperatures are evaluated from Eqs. (42) and (43).

## 4. Validations of the Rosseland approximation

Smirnova et al. [23] numerically studied the compressible fluid flow and heat transfer within the solar receiver with silicon carbide monolithic honeycombs. In their paper, the following input data were collected to obtain the analytic solutions based on the present local thermal non-equilibrium model:

However, it should be noticed that the porosity of the silicon carbide monolithic honeycombs is not available in Smirnova et al. [23], its value was estimated to be

As for possible instabilities, the fifth-order characteristic Eq. (61) based on the P1 model should be examined carefully. Figure 2 shows the residual of the fifth-order characteristic equation

The third-order characteristic Eq. (38) based on the Rosseland approximation, on the other hand, yields only one positive root

In Figure 4, the present analytic solutions are compared against the large-scale FEM numerical calculations based on COMSOL, reported by Smirnova et al. [23]. It should be mentioned that the direct numerical integrations of Eqs. (20)–(22) were also carried out using the finite volume method code, SAINTS [12]. As the convergence criteria, the residuals of all equations are less than 10^{−5}. It can be clearly seen that the air temperature increases as receiving heat from the monolithic receiver. Eventually, these two phases reach local thermal equilibrium near the exit. Both sets of solutions agree very well with each other, indicating the validity of the present local thermal non-equilibrium model.

## 5. Applications to silicon carbide ceramic foam volumetric receiver

In order to overcome the problems associated with thermal spots and flow instabilities, we would like to study fluid flow and heat transfer characteristics in silicon carbide ceramic foams based on the analytical expressions of pressure and temperature fields within a solar volumetric receiver. The performance of the receiver may be assessed in terms of the receiver efficiency

All other parameters are evaluated using Eqs. (17), (18) and from Eq. (23) to (28).

As shown in Figure 5, it is interesting to note that *PP*, which means that the pore diameter

As indicated in Eq. (46), it can be easily deduced that *PP* and *PP*, which results in that the amount of heat carried by the air, *PP* range, in which

or

since

Thus, Eq. (46) may be written for the case in which the sharp rise in

which, for given PP, gives the minimum value of the pore diameter

For *PP* = 300, 500 and 1000 W/m^{2} studied here, Eq. (69) gives

In order to achieve local thermal equilibrium for the two phases within the receiver, the length of the receiver is assumed to be sufficiently long in the present study. In view of minimizing the required pumping power, however, it is noticeable that shorter length is better, as clearly seen from Eq. (46). Hence, a minimum length required to approach local thermal equilibrium may be chosen to design a receiver, which would guarantee both maximum receiver efficiency and minimum pumping power. Therefore, we may roughly set the optimal receiver length as

such that

Eq. (71) together with Eq. (69) provides useful information for designing a volumetric solar receiver of silicon carbide ceramic foam.

## 6. Conclusions

For the first time, the complete set of analytical solutions, which fully considers the combined effects of turbulence, tortuosity, thermal dispersion, compressibility on the convective, conductive and radiative heat transfer within a ceramic foam receiver, is presented based on the two-energy equation model of porous media. Both the Rosseland approximation and the P1 model are applied to account for the radiative heat transfer through the solar receiver, while the low Mach approximation is exploited to investigate the compressible flow through the receiver. Based on the P1 model, two positive roots were found from the characteristic equations of the fifth-order differential equation, indicating possible occurrence of hydrodynamic and thermal instabilities. However, it has been found that the Rosseland approximation for this case, despite its failure near the inlet boundary, is fairly accurate and may well be used for quick estimations and further analysis. Due to their advantages, such as high thermal conductivity and fluid mixing, silicon carbide ceramic foams are considered as a possible candidate for the receiver, which can overcome the problems associated with thermal spots and flow instabilities. The results show that the pore diameter must be larger than its critical value to achieve high receiver efficiency. As a result, there exists an optimal pore diameter for achieving the maximum receiver efficiency under the equal pumping power. The optimal pore diameter yielding the maximum receiver efficiency may be found around the critical value given by Eq. (71). A simple relation is derived for determining the length of the volumetric solar receivers of silicon carbide ceramic foam.

## Conflict of interest

The authors declare no conflict of interest.

## Nomenclature

A | surface area (m2) |

Aint | interfacial surface area between the fluid and solid (m2) |

b | inertial coefficient (1/m) |

c | specific heat (J/kg K) |

cp | specific heat at constant pressure (J/kg K) |

dm | pore diameter (m) |

G | mass flux (kg/m2 s) |

h | specific enthalpy (J/kg) |

hv | volumetric heat transfer coefficient (W/m3K) |

I0 | intensity of radiation (W/m2) |

k | thermal conductivity (W/m K) |

K | permeability (m2) |

L | receiver length (m) |

nj | normal unit vector from the fluid side to solid side (−) |

PP | pumping power per unit frontal area (W/m2) |

Pr | Prandtl number (−) |

q | heat flux (W/m2) |

R | gas constant (J/kg K) |

T | temperature (K) |

ui | velocity vector (m/s) |

V | representative elementary volume (m3) |

xi | Cartesian coordinates (m) |

x | axial coordinate (m) |

β | mean extinction coefficient (1/m) |

γ | dimensionless parameter (−) |

ε | porosity (−) |

ε∗ | effective porosity (−) |

ξ | incidence angle (rad) |

η | receiver efficiency (−) |

λ | characteristic coefficient (1/m) |

μ | viscosity (Pa s) |

ν | kinematic viscosity (m2/s) |

ρ | density (kg/m2) |

σ | Stephan-Boltzmann constant (W/m2K4) |

κ | absorption coefficient (1/m) |

τij | stress tensor (Pa) |

Special symbols | |

φ˜ | deviation from intrinsic average |

ϕ | Darcian average |

ϕf,s | intrinsic average |

Subscripts and superscripts | |

dis | dispersion |

eq | equilibrium |

f | fluid |

s | solid |

stag | stagnation |

0 | reference |