Search NASA⌕ Search

SEARCH · Search NASA

Results for “experimental comparison”

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

Modeling root system growth around obstacles

State-of-the-Art models of Root System Architecture (RSA) do not allow simulating root growth around rigid obstacles. Yet, the presence of obstacles can be highly disruptive to the root system. We grew wheat seedlings in sealed petri dishes without obstacle and in custom 3D-printed rhizoboxes containing obstacles. Time-lapse photography was used to reconstruct the wheat root morphology network. We used the reconstructed wheat root network without obstacle to calibrate an RSA model implemented in the R-SWMS software. The root network with obstacles allowed calibrating the parameters of a new function that models the influence of rigid obstacles on wheat root growth. Experimental results show that the presence of a rigid obstacle does not affect the growth rate of the wheat root axes, but that it does influence the root trajectory after the main axis has passed the obstacle. The growth recovery time, i.e. the time for the main root axis to recover its geotropism-driven growth, is proportional to the time during which the main axis grows along the obstacle. Qualitative and quantitative comparisons between experimental and numerical results show that the proposed model successfully simulates wheat RSA growth around obstacles. Our results suggest that wheat roots follow patterns that could inspire the design of adaptive engineering flow networks.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Burning velocities of R-32/O 2 /N 2 mixtures: Experimental measurements and development of a validated detailed chemical kinetic model

This work entails characterizing the flammability of the refrigerant R-32 (CH 2 F 2 ) by both experimental measurements and modeling. Burning velocities S u were measured using a constant-volume spherical-flame method for R-32/O2/N 2 mixtures with O 2 /N 2 ratios ranging from 21% (synthetic air) to 40%, pressures of (1 to 3) bar, and equivalence ratios $\phi$ of (0.8 to 1.3). Based on a critical assessment of available data, and extended by our own calculations, a detailed chemical kinetic model was developed and key reactions determined using reaction path and sensitivity analyses. Initiation and combustion were identified as distinct kinetic regimes and burning velocities were found to be controlled by two primary reactions: unimolecular decomposition of CH 2 F 2 → CHF + HF and the subsequent reaction, CHF + O 2 → CHFO + O, the latter reaction initiating the radical chain propagating and branching by producing O atoms. Sensitive rate constants in the kinetic model were critically adjusted within their uncertainties and current knowledge bounds to best fit the experimental burning velocities. We found that rate constants in the model could be adjusted to match a given experimental S u for specific conditions (O 2 loading, P, T, $\phi$). This, however, then fixes predicted burning velocities for other all conditions within (3 to 4)% if physically realistic rate parameters are maintained. Thus, the entire set of experimental data is fit, not just to particular conditions. Relative random uncertainties in the experimental Su measurements were (4 to 6)%, but assumptions made for thermal radiation lost by the burned gas in the spherical-flame experiments add an additional systematic uncertainty. Systematic differences between the limiting cases of adiabatic (no thermal radiation lost) and optically-thin (all thermal radiation lost) varied significantly with conditions and ranged from (4 to 30)% at high to low velocities, respectively, translating into uncertainties of (2 to 15)% considering the average of two limiting cases. Comparison of experimental and kinetically modeled Su values suggests that the burned gas tends towards the optically-thin limit at the lowest pressures and fuel loadings and toward the adiabatic limit at the highest pressures and loadings. We tested and found support for this conclusion with a detailed analysis as a function of all the conditions (T, P, % O 2 , $\phi$). This behavior appears to transition from optically-thin to adiabatic as the density of the initial fuel increases, which results in increased CO 2 in the burned gas and thus increased absorption of the thermal radiation (consistent with the Beer-Lambert Law). The validated detailed model based on evaluated kinetics is shown to accurately predict burning velocities for R-32 O 2 /N 2 mixtures over a wide range of conditions and provides a reliable basis for extrapolation to other conditions.

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

Phonon Spectroscopy in Antimony and Tellurium Oxides

α-Sb 2 O 3 (senarmontite), β-Sb 2 O 3 (valentinite), and α-TeO 2 (paratellurite) are compounds with pronounced stereochemically active Sb and Te lone pairs. The vibrational and lattice properties of each have been previously studied but often lead to incomplete or unreliable results due to modes being inactive in infrared or Raman spectroscopy. Here in this paper, we present a study of the relationship between bonding and lattice dynamics of these compounds. Mössbauer spectroscopy is used to study the structure of Sb in α-Sb 2 O 3 and β-Sb 2 O 3 , whereas the vibrational modes of Sb and Te for each oxide are investigated using nuclear inelastic scattering, and further information on O vibrational modes is obtained using inelastic neutron scattering. Additionally, vibrational frequencies obtained by density functional theory (DFT) calculations are compared with experimental results in order to assess the validity of the utilized functional. Good agreement was found between DFT-calculated and experimental density of phonon states with a 7% scaling factor. The Sb–O–Sb wagging mode of α-Sb 2 O 3 whose frequency was not clear in most previous studies is experimentally observed for the first time at ~340 cm -1 . Softer lattice vibrational modes occur in orthorhombic β-Sb 2 O 3 compared to cubic α-Sb 2 O 3 , indicating that the antimony bonds are weakened upon transforming from the molecular α phase to the layer-chained β structure. The resulting vibrational entropy increase of 0.45 ± 0.1 kB/Sb 2 O 3 at 880 K accounts for about half of the α–β transition entropy. The comparison of experimental and theoretical approaches presented here provides a detailed picture of the lattice dynamics in these oxides beyond the zone center and shows that the accuracy of DFT is sufficient for future calculations of similar material structures.

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

Nonadiabatic Derivative Couplings through Multiple Franck–Condon Modes Dictate the Energy Gap Law for Near and Short-Wave Infrared Dye Molecules

Near infrared (NIR, 700–1000 nm) and short-wave infrared (SWIR, 1000–2000 nm) dye molecules exhibit significant nonradiative decay rates from the first singlet excited state to the ground state. While these trends can be empirically explained by a simple energy gap law, detailed mechanisms of nearly universal behavior have remained unsettled for many cases. Theoretical and experimental results for two representative NIR/SWIR dye molecules reported here clarify the key mechanism for the observed energy gap law behavior. It is shown that the first derivative nonadiabatic coupling terms serve as major coupling pathways for nonadiabatic decay processes from the first excited singlet state to the ground state for these NIR and SWIR dye molecules and that vibrational modes other than the highest frequency modes also make significant contributions to the rate. Finally, this assessment is corroborated by further theoretical comparison with possible alternative mechanisms of intersystem crossing to triplet states and also by comparison with experimental data for deuterated molecules.

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

State of the art, gaps, and prospects in fusion materials theory and modelling

Advancing the theory and simulation of materials for fusion applications remains a key component of global roadmaps aimed at delivering much-needed fusion power. Especially as the drive for commercial application increases, prototypes must be designed against radiation damage before the relevant experimental data can be collected and cost reductions that are possible by testing materials in silico become even more important. Here, we summarise the state of the art as it emerged during the 7 th Fusion Materials Theory & Modelling Workshop that took place in 2024, with the aim to highlight present gaps and future directions for the fusion materials modelling community. Of particular interest were the effects of transmutations, chemical complexity with the development of novel alloys and interatomic potentials, advancements in modelling high-dose microstructures, comparison with experimental data and multiscale models for structural assessment relying on high-performance computing and virtual reality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A Rapid Compression Machine Study of 2-Phenylethanol Autoignition at Low-To-Intermediate Temperatures

To meet the increasing anti-knocking quality demand of boosted spark-ignition engines, fuel additives are considered an effective approach to tailor fuel properties for satisfying the performance requirements. Thus, screening/developing bio-derived fuel additives that are best-suited for advanced spark-ignition engines has become a significant task. 2-Phenylethanol (2-PE) is an attractive candidate that features high research octane number, high octane sensitivity, low vapor pressure, and high energy density. Recognizing that the low temperature autoignition chemistry of 2-PE is not well understood and the need for fundamental experimental data at engine-relevant conditions, rapid compression machine (RCM) experiments are therefore conducted herein to measure ignition delay times (IDTs) of 2-PE in air over a wide range of conditions to fill this fundamental void. These newly acquired IDT data at low-to-intermediated temperatures, equivalence ratios of 0.35–1.5, and compressed pressures of 10–40 bar are then used to validate the 2-PE model developed by Shankar et al. (2017). It is found that this literature model greatly overpredicts the current RCM data. The comparison of experimental and simulated results also provides insights into 2-PE autoignition behaviors at varying conditions. Further chemical kinetic analyses demonstrate that the absence of the O2-addition pathway of β-R. radical in the 2-PE model of Shankar et al. (2017) could account for the model discrepancies observed at low-to-intermediated temperatures.

Fang, Ruozhou (ORCID:0000000348315100)↗
DOE OSTI2020

Achieving a robust grassy-ELM operation regime in CFETR

We have identified a robust grassy-edge localized mode (ELM) operation regime for future tokamak reactors. The regime exists within a pedestal top electron collisionality (ν*) window at high global poloidal beta (β p ). The existence of an upper ν* limit for grassy-ELMs is consistent with results previously reported in experiments (Oyama et al 2010 Nucl. Fusion 50 064014), while the existence of a lower ν* limit has not been reported previously. Using EPED and BOUT++, a theoretical model that quantitatively explains the physics of the grassy-ELMs within the window, which distinguishes them from the small mixed-ELMs at lower ν*, is presented for the first time. A peeling-ballooning stability boundary is obtained by scanning the operating density space. The change in density corresponds to a change in ν* that affects the pedestal bootstrap current. High β p leads to a strong Shafranov shift, which affects the flux surface averaged pressure drive. The two effects combine to create a peeling-dominated window in intermediate ν* buffered by ballooning-dominated regimes. Only the peeling-dominated regime shows a cyclic behavior in the perturbed pressure during the nonlinear simulation of an ELM crash, reminiscent of grassy-ELM dynamics. Similarly, the energy released across the separatrix is demonstrated to be significantly smaller. The quick recovery of the ELM crash is explainable by the rapid rise of a low n kink-peeling instability when the pedestal current I ped exceeds a threshold at high β p . It minimizes the excursion beyond marginal stability and is absent in the ballooning-dominated regime. Comparison with recent experiments over a range of β p and ν* strongly supports the physical picture proposed by the modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A knowledge gap analysis for transient CHF prediction within RELAP5-3D

Here we compare several experimental pool and flow boiling critical heat flux experiments with corresponding computational models developed with RELAP5-3D. Encompassed within this paper are several sensitivity analyses for the boiling heat transfer correlations, critical heat flux (CHF) prediction within RELAP5-3D and the thermal properties of the cladding for experiments in the Transient Reactor Test Loop (TRTL) and Transient Reactor Test (TREAT) facilities. The sensitivity analyses were consistent in determining that there is a large contribution from the critical heat flux multiplier which in this application, accounts for the effects of power transient critical heat flux. Comparison to experimental results of both the TREAT CHF-SERTTA experimental campaign and a University of New Mexico (UNM) pool boiling experimental campaign demonstrate both the under-prediction of transient CHF and nucleate boiling heat transfer. To address the gaps in knowledge regarding the transient heating effect for both boiling heat transfer and CHF prediction, we present an experimental test matrix for the TRTL facility that will isolate the power transient effect under PWR relevant operational conditions. The test matrix will incorporate 6 power pulses with a constant energy deposition and a timescale that spans over several orders of magnitude, including a steady state test with an exponential period of 10 s.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Experiment-modeling studies comparing energy dissipation in the DIII-D SAS and SAS-VW divertors

Recent DIII-D experiments on Small Angle Slot (SAS) divertors have confirmed that a combination of divertor closure and target shaping can enhance cooling across the divertor target and increase energy dissipation, but with significant dependence on B T (toroidal magnetic field) direction. In these novel divertors, the roles of closure, target shaping, drifts, and scale lengths are all interconnected in optimizing dissipation, with the separatrix electron density n eSEP being the key parameter associated with the level of dissipation/detachment. After modifying the original flat-targeted graphite SAS to include a V shape with a tungsten coating on the outer side of the divertor (SAS-VW), matched series of discharges were run to compare to detailed SOLPS-ITER modeling. Experimentally, when run as designed with the outer strike point at the slot vertex, SAS-VW requires nearly identical n eSEP for detachment as the original SAS, with little difference in dissipation for the new geometry. This is in contrast to (1) earlier modeling predictions that a small change of the SAS geometry to a V shape should enhance dissipation at the same n eSEP for magnetic configurations having better H-mode access (ion B × ∇B drift directed into the divertor), and (2) despite the achievement of significantly higher (2-7x) neutral pressures and compression in the SAS-VW slot. Comparisons of experimental density scans to the most recent SOLPS-ITER modeling with ExB drifts show reasonable agreement for dissipation/detachment onset when using separatrix density as the independent parameter. In order to help understand the discrepancy in modeled vs actual performance for the new configuration, additional measurements varying gas injection location and impurity injection were undertaken. In-slot D 2 gas fueling is more effective (5–22 %) in promoting detachment, in accord with modeling. In-slot impurity injection (N 2 or Ne) can yield 30 % lower core Z eff and 15 % less confinement degradation after detachment compared to main chamber puffing, as well as relatively lower tungsten leakage from the divertor. Modeling can also reproduce the improved detachment seen as the strike point moves inboard of the slot vertex. While we can explain the effects of the most important parameters causing energy dissipation in these slot divertors, it remains that many aspects of their behavior cannot be accurately modeled using state-of-art codes such as SOLPS-ITER. This is of concern for future model-driven designs utilizing similar V-shaped geometries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

FY24 Progress Report on Viscosity and Thermal Conductivity Measurements of Nuclear Industry Relevant Chloride Salts: An Experimental and Computational Study

As presented in this report, experimental and computational techniques were performed to assess the viscosity and thermal conductivity of key alkali and actinide chloride mixtures for molten salt reactor developers. These mixtures were pure LiCl, NaCl-KCl, LiCl-NaCl, LiCl-KCl, LiCl-NaCl-KCl, and NaCl-UCl 3 . Experimental measurements of viscosity were performed with a rolling ball viscometer, whereas experimental measurements of thermal conductivity were performed with a variable gap apparatus. Additional benchmarking work was performed using both property measurement systems to prepare for x-ray radiography in stainless-steel crucibles for viscosity and to ensure that calibration methods were accurate for thermal conductivity before assessing the NaCl-UCl 3 system. Validation data for the NaCl-UCl 3 in literature are minimal. Details on the calibration methods, salt measurement processes, and sources of error and uncertainty are discussed in detail for both property measurements. The computational methods described herein involved ab-initio molecular dynamics (AIMD) calculations using CP2K. The calculations were performed for the LiCl-KCl-NaCl and NaCl-UCl 3 systems. These calculations not only provided thermophysical property estimations for comparison to experimental data, but they also allowed for the determination of diffusion coefficients, coordination numbers, and radial distribution functions to provide insight into ion mobility and local coordination environments, which is linked to macroscopic property trends.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Evaluation of Reduced Order Model for HT-9 Creep and Modifications to Current HT-9 Creep Model in BISON

In order to leverage existing reduced order models (ROM) for the modeling of LANL-developed HT-9 cladding constitutive behavior, mechanistic-based models have been implemented into BISON. It is posited that the increased fidelity of the ROM will show improvement over engineering-scale models in terms of comparison against experimental measurements. INL will assist LANL in implementing an HT9 ROM in BISON. Once the HT-9 mechanistic constitutive model can be leveraged in nuclear performance simulations, the numerical results of fast reactor models will be analyzed and compared against the Fuels Irradiation & Physics Database (FIPD) and separate effects tests to attempt a validation process. INL and LANL will work together to perform any necessary improvement that is identified during the implementation and use processes. Successful completion of the milestone will enable the general use of validated mechanistic ROMs for HT-9 cladding in BISON. This report provides a brief introduction to the Los Alamos Reduced Order Model Applied to Nonlinear Constitutive Equations (LAROMANCE) code, mechanical testing used to illustrate its creep and plastic deformation predictions for HT-9 cladding, and a comparison with existing models implemented into BISON. Current models have been modified to more appropriately account for primary thermal creep and results for corresponding testing are included in this report.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Performance Characterization of a Natural Gas–Air Rotating Detonation Engine

An experimental study of a rotating detonation engine (RDE) operating with natural gas and air at elevated chamber pressures and air preheat temperatures was conducted to quantify its performance at conditions representative of land-based power generation gas turbine engines. Here, the thrust produced by the combustor was measured to characterize its work output potential. High-frequency pressure transducers and broadband chemiluminescence measurements of the flame provided information about the wave structure and dynamics. Analysis of common performance metrics demonstrated the necessity of normalizing any RDE performance parameter by the driving system potential, typically the reactant manifold pressure. Application of a thermodynamic performance model to a generic RDE identified the area ratio between the RDE exhaust and injection throats as the primary parameter affecting delivered pressure gain. The model was further applied to draw comparison with experimental measurements of net pressure gain for identical flow conditions. Only one of the two tested injector configurations followed the predicted trends, suggesting that performance of the second was governed by physical processes other than the reactant thermodynamics. Although an absolute pressure gain was not demonstrated, it is promising that the natural gas–air RDE delivered up to 90% of the theoretical performance.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2021

Enabling Particulate Materials Processing Science for High-Consequence, Small-Lot Precision Manufacturing

This Laboratory Directed Research and Development project developed and applied closely coupled experimental and computational tools to investigate powder compaction across multiple length scales. The primary motivation for this work is to provide connections between powder feedstock characteristics, processing conditions, and powder pellet properties in the context of powder-based energetic components manufacturing. We have focused our efforts on multicrystalline cellulose, a molecular crystalline surrogate material that is mechanically similar to several energetic materials of interest, but provides several advantages for fundamental investigations. We report extensive experimental characterization ranging in length scale from nanometers to macroscopic, bulk behavior. Experiments included nanoindentation of well-controlled, micron-scale pillar geometries milled into the surface of individual particles, single-particle crushing experiments, in-situ optical and computed tomography imaging of the compaction of multiple particles in different geometries, and bulk powder compaction. In order to capture the large plastic deformation and fracture of particles in computational models, we have advanced two distinct meshfree Lagrangian simulation techniques: 1.) bonded particle methods, which extend existing discrete element method capabilities in the Sandia-developed , open-source LAMMPS code to capture particle deformation and fracture and 2.) extensions of peridynamics for application to mesoscale powder compaction, including a novel material model that includes plasticity and creep. We have demonstrated both methods for simulations of single-particle crushing as well as mesoscale multi-particle compaction, with favorable comparisons to experimental data. We have used small-scale, mechanical characterization data to inform material models, and in-situ imaging of mesoscale particle structures to provide initial conditions for simulations. Both mesostructure porosity characteristics and overall stress-strain behavior were found to be in good agreement between simulations and experiments. We have thus demonstrated a novel multi-scale, closely coupled experimental and computational approach to the study of powder compaction. This enables a wide range of possible investigations into feedstock-process-structure relationships in powder-based materials, with immediate applications in energetic component manufacturing, as well as other particle-based components and processes.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Measurement of flow properties coupled to experimental and numerical analyses of dense, granular flows for solar thermal energy storage

Granular flows of sintered bauxite proppants were researched along an inclined plane for solar thermal energy storage applications. Granular flow properties needed to drive numeric granular models were measured for improved numerical model predictions for Carbobead CP 50/140, 40/100, and 30/60 particles. Particle shape and size distributions were determined by coupling optical microscopy to an in-house image processing algorithm. The impulse excitation technique was used to measure elastic and shear moduli, and compute Poisson’s ratio. The coefficient of static sliding friction was measured using the slip-stick method, and the static rolling friction was determined from measured shear on particles positioned between two hot-pressed plates. The coefficient of restitution was measured by dropping particles on a surface and determining the kinetic energy before and after impact with the surface using high resolution particle tracking velocimetry. Particle size did not significantly impact the coefficients of restitution and static rolling friction, however, particle shape distribution resulted in a large variation in measurements. An inclined flow experiment was performed to characterize granular flows of Carbobead CP 30/60 particles using particle image velocimetry. Numerical models of the experiment using discrete element method were generated with the measured mechanical properties as inputs for comparison with experimental results. A constant directional torque rolling friction model best predicted bulk granular flow behavior. Furthermore, good agreement between the model and experiment was achieved at ambient, steady state conditions, with average velocity differences <10%.

14 SOLAR ENERGY↗
DOE OSTI2022

Deciphering Distinct Overpotential-Dependent Pathways for Electrochemical CO 2 Reduction Catalyzed by an Iron–Terpyridine Complex

[Fe(tpyPY2Me)] 2+ ([Fe] 2+ ) is a homogeneous electrocatalyst for converting CO 2 into CO featuring low overpotentials of <100 mV, near-unity selectivity, and high activity with turnover frequencies faster than 100 000 s -1 . To identify the origins of its exceptional performance and inform future catalyst design, we report a combined computational and experimental study that establishes two distinct mechanistic pathways for electrochemical CO 2 reduction catalyzed by [Fe] 2+ as a function of applied overpotential. Electrochemical data shows the formation of two catalytic regimes at low (η TOF/2 of 160 mV) and high (η TOF/2 of 590 mV) overpotential plateaus. We propose that at low overpotentials [Fe] 2+ undergoes a two-electron reduction, two-proton-transfer mechanism (electrochemical-electrochemical-chemical-chemical, EECC), where turnover occurs through the dicationic iron complex, [Fe] 2+ . Computational analysis supports the importance of the singlet ground-state electronic structure for CO 2 binding and that the rate-limiting step is the second protonation in this low-overpotential regime. When more negative potentials are applied, an additional electron-transfer event occurs through either a stepwise or proton-coupled electron-transfer (PCET) pathway, enabling catalytic turnover from the monocationic iron complex ([Fe] + ) via an electrochemical-chemical-electrochemical-chemical (ECEC) mechanism. Comparison of experimental kinetic data obtained from variable controlled potential electrolysis (CPE) experiments with direct product detection with calculated rates obtained from the energetic span model supports the PCET pathway as the most likely mechanism. Moreover, we build upon this mechanistic understanding to propose the design of an improved ligand framework that is predicted to stabilize the key transition states identified in our study and explore their electronic structures using an energy decomposition analysis. Finally, taken together, this work highlights the value of synergistic computational/experimental approaches to decipher mechanisms of new electrocatalysts and direct the rational design of improved platforms.

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

Phase field modeling of shearing processes of a dual-lobed γ"|γ'|γ" coprecipitate

The Ni-based superalloy IN718 is one of the most widely used commercial alloy in the aerospace industry since its development in 1960s. The excellent mechanical properties have been attributed in a great deal to strengthening by coherent γ' and γ" precipitates. The deformation mechanisms of these two phases have been well studied individually. Recent experimental characterization has shown coprecipitates of these two phases in a variety of morphologies and it was argued that these coprecipitates may lead to improved strengthening as compared to their monolithic counterparts. However, the deformation mechanisms of these coprecipitates are still not well understood. In this study, we performed microscopic phase field simulations, with generalized-stacking-fault (GSF) energy surfaces from ab initio calculations as inputs, to systematically study the shearing processes of a periodical array of dual-lobed coprecipitates as well as monolithic precipitates. We found that the coupling between the γ' and γ" phases in the coprecipitates forces dislocations to take high energy shearing pathways in both phases that would not occur if they were in monolithic forms. The coupling also creates stacking fault configurations in the coprecipitates that require high stress to form. Thus, the presence of coprecipitates in general should increase the resistance to dislocation shearing and lead to higher strength levels. Various fault configurations observed during the shearing process are documented as a reference for future comparison with experimental observations. The link between stacking fault shearing and microtwinning is also discussed. Here, the mechanisms analyzed in this study deepens our understanding of coprecipitation effects on alloy strength and may form a cornerstone for multi-precipitate strengthened alloy design strategies.

36 MATERIALS SCIENCE↗
DOE OSTI2025

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Efficient end-to-end simulation of time-dependent coherent X-ray scattering experiments

Physical optics simulations for beamlines and experiments allow users to test experiment feasibility and optimize beamline settings ahead of beam time in order to optimize valuable beam time at synchrotron light sources like NSLS-II. Further, such simulations also help to develop and test experimental data processing methods and software in advance. The Synchrotron Radiation Workshop ( SRW ) software package supports such complex simulations. We demonstrate how recent developments in SRW significantly improve the efficiency of physical optics simulations, such as end-to-end simulations of time-dependent X-ray photon correlation spectroscopy experiments with partially coherent undulator radiation (UR). The molecular dynamics simulation code LAMMPS was chosen to model the sample: a solution of silica nanoparticles in water at room temperature. Real-space distributions of nanoparticles produced by LAMMPS were imported into SRW and used to simulate scattering patterns of partially coherent hard X-ray UR from such a sample at the detector. The partially coherent UR illuminating the sample can be represented by a set of orthogonal coherent modes obtained by simulation of emission and propagation of this radiation through the coherent hard X-ray (CHX) scattering beamline followed by a coherent-mode decomposition. GPU acceleration is added for several key functions of SRW used in propagation from sample to detector, further improving the speed of the calculations. The accuracy of this simulation is benchmarked by comparison with experimental data.

36 MATERIALS SCIENCE↗