A Comparison Between the Monte Carlo Ray Trace and the FLUENT Discrete Ordinates Methods for Treating Solar Input to a Particle Receiver

Author(s):  
Ahmet Murat Mecit ◽  
Fletcher Miller

A new type of high temperature solar receiver for Brayton Cycle power towers is being designed and built in the Combustion and Solar Energy Laboratory at San Diego State University under a DOE Sunshot Award. The Small Particle Solar Receiver is a pressurized vessel with a window to admit concentrated solar radiation that utilizes a gas-particle suspension for absorption and heat transfer. As the particles absorb the radiation that enters the receiver through the window, the carrier fluid (air in this case) heats which oxidizes the particles and the flow leaves the receiver as a clear gas stream. After passing through an in-line combustor if needed, this hot gas is used to power a turbine to generate electricity. The numerical modelling of the receiver is broken into three main pieces: Monte Carlo Ray Trace (MCRT) method (written in FORTRAN), ANSYS Fluent (CFD), and the User Defined function (written in C code) for oxidation. Each piece has its advantages, disadvantages, and limitations and the three pieces are coupled to finalize the calculation. While we have successfully demonstrated this approach to obtaining the velocity and temperature fields, one big challenge to this method is that the definition of the geometry is a time consuming programming task when using MCRT. On the other hand, arbitrary geometries can be easily modelled by Computational Fluid Dynamics (CFD) codes such as FLUENT. The goal of this study is to limit the use of MCRT method to determining the appropriate input boundary condition on the outside of the window of the receiver and to use the built-in Discrete Ordinates (DO) method for all the radiation internal to the receiver and leaving the receiver due to emission. To reach the goal, this paper focuses on the DO method implemented within FLUENT. An earlier study on this subject is based and advanced. Appropriate radiation input for the DO method is extensively discussed. MIRVAL is used to simulate the heliostat field and VEGAS is used to simulate a lab-scale solar simulator; both of these codes utilize the MCRT method and provide intensity information on a surface. Output from these codes is discretized into DO parameters allowing the solution to proceed in FLUENT. Suitable benchmarks in FLUENT are used in a cylindrical geometry representing the receiver for the comparison and validation. This method will allow FLUENT to be used for a variety of problems involving concentrated solar energy.

Author(s):  
Eugene Cho ◽  
Fletcher Miller

The Combustion and Solar Energy Laboratory (C&SEL) at San Diego State University is developing a Small Particle Heat Exchange Receiver (SPHER) to absorb and transfer heat from concentrated solar radiation to a working fluid for a gas turbine. The SPHER is to be used with a Concentrated Solar Power (CSP) system where a heliostat field highly concentrates solar radiation on the optical aperture of the SPHER. The solar radiation is volumetrically absorbed by a unique carbon nanoparticle gas mixture within the cavity of the SPHER. This research focuses on comparing a Computational Fluid Dynamics (CFD) model using the ANSYS FLUENT Discrete Ordinates (DO) Model and a program developed by the C&SEL which uses a Monte Carlo Ray Trace (MCRT) method to calculate the spatial and directional distribution of radiation for an idealized solar receiver geometry. Previous research at the C&SEL has shown successful implementation of the MCRT method to calculate the spatial and directional distribution of radiation for an idealized solar receiver geometry. The MCRT method is highly accurate and will serve as the benchmark solution for this research. However the MCRT code takes several days to run, is inflexible to geometry changes, and is cumbersome to implement as the MCRT code needs to be rewritten for each new receiver geometry being considered. These factors necessitate the need to find an alternate method that accurately calculates the spatial and directional distribution of radiation for a solar receiver and can be efficiently implemented for various receiver geometries being studied. The Discrete Ordinates method is a new method for solving the Radiative Transport Equations (RTE) using a FORTRAN program, developed by the C&SEL, and the ANSYS FLUENT Discrete Ordinates model for calculating the RTE. The FORTRAN program calculates the proper inlet radiation boundary conditions that ANSYS FLUENT uses for calculating the RTE. The methodology used for determining the correct CFD mesh, radiative boundary conditions, optimal number of DO theta and phi discretization, as well as the optical properties of the working fluid will be presented in this paper. The main focus of this research is to compare two different methods for solving the Radiative Transport Equations within the idealized SPHER. The solution data for several cases using the previous coupled MCRT method and the ANSYS FLUENT Discrete Ordinates method is presented for both a collimated and diffuse gray radiation approximation. The case studies focus on researching how the MCRT method and Discrete Ordinates method differ when comparing critical receiver parameters such as the mean outlet temperature, wall temperature profile, outlet tube temperature profile, and total receiver efficiency while keeping the total inlet radiation flux of 5 MW and inlet mass flow rate of 5 kg/s constant. This research also presents a study on the optimal Discrete Ordinates angular discretization, as well as a study to determine the solution’s dependence on the number of inlet boundary conditions imposed on the window.


1990 ◽  
Vol 112 (3) ◽  
pp. 169-176 ◽  
Author(s):  
Douglas R. Adkins

Heat pipes are being developed to transfer solar energy from the focal point of a parabolic dish concentrator to the working fluid of Stirling engines. With these receivers, concentrated solar energy that is absorbed on the concave surface of a dome is removed by the evaporation of liquid sodium on the convex side of the dome. Vaporized sodium then condenses on an engine’s heater tubes and transfers energy to the working fluid of the engine. The condensed sodium returns to the absorber surface where it is redistributed across the dome by the capillary action of a wick. Issues concerning the flow of sodium in a heat-pipe solar receiver are investigated in this paper. A comparison is made between various wick options, and general issues concerning the design of heat-pipe receivers are also discussed.


Author(s):  
Jo¨rg Petrasch

A free and open source Monte Carlo ray-tracing program for concentrating solar energy research and development is presented. The program uses non energy partitioning Monte Carlo methods to model radiative exchange between arbitrarily arranged surfaces. Surface models include concentrating geometries, such as spherical, parabolic, and elliptical concentrators as well as compound parabolic concentrators. The program’s modular design allows implementation of additional surface and source models. The program has been thoroughly tested and experimentally validated. It has been used to model several concentrating devices including PSI’s high flux solar furnace and ETH’s high flux solar simulator. Furthermore, it has been used to design PSI’s high flux solar simulator and UFL’s high flux solar simulator. The code is particularly suited to provide radiative boundary conditions for numerical models of high temperature solar receivers and solar thermochemical reactors.


2021 ◽  
Vol 5 (2) ◽  
pp. 16
Author(s):  
Isabel Padilla ◽  
Maximina Romero ◽  
José I. Robla ◽  
Aurora López-Delgado

In this work, concentrated solar energy (CSE) was applied to an energy-intensive process such as the vitrification of waste with the aim of manufacturing glasses. Different types of waste were used as raw materials: a hazardous waste from the aluminum industry as aluminum source; two residues from the food industry (eggshell and mussel shell) and dolomite ore as calcium source; quartz sand was also employed as glass network former. The use of CSE allowed obtaining glasses in the SiO2-Al2O3-CaO system at exposure time as short as 15 min. The raw materials, their mixtures, and the resulting glasses were characterized by means of X-ray fluorescence, X-ray diffraction, and differential thermal analysis. The feasibility of combining a renewable energy, as solar energy and different waste for the manufacture of glasses, would highly contribute to circular economy and environmental sustainability.


Author(s):  
S. Kh. Suleimanov ◽  
V. G. Babashov ◽  
M. U. Dzhanklich ◽  
V. G. Dyskin ◽  
M. I. Daskovskii ◽  
...  

Catalysts ◽  
2021 ◽  
Vol 11 (6) ◽  
pp. 723
Author(s):  
Mahesh Muraleedharan Nair ◽  
Stéphane Abanades

The CeO2/CeO2−δ redox system occupies a unique position as an oxygen carrier in chemical looping processes for producing solar fuels, using concentrated solar energy. The two-step thermochemical ceria-based cycle for the production of synthesis gas from methane and solar energy, followed by CO2 splitting, was considered in this work. This topic concerns one of the emerging and most promising processes for the recycling and valorization of anthropogenic greenhouse gas emissions. The development of redox-active catalysts with enhanced efficiency for solar thermochemical fuel production and CO2 conversion is a highly demanding and challenging topic. The determination of redox reaction kinetics is crucial for process design and optimization. In this study, the solid-state redox kinetics of CeO2 in the two-step process with CH4 as the reducing agent and CO2 as the oxidizing agent was investigated in an original prototype solar thermogravimetric reactor equipped with a parabolic dish solar concentrator. In particular, the ceria reduction and re-oxidation reactions were carried out under isothermal conditions. Several solid-state kinetic models based on reaction order, nucleation, shrinking core, and diffusion were utilized for deducing the reaction mechanisms. It was observed that both ceria reduction with CH4 and re-oxidation with CO2 were best represented by a 2D nucleation and nuclei growth model under the applied conditions. The kinetic models exhibiting the best agreement with the experimental reaction data were used to estimate the kinetic parameters. The values of apparent activation energies (~80 kJ·mol−1 for reduction and ~10 kJ·mol−1 for re-oxidation) and pre-exponential factors (~2–9 s−1 for reduction and ~123–253 s−1 for re-oxidation) were obtained from the Arrhenius plots.


Sign in / Sign up

Export Citation Format

Share Document