Abstract

This study provides information to design heterogeneous photocatalytic solar reactors with flat plate geometry used in treatment of effluents and conversion of biomass to hydrogen. The concept of boundary layer of photon absorption taking into account the efficient absorption of radiant energy was introduced; this concept can be understood as the reactor thickness measured from the irradiated surface where 99% of total energy is absorbed. Its thickness and the volumetric rate of photons absorption (VRPA) were used as design parameters to determine (i) reactor thickness, (ii) maximum absorbed radiant energy, and (iii) the optimal catalyst concentration. Six different commercial brands of titanium dioxide were studied: Evonik-Degussa P-25, Aldrich, Merck, Hombikat, Fluka, and Fisher. The local volumetric rate of photon absorption (LVRPA) inside the reactor was described using six-flux absorption-scattering model (SFM) applied to solar radiation. The radiation field and the boundary layer thickness of photon absorption were simulated with absorption and dispersion effects of catalysts in water at different catalyst loadings. The relationship between catalyst loading and reactor thickness that maximizes the absorption of radiant energy was obtained for each catalyst by apparent optical thickness. The optimum concentration of photocatalyst Degussa P-25 was 0.2 g/l in 0.86 cm of thickness, and for photocatalyst Aldrich it was 0.3 g/l in 0.80 cm of thickness.

1. Introduction

Heterogeneous photocatalysis based on TiO2 and modified photocatalysts is widely used in energetic and environmental applications, including water and air purification systems [1, 2], self-cleaning surfaces [3], wettability patterns, and printing applications [4], and as a clean production route of hydrogen fuel [5]. The first step of photocatalytic reaction mechanism [6] is the radiant energy absorption higher than the band gap of the semiconductor to generate electron-hole pairs on its surface, which can initiate oxidation and reduction reactions.

A clean route for hydrogen production has been proposed by photogenerated electrons at the conduction band from reduction of water in the absence of oxygen and oxidation of waste biomass using solar energy [7].

Regarding the mineralization of organic contaminants, Turchi and Ollis proposed a path based on the generation of oxidizing species (Table 1) [6]; adsorbed water on the TiO2’s surface (reactions 2a and 2b) reacts with holes to generate hydroxyl radicals (reactions 6a and 6b). The attack of the hydroxyl radical (reactions 9–12) on a C–H bond is the prevailing step for the oxidative breakdown, which is responsible for the mineralization of a wide range of pollutants (e.g., cyanide, arsenic, halogenated hydrocarbons, pesticides, endocrine disruptors, estrogens, inorganics, aromatics, emerging pollutants, and many others) [814]. More details of the mechanism and their considerations are in [15].

The above mentioned applications present a trend towards the use of solar energy as a sustainable source of energy, development of efficient reactors geometries, and scale-up of reactors. Therefore, rigorous methodologies have been developed for design and optimization of photoreactors based on scientific knowledge of the photocatalytic process [8, 9, 14, 1618], allowing a more economical and technically feasible process.

The rigorous description of photoreactor requires (i) radiation field, (ii) hydrodynamic field, and (iii) mass balance with kinetic expression whose parameters should be independent of photon absorption. The study of the radiation field provides optimal design information such as reactor geometry, catalyst selection, catalyst loading, and dimensions of the reactor [19, 20].

The variable most employed to quantify the spatial distribution of radiant energy absorbed within the reactor is the local volumetric rate of photon absorption (LVRPA); it depends on the geometry, radiation source, loading, and type of photocatalyst, in some cases on the pollutant if this presents absorption of radiant energy. A rigorous approach to determine the LVRPA is the solution of the radiative transfer equation (RTE), by the discrete ordinate method (DOM), which requires spatial, directional, and spectral discretizations inside the reactor [8, 20, 21].

A useful approach is the six-flux absorption-scattering model (SFM); it has been applied to different geometries: compound parabolic collector (CPC) photoreactors [10, 19], annular reactors [12, 13, 2224], and flat plate reactors [25]. This model calculates the LVRPA through a system of algebraic equations, where the scattering event can occur only along the six directions of the Cartesian coordinates [25]; it reduces the mathematical complexity to an analytical system, thus saving computational time and effort.

Flat plate reactors are scalable, and these can be used with solar radiation, so they are very attractive and also provide an excellent configuration for efficient excitation of the semiconductor photocatalyst TiO2 [26]; their modeling requires a complex analysis of the radiation field inside the photoreactor [27]. In water treatment, for a dilute pollutant is considered the photon absorption only by the solid photocatalyst particles [28, 29], therefore, the pollutant concentration and the radiation field are decoupled [20, 25] allowing us to separately study the effect of photon absorption of the water-catalyst solution.

In this study, design parameters of a flat plate photoreactor with photocatalysts based on TiO2 and solar radiation were analyzed in terms of absorption of photons. For mathematical convenience, the variable was employed; it is estimated by SFM model and does not require emission model. Also, a new design parameter for selection of the optimum reactor thickness called boundary layer thickness of photon absorption is presented. The information presented might be used in kinetic models of pollutants degradation and hydrogen production.

2. Methodology

2.1. Model of Boundary Layer of Photon Absorption

This paper proposes a new concept for heterogeneous photocatalytic reactor design, which allows us to determine the best thickness. Figure 1 shows a schematic representation of a double flat plate reactor with thick , located perpendicular to the source. Additionally, the figure presents a profile of LVRPA along the -coordinate, with maximum value at the interface surface-flow () and minimum value at reactor bottom ().

The region where there is a gradient of energy absorption has been called “boundary layer of photon absorption,” and its thickness is defined as -value which satisfies the following ratio: This definition of boundary layer thickness for transport phenomena and absorption of radiant energy is analogous to the definition of thickness of hydrodynamic, thermal, and concentration boundary layer on a flat plate [30]. The can be understood as the reactor thickness measured from the irradiated surface where 99% of total energy is absorbed. For this reason, boundary layer thickness becomes a design parameter for selection of optimum thickness of flat plate reactors. Oversizing of reactor thickness can lead to the presence of a dark sublayer, where there is no considerable photon absorption.

This design parameter is similar to the apparent optical thickness used for sizing of annular reactors and CPCs [19, 31]. The uses the analogous concept of boundary layer thickness to obtain the optimal light absorption. depends of catalyst (type and loading) and it can be used as a geometrical parameter equal to the optimum reactor thickness (see Figure 1).

2.2. Volumetric Rate of Photon Absorption VRPA

VRPA is defined as an average value of LVRPA in the whole volume. For a flat plate reactor the VRPA is expressed as [20] VRPA is a design parameter established and validated to determine catalyst concentrations in photocatalytic processes [20].

VRPA and boundary layer thickness presented above are design parameters of photoreactors and their calculations require a model of radiant field, which quantify LVRPA within the reactor.

2.3. Radiant Field Model

The radiation field of a reactor with incident solar radiation was modeled by six-flux absorption-scattering model (SFM) [10, 12, 14, 19, 24, 25]. LVRPA divided by flux of incident radiant energy in reactor expressed by this model is where is the coordinate of photon, in a flat plate reactor, is the corrected scattering albedo, corresponds to the corrected value of mean free path of photons in the suspension, and is a SFM parameter.

The parameters of (3) are function of average values of the scattering and absorption coefficients of catalyst, which are calculated as well: where and correspond to 280 y 395 nm, respectively, given by the interception of solar emission spectrum and TiO2 absorption spectrum; is the spectrum of solar emission power [32]; and are the spectrums of scattering and absorption coefficients, respectively, per mass unit of catalyst, reported in [33] using a phase function with diffuse reflectance for different commercial brands of TiO2.

Extinction coefficient is the sum of the absorption and scattering average coefficients . This parameter is used to calculate the scattering albedo which determines the fraction of dispersed energy.

Scattering albedo requires a correction to the model of six directions [25], through probabilities of forward, backward, and sideways scattering (, , ): where , , values are 0.11, 0.71, and 0.045 for a phase function of great sphere with diffuse reflectance [10]. Mean free path of photons in the suspension and its corrected value , which depend on catalyst concentration, are calculated as well [10]: Finally, the parameter is calculated through the following equation: where is apparent optimal thickness defined as and is the optical thickness, for a flat plate reactor [25, 26].

To obtain the new parameter, boundary layer thickness of photon absorption is proposed to minimize an objective function (10), which is obtained from the definition (1), using the Nelder-Mead optimization algorithm: where and correspond to the SFM evaluated in equal to zero and , respectively.

3. Results and Discussion

3.1. Optical Properties under Solar Radiation

The SFM parameters that are independent of the catalyst concentration are shown in Table 2 for six different commercial brands of TiO2 and solar radiation using (4)–(6). Parameters of photocatalytic activity are different due to the percentage of anatase, rutile, and brookite present in each trademark.

The values reported here disagree with those reported by [10] for Degussa under the same radiation condition ( and ) due to the use of different phase functions. Colina-Márquez et al. [10] applied optical properties described by Henyey and Greenstein phase function [34]; in contrast, this study employed optical properties of catalyst and scattering probabilities , , from a phase function model of diffuse reflectance [22, 25].

3.2. LVRPA

Figure 2 shows the profiles (3) as a function of -coordinate inside a flat plate photoreactor with solar radiation for catalyst Degussa and Aldrich for different catalyst concentrations. The optical properties presented in Table 2 were used. A reactor thickness equal to 1 cm was considered, because it is enough for photon absorption [20, 35]. Therefore, the dimensionless optical thickness varies from 4 to 22 depending on loading and type of catalyst.

The at irradiated surface of the reactor increases with catalyst concentration due to the amount of surface-exposed catalyst and back-scattering energy absorption from internal layer inside the reactor. These results agree with the literature [20]. In addition, it was found that the distribution of LVPRA is similar when the apparent optical thickness is used in annular and CPCs reactors [31].

The change of with respect to near the irradiated surface is associated with extinction coefficients; the total energy available in the system is quickly extinguished by increases in catalyst concentration and high extinction coefficients .

Finally, the distance from the surface () where the radiant energy is absorbed (i.e., boundary layer thickness) depends on concentration and type of catalyst. For catalyst loading greater than 0.3 g/l, the at the bottom of the reactor is negligible.

To validate the SFM in a flat plate reactor, the model was run against rigorous solution of RTE proposed by [9, 36]. The reported reactor had 1.2 cm of thickness, UV lamps as radiation source, Aldrich and Degussa P-25 as photocatalysts. The SFM parameters were calculated to the experimental conditions reported ( Einstein cm−2 s−1, catalyst Aldrich,  cm2 g−1, and  cm2 g−1) with a local percent error of LVRPA lower than 12, and 2.4 for VRPA.

3.3. Design of Reactors

Figure 3 shows the profiles (2) in solid lines as a function of catalyst loading inside a flat plate reactor, for six different commercial brands of TiO2 catalysts. The dashed lines represent the values of that satisfy (1).

For all catalysts, a catalyst loading less than 0.2 g/l has low rate of photon absorption; from 0.2 to 0.4 g/l there are high values of and thus the best operating conditions of the system are within this interval; loadings greater than 0.4 g/l reach some form of saturation; therefore, there is no justification for extra consumption of catalyst.

From profile in Figure 3 and Table 2, it is established that the absorption of radiant energy is faster for catalysts with a low coefficient of scattering albedo . In other words, high absorption coefficients and low scattering coefficients have greater absorption of photons, increasing the mathematical relationship . As a result, the absorption presents the following order: Fluka > Fischer > Merck > Aldrich > Degussa P-25 > Hombikat.

The catalyst Aldrich has a maximum VRPA 7% higher than the one presented by the catalyst Degussa. A similar trend was reported for polychromatic radiation using UV lamps, where the catalyst Aldrich is 19% more efficient than Degussa P-25 [36].

For all catalysts and loadings lower than 0.1 g/l, the absorption of photons occurs in the entire thickness of the reactor (1 cm), but the VRPA has small values. Catalyst loadings between 0.1 and 0.2 g/l employ the entire thickness of the reactor with approximately zero values of the LVRPA at the bottom of the reactor and begin to show a boundary layer of photon absorption thinner than 1 cm, due to some form of saturation between photon transport and the catalyst amount to absorb them.

For catalyst loading higher than 0.2 g/l, different boundary layer thicknesses of absorption can be seen. Catalysts with higher coefficient of photon extinction have a thinner boundary layer because the energy is quickly extinguished as it travels through the fluid.

However, Hombikat catalyst has the largest boundary layer of photon absorption, despite having an intermediate value of extinction coefficient; this is due to its low absorption coefficient and its high coefficient of scattering albedo ; most of the photons are dispersed in a large layer of fluid phase with a photon absorption rate slower than other catalysts.

Profile of boundary layer thickness of photon absorption represents the best conditions of design. A point above this curve presents a dark layer inside reactor, while a point below it presents VRPA less than the maximum possible.

3.4. Optimal Configuration of Reactors

Figure 3 provides useful information to define the design thickness of a flat plate reactor and type and loading of catalyst. For applications in very thin reactors (e.g., thin films of continuous fluid by sliding on an inclined flat plate), the catalyst with the best performance is Degussa P-25 followed by Aldrich; the catalyst selection is based on reactor thickness. On the other hand, for applications in very thick reactors, the catalyst with the best performance is Fluka followed by Fischer; the catalyst selection is based on the greater amount of energy absorbed.

To design a photoreactor using a specific catalyst, for example, Degussa P-25 (point line), the profile in Figure 3 recommends 0.2 g/l of catalyst to take the maximum energy absorbed with the minimum consumption of catalyst. The profile in the same figure suggests a design thickness of reactor equal to or slightly higher than 0.86 cm and take advantage of the entire volume within the reactor for the photon absorption. Table 3 was reported using the procedure above.

The optimum catalyst loading was calculated for catalysts water system; this increases in systems where there is absorption of radiant energy by substrate. Other authors have reported optimum values of catalyst loading similar to those reported in this paper. Brandi et al. report a loading close to 0.3 g/l based on the simulation of photon absorption for Aldrich and Degussa catalysts with radiation from UV lamps in a flat plate reactor [20, 36].

Experimental researches developed at the solar platform from Almeria obtained an optimal loading of Degussa P-25 in a plate flat solar reactor equal to 0.2 g/l [3739]; this concentration is the same obtained by simulation in this work.

Different geometries are not comparable; however, the reported values of optimum catalyst loading have proved to be very similar in different reactors. An example of this result is the work made by Colina-Márquez and coworkers. They found optimum catalyst loadings from 0.17 to 0.4 g/l for catalysts with scattering albedo between 0.75 and 0.95 in a CPC solar reactor and catalyst loadings between 0.2 and 0.6 g/l for a tubular reactor [19]. Particularly, they obtained 0.21 g/l as optimal loading for Degussa P-25.

For thicknesses not near 1 cm, it is recommended to use the design parameter called apparent optical thickness , where ; the parameter was proposed by [19] and modified here for its application to flat plate reactor. At conditions reported in Table 2, the values obtained were 4.39, 4.08, 4.15, 3.77, 3.99, and 4.10 to Aldrich, Degussa, Merck, Hombikat, Fisher, and Fluka, respectively. The best design conditions for different catalyst loadings and reactor thicknesses are obtained by keeping constant.

Reactor design information presented in this paper is useful to determine the optimum catalyst loading, the optimal thickness of reactor, and the radiant energy absorbed. In addition, it eliminates the need of statistical analysis of experiment design which involves considerable consumption of time and resources [19]. New photocatalysts can be analyzed using the methodology presented here by varying absorption and scattering coefficients.

4. Conclusions

This study has presented a new parameter for selection of thickness of a flat plate reactor called boundary layer thickness of photon absorption , which represents an effective thickness of absorption. Its use allows the elimination of dark zones within the reactor and absorption of the largest amount of photons. Also, together with the VRPA, it provides sufficient criteria for sizing of reactor thickness, the choice of catalyst, and its loading. Its use is applicable to the scale-up of solar reactors.

The parameter was evaluated for TiO2 semiconductors Degussa P-25, Aldrich, Merck, Hombikat, Fluka, and Fisher, as a function of reactor thickness, catalyst loading, and the maximum amount of energy absorbed getting optimal operating conditions in the catalyst loading.

Nomenclature

SFM parameter, dimensionless
:SFM parameter, dimensionless
:Catalyst concentration,
:Objective function,
:Thickness, cm
:Emission spectrum of radiation source, Einstein
:Flux of radiation at surface, Einstein
Local volumetric rate of photon absorption, Einstein
:Probability of backward scattering, dimensionless
:Probability of forward scattering, dimensionless
:Probability of side scattering, dimensionless
:Auxiliary coordinate in the photon flux direction, cm
:Cartesian coordinate, cm
:Cartesian coordinate, cm
:Cartesian coordinate, cm.
Greek Letters
:Specific extinction coefficient,
:SFM parameter, dimensionless
:Boundary layer thickness of photon absorption, cm
:Specific absorption coefficient,
:Radiation wavelength, nm
:Mean free path of photons, cm
:Extinction length, cm
:Specific scattering coefficient,
:Optical thickness, dimensionless
:Apparent optical thickness, dimensionless
:Scattering albedo, dimensionless
:Scattering albedo corrected, adimensional.
Subscripts
corr: Corrected
max: Maximum
min: Minimum
: Reactor
: Irradiated surface
: Bottom of reactor.
:Spectral.
Superscripts
:Specific.
Special Symbols
:Average value.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors are grateful to Universidad del Valle for the financial support to produce this work and to Colciencias for the financial support during their Ph.D. studies. M. A. Mueses thanks Universidad de Cartagena for the financial support. Finally, the authors thank Fernando Otálvaro for the English revision.