Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “REACTOR LATTICES”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

Study of iron structure stability in high temperature molten lead-bismuth eutectic with oxygen injection using molecular dynamics simulation

Corrosion of structural materials in high temperature molten lead-bismuth eutectic is a major problem for design of PbBi cooled reactor. One technique to inhibit corrosion process is to inject oxygen into coolant. In this paper we study and focus on a way of inhibiting the corrosion of iron using molecular dynamics method. For the simulation results we concluded that effective corrosion inhibition of iron may be achieved by injection 0.0532 wt% to 0.1156 wt% oxygen into liquid lead-bismuth. At this oxygen concentration the structure of iron material will be maintained at about 70% in bcc crystal structure during interaction with liquid metal.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Time-Resolved PIV/PTV Measurements on Interior Subchannels of a Wire-Wrapped 61-pin Hexagonal Fuel Bundle

Liquid metal fast reactor fuel bundles using sodium as a coolant typically utilize a tightly packed triangular lattice of fuel pins enclosed in a hexagonal duct. Often, the fuel pins are helically wrapped with a wire spacer. The primary function of the wire spacers is to maintain a gap between adjacent fuel pins. They also enhance sub-channel mixing to increase convective heat transfer and mitigate vortex-induced vibration. Texas A and M University has conducted isothermal flow experiments in a wire-wrapped 61-pin hexagonal fuel bundle to support the research on advanced nuclear fuel development sponsored by the Department of Energy (DOE). These tests produced high spatial and temporal resolution measurements of the velocity field and pressure at different locations in the fuel bundle. The experimental facility employs matched-index-of-refraction (MIR) techniques and laser diagnostic velocity measurement techniques such as particle image velocimetry (PIV), particle tracking velocimetry (PTV), and laser Doppler velocimetry (LDV). The result of this experimental activity will be an experimental database of pressure and velocity field measurements that will be suitable for validating computational fluid dynamics (CFD) methods, such as Reynolds-averaged Navier-Stokes (RANS) and highly resolved large eddy simulation (LES). PIV results on the set of edge subchannels near the hexagonal duct have been presented in another paper. This paper presents PIV results on two sets of interior subchannels, defined as Interior- 1 and Center-2, at a Reynolds number (Re) of 19,000 in the turbulent regime.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Particle Image Velocimetry Measurements in a Wire-Wrapped 61-Pin Hexagonal Fuel Bundle

Liquid metal fast reactor fuel assemblies using sodium as a coolant typically utilize a tightly packed triangular lattice of fuel pins enclosed in a hexagonal duct. The fuel pins are typically wrapped with a helical wire spacer. The primary function of the wire spacers is to maintain a gap between neighboring fuel pins. They also mitigate vortex-induced vibration and enhance sub-channel mixing to increase convective heat transfer. Many experiments have been performed to study thermal-hydraulic characteristics, primarily pressure drops, of the wire-wrapped fuel bundles. Novendstern developed a semi-empirical model to predict pressure losses in a hexagonal array of wire-wrapped pins in the turbulent flow regime. Rehme proposed pressure drop correlations for a wide range of geometrical parameters, such as pitch-to-diameter ratio P/D from 1.125 to 1.417, lead-to-wire diameter ratio H / dm from 6 to 45, and number of pins n from 7 to 61. Cheng and Todreas introduced their hydrodynamic models for subchannel friction factors and mixing parameters in a 37-pin fuel bundle, which covered the laminar, transition and turbulent flow regimes. Computational fluid dynamics (CFD) offers an opportunity to validate the above correlations and investigate the flow characteristics in the subchannels. Several numerical studies of the wire-wrapped fuel bundles have been performed and most of them are Reynolds-Averaged Navier-Stokes (RANS) based simulations with various turbulence models. Ahmad and Kim performed RANS simulations of 7- and 19-pin bundles using the k-omega SST turbulent model, while Gajapathy et al. performed RANS simulations on 7-pin bundles using the k-epsilons model. Recently, Fischer et al. carried out large-eddy simulation (LES) calculations on wire-wrapped fuel bundles to obtain higher-resolution numerical solutions of subchannel behavior. The geometrical complexity of the hexagonal wire-wrapped fuel bundle has limited experimental activities to pressure drop measurements. A database of pressure drops and flow-field measurements in a 61-pin wire-wrapped hexagonal fuel bundle is not available to benchmark the existing correlations and validate the CFD calculations. Texas A and M University is conducting isothermal flow experiments in a wire-wrapped 61-pin hexagonal fuel bundle to support the research on advanced nuclear fuel development sponsored by the Department of Energy (DOE). The purpose of these tests is to perform high spatial and temporal resolution measurements of the flow and pressure measurements at different locations in the fuel bundle. The experimental facility employs matched-index-of-refraction (MIR) techniques and laser diagnostic velocity measurement techniques such as particle image velocimetry (PIV) and laser Doppler velocimetry (LDV). The result of this experimental activity will be an experimental database of pressure and flow-field measurements that will be suitable for validating CFD computer codes. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Vibrational modes and quantum zero-point energy of hydrogen in ZrH 0.0155 and ZrH 2

We report on an inelastic neutron scattering study of the proton dynamics in ZrH 0.0155 and ε-ZrH 2 . In particular, we present measurements of the incoherent dynamic structure factor, generalized vibrational density of states, and proton momentum distribution of these two materials. Our results are generally consistent with theoretical predictions of Elsässer et al. [Mat. Res. Soc. Symp. Proc. 453 221–226 (1997)]. They argued that the effective Born-Oppenheimer potential experienced by the hydrogen atoms in ε-ZrH 2 is nearly isotropic and harmonic at energies below 0.3 eV, but becomes anisotropic and anharmonic for higher energies. At low temperatures, the proton momentum distribution is dominated by the quantum-mechanical ground state of the protons. We find that it assumes a Gaussian shape, consistent with the concept that the potential surface is approximately harmonic for small displacements of the hydrogen atoms. However, the anharmonicity of the potential becomes readily apparent in the excited states of the hydrogen atoms, as the harmonic approximation breaks down in the description of the multiphonon bands.

36 MATERIALS SCIENCE↗

IER-539 CED-1: Preliminary Design of a New Horizontal Split Table

This report presents the preliminary design (CED-1) of IER-539, focused on a new Horizontal Split Table (HST) critical assembly machine. The Department of Energy (DOE) National Nuclear Criticality Safety Program (NCSP) funds critical experiment R&D and its Mission and Vision document identifies a long-standing technical gap of a general purpose HST. HSTs have historically been used for experiments that are impractical to assemble on a vertical lift machine (VLM), usually due to their large physical size or weight. While many HSTs were historically employed in the US (including at LLNL), there are currently no surviving HSTs in the US. While there are two critical experiments facilities remaining in the US, the National Critical Experiments Research Center (NCERC) operated by LANL and the Critical Experiments Facility (CX) operated by SNL, these facilities are not equipped with a general purpose HST. There are a number of experiment designs that require an HST, including needs for Defense Programs and Nuclear Energy. Lack of an HST constraints the national ability to perform large and heavy experiments. Large footprints are needed to successfully create mockup experiments for advanced reactors (mock-up cold critical reactor cores, which is a urgent need for new reactor designs using novel materials and longer cycle lengths), to create solution experiments in lattices, or to test nuclear data of heavy metals. A new HST is a vital piece of criticality experiment infrastructure that would allow for the conduct of integral critical experiments to address needs in the nuclear criticality safety, nuclear data, and other DOE stakeholder communities.

42 ENGINEERING↗

Improved Modelling for the Assembly Discontinuity Factors in the APEC Method

Modern reactor analysis of light water reactor depends largely on a two-step procedure based on the simplified equivalence theory (SET). Consequently, the resulting homogenized two-group cross-section (XS) and associated discontinuity factor (DF) are subject to noticeable, intrinsic uncertainties. In our previous study, we introduced the albedo-corrected parameterized equivalence constants (APEC) method, which is a new leakage correction method applicable to the fuel assembly homogenized two-group XS of the standard two-step procedure. In the APEC method, fuel assembly two-group XS is expressed as a simple polynomial function of assembly-wise current-to-flux ratio (CFR), and it was demonstrated that two-group XSs can be well corrected during a nodal calculation and the nodal accuracy is improved a lot. However, since the conventional assembly discontinuity factor (ADF) is still used in the current APEC method, a small but non-negligible error is unavoidable. Therefore, it is necessary to include improved DF model in the APEC method for a higher accuracy. The functionalized discontinuity factor (FDF) using the definition of CFR was investigated by authors before the development of the APEC method, where it was motivated by the boundary perturbation theory by F. Rahnema. However, in our previous study on FDF, lattice calculations were performed only with one-dimensional two-group diffusion equation and could not be applied to actual core analysis. In this paper, feasibility study of DF functionalization is performed based on a 2-D method of characteristics (MOC) lattice code, DeCART2D.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Effect of Cellular Lattices, Heat-Exchange Intensifiers, and Coolant Mixing on the Critical Power of VVER FA

The results of one stage in the development of methods for increasing the critical power of water-cooled nuclear reactors – the capacity of a reactor installation with sharp degradation of the heat exchange characteristics of the coolant (crisis of heat exchange) – are presented. This stage included experimental determination of the effect on the crisis of heat exchange of arranging heat-exchange intensifiers in FA and transverse mixing of the coolant flow, for which lattices assembled from properly shaped cellular elements where used. The FA models consisted of an assembly of 19 fuel-element simulators with a non-uniform distribution of heat emission along the length. Freon-12 was used as the coolant. The models of the FA differed by the type of intensifier lattices and the number of FA arranged along the length. The efficacy of the influence of lattices on the critical power as a result of the intensification of heat exchange and transverse mixing of the coolant was determined separately. The optimal rigging of FA with lattices was determined experimentally.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Minor actinide transmutation in thorium and uranium matrices in heavy water moderated reactors

The irradiation of Th{sup 232} breeds fewer of the problematic minor actinides (Np, Am, Cm) than the irradiation of U{sup 238}. This characteristic makes thorium an attractive potential matrix for the transmutation of these minor actinides, as these species can be transmuted without the creation of new actinides as is the case with a uranium fuel matrix. Minor actinides are the main contributors to long term decay heat and radiotoxicity of spent fuel, so reducing their concentration can greatly increase the capacity of a long term deep geological repository. Mixing minor actinides with thorium, three times more common in the Earth's crust than natural uranium, has the additional advantage of improving the sustainability of the fuel cycle. In this work, lattice cell calculations have been performed to determine the results of transmuting minor actinides from light water reactor spent fuel in a thorium matrix. 15-year-cooled group-extracted transuranic elements (Np, Pu, Am, Cm) from light water reactor (LWR) spent fuel were used as the fissile component in a thorium-based fuel in a heavy water moderated reactor (HWR). The minor actinide (MA) transmutation rates, spent fuel activity, decay heat and radiotoxicity, are compared with those obtained when the MA were mixed instead with natural uranium and taken to the same burnup. Each bundle contained a central pin containing a burnable neutron absorber whose initial concentration was adjusted to have the same reactivity response (in units of the delayed neutron fraction β) for coolant voiding as standard NU fuel. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Preparation of nuclear libraries with deterministic and stochastic methods for LWR reflectors

The explicit reflector methodology is used in the system of codes CASMO-5 / SIMULATE-3 to include the reflector around the active core into the computational region and avoid adopting any ad-hoc or experimental albedo coefficients as boundary conditions. However, to complete the core calculation, a set of cross sections and discontinuity factors is needed for the reflector nodes and the accuracy of these nuclear parameters influences the final results, in particular along the peripheral regions of the core. In this paper the explicit reflector methodology of CASMO-5 is adopted to evaluate the few-group cross sections and discontinuity factors of the different reflector cases, based on the design of Generation II and III LWR reactors. In addition, in the perspective of using Monte-Carlo codes as a complementary option for lattice calculations of reflector configurations, the stochastic SERPENT code is also included as part of this benchmark. With the latter, the impact of applying 2-D reflector models with homogenized materials instead of explicit representation of the real geometrical structures is moreover evaluated and shown to be limited. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Adjusted Graphite Thermal Scattering Data Applied to Fluoride-Salt-Cooled Reactors

Thermal scattering data is a key component of fission reactor simulation. While most materials can be reasonably modeled with a free gas approximation, the moderators in thermal systems require a more rigorous treatment to sufficiently capture their specific neutron interaction physics. Graphite does not typically feature as a principle moderator in power reactors in service today, but remains an attractive prospect for a number of new reactor concepts currently in development. In comparison to light water, a material nearly ubiquitous in thermal reactors, graphite offers a better neutron economy by way of a small absorption cross section, high thermodynamic stability and can be incorporated into the structure of the core itself. Graphite in its purest form is a crystalline material composed of stacked layers of graphene, a hexagonal planar lattice of carbon atoms. In-plane, each carbon atom is covalently bonded to 3 neighbouring carbon atoms, creating a strong bond, while between planes only a weak van der Waals force exists. Pure graphite has a density of 2.26 g/cm{sup 3}. However, reactor-grade contains a significant amount of carbon binder, leading to a porous bulk medium with randomly oriented micro crystals. Densities of 1.5 - 1.9 g/cm{sup 3} are typical for reactor-grade graphite, corresponding to porosities on the order of 10-30%. This deviation from pure crystalline graphite should be accounted for in the thermal scattering data. When an incident neutron is in or near the thermal region, its de Broglie wavelength is the same order as the molecular bond lengths, leading to elastic scattering. At similar energies, the dynamics of the moderator lattice (in the case of solid moderators) can cause inelastic scattering. For these moderators, special cross sections should be calculated and used in reactor simulations. To generate this thermal scattering data, parameters pertaining to the properties and structure of the material in question are fed into a nuclear data processing code, such as NJOY. Within NJOY is a module called LEAPR which is used to generate the thermal scattering law S(α,β), which is in turn is fed into another module THERMR which is used to generate thermal cross sections that factor in the molecular-scale scattering phenomena. The core of the thermal scattering law is the phonon spectrum ρ(β) which describes the modes of excitation of the material lattice. An investigation into the phonon frequency spectrum was conducted at North Carolina State University (NCSU) for reactor-grade graphite which accounts for porosity. This phonon spectrum was calculated using the ab initio density functional theory (DFT) code VASP. The inelastic cross section calculated using this spectrum was found to have good agreement with relevant experimental data. The resulting spectra differed from the benchmark spectrum produced by Young and Koppel, the same data used in current nuclear data evaluations. The latter uses the lattice dynamics of perfect crystalline graphite with four force constants fitted to the model from thermodynamics data from reactor-grade (porous) graphite, which yields a reasonable approximation to the phonon spectrum of reactor-grade graphite.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Nuclear reactions from lattice QCD

In this study, one of the overarching goals of nuclear physics is to rigorously compute properties of hadronic systems directly from the fundamental theory of strong interactions, Quantum Chromodynamics (QCD). In particular, the hope is to perform reliable calculations of nuclear reactions which will impact our understanding of environments that occur during big bang nucleosynthesis, the evolution of stars and supernovae, and within nuclear reactors and high energy/density facilities. Such calculations, being truly ab initio, would include all two-nucleon and three- nucleon (and higher) interactions in a consistent manner. Currently, lattice QCD provides the only reliable option for performing calculations of some of the low-energy hadronic observables. With the aim of bridging the gap between lattice QCD and nuclear many-body physics, the Institute for Nuclear Theory held a workshop on Nuclear Reactions from Lattice QCD on March 2013. In this review article, we report on the topics discussed in this workshop and the path planned to move forward in the upcoming years.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Demonstration of Pronghorn’s Subchannel Code Modeling of Liquid-Metal Reactors and Validation in Normal Operation Conditions and Blockage Scenarios

Pronghorn-SC is a subchannel code within the Multiphysics Object-Oriented Simulation Environment (MOOSE). Initially designed to simulate flows in water-cooled, square lattice, subchannel assemblies, Pronghorn-SC has been expanded to simulate liquid-metal-cooled flows in triangular lattices, hexagonal subchannel assemblies. For this purpose, the algorithm of Pronghorn-SC was adapted to solve the subchannel equations as they are applicable to a hexagonal wire-wrapped sodium-cooled fast reactor. Cheng–Todreas models for pressure drop and cross-flow models were adopted and a coolant heat conduction term was added. To solve these equations, an improved implicit algorithm was developed robust enough to deal with the numerical issues, associated with low flow and recirculation phenomena. To confirm the prediction capability of Pronghorn-SC, calculations and comparisons with available experimental data of 19- and 37-pin assemblies were performed, as well as other subchannel codes. Finally, a flow blockage modeling feature was added. This capability was validated for both water-cooled square sub-assemblies and sodium-cooled hexagonal sub-assemblies, using experimental data of partially and fully blocked cases.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Chapter 5: Metamorphic Growth and Multijunction III-V Solar Cells

Metal organic vapor phase epitaxy (MOVPE) has proved to be the primary materials growth technique for low-cost, high-volume III-V multijunction (MJ) solar cells. This chapter shows the typical applications that drive volume production for MJ solar cells using MOVPE. The two commonly used structures for lattice-mismatched MJ solar cells are the upright metamorphic (UMM) structure and the inverted metamorphic (IMM) structure. The chapter discusses the key aspects of reactor technology pertinent to UMM and IMM solar cell growth. The main commercial low-pressure, As/P MOVPE reactor technology are vertical chambers where wafers are held on a planar rotating graphite susceptor, acting to aid the laminar flow or subplanetary rotation used to planarize the growth. The in situ monitoring becomes a critical component when considering lattice-mismatched growth, in particular the measurement of wafer curvature. One may relate the wafer curvature to the net strain at growth temperature through the use of Stoney's equation.

14 SOLAR ENERGY↗

Thermal plasma processed ferro-magnetically ordered face-centered cubic iron at room temperature

Here, we report tailor made phase of iron nanoparticles using homogeneous gas phase condensation process via thermal plasma route. It was observed that crystal lattice of nano-crystalline iron changes as a function of operating parameters of the plasma reactor. In the present investigation iron nanoparticles have been synthesized in presence of argon at operating pressures of 125–1000 Torr and fixed plasma input DC power of 6 kW. It was possible to obtain pure fcc, pure bcc as well as the mixed phases for iron nanoparticles in powder form as a function of operating pressure. The as synthesized product was characterized for understanding the structural and magnetic properties by using X-ray diffraction, vibrating sample magnetometer, and Mössbauer spectroscopy. The data reveal that fcc phase is ferromagnetically ordered with high spin state, which is unusual whereas bcc phase is found to be ferromagnetic as usual. Finally, the structural and magnetic properties are co-related.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Hydrogen trapping energetics at BCC iron-helium interfaces

Density functional theory (DFT) calculations have been performed to assess the trapping and segregation strength of hydrogen (H) to helium (He) interfaces in BCC (Fe) as a function of He density and surface orientation. The He density ranges from 0 to 2 He/V which corresponds to equilibrium bubbles that are expected to form in the structural components in a fusion reactor environment. The DFT calculations consist of a slab of BCC metal and an initial lattice of FCC He. Here, the report on the binding energy of H to these interfaces as well as the migration barriers into and away from the surface, which together provide information on the de-trapping energy for H at He bubbles. It was found that the binding energy of H to He-Fe interfaces decreases with increasing He density. Migration barriers are also modified due to the He density, which are sensitive to the surface orientation as well as the diffusion pathway. These results provide valuable first-principles energetics that are necessary for higher-order, mesoscale models, and provide insight into the extent to which He bubbles may trap tritium in fusion structural materials.

42 ENGINEERING↗

Development of ORIGEN Libraries for Mixed Oxide (MOX) Fuel Assembly Designs

In this research, ORIGEN cross section libraries for reactor-grade mixed oxide (MOX) fuel assembly designs have been developed to provide fast and accurate depletion calculations to predict nuclide inventories, radiation sources and thermal decay heat information needed in safety evaluations and safeguards verification measurements of spent nuclear fuel. These ORIGEN libraries are generated using two-dimensional lattice physics assembly models that include enrichment zoning and cross section data based on ENDF/B-VII.0 evaluations. Using the SCALE depletion sequence, burnup-dependent cross sections are created for selected commercial reactor assembly designs and a representative range of reactor operating conditions, fuel enrichments, and fuel burnup. The burnup dependent cross sections are then interpolated to provide problem-dependent cross sections for ORIGEN, avoiding the need for time-consuming lattice physics calculations. The ORIGEN libraries for MOX assembly designs are validated against destructive radiochemical assay measurements of MOX fuel from the MALIBU international experimental program. This program included measurements of MOX fuel from a 15 × 15 pressurized water reactor assembly and a 9 × 9 boiling water reactor assembly. The ORIGEN MOX libraries are also compared against detailed assembly calculations from the Phase IV-B numerical MOX fuel burnup credit benchmark coordinated by the Nuclear Energy Agency within the Organization for Economic Cooperation and Development. Finally, the nuclide compositions calculated by ORIGEN using the MOX libraries are shown to be in good agreement with other physics codes and with experimental data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Issue Summary of INL Phase IV Transient Results for IAEA CRP on HTGR UAM Benchmark

This report details the Parallel and Highly Innovative Simulation for Idaho National Laboratory (INL) Code System (PHISICS)/Reactor Excursions and Leak Analysis Program (RELAP5)-3D results obtained for the transient core exercises defined for Phase IV of the International Atomic Energy Agency (IAEA) Coordinated Research Project (CRP) on high-temperature gas cooled reactor (HTGR) uncertainty analysis in modeling (UAM). The Phase III models and results are linked to the earlier Standardized Computer Analyses for Licensing Evaluation (SCALE)/Sampler/New ESC-based Weighting Transport (NEWT) data generated for the lattice physics (lattice) stage Phase I of the CRP. The focus of this report is the Uncertainty/Sensitivity Assessment (U/SA) of the prismatic modular high-temperature gas cooled reactor (MHTGR)-350 design, and specifically for Exercises IV-1 and IV-2 of the benchmark: the Control Rod Withdrawal (CRW) and Pressurised Loss of Cooling (PLOFC) events. The statistical U/SA methodology is implemented and demonstrated using the RAVEN code, based on perturbed cross-section libraries obtained from the SCALE/Sampler sequence. Uncertainties in nuclear data (cross-sections and the average number of neutrons produced per fission, 235U[¯v ]) lead to standard deviations (uncertainties of one s) of approximately 0.5% in the core eigenvalues of the MHTGR-350 and core models. For the coupled neutronics/thermal fluid model, local power density uncertainties up to 3.6% were observed in the colder regions of the core, while the local maximum fuel temperature uncertainties reached 1.5% for the models that included thermal fluid uncertainties. The addition of thermal fluid uncertainties dominated the impacts of nuclear data uncertainties in all cases. The main contributors to uncertainties in the power density and fuel temperatures during the transients were uncertainties in the reactor operating conditions (total power, inlet mass flow rate and inlet gas temperature). Variations in the bypass flows did not have significant impact on any of the output variables. For the nuclear data uncertainties it was found that the 235U(¯v ) / 235U(¯v ) covariance produced the largest sensitivities in terms of its impact on the eigenvalue and peak reactor power. It was also observed that the impact of any nuclear data uncertainties on the maximum fuel temperature was much less significant that the impact on eigenvalue and power. Another important finding was that although the use of eight or more energy groups is recommended for best-estimate HTGR simulation, two-group models produced acceptable uncertainty and sensitivity results for most FOMs. Since the statistical U/SA methodology is computationally expensive, and most transient solver requirements will scale directly with the number of energy groups, two energy groups could be used by HTGR developers during the early stages of design when larger uncertainty margins can be tolerated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Solar Thermochemical Fuel Production via a Novel Low Pressure, Magnetically Stabilized, Non-volatile Iron Oxide Looping Process (University of Florida)

The UF Solar Fuels team has developed a modular reactor at the 10kW scale for the thermochemical conversion of H 2 O and CO 2 to H 2 and CO, respectively, using concentrated solar energy to drive high-temperature redox reactions. Various materials, including ferrites, doped and mixed metal ferrites, and ceria, were investigated and reported on. Ceria was chosen as the reactive material, with thermogravimetric and bench-scale analysis concurring with literature review and pointing towards reasonable kinetics and stability at temperatures on the order of 1500°C. A combined radiation, conduction, and species reaction/mass transport model was developed, utilizing Monte Carlo ray tracing, Lattice Boltzmann, and random walk particle tracking for the respective components to direct development and optimization of the reactor. With experimental data used for reaction rates, the model showed a path to efficiencies >20%, which could be market competitive with current PV-hydrolysis hydrogen generation systems. Economic analysis of a larger-scale plant (100kW modular system with cost performed at 500 modules) has been performed, assuming 0.5% and 10% internal rate of return. Without compression costs, H 2 production cost using this is technology at 12 USD/kg H 2 and 17 USD/kg H 2 for realized efficiencies of 20% and 15%, respectively.

14 SOLAR ENERGY↗