Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Solid phase processing science”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Influence of physical and chemical hydrology on bioremediation of a U-contaminated aquifer informed by reactive transport modeling incorporating 238 U/ 235 U ratios

Microbially-catalyzed reductive immobilization of aqueous uranium (U) as a solid phase has been proposed as a U remediation technique. Both laboratory and field experiments have demonstrated that this reduction reaction alters the 238 U/ 235 U ratio, producing a 238 U-enriched U(IV) solid. In contrast, other major U reactive transport processes fractionate these isotopes much less. This suggests the potential to quantify the extent of bioreduction occurring in groundwater containing U using the 238 U/ 235 U ratio, which would substantially improve upon current practices largely relying on U concentration measurements alone. Current reactive transport models for uranium dynamics include only concentration measurements, which are strongly influenced by highly coupled reactive transformation and aqueous transport processes. The complex physical and chemical behavior of U potentially compromises such quantitative analysis of U storage and release. Here we report the first numerical reactive transport model which explicitly incorporates variations in the 238 U/ 235 U ratio of U and demonstrates improved interpretation of the principal chemical reactions and groundwater transport processes affecting the subsurface mobility and distribution of this widespread contaminant. Recent U bioreduction studies performed in a contaminated aquifer in Rifle, Colorado, USA applied Rayleigh distillation models to interpret U stable isotope fractionation observed as a result of acetate amendment. These simplified models were unable to resolve the spatiotemporal pattern of U isotope fractionation recorded in the aqueous solutes. Here, we employ the multi-component, isotope-enabled CrunchTope reactive transport software to interpret these measured U isotope ratios, and demonstrate accurate reproduction of observed trends in both geochemistry and 238 U/ 235 U ratios for two consecutive years of field experiments. Overall, our results indicate that accurately modeling both the U concentration and isotope ratio distributions greatly constrains the parameter space of the model. We find that the transport properties of U in the Rifle aquifer are governed by the presence of low-permeability regions, which the isotopes are uniquely sensitive to. When U reduction is spatially constrained by these low-permeability regions, the shift in the 238 U/ 235 U ratio becomes more muted. Accurate modeling of observed U isotope ratios thus provides a powerful means to better understand bioremediation, and the current study serves to advance the application of this novel method.

238U/235U↗
DOE OSTI2022

High Performance SOFCs with a Superior Stability for Reliable and Durable Power Systems

Next generation of fuel cells, electrolyzers, and batteries requires higher power, faster kinetics, and larger energy density, which necessitate the use of compositionally complex oxides to achieve multifunctionalities and activity. These compositionally complex oxides may change their phases and structures during an electrochemical process – a so-called “electrochemically driven phase transformation”. The origin for such a phase change has remained obscure. More importantly, there is a need to develop high performance solid oxide fuel cells with an enhanced stability. In this work, the La 0.6 Sr 0.4 Co 0.2 Fe 0.8 O 3 (LSCF) cathode surface is modified by infiltration of Pr 6 O 11 and the power density at 0.8V and 750 °C is improved by 21%. In addition, by replacing the traditional barrier layer Gd 0.2 Ce 0.8 O 1.9 with mixed conducting Pr 0.1 Gd 0.1 Ce 0.8 O 1.9 , the power density increases by 38%. The different mechanism of promotions was investigated by electrochemical impedance spectroscopy. The ohmic resistance is dramatically reduced by applying the PGCO interlayer, and the distribution of relaxation time was used to analyze the mechanism for which the polarization resistance was decreased attributing to the mixed conduction nature in PrO x . An increase of power density at 0.8 V of 0.358 W/cm 2 (71%) is achieved with the implementation of both surface modification and buffer layer engineering. An experimental study and a theoretical analysis were then carried out on phase evolution in praseodymium nickelates. Nickelate-based electrodes show up to 60× greater phase transformation during operation when compared to thermally annealed ones. Theoretical analysis suggests that the presence of a reduced oxygen partial pressure at the interface between the oxygen electrode and the electrolyte is the origin for the phase change in an oxygen electrode. Guided by the theory, an addition of the electronic conduction in the interface layer leads to the significant suppression of phase change, while improving cell performance and performance stability. When an oxygen electrode is under polarization, the oxygen partial pressure at the interface between the oxygen electrode and the electrolyte is lower than that of incoming oxidant. Under a high polarization, the environment at the aforementioned interface may lead to phase transformation of the oxygen electrode. The local oxygen partial pressure is determined by the transport properties at the interfaces. An addition of the electronic conduction in the interface layer, for instance using (Pr,Gd)-doped ceria to replace Gd-doped ceria, results in improved cell performance and performance stability, while the phase transformation is significantly suppressed. This work provides a fundamental understanding of the origin for phase transformation in oxygen electrodes during operation and use this knowledge to develop a high-performance electrode that exhibits improved performance stability.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2024

Recovery of Oxpure 612C-50 Coconut Carbon Particles Using Bump Arrays

Cyanide is used to leach gold from crushed ore (solid matrix) into a gold cyanide solution. The gold is extracted from the cyanide solution by adsorption onto activated carbon. The gold extraction process occurs when the activated coconut carbon is placed into tanks that contain the gold cyanide solution either in a batch process or into a continuous flow circuit. The coconut carbon is then removed from the solution for gold recovery. The use of the mesofluidic separation system (a deterministic lateral displacement system) can significantly reduce costs and waste generated from the process of removing the coconut carbon from the solution. In the coconut carbon recovery process used at many gold mines in Nevada, the gold is attached to carbon particles that are still in the cyanide liquor. The mesofluidic system can rapidly remove the gold bearing carbon particles from the cyanide liquor quickly and with no operating costs. The benefits include: • Reduce operational costs related to filtration and hydro-cyclones • Reduced acid usage for the elution gold stripping phase as only 25% of the fluids will follow the carbon particles to the final express lane when two mesofluidic systems are used in series This alternative particle removal techniques would be advantageous. A promising technique for removing larger particles from slurries is mesofluidic separation, similar to deterministic lateral displacement arrays or “bump” arrays [1] but operates at much larger flow rates. As described by Pease, et al. [2], “Bump arrays in deterministic lateral displacement devices separate large particles from small particles using arrays of staggered posts. Large particles, defined as those with radii larger than the distance between the edge of a post and the stagnation streamline from the next downstream post, must bump toward one side of the device, whereas particles smaller than this distance slalom from entrance to exit without net lateral displacement.” Unlike filters or sieves, the posts in the arrays that cause separation do not block or occlude particle flow but work because particles go around the posts. Unlike hydrocyclones separation, separation is not driven by particle density, because gravitational forces are unimportant to mesofluidic separation. Burns, et al., and Pease, et al., have shown that particles may be separated from complex suspensions, under turbulent flow conditions, and at industrially important flow rates [3-5]. They have also recently shown that mesofluidic devices may be arranged in series to increase separation performance [6]. In this paper, we evaluate the mesofluidic system for the separation of commercial OxPure GR 612 Charred Coconut shell particles as a proof of principle test for the rapid separation using the 1500 micron cut mesofluidic separator. We first describe the experimental system and conditions. We then present the experimental results. [1] Huang, L.R., E.C. Cox, R.H. Austin, J.C. Sturm. 2004. Continuous particle separation through deterministic lateral displacement. Science, 304, 987–990. https://doi.org/10.1126/science.1094567. [2] Pease L.F., J.A. Bamberger, C.A. Burns, and M.J. Minette. 2021. Large Particle Separation from Non-Newtonian Slurries using Bump Arrays. In Proceedings of the ASME 2021 Fluids Engineering Division Summer Meeting FEDSM2021-65904, V003T08A023. New York, New York: ASME. doi:10.1115/FEDSM2021-65904 [3] Burns C.A., T.G. Veldman, J. Serkowski, R.C. Daniel, X.-Y. Yu, M.J. Minette, L.F. Pease. 2021. Mesofluidic separation versus dead-end filtration. Separation and Purification Technology, 254, 117256. https://doi.org/10.1016/j.seppur.2020.117256. [4] Pease L.F., J.E. Serkowski, T.G. Veldman, J. Williams, X.-Y. Yu, M.J. Minette, J.A. Bamberger, C.A. Burns. 2021. Can Bump Arrays Separate Particles from Turbulent Flows?. In Proceedings of the ASME 2021 Fluids Engineering Division Summer Meeting FEDSM2021-67696, V003T08A024. New York, New York: ASME. doi:10.11

mesofluidic separation, slurry, bump array↗
DOE OSTI2025

Multi-physics Topology OPtimization and Additive Manufacturing for High-temperature Heat Exchangers

This research significantly advances the understanding of high-temperature heat exchanger design through an integrated approach that combines topology optimization (TO), triply periodic minimal surface (TPMS) structures, additive manufacturing (AM) and thermohydraulic testing. Each of these components contributes uniquely to a unified, high-performance design, fabrication and testing workflow. Topology optimization serves as the foundation of the design methodology by providing a systematic way to determine the most effective material layout for separating hot and cold fluids while maximizing thermal performance. The researchers introduced a novel three-material optimization framework using two density fields to represent hot fluid, cold fluid, and solid domains. This approach enables automated discovery of optimal shapes and flow paths that cannot be intuitively designed, especially under constraints imposed by manufacturing technologies. Furthermore, constraints such as minimal wall thickness and overhang angles were embedded into the optimization process, ensuring that resulting designs are not only thermally efficient but also manufacturable using modern additive techniques. In parallel, the study delves into the use of Gyroid-based TPMS geometries for constructing the core of the heat exchanger. TPMS structures are known for their high surface area, excellent fluid mixing capabilities, and minimal pressure drop characteristics. The researchers applied a data-driven modeling framework using Heteroscedastic Sparse Gaussian Process Regression (HSGPR) combined with genetic algorithms. This allowed for the rapid evaluation and optimization of key geometric parameters such as frequency, iso-value, and phase shift. The result was a set of Gyroid structures tailored for high heat transfer and low flow resistance, demonstrating clear improvements over conventional straight-channel designs. After the designing process, additive manufacturing played a critical role by turning these highly complex, optimized geometries into physical components. Utilizing Laser Powder Bed Fusion (LPBF) with Haynes 282, the study demonstrated the feasibility of fabricating these heat exchangers at high precision. Post-processing methods, including dilation-erosion operations, were applied to ensure local features adhered to self-supporting constraints. The fabricated structures were then subjected to thermohydraulic testing under conditions representative of supercritical CO 2 Brayton cycles, validating the predicted performance and confirming the viability of the full design-to-fabrication pipeline. Finally, thermohydraulic testing across the above studies served as a crucial experimental validation of advanced heat exchanger. Under consistent high-temperature and high-pressure conditions using supercritical CO 2 , the testing demonstrated that both TO and Gyroid-based TPMS designs significantly outperformed conventional straight-channel HXs. The TO design achieved a 115% increase in UA and NTU and a 27.6% boost in gravimetric power density, while the data-driven optimized Gyroid design delivered a 166% increase in UA and NTU and improved effectiveness from 68.7% to 86.1%. These results validate the simulation models, confirm the manufacturability of complex geometries under AM constraints, and provide key insights into design-performance trade-offs, thereby advancing the development of high-efficiency, compact heat exchangers for extreme environments.

36 MATERIALS SCIENCE↗
DOE OSTI2022

High-Temperature 2D Optical Relaxation Visualizes Enhanced Oxygen Exchange Kinetics at Metal-Mixed Conducting Oxide Interfaces

Solid-state heterointerfaces are of interest for emergent local behavior that is distinct from either bulk parent compound. One technologically relevant example is the case of mixed ionic/electronic conductor (MIEC)–metal interfaces, which play an important role in electrochemistry. Metal–MIEC composite electrodes can demonstrate improved catalytic activity vs single-phase MIECs, improving fuel cell efficiency. Similarly, MIEC surface reaction kinetics are often evaluated using techniques that place metal current collectors in contact with the surface under evaluation, potentially altering the response vs the native surface. Here, techniques enabling direct and local in situ observation of the behavior at and around such heterointerfaces are needed. Here, we develop a spatially resolved optical transmission relaxation (2D-OTR) method providing continuous evaluation of local, high-temperature, controlled atmosphere defect kinetics across a ~1 cm2 sample area simultaneously in a contact-free manner. We apply it to observe the spatial variance of oxygen incorporation and evolution rates at ~525–620 °C, in response to step changes in oxygen partial pressure, on MIEC SrTi 0.65 Fe 0.35 O 3–x films as a function of distance from porous Pt and Au layers. Using this model geometry, we find significant enhancements in kinetics adjacent to the metals that decay over a few millimeter distance. To extract kinetic parameters, we fit the short-term optical data (initial portion of relaxations) with an exponential decay function appropriate for surface-exchange-limited kinetics, yielding apparent surface exchange coefficients (k chem ) with spatial resolution, decreasing with distance from the metal. To understand the kinetic processes governing the complete (long-term) optical relaxations, we performed COMSOL simulations, which demonstrated that a combination of laterally varying k chem and in-plane diffusion controls the observed kinetics over the full time range. Further support for spatially varying kchem comes from demonstrations of changing surface and bulk chemistry vs distance from the metal–MIEC interface, by X-ray photoelectron and optical absorption spectroscopies, respectively. Although microporous Pt and Au are not excellent electrodes in isolation, both metals exert a synergistic effect on the oxygen surface exchange rate in the presence of the mixed conducting film.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Development of Corrosion- and Erosion-Resistant Coatings for Advanced Ultra-Supercritical Materials

This final report summarized the research efforts and major findings of the Phase I Project “Development of Corrosion- and Erosion-Resistant Coatings for Advanced Ultra-Supercritical Materials”, for the period of October 1, 2019 – Sept. 30, 2021. This project is a collaborative endeavor between Tennessee Tech University, Purdue University, Oak Ridge National Laboratory, Siemens Corporation, and Eastern Plating, LLC, aiming at improving the durability and lifetime of high-pressure (HP) steam turbine blades in advanced ultra-supercritical (A-USC) coal-fired power plants through the development of corrosion/erosion-resistant coatings manufactured via a low-cost electrolytic codeposition process. While Tribaloy alloy T-400C was identified by the U.S. A-USC Materials Consortium as a promising coating composition, further composition optimization is needed to enhance its corrosion and erosion resistance for protecting the A-USC Ni-base turbine components. An integrated computational and experimental approach was employed to optimize coating composition/microstructure and processing parameters. In order to identify candidate coating compositions that could offer balanced properties, thermodynamic calculations were performed to explore the γ+Laves composition space in the Co-Ni-Cr-Mo-Si system with different alloying additions at 600-800 °C. Guided by the calculation results, experimental assessment of selected alloys led to the development of a new generation of Tribaloy compositions with the optimal levels of Cr, Mo and Si, reactive element (e.g., 0.4-0.6 wt.% Y) and other alloying additions. The low-cost and non-line-of-sight electro-codeposition process was employed to deposit a Ni(Co)-CrMoSiY composite coating on commercial Haynes 282 (H282) Ni-base alloy. A diffusion treatment was subsequently applied to convert the composite to the Tribaloy-type coating. Both the codeposition parameters and heat treatment conditions were varied to achieve the desired coating composition, microstructure and phase constituents. In addition, since additive manufacturing (AM) may be an alternative cost-saving option for potential A-USC turbine repair, laser direct deposition was explored to fabricate the H282 alloy with minimal defects. The electro-codeposited Tribaloy coating was also applied to the AM H282 substrate to demonstrate the viability of the coating process in improving the surface finish of AM alloys. Both high-temperature oxidation performance and solid particle erosion (SPE) resistance of model alloys and electro-codeposited coatings were evaluated. About 25 model alloys with various Cr/Mo ratios and reactive element levels, as well as partial substitution of Mo with Nb were evaluated with regard to their oxidation resistance in both air and pure steam at 760-800°C. Compositions based on Ni(Co)-20Cr-18Mo-2.6Si-0.6Y (wt.%) showed significantly improved oxidation resistance over the baseline T400-C. Based on the alloy development results, three generations of new Tribaloy coatings (Gen-1, Gen-2, and Gen-3) with various Mo/Cr contents and Y levels were prepared via electro-codeposition and their microstructure/performance were evaluated. Outstanding air and steam oxidation resistance was achieved for the Gen-2 and Gen-3 coatings. Furthermore, while the SPE resistance of the coatings depended on many factors such as temperature, environment, erodent, velocity and impact angle, the coated samples exhibited similar or better SPE resistance compared to the H282 substrate when magnetite was used as erodent (which is a realistic erodent in A-USC steam turbines). The new Tribaloy coatings also had good long-term compatibility with the H282 alloy substrate. The two large-sized rotating barrels were designed, constructed, and employed to coat dummy HP blades. Uniform coating thickness and microstructure were achieved at various blade locations. Also, a preliminary techno-economic analysis of the proposed coating process was conducted to quantify the cost-effectiveness and to assess the commercial viability of the corrosion- and erosion-resistant coatings. It is estimated that a cost reduction of ~30% could be achieved with the electro-codeposition coating process over the state-of-the-art high velocity oxygen fuel (HVOF) thermal spray. Compared to the leading Tribaloy coating technologies such as HVOF and plasma transfer alloying, electro-codeposition based process has advantages such as low-cost process equipment, uniform deposition even for complex shapes, low levels of contaminants/porosities, reduction of powder waste, and potentially better surface finish and longer turbine service life. The Phase 1 study has demonstrated that it is feasible to develop an electro-codeposited Tribaloy coating with balanced corrosion and erosion properties, even though additional research efforts such as further coating process scale-up and longer-term performance evaluation under realistic A-USC conditions are clearly needed.

20 FOSSIL-FUELED POWER PLANTS↗
NASA NTRSDate not supplied

Venus Global Reference Atmospheric Model (Venus-GRAM) Upgrades

Introduction: The Venus Global Reference Atmospheric Model (Venus-GRAM) is one of the most widely used engineering models of Venus’ atmosphere. The Venus-GRAM upgrades are being developed by NASA Marshall Space Flight Center (MSFC) and NASA Langley Research Center (LaRC). This presentation will provide details regarding the upgrades that have been made to Venus-GRAM and the ongoing objectives, tasks, and milestones related to the GRAM upgrades funded by the NASA Science Mission Directorate (SMD). Venus-GRAM: Venus-GRAM is an engineering-oriented atmospheric model that estimates mean values and statistical variations of the atmospheric properties of Venus. Venus-GRAM provides mean values and variability for any point in the atmosphere as well as seasonal, geographic, and altitude variations. Venus-GRAM outputs include atmospheric density, temperature, pressure, winds, and chemical composition along a user-defined path. It is extensively used by the engineering community because of its ability to create realistic dispersions. GRAMs have been integrated into high fidelity flight dynamic simulations of launch, entry, descent and landing (EDL), aerobraking and aerocapture. GRAMs are currently available for Earth, Mars, Venus, Neptune, Titan, and Uranus. The lower atmosphere model in Venus-GRAM (up to 250 km) is based on the Venus International Reference Atmosphere (VIRA) [1]. The Venus-GRAM thermosphere (250 to 1000 km) is based on a MSFC-developed model [2] which assumes an isothermal temperature profile initialized using VIRA conditions at 250 km [3]. The VIRA version included in Venus-GRAM includes Pioneer Venus Orbiter and Probe data as well as Venera probe data, but it does not include a solid planet model or a high-resolution gravity model [4]. Venus-GRAM Upgrade Status: Code Modernization. Venus-GRAM has been rearchitected from Fortran to a common object-oriented C++ framework called the GRAM Suite. This new architecture creates a common GRAM library of data models and utilities. The first C++ release of the rearchitected Venus-GRAM is a straight conversion from the latest Fortran version. Model Upgrades. The focus of the model upgrade task is to improve the atmosphere models in the existing GRAMs and to establish a foundation for developing GRAMs for additional destinations. The GRAM ephemeris has been upgraded to the NASA Navigation and Ancillary Information Facility (NAIF) SPICE toolkit (version N0066). The calculation of the speed of sound has also been improved in the GRAMs. In FY20, the GRAM project established contracts to improve the model data within Venus-GRAM. Hampton University is developing an empirical global model for Venus. The University of Wisconsin is reanalyzing the Venus Express radio occultation observations and analyzing the Akatsuki thermal imaging data. Upgraded Venus-GRAM Release. GRAM Suite Version 1.3 will be released in September 2021 and will contain the rearchitected Venus-GRAM, including the common GRAM framework and planet–specific code. A User Guide and Programmer’s Manual are released with all GRAMs. Conclusions: GRAMs are frequently used toolsets and vital in assessing effects of atmospheres on interplanetary spacecraft during the program life cycle process. Releases of the GRAM Suite, upgrades of the existing planetary GRAMs, and development of new planetary GRAMs are ongoing. Venus-GRAM atmosphere model upgrades will be included in the next phase of GRAM tasks. References: [1] Kliore, A. J. et al. (1985) ASR, 5, 11, 1-304. [2] Justh, H. L. et al. (2006) AIAA/AAS Astrodynamics Specialist Conference & Exhibit, Abstract AIAA-2006-6394. [3] Guide to Reference and Standard Atmosphere Models, BSR/AIAA G-003-2010. [4] Limaye, S. S. (2012), LPSC VEXAG Townhall Meeting. Acknowledgments: The authors gratefully acknowledge support from the NASA SMD.

atmospheric models↗
DOE OSTI2020

BWR Spent Nuclear Fuel Acquisition and Testing to Support DOE-NE High Burnup Spent Fuel Data Project

The Office of Spent Fuel and Waste Disposition (SFWD) within the US Department of Energy (DOE) Office of Nuclear Energy (NE) established the Spent Fuel and Waste Science and Technology (SFWST) campaign to conduct research and development (R&D) activities related to the storage, transportation, and disposal of spent nuclear fuel (SNF) and high-level radioactive waste. The SFWST program was created within SFWD to address issues of extended or long-term SNF storage and transportation. Some near-term objectives of SFWST are to use a science-based, engineering-driven approach to: Support the enhancement of the technical bases to support the continued safe and secure dry storage of SNF for extended periods; Support the enhancement of the technical bases for retrieving SNF after extended dry storage; Support the enhancement of the technical bases for transporting high burnup (HBU) fuel and transporting low burnup fuel and HBU fuel after dry storage DOE-NE, in partnership with the Electric Power Research Institute, developed the High Burnup Spent Fuel Data Project to perform a large-scale demonstration and laboratory-scale testing of HBU pressurized water reactor (PWR) fuels (exceeding 45 gigawatt-days per metric ton of uranium [GWd/MTU]). Under this project, 25 sister rods—which are rods that have the same design, power histories, and other characteristics—were removed from assemblies at the North Anna Nuclear Power Station and sent to Oak Ridge National Laboratory (ORNL) in January 2016. ORNL performed detailed nondestructive examination (NDE) on all 25 rods. The NDE consisted of visual examinations, gamma and neutron scanning, profilometry and rod length measurements, and eddy current examinations. After completing the NDE, 10 of the sister rods were delivered to Pacific Northwest National Laboratory (PNNL) in a NAC International, Inc. legal-weight truck cask in September 2018 for destructive examination (DE). To date, SFWD work has focused on the PWR fuel that is part of the Sister Rod Test program. No boiling water reactor (BWR) fuel has been tested in the program, and the data needs that were identified for the PWR fuel have not been collected for BWR fuel. The goal to obtain six to nine BWR rods and test them at ORNL will support closing this important data gap. BWR fuel comprises approximately 56% of the total fuel assemblies currently in storage at nuclear power plants in the United States. BWR nuclear fuel and cladding designs and manufacturing are significantly different from PWRs. Differences include the following: BWR fuel pellets are larger than PWR pellets; Variations of Zircaloy-2 (including liners) are used instead of the Zircaloy-4 cladding materials used in PWRs; Clad manufacturing and stress-relief processes are different between PWRs and BWRs; The fuel rod dimensions are different because larger rod diameters and thicker cladding are used in BWRs; BWR fuel typically has lower internal rod pressures and sees vastly different operating conditions than PWR fuel (i.e., two-phase flow); BWR assemblies are “canned,” meaning each assembly is surrounded by a metal fuel channel; BWR cladding is often composed of an inner pure Zr liner that has widely different mechanical properties than the Zircaloy-2 alloy and exhibits a stronger affinity for hydrogen; The BWR SNF generally has more total hydrogen in the cladding/liner than typical PWR fuel; The construction of the PWR and BWR assemblies is vastly different; BWR rods are solidly attached to the assembly nozzles and experience a much different vibration and shock load than PWR rods, which are “floating” within a grid system attached to guide tubes, and the rods sit loosely on the bottom end plates. These numerous differences will affect the way the BWR SNF responds under dry storage preparation processes (e.g., vacuum drying) and during transportation. The results collected in the PWR experimental program must be compared with a subset of similar data collected on BWR SNF to establish a technical basis for whether the larger PWR database is sufficient to bound the BWR SNF end-of-life conditions as is currently assumed for several fuel/clad properties. Changes that occur in both fuel types at HBU could exacerbate any mechanical property differences. As the fuel burnup increases, several changes occur that might affect the performance of the fuel, cladding, and assembly hardware in storage and transportation. These changes include increased cladding corrosion layer thickness, increased cladding hydrogen content, increased cladding creep strains, increased fission gas release, and the formation of the HBU structure at the surface of the fuel pellets. The Nuclear Regulatory Commission (NRC) limits the current maximum rod-averaged burnup to 62 GWd/MTU due to these changes and the lack of data at higher burnups.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Task Sharing of Proton Incorporation in Vertically Aligned Nanocomposite Triple Conductors: Growth, Structure, and Surface Exchange Kinetics

As protonic ceramic electrolysis cells emerge for efficient H 2 production, there is a need to develop air electrode materials enabling fast, durable steam splitting and proton incorporation. Single-phase triple conductors may fail to satisfy the myriad performance/stability requirements, and their critical charge-carriers (holes, oxygen vacancies, and protons) are in competition, limiting their concentrations. Instead, we propose task-sharing, vertically aligned nanocomposites (VANs), comprising a proton conductor (BaZr 0.9 Y 0.1 O 3-δ ) and a redox-active mixed ionic electronic conductor (Ce 0.9 Pr 0.1 O 2-δ ), that may enable rapid proton surface exchange at the solid–gas interface and transport along the solid–solid heterointerfaces. We grew VANs by pulsed laser deposition and investigated the interplay between their processing conditions, structure, and proton and oxygen surface exchange kinetics. We varied the substrate temperature, laser repetition rate, laser fluence, and processing oxygen pressure. The crystallinity and phases were characterized by grazing-incidence X-ray diffraction, and the strain and structural order as a function of depth were evaluated by angle-dependent synchrotron X-ray pair distribution function analysis. To evaluate the potential for interdiffusion, the formation energies of substitutional defects were simulated with density functional theory. Corresponding structural analysis and elemental mapping were performed by scanning/transmission electron microscopy, energy-dispersive X-ray spectroscopy, and electron energy-loss spectroscopy, indicating distinct nanoscale compositional regions with a hierarchical structure embedded in individual VANs columns and minimal interdiffusion across a bilayer film. Proton and oxygen surface exchange coefficients (k H , k O ) and polarization resistances were evaluated by electrical and optical relaxations and impedance spectroscopy of VAN-incorporated protonic ceramic electrochemical cells, respectively, at 400–500 °C, demonstrating values comparable to some of the best-known triple and mixed conductors.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Cryogenic Carbon Capture™ (CCC) Status Report

The Cryogenic Carbon Capture™ (CCC) process separates CO2 from light gases in essentially any continuous process. CCC cools the gases to the frost or desublimation point of CO2 (-100 to -135 °C), separates and pressurizes the solids, and warms all streams to produce a CO2-depleted stream at ambient pressure and a pure (99+%) pressurized liquid CO2 stream typically to about 150 bar, both at ambient temperature. The process also recovers all gas moisture and most gas impurities less volatile than CO2 (NOx, SOx, Hg, PM, UHC, CCC, etc.) in separable streams. CCC nearly eliminates refrigeration energy for sensible temperature changes through heat integration. CCC does require energy to change the CO2 phase from a mixed vapor to a pressurized fluid, which represents the minimum energy required of any process for this separation. CCC uses additional energy for turbomachinery inefficiencies, heat losses, moisture removal and overall process pressure drop. Aside from these real-world energy demands, CCC operates near the minimum energy required to perform this gas separation by minimizing stream recycling. CCC compresses CO2 as a liquid, which is one of several reasons it costs about about half as much and consumes about half as much energy as an amine process when using flue gases with about 15% CO2. The process also has several major additional advantages, including (a) it is a bolt-on retrofit technology that does not need steam or any modification of existing equipment, (b) it recovers water and nearly all pollutants in addition to CO2 from the flue gas, (c) it enables highly efficient and cost effective energy storage at grid scale and on time scales of minutes, (d) it enables NG storage if the energy storage option is used, and (d) it has a small footprint and is minimally disruptive to existing plants, requiring only electrical power and a gas source to operate. Sustainable Energy Solutions (SES) has scaled this technology through several levels, the largest of which captures nominally 1 tonne of CO2/day and is called the skid system. Skid system field tests include utility-scale power plants, cement plants, heating plants, and other utility or industrial sites that burn natural gas, biomass, coal, shredded tires, municipal waste, and combinations of these fuels. These field tests produced 95-99% CO2 capture with CO2 purities of 99+% and initial CO2 contents that range from 4 to 28%. SES currently seeks to scale the system to merchant scale (10-80 tonnes of CO2 per day). In the process of doing so, SES has demonstrated the potential for CCC to contribute to energy storage and direct air capture in innovative and cost-effective ways. This presentation discusses the overall process and highlights results from field and in-house tests. These include (a) measured CO2 capture rates and operating conditions from in-house and field tests, and (b) predicted utility-scale costs and energy demands. This discussion also includes the application of the CCC technology to energy storage and direct air capture.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2025

Rapid Bayesian High Entropy Alloy Designs Fabricated via Wire Arc Additive Manufacturing

Purpose: This project seeks to demonstrate a new high-throughput (rapid) alloy design technique applied to creating new high entropy alloys (HEAs) for extreme environments. High entropy alloys shift the design paradigm from being focused on a single principal element (e.g. nickel-based alloys) to target alloys that include high atomic fractions (X >10%) of multiple elements. These HEA materials can exhibit sluggish diffusion and enhanced corrosion resistance, ideal for potential applications in advanced ultra supercritical (A-USC) steam cycles for power generation. Scope: The addition of multiple elements in high atomic fractions creates an enormous design space that cannot easily be investigated by traditional material design strategies such as designed of experiments (DOE). This project utilizes a Bayesian machine learning algorithm that has been modified to work with calculation of phase diagrams (CALPHAD) software. This Bayesian algorithm reduces manual inputs and increase the likelihood of achieving an optimal solution. Compositional inputs to this algorithm will be assessed using existing material property models for high temperature strength and corrosion resistance. The target for alloy performance will be a 15% (~100 ⁰C) increase in allowable service temperature beyond heat-resistant stainless steels while maintaining or improving alloy cost and corrosion resistance. Haynes 230 was selected as a baseline, which is 57 wt% Ni with 22 wt% Cr 14 wt% W, and 2 wt% Mo as solid solution strengtheners. In addition to rapid design via Bayesian machine learning, the alloys were rapidly fabricated using a multi-wire arc additive manufacturing (mWAAM) technique which allows for precise control of alloy composition and assessing of alloy design “windows” to study composition effects. Build speeds for wire-arc additive processes are among the highest for additive technologies enabling rapid and reliable sample fabrication when compared to conventional methods such as arc button melting. The mWAAM samples will be rapidly characterized via instrumented indentation for room temperature modulus and strength and for elevated temperature strength via hot hardness tests. After being screened with hardness testing, potential alloys will be further evaluated with conventional microscopy techniques including scanning electron microscopy (SEM) and transmission electron microscopy (TEM) to assess agreement with modeling results. The most promising compositions will also be evaluated by printing full sized tensile specimens for mechanical behavior tests at elevated temperatures. Results: Bayesian machine learning of a single performance function was initially used to optimize five performance metrics: 1) single phase stability, 2) yield strength, 3) creep resistance (low diffusion coefficient), 4) freezing range (weldability), and 5) material cost. The single performance function was suboptimal as assumptions had to be made about the results while formulating the optimization. A goal-oriented Bayesian optimization strategy (Hanaoka, 2021) was implemented with CALPHAD for use with the five metrics above. This multi-objective Bayesian optimization (MOBO) enabled the design of NiCrCoFe alloys with V and W additions. A base composition of NiCoCr was selected as Ni provides a stable FCC matrix, Cr aids corrosion/oxidation resistance, and Co is a solid-solutions strengthener that also improves creep by increasing the activation energy. Fe helps reduce diffusion coefficients and cost. Finally, V and W were selected for their reasonable solubility and high atomic misfit to aid in solid solution strengthening. Cracking of the mWAAM specimens was an early issue, and the Easton solidification cracking model (Easton et al., 2014a) was selected for addition to the MOBO function. High performing alloys fabricated by mWAAM included Ni 28 Cr 25 Co 26 Fe 15 V 8 and Ni 62 Cr 18 Co 1 Fe 3 W 15 . It was observed that even after adapting the mWAAM process for W, the W did not fully dissolve. To fully evaluate the Ni 62 Cr 18 Co 1 Fe 3 W 15 composition, a cored wire (80-20 NiCr sheath/powder core) was manufactured and printed via WAAM, and HIP’ing was utilized to homogenize and densify the printed alloy. The V and W alloys produced met metrics 1 (solid solution), 4 (solidification cracking), and 5 (cost). However, an unmodeled mechanism of thermal stress cracking was identified in the WAAM produced materials, perhaps exacerbated by the lack of grain boundary strengthening elements (B, C). Conclusions & Recommendations: A high-throughput (rapid) alloy design technique was applied to designing and manufacturing new high entropy alloys (HEAs) for extreme environments utilizing MOBO and mWAAM. The developed process was rapid and effective in addressing the mechanisms included in the model. The lack of grain boundary strengthening element additions (e.g., B, C) was a simplification that likely produced thermal stress cracking that turned into a large part of the investigation. Additions on the order of 0.005 wt% B and 0.05 wt% C likely would have minimized thermal stress grain boundary cracking. Overall, the high throughput design strategy is promising for rapid design of metrics-driven alloys for advanced ultra supercritical (A-USC) steam cycles for power generation. The MOBO and mWAAM process could be commercialized to accelerate metrics-driven alloy design. In addition, the cored-wire process utilized for scale-up is a promising high-volume process for WAAM alloy development and scale-up.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Scalable Nano-Scaffold SOFC Anode Architecture Enabling Direct Hydrocarbon Utilization

This project is based on WVU’s pending patents, technology and aims to design and modify the internal surfaces of the Ni/YSZ anode from currently commercially viable Solid Oxide Fuel Cells (SOFCs) using the additive manufacturing process of Atomic Layer Deposition (ALD). The surface architecture/scaffold added onto the internal surface of the anode possesses an engineered nanostructure but it features only commonly-used oxide conductors and electro-catalyst materials. The surface layer possesses a minimum thickness of ~2-40 nm and is solely designed to control the surface reforming reactions and to increase catalytic activity. Three-dimensional (3D) nano scaffold architectures with the noble metal nano-catalyst, low-cost bimetallic catalytic alloys, and nano-scale ionic conducting oxide fully compatible with the state-of-the-art Ni/YSZ anode, were applied to the internal surface of the entire porous SOFC anode using ALD. In the present work, the surface scaffold architecture is essentially multi-functional at the nano-scale, facilitated by the multiple heterostructured interfaces. It will significantly enhance the power density and cell durability for direct hydrocarbon utilization by (1) increasing the number of electrochemical reaction sites to enhance the hydrogen/hydrocarbon oxidation reactions; (2) reducing carbon formation; (3) mitigating the coarsening of backbone Ni phase and the oxidation attack of Ni from oxidants (e.g., H2O, CO2); and (4) promoting the internal reforming capabilities, especially for natural gas applications. ALD is employed to generate stable anode surface architectures that are uniform, precisely controllable at the atomic scale, and accurately repeatable for processing. The engineered anode surface nano-scaffold architecture was cataloged and analyzed using High-resolution Transmission Electron Microscopy (TEM), and cell power/durability performance assessed via comprehensive electrochemical performance testing with commercial specimens and relevant environments using hydrocarbon fuels. To the best of our knowledge, this project is the First Report on ALD of Ni/YSZ. The actual achievement of this Project includes (1). Successful demonstration of 7 types of ALD layers on Ni/YSZ anode, including Co, Ni, Mn, Pt, Ru, ZrOx and multi-functional nano-composite. (1). Conformal coating and subsequently spontaneously pinning the discrete nano-catalyst, including the precious metal nano-catalyst and the Ni and Co catalysts, on the YSZ surface upon the electrochemical operation in the reduced atmosphere. Those nano-catalysts on the ionic-conducting YSZ provided excellent sites for promoting internal reforming; (2). Demonstrated ALD coating increased both catalytic activity and conductivity of Ni/YSZ. Conformal coating provided dopants and introduced additional electrical conducting pathways on the YSZ ionic conductor. The doped surface layer of YSZ with mixed conductivity thus further introduces the active triple phase boundaries adjacent to the ALD-coated nano-catalysts such as Pt, Co, and Ni that are pinned on the YSZ surface. The nano-composite ALD coating on Ni/YSZ anode has significantly increased cell durability; and (3). ALD coating of Ni/YSZ anode increased the power density of the entire cell by 300%. For a long time, the SOFC performance, such as the power density, was deemed hindered by the cathode. The sluggish oxygen reduction reaction (ORR) in the cathode was deemed as hindering the power density of the SOFCs. For the anode-supported commercial SOFCs, the cell performance is considered to be limited by the cathode's performance. For the first time in the field of SOFC, this project has demonstrated that (1). the performance of commercial SOFCs can be further increased by the ALD coating on Ni/YSZ anode backbone. (2). ALD coating on Ni/YSZ fuel electrodes results in the enhancement of power density, and increased reliability, robustness, and endurance of SOFCs, for their application using both hydrogen and hydrocarbon fuels over the entire operating temperature range of 650-800ºC for the inherently functional commercial cells. (3). ALD coating provides alternative approaches of exsolutions for introducing the stable catalyst onto the internal surface of the Ni/YSZ electrode. ALD coating could be much more versatile than exsolution in employing the catalysts with various chemistries onto the various backbones. (4). Due to the negligible amount of ALD materials coated onto the internal surface of the porous cathode of the as-fabricated cells, a peak power density increase up to 300 % induced by ALD coating was simultaneously achieved in terms of both power density and specific power. (5). The ALD coating developed through this project was applied to both the SOFC and Solid Oxide Electrolysis Cells (SOEC). SOEC’s face a similar but more demanding need to improve the fuel electrode's performance. It opens further research directions for electrocatalytic surface nanoionics with a wide range of chemistry. It will revolutionize our ability to render the formation of a nanostructured electrode that has been constantly pursued yet barely achieved for practical SOFC/SOEC applications. The research is also immediately transformative since both the preliminary data and the proposed work are on the direct implantation of nanoionics into the state-of-the-art inherently functional SOCs. It represents an immediate impact on the commercial sectors in SOC technology since the applied ALD processing is computer-controlled ALD coating using the commercial ALD systems, and it is scalable to both the single cells and SOC stacks.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Thermodynamically consistent incorporation of the Langmuir adsorption model into compressible fluctuating hydrodynamics

For a gas–solid interfacial system where chemical species undergo reversible adsorption, we develop a mesoscopic stochastic modeling method that simulates both gas-phase hydrodynamics and surface coverage dynamics by coupling the Langmuir adsorption model with compressible fluctuating hydrodynamics. To this end, we derive a thermodynamically consistent mass–energy update scheme that accounts for how the mass and energy variables in the gas and surface subsystems should be updated according to the changes in the number of molecules of each species in each subsystem due to adsorption and desorption events. By performing a stochastic analysis for the ideal Langmuir model and the full hydrodynamic system, we analytically confirm that our mass–energy update scheme captures thermodynamic equilibrium predicted by equilibrium statistical mechanics. We find that an internal energy correction term is needed, which is attributed to the difference in the mean kinetic energy of gas molecules colliding with the surface from that computed from the Maxwell–Boltzmann distribution. By performing an equilibrium simulation study for an ideal gas mixture of CO and Ar, with CO undergoing reversible adsorption, we validate our overall simulation method and implementation.

Adsorption↗
DOE OSTI2020

Vadose Zone Model for SALDS for Composite Analysis

The objectives of the vadose modeling for the updated Hanford Site composite analysis (CA) are to simulate the flow and transport of water and radionuclide releases from the surface to the water table and to provide radionuclide transfer rates for the plateau to river (P2R) model, version 8.3 (CP-57037, Model Package Report: Plateau to River Groundwater Model, Version 8.3). Water additions include natural recharge and water discharged to the ground as a result of industrial processes associated with Hanford Site operations. Contaminant sources include radionuclides in water discharged to the ground during operations and radionuclides disposed “dry” in solid waste burial grounds or other means. The following 16 radionuclides were selected for this modeling effort: carbon-14 (C-14), chlorine-36 (Cl-36), tritium (H-3), iodine-129 (I-129), neptunium-237 (Np-237), rhenium-187 (Re-187), strontium-90 (Sr-90), technetium-99 (Tc-99), uranium-232 (U-232), uranium-233 (U-233), uranium-234 (U-234), uranium-235 (U-235), uranium-236 (U-236), uranium-238 (U-238), radium-226 (Ra-226), and thorium-230 (Th-230). The simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070. The parallel version of the Subsurface Transport Over Multiple Phases (STOMP1) simulator, officially named the exascale Subsurface Transport Over Multiple Phases (eSTOMP) is used to simulate flow and transport for the vadose models. The documentation for the STOMP code is comprehensive. The theoretical and numerical approaches applied in the STOMP code are documented in a published theory guide (PNNL-12030, STOMP Subsurface Transport Over Multiple Phases Version 2.0 Theory Guide). The code has undergone a rigorous verification procedure against analytical solutions, laboratory-scale experiments, and field-scale demonstrations. The application guide (PNNL-11216, STOMP Subsurface Transport Over Multiple Phases Application Guide) provides instructive examples in the application of the code to classical groundwater problems. The user’s guide (PNNL-15782, STOMP: Subsurface Transport Over Multiple Phases Version 4.0: User’s Guide) describes the general use, input file formatting, compilation, and execution of the code. The primary output of the vadose zone modeling is radionuclide transfer rates to the groundwater for input into the P2R model. The rates will be summed over the 100 by 100 m P2R grid cells that fall within the vadose zone model source domain. The Hanford Site Central Plateau was subdivided into 26 individual vadose zone models, with 13 in the 200 East Area and 13 in the 200 West Area. Waste sites that have a completed performance assessment (PA) or past-leak analysis were not included as sources of radionuclides. Instead the vadose zone to groundwater transfer rates of the Environmental Restoration Disposal Facility, Integrated Disposal Facility, US Ecology, and Waste Management Area C (WMA C) PAs and the past-leak analysis for WMA C were used as direct input to the P2R model. Each of the vadose zone models is documented in separate environmental calculation files (ECFs). This ECF describes the State-Approved Land Disposal Site (SALDS) model. The scope of this ECF is to document the development and results of the SALDS vadose zone model. CP-63515, Model Package Report: Central Plateau Vadose Zone Models, describes the approach, assumptions, process of determining the number of models required and domain of each model, input data, and processing common to all the models.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Meshfree simulation and prediction of recrystallized grain size in friction stir processed 316L stainless steel

Friction stir processing (FSP) is a promising solid-phase microstructural modification technique that can repair and enhance damaged stainless steel surfaces exposed to harsh environments. The quality of the repaired material is closely correlated to the recrystallized grain size in the stir zone (SZ), which is influenced by the thermomechanical conditions dictated by FSP process parameters. Thus, establishing a reliable relationship between these parameters and recrystallized grain size in the SZ is crucial for optimizing repair quality. However, existing experimental approaches often rely on indirect temperatures measured far from the SZ, along with rough strain rate estimations, which are imprecise and time-consuming. Meanwhile, existing mesh-based modeling methods usually face numerical challenges when dealing with the large material deformations inherent in FSP. Here, to address these issues, this study introduces a meshfree process model for FSP based on the smoothed particle hydrodynamics (SPH) method, aimed at predicting process conditions under different parameters. The model is validated using experimental data from 11 combinations of tool traverse and rotation speeds on 316 L stainless steel. Correlations between process parameters, material flow, temperature, strain, strain rate, and recrystallized grain size are revealed through SPH simulations and electron backscatter diffraction (EBSD) imaging. The results show that in situ SZ temperatures range from 1071 to 1322°C, which exceed the tool temperature by over 300°C. Furthermore, SZ temperature, strain rate, and grain size increase monotonically with higher tool temperature and faster traverse speed. A relationship is then established between the model-predicted Zener-Hollomon parameter and the recrystallized grain size based on EBSD data, expressed as ln(d) = -0.364 ln(Z) + 14.673. Finally, this relationship exhibits satisfactory accuracy with errors of less than 26.9% in predicting grain sizes at various SZ locations, which offers valuable insights for optimizing FSP repair processes for 316 L stainless steel.

316L stainless steel↗
DOE OSTI2021

Consistent Controls on Trace Metal Micronutrient Speciation in Wetland Soils and Stream Sediments

Trace metal are essential for microbially-mediated biogeochemical processes occurring in anoxic wetland soils and stream bed sediments, but low availability of these elements may inhibit anaerobic element cycling and transformations. Solid-phase speciation is likely a critical control on trace metal availability but has seen limited study in anoxic systems having concentrations similar to geological background levels, where metal limitations may be most prevalent. We have investigated trace metal concentrations and solid-phase speciation in three freshwater subsurface aquatic systems: marsh wetland soils, riparian wetland soils, and the sediments of a streambed. These systems displayed low solid-phase trace metal concentrations, generally at or below geological background levels, which generally followed the trend Zn > Cu ≈ Ni > Co and showed no correlation with major element compositions. All soils and sediments were dominated by quartz but varied in clay mineralogy as well as the organic matter, total sulfur, and total iron contents. X-ray absorption near-edge structure (XANES) spectroscopy shows that sulfur speciation in both wetlands is dominated by organic sulfur. Elemental sulfur and iron sulfides together made up <25% of the sulfur in the wetland soils, but the distribution between inorganic and organic forms was reversed in the stream sediments. Ferrous and ferric iron in clay minerals were common species identified by both XANES and extended X-ray absorption fine structure (EXAFS) spectroscopies at all sites. Iron(III) oxides were substantial components in all but the marsh wetland soils. Quantitative analysis of copper, nickel, and zinc XANES spectra revealed similar metal speciation across all sites. Copper speciation was dominated by sulfides, adsorbed species, and minor amounts of copper bound to organic matter; no metallic copper was detected. Nickel speciation also varied little and was dominated by nickel in clay mineral octahedral sheets and nickel sulfide, with adsorbed species also present. Zinc speciation was slightly more varied, with the marsh wetland soils and stream bed sediments containing adsorbed species, zinc associated with clay mineral structures, and zinc bound to reduced sulfur groups on organic matter, whereas the riparian wetland soils lacked clay-associated zinc but contained zinc sulfide. Trace metals bound to reduced sulfur occurred at every site, with a greater sulfur-bound fraction for copper. The fractional abundance of sulfur-bound species showed no relationship with soil or sediment total sulfur content, which varied by two orders of magnitude. More broadly, the observations in this study suggest that trace metal speciation in freshwater wetland soils and stream sediments is consistently dominated by a small set of recurring components which are distinct for each metal. Furthermore, this may represent a general geochemical phenomenon in anoxic soils and sediments containing trace metals at background concentrations (as low as 3 µg g -1 ) that was not predicted from systems that are contaminated with or naturally-enriched in copper, nickel, or zinc.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Design Study for Bio-Oil Production from Biomass Using a Dual Fluidized-Bed Reactor

To evaluate the re-design and reconfiguration of a dual-fluidized bed (DFB) gasification system into a recirculating pyrolysis reactor, Computation Fluid Dynamic (CFD) simulations of the system were conducted. The Barracuda Virtual Reactor® computational particle fluid dynamic code was used to perform simulations of the pyrolysis process. Modeling of the chemical reaction kinetics for both gas phase and solid particle phase were included. The recirculating pyrolysis reactor shown in Fig. 1a is based on a bubbling-bed biomass pyrolyzer and a riser combustor to convert the remaining char. The operational differences in the re-configuration of the DFB gasification system into a recirculating pyrolysis system for the production of bio-oil are (1) replacement of a low-surface area inert bed material with a high-surface-area bed material that has acidic properties to provide catalytic activity for the production of bio-oil with reduced oxygen content, (2) lower temperature and residence time for bio-oil production from pyrolysis, (3) replacement of the fluidization gas in the bubbling bed pyrolyzer from steam to nitrogen, and (4) the reduction of pyrolyzer freeboard volume. The bed material used for catalytic pyrolysis is Sasol 300 (300-micron dia., bulk density 0.94 kg/l, and surface area 130 m2/g) and is a theta-alumina with mild acidity. This is in comparison with previous standard bed material Carbo HSP (430-micron dia., bulk density 2.01 kg/l, and surface area 0.03 m2/g) used for gasification. For bio-oil production, pyrolysis in the bubbling bed requires temperatures in the range of 550 C in comparison with gasification temperatures near 850 C. To attain this lower temperature requires management of the energy mass balances, with control of the bed material recirculation rate between bubbling bed pyrolyzer and riser combustor, the introduction of a nitrogen purge in the pyrolyzer, and adjusting the pressure balance between the two vessels. To extract bio-oil from the pyrolysis reactor with a snorkel, two different freeboard configurations were evaluated. In Fig. 1 b the existing high freeboard configuration is shown and in Fig. 1 c the reduced freeboard design is presented. The introduction of a nitrogen purge for the high free board configuration provided the highest bio-oil production from the CFD simulations. To decrease the bed-material circulation rate, primary and secondary air on the combustor side were reduced, and the pressure on the combustor was slightly increased. A portion of the biomass and bio-oil was observed to be transported to the combustor, leading to a smaller pyrolysis yield. The control of the pyrolyzer temperature is performed by controlling the circulation rate. Good fluidization of the bubbling bed and cascade PID control are required to keep the temperature from oscillating due to large time delay experienced when changing primary and secondary air. (a) (b) (c) Figure 1. (a) Dual Fluidized Bed pyrolysis system configuration, (b) mole fraction of nonpolar biooil in high freeboard configuration, and (c) mole fraction of nonpolar bio-oil in low freeboard configuration with nitrogen purge introduced in both configurations.

BASIC BIOLOGICAL SCIENCES,BIOMASS FUELS↗
DOE OSTI2021

Machine learning to discover mineral trapping signatures due to CO 2 injection

Mineral trapping is pursued as a geological CO 2 sequestration (GCS) mechanism because it permanently stores CO 2 in solid phases or minerals. However, CO 2 mineral-trapping mechanisms are poorly understood due to (1) lack of sufficient field and laboratory data characterizing these complex processes, and (2) challenges to develop site-specific reactive-transport models coupling fluid flow and geochemical reactions occurring at various temporal (from milliseconds to years) and spatial (from pore (millimeters) to field (kilometers)) scales. Reactive transport with additional complexities such as heterogeneity can make the simulation outputs even more difficult to interpret because of complex nonlinearity and multi-scale interdependencies. Furthermore, the values of model outputs such as concentrations can vary by several orders of magnitude, making it harder to correlate and characterize the impact of the variables via traditional data interpretation techniques such as exploratory data analyses. Recently, machine learning (ML) has shown promise in feature discovery and in highlighting hidden mechanisms that cannot be obtained by existing data-analytics and statistical methods. In this study, we applied an unsupervised ML approach, non-negative matrix factorization with custom -means clustering (NMF) to the data generated by reactive-transport simulations of GCS. The reactive-transport data consisted of 19 attributes, including four physio-chemical variables (pH, porosity, aqueous CO 2 , and sequestered CO 2 ), six chemical species (K + , Na + , HCO, Ca 2+ , Mg 2+ , Fe 2+ ), and four carbonate minerals (calcite, dolomite, siderite, and ankerite), a feldspar mineral (albite), and four clay minerals (illite, clinochlore, kaolinite, and smectite) over a period of 200 years of simulation time. Furthermore, the simulation data used was for Morrow B sandstone at the Farnsworth hydrocarbon unit in Texas. Data are sampled at two locations within the model domain: (1) at the injection well and (2) 200 m west of the injection well. The injection was performed for a period of 10 years. Using NMF, we estimated the temporal interdependencies among the 19 attributes over a span of 200 years. We found that NMF was able to identify four reaction stages and their dominant attributes; these cannot be directly discerned through traditional visualization (e.g., line plots, Pareto analysis, Glyph-based visualization methods) or exploratory data analysis tools of the simulation data. The four stages were: reactions in the injection phase followed by short-, mid-, and long-term reactions. The NMF analysis also revealed that 10 among the 19 attributes are dominant. These dominant attributes for mineral trapping include calcite, dolomite at injection well, siderite at 200 m away from the injection well, clinochlore, kaolinite, Na + , K + , Ca 2+ , Mg 2+ , pH, and aqeuous CO 2 . Finally, at late times (65–200 years), our results showed that calcite plays a major role in mineral trapping with insignificant contribution from siderite, ankerite, and clay minerals. These findings make the proposed unsupervised ML-model attractive for reactive-transport sensing towards real-time GCS monitoring.

54 ENVIRONMENTAL SCIENCES↗