Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fuel elements”

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.

1,187 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Enhancement of Carbon Capture Reactor Performance (Final Technical Report)

Significant challenges are still present in post-combustion CO 2 capture and new technologies and advanced components are needed to significantly advance the deployment of CO 2 capture for natural gas combined cycle (NGCC) plants. Critical elements of CO 2 capture that still need to be addressed include how to increase CO 2 mass transfer in the absorber column with liquid to gas ratios of <1.2, while reducing the size of the absorber column to reduce capital costs. Research involving chemical mechanism with design, synthesis, and assembly of materials with targeted functionally were combined with advanced additive manufacturing techniques towards development of enhanced CO 2 capture reactors that can lead to safe, reliable, and low-cost carbon capture technologies. The objective of the project was to develop and test novel carbon capture materials and reactor components that contribute to increased CO 2 mass transfer through increased turbulent gas-liquid interface and improved solvent wetting within the absorber. A technoeconomic analysis (TEA) was completed showing how the proposed technology decreases capital costs by reducing the size of the absorber column and the amount of packing required for high CO 2 capture rates. A technology maturation plan (TMP) was also developed to describe the current technology readiness levels (TRL) and outline additional research and development (R&D) needed to further develop these advanced components for NGCC CO 2 capture plants. The successful completion of this project has shown a pathway to reduce the absorber size and associated construction costs of post-combustion NGCC CO 2 capture systems at 97% capture and promote the utilization of abundant natural gas for production of reliable electricity.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2023

Theory and Simulation of Metal–Insulator–Semiconductor (MIS) Photoelectrodes

A metal–insulator–semiconductor (MIS) structure is an attractive photoelectrode-catalyst architecture for promoting photoelectrochemical reactions, such as the formation of H 2 by proton reduction. The metal catalyzes the generation of H 2 using electrons generated by photon absorption and charge separation in the semiconductor. The insulator layer between the metal and the semiconductor protects the latter element from photo-corrosion and, also, significantly impacts the photovoltage at the metal surface. Understanding how the insulator layer determines the photovoltage and what properties lead to high photovoltages is critical to the development of MIS structures for solar-to-chemical energy conversion. Herein, we present a continuum model for charge-carrier transport from the semiconductor to the metal with an emphasis on mechanisms of charge transport across the insulator. The polarization curves and photovoltages predicted by this model for a Pt/HfO 2 /p-Si MIS structure at different HfO 2 thicknesses agree well with experimentally measured data. The simulations reveal how insulator properties (i.e., thickness and band structure) affect band bending near the semiconductor/insulator interface and how tuning them can lead to operation closer to the maximally attainable photovoltage, the flat-band potential. This phenomenon is understood by considering the change in tunneling resistance with insulator properties. The model shows that the best MIS performance is attained with highly symmetric semiconductor/insulator band offsets (e.g., BeO, MgO, SiO 2 , HfO 2 , or ZrO 2 deposited on Si) and a low to moderate insulator thickness (e.g., between 0.8 and 1.5 nm). Beyond 1.5 nm, the density of filled interfacial trap sites is high and significantly limits the photovoltage and the solar-to-chemical conversion rate. These conclusions are true for photocathodes and photoanodes. This understanding provides critical insight into the phenomena enhancing and limiting photoelectrode performance and how this phenomenon is influenced by insulator properties. The study gives guidance toward the development of next-generation insulators for MIS structures that achieve high performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Energy-Storing Cryogenic Carbon Capture™ for Utility and Industrial-scale Processes (Final Report)

This project leverages work completed under funding from previous U.S. Department of Energy (DOE) projects and other sources of funding around the Cryogenic Carbon Capture Process™ (CCC). The purpose of this project is to further develop a perturbation of the CCC process that stores energy called the Cryogenic Carbon Capture™ Energy Storing (CCC-ES) process. The successful completion of the tasks in this project has prepared the CCC-ES process for scaling and developed initial plans and a feasibility analysis for an engineering-scale or small-commercial scale system that can store at least 10 MWh of energy. The specific areas of work in this project are: (1) Develop a conceptual, site-specific integrated process flow diagram: Sustainable Energy Solutions (SES) developed a conceptual, site-specific plan including a process flow diagram for integrating CCC-ES at a specific location. The proposed system was based on a typical engineering-scale CCC system designed to capture 30-50 metric tonnes of CO 2 per day (TPD), with modifications and additional changes and equipment put in place to add the energy storage element to the system. The 30-50 TPD scale corresponds to the required 10 MWh of stored energy. This task includes updating simulation software for transient analysis. SES also coordinated with the Jim Bridger Power Plant in Point of Rocks, Wyoming to get site-specific data. Using these data, SES completed site-specific preliminary process flow diagrams, estimated capital and operating costs for the CCC-ES system, and used the upgraded software package to estimate the primary figures of merit. (2) Complete a technoeconomic analysis (TEA): SES completed a full-scale TEA of the CCC-ES system to meet specific success criteria. This TEA used process information from the host site selected above. (3) Technology and data gap assessment for the CCC-ES process: As part of the technology gap assessment, SES performed an assessment of the best alternative energy storage solutions including figures of merit and operational strengths and weaknesses. SES then compared the proposed system with a discussion of how it addresses the weaknesses of alternative systems, showing for example a cost of stored energy at less than $50/MWh and with a roundtrip efficiency of greater than 95%. This assessment includes technical and other risk analyses and technology gaps with research and development needs to commercialize the process by 2030. Finally, work includes a commercialization roadmap/development pathway outlining the major milestones to move to a commercially operating CCC-ES system. (4) Phase II pre-FEED project plan: SES completed a project plan for the Phase II pre-FEED analysis using the same host site as in the Phase I conceptual study.

03 NATURAL GAS↗
DOE OSTI2026

Advected carbon younger than the sediment fuels microbial metabolism in a pumped deep aquifer

Ever deeper wells are drilled worldwide to pump potable groundwater. Recent studies argue that overpumping compresses clays and releases reactive dissolved organic carbon (DOC), which in turn drives a series of microbial reactions that affect groundwater potability including arsenic concentrations. Here, we use a novel method to measure the radiocarbon ages of microbial RNA to determine the source of reactive or metabolizable DOC and argue against the compression of clays as the sole source of carbon. We show that microbial RNA (5,230; 5,550; 6,250 yr; n = 3 wells), DOC (280-10,800 yr; n = 13), and methane (modern to 6,240 yr; n = 3), from an overpumped deep aquifer in Bangladesh are much younger than the overlying clay layers deposited during the Pleistocene over 12,000 years ago. Mass-balance indicates that at least half of the carbon incorporated into RNA has to come from reactive DOC or methane that is advected downward via vertical recharge. This metabolism of advected organic carbon could have implications for the quality of water pumped from deep aquifers.

Biological and medical sciences↗
DOE OSTI2023

A diffusion–deformation model with damage for polymer undergoing rapid decompression failure

Green hydrogen is emerging as one of the most promising alternatives to replace fossil fuels. While hydrogen gas has a good energy density by weight, its poor energy density by volume requires it to be stored under high pressure for commercial use. The hydrogen infrastructure developed to handle this high pressure hydrogen contains multiple components made of polymers as sealing agents. Although polymers do not react chemically with hydrogen gas, they undergo a mechanical failure when high pressure hydrogen gas is suddenly depressurized. This phenomenon known as 'rapid decompression failure' occurs due to the diffusion of hydrogen through polymer and getting trapped inside pre-existing cavities or voids. Here, in this paper, a continuum mechanics-based coupled diffusion-deformation-damage model was developed to predict the hydrogen distribution, stress distribution, and damage propagation inside the polymer while it undergoes rapid decompression failure. The polymer was modeled as a hyperelastic material because it represents the nonlinear material response observed in uniaxial tensile tests perfectly. The effects of hydrogen diffusivity, pre-existing cavity size, cavity location, applied hydrogen pressure, and depressurization rate on damage initiation were studied. It was found that the coefficient of diffusion plays an important role in damage initiation and damage was mostly concentrated in the inside areas rather than near the surface. Experiments were conducted with EPDM polymer which agreed well with the predicted trends using the given model. The effect of adding carbon black and silica filler particles and plasticizer to the pure EPDM polymer was also studied. It was found that damage during RDF decreases with the addition of fillers, but increases with the addition of the plasticizer. Finally, the damage evolution in the presence of two cavities was also studied, and was found that the interaction of stress fields around the cavities alters the damage occurring during RDF.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Computational Fluid Dynamics Analysis of the Molten Salt Tritium Transport Experiment Test Section

Tritium, a radionuclide produced through neutron capture by lithium and other elements (beryllium and fluoride) in molten salts, presents unique challenges to radionuclide release. This is true for both fusion energy breeder blankets and molten salt fission reactors. The fundamental understanding of tritium transport is crucial to the safe design and operation of these reactors. The Molten Salt Tritium Transport Experiment (MSTTE), currently under construction at Idaho National Laboratory, aims to investigate tritium transport phenomena using a forced-convection fluoride salt loop. This loop is designed to study various transport mechanisms, such as permeation through metals and gas-liquid interactions, and is intended to support future research on tritium extraction units. A critical aspect of the MSTTE loop design is ensuring a fully developed velocity profile before the fluid reaches the permeation test section where measurements are made. This study employs computational fluid dynamics to model the salt flow behavior within the MSTTE permeation test section. A realizable k-ε turbulent model with enhanced wall treatment is used to simulate the single-phase, vertical upward flow of molten salt FLiNaK under isothermal conditions. The simulation results indicated flow distortion and underdeveloped profiles at all planned flow rates within the test section due to the 85-deg sharp bend. To address this issue, a reduced diameter with a reducer and expander and a flow conditioner are investigated to achieve fully developed flow. The analysis showed that the flow conditioner successfully corrected the flow profile, achieving fully developed behavior at a flow rate of 50 liters per minute (LPM). In conclusion, this research enhances our understanding of flow dynamics in molten salt systems and contributes to optimizing tritium transport control technologies.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Development and Validation of SAM Multi-dimensional Flow Model for Thermal Mixing and Stratification Modeling

Thermal mixing and stratification in large pools or enclosures are very important phenomena that are critical to nuclear reactor safety. Because of the wide ranges of time and length scales associated with such phenomena, accurate modeling and simulation of thermal mixing and stratification remain as the key unresolved, challenging problems for reactor transient analyses. In traditional system analysis codes, simplified zero-dimensional (0-D) models are widely used for their high numerical efficiency, but they generally suffer from very limited prediction accuracies or range of applicability. Like traditional system analysis codes, the current version of the SAM code has implemented such simplified 0-D and one-dimensional mixing models. On the other hand, high-resolution Computational Fluid Dynamics (CFD) tools are often used to model complex thermal mixing and stratification phenomena. They are, however, generally numerically expensive, and they require large amounts of computational resources. It is therefore desirable to implement advanced and efficient thermal mixing and stratification modeling capabilities embedded in a system analysis code. This approach will improve the accuracy of reactor safety analyses when thermal mixing and stratification are involved, and also avoid using the large computational resources needed for high-resolution CFD analysis. Currently, with the support of the U.S. DOE Office of Nuclear Energy’s Nuclear Energy Advanced Modeling and Simulation program, an effort has been launched to develop and implement a multi-dimensional flow model in the system analysis code SAM, and demonstrate its applications to model thermal mixing and stratification phenomena in large enclosures. The main outcomes of this research and development activity are summarized in this report, which presents an attempt to include a built-in advanced multi-dimensional flow model in a system analysis code with the focus on overcoming the simulation challenges of thermal mixing and stratification phenomena. In this report, we start with the introduction of existing SAM code capabilities to simulate thermal mixing and stratification phenomena, which is followed by a short summary of the multi-dimensional model implemented in the SAM code, including both the physical model and the Finite Element Method code implementation. In this study, two options were implemented in the SAM code to model turbulent flows: a relatively simple built-in turbulence model and an interface to accept externally computed turbulent viscosities (e.g., from a high- resolution CFD simulation). Code validation studies on this newly added capability were then carried out to compare SAM simulation results with experimental data from the SUPERCAVNA facility, which was designed to study the complex flow recirculation and thermal stratification phenomena relevant to sodium fast reactor designs. In this study, one transient and two steady-state test cases were used for code validation. Different approaches have been used to model the complex turbulence flow fields in the SUPERCAVNA facility. A highly simplified zero-equation turbulence model was first used, but it was determined that it is too simple to capture the complex turbulence flow fields in these test cases. Subsequently, the code validation continued with the use of turbulent viscosity data from high-resolution STAR-CCM+ CFD simulations to improve the accuracy of the results. Using this approach, the SAM simulation results showed very good agreement with both the SUPERCAVNA experimental data and STAR-CCM+ simulation results. In this report, we demonstrate the development, implementation, and successful validation of a multi-dimensional flow model in the SAM code, which aims to improve the simulation accuracy for complex thermal mixing and stratification phenomena. Lessons have also been learned, including that in cases where the flow fields are not well predicted by the zero-equation model, the iinput of turbulent viscosities from an external source can enhance the overall predictive capabilities needed to accurately capture complex thermal-hydraulic phenomena. Therefore, future research will be needed to further improve the code’s capabilities, such as by developing a more efficient and robust approach to capture the turbulence effects in the SAM code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Nuclear Materials Packaging, Transportation, and Systems Analysis Group Software Quality Assurance Plan: ANSYS Mechanical Finite Element Analysis Software Version 2023R1

ANSYS Inc. develops and markets engineering simulation software and services used in the aerospace, automotive, manufacturing, electronics, biomedical, energy, defense, and many other industries. ANSYS is dedicated to engineering simulation and is the world’s leading software provider. ANSYS was founded in 1970 and is headquartered in Canonsburg, Pennsylvania. ANSYS provides an engineering analysis tool combining structural, thermal, computational fluid dynamics, acoustic, and electromagnetic simulation capabilities. ANSYS has two main programs, which use the same solvers: (1) Mechanical APDL (ANSYS Design Parametric Language), a Fortran-based coding platform, and (2) ANSYS Workbench, which uses a graphical user interface to aid in finite element analysis implementation. This plan covers both APDL and Workbench. The ANSYS computer program is a large-scale, multipurpose finite element program that can be used to solve several classes of engineering analyses. The analysis capabilities of ANSYS include the ability to solve static and dynamic structural analyses, steady-state and transient heat transfer problems, mode-frequency and buckling eigenvalue problems, static or time-varying magnetic analyses, and various types of field and coupled-field applications. The program contains many special features that allow nonlinearities or secondary effects such as plasticity, large strain, hyperelasticity, creep, swelling, large deflections, contact, stress stiffening, temperature dependency, material anisotropy, and radiation to be included in the solution. As ANSYS has been developed, other special capabilities such as substructuring, submodeling, random vibration, kinetostatics, kinetodynamics, free convection fluid analysis, acoustics, magnetics, piezoelectrics, coupled-field analysis, and design optimization have been added to the program. These capabilities contribute further to making ANSYS a multipurpose analysis tool for varied engineering disciplines. The ANSYS program has been in commercial use for over 50 years and has been used extensively in the aerospace, automotive, construction, electronic, energy services, manufacturing, nuclear, plastics, oil, and steel industries. Additionally, many consulting firms and hundreds of universities have used ANSYS for analysis, research, and educational purposes. ANSYS is recognized worldwide as one of the most widely used and capable programs of its type. Ansys design analysis software is the first created within a quality system with ISO 9001 certification, the internationally accepted quality standard. Product development, testing, maintenance and support processes also meet the United States Nuclear Regulatory Commission's quality requirements, as they have for nearly four decades. The Quality Assurance Service Agreement is suitable for the customers working in the nuclear industry who need to meet specific federal regulations including 10CRF50 Appendix B and provisions of 10CFR21. ANSYS has retained its original International Organization for Standardization (ISO) 9001 accreditation certificate since1995-05-04, It’s current certificate is valid until 2027-05-29.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

2023 Critical Materials Strategy

The global effort to curb carbon emissions is accelerating demand for clean energy technologies and the materials they rely on. Demand for these materials will only continue to grow, especially as some nations aim to achieve net zero emissions by 2050. While some major materials like steel, copper, and aluminum are already powering the fossil fuel economy, others are more minor materials with potential supply risks. These risks could jeopardize the ability to reduce greenhouse gas emissions within the desirable timeframe to avoid significant climate change. In some cases, it may be necessary to take action to improve the resilience of material supply chains and mitigate supply risks. Understanding the importance of individual materials to clean energy and the supply risks associated with them is necessary to identify which materials may serve as potential roadblocks to a clean energy future. The U.S. Department of Energy (DOE) issued a series of 13 supply chain deep dive assessment reports on various energy technologies in 2022 in response to President Biden’s Executive Order on America’s Supply Chains (E.O. 14017). These reports emphasized that supply chain bottlenecks can occur at any stage of the value chain from mining and refining to component and even sub-system manufacturing. The bottlenecks are a combination of factors such as material availability, equipment availability, work force availability and quality, logistics, regulatory framework, and market conditions. These bottlenecks were worsened during the global Covid-19 pandemic. Its lingering impacts have hindered capacity expansion for material supply chains and prevented product lead-time recovery. One approach to reduce supply chain risks for the United States is to have a strong domestic manufacturing sector with a diverse set of producers. Boosting responsible domestic production would require leveraging the latest science not only in material extraction but also in developing substitutes, recycling, reuse, and remanufacturing. This report is an updated analysis of previous Critical Materials Strategy (CMS) reports published by the DOE in 2010, 2011, and 2019 based on national and global priorities, technology advancement, and technology adoption trends. Like the CMS reports, this analysis presents the results of a formal material criticality assessment to identify which materials are critical to the continued deployment of clean energy technologies globally. The analysis in this report leveraged the DOE supply chain deep dive assessments to develop the initial list of materials to evaluate. This DOE Critical Materials Assessment (CMA) is conducted independently of criticality assessments performed by other U.S. government agencies, such as that conducted by the U.S. Geological Survey (USGS). This analysis complements the USGS critical minerals determination in three aspects. First, the DOE assessment is performed from a global perspective, while the USGS analysis focusses on the importance of minerals to the U.S. economy. Second, this report focuses on the importance of materials to clean energy technologies, rather than to the economy in general. Lastly, this study is forward looking to 2035 based on clean energy deployment scenarios, whereas the USGS assessment is retrospective. Materials evaluated in this report that do not appear in the USGS Critical Minerals List include copper, uranium, electrical steel, and SiC. A draft version of this report received ~80 public comments related to supporting data and methodological improvement. Those comments have been incorporated as much as possible where appropriate. Highlights of findings from this 2023 CMA include: Rare earth materials (neodymium, praseodymium, dysprosium, and terbium) used in magnets in electric vehicle (EV) motors and wind turbine generators continue to be critical. While dysprosium (Dy) and terbium (Tb) are both heavy rare earth elements that serve the same function in magnets, the criticality of Tb is slightly lower than that for Dy in the short term due to the widespread use of Dy in high-grade magnets and Tb’s present role as a substitute. Similarly, praseodymium (Pr) is critical in the medium term but only near critical in the short term because it is more substitutable in magnets than neodymium (Nd); Materials used in batteries for EVs and stationary storage are now considered to be critical. While cobalt (Co) was found to be critical in this and previous reports, lithium (Li) becomes critical in the medium term due to its broader use in various battery chemistries and the rampant growth of the EV industry. Natural graphite is a new addition in this assessment and is also found to be critical; Platinum group metals used in hydrogen electrolyzers, such as platinum (Pr) and iridium (Ir), are critical due to an increased focus on hydrogen technologies to achieve net zero carbon emissions, while those used in catalytic converters, such as rhodium (Rh) and palladium (Pd), were screened out due to the decreased importance of catalytic converters in the medium term; Gallium (Ga) continues to be critical due to its use in light-emitting diodes (LEDs). In addition, the use of Ga has increased in magnet manufacturing and in semiconductor in forms such as gallium arsenide (GaAs) or gallium nitride (GaN); Major materials like Aluminum (Al), copper (Cu), nickel (Ni), and silicon (Si) move from noncritical in the short term to near critical in the medium term due to their importance in electrification; Electrical steel is near critical due to its use in transformers for the grid and electric motors in EVs.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

FY23 Progress Report on Viscosity and Thermal Conductivity Measurements of Molten Salts

As presented in this report, thermal conductivity and viscosity measurements were performed on key chloride pseudo-binary molten salt systems of relevance to molten salt reactor developers. Thermal conductivity measurements were conducted with a variable gap technique, in which a known heat flux is driven across a molten salt specimen and the temperature difference is measured, allowing calculation of the thermal conductivity. This is achieved by establishing a small gap between the bottom of a cylindrical inner containment, which houses electrical heating elements, and an outer containment, which houses cooling channels; the gap size can be varied by compression of a formed bellows. A new calibration scheme was developed herein, involving a correction factor to the heat flux based on He measurements at various temperatures. Furthermore, the data processing methodology was improved to minimize the impact of radiative heat transfer in the calculation of salt specimen thermal conductivity from the temperature difference measurements. Viscosity measurements were conducted with a rolling ball viscometer, in which a ball rolls some known distance in an angled tubular crucible, and the terminal velocity can be used to calculate the viscosity of the salt. The measurement can be performed in a quartz crucible, with which a standard camera can be used to track the ball, or in a metal crucible, with which x-ray radiography is required to track the ball. A new custom x-ray system was made and dedicated to the rolling ball viscometer to enable high throughput automated measurements with salts which require containment with metal. Both systems have been integrated with a ventilation stack which allows for off-gassing of radioactive material, enabling future measurements with U-bearing salts. The thermal conductivity measurements performed herein were with NaCl-KCl (44 mol% NaCl). This salt system was measured in the previous fiscal year, however the thermal conductivity values obtained were comparatively low, and so the system was remeasured with the aforementioned calibration scheme and improved post-processing techniques. The newly obtained thermal conductivity values for NaCl-KCl (44 mol% NaCl) indicate good agreement with kinetic theory and ab-initio models (within 5–10 %). The viscosity measurements performed herein were with three different compositions of NaCl-KCl: 75, 50, and 25 mol% NaCl; a new calibration scheme was applied to account for variable flow effects in the laminar regime. The results show reasonable agreement with literature (5–20 %, depending on the temperature, composition, and study); however, literature values are likely higher than true values based on pure end-member measurements performed in the comparative studies. The measurements conducted herein do show a trend such that viscosity increases with increasing NaCl concentration, which agrees well with one of two comparative studies. In general, the measurements conducted herein gives confidence in the capability to use these systems to accurately measure thermal conductivity and viscosity of actinide-bearing salts within the next fiscal year.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Comparison of the Chemical Composition of Liquids from the Pyrolysis and Hydrothermal Liquefaction of Lignocellulosic Materials

Major differences in thermal stability and hydrotreatment behavior of HTL and pyrolysis oils have been reported in the literature. However, little is known about the variations in the chemical composition of these oils that could explain such differences. Two commercial wood pyrolysis oils (Pyrovac and BTG), and their water-soluble (WS) and water-insoluble fractions (WIS) were analyzed and compared with the aqueous (WS WD-57 ) and oily (WIS WS-57 ) fractions obtained from hydrothermal liquefaction (HTL) of Douglas-fir. The samples were characterized by GC/MS, Karl Fischer titration, carbonyl content, total acid number, elemental composition, calorific value, proximate analysis, Fourier Transform Infrared Spectroscopy (FTIR), Folin-Ciocalteu (FC), and UV fluorescence. All the fractions were also analyzed by Fourier Transform Ion Cyclotron Resonance Mass Spectroscopy (FT-ICR-MS) and by Electrospray Ionization (EI). The most prevalent class of compounds in the water insoluble phases were phenols derived from lignin. Water-soluble phases contain mostly the oxygenated compounds derived from cellulose and hemicellulose and were richer in carbonyl functional groups. The water content of the resulting aqueous phases were between: 65 (WS BTG ) and 96 (WS Pyrovac ) wt. %. The bio-oil from BTG has higher water content and lower HHV, compared to Pyrovac oil. The GC/MS results of BTG oil show the presence of a more prominent acetic acid peak and higher TAN number than the Pyrovac oil. The GC/MS of Pyrovac oil showed more obvious mono-phenol peaks. The quantification of this family by Folin-Ciocalteu method confirmed higher content of monophenolic compounds compared with the BTG oil. The lower thermal stability of pyrolysis oils compared with HTL biocrudes can be partially explained by the fact that pyrolysis oils (BTG and Pyrovac) contain carbohydrates while HTL biocrude (WIS WD-57 ) doesn’t. Thus, we decided to further investigate the chemical differences between the phenolic rich fractions insoluble in water and the holocellulose derived compounds soluble in water. Even after water extraction, the acid content of the water insoluble fraction from BTG (WISBTG) was higher than the acid content of the water insoluble fraction obtained by HTL (WISWSD-57). Likewise, the acid content of the aqueous phases derived from pyrolysis oils (WS Pyrovac , WS BTG ) was also higher than for the aqueous phase obtained by HTL (WSWSD-57). This result is in part due to the use of bases in the HTL process that neutralizes the acid formed in that process. Moreover, the starting feedstock may also influence the differences between the oils. Although, the UV-Fluorescence spectra, ICR-MS and the EI analyses showed some minor differences in the molecular weight and chemical make-up of the oligomers soluble and insoluble in water from pyrolysis and HTL; the differences observed were not large enough to justify the differences in behavior between these oils reported in the literature. Our results suggest that the differences observed between HTL biocrudes and pyrolysis oils are likely partially due to the presence of holocellulose derived products in the pyrolysis oils and higher acid contents.

09 BIOMASS FUELS↗
DOE OSTI2026

Biomass Demineralization: A Critical Need for Future Biorefineries

Biomass contains up to 14 essential elements that serve as nutrients for plant growth and development, including photosynthesis and enzyme functionalities. These elements in different chemical forms (e.g., minerals) constitute the inorganic fraction of biomass and can cause operational issues in thermal and biochemical biomass conversion technologies. In biomass gasification and pyrolysis processes, for instance, inorganics can cause fouling, tar formation, and corrosion. In catalytic and biochemical processes, inorganics can poison catalysts, alter biochemical pathways, and modify product yields and selectivity. This review provides an overview of the inorganic content in biomass feedstocks, the critical role inorganics play in plant biochemistry, the effect that inorganics have in various thermochemical and biochemical biomass conversion technologies, and different approaches to remove them from biomass. We provide recommendations for future research, focusing on developing technologies to effectively remove inorganics from biomass, using the inorganics to improve soil quality and for alternative applications, and designing biorefineries to convert demineralized biomass obtained from diverse sources.

09 - BIOMASS FUELS↗
DOE OSTI2025

Quantifying Climate Change Effects of Bioenergy and BECCS: Critical Considerations and Guidance on Methodology

Bioenergy is a critical element in many national and international climate change mitigation efforts, including as a carbon dioxide removal strategy combined with the capture and durable geological storage of flue gas emissions (BECCS). However, divergent results on the effectiveness of bioenergy as a climate change mitigation measure are reported in the scientific literature. Climate impacts of bioenergy depend on case-specific factors, primarily biophysical features of the biomass production system, and the design and efficiency of conversion and capture processes. Estimates of climate impacts are also strongly affected by methodological choices and assumptions, and much of the divergence between studies derives from differences in the assumed alternate use of the land or feedstock, the alternate energy source and the system boundaries applied. We present a methodology to support robust estimates of the climate change effects of bioenergy systems, updating the standard methodology developed by the International Energy Agency's Technology Collaboration Program on Bioenergy. We provide guidance on the key choices including the reference land use and energy system that bioenergy is assumed to displace, spatial and temporal system boundaries, co-product handling, climate forcers considered, metrics applied and time horizon of impact assessment. Researchers should consider the whole bioenergy system including all life cycle stages, and choose system boundaries, reference systems and treatment of co-products that are consistent with the intended application of the results. The assessment should be normalised to a functional unit that can be compared with other systems delivering an equivalent quantity of the same function. All significant climate forcers should be included, and climate effects should be quantified using appropriate impact assessment methods that distinguish the impact of time. Consistency in methodology and interpretation will facilitate comparison between studies of different bioenergy systems.

09 BIOMASS FUELS↗
DOE OSTI2021

MFIX DEM Enhancement for Industry-Relevant Flows (Final Report)

The overall goal of this two-phase project is to implement performance improvements of the Multiphase Flow with Interphase Exchanges (MFIX) Discrete Element Model (DEM) code that enable a transformative shift for industrial use. Prior to this effort, the largest simulations performed using MFIX are O(10 7 ) particles. This falls short of the O(10 9 ) particle simulations that must be completed on a timescale of days or weeks (vs. months or years) to enable simulations with physically-relevant domain sizes to be incorporated into industrial design cycles within five years. This was accomplished by tailoring best-in-class practices to bear on the unique challenges posed by the MFIX-DEM algorithm and code base. Scientific simulations (e.g., in cosmology, turbulent combustion) routinely use massively parallel computing to update far more particles in short wall clock times. Results from Phase 1 (1.5 years in duration) indicated significant gains in speed were possible for a wide range of benchmark cases. Moreover, a survey sent to >35 companies indicates that the timing is ideal for such an enhanced tool, with >80% of the respondents indicating that DEM is already value-added or will be within the next 5 years, and >70% of the respondents indicating that improved speed is the top computational priority. In Phase 2 (3.5 years in duration), the two major barriers that hinder industry from effectively using multiphase Computational Fluid Dynamics (CFD) to cut costs and improve performance, namely computational overhead and confidence in predictions, continued to be addressed. Regarding the former, the results from Phase 1 to guide the effort, with enhancements focused on an improved time-stepping algorithm and particle sorting. Four target problems of 1 billion particles each and increasing complexity were identified: homogeneous cooling, tumbler with continuous particle size distribution, discharge from a rectangular hopper and a cylindrical riser. Each of these were successfully simulated for relevant time scales (on order of seconds) using less than 24 hours of wall clock time. These represent the first 1-billion particle DEM simulations performed with MFIX, namely using the MFIX-Exa code. This code is currently under development at NETL in collaboration with Lawrence Berkeley National Laboratory. Regarding the second barrier on predictive uncertainty, experiments from Phase 1 (interacting nozzles - hydrodynamics only) and Phase 2 (very small-scale segregation experiments) were used to demonstrate the ability of two simplified approaches to uncertainty quantification (UQ). By limiting the number of particles, UQ based on the simplified treatment was compared to standard UQ, which was shown to have much higher computational demands. Experiments were also performed on a pilot-scale stripper unit to provide validation data for future CFD-DEM simulations and UQ.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2021

Striving to translate shale physics across ten orders of magnitude: What have we learned?

Shales will play an important role in the successful transition of energy from fossil-based resources to renewables in the coming decades. Aside from being a significant source of low-carbon intensity fuels, like natural gas, they also serve as geologic seals of subsurface formations that may be used to isolate nuclear waste, sequester CO 2 , or store intermittent energy (e.g., solar hydrogen). Despite their importance, shales pose significant engineering and environmental challenges due to their nanoporous structure and extreme heterogeneity that spans at least ~10 orders of magnitude in spatial scale. Two challenges inhibit a system-level understanding: (1) the physics of fluid flow and phase behavior in shales are poorly understood due to the dominant molecular interactions between minerals and fluids under confinement, and (2) the apparent lack of scale separation that prevents a reliable (closed) description of the physics at any single scale of observation. In this review, we focus on the latter issue and discuss scale translation, which in its broadest sense is transforming data or simulations from one spatiotemporal scale to another. While effective scale translation is not exclusive to shales, but all geologic porous media, the need for it is especially acute in shales given their high degree of heterogeneity. Classical theories like homogenization, while indispensable, fail when scales are not separated. Other methods, like numerical upscaling, scale-translate in only one direction: small to large, but not the reverse, called downscaling. However, the confluence of advances in three areas are bringing challenging problems such as shales within reach: increased computational power and scalable algorithms; high-resolution imaging and multi-modal data acquisition; and machine learning to process massive amounts of data. While these advances equip geoscientists with a wide array of experimental and computational tools, no individual tool can probe the entire gamut of heterogeneity in shales. Their effective use, therefore, requires an ability to bridge between various data types obtained at different scales. The aim of this review is to present a coherent account of computational and experimental methods that may be used to achieve just that, i.e., to perform scale translation. We provide a broader definition of scale translation, one that transcends classical homogenization and upscaling methods, but is consistent with them and accommodates notions like downscaling and data translation. After a brief introduction to homogenization, we review hybrid methods, numerical upscaling and its recent extensions, multiscale computing, high-resolution imaging, and machine learning. We place particular emphasis on multiscale computing and propose an algorithmic framework to bridge between the pore (micro) and Darcy (macro) scales. Throughout the paper, we draw comparisons between the various methods and highlight their (often hidden) similarities, differences, benefits, and pitfalls. We finally conclude with two case studies on shales that exemplify some of the methods presented.

58 GEOSCIENCES↗
DOE OSTI2022

Thermocouple Testing in Support of the AGR-5/6/7 Experiment

This report documents thermocouple testing performed in IRC Lab C-15 over a period of seven years. This testing supported selection and characterization of the thermocouple set used in the AGR-5/6/7 experiment. The following summary was taken directly from the report. Temperature measurement is a challenging aspect of very high temperature irradiation experiments because commonly used high-temperature commercial thermocouples such as platinum-rhodium (Types S, R, and B) and tungsten-rhenium (Type C), suffer dramatic drift because of neutron-induced transmutation. As a result, these types of thermocouples, which are used routinely for industrial temperature measurements outside of reactors, are used only in very special circumstances for reactor experiments. Conversely, because of their low neutron cross-sections, Type N thermocouples are affected to only a limited extent by neutron irradiation. However, the use of these nickel-based thermocouples is limited when the temperature exceeds 1050°C due to drift arising from minor alloying elements migrating from the thermocouple's metal sheath to the thermoelements. This change in the composition of the thermo-elements results in significant decalibration of the signal. The issues described above were recognized during the early planning stages of the final AGR experiment (designated AGR-5/6/7), and a thermocouple furnace testing program was performed over a seven-year period (2014-2019, 2021) to first select and then characterize the best thermocouple set for the high temperature regions of the AGR-5/6/7 irradiation experiment. The calculated temperature range of the AGR-5/6/7 experiment was 600–1500°C. For temperatures below 1000°C standard Type N thermocouples were deemed adequate. The furnace testing campaign identified two thermocouple types suitable for measuring temperatures above 1000°C, a Mo/Nb thermocouple developed at INL called HTIR-TC, and a Type N thermocouple developed by Cambridge University (called herein Cambridge Type N), which featured a custom high nickel alloy sheath. One of the original goals of the furnace testing program was to identify a thermocouple capable of low drift operation near the peak temperature expected in AGR-5/6/7, i.e., about 1400°C. The HTIR-TC design appeared promising in this regard, however a manufacturing difficulty proved to be a barrier and instead the furnace testing focused on drift performance at 1250°C. The manufacturing difficulty was that the Nb sheaths of the HTIR-TCs experienced extreme embrittlement when heat treated at 1600°C or greater. Heat treatment is needed to stabilize the emf output of this TC type, and the higher the heat treatment temperature the higher the peak temperature of stable operation. Because of the sheath embrittlement the heat treatment temperature had to be lowered to 1450°C resulting in a stable operating temperature of about 1250°C. One of the successes of the furnace testing program was identification of a shortcoming in the heat treatment procedure that had been traditionally used in the production of HTIR-TCs. The shortcoming was that the entire heated length of the HTIR-TC sensor was not being heat treated, but rather only the part of the sensor expected to experience temperatures above 1000°C. The problem manifested itself when the thermocouples were removed from the heat treat furnace and placed in another furnace with a different geometry, their indicated temperatures would be widely scattered, but mostly in the negative direction. The solution was to heat treat the entire heated length of the sensor. Since the deepest immersion depth in the AGR-5/6/7 experiment was about 40 inches, a heat treatment length of 48 inches was used. After this change was implemented, thermocouples which were moved into a new environment with a different temperature profile (i.e., a different furnace), produced accurate temperature measurements. Although assembly of the AGR-5/6/7 experiment was completed in September of 2017 (and irradiation begun in 2018), furnace testing of thermocouples continued in 2018 and 2019. The main purpose of this testing was to establish very long-term drift characteristics of the HTIR and Cambridge Type N thermocouples installed in the experiment. Representative thermocouples from the same lots as those installed in the AGR-5/6/7 experiment were used. Additionally, thermocouples of different designs, (particularly variations on the HTIR-TC design) were "piggy-backed" on this testing program to provide insights for instrumenting future very high temperature irradiation experiments. This two-year testing program demonstrated that HTIR-TCs and Cambridge Type N TCs could operate at 1250°C for up to 10,000 hrs (and in some cases longer) while experiencing negative drifts on the order of 2-4°C/1000 hrs. This performance was considered acceptable given the extreme operating environment the sensors faced.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Experimental and modeling studies of PR and ND oxalate solubility to high ionic strengths: Insight into actinide(III) oxalates

Actinide oxalates are chemical compounds important to nuclear industry, ranging from actinide separation in waste reprocessing, to production of specialty actinides, and to disposal of high level nuclear waste (HLW) and spent nuclear fuel (SNF). In this study, the solubility constants for Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O by performing solubility experiments in HNO 3 and mixtures of HNO 3 and H 2 C 2 O 4 at 23.0 ± 0.2 °C have been determined. The targeted starting materials, Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O, were successfully synthesized at room temperature using PrCl 3 , NdCl 3 and oxalic acid as the source metrials. Then, we utilized the targeted solubility-controlling phases to conduct solubility measurements. There was no phase change over the entire periods of experiments, demonstrating that Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O were the solubility-controlling phases in our respective experiments. Based on our experimental data, we have developed a thermodynamic model for Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O in the mixtures of HNO 3 and H 2 C 2 O 4 to high ionic strengths. The model for Pr 2 (C 2 O 4 ) 3 •10H 2 O reproduces well the reported experimental data for Pu 2 (C 2 O 4 ) 3 •10H 2 O, which are not utilized for the model development, demonstrating that Pr(III) is an excellent analog for Pu(III). Similarly, the model for Nd 2 (C 2 O 4 ) 3 •10H 2 O reproduces the solubility of Am 2 (C 2 O 4 ) 3 •10H 2 O and Cm 2 (C 2 O 4 ) 3 •10H 2 O. The Pitzer model was used for the calculation of activity coefficients. Based on the published, well established model for dissociation constants for oxalic acid and stability constants for actinide-oxalate complexes [i.e., AmC 2 O 4 + , and Am(C 2 O 4 ) 2 - ] to high ionic strengths, we have obtained the solubility constants (log 10 K 0 ) for the following reactions at 25 °C, Pr 2 (C 2 O 4 ) 3 •10H 2 O ⇌ 2Pr 3+ + 3C 2 O 4 2- + 10H 2 O(l). Nd 2 (C 2 O 4 ) 3 •10H 2 O ⇌ 2Nd 3+ + 3C 2 O 4 2- + 10H 2 O(l). to be -30.82 ± 0.30 (2σ), and - 31.14 ± 0.35 (2σ), respectively. These values for can be directly applied to Pu 2 (C 2 O 4 ) 3 •10H 2 O, Am 2 (C 2 O 4 ) 3 •10H 2 O and Cm 2 (C 2 O 4 ) 3 •10H 2 O. The model established for actinide oxalates by this study provides the needed knowledge with regard to solubilities of actinide/REE oxalates at various ionic strengths, and is expected to find applications in many fields, including the geological disposal of nuclear waste and the mobility of REE under the surface conditions, as Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O can be regarded as the pure Pr and Nd end-members of deveroite, a recently discovered natural REE oxalate with the following stoichiometry, (Ce 1.01 Nd 0.33 La 0.32 Pr 0.11 Y 0.11 Sm 0.01 Pb 0.04 U 0.03 Th 0.01 Ca 0.04 ) 2.01 (C 2 O 4 ) 2.99 •9.99H 2 O. Regarding its importance in the geological disposal of nuclear waste, Am 2 (C 2 O 4 ) 3 •10H 2 O/Pu 2 (C 2 O 4 ) 3 •10H 2 O/Cm 2 (C 2 O 4 ) 3 •10H 2 O can be the source-term phase for actinides, as demonstrated by the instance in the disposal in clay/shale formations. This is exemplified by the stability of Am 2 (C 2 O 4 ) 3 •10H 2 O in comparison with Am(OH) 3 (am), Am(OH) 3 (s) and AmCO 3 (OH)(s) under the relevant geological repository conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗