Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Bulk behavior”

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 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Bioblendstocks to Optimize Mixing Controlled Compression Ignition (MCCI) Engines

In this project, a team of researchers from the University of Massachusetts Lowell, the University of Maine, and Mainstream Engineering developed an integrated process for the product of bioblendstocks to optimize mixing controlled compression ignition (MCCI) engines. The objective was to improve the energy density, sooting propensity, and cetane number of base diesel fuel while maintaining cold weather behavior. The process converts woody biomass (e.g. sawmill residues) into bio-oil through selective fast pyrolysis; the bio-oil is then selectively upgraded to form selectively oxygenated, minimally-branched hydrocarbons using non-noble metal catalysts in combination with metal-catalyzed hydrogenation. Advanced predictive models, in conjunction with existing property databases, and experimental testing are used to evaluate overall bioblendstock properties and their impact on base diesel fuel. An iterative, targeted upgrading approach was implemented to optimize the proposed bioblendstock’s properties. Assessment methodologies included techno-economic analysis, life-cycle assessment, property testing, and engine testing. Ultimately, the project team successfully produced a viable bioblendstock while identifying critical process points related to scale-up efforts. It was found that producing pyrolysis oils at 500 degrees C and with pine particle sizes of 1-2 mm led to bio-oil with a higher yield (of approximately 45 wt%) and rich amounts of aromatic alcohols. The resultant pyrolysis oil was then upgraded using a sequence of mild hydrotreating, followed by catalytic etherification and esterification, followed by another final mild hydrotreating to produce a blendstock containing saturated species with a limited, but non-zero, amount of oxygen. The aromatic alcohols produced by pyrolysis were especially helpful in this regard, as the resulting bicycloethers and derivatives exhibited high cetane numbers. While most bulk properties of the bioblendstock met or exceeded targeted thresholds, viscosity and cloud point notably fell outside the expected range; this could be addressed by blending limits and/or through the use of additives that are commonplace in current refinding practices. Identification of a bioblendstock that can be produced economically at scale while improving the performance and emissions characteristics of internal combustion engines positively affects the economy by boosting domestic fuel production and the environment by decreasing harmful emissions and increasing efficiency.

09 BIOMASS FUELS↗
DOE OSTI2023

Improved melt model for power flow

Accelerators that drive z-pinch experiments transport current densities in excess of 1 MA/cm 2 in order to melt or ionize the target and implode it on axis. These high current densities stress the transmission lines upstream from the target, where rapid electrode heating causes plasma formation, melt, and possibly vaporization. These plasmas negatively impact accelerator efficiency by diverting some portion of the current away from the target, referred to as “current loss”. Simulations that are able to reproduce this behavior may be applied to improving the efficiency of existing accelerators and to designing systems operating at ever higher current densities. The relativistic particle-in-cell code CHICAGO ® is the primary code for modeling power flow on Sandia National Laboratories’ Z accelerator. We report here on new algorithms that incorporate vaporization and melt into the standard power-flow simulation framework. Taking a hybrid approach, the CHICAGO® kinetic/multi-fluid treatment has been expanded to include vaporization while the quasi-neutral equation-of-motion has been updated for melt at high current-densities. For vaporization, a new one-dimensional substrate model provides a more accurate calculation of electrode thermal, mass, and magnetic field diffusion as well as a means of emitting absorbed contaminants and vaporized metal ions. A quasi-fluid model has been implemented expressly to mimic the motion of imploding liners for accurate inductance histories. For melt, a multi-ion Hall-MHD option has been implemented and benchmarked against Alegra MHD. This new model is described with sufficient detail to reproduce these algorithms in any hybrid kinetic code. Physics results from the new code are also presented. A CHICAGO ® Hall-MHD simulation of a radial transmission line demonstrates that Hall physics, not included in Alegra, has no significant impact on the diffusion of electrode material. When surface contaminant desorption is mocked in as a hydrogen surface plasma, both the surface and bulk-material plasmas largely compress under the influence of the j × B force. Similar results are seen in Alegra, which also shows magnetic and material diffusion scaling with peak current. Test vaporization simulations using MagLIF and a power-flow experimental geometry show Fe + ions diffuse only a few hundred µm from the electrodes, so present models of Z power flow remain valid.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Deformation and Transformation Textures in the NaMgF3 Neighborite—Post-Perovskite System

The D″ region of the lower mantle, which lies just above the core–mantle boundary, is distinct from the bulk of the lower mantle in that it exhibits complex seismic heterogeneity and seismic anisotropy. Seismic anisotropy in this region is likely to be largely due to the deformation-induced texture (crystallographic preferred orientation) development of the constituent mineral phases. Thus, seismic anisotropy can provide a marker for deformation processes occurring in this dynamic region of the Earth. Post-perovskite-structured (Mg,Fe)SiO3 is believed to be the dominant mineral phase in many regions of the D”. As such, understanding deformation mechanisms and texture development in post-perovskite is important for the interpretation of observed seismic anisotropy. Here, we report on high-pressure diamond anvil cell deformation experiments on NaMgF3 neighborite (perovskite structure) and post-perovskite. During deformation, neighborite develops a 100 texture, as has been previously observed, both in NaMgF3 and MgSiO3 perovskite. Upon transformation to the post-perovskite phase, an initial texture of {130} at high angles to compression is observed, indicating that the {100} planes of perovskite become the ~{130} planes of post-perovskite. Further compression results in the development of a shoulder towards (001) in the inverse pole figure. Plasticity modeling using the elasto-viscoplastic self-consistent code shows this texture evolution to be most consistent with deformation on (001)[100] with some contribution of glide on (100)[010] and (001)<110> in NaMgF3 post-perovskite. The transformation and deformation mechanisms observed in this study in the NaMgF3 system are consistent with the behavior generally observed in other perovskite–post-perovskite systems, including the MgSiO3 system. This shows that NaMgF3 is a good analog for the mantle bridgmanite and MgSiO3 post-perovskite.

Geochemistry & Geophysics↗
DOE OSTI2020

Influence of Vegetation Canopy Architecture on Dynamics of Spatiotemporal Light Distribution and its Effect on Land-Biosphere-Atmosphere Fluxes of Energy, Carbon, and Water

The dynamic spatiotemporal distribution of light and micro climatic variables (wind speed, atmospheric CO 2 concentration, humidity, pressure, and air temperature) within the plant canopy as determined by the canopy architecture play a critical role in land-biosphere-atmosphere (LBA) photosynthetic carbon uptake, transpiratory water loss and energy partitioning. Vegetation models that include effects of sun and shade flecks, leaf temperature and energy balance partitioning between sensible and latent heat fluxes, and reduced turbulent mixing in dense vegetation canopies can accurately predict canopy performance, particularly its shaded leaves, which can constitute 80% of leaf area and contribute up to 50% of carbon update and 65% of water loss. However, current vegetation models with their bulk characterization of shaded leaves poorly represent the complex three dimensional (3D) canopy architecture are insufficient to capture the spatiotemporal dynamics of shaded leaves. Here, we propose to employ an x-ray computed tomography (XCT) imaging technique in combination with an automated, high throughput, image processing pipeline (3DArch, to be developed) to obtain an explicit 3D triangulated surface mesh of a plant canopy (Fig 1). 3DArch would provide a much needed technological jump in 3D LBA modeling where computer generated artificial 3D models suffer from lack of high resolution detailed canopy architecture parameters. The 3D canopy surface mesh output from 3DArch is coupled with an explicit 3D forward ray tracing model (fastTracer) to obtain a high resolution representation of the dynamic spatiotemporal light environment within the canopy over the course of a day (Fig 1). The results from the light attenuation model are then coupled with a process based multi-scale bio-physical canopy and leaf model of photosynthesis (PlantGro) to get fluxes of carbon, water and energy at the leaf and canopy scale, in addition to estimating whole plant carbon uptake, and water loss. Our proposed model framework would be the first of its kind to provide a breakthrough capability in the field of process-based plant modeling that bridges the key technological gap in shaded leaf modeling through explicit representation of 3D canopy architecture. Computer generated artificial 3D model canopy simulations on grasses show that the canopy light distribution does not follow an exponential behavior as assumed in most 1D vegetation models, but rather exhibits a heterogeneous behavior, with a small but significant portion of interior canopy leaves experiencing peak light levels due to sun flecks while the rest of the canopy receives substantially lower light levels. These results indicate a bias in traditional models that overestimates shaded leaf contribution towards whole canopy photosynthesis by up to 15% and underestimate transpiration by up to 12% with an associated increase in sensible heat fluxes since photosynthesis, transpiration and energy partitioning are tightly coupled processes. However, these test results need to be validated using the 3D architecture obtained from a real plant imaged in 3D using 3DArch, rather than an artificial model plant. Our proposed model framework would allow us to identify inefficiencies in canopy light distribution and prescribe optimal canopy architectures (leaf optical properties, leaf area, leaf angles, row spacing, planting density, etc.) to traditional and genetic engineering based crop breeders, for maximizing target functions such as light use efficiency, water use efficiency, drought tolerance, crop yield, crop resilience etc., with important applications in biofuel crops such as sorghum, and switchgrass. Our model can also be applied for designing optimal canopy architecture that increase surface albedo as a mitigation strategy for global warming and projections of crop yield for future food and water security under climate change.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Surface stability and H adsorption and diffusion near surfaces of W borides: a first-principles study

Abstract Understanding the behavior of tungsten boride (W x B y ) surfaces in a fusion reactor environment is an important topic since boronization is a common wall conditioning method used in fusion Tokamaks. We report the results of density functional theory calculations that investigate the surface stability of W x B y (tetragonal I4 1 /amd -WB, hexagonal P 6 3 / mmc -WB 2 and tetragonal I4/m- W 2 B) with low-index orientations, as well as hydrogen (H) energetics near W x B y surfaces. For single element terminated W x B y surfaces, B terminated surfaces are more energetically stable than W terminated as a result of significant reconstruction of B. The H surface adsorption energy and activation energy of H diffusion penetration below W x B y surfaces are mainly related to the outer termination. Specifically, the WB(001) surface terminated with two B layers, referred to as WB(001)-T BB , has higher H adsorption affinity and lower H diffusivity on this surface than other terminations, which is controlled by the significant charge transfer from B to H. However, B atoms on the WB 2 (0001)-T BB surface decrease both H adsorption and diffusivity on the surface, but enhance H diffusion below the surface in comparison to W terminated WB 2 (0001) surface. H would be trapped and diffuse within atomic surface gaps on the WB 2 ( 2 1 ˉ 1 ˉ 0 ) surface, while H below the surface layer would jump along the [0001] direction rather than diffuse into bulk. The surface diffusion activation energy of H on the W 2 B(001) surface slightly varies with terminations. Once H crosses the surface layer of W 2 B(001) with either termination, it prefers to diffuse into the bulk, or back towards the surface, rather than move parallel to the surface. Interestingly, WB 2 (0001) and WB 2 ( 2 1 ˉ 1 ˉ 0 ) surfaces will have relatively higher H retention than the other W x B y surfaces evaluated in this work.

Yang, L. (ORCID:0000000322166071)↗
DOE OSTI2021

Experimental Investigation of Elastodynamic Nonlinear Response of Dry Intact, Fractured and Saturated Rock

Nonlinear elastodynamic response of fractured rocks carries crucial information on fracture features that can be exploited to forecast flow properties, friction constitutive behavior and poromechanical response. Well-controlled laboratory experiments are designed to measure the nonlinear elastodynamic response of Westerly granite in three states: dry intact, dry fractured and saturated fractured. We study the effect of fracturing and saturation in modifying the elastodynamic response of the rock. Each sample is tested at a normal stress level of 15 MPa. We measure the elastodynamic response of an intact L-shaped sample of Westerly Granite subjected to normal stress oscillations of prescribed amplitudes (0.2–1 MPa) and frequencies (0.1, 1, 10 Hz). Ultrasonic waves transmitted across the sample are used to monitor the evolution of wave velocity before, during and after dynamic stressing. Furthermore, the nonlinearity of the elastodynamic response is measured in terms of: (1) the offset in normalized wave velocity; (2) the amplitude of wave velocity fluctuation during the oscillations; and (3) recovery rate of the wave velocity post-oscillation. We observe that the three nonlinearity parameters show a similar trend. Irrespective of the parameter, the nonlinearity measures higher for sample in dry-intact condition than that for dry-fractured and the saturated-fractured sample exhibits smaller nonlinearity than the dry-fractured sample. As expected, the saturated sample exhibits less nonlinearity than the dry intact and fractured samples due to the presence of interstitial fluid and the resulting increased interface stiffness. Conversely, the dry intact rock shows a higher nonlinearity than the dry fractured. We use numerical simulations to show that the presence of fracture significantly alters the strain distribution across the bulk of the sample and only the contacting asperities are highly strained, thus resulting in a decrease in the measured elastodynamic nonlinearity.

58 GEOSCIENCES↗
DOE OSTI2025

Epoxide Alcoholysis over M-BEA Zeolites: Effects of Alcohol Chain Length on Rates and Regioselectivities

The structures of nucleophilic reactants affect their coordination behavior among solvent molecules and kinetics of reactions with surface intermediates within the confines of fluidfilled pores of zeolites and other microporous materials. Consequently, rates and regioselectivities of diverse chemistries may depend sensitively on nucleophile identity in manners not observed for classic fluid phase reactions. Here, we examine the impact of varying the primary alcohol (ROH) chain length on the kinetics of 1,2-epoxybutane (C 4 H 8 O) ring-opening within Brønsted (Al-BEA) and Lewis acid (Zr-BEA) zeolites. Turnover rates increase by factors of ∼6 (Al-BEA) and 4-fold (Zr-BEA) between methanol and 1-hexanol, yet the reaction mechanisms remain comparable. Despite modest rate differences, apparent activation enthalpies calculated from rates and activities of solvated reactants decrease linearly by 12 (Al-BEA) to 33 kJ mol −1 (Zr- BEA) with increased proton affinity, which suggests bond formation energies for the nucleophile strongly influence rate increases. The molecular interpretation of these trends demonstrates, however, that the solvation of ring-opening transition states by zeolite pore structures and solvent molecules also governs rates. The impact of local solvating interactions appears most directly as changes in regioselectivities, which tend to enhance terminal alcohol formation with increasing ROH chain length. Regioselectivities largely do not vary with differences in fluid composition for a given ROH. The addition of H 2 O increases the number of hydrogen bonds among reactive species, and trends in regioselectivities imply that the decreased hydrogen bonding ability of longer chain ROH, and not the nucleophile strength or steric bulk, determines the regioselectivities of the resulting products. This work provides direct experimental evidence that nucleophilicity and hydrogen bonding influence reaction barriers and regioselectivities in zeolitecatalyzed epoxide ring-opening, offering pathways to better control reaction kinetics.

acidic zeolites↗
DOE OSTI2025

Modeling active nematics via the nematic locking principle

Active nematic systems consist of rod-like internally driven subunits that interact with one another to form large-scale coherent flows. They are important examples of far-from-equilibrium fluids, which exhibit a wealth of nonlinear behavior. This includes active turbulence, in which topological defects in the nematic order braid around one another in a chaotic fashion. One of the most studied examples of active nematics consists of a dense two-dimensional layer of microtubules, crosslinked by kinesin molecular motors that inject extensile deformations into the fluid. Though numerous theoretical studies have modeled microtubule-based active nematics, no consensus has emerged on how to fully and quantitatively capture the features of the experimental system. Here, to better understand the theoretical foundations for modeling this system, we propose a fundamental principle we call the nematic locking principle—individual microtubules cannot rotate without all neighboring microtubules also rotating. Physically, this is justified by the high density of the microtubules, their elongated nature, and their corresponding steric interactions. We assert that the nematic locking principle holds throughout the majority of the material, but breaks down in the neighborhood of topological defects and other regions of low density. We derive the most general nematic transport equation consistent with this principle and also derive the most general term that violates it, introducing fracturing into the material. We then examine the standard Beris–Edwards approach, commonly used to model this system, and show that it violates the nematic locking principle throughout the majority of the material due to fracturing. We then propose a modification to the Beris–Edwards model that enforces nematic locking nearly everywhere. This modification shuts off fracturing except in regions where the order parameter (a proxy for density) is reduced. In these regions fracturing is turned on. The resulting simulations in turn show strong nematic locking throughout the bulk of the material, with narrow bands of fracturing, consistent with experimental observation. One additional advantage of enforcing nematic locking is that nontrivial stationary state solutions, common in Beris–Edwards simulations but not seen in experiments, are eliminated.

Mitchell, Kevin A. [Univ. of California, Merced, C↗
DOE OSTI2020

Physical and Thermomechanical Properties of Yttrium Hydride from Large Scale Bulk Metal Hydriding Furnace

Given the superior thermal stability and highly attainable hydrogen density, yttrium hydride is an excellent high-temperature moderator material in advanced thermal neutron spectrum reactors that require small core volumes. Yttrium hydride has been selected as the moderator material for the Transformational Challenge Reactor, which was launched at Oak Ridge National Laboratory (ORNL) in 2019. However, fabrication of large-scale crack-free yttrium hydride is challenging and very limited efforts have been committed to the characterization of bulk yttrium hydride in response to the need to establish a complete database of the thermomechanical properties of YHx. In this report, the challenges associated with fabricating large-scale crack-free yttrium hydride are discussed herein. In response to those challenges, a hydriding system was designed and constructed at ORNL and was used to successfully fabricate crack-free yttrium hydride in complex geometries at large scales. This was accomplished by precisely controlling the hydrogen’s partial pressure and the retort temperature, which was informed by the well-established thermodynamic properties of the binary H-Y system. Hydrogen content in as-fabricated hydride was determined by the weight change method and vacuum hot extraction technique, complemented by the X-ray diffraction (XRD). In addition, significant efforts are being dedicated to establishing a complete database of the thermomechanical properties of as-fabricated yttrium hydride. In FY2020, we investigated the thermophysical properties of yttrium hydrides as a function of temperature (room temperature to 700°C) and hydrogen concentration (H/Y ratio ranges from 1.52 to 1.93). The results indicate that at the temperatures below 300 °C, the hydrogen content did not have a significant influence on the thermal expansion, while the specific heat capacity, the thermal diffusivity, and the calculated thermal conductivity were slightly higher for the higher H/Y ratio. Between 300°C and 700 °C, a reversible second-order endothermic transition in all measured thermal properties was observed. It was also found that the onset temperatures of the observed transition varied, with the composition having inverse dependence on the hydrogen content. An attempt was made to explain the behavior of the thermophysical properties at higher temperatures by considering the order– disorder transition as a result of hydrogen redistribution. In addition, nanoindentation was employed to determine the elastic modulus and hardness and to capture the crystal orientation dependence of these parameters. Vickers hardness was also reported. The final section of the report introduces ongoing neutron irradiation campaign of yttrium hydride.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Exploring the hardness and high-pressure behavior of osmium and ruthenium-doped rhenium diboride solid solutions

Rhenium diboride (ReB 2 ) exhibits high differential strain due to its puckered boron sheets that impede shear deformation. Here, we demonstrate the use of solid solution formation to enhance the Vickers hardness and differential strain of ReB 2 . ReB 2 -structured solid solutions (Re 0.98 Os 0.02 B 2 and Re 0.98 Ru 0.02 B 2 , noted as “ReOsB 2 ” and “ReRuB 2 ”) were synthesized via arc-melting from the pure elements. In-situ high-pressure radial x-ray diffraction was performed in the diamond anvil cell to study the incompressibility and lattice strain of ReOsB 2 and ReRuB 2 up to ~56 GPa. Both solid solutions exhibit higher incompressibility and differential strain than pure ReB 2 . However, while all lattice planes are strengthened by doping osmium (Os) into the ReB 2 structure, only the weakest ReB 2 lattice plane is enhanced with ruthenium (Ru). These results are in agreement with the Vickers hardness measurements of the two systems, where higher hardness was observed in ReOsB 2 . The combination of high-pressure studies with experimentally observed hardness data provides lattice specific information about the strengthening mechanisms behind the intrinsic hardness enhancement of the ReB 2 system.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modification of Ni-20Cr corrosion dealloying behavior in molten fluorides via cold work induced plastic deformation

The corrosion dealloying behavior of cold-worked (CW) Ni20Cr alloy (wt%) was studied in molten LiF-NaF-KF (or FLiNaK) salts at 600 °C, equal to a homologous temperature (TH) of 0.52. Alloys were cold-rolled to achieve reductions of thickness of 10%, 30%, and 50% introducing plastic deformation and a high density of dislocations. Potentiostatic holds (Eapplied) were applied in two different electrode potential regimes. At 1.75VK+/K$${\text{V}}_{{\text{K}}^{+}/{\text{K}}}$$, Cr dealloying to Cr(II) and Cr(III) is predominant, while at 1.90 VK+/K$${\text{V}}_{{\text{K}}^{+}/{\text{K}}}$$, both Ni and Cr are oxidized in molten FLiNaK at 600 °C. In these potential regimes, dealloyed Ni20Cr displayed bicontinuous porosity within the grain interior and at grain boundaries, driven by the high driving force for Cr dissolution and sustained by defect mediated outward solid state diffusion of Cr in parallel with surface diffusion of Ni. The bicontinuous porous structure developed was observed to undergo further coarsening and densification of the Ni-rich ligaments at a higher electrode potential. The main effect of CW observed is the introduction of plastic deformation and dislocation substructures that serve as short-circuit paths for Cr solid state diffusion to surfaces exposed to FLiNaK. This modified the evolution of the bicontinuous porous structure which increased with CW substantially. Kinetic analysis reveals that the Cr dealloying at +1.75 VK+/K$${\text{V}}_{{\text{K}}^{+}/{\text{K}}}$$ and 1.90 VK+/K$${\text{V}}_{{\text{K}}^{+}/{\text{K}}}$$ is initially charge transfer controlled, except for the 50% CW condition at 1.90 VK+/K$${\text{V}}_{{\text{K}}^{+}/{\text{K}}}$$, where the process becomes limited by slow Cr defect mediated bulk diffusion. The rate determining factors are explored and compared to experimental results.

Chan, Ho Lun↗
DOE OSTI2023

Metal forming and working of stabilized nanocrystalline Cu-Ta for electrical contacts

The commercialization of nanocrystalline metals and alloys is currently entering a renaissance period. Many of the processing and consolidation challenges that have haunted them are now more fully understood, opening the doors for stabilized nanocrystalline metals to be produced on a bulk scale. While challenges remain, the increased volume at which these materials are being supplied is for the first time allowing for investigations into more traditional methods of metal working, such as extruding, rolling, forming, and forging. Recently, the manufacturing science has been developed to allow nanocrystalline Cu-Ta alloys to progress to this point. This article therefore builds upon the last decade of evolutionary progression within the family of stabilized nanocrystalline Cu-Ta alloys by presenting some of the first findings related to scaled powder synthesis, production of billets, thin sheets, and foils. Here, the mechanical performance and physical properties relevant to forming electrical contacts and pins including, tensile, J-integral fracture toughness, Charpy Impact, bi-axial tension, and conductivity are reported. This introductory investigation into forming such a novel material, provides evidence to these alloys potential at bridging the gap between being a scientific curiosity to that of a real engineering material.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Model-based Hierarchical Reinforcement Learning for Improved Physical Security Design: A Prototype

Prior work in FY24 developed an adversarial AI agent aid in path analysis of physical protection systems. This agent, trained using a model-based reinforcement learning algorithm, was able to successfully learn the most vulnerable path in facilities. It was able to extend the current state of practice for physical protection design by exhibiting dynamic behavior based on current environmental conditions. Whereas PathTrace largely performs a static, graph-based analysis, the AI agent was able to make decisions based on relative position in the facility, current conditions (was the adversarial agnet discovered?), and proximity to secondary targets. The agent demonstrated some novel capabilities, but had limitations that need to be resolved before it can be used for production purposes. For example, the adversarial agent generalizes poorly and takes a relatively long time to train. Nonetheless, there is still considerable promise for developing the adversarial agent further in order to explore even richer, more dynamic behaviors (e.g., adversary motivations, environmental debris, and more). This work considers a complementary idea; development of a planning agent. The planning agent is envisioned as an auto-complete-like tool that can help accelerate security system design by human experts. The agent would respect existing barriers and sensors placed by a human expert while offering cost-effective suggestions (i.e., implicitly balancing effectiveness with cost) to improve the design. The goal is for this agent to be part of an expert’s toolbox, not to totally upend the current state-of-practice, or to displace human experts. The ultimate goal would be concurrent training of both the adversarial and planning agent together, to learn entirely through self-play. This would represent an entirely new way of performing system deign. We selected a hierarchical, model-based reinforcement learning algorithm to serve as the planning agent. This is an extension of concepts used in the prior FY24 adversarial agent work. There, we had a single agent acting an environment. Here, we have two different sub-agents (policies), working together, to form a complete agent. There is a manager policy, which can select abstract goals on slower time scales, and a worker, which performs primitive actions to reach goals selected by the manager. It is worth noting that this class of algorithm is challenging to work with. From our understanding, our work is one of the first successful uses of model-based reinforcement learning (MBRL) in nuclear energy1 , and likely the first hierarchical model-based reinforcement learning application in nuclear energy. Further, this work is one of the first known attempts to apply AI to perform a design tasks in nuclear energy. Consequently, there were significant implementation challenges and the bulk of the work was focused on successful implementation and algorithm design. The results presented here are very low technology readiness level as a consequence of the lack of related literature, but still represent a significant step forward in the pursuit of applied AI for design.

42 ENGINEERING↗
DOE OSTI2026

Determining the Solubility Behavior of Kogarkoite in Simulated Nuclear Waste

Kogarkoite (Na 3 FSO 4 ) is a sparingly soluble fluoride–sulfate double salt that has been identified in high level nuclear waste sludge at the Hanford Site and, more recently, in sludge batch compilation samples at the Savannah River Site (SRS). Due to its complex dissolution behavior, which exhibits an inverse dependence on sodium ion activity, the presence of this mineral poses significant challenges to waste retrieval and processing. Incomplete dissolution during sludge washing can lead to the retention of fluoride and sulfate in the high-level waste feed, potentially causing the formation of corrosive, immiscible molten salt layers, known as "glass gall,” in vitrification melters. Current efforts to optimize flowsheet parameters and wash-water volumes are hindered by the absence of a commercially available, certified reference material, which prevents the accurate calibration of analytical methods and the verification of dissolution kinetics. To address this critical gap, this research focuses on the laboratory synthesis of pure Kogarkoite to serve as a standard for comprehensive solubility and washing performance testing. A coupled synthesis and simulant campaign was executed using an evaporative crystallization protocol designed to replicate the dynamic concentration effects observed in tank farm operations. Thirteen simulant matrices were prepared by dissolving systematically varied ratios of sodium fluoride (NaF) and sodium sulfate (Na 2 SO 4 ) in deionized water under three distinct caustic regimes: 0.0 g (control), 4.0 g (~1 M), and 12.0 g (~3 M) sodium hydroxide (NaOH). While thermodynamic equilibrium models suggest that high-caustic environments should favor the stability of the double salt7, results from this evaporative study at 25 0 C revealed a distinct kinetic divergence. Simulants with high hydroxide loading predominantly yielded large, blocky crystals of sodium sulfate decahydrate (Na 2 SO 4 .10H 2 O). Successful synthesis of pure Kogarkoite was achieved exclusively in specific NaOH-free compositional windows, where the precipitate manifested as fine, opaque granular aggregates. Ion chromatography (IC) analysis confirmed phase purity through the simultaneous stoichiometric depletion of both fluoride and sulfate from the supernatant. This successful synthesis establishes a reproducible route to generate bulk Kogarkoite, enabling the subsequent phase of quantitative dissolution testing using inhibited water to optimize sludge-batch assembly.

Sarker, Md Sharif [Florida International Univ. (FI↗
DOE OSTI2023

Demonstration of Scaled-Production of Rare Earth Oxides and Critical Materials from U. S. Coal-Based Sources (Final Report)

The project objective was to demonstrate scaled production of high purity rare earth oxides (REO), nominally exceeding 90% grade, from coal refuse sources using innovative technologies that reduce cost and improve environmental outcomes relative to traditional rare earth processing technologies. The project utilized a critical material pilot plant constructed and tested as part of a previous U.S. Department of Energy project. Target performance criteria was a 50% reduction in production cost based on previous optimum values, 150% increase in recovery and greater than 2% concentrates of rare earth oxides, cobalt and manganese. Concentrate production goals were to produce a rare earth mix oxide product at a rate of 200 grams per day having a minimum purity of 50% as well as products of cobalt and manganese having a minimum purity of 2%. A previous pilot plant investigation identified acid cost as the major contributor to an operating cost that made the recovery of rare earth and other critical metals from bituminous coal sources economically challenging. As such, acid cost reduction was a major target using bio-oxidation reactors to produce sulfuric acid from naturally occurring coal pyrite. Based on laboratory data, a bio-oxidation circuit was designed for the pilot plant to produce 7.5 l/min of acid using two 11-m3 (3000 gallon) reactors equipped with 40 hp aerators for air dispersion. The pilot-scale tests revealed that acid concentration equivalent to as high as 1.0 M sulfuric acid could be continuously produced. However, the bio-acid contained exceptionally high iron concentrations, which complicated downstream processing of the pregnant leach solution (PLS). A TEA of the bio-oxidation circuit showed that production cost was approximately $0.13 per kg acid equivalent if produced using a four-day retention time in the reactors. This value represents a 48% reduction from that of purchased bulk sulfuric acid. Calcination (or roasting) studies were conducted on coarse refuse from West Kentucky No. 13 and Fire Clay coal seam sources. The test results revealed the potential to increase recovery by nearly 100% using temperatures between 500°C to 700°C with light REE recovery value being the most improved. Acid baking of the calcined products using sulfuric acid at 250°C increased heavy rare earth recovery from around 40% to 80% while decreasing the acid requirements by over 50%. The existing pilot plant was upgraded for the pilot scale demonstration of REE and CM recovery. The primary feedstocks were West Kentucky No.13 and heap leach pregnant leach solution (PLS) while a secondary feedstock was a lignite waste material from a construction sand operation. The pilot scale operation started with PLS generation through leaching followed by iron and aluminum removal, nominally at 3.3 and 4.5 pH, respectively. Leaching lixiviants used for the test were industrial grade sulfuric acid or bio-acid generated at the pilot scale facility. For most of the tests, the solid feed rate was 200 lb/hr whereas lixiviant was added at 2 gpm to provide an optimal residence time of 45 minutes. Following the contaminant removal step, several different process schematics were tested with the goal of maximizing REE recovery and purity. In the first test, direct processing of aluminum precipitation raffinate for REE recovery using oxalic acid at pH 1.5 was investigated. Overall REE recovery was approximately 45%. Unfortunately, elevated calcium content in the PLS significantly impacted the RE-Oxide product grade. Similarly, high calcium content decreased both the product purity and grades of CM products. As such, a new flowsheet was tested with the same initial process schematic but different precipitation stages for REEs and CMs at pH 6.0 and 9.0, respectively. It was noted that the overlapping precipitation behavior of Co, Ni and Zn with REEs limited the applicability of this process schematic. While this change increased the RE-Oxide grade from 36% in the first test to 87%, the loss of critical metals to the REE cake and bypass of the REEs to the CM cake significantly impacted the recovery of both REEs and CMs. Therefore, the modified process flowsheet combined oxalic acid precipitation stage raffinate and redissolved CM cake filtrate to maximize both the recovery and purity of the products. Consequently, a RE-Oxide product with 85% purity and CM cakes with over 19% Co, 38% Ni, 14% Zn and 9% Mn content were generated with significantly higher recoveries. While the modified process flowsheet improved recoveries and grades, elemental losses observed in separate precipitation and redissolution losses inspired the adaptation of a single precipitation stage at pH 9.0 for both REEs and CMs. This change was anticipated to maximize the REE recovery while minimizing the costs associated with separated redissolution and processing stages. As expected, REE recovery in this new circuit arrangement was over 56% with a product grade of over 87% RE-Oxide content. Similarly, Co, Ni, Mn, and Zn recoveries of 54%, 40%, 67%, and 66%, respectively, were achieved. Unfortunately, the elevated calcium content present in the solution due to its precipitation at pH 9.0 caused a decrease in the CM cake quality. Therefore, the final process flowsheet involved the addition of calcium oxalate precipitation following the oxalic acid precipitation stage, which effectively eliminated calcium contamination of the CM products. Finally, the pilot scale experiments conducted using bio-acid achieved comparable REE leaching recoveries to the conventional sulfuric acid leaching. Elevated iron concentration in the solution caused the co-precipitation of REEs with the iron and aluminum cake, resulting in the REE losses. Furthermore, elevated iron content bypassing the iron and aluminum precipitation stages contaminated the metal sulfide and manganese cake, respectively. A techno-economic analysis was performed based on a commercial facility capable of treating 500 tph of coal-based material. The production cost for West Kentucky No. 13 coarse refuse material ranged from approximately $500-$700/kg of total rare earth oxide whereas the lignite source had significantly lower production cost of $100-$300/kg. The significant difference in cost was due to the easier leaching characteristics of the lignite material and the higher feed concentrations. All process scenarios resulted in a negative net present value (NPV). For the lignite feedstock, laboratory REE leach recovery values were about 30% higher than the pilot plant data. Using the lab leach results, a positive net present value was achieved and the production cost decreased from $100-$300 $/kg to less than $150/kg of total REO.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Magnetic and dielectric property control in the multivalent nanoscale perovskite Eu 0.5 Ba 0.5 TiO 3

We report nanoscale Eu 0.5 Ba 0.5 TiO 3 , a multiferroic in the bulk and candidate in the search to quantify the electric dipole moment of the electron. Eu 0.5 Ba 0.5 TiO 3 , in the form of nanoparticles and other nanostructures is interesting for nanocomposite integration, biomedical imaging and fundamental research, based upon the prospect of polarizability, f-orbital magnetism and tunable optical/radio luminescence. We developed a [non-hydrolytic]sol–[H 2 O-activated]gel route, derived from in-house metallic Ba (s) /Eu (s) alkoxide precursors and Ti{(OCH(CH 3 ) 2 } 4 . Two distinct nanoscale compounds of Ba:Ti:Eu with the parent perovskite crystal structure were produced, with variable dielectric, magnetic and optical properties, based on altering the oxidizing/reducing conditions. Eu 0.5 Ba 0.5 TiO 3 prepared under air/O 2 atmospheres produced a spherical core–shell nanostructure (30–35 nm), with perovskite Eu 0.5 Ba 0.5 TiO 3 nanocrystal core-insulating oxide shell layer (~3 nm), presumed a pre-pyrochlore layer abundant with Eu 3+ . Fluorescence spectroscopy shows a high intensity 5 D 0 → 7 F 2 transition at 622 nm and strong red fluorescence. The core/shell structure demonstrated excellent capacitive properties: assembly into dielectric thin films gave low conductivity (2133 GΩ mm -1 ) and an extremely stable, low loss permittivity of ε eff ~25 over a wide frequency range (tan δ < 0.01, 100 kHz–2 MHz). Eu 0.5 Ba 0.5 TiO 3 prepared under H 2 /argon produced more irregular shaped nanocrystals (20–25) nm, with a thin film permittivity around 4 times greater ( ε eff 101, tan δ < 0.05, 10 kHz–2 MHz, σ ~59.54 kΩ mm -1 ). Field-cooled magnetization values of 0.025 emu g -1 for EBTO-Air and 0.84 emu g -1 for EBTO-Argon were observed. X-ray photoelectron spectroscopy analysis reveals a complex interplay of Eu II/III /Ti III/IV configurations which contribute to the observed ferroic and fluorescence behavior.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Assessment of Sodium Thermal Stratification Models Utilizing the TSTF Benchmark

As a result of certain transient scenarios, a thermally stratified layer of liquid sodium can develop in the bulk coolant volumes of a sodium-cooled fast reactor (SFR). In addition to the effects a stratification layer has on the temperature of the heat transport system, a stratification layer can also influence the transition to and establishment of natural circulation flow, which plays an important role in passive cooling and the inherent safety of a pool-type SFR. Therefore, the ability to accurately capture thermal stratification phenomena is important when demonstrating the safety basis of a pool-type SFR during transient sequences. The present work assesses various computational models with different fidelities in their ability to predict thermal stratification in the upper plenum of an SFR. Each computational model will be assessed using the data generated at the Thermal Stratification Test Facility (TSTF) located at the University of Wisconsin-Madison. Using measured flow rate and inlet temperature data, the measured temperature distributions of the tests are compared to the predictions of the lumped volume-based models in SAS4A/SASSYS-1, a 1D-based model in SAM, and a 3-D computational fluid dynamics (CFD) model using STAR-CCM+. The relative performance of the various computational methods is assessed with respect to key metrics such as bulk coolant temperature distribution and plenum exit temperature. A total of eight tests are analyzed, covering different combinations of flow rates (3 and 10 GPM) and upper internal structure (UIS) configurations (none, solid, porous, and open) The perfect mixing model of SAS4A/SASSYS-1 provides the highest accuracy when the flow rate is high and there is no UIS in the test vessel, as high flow rate injection promotes thermal mixing of the sodium in the test vessel. For most of the analyzed tests, the stratified volume model of SAS4A/SASSYS-1 is able to predict the delay in the outlet temperature drop and temperature distribution in the test vessel by a small number of layers to represent thermal stratification. However, the stratified volume model can only simulate a maximum of three temperature layers within a volume and when a layer approaches the elevation of the outlet, the predicted outlet temperature can demonstrate rapid, non-physical changes. The 1-D axial mixing model of SAM provides results that agree reasonably well with the measured data in the prediction of the temporal evolution of the outlet temperature with the exception of the case with a high flow rate and no UIS. The SAM 1-D model has a similar level of accuracy to CFD results when it comes to predicting the outlet temperature. CFD shows overall good agreement in predicting the temperature distribution in the test vessel and outlet temperature. As CFD can model the test vessel geometry in detail, it performs well in the cases of complex geometries such as tests that included a UIS and internal flow through the UIS resulting in active mixing of the coolant in the test vessel. Each of the models discussed in the present work has the potential to be useful during the various stages of reactor design, analysis, and licensing. The lumped-volume approach can be applied for fast turnaround safety calculations to obtain overall reactor behavior during transients. The 1-D models provide improved accuracy when stratification is expected for a relatively low increase in the computational cost. The CFD model can be utilized for confirmatory analysis of the 1-D model, when experimental measurements are not available.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗