Exploring irradiated granular flows with rapid heating for concentrated solar thermal energy collection and storage
Not Available
Engineering topics
Publications and source records attributed to Loutzenhiser, Peter G..
Not Available
Utilizing thermochemical reactions for thermal energy storage and solar fuel production has been an emerging research topic. Thermal transport properties of the materials are an important parameter that can determine the kinetics and efficiency of thermochemical reactions. With the increasing number of new thermochemical materials (TCMs); however, there is a lack of reliable techniques to monitor the thermal transport property of the materials and their changes as a function of reactions in real time. In this work, we report the in situ monitoring of thermochemical reactions using modulated photothermal radiometry (MPR). The thermal conductivities of two TCMs, namely, calcium hydroxide (Ca(OH) 2 ) and Ba 0.15 Sr 0.85 FeO 3–δ (BSF1585), were measured as a function of temperature and time using the MPR technique. The measured thermal conductivities were correlated to the reaction. The work has two significant contributions to the research communities. First, it provides a non-invasive diagnostic tool for monitoring the thermal transport properties of TCMs that can potentially be a high-throughput measurement technique conducive to optimizing TCMs, reactors, and related thermal systems. Second, for TCMs that show observable changes in thermal transport properties, a correlation between the measured thermal conductivity and the conversion fraction of the reaction can be established for monitoring the reaction kinetics based on thermal characterization.
A two-step solar thermochemical looping cycle based on Co 3 Mo 3 N/Co 6 Mo 6 N reduction/nitridation reactions offers a pathway for green NH 3 production that utilizes concentrated solar irradiation, H 2 O, and air as feedstocks. The NH 3 production cycle steps both derive process heat from concentrated solar irradiation and encompass 1) the reduction of Co 3 Mo 3 N in H 2 to Co 6 Mo 6 N and NH 3 ; and 2) nitridation of Co 6 Mo 6 N to Co 3 Mo 3 N with N 2 . Co 3 Mo 3 N reduction/nitridation reactions are examined at different H 2 and/or N 2 partial pressures and temperatures. NH 3 production is quantified in situ using liquid conductivity measurements coupled with mass spectrometry (MS). Solid-state characterization is performed to identify a surface oxygen layer that necessitates the addition of H 2 during cycling to prevent surface oxidation by trace amounts of O 2 . H 2 concentrations of > 5% H 2 /Ar and temperatures >500 °C are required to reduce Co 3 Mo 3 N to Co 6 Mo 6 N and form NH 3 at 1 bar. Complete regeneration of Co 3 Mo 3 N from Co 6 Mo 6 N is achieved at conditions of 700 °C under 25–75% H 2 /N 2 . H 2 pressure-swings are observed to increase NH 3 production during Co 3 Mo 3 N reduction. In conclusion, the results represent the first comprehensive characterization of and definitive non-catalytic production of NH 3 via chemical looping with metal nitrides and provide insights for technology development.
Bauxite and silica particles have gained increasing attention for applications in the field of concentrated solar power. In this work, a Monte Carlo ray-tracing simulation is performed to predict the radiative properties (absorptance, reflectance, and transmittance) of packed beds with mixed bauxite and silica spherical particles at wavelengths of 0.5, 2.6, 9.0, and 9.35 µm. These wavelengths are representative for the visible, near-infrared, and the mid-infrared regions that are important for solar and thermal radiation, respectively. A repeating unit-column approach is used to mathematically represent the particle bed. The effects of particle mixing ratio, volume fraction, and wavelength on the predicted radiative properties are examined. The obtained radiative properties are inputted to an inverse method to retrieve the effective absorption and scattering coefficients as well as the scattering albedo, which may be later used in a continuous-scale radiative heat transfer analysis. Furthermore, the independent scattering model is used to obtain the absorption coefficient and scattering albedo based on the absorption and scattering cross sections predicted by a Monte Carlo algorithm for a single particle. It is shown that the independent scattering model underpredicts the scattering coefficient for opaque particles but overpredicts the scattering coefficient for semitransparent particles for sufficiently high particle volume fractions. Here, the radiative properties calculated from the independent scattering model are compared to the full Monte Carlo simulation of the particle bed to examine the influence of particle mixing on dependent scattering.
High-temperature, dense granular flows along an inclined plane were considered for solar thermal energy transport and storage with sintered bauxite particles. A series of experiments was performed for particle inlet temperatures of ~ 200, 400, 600, and 800 °C to understand the mechanisms of granular flows at extreme temperatures. Mass flow rates were measured using a load cell and free-surface velocities were measured and computed using particle image velocimetry. Surface temperatures were measured using infrared cameras. A significant decrease in steady-state particle mass flow rate was observed with increasing temperature due to changing flow properties. A decrease in bulk particle free-surface velocities was observed at higher temperatures. Free-surface velocity measurement error between experiments were within 20% of the average. The particle surface temperatures decreased from inlet to outlet with larger gradients at higher temperatures observed due to increasing convection and radiative heat losses. Here, a decrease in temperature was observed along the side walls due to a decrease in particle velocities.
Bauxite and silica particles are candidate materials for solar thermal energy storage at high temperatures. The temperature-dependent emittance of packed beds with bauxite and silica particles was measured using a newly upgraded emissometer at wavelengths 2 μm ≤ λ ≤ 16 μm and temperatures up to ~730 K. The room-temperature emittance was obtained from the measured directional-hemispherical reflectance. A fused silica disc was used to test the emissometer by comparing the measured spectral emittance with the calculated emittance from a fitted Lorentz oscillator model. For the polycrystalline silica particles and the fused silica disc, the measured emittance increases with temperature in the mid-infrared region. The underlying mechanism is interpreted as the temperature-dependent damping coefficient in the Lorentz oscillator model. Two types of bauxite particles with different compositions and sizes were investigated. For λ > 10 μm, the measured emittance at elevated temperatures is higher than that at room temperature. In the region 2 μm < λ < 6 μm, the temperature dependence varies for different types of particles. The total emittance of bauxite particle beds was calculated by spectral integration using Planck’s distribution at the prescribed temperature. As a result, the calculated total emittance is between 0.89 and 0.96, but it does not change monotonically with temperature.
CaAl 0.2 Mn 0.8 O 3-δ (CAM28) and CaTi 0.2 Mn 0.8 O 3-δ (CTM28) are perovskite metal oxides developed for high-temperature thermochemical energy storage (TCES) applications, e.g., in support of air Brayton power generation. Previous reports for these compounds focus on the equilibrium non-stoichiometry (δ) as a function of temperature and oxygen partial pressure (pO 2 ) and the endotherm (or exotherm) accompanying changes in δ resulting from thermal reduction (or re-oxidation). Herein, we report results for elemental substitution and doping (Al, Co, Fe, La, Sr, Ti, Y, Zn, and Zr) of calcium manganites (CM) that establish the preference for CAM28 and CTM28. Techniques employed include conventional (screening and equilibrium) and ballistically heated multi-cycle thermogravimetric analysis (TGA), conventional and high temperature ( in-situ ) X-ray diffraction (XRD), and differential scanning calorimetry (DSC). Forward-looking results for A-site Y-doped materials, e.g., Ca 0.9 Y 0.1 MnO 3-δ (CYM910), establish a route to increasing the reduction enthalpy relative to CAM28 and CTM28, albeit at the expense of increased reduction temperatures and raw materials costs. A thermodynamic model presented for CAM28, but extendable to related materials, provides values for the reaction enthalpy and extent of reduction as a function of temperature and oxygen partial pressure for use in design efforts. Taken as a whole, the results support the choice of Al-doped CaMnO 3-δ as a low-cost material for TCES in a high temperature air Brayton application, but point the way to achieving higher stored energy densities that could lead to overall cost savings.
The focus of this work is on the measurement and analysis of the radiative properties of polycrystalline SiO 2 particle beds with various layer thicknesses. The particles are polydispersed with average diameters of 222, 150, and 40 μm . The spectral, directional–hemispherical reflectance and transmittance of the particle bed are measured at wavelengths from 0.4 to 1.8 μm using a monochromator, and the reflectance measurement is extended to 15 μm using a Fourier-transform infrared spectrometer. Particles are closely packed between two transparent windows for measuring the radiative properties. In the visible and near-infrared region up to 1.8 μm, the inverse adding–doubling method yields the effective absorption and scattering coefficients. The results suggest that short wavelength absorption needs to be included in modeling the behavior of particle beds due to multiple scattering. A discrete-scale Monte Carlo ray-tracing method is developed to model the radiative properties by assuming monodispersed spherical particles, and the simulated results compare well with measurements. The effective absorption and scattering coefficients of the particle beds obtained from the independent scattering theory are compared to those from the inverse method. As a result, the impact of dependent scattering on the packed beds is observed for smaller-sized particles.
Submillimeter solid particles have been considered as thermal storage media for concentrated solar power applications. Knowledge of the scattering and absorption properties of individual particles is crucial for modeling the radiative heat transfer of the particle bed. In this work, a laser scatterometer is used to measure the single-particle scattering properties at a wavelength of 635 nm by using two configurations: (1) a falling particle curtain and (2) a taped particle layer. Because a one-particle nominal thickness is formed with area fractions of 5-55% depending on the configuration, multiple scattering is minimized and hence the single scattering phase function, averaged over all illuminated particles, is directly measured. Bauxite-based ceramic particles that are strongly absorbing in the solar spectrum and silica particles that are nonabsorbing in the visible and near-infrared are investigated. The directional-hemispherical reflectance and transmittance of the taped particles are also measured to deduce the forward and backward scattering efficiency factors and the absorption efficiency factors. Only weak wavelength dependence is observed in the measured region from 380 nm to 1020 nm. Furthermore, the scattering phase functions of all bauxite-based particles with varying sizes and compositions are very similar and can be fitted to a Henyey-Greenstein phase function with an asymmetry factor g = -0.20. For the silica particles, forward scattering dominates and g = 0.45 yields the best fit. A Monte Carlo method is developed to model the particle scattering characteristics, and reasonable agreements between the modeling and experimental results are observed by introducing a specularity parameter.
Two-step solar thermochemical cycles based on reversible reactions of SrFeO 3–δ and (Ba,La) 0.15 Sr 0.85 FeO 3–δ perovskites were considered for air separation. The cycle steps encompass (1) the thermal reduction of SrFeO 3–δ or (Ba,La) 0.15 Sr 0.85 FeO 3–δ perovskites driven by concentrated solar irradiation and (2) oxidation in air to remove O 2 and produce N 2 . Rate limiting mechanisms were examined for both reactions using a combination of isothermal and non-isothermal thermogravimetry for temperature-swings between 673 and 1373 K, heating rates of 10, 20, and 50 K min –1 , and O 2 pressure-swings between 20% O 2 /Ar and 100% Ar at atmospheric pressure. Evolved O 2 and associated lag due to transport behavior were measured with gas chromatography and used with measured sample temperatures to predict equilibrium compositions from a compound energy formalism thermodynamic model. Measured and predicted chemical equilibrium changes in deviation from stoichiometry were compared. Rapid chemical kinetics were observed as the samples equilibrated rapidly for all conditions, indicative that heat and mass transfer were the rate limiting mechanisms. Furthermore, the effects of bulk diffusion (or gas diffusion through the bed or pellet) were examined using pelletized and loose powdered samples and determined to have no discernable impact.
A two-step solar thermochemical cycle was considered for air separation to produce N2 based on (Ba,La) x Sr 1-x FeO 3-δ perovskite reduction/oxidation (redox) reactions for A-site fractions of 0 ≤ x ≤ 0.2. The cycle steps encompassed (1) thermal reduction and O 2 release via concentrated solar input and (2) re-oxidation with air to uptake O 2 and produce high-purity N 2 . Thermogravimetry at temperatures between 400 and 1100 °C in atmospheres of 0.005 to 90% O 2 /Ar at 1 bar was performed to measure equilibrium nonstoichiometries. The compound energy formalism was applied to model redox thermodynamics for both Ba 2+ and La 3+ substitution. Non-linear regression was used to determine the empirical parameters based on the thermogravimetric measurements. The model was used to define partial molar reaction enthalpies and entropies and predicted equilibrium oxygen nonstoichiometry as functions of oxide stoichiometry, site fraction, temperature, and O 2 partial pressure. The thermodynamic analysis showed the materials are appealing for air separation at temperatures below 800 °C.
Abstract not provided.
High temperature particle flow properties necessary to predict granular flow behavior for solar thermal energy storage applications were measured and calculated for Carbobead CP 30/60 up to 800 °C. Here, the measured properties included elastic and shear moduli, particle-particle coefficients of static sliding and rolling friction, and particle-particle coefficients of restitution. Poisson’s ratio was calculated with elastic and shear moduli. The flow properties were used as inputs for a numerical model using the discrete element method to examine granular flows along an inclined plane at high temperature. The flow behavior was strongly influenced by the coefficients of static friction, which impacted the particle residence time, shear effects from the side walls, and particle flow mass flux. An 8.7%, 15.6%, and 8.5% increase and 37.9% decrease in steady state mass flow rate was observed for 200 °C, 400 °C, 600 °C, and 800 °C, respectively, when compared to room temperature simulations. A 52%, 59%, and 33% decrease in the time to reach steady state was observed for 200 °C, 400 °C, and 600 °C, respectively, while a 53% increase in time was observed for 800 °C. A significant delay in the flow development at 800 °C was observed due to significantly higher frictional forces.
An A‑ and B‑site substitutional study of SrFeO3−δ perovskites (A’xA1−xB’yB1−yO3−δ, where A = Sr and B = Fe) was performed for a two‑step solar thermochemical air separation cycle. The cycle steps encompass (1) the thermal reduction of A’xSr1−xB’yFe1−yO3−δ driven by concentrated solar irradiation and (2) the oxidation of A’xSr1−xB’yFe1−yO3−δ in air to remove O2, leaving N2. The oxidized A’xSr1−xB’yFe1−yO3−δ is recycled back to the first step to complete the cycle, resulting in the separation of N2 from air and concentrated solar irradiation. A-site substitution fractions between 0 ≤ x ≤ 0.2 were examined for A’ = Ba, Ca, and La. B-site substitution fractions between 0 ≤ y ≤ 0.2 were examined for B’ = Cr, Cu, Co, and Mn. Samples were prepared with a modified Pechini method and characterized with X-ray diffractometry. The mass changes and deviations from stoichiometry were evaluated with thermogravimetry in three screenings with temperature- and O2 pressure-swings between 573 and 1473 K and 20% O2/Ar and 100% Ar at 1 bar, respectively. A’ = Ba or La and B’ = Co resulted in the most improved redox capacities amongst temperature- and O2 pressure-swing experiments.
Abstract not provided.
A novel method for pairing surface irradiation and volumetric absorption from Monte Carlo ray tracing to computational heat transfer models is presented. The method is well-suited to directionally and spatially complex concentrated radiative inputs (e.g., solar receivers and reactors). The method employs a generalized algorithm for directly mapping absorbed rays from a Monte Carlo ray tracing model to boundary or volumetric source terms in the computational mesh. The algorithm is compatible with unstructured, two and three-dimensional meshes with varying element shapes. Four case studies were performed on a directly irradiated, windowed solar thermochemical reactor model to validate the method. The method was shown to conserve energy and preserve spatial variation when mapping rays from a Monte Carlo ray tracing model to a computational heat transfer model in ansys fluent.