SEARCH · Engineering Papers

Results for “time dependent 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.

SourceResearch connectionFrom yearThrough year

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

At least 757 records · Page 42

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Enhanced Diffusive Transport in Fluctuating Porous Media

Mass transport within porous structures is a ubiquitous process in biological, geological, and technological systems. Despite the importance of these phenomena, there is no comprehensive theory that describes the complex and diverse transport behavior within porous environments. While the porous matrix itself is generally considered a static and passive participant, many porous environments are in fact dynamic, with fluctuating walls, pores that open and close, and dynamically changing cross-links. While diffusion has been measured in fluctuating structures, notably in model biological systems, it is rarely possible to isolate the effect of fluctuations because of the absence of control experiments involving an identical static counterpart, and it is generally impossible to observe the dynamics of the structure. In this paper, we present a direct comparison of the diffusion of nanoparticles of various sizes within a trackable, fluctuating porous matrix and a geometrically equivalent static matrix, in conditions spanning a range of regimes from obstructed to highly confined. The experimental system comprised a close-packed layer of colloidal spheres that were either immobilized to a planar surface or allowed to fluctuate locally, within the space defined by their nearest neighbors. Interestingly, the effective long-time diffusion coefficient was approximately 35–65% greater in the fluctuating porous matrix than in the static one (depending on the size of the nanoparticle probes), regardless of the geometric regime. This was explained by considering the enhancing effects of matrix fluctuations on the short-time diffusion coefficient and cooperative “gate-opening” motions of matrix particles and nanoparticle probes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Fuel reid vapor pressure level and ethanol content on stochastic preignition, effects at steady and unsteady engine operation

The present work investigates relations between fuel Reid vapor pressure (RVP) and biofuel (ethanol) content on stochastic preignition (SPI) at both sustained steady-state engine operation and following load transients. This work stems from in-field observations that automotive original equipment manufacturers have observed consistent seasonal increases in United States customer drivability complaints and warranty claims during September and October where SPI is suspected to be responsible. The seasonal timing of these events coincides with the United States seasonal fuel property changeover initiating on September 15 each year, where fuel RVP increases. To explore potential linkage between fuel RVP and SPI the present study employs engine SPI experiments coupled with laboratory spray measurements of fuels with RVPs of 8, 12, and 16 psi in both E10 (10% ethanol) and E25 (25% ethanol) fuels. Engine results are partitioned into fuel RVP and ethanol content effects on SPI in steady-state, sustained high-load engine operation and unsteady-state low- to high-load transitions, where off-engine spray vessel patternation and tip penetration results help to elucidate the observed fuel effects on SPI. A boosted direct-injected, spark-ignition engine was fueled with three market relevant E10 and E25 fuels with RVPs of 8, 12, and 16 to characterize the interplay between winter fuels and abnormal combustion behavior. The steady-state work shows that for high-load, steady-state engine operation, SPI is directly linked to fuel retention, which was found to be dependent on fuel distillation. The unsteady-state engine operation work shows that following low-to high-load transitions, SPI can occur from a memory of fuel property effects at low-load operation. Specifically, the fuel RVP effect on fuel spray collapse at low loads was found to correlate with SPI with a more than 95% confidence interval following low- to high-engine-load transitions. Results suggest that fuel-wall impingement at low-load operation could carry over into high-load transitions and generate SPI events following low- to high-load transitions.

09 BIOMASS FUELS↗

DOE OSTI2025

Monitoring Fracture Hydromechanical Evolution in the Lab and Field Using Unsupervised Metric Learning

Fractures evolve in time through thermal‐hydraulic‐mechanical‐chemical (THMC) processes that alter their long‐range hydraulic transport properties and modify subsurface behavior and activities. The location of subsurface fractures makes it necessary to use remote sensing techniques such as passive or active seismic monitoring for fracture characterization. In this paper, we develop a machine learning approach to monitor the evolution of fracture properties using passive seismic sources in a laboratory setting and using active seismic monitoring from the Sanford Underground Research Facility in Lead, South Dakota, at a depth of 1.25 km in amphibolite rock during stimulation of natural fractures as well as during induced fracturing. The unsupervised metric learning technique applies tandem neural networks (twin (Siamese) or triplet) with contrastive loss and adaptive margins to track slowly varying systems for which class or similarity labels are not available. The approach adopts locality‐sensitive hashing to divide time‐ordered contiguous data into an arbitrary number of pseudo‐classes. Contrastive‐loss training with many hash bins generates an evolving latent‐space trajectory. This approach enables unsupervised metric learning for seismic data stacks under the condition of contiguous state sampling and slowly varying fracture properties. The displacement discontinuity theory provides a mechanistic foundation for the fracture‐dependent trajectories that are related to relaxation of fractures with time‐dependent specific stiffness responding to changes in stress or fluid saturation.

02 PETROLEUM↗

DOE OSTI2022

Characterization of single-shot attosecond pulses with angular streaking photoelectron spectra

Most of the traditional attosecond pulse retrieval algorithms are based on a so-called attosecond streak camera technique, in which the momentum of the electron is shifted by an amount depending on the relative time delay between the attosecond pulse and the streaking infrared pulse. Thus, temporal information of the attosecond pulse is encoded in the amount of momentum shift in the streaked photoelectron momentum spectrogram S(p,τ), where p is the momentum of the electron along the polarization direction and τ is the time delay. An iterative algorithm is then employed to reconstruct the attosecond pulse from the streaking spectrogram. This method, however, cannot be applied to attosecond pulses generated from free-electron x-ray lasers where each single shot is different and stochastic in time. However, using a circularly polarized infrared laser as the streaking field, a two (or three)-dimensional angular streaking electron spectrum can be used to retrieve attosecond pulses for each shot, as well as the time delay with respect to the circularly polarized IR field. Here we show that a retrieval algorithm previously developed for the traditional streaking spectrogram can be modified to efficiently characterize single-shot attosecond pulses. The methods have been applied to retrieve 188 single shots from recent experiments. We analyze the statistical behavior of these 188 pulses in terms of pulse duration, bandwidth, pulse peak energy, and time delay with respect to the IR field. Furthermore, the retrieval algorithm is efficient and can be easily used to characterize a large number of shots in future experiments for attosecond pulses at free-electron x-ray laser facilities.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

Numerically Testing Conceptual Models of the Utah FORGE Reservoir Using July 2023 Circulation Test Data

Over the past several years, many new data sets have become available regarding the characterization of the Utah FORGE reservoir. These include, but are not limited to, the stimulation of Well 16A, the drilling and completion of Well 16B, and interwell circulation confirmatory testing. As part of the characterization efforts, conceptual models of the reservoir are re-examined as new data become available. As part of the planning for FORGE activities, numerical models are often used to predict the reservoir response to the planned testing. Stochastic methods are often employed to bound uncertainty and allow for evaluation of comprehensive ranges of key reservoir parameters. For the most recent interwell circulation confirmatory testing (July 2023), a priori numerical model predictions did bound the observed behavior (Xinj et al., 2023), but key deviations from expected behavior prompted the FORGE team to reevaluate our conceptual model of the reservoir. In early October 2023, key members of the development, testing, and monitoring teams met for 2 days to review newly collected data and discuss ‘interesting’ or ‘key’ observations. From these discussions, 15 Key Observations were documented, with several significant ones being that the discrete fracture network developed from the 16A stimulation data may not be appropriate and that the early time pressure data obtained during the summer 2023 reservoir testing were best described using radial solutions. In July 2023, two campaigns of interwell confirmatory testing were conducted, the first set of tests occurred on July 4-5 and the second set on July 18-19. The second set of circulation tests conducted at the Utah FORGE site between the injection well 16A(78)-32 and production well 16B(78)-32 on July 18 and 19, 2023 are used to calibrate material properties in a thermal-hydraulic-mechanical (THM) simulation of the discrete fracture network connecting the wells. The spatially and temporally varying reservoir properties are calibrated to match the time dependent pressure and production profiles from the circulation tests. In future work, this calibrated model will be coupled to the native state THM model of the FORGE reservoir to predict surface deformation and strains resulting from pumping schedules.

58 GEOSCIENCES↗

DOE OSTI2024

Numerically Testing Conceptual Models of the Utah FORGE Reservoir Using July 2024 Circulation Test Data

15 GEOTHERMAL ENERGY↗

DOE OSTI2024

Nature of polariton transport in a Fabry-Perot cavity

Fabry-Perot microcavities can strongly enhance interactions between light and molecules, leading to the formation of hybrid light-matter states known as polaritons. Polaritons possess much smaller effective masses and much larger group velocities when the molecules are resonant with cavity modes that have finite (nonzero) in-plane wave vectors, giving rise to the possibilities of long-range and ultrafast ballistic transport. In this paper, we present the results of numerical simulations of the ultrafast ballistic transport phenomenon in real space and time during and after initialization with a short, spatially localized pulse. Here, we address this problem with two approaches: the standard transfer-matrix method (TMM) for planar structures and a second simulation based on the numerical solution of the Maxwell-Bloch equations that can be used for general configurations (complex metasurfaces, for instance) and excitation modes. The agreement between the TMM and the full numerical calculation with the Maxwell-Bloch equations when applied to multilayer planar configurations provides proof of the validity of both approaches for the present analyses. Overall, we find that the transport of the molecular excitons inside the cavity synchronize with the evolution of the enhanced electromagnetic field inside the cavity. Moreover, the synchronized transport rate is in good agreement with the group velocities predicted from a calculated dispersion relation across a wide range of frequencies. Finally, we relate the group velocity to the Hopfield coefficient obtained from quantum modeling and suggest that the dependence of light-matter coupling on the in-plane wave vector can be an important but overlooked factor for understanding the transport behavior of polaritons. These simulations provide an intuitive tool for understanding the collective motion of light and excitons and helps us to better understand how experimental observations of polaritons should be interpreted.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

High Performance Equilibrium Solvers for Integrated Magnetic Fusion Simulations

With the funding provided by this award, we developed numerical codes for the study of magnetically confined plasmas for fusion applications. Accordingly, our work can be divided into two separate categories: 1) the design and analysis of novel numerical methods providing high accuracy and high efficiency; 2) the study of the equilibrium and stability of magnetically confined plasmas with some of these numerical codes, as well as the study of the nature of the turbulent behavior which may arise in the presence of instabilities. We first developed new numerical schemes based on integral equation methods for the computation of steady-state magnetic configurations in fusion experiments, providing high accuracy for the magnetic field and its derivatives, which are required for stability and turbulence calculations. We employed different integral formulations depending on the application of interest: axisymmetric or non-axisymmetric equilibria, force-free or magnetohydrodynamic equilibria, fixed-boundary equilibria or free-boundary equilibria. While efficient, these methods do not yet apply to plasma boundaries which are not smooth, a situation which is fairly common in magnetic confinement experiments. To address this temporary weakness, we also constructed a new steady-state solver based on the Hybridizable Discontinuous Galerkin (HDG) method, which provides full geometric flexibility. In addition to these numerical tools focused on steady-states, we also contributed to the improvement of the speed and accuracy of codes simulating the plasma dynamics of fusion plasmas, by developing a novel velocity space representation for the efficient solution of kinetic equations, which most accurately describe the time evolution of hot plasmas in fusion experiments. Using the tools discussed above, we studied several questions pertaining to the equilibrium and stability of magnetically confined plasmas. In particular, we derived a new simple model for axisymmetric devices called tokamaks, to predict how elongated a fusion plasma can be before it becomes unstable and collapses. We also looked at the effect of the shape of the outer plasma surface on key properties of the steady-state magnetic configurations, and how these properties impact turbulence in fusion plasmas, and the corresponding transport of momentum. Likewise, we studied the role of large localized flows on the steady-state magnetic configurations, and how they may influence plasma stability and turbulence. Non-axisymmetric steady-state magnetic configurations are inherently more complex than axisymmetric steady-state configurations, and the subject of ongoing controversies regarding the regularity of the equations determining such steady-states, and their solutions. Implementing an existing NYU code in a new geometry, we studied the nature of the singularity of the solutions observed in the code, and methods to eliminate them. Our main conclusion is that by appropriately tailoring the plasma boundary, it is possible to eliminate the singularities otherwise appearing in our simulations, and to obtain steady-states which appear to be smooth. To gain further insights on incompletely understood turbulence phenomena, we proposed a new reduced model capturing most of these phenomena, which is simple enough to not require expensive numerical simulations on massive supercomputers to investigate them. We demonstrated the strong similarity between our simulations and published results obtained from computationally expensive simulations, and plan to rely on our reduced model to identify the key mechanisms determining the evolution and strength turbulent driven transport in fusion plasmas. Finally, we proposed a new framework for tokamak reactor design studies, enabling us to consider the relative merits of steady-state versus pulsed fusion reactors. We found that pulsed fusion reactors may benefit most from recent advances in magnet technology, and the availability of very high field magnets. As such, they may become more desirable than steady-state tokamak reactors for cost efficient electricity generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

The Status and Ambitions of the US Heavy Element Program

The aim of this whitepaper is to highlight the current capabilities and priorities of the US Heavy Element community and to provide the framework for a coordinated advancement of nuclear science from these studies. This is an organized effort to reflect on what has been achieved in the field given the recommendations and initiatives of the 2015 NSAC Long Range Plan, and on what can be realized in the next decade given current and possibly expanded investments. Current investments have positioned the US community to be among the world leaders in studies of the nuclear and chemical properties of the heaviest elements. These include studies of reaction mechanisms, moving us ever closer to the “island of stability”, in spectroscopy, allowing us to better understand nuclear structure at these extreme proton numbers, in chemical behavior, looking to determine how these elements should be placed on the Periodic Table, in performing the first measurements where isotopes are directly identified by their mass numbers, and in laying the foundation for a potential US-led new element discovery experiment. At present, the highest priority of the US Heavy Element community is to capitalize on the current investments by supporting the operations of US facilities at optimal values. These facilities include the Argonne Tandem Linac Accelerator System at Argonne National Laboratory and other Department of Energy facilities such as the 88-Inch Cyclotron at Lawrence Berkeley National Laboratory, which has a dedicated superheavy element program, as well as university laboratories, including Texas A&M University. The High Flux Isotope Reactor at Oak Ridge National Laboratory is crucial to providing the radioactive isotopes required for heavy element science targets. This facility should be supported to provide the actinide materials that are essential for US-based science. Production of stable, rare isotopes for beam material, including 48 Ca, 50 Ti, 54 Cr and 58 Fe, at the Stable Isotope Production and Research Center is critical to continued research in heavy element science and should be a priority. The continued development of targets for heavy element science and retaining US-based expertise is critical for the heavy element community. This is an area that is currently under pressure. For example, the target laboratory at Argonne National Laboratory serves a broad community and is currently under threat due to loss of critical personnel. The skills needed to make targets for nuclear science and develop new targetry methods need to be supported long term at Argonne, Oak Ridge and Lawrence Livermore National Laboratories as well as maintaining the programs at Oregon State University and San José State University as vital pipelines for training students. Advances in theory are the foundation to understand how nuclei behave and to predict those behaviors in new circumstances. Progress in these studies will necessitate continued and new investment and access to high-performance computing. The future health of the heavy element field is dependent on the continuous support of talented early-career professionals at all levels. It is critical that opportunities continue to be created for the next generation to become established in heavy element research so that we can ensure the field is attracting and retaining the best minds for continued success. It is also clear that to ensure diversity of ideas, perspectives and techniques, we need to recruit diverse personnel that are trained at the best facilities. The heavy element community is in support of continued investment to programs with initiatives in diversity, equity, and inclusion. Looking to the next decade of research, support needs to maintain and grow US leadership in heavy element science. Specifically, new investments in state-of-the-art instrumentation will be essential to scientific development of the field and in expanding scientific knowledge. Advances in the next generation of electron cyclotron resonance ion sources, multi-reflection time-of-flight devices, laser spectroscopy, trapping methods and next generation alpha and gamma spectroscopy systems should be prioritized.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2024

A prediction model of failure threshold for shear deformation in a Zr-based bulk metallic glass

The failure of bulk metallic glasses (BMGs) during plastic deformation at room temperature is abrupt and instantaneous, while the analysis of precursor information based on avalanche events helps predict catastrophic failure. An acoustic emission (AE) signal can provide accurate precursor information for material failure, due to its sensitive and high fast calculation ability. In the current study, AE monitoring tests are carried out during uniaxial compression tests of BMGs at different strain rates. The AE experimental failure threshold, E max , is proposed on the basis of AE cumulative energy, which reflects the intensity of damage evolution at different loading conditions. Compared with the critical shear band velocity (CSBV) associated with stick-slip dynamics of serrated flow, E max is a more sensitive failure parameter since it is connected with the local microscopic changes that occur during the material response process. Here, the E max is obtained prior to reaching the CSBV since the calculation of these two avalanches analysis focuses on the different stages of shear band growth. In particular, AE events are related to the “dry” friction process in the first stage, however, the CSBV is responsible for the “viscous” glide in the second stage. Therefore, Emax is not affected by the complex interactions between the shear bands during the stick-slip process. The maximum avalanche of serrated flow, S max , is proposed as the experimental failure threshold, which depends on the applied strain rate as S max ~ $\dot{ε}$ –λ . According to the relationship of E max and S max , the theoretical failure threshold, E max , follows a criterion E max = 2545$\dot{ε}$ –λ - 4468, where λ is equivalent to 0.15 for this work. Finally, combining the different calculations and AE measurements, this model gives new insights to predict the deformation failure behavior of Zr-based BMGs.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Suppression of antiferromagnetic order and strong ferromagnetic spin fluctuations in Ca(Co 1-x Ni x ) 2-y As 2 single crystals

CaCo 2–y As 2 is a unique itinerant system having strong magnetic frustration. Here, we report the effect of electron doping on the physical properties resulting from Ni substitutions for Co. The single crystals of Ca(Co 1–x Ni x ) 2–y As 2 were characterized by single-crystal x-ray diffraction, energy-dispersive x-ray spectroscopy, magnetization M versus temperature T, magnetic field H, time t, and heat capacity C p (H, T) measurements. The A-type antiferromagnetic (AFM) transition temperature T N = 52 K for x = 0 decreases to 22 K with only 3% Ni substitution and is completely suppressed for x > 0.16. For 0.11 ≤ x ≤ 0.52 strong ferromagnetic (FM) fluctuations develop as revealed by magnetic susceptibility χ(T) = M(T) / H measurements. For x = 0.11 and 0.16 competing AFM and FM interactions result in a reentrant spin-glass behavior below T N , as evidenced by the observations of thermomagnetic hysteresis and magnetic relaxation. Enhanced FM fluctuations are also found for the x = 0.21 and 0.31 crystals, where χ c increases significantly at low T. A large χ anisotropy in these compositions where χ c is up to a factor of two larger than χ ab suggests that the FM spin fluctuations are quasi-1D in nature. Weak ferromagnetic contributions to the magnetization are found at T = 2 K for x = 0.11 –0.31. Heat-capacity Cp(T) measurements reveal the presence of FM quantum spin fluctuations for 0.11 ≤ x ≤ 0.52, where a logarithmic T dependence of Cp(T) / T is observed at low T. Here, the suppression of AFM order by the development of strong FM fluctuations in Ca(Co 1–x Ni x ) 2–y As 2 crystals suggests the presence of a FM quantum-critical point at x ≈ 0.20. Our density-functional theory (DFT) calculations confirm that FM fluctuations are enhanced by Ni substitutions for Co in CaCo 2–y As 2 . The Sommerfeld electronic heat-capacity coefficient is enhanced for x = 0, 0.21, and 0.42 by about a factor of two compared to DFT calculations of the density of states (DOS) at the Fermi energy, suggesting an enhancement of the DOS from electron-phonon and/or electron-electron interactions. The crystals with x > 0.52 do not exhibit FM spin fluctuations or magnetic order at T ≥ 1.8 K, which was found from the DFT calculations to arise from a Stoner transition. Superconductivity is not observed above 1.8 K for any of the compositions. Neutron-diffraction studies of crystals with x = 0.11 and 0.16 in the crossover regime (0.1 ≲ x ≲ 0.2) show no evidence of A-type ordering as observed in the parent compound with x = 0. Furthermore, no other common magnetic structures, such as ferromagnetic (FM), helical stacking of in-plane FM layers, or in-plane AFM structure, are found with an ordered moment greater than the uncertainty of 0.05 μ B per transition-metal atom.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Quartz Dissolution Effects on Flow Channelization and Transport Behavior in Three‐Dimensional Fracture Networks

We perform a set of reactive transport simulations in three-dimensional fracture networks to characterize the impact of geochemical reactions on flow channelization. Flow channelization, a frequently observed phenomenon in porous and fractured subsurface rock formations, results from the spatially variable hydraulic resistance offered by a geological structure. In addition to geo-structural features such as network connectivity, geometry, and hydraulic resistance, geochemical reactions, for example, dissolution and precipitation, can dynamically inhibit or enhance flow channelization. These geochemical processes can change the fracture permeability leading to increased flow channelization, which are localized connected regions of high volumetric flow rates that are seemingly ubiquitous in the subsurface. In our simulations, fractures partially filled with quartz are gradually dissolved until quasi-steady state conditions are obtained. We compare the flow field's initial unreacted and final dissolved states in terms of flow and transport observations. We observe that the dissolved fracture networks provide less resistance to flow and exhibit increased flow channelization when compared to their unreacted counterparts. However, there is substantial variability in the magnitude of these changes which implies that the channelization strongly depends on the network structure. In turn, we identify the interplay between the particular network structure and the impact of geochemical dissolution on flow channelization. The presented results indicate that geological systems that have been weathering or reactive for longer times in older landscapes are likely to have increased flow channelization compared to their equivalent but younger counterparts, which implies a time dependence on flow channelization in fractured media.

Hyman, Jeffrey D.↗

DOE OSTI2024

Measuring Thread Timing to Assess the Feasibility of Early-Bird Message Delivery Across Systems and Scales

Early-bird communication is a communication/computation overlap technique that leverages fine-grained communication to improve application run-time. Communication is divided such that each individual thread can initiate transmission of its portion of the data upon completion rather than waiting for a dedicated communication phase. The benefit of early-bird communication depends on the completion timing of the individual threads: On the one hand, if all threads are complete at nearly the same time, the overheads of sending multiple messages will accumulate, leading to performance that is worse than if a single message had been sent. On the other hand, if thread completions are spread out in time, those that complete earlier can send data while others continue working, leading to performance that is better than if a single message had been sent. The challenge is that the completion times are currently unknown and can vary based on application, problem size, system software, and underlying hardware. In this paper, we address this lacuna by measuring and evaluating the potential overlap afforded by early-bird communication for a selection of proxy applications. These measurements help us understand whether a given application could benefit from early-bird communication. Here, we present our technique for gathering this data and evaluate data collected from three proxy applications: MiniFE, MiniMD, and MiniQMC. Each application is run on three systems with distinct CPU architectures and strong scales across three run sizes. To characterize the behavior of these workloads, we study the trends of thread timings at both a macro level, across all threads across all runs of an application, and a micro level, that is, within a single process of a single run. We observe that our tested applications exhibit significantly different thread arrival distributions. The machine used had a significant impact, with the window of potential overlap varying by as much as an order of magnitude.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Time-temperature-stress superposition of PBX 9502 compressive creep data for lifetime predictions

The plastic-bonded explosive (PBX) 9502 is a viscoelastic, high-solids loaded polymer bound composite comprised of 95 weight% (wt%) TATB explosive crystals and 5 wt% FK-800 polymer binder. The uniaxial quasi-static mechanical properties have been studied and characterized extensively over 25+ years, including creep. Creep is the strain evolution over time as a result of an asymmetric load typically at temperatures above ambient. Creep measurements to failure can take excessively long periods of time depending on various factors including applied stress, temperature, and loading rate that render a single test to failure impossible. Using the principles of time-temperature-stress superposition (TTSSP), short-term creep tests at various applied stresses and temperatures can be used for long-term creep predictions. The TTSSP model shifts short-term test data using two factors, one for temperature and stress, creating a master curve at a reference test condition. The master curve can be shifted to predict the creep response at any desired test condition. PBX 9502 compression specimens were tested at three temperatures and various applied stresses. This work discusses the TTSSP analysis of PBX 9502 compressive creep data to create a single comprehensive master curve.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Using machine-learning to understand complex microstructural effects on the mechanical behavior of Ti-6Al-4V alloys

Structural materials properties are highly dependent on their microstructure. Their microstructure is in turn affected by multiple fabrication and thermo-mechanical treatment parameters, all of which conform a highly-dimensional parametric space with often hidden correlations that are difficult to extract by experimentation alone. This is particularly true for alloys of the dual-phase Ti-6Al-4V family, with their greatly complex and rich microstructures, which combine several intrinsic length scales associated with multiple grain and subgrain structures, grains with different crystal lattices (α and β phases), and complex chemistry. In this paper we use a comprehensive set of machine learning techniques to develop predictive tools relating the yield strength and hardening rate of these alloys to a set of input parameters covering extensive ranges. The data generator is a finite-element crystal plasticity model for polycrystal deformation that takes into account slip anisotropy and employs standard dislocation evolution models for the α and β phases of Ti-based alloys. Our dataset includes over two thousand independent simulations and is used to train the machine learning models, which are then used to establish correlations between microstructural parameters and the alloys’ mechanical response. Our results point to the most influential parameters affecting yield strength and hardening rate, information that can then be used to guide experimental synthesis and characterization efforts to save time and resources.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Pressure-Stabilized MnSb2 with Complex Incommensurate Magnetic Order

Marcasite-type compounds have been proposed as promising hosts of exotic magnetic quantum states, yet experimental realizations in stoichiometric, disorder-free systems remain limited. Here, we report the high-pressure stabilization and magnetic characterization of MnSb2, a marcasite-type compound that is thermodynamically metastable under ambient pressure. Single crystals were synthesized using a cubic multianvil press at 3.3 GPa and 490 °C for 24 h, and powder and single-crystal X-ray diffraction confirm the orthorhombic Pnnm structure. These crystals are stable at ambient pressure for a long time up to between 450 and 500 K. Heat-capacity measurements reveal phase transitions at approximately T0 ∼ 118 K and T1 ∼ 220 K. Neutron diffraction uncovers an unconventional magnetic state below T1 ∼ 220 K. Magnetic powder neutron diffraction refinements reveal possible multiple magnetic configurations that provide comparably acceptable fits to the experimental data. While most solutions are consistent with a spin-density-wave (SDW) description, helical models systematically yield inferior agreement factors. Across a broad range of models, the Mn ordered moment reaches a maximum value of approximately 2 μB and remains predominantly collinear, with minimal canting along the c-axis. At 200 K, the magnetic propagation vector is q = (0, 0.3975, 0.3783); upon cooling, the b component increases toward 0.5, reflecting a temperature-dependent evolution of the modulation. The need for modification of the magnetic model between high and low temperatures further highlights the complex and strongly temperature-dependent nature of the magnetic order in this system. These results establish MnSb2 as a pressure-stabilized marcasite magnet with a tunable, complex magnetic state and a compelling stoichiometric platform for exploring unconventional magnetic behavior, including potential altermagnetism.

Xu, Mingyu [Michigan State University , , , ,; Iow↗

DOE OSTI2025

Modeling and Simulation of Fuel Dispersal During the Loss-of-Coolant Accident

This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

Oil price states and drivers: An analysis of the second-month spot-futures price differential

Oil remains a dominant component of global energy use, and its price, characterized by frequent changes and an ever-present potential for large swings, continues to be a focus of industry participants, policymakers and analysts attention. Here, this study examines the behavior of future spot oil prices using a continuous-time hidden Markov model (HMM) and daily price data spanning years 2007 to 2024. We identify six states in the second-month WTI spot-futures price differential and assess the roles of eleven futures price, quantity, financial market, and geopolitical/volatility variables in each state. The model yields several insights into the workings of the oil market and the relative roles of these drivers. We find support for several theoretical and empirical findings in the oil market literature, including the role of inventory, volatility/risk, and market responses to contango/backwardation in futures markets. A novel finding is that “normal contango” conditions represent a significant portion of next-day states in our in-sample data. Under the most volatile normal contango state, many of the oil market drivers differ markedly in coefficient signs and magnitudes from those in other states. The resulting model also performed well out-of-sample and would, in addition to understanding the impact of market drivers, be useful for short-term forecasting. Overall, the findings highlight the highly non-linear, regime-dependent interactions of the oil price and its drivers, emphasizing the importance of detailed information to market stakeholders.

Oladosu, Gbadebo A. [Oak Ridge National Laboratory↗