Next Article in Journal
A Machine-Learning-Based Approach to Critical Geometrical Feature Identification and Segmentation in Additive Manufacturing
Previous Article in Journal
Application of Blockchain Technology in Environmental Health: Literature Review and Prospect of Visualization Based on CiteSpace
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solar Energy Management Using a V-Groove: An Approach Based on a Multiple Optical Path Algorithm

1
Laboratoire Énergétique Mécanique Électromagnétisme EA 4416, Université Paris Nanterre, 92000 Nanterre, France
2
ICube UMR-7357, CNRS, University of Strasbourg, INSA, 67400 Strasbourg-Illkirch, France
3
SDM Research Group, International University of Beirut BIU, Beirut 146404, Lebanon
4
School of Engineering, Lebanese International University, Bekaa 1803, Lebanon
*
Author to whom correspondence should be addressed.
Technologies 2022, 10(5), 101; https://doi.org/10.3390/technologies10050101
Submission received: 16 July 2022 / Revised: 17 August 2022 / Accepted: 2 September 2022 / Published: 12 September 2022

Abstract

:
Angular and spectral separations of thermal radiation have become a key challenge in solar concentration or thermal management of sources radiating at extremely high or low temperatures. Reflections obtained from electromagnetic theory in an open cavity geometry increase the emission and absorption compared to a flat surface due to the cavity effect. In this paper, in order to obtain the directional emission of geometric surfaces (V-Grooves) using ray tracing and studying the propagation of light, a new algorithm is developed. The numerical simulations take into account the materials properties of both facets of the V-shape, thus simulating an original asymmetric and a multilayer V-shape and providing a very interesting directive thermal emission behavior. We evaluated the emission behavior from the reflection and emission coefficients of different rays at different angles for different parameters (materials properties, wavelength, and geometry). The simulations of a V-groove showed that due to the different reflections occurring inside an aluminum V-cavity with an aperture angle of 28°, the emissivity was well enhanced by 86% in the normal direction compared to a flat surface made of the same material. Moreover, using the original asymmetric opposite-sided materials (Al and Si) in a V- groove, it was possible to separate and control the emission by focusing the radiation or directing different spectral bands in different directions.

1. Introduction

The thermal management of energy systems represents a fundamental scientific issue in the field of sustainable development [1], in particular to increase the eco-efficiency in the various fields of use such as aeronautics [1], solar energy [2], or building glazing. Controlling the radiative properties of the surface enables optimization of thermal management systems, and offers an important perspective for improving collection performance by enhancing the emission or the reflection of the surface over specific wavelength range. The thermal radiation depends mainly on the nature of the radiating surface. Actually, the emission of gray surfaces is evaluated by the energy emitted by a black body (Plank’s law) [3] corrected by its emissive factor ( ϵ ). Its emissivity depends on the wavelength ( λ ), temperature ( Τ ), angle of incidence ( θ ), and roughness value. For other rough surfaces due to surface slopes such as periodic surfaces (periodic V-grooves, semicircular gratings, photonic crystals, etc.) or multiple thin films [4], the emission geometry and scattering properties of the medium change and affect both the directional and spectral properties of emitted thermal radiation [3,4,5,6,7]. V-groove geometries are found in a variety of mechanical designs such as deployable mechanisms [8], solar absorbers [8,9], emitters [8], solar cells [9], and thermal radiators [8]. For these complex surfaces, the emissivity depends on the geometric properties of the surface. When radiant energy is emitted from the surface of a cavity or enters a cavity, the resulting multiple reflections increase the absorption and thus, according to the Kirchhoff’s law of thermal radiation, increase the emission of the open cavity compared to a flat surface due to the cavity effect [8]. Sparrow and Lin [10,11] initially studied the apparent radiative properties of V-grooved surface. Later, Kanayama [12] studied the apparent directional emittances of V and circular grooved surfaces. He experimentally validated the theoretical calculations by testing an aluminum V-groove for several geometric opening angles. Recent new studies published by Rydge et al. [8,13,14] have considered the total hemispherical apparent radiative properties of the infinite V-groove. These studies discuss the methodology for calculating the apparent emissivity for an isothermal cavity and present some numerical results of the apparent emissivity of a groove as a function of the intrinsic emissivity for different opening angles using a Monte Carlo ray tracing algorithm. The Monte Carlo ray tracing algorithm is a direct numerical method that can be used to quantify the apparent emissivity and apparent absorptivity of arbitrarily shaped cavities [3,15,16]. It is worth noted that surfaces are made directionally selective by altering mechanically their finishing on a microscopic (microgrooves) or macroscopic scale. For example, large V grooves-compared to the wavelengths of the radiation-tend to reflect the incoming radiation several times for near-normal incidence [17]. The number of reflections decreases with the increase of the incidence angle.
In the presented works, the apparent emissivity of V-groove shapes is evaluated using geometric methods that take into account the intrinsic emissivity of the radiating cavity surface, but without considering the directional and specular dependence of the reflections produced by the groove facets. The novelty of this paper lies in the new approach that allows to study new arrangements such as a V-shape with two facets made of different materials where each reflection or absorption by a facet of the groove depends on both the complex optical index of the material and the angle of incidence. A new model considering the calculations of the reflection and emission coefficient of a V-groove using ray tracing is presented.

2. Theoretical and Numerical Methods

A multi-path ray tracing method to calculate and find the apparent directional emissivity of a V-groove shape is developed. The calculations of the coefficients of reflection and specular absorption of the various rays as a function of angles of incidence are made using the Fresnel formulas between a homogeneous and isotropic medium [3]. In this study, only perfect smooth conductive or dielectric surfaces are considered, i.e., diffraction and diffusion are considered negligible. In this study, we were interested in determining the radiative behavior of V-shape cavity emitters whose facets are considered perfectly smooth producing a specular reflection. In this case, the reflection and absorption coefficients are described by the Fresnel equations deduced from Maxwell’s electromagnetic theory. The optical properties of a material are linked to its refractive index, real index if it is a dielectric without loss, complex index if it exhibits absorption or conduction [17]. The complex refractive index (Equation (1)) of a conducting material depends on the wavelength ( λ ) and the angle of incidence ( θ ):
n ˜ ( λ ) = n ( λ ) + i k ( λ )
All optical and radiative properties-reflectivity, absorptivity and emissivity (resp. ρ ,   α ,   ε )-can be calculated using Fresnel’s formulas, for a given wavelength, temperature and angle of incidence. We present in Figure 1 the variation of the real and imaginary parts of the complex refractive index of gold, platinum, silicon and aluminum [18,19]. In order to calculate the optical properties for a given temperature, it is necessary to know the evolution of the real and imaginary parts of the complex refractive index for this temperature.
The amount of light reflected at an interface depends on the polarization (perpendicular or parallel to the incident plane wave), the angle of incidence and the refractive indices of the two media. For the simulations, an interface between a perfect dielectric and an absorbing medium is considered. In this case, the plane wave is incident from a perfect dielectric (air) on an absorbing medium. By solving the following system of Equations (2)–(4), the two reflection coefficients of perpendicular and parallel polarizations [4] are calculated:
ρ ( λ , θ ) = ( n 1 cos θ p ) 2 + q 2 ( n 1 cos θ + p ) 2 + q 2 ρ ( λ , θ ) = ( p n 1 sin θ tan θ ) 2 + q 2 ( p + n 1 sin θ tan θ ) 2 + q 2 ρ ( λ , θ )
p 2 = 1 2 [ ( n 2 2 k 2 2 n 1 2 sin 2 θ ) 2 + 4 n 2 2 k 2 2 + ( n 2 2 k 2 2 n 1 2 sin 2 θ ) ]
q 2 = 1 2 [ ( n 2 2 k 2 2 n 1 2 sin 2 θ ) 2 + 4 n 2 2 k 2 2 ( n 2 2 k 2 2 n 1 2 sin 2 θ ) ]
ρ ( λ , θ )   and ρ ( λ , θ ) are the energy reflection coefficients in the polarized directions, parallel and perpendicular to the incident plane, respectively [4]. These two coefficients are defined as the square of the ratio of the parallel or perpendicular component of the electric field by the incident electric field. In this study, n 2   and k 2   are the refractive and extinction coefficients of the complex refractive index of the absorbing material n ˜ 2 ( λ ) and n 1 = 1   is the real refractive index of the incident medium (air). Since thermal radiation is an unpolarized wave, then the average unpolarized reflectivity ( ρ ) and emissivity ( ε ) are given by the following Equations (5) and (6):
ρ ( λ , θ ) = ρ ( λ , θ ) + ρ ( λ , θ ) 2
ε ( λ , θ ) = 1 ρ ( λ , θ )
Indeed, for thermal emission, emissivity is analogous to an unpolarized electromagnetic coefficient [3].
The radiative heat transfer analysis of a V-groove geometry (or any other cavity geometry) is much more complicated than that of a normal flat surface due to the multiple reflections that occur within the cavity walls. The sides of the V-groove are assumed as two flat surfaces forming an angle ( ϕ ). The parameters are depicted in Figure 2. The apparent radiative properties of a cavity are function of the intrinsic radiative properties of the cavity walls, the type of reflection (specular or diffuse), the wavelength ( λ ), the overall shape and geometry of the cavity ( h ,   t , ϕ ), and the cavity material (n, k) [17]. They are affected by the temperature (T) of the cavity walls [8]. In this study, materials are considered at room temperature.
The apparent directional emittance is calculated based on the following assumptions:
(1)
The dimensions of the groove surface are much larger than the wavelength of the radiation [12].
(2)
The inner surface of the groove is optically smooth, and its material has a specific property of the specular reflection.
Figure 3 presents the ray tracing diagram of a V-cavity showing the minimum distance dx at which the algorithm works, verifying the two important assumptions above. The simulations are performed in the solar spectral range [ 0.1   μ m < λ < 2   μ m ], for a maximum wavelength of   λ m a x = 2   μ m , dx min = 10   λ m a x .
Homogeneous, isotropic materials with a perfectly defined chemical composition are considered. The separation surface between the media is considered optically smooth. The energy propagation medium between the reflecting and collecting surfaces is not diffusive. These conditions are ideal; they allow us to determine the influence of the optical properties of complex and real index materials on the radiative properties. Several functions were implemented to determine the intersections, angles, and reflection coefficients of the reflected rays inside the V- cavity; the sides of the V-cavity are treated as two flat surfaces forming an angle ϕ . The sides of the cavity are discretized into different facets. Simulations of the emitted and reflected rays inside the groove are done by choosing a large width with respect to the wavelength (a minimum ratio d x / λ = 10 for the study case considering infrared wavelength of 10 µm).
The groove is illuminated by a wavelength λ , beams of rays are projected onto the cavity for different angles and in all directions (diffuse irradiation) [8]. Equation (8) explains how the emission components for rays of different reflection numbers were calculated using our bi-directional tracing algorithm. For each ray, the reflection coefficient was evaluated by considering the angle of incidence ( β ) for all reflections within the cavity, as shown in Figure 4. It is noted that for each number of reflections (the beam of rays that are emitted and having the same number of reflections inside the cavity), there is a different emission behavior. From the calculation of the maximum number of reflections inside the cavity, the variation of the maximum number of reflections as a function of the opening angle of the cavity and the incident rays on the cavity are studied. The maximum number of reflections (nmax) can also be calculated using the following equation [8]:
  n m a x = [ ( π γ ) ϕ + 1 2 ]
where γ   is the incident angle, and ϕ is the cavity opening angle. The maximum number of reflections should be rounded to the nearest integer value. If the initial evaluation is exactly a whole number, then the number of reflections must be increased by one. The apparent directional emittance of the radiation is calculated using Equation (8) [12]:
ε a β = ε β 1 + ( 1 ε β 1 ) ε β 2 + ( 1 ε β 1 ) ( 1 ε β 2 ) ε β 3 + + ( 1 ε β 1 ) ( 1 ε β 2 ) ( 1 ε β 3 ) ( 1 ε β i 1 ) ε β i
where:
i = 1 ,   2 ,   3 ,   , N   ( number   of   reflections ) ; β : is the angle of emission or reflection; ε β 1 : is the emission component (rays emitted without any reflections); ( 1 ε β 1 ) ε β 2 : is for one-reflection component (rays emitted with a single reflection inside the groove). ( 1 ε β 1 ) ( 1 ε β 2 ) ε β 3 : is for a two-reflections component (rays emitted with two successive reflections inside the groove); and ( 1 ε β 1 ) ( 1 ε β 2 ) ( 1 ε β 3 ) ( 1 ε β i 1 ) ε β i : is for (i − 1) component of the reflections (rays emitted with i − 1 reflections), where (i − 1) is the maximum number of reflections that occurred inside the groove at a specific opening angle.

3. Multi-Layer V-Groove Design

A multilayer groove is modelled using the algorithm described in [20] combined with our optical algorithm. Figure 5 depicts the design. It consists of 6 layers of thickness λ/4 ni; and we have chosen to simulate by taking the average of the reflection and transmission on wavelengths in the visible and the near infrared domain, separately.

4. Results and Discussions

In the following, we present and discuss various results obtained using the ray tracing method [21]. The calculations and simulations were performed in the spectral band of solar radiation ( 0.1   μ m < λ < 2   μ m ). The reflectivity of several materials suitable for the design of a V-groove solar concentrator was evaluated under normal incidence. The behaviors show their spectral dependence on different parts of the electromagnetic spectrum (visible to near infrared NIR). Figure 6 presents the evolution of the spectral reflectivity of several materials at normal incidence. Results show the responses of different materials on the solar spectrum. As can be seen, for instance, gold has a lower reflectivity in the ultraviolet and in a small part of the visible region (~0.38 of λ [ 0.2   μ m   t o   0.5   μ m ] ). However, it tends to have a high reflectivity in other regions. Both gold and aluminum have perfect reflectivity over the infrared part of the solar spectrum about 0.95 for λ [ 1   μ m   t o   2.2   μ m ] , while Si and Pt have about 0.3 and 0.7, respectively on the same spectrum.
The optical properties of these materials are well known. Since they have different behaviors in the solar spectrum, they are good candidates to simulate the angular and spectral change of the emission of a V-groove surface made of these materials. It is also very interesting to find the change in spectral and directional emission using V cavity formed by two of these materials. The reflection coefficients ρ ( θ ) as a function of the angle of incidence were calculated using Fresnel equations. Based on the spectral distribution of the different materials in the solar spectrum, we chose to divide the spectrum into two different parts [ 0.1   μ m < λ 1 < 0.5   μ m ]   &   [ 0.5   μ m < λ 2 < 2   μ m ] , as there is a noticeable difference in reflectivity between these two regions, as shown in Figure 5. The directional and spectral dependence are well known. For instance, the reflectivity of gold in the first spectral range ( λ 1 ) increases from 0.4 to one, when the angle of incidence increases from 0° to 90°, while in the second spectral range ( λ 2 ), gold shows high reflectivity close to unity for all angles of incidence. The large variation in reflection coefficients was recorded at grazing angles in the range [70 ° to 90 ° ] for all materials. Therefore, it was important to take this dependence and variation into account in the calculations of the different reflections in our V-groove model.

4.1. Influence of the Angle of Incidence and the Opening Angle

The study of the influence of the angle of incidence and the angle of the opening groove on the number of reflections have been performed. A parametric study on V-cavity has been done. At a first stage, the maximum number of reflections as a function of incident angles for different opening angles of the V-groove was calculated using Equation (7) and compared to our ray-tracing program. In this way, the maximum number of reflections that occur inside the cavity was recorded for different angles of the cavity. Figure 7 reveals the evolution of the total number of reflections calculated. It was noted that the largest number of reflections occurs for the smallest opening angle. Furthermore, for a specific aperture angle the number of reflections that occur inside the cavity gradually decreases as the angle of incidence increases. For example, at ϕ = 30 ° , a maximum of 6 reflections occur for rays of incidence angle 0 ° < γ < 25 ° , and a maximum of 5 reflections occur for rays of incidence angle 25 ° < γ < 55 ° , and when the incident angles increase to 55 ° < γ < 85 ° , a maximum of 4 reflections occur inside the cavity.
Figure 8 was obtained by recording the maximum total number of reflections that occur as a function of the opening angle, and for different emission angles inside the cavity. Obviously, it has been observed that as the opening angle increases, the number of reflections inside the groove gradually decreases. For instance, for a specific angle of 30°, a maximum of six reflections occurred within the groove, whereas when the angle increased to 100°, only two reflections were encountered.

4.2. Total Number of Reflections and Light-Matter Interaction

The components of the emission by considering the total number of reflections and the light-matter interaction are here presented. The total directional emissivity of the cavity depends on the number of reflections that occur within the groove. The emission components of the rays of different numbers of reflections have been calculated and show the influence of the number of reflections on the emission diagram for a typical V-groove. Figure 9 shows the emission components of an aluminum V-groove for different numbers of reflections. Data for all rays of a single reflection inside the groove are extracted and considered to compute the polar diagram of the thermal emission. The same is done for the rays present in a second reflection, third, and so on. Then, the data of all the different reflection beams are collected alone. Finally, all data collected for a small solid angle are summarized. The maximum value of the emissivity ( 0.19) occurs for rays of 6 reflections, and it decreases to a minimum value of 0.022 for rays with only one reflection which means that about 88.42% of the energy was absorbed by the cavity due to the 6 consecutive reflections.

4.3. Total Directional Apparent Emissivity

In the following, the total directional apparent emissivity is evaluated. The apparent directional emittance was calculated by summing the emission components of the different reflections as presented in Figure 10 for a small solid angle. In order to validate our method, we compared our results to those of Kanayama [12] and found a very good agreement. The directional variation of emission and reflection from different angles was accurately calculated and using ray tracing, all components of the reflection and emission were recorded as a function of angle of incidence. Finally, the emission components of the different reflections were summed in all directions to obtain the apparent emissivity of a groove with an opening angle . The emissivity was enhanced from 0.022 to 0.28 in the normal direction by using an aluminum V-groove surface with a 30° opening angle compared to a normal flat surface.

4.4. Influence of the Opening Angle on the Thermal Emission Diagram

The influence of the opening angle on the thermal emission diagram is here studied. Several simulations of the aluminum V-groove were performed at different opening angles. The results from three different angles are presented in Figure 11. It can be clearly seen that for the smallest angle between the sides of the groove, it tends to focus more energy in the normal direction.
Figure 12 shows the results of the average apparent emissivity or absorptivity ( ε a = α a ) of V-groove subjected to diffuse irradiation at different groove angles ( ϕ ), calculated for different intrinsic surface emissivities ( ϵ = 0.05 ,   0.2 ,   0.4 ,   0.6 ) [13], and compared to the results obtained by our ray tracing algorithm. The emissivity of each side surface of the groove is assumed to be constant in all directions.

4.5. Comparison of Spectral Reflectivity

A comparison of the spectral reflectivity between symmetric V-cavities of different aperture angles and a flat surface is done in this section. After validation of the model, different simulations on a V-cavity with different aperture angles were performed, in order to study the spectral distribution of the reflectivity of V-cavities at different aperture angles. Figure 13 depicts the spectral reflectivity results in the normal direction ( θ = 0 ° ), of V-shaped surfaces at different angles and a flat surface, for two different study materials (Au and Si), respectively. By comparing the results in Figure 12, we have highlighted that V-shaped surfaces always have a lower reflectivity than a flat surface, especially for small opening angles. In addition, for spectral regions of lower reflectivity, there is a greater decrease or decline in the reflective coefficient. For instance, the reflectivity of gold in the spectral band where λ < 0.5   μ m is 0.4 for a flat surface, whereas for a groove surface it gradually decreases from 0.33 to 0 as the angle of the groove decreases from 100° to 30°. While, in the other NIR spectral bands λ > 1   μ m , there is a slight decrease in the reflection coefficient, about 0.075 drop. This decrease is due to the consecutive reflections that occur inside the cavity, because for small angles of the groove, the number of consecutive reflections is high. Similarly, for spectral regions with low reflection coefficients, the absorption is higher than that of high reflection coefficients. Thus, for a higher number of consecutive reflections, there is a sharp increase in absorption, and this is the reason of a sharp drop observed in the reflection coefficient in these regions, in contrast to the other high reflection coefficient region.

4.6. Influence of Parameters on the Reflection and Directional Emission Diagram

The influence of spectral distribution and opening angle on the directional reflection and emission diagram of symmetrical V-shaped cavities has also studied. The results for Au and Si are presented in Figure 14 and Figure 15. Each Figure presents four polar diagrams; (a) and (b) are for the directional reflectivity in the two spectral regions [ 0.1   μ m < λ 1 < 0.5   μ m ]   and   [ 0.1   μ m < λ 2 < 2   μ m ] respectively, (c) and (d) are for directional emissivity in the two spectral regions   λ 1   and   λ 2 , respectively. Each material has different directional behavior for the different spectral bands. Comparing the results of the V-shaped surface and the flat surface, it can be seen that the use of a V- shaped surface leads to a decrease in directional reflectivity in all directions and in all spectral bands, especially in the normal direction with small aperture angles.

4.7. Design of Asymmetric Materials

Directional thermal emission can be controlled by using an asymmetric V-groove materials design. The structure is depicted in Figure 16. Our model was used to simulate the radiative behavior of an asymmetric sided V-groove cavity, composed of the good conductor (Au) and the dielectric (Si), and presented in Figure 16, in the spectral of wavelength range between 0.1 and 2 µm. These two materials are available in micro technology and could be easily implemented in more elaborate structures, in particular to define optimized assemblies to separate angularly the wavelengths, infrared and visible range for example.
For Figure 17, calculations were performed on two parts of the solar spectrum [ 0.1   μ m < λ 1 < 0.5   μ m ]   a n d   [ 0.5   μ m < λ 2 < 2   μ m ] , and for different cavity opening angles ϕ = [ 30 ° ,   45 ° ,   60 ° ] . As expected, they show an asymmetric emission in the second wavelength band   λ 2 , which is distributed in different directions depending on the opening angle. Compared to the directional emission distribution of the two spectral bands, we see that there is an angular separation for these bands. For instance, for an opening angle of 30°, there is 60% absorption and 40% reflection for the first spectral band λ 1 in the angular range θ = [ 90 ° 75 ° ] , while there is no absorption for the second spectral band   λ 2 , and about 100% of this band reflected in the same angular range. This means that this surface is able to absorb 60% of the first spectral band and reflect all the other bands in a specific angular direction. These combinations of material in a V-shaped surface are able to separate angularly two different spectral bands by absorption and reflection.

4.8. Multilayers on the Facets of a V-Groove

In order to well separate the solar spectrum, we used the model shown in Figure 5. We have alternated the stacking of layers with refractive indices of 2.36 and 1.38 on a glass substrate (n = 1.52) and we have calculated separately for the main regions of the solar spectrum, for an opening angle equal to 30°. Figure 18 reveals that the angular separation for the reflection between the visible range (λ between 0.5 and 0.75 µm) and the near infrared range (1.15 µm > λ > 0.75 µm) is excellent. Indeed, the reflection is about 90% in the direction for the visible range (between 75° and 90°) and about 90% in the direction between 0° and 25° for the near infrared range. It is then possible, for example, to redirect the solar energy to the right for photovoltaic cells (visible) and to the left for water heating (near infrared).

5. Conclusions

In this paper, we have discussed and presented a way to calculate the reflection and emission diagrams of V-grooved surfaces, using an algorithm dealing with the different reflections that occur in a V-groove for different angles, materials and at different wavelengths. A Matlab program has been developed for this purpose. It allows to model the emission of different rough surfaces (sloped surfaces) and complex periodic structured surfaces simulated for different dependent parameters (geometrical aperture, materials combinations with respect to the spectral bandwidth separation). It is very interesting to study the emission from these surfaces, since they could play an important role in thermal management by controlling the thermal emission in a directional or spectral. Simulations of a V-groove showed that due to the different reflections occurring inside an aluminum V-cavity with an aperture angle of 28°, the emissivity was well enhanced by 86% in the normal direction compared to a flat surface made of the same material. In addition, it was demonstrated that by using original asymmetric opposite-sided materials (Al and Si) in a V- groove, it was possible to separate and control the emission by focusing the radiation or directing different spectral bands in different directions depending on the opening angle of the cavity.
Moreover, with the algorithm, we are able to simulate different wavelength ranges (Visible, and near infrared), and different complex surface designs (asymmetric sides, circular cavities, and multiple thin films), which will be crucial to find, compare, and then design spectrally selective surfaces. That is useful for directional and spectral control of thermal emission, and for improving the collection performance of solar electromagnetic radiation by separating different parts of the solar spectrum in different direction. The separation of the solar spectrum, using multilayer stacks on the both sides of the V-groove, is evidenced in this paper.

Author Contributions

Conceptualization, F.K., B.S., G.G. and M.H.; methodology, F.K., B.S. and G.G..; validation, F.K., B.S., G.G., M.H. and P.T.; formal analysis, F.K. and G.G.; writing—original draft preparation, F.K.; writing—review and editing, B.S., G.G., M.H. and P.T.; supervision, B.S., G.G. and M.H.; All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hill, S.A.; Kostyk, C.; Motil, B.; Notardonato, W.; Rickman, S.; Swanson, T. Draft Thermal Management Systems Roadmap: Technology Area 14. Available online: https://www.nasa.gov/pdf/501320main_TA14-Thermal-DRAFT-Nov2010-A.pdf (accessed on 1 February 2020).
  2. Alonso-Álvarez, D.; Augusto, A.; Pearce, P.; Llin, L.F.; Mellor, A.; Bowden, S.; Ekins-Daukes, N. Thermal emissivity of silicon heterojunction solar cells. Sol. Energy Mater. Sol. Cells 2019, 201, 110051. [Google Scholar] [CrossRef]
  3. Modest, M.F. Radiative Heat Transfer, 2nd ed.; Academic Press: Cambridge, MA, USA, 2003. [Google Scholar]
  4. Thakur, A.; Raman, R. Thermal emissivities of films on substrates. Appl. Energy 1983, 15, 1–13. [Google Scholar] [CrossRef]
  5. Prokhorov, A.V.; Hanssen, L.M.; Mekhontsev, S.N. Radiation Properties of IR Calibrators with V-grooved Surfaces. In Proceedings of SPIE, Proceedings of the Thermosense XXVIII, Kissimmee, FL, USA, 17–20 April 2006; SPIE: Bellingham, DC, USA, 2006; Volume 6205. [Google Scholar]
  6. Sai, H.; Yugami, H.; Akiyama, Y.; Kanamori, Y.; Hane, K. Spectral control of thermal emission by periodic microstructured surfaces in the near-infrared region. J. Opt. Soc. Am. A Opt. Image Sci. Vis. 2001, 18, 1471. [Google Scholar] [CrossRef] [PubMed]
  7. Hollands, K.G.T. Directional selectivity, emittance, and absorptance properties of vee corrugated specular surfaces. Sol. Energy 1963, 7, 108–116. [Google Scholar] [CrossRef]
  8. Mulford, R.B.; Collins, N.S.; Farnsworth, M.S.; Jones, M.R.; Iverson, B.D. Total hemispherical apparent radiative properties of the infinite V-groove with specular reflection. Int. J. Heat Mass Transf. 2018, 124, 168–176. [Google Scholar] [CrossRef]
  9. Bailey, S.; Fatemi, N.; Landis, G.A.; Brinker, D.; Faur, M. Application of v-groove technology to InP solar cells. In Proceedings of the International Conference on Indium Phosphide and Related Materials, Denver, CO, USA, 23–25 April 1990; pp. 73–79. [Google Scholar]
  10. Sparrow, E.M.; Lin, S.H. Absorption of thermal radiation in a v-groove cavity. Int. J. Heat Mass Transf. 1962, 5, 1111–1115. [Google Scholar] [CrossRef]
  11. Sparrow, E.M.; Lin, S.H. Radiation heat transfer at a surface having both specular and diffuse reflectance components. Int. J. Heat Mass Transf. 1965, 8, 769–779. [Google Scholar] [CrossRef]
  12. Kanayama, K. Apparent Directional Emittances of the Random Rough Surfaces of Metals and Nonmetals. Trans. Jpn. Soc. Mech. Eng. 1973, 39, 1635–1643. [Google Scholar] [CrossRef]
  13. Mulford, R.B.; Collins, N.S.; Farnsworth, M.S.; Jones, M.R.; Iverson, B.D. Total hemispherical apparent radiative properties of the infinite V-groove with specular reflection. J. Thermophys. Heat Trans. 2018, 32, 1109–1112. [Google Scholar] [CrossRef]
  14. Prokhorov, A.V. Monte Carlo method in optical radiometry. Metrologia 1998, 35, 465–471. [Google Scholar] [CrossRef] [Green Version]
  15. Sapritsky, V.I.; Prokhorov, A.V. Calculation of the effective emissivities of specular-diffuse cavities by the Monte Carlo method. Metrologia 1992, 29, 9–14. [Google Scholar] [CrossRef]
  16. Prokhorov, A.V.; Hanssen, L.M.; Mekhontsev, S.N. Calculation of the radiation characteristics of blackbody radiation sources. Exp. Methods Phys. Sci. 2009, 42, 181–240. [Google Scholar]
  17. Raman, R.; Thakur, A. Thermal emissivity of materials. Appl. Energy 1982, 12, 205–220. [Google Scholar] [CrossRef]
  18. Polyanskiy, M.N. Refractive Index Database. Available online: https://refractiveindex.info (accessed on 14 January 2020).
  19. Palik, E.D. Handbook of Optical Constants of Solids; Academic Press: Orlando, FL, USA, 1985. [Google Scholar]
  20. Petit, R. Ondes Électromagnétiques En Radioélectricité Et En Optique; Masson: Paris, France, 1989; Volume 2. [Google Scholar]
  21. Kawtharani, F. Spectral and Angular Study of Thermal Emission of Cavities: Prospective for Solar Energy & Thermal Management. Ph.D. Thesis, University of Paris-Nanterre, Nanterre, France, July 2022. [Google Scholar]
Figure 1. Refractive and extinction coefficients of some materials as function of wavelength.
Figure 1. Refractive and extinction coefficients of some materials as function of wavelength.
Technologies 10 00101 g001
Figure 2. V groove parameters.
Figure 2. V groove parameters.
Technologies 10 00101 g002
Figure 3. Ray tracing diagram showing the minimum width in the V-cavity.
Figure 3. Ray tracing diagram showing the minimum width in the V-cavity.
Technologies 10 00101 g003
Figure 4. Geometric representation of the V-groove simulation.
Figure 4. Geometric representation of the V-groove simulation.
Technologies 10 00101 g004
Figure 5. Representation of a V-shaped cavity whose sides are made of thin layer materials (the figure on the right corresponds to a stack of 6 layers).
Figure 5. Representation of a V-shaped cavity whose sides are made of thin layer materials (the figure on the right corresponds to a stack of 6 layers).
Technologies 10 00101 g005
Figure 6. Unpolarized reflectivity versus wavelength for several conductors and dielectric materials calculated for normal incidence.
Figure 6. Unpolarized reflectivity versus wavelength for several conductors and dielectric materials calculated for normal incidence.
Technologies 10 00101 g006
Figure 7. Total number of reflections occurring inside the cavity as a function of the angle of incidence for different opening angles.
Figure 7. Total number of reflections occurring inside the cavity as a function of the angle of incidence for different opening angles.
Technologies 10 00101 g007
Figure 8. Maximum number of reflections as a function of the opening angle of the V-groove ϕ as computed from Equation (7) (—) and our ray tracing algorithm (●).
Figure 8. Maximum number of reflections as a function of the opening angle of the V-groove ϕ as computed from Equation (7) (—) and our ray tracing algorithm (●).
Technologies 10 00101 g008
Figure 9. Polar diagrams of the NIR emissivity computed using ray tracing of each of the five emission components through the sides of the groove composed of aluminum for different numbers of reflections at an opening angle 30°.
Figure 9. Polar diagrams of the NIR emissivity computed using ray tracing of each of the five emission components through the sides of the groove composed of aluminum for different numbers of reflections at an opening angle 30°.
Technologies 10 00101 g009
Figure 10. Total apparent directional emissivity of a V-groove and a flat aluminum surface in the infrared range compared to [12].
Figure 10. Total apparent directional emissivity of a V-groove and a flat aluminum surface in the infrared range compared to [12].
Technologies 10 00101 g010
Figure 11. Apparent directional emissivity of an aluminum V-groove at different opening angles compared to the calculated data [12].
Figure 11. Apparent directional emissivity of an aluminum V-groove at different opening angles compared to the calculated data [12].
Technologies 10 00101 g011
Figure 12. Average apparent emissivity or absorptivity calculated for diffuse irradiation [8].
Figure 12. Average apparent emissivity or absorptivity calculated for diffuse irradiation [8].
Technologies 10 00101 g012
Figure 13. Spectral reflectivity of V-cavity at different aperture angles ( ϕ = 30 ° ,   45 ° ,   60 ° ,   80 ° ,   100 ° ) and flat surface, formed for the two chosen materials (Au and Si) in the normal direction ( θ = 0 ° ).
Figure 13. Spectral reflectivity of V-cavity at different aperture angles ( ϕ = 30 ° ,   45 ° ,   60 ° ,   80 ° ,   100 ° ) and flat surface, formed for the two chosen materials (Au and Si) in the normal direction ( θ = 0 ° ).
Technologies 10 00101 g013
Figure 14. (a,b): Directional reflectivity; (c,d): apparent emissivity of a V-cavity and flat surface formed of gold in two different spectral ranges [ λ 1 , λ 2 ], at different opening angles ϕ = [30°, 45°, 60°].
Figure 14. (a,b): Directional reflectivity; (c,d): apparent emissivity of a V-cavity and flat surface formed of gold in two different spectral ranges [ λ 1 , λ 2 ], at different opening angles ϕ = [30°, 45°, 60°].
Technologies 10 00101 g014
Figure 15. (a,b): Directional reflectivity; (c,d) apparent emissivity of a V-cavity and flat surface formed of silicon in two different spectral ranges [ λ 1 , λ 2 ], at different opening angles ϕ = [30°, 45°, 60°].
Figure 15. (a,b): Directional reflectivity; (c,d) apparent emissivity of a V-cavity and flat surface formed of silicon in two different spectral ranges [ λ 1 , λ 2 ], at different opening angles ϕ = [30°, 45°, 60°].
Technologies 10 00101 g015
Figure 16. V cavity of asymmetric sides formed of Au & Si.
Figure 16. V cavity of asymmetric sides formed of Au & Si.
Technologies 10 00101 g016
Figure 17. Polar diagrams for the distribution of directional emission and reflection from the sides of a V-cavity composed of two different materials (Au and Si).
Figure 17. Polar diagrams for the distribution of directional emission and reflection from the sides of a V-cavity composed of two different materials (Au and Si).
Technologies 10 00101 g017
Figure 18. Directional reflection and transmission for three different spectral regions for five double multilayers (10 Layers). ϕ = 30 ° .
Figure 18. Directional reflection and transmission for three different spectral regions for five double multilayers (10 Layers). ϕ = 30 ° .
Technologies 10 00101 g018
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kawtharani, F.; Serio, B.; Guida, G.; Twardowski, P.; Hammoud, M. Solar Energy Management Using a V-Groove: An Approach Based on a Multiple Optical Path Algorithm. Technologies 2022, 10, 101. https://doi.org/10.3390/technologies10050101

AMA Style

Kawtharani F, Serio B, Guida G, Twardowski P, Hammoud M. Solar Energy Management Using a V-Groove: An Approach Based on a Multiple Optical Path Algorithm. Technologies. 2022; 10(5):101. https://doi.org/10.3390/technologies10050101

Chicago/Turabian Style

Kawtharani, Fadel, Bruno Serio, Geraldine Guida, Patrice Twardowski, and Mohammad Hammoud. 2022. "Solar Energy Management Using a V-Groove: An Approach Based on a Multiple Optical Path Algorithm" Technologies 10, no. 5: 101. https://doi.org/10.3390/technologies10050101

APA Style

Kawtharani, F., Serio, B., Guida, G., Twardowski, P., & Hammoud, M. (2022). Solar Energy Management Using a V-Groove: An Approach Based on a Multiple Optical Path Algorithm. Technologies, 10(5), 101. https://doi.org/10.3390/technologies10050101

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop