Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reference Model 2”

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 OSTI2024

DECOVALEX-2023: Task F1 Final Report

DECOVALEX-2023 Task F is a comparison of models and methods for post-closure performance assessment (PA) of a deep geologic repository for radioactive waste. The general aims of Task F are to build confidence in the models, methods, and software used for PA and to stimulate additional research and development in PA methodologies. The task objectives are to motivate development of PA modelling skills and capabilities, to examine the influence of model choices on calculated repository performance, and to compare the uncertainties introduced by model choices to other sources of uncertainty. Task F involves no actual experiment or site. It is a PA modelling exercise that requires the conceptual development of hypothetical repository designs and geologic settings. Because three of the teams were interested in salt and the rest of the teams were interested in crystalline rock, Task F was split into two branches: Task F1 for crystalline rock and Task F2 for salt. This report is for Task F1, crystalline rock. Teams from seven countries (Canada, Czech Republic, Germany, Korea, Sweden, Taiwan, and United States) participated in Task F1. The teams worked together to define the features, events, and processes of the reference case repository and established a set of performance measures. In addition, they defined a set of benchmark problems designed to test and compare modelling capabilities for fracture flow and transport at different scales. The repository design and benchmark problems are documented in a Task Specification that evolved over time as the group honed the specifications. The benchmark problems verified that each team can aptly model flow and transport in fractured media in 1-, 2-, and 3-dimensions. Two general approaches were used for the 3-dimensional benchmarks: discrete fracture network (DFN) and equivalent continuous porous medium (ECPM). DFN modelling involves explicit meshing of each fracture while ECPM modelling aims to capture the effective porosity and directional permeability of each cell in a space-filling mesh as affected by intersecting fractures. In some models, a combination of the two is used, i.e., DFN for large known fractures and ECPM for the rest of the domain. Transport is solved by using either the advection-dispersion equation or particle tracking. Although some variation is observed among model breakthrough curves in the benchmark problems, there is strong agreement in breakthrough behaviour up to at least the 75 th percentile for all benchmarks. At the 90 th percentile, breakthrough results show larger differences, suggesting several models retain substantially higher fractions of tracer in regions of slower moving water. In addition to the flow and transport benchmarks, several teams completed the source term benchmark, verifying capabilities for modelling radionuclide decay and ingrowth, waste package breach, instant release fractions, fuel matrix degradation rates, and radionuclide solubility limitations. The reference case is conceptualized as a generic spent fuel repository at a depth of 450 m in fractured crystalline rock. The repository has 50 parallel backfilled drifts, each with 50 deposition holes 6 m apart. Each deposition hole contains a 4-PWR waste package and bentonite buffer. The rock domain is 5 km in length, 2 km in width, and 1 km in depth. It has 6 deterministic fractured deformation zones and a multitude of stochastic fractures. Teams generally used the ECPM approach for the entire rock or a hybrid approach in which the deterministic fracture zones are modelled with a DFN and the rest of the rock is modelled by ECPM. Of the reference case problems specified, only the results of the initial reference case problem are compared in this report. The initial problem focuses on transport from the deposition holes to the surface, i.e., it neglects waste package performance. Tracers are released at all waste package locations at time zero and tracked for their releases to the near field and ground surface. The water fluxes calculated at the ground surface entry and exit regions of the domain are similar for all models except for two that have considerably lower fluxes. For tracer transport, large differences are observed among models in the magnitude of tracer transported. Much of the difference appears to be due to how the repository is implemented and hence the different degrees of repository simplification. Models that exclude the drifts, buffer, and backfill from the domain tend to show greater release of tracers and radionuclides from the repository. The initial study presented here indicates that major differences in modelling important processes within the repository (e.g., diffusion through buffer and backfill) can produce broadly different release and transport results, especially when those processes are excluded. Even for the models that included all specified features, events, and processes, the results show significant differences and demonstrate the importance of examining multiple modelling approaches in performance assessment. The differences in results observed in this study are expected to motivate teams to either increase complexity in future versions of the reference case models or to improve methods to account for the effects of simplified features and processes. Either way, future improvements in these models are expected to produce results that more closely agree.

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

Enhanced HLW Glass Property-Composition Models – Phase 2

As a continuation of the Phase 1 HLW model development work as part of a phased approach to enhance and expand HLW glass property-composition models, additional data were collected for the HLW glass dataset and used to revise the models developed during Phase 1. These data are for glasses that are collectively referred to as ORP-2014 glasses. The sources of additional data for the Phase 2 work included: (i) actively designed glasses formulated to support various glass studies at VSL since the completion of Phase 1; and (ii) glasses from two statistically designed matrices to supplement the existing HLW compositional space. The Phase 2 matrix design was intended to augment and improve the coverage in the high aluminum region.

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

A model to evaluate 100-year energy mix scenarios to facilitate deep decarbonization in the southeastern United States

The Southeast United States consumes approximately one billion megawatt-hours of electricity annually; roughly two-thirds from carbon dioxide (CO 2 ) emitting sources. The balance is produced by non-CO 2 emitting sources: nuclear power, hydroelectric power, and other renewables. Approximately 40% of the total CO 2 emissions come from the electric grid. The CO 2 emitting sources, coal, natural gas, and petroleum, produce approximately 372 million metric tons of CO 2 annually. The rest is divided between the transportation sector (36%), the industrial sector (20%), the residential sector (3%), and the commercial sector (2%). An Energy Mix Modeling Analysis (EMMA) tool was developed to evaluate 100-year energy mix strategies to reduce CO 2 emissions in the southeast. Current energy sector data was gathered and used to establish a 2016 reference baseline. The spreadsheet-based calculation runs 100-year scenarios based on current nuclear plant expiration dates, assumed electrical demand changes from the grid, assumed renewable power increases and efficiency gains, and assumed rates of reducing coal generation and deployment of new nuclear reactors. Within the model, natural gas electrical generation is calculated to meet any demand not met by other sources. Thus, natural gas is viewed as a transitional energy source that produces less CO 2 than coal until non-CO 2 emitting sources can be brought online. The annual production of CO 2 and spent nuclear fuel and the natural gas consumed are calculated and summed. A progression of eight preliminary scenarios show that nuclear power can substantially reduce or eliminate demand for natural gas within 100 years if it is added at a rate of only 1000 MWe per year. Any increases in renewable energy or efficiency gains can offset the need for nuclear power. However, using nuclear power to reduce CO 2 will result in significantly more spent fuel. More efficient advanced reactors can only marginally reduce the amount of spent fuel generated in the next 100 years if they are assumed to be available beginning around 2040. Thus closing the nuclear fuel cycle to reduce nuclear spent fuel inventories should be considered. Future work includes the incorporation of economic features into the model and the extension of the evaluation to the industrial sector. It will also be necessary to identify suitable sites for additional reactors.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2015

Unified continuum damage model for matrix cracking in composite rotor blades

This paper deals with modeling of the first damage mode, matrix micro-cracking, in helicopter rotor/wind turbine blades and how this effects the overall cross-sectional stiffness. The helicopter/wind turbine rotor system operates in a highly dynamic and unsteady environment leading to severe vibratory loads present in the system. Repeated exposure to this loading condition can induce damage in the composite rotor blades. These rotor/turbine blades are generally made of fiber-reinforced laminated composites and exhibit various competing modes of damage such as matrix micro-cracking, delamination, and fiber breakage. There is a need to study the behavior of the composite rotor system under various key damage modes in composite materials for developing Structural Health Monitoring (SHM) system. Each blade is modeled as a beam based on geometrically non-linear 3-D elasticity theory. Each blade thus splits into 2-D analyzes of cross-sections and non-linear 1-D analyzes along the beam reference curves. Two different tools are used here for complete 3-D analysis: VABS for 2-D cross-sectional analysis and GEBT for 1-D beam analysis. The physically-based failure models for matrix in compression and tension loading are used in the present work. Matrix cracking is detected using two failure criterion: Matrix Failure in Compression and Matrix Failure in Tension which are based on the recovered field. A strain variable is set which drives the damage variable for matrix cracking and this damage variable is used to estimate the reduced cross-sectional stiffness. The matrix micro-cracking is performed in two different approaches: (i) Element-wise, and (ii) Node-wise. The procedure presented in this paper is implemented in VABS as matrix micro-cracking modeling module. Three examples are presented to investigate the matrix failure model which illustrate the effect of matrix cracking on cross-sectional stiffness by varying the applied cyclic load.

36 MATERIALS SCIENCE↗
DOE OSTI2016

TU-FG-BRB-03: Basis Vector Model Based Method for Proton Stopping Power Estimation From Experimental Dual Energy CT Data

Purpose: This work aims at reducing the uncertainty in proton stopping power (SP) estimation by a novel combination of a linear, separable basis vector model (BVM) for stopping power calculation (Med Phys 43:600) and a statistical, model-based dual-energy CT (DECT) image reconstruction algorithm (TMI 35:685). The method was applied to experimental data. Methods: BVM assumes the photon attenuation coefficients, electron densities, and mean excitation energies (I-values) of unknown materials can be approximated by a combination of the corresponding quantities of two reference materials. The DECT projection data for a phantom with 5 different known materials was collected on a Philips Brilliance scanner using two scans at 90 kVp and 140 kVp. The line integral alternating minimization (LIAM) algorithm was used to recover the two BVM coefficient images using the measured source spectra. The proton stopping powers are then estimated from the Bethe-Bloch equation using electron densities and I-values derived from the BVM coefficients. The proton stopping powers and proton ranges for the phantom materials estimated via our BVM based DECT method are compared to ICRU reference values and a post-processing DECT analysis (Yang PMB 55:1343) applied to vendorreconstructed images using the Torikoshi parametric fit model (tPFM). Results: For the phantom materials, the average stopping power estimations for 175 MeV protons derived from our method are within 1% of the ICRU reference values (except for Teflon with a 1.48% error), with an average standard deviation of 0.46% over pixels. The resultant proton ranges agree with the reference values within 2 mm. Conclusion: Our principled DECT iterative reconstruction algorithm, incorporating optimal beam hardening and scatter corrections, in conjunction with a simple linear BVM model, achieves more accurate and robust proton stopping power maps than the post-processing, nonlinear tPFM based DECT analysis applied to conventional reconstructions of low and high energy scans. Funding Support: NIH R01CA 75371; NCI grant R01 CA 149305.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Light right-handed smuons at the LHC: Natural dark matter and g μ − 2 in an unexplored realm of the phenomenological MSSM

We introduce an in-depth study of an unprobed phenomenological minimal supersymmetric Standard Model region offering a natural solution to dark matter. Since the 1980s this region has been referred to as the “bulk,” consisting of sub 200 GeV neutralinos and right-handed smuons. The bulk satisfies recent muon g μ − 2 measurements and sustains consistency with all presently operating SUSY experiments and LHC constraints. Initial ingress into the bulk will arrive soon via the LUX-ZEPLIN 1000-day experiment. The ATLAS Collaboration at the LHC has confirmed that observation of these light right-handed smuon events can occur at the ongoing LHC Run 3 and forthcoming High-Luminosity LHC. Moreover, the future FCC-ee and CEPC circular colliders should handily observe the events. Published by the American Physical Society 2025

Yin, Xiangwei (ORCID:0000000193823206)↗
DOE OSTI2026

Modeling CO 2 flow through faulted/fractured reservoirs using tEDFM in corner-point grids

The interest in underground CO 2 storage has increased significantly over the last decade because of the rising concern about global warming due to the growing levels of greenhouse gases in the atmosphere. Considering that CO 2 accounts for 80% of these greenhouse gases, carbon capture, utilization, and storage (CCUS) is regarded as one of the most direct approaches to achieving the net zero carbon target. Although CO 2 storage in deep saline aquifers and depleted gas reservoirs has been studied extensively, most studies use commercial simulators that model faults/fractures by simply modifying the transmissibility in the direction perpendicular to the fault surfaces. Here, this work shows that this simplistic approach ignores the accelerated flow in the directions parallel to the fault plane, leading to significantly higher leakage along the fault surface. To accurately model the flow of CO 2 in faulted reservoirs, we present the first transient embedded discrete fracture model for corner-point grids (tEDFM-CPG). By comparing the results of the tEDFM-CPG to high-resolution reference solutions, we show that this approach is accurate and efficient at predicting CO 2 flow in faulted/fractured reservoirs. Finally, this work presents the use of mixed reality (MR) to efficiently observe CO 2 gas migration in the interior of these corner-point grid systems.

25 ENERGY STORAGE↗
DOE OSTI2022

Quasielastic least-squares reverse time migration of PS reflections

Conventional analysis of amplitude variation with offset for elastic PS reflections is based on analytical reflection coefficients in a layered medium, and wave-equation-based PS migration is mainly used to produce a structural image. To overcome this problem, we have developed a least-squares reverse time migration (LSRTM) method for elastic PS reflections based on a quasielastic wave equation. Here, the quasielastic wave equation can accurately model PS reflections with elastic amplitudes under the first-order Born approximation. Our LSRTM method inverts for perturbations of the S-wave velocity and density by minimizing the L 2 norm of the difference between recorded and predicted PS reflections modeled using a quasielastic wave equation. We refer to our new method as quasielastic LSRTM of PS reflections. Numerical tests on synthetic and field data indicate that our method can properly handle the amplitudes of elastic PS reflections and provides an accurate estimate of the perturbations of S-wave velocity and density. Extending the method to the 3D case is not straightforward and might require incorporating certain PS data processing techniques into the inversion itself.

58 GEOSCIENCES↗
DOE OSTI2016

Marginalizing instrument systematics in HST WFC3 transit light curves

Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) infrared observations at 1.1–1.7 μm probe primarily the H{sub 2}O absorption band at 1.4 μm, and have provided low-resolution transmission spectra for a wide range of exoplanets. We present the application of marginalization based on Gibson to analyze exoplanet transit light curves obtained from HST WFC3 to better determine important transit parameters such as R{sub p}/R{sub *}, which are important for accurate detections of H{sub 2}O. We approximate the evidence, often referred to as the marginal likelihood, for a grid of systematic models using the Akaike Information Criterion. We then calculate the evidence-based weight assigned to each systematic model and use the information from all tested models to calculate the final marginalized transit parameters for both the band-integrated and spectroscopic light curves to construct the transmission spectrum. We find that a majority of the highest weight models contain a correction for a linear trend in time as well as corrections related to HST orbital phase. We additionally test the dependence on the shift in spectral wavelength position over the course of the observations and find that spectroscopic wavelength shifts δ{sub λ}(λ) best describe the associated systematic in the spectroscopic light curves for most targets while fast scan rate observations of bright targets require an additional level of processing to produce a robust transmission spectrum. The use of marginalization allows for transparent interpretation and understanding of the instrument and the impact of each systematic evaluated statistically for each data set, expanding the ability to make true and comprehensive comparisons between exoplanet atmospheres.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Dynamic Wind Loading on CSP Collectors Caused by Turbulent Wind Fluctuations: Insights from a 2-Year Field Campaign

Concentrating Solar Power (CSP) is a promising solar technology for electricity generation with thermal energy storage and with the additional benefit of industrial heat production. Wind loading on CSP collector structures, such as parabolic troughs or heliostats, is one of the primary drivers of their structural design costs. In particular, dynamic wind loading is a major source of uncertainty in the collector design process, which heavily relies on wind tunnel testing. In the field, the turbulent nature of the incoming wind creates fluctuating loads (support structure loads and resulting mirror deflections) on the collectors, with impacts on fatigue lifetime and optical performance. As is well known, wind tunnel tests cannot entirely reproduce the complex turbulent wind conditions typically observed at full-scale plants. To shed light on this topic, NREL initiated a field campaign at the operational Nevada Solar One (NSO) powerplant that uses parabolic troughs as solar collectors. The aim of the project is a detailed characterization of prevailing wind and turbulence conditions and resulting operational loads on parabolic troughs. We use the published 2-year dataset of high-resolution combined wind and structural loads measurements [1] to characterize the dynamic structural wind response. For quantifying dynamic wind loading, we apply the concept of admittance functions, which are spectral transfer functions that couple the turbulent wind to resulting structural loads (aerodynamic admittance), and to the structural response (mechanical admittance). In practice, aerodynamic admittance describes which turbulent eddy sizes are effective in creating structural loads. The mechanical admittance describes in which frequency ranges these loads are reinforced or dampened by the structure. While these functions are an established concept in civil engineering, their recent application to a single full-scale heliostat [2] proved their broader applicability to CSP collectors. Here, we present a characterization of admittance functions for full-scale parabolic trough collectors and show how wind characteristics (mean wind speed and direction, turbulent kinetic energy, turbulent length scales), the sun-tracking trough angle, and row position alter the admittance functions. Further, we study to which extent the admittance functions are universal for a specific trough geometry and how our findings compare to reported heliostat results. References [1] https://data.openei.org/submissions/5938. [2] Blume, K., Roger, M., and Pitz-Paal, R. 2023b. "Simplified analytical model to describe wind loads and wind-induced tracking deviations of heliostats." Solar Energy, 256, 96-109. https://doi.org/10.1016/j.solener.2023.03.055.

admittance functions↗
DOE OSTI2024

Consequences of plane-strain and plane-stress assumptions in fully coupled chemo-mechanical Li-ion battery models

In Li-ion battery research, it is common to simulate chemo-mechanical phenomena in reduced dimensions (e.g., 2-D) as opposed to fully resolve these complex physics in 3-D. It is common to assume either (1) the out-of-plane strain is negligible (commonly referred to as plane-strain), or (2) the out-of-plane stress is negligible (commonly referred to as plane-stress). However, there is typically little consideration as to the quantitative consequences of these approximations. Furthermore, the influence of these out-of-plane assumptions can be compounded and convoluted when chemo-mechanics models implement so-called “fully coupled” formulations, where the local species concentrations influence the stress-state and the stress-state influences the local species fluxes. Here, the present manuscript explores the implications of using plane-stress and plane-strain assumptions in 2-D as compared to simulating a full 3-D electrode particle. This comparative study includes simulating both isotropic and anisotropic particle intercalation where the particles can be surrounded by either a liquid or solid electrolyte. Additionally, common Li-ion battery-model metrics such as the state-of-stress, intercalation fraction distribution, and specific capacity are compared, while also considering the effects of particle size and C-rate. As alternatives to the pure plane-strain and plane stress approximations, two modified plane-strain assumptions are found to better approximate the fully coupled chemo-mechanical 3-D behavior.

25 ENERGY STORAGE↗
DOE OSTI2022

Linear relationships for modeling CO 2 absorption in aqueous alkanolamine solutions in a thermodynamically consistent way

A thermodynamically consistent model for the carbon dioxide (CO 2 ) absorption in aqueous alkanolamine system is of great importance in the research and development of a CO 2 capture process. To facilitate the development of thermodynamic models, linear Gibbs free energy, enthalpy, and heat capacity relationships using well-known amines as reference are used to correlate the standard reference state properties of ionic species with those of molecular species in the electrolyte system, which has been approved to provide a reliable and consistent way to estimate required parameters when there is minimal or no appropriate experimental data available. Here, the proposed relationships have been applied to the development of an electrolyte nonrandom two liquid (NRTL) activity coefficient model for CO 2 absorption in aqueous 1-amino-2-propanol (A2P) solution, as an example to demonstration the methodology. With limited vapor–liquid equilibrium data and other thermodynamic properties, the parameters in the electrolyte NRTL model are identified with good accuracy.

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

Automated ICRF heating surrogate modeling via machine learning

This work introduces automated machine learning workflows that address critical bottlenecks in surrogate model development for Ion Cyclotron Range of Frequencies (ICRF) heating applications. The automated framework includes data analysis tools that transform raw datasets into actionable insights in seconds, replacing weeks of manual exploratory effort and ensuring consistent, reproducible dataset characterization. By integrating advanced hyperparameter optimization (HPO) methods including Bayesian optimization via BoTorch and Tree-structured Parzen Estimators (TPE), the framework significantly reduces model development time from weeks to hours, decreasing computational cost and required expertise, while enabling high-accuracy surrogate models. Compared to traditional hyperparameter scanning (HPS) techniques such as methodical, randomized, and grid searches, HPO methods achieve superior convergence and predictive performance, even when compared to already well-tuned reference models. On NSTX High Harmonic Fast Wave (HHFW) heating datasets, both Random Forest Regressor (RFR) and neural network surrogates demonstrate improved accuracy, achieving R 2 values beyond 0.97 and 0.98, respectively. The results show that while HPO gains are modest for robust architectures like RFR, they become essential for more sensitive models such as neural networks, highlighting the trade-offs across optimization strategies. Through automated workflows that eliminate manual hyperparameter tuning and require minimal ML expertise, this work enables widespread adoption of high-fidelity surrogate models across the fusion community for real-time plasma control, uncertainty quantification, rapid experimental scenario development, and integrated system optimization.

Sanchez-Villar, Alvaro [Princeton Plasma Physics L↗
DOE OSTI2015

On the need for the light elements primary process (LEPP)

Extant chemical evolution models underestimate the galactic production of Sr, Y, and Zr as well as the solar system abundances of s-only isotopes with 90 < A < 130. To solve this problem, an additional (unknown) process has been invoked, the so-called light element primary process (LEPP). In this paper we investigate possible alternative solutions. Based on Full Network Stellar evolutionary calculations, we investigate the effects on the solar system s-only distribution induced by the inclusion of some commonly ignored physical processes (e.g., rotation) or by the variation of the treatment of convective overshoot, mass loss, and the efficiency of nuclear processes. Our main findings are (1) at the epoch of the formation of the solar system, our reference model produces supersolar abundances for the whole s-only distribution, even in the range 90 < A < 130, (2) within errors, the s-only distribution relative to {sup 150}Sm is flat, (3) the s-process contribution of the less massive AGB stars (M < 1.5 M {sub ⊙}) as well as of the more massive ones (M > 4.0 M {sub ⊙}) are negligible, (4) the inclusion of rotation implies a downward shift of the whole distribution with a higher efficiency for the heavy s-only isotopes, leading to a flatter s-only distribution, (5) different prescriptions on convection or mass loss produce nearly rigid shifts of the whole distribution. In summary, a variation of the standard paradigm of AGB nucleosynthesis would allow reconciliation of model predictions with solar system s-only abundances. Nonetheless, the LEPP cannot be definitely ruled out because of the uncertainties still affecting stellar and galactic chemical evolution models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Testing Lithium-Ion Battery with the Internal Reference Electrode: An Insight into the Blocking Effect

Long-term stability of reference electrode is a major issue in lithium-ion batteries. A composite Ni-Li reference electrode (RE) with a longer lifetime than traditional Cu-wire Li reference electrodes is proposed in this study. However, the larger size of the Ni-Li RE creates a blocking effect when positioned internally in the battery. The blocking effect manifests as an abnormal feature in the measured electrode potential. A two-dimensional (2D) finite element model of the cell with an internal RE is constructed to investigate the mechanism of the blocking effect. The internal RE blocks the path of local current flow in the electrolyte, and decreases the electrochemical reaction rate of particles in the anode and cathode, leading to measurement errors. Characteristic lengths of RE and applied currents were examined as parameters in the finite element model to determine their effect on the magnitude of error. A useful guide is provided for proper selection of size of RE and applied current in dynamic tests of lithium-ion batteries.

2-D finite element model↗
DOE OSTI2014

Accurate prediction of diradical chemistry from a single-reference density-matrix method: Model application to the bicyclobutane to gauche-1,3-butadiene isomerization

Multireference correlation in diradical molecules can be captured by a single-reference 2-electron reduced-density-matrix (2-RDM) calculation with only single and double excitations in the 2-RDM parametrization. The 2-RDM parametrization is determined by N-representability conditions that are non-perturbative in their treatment of the electron correlation. Conventional single-reference wave function methods cannot describe the entanglement within diradical molecules without employing triple- and potentially even higher-order excitations of the mean-field determinant. In the isomerization of bicyclobutane to gauche-1,3-butadiene the parametric 2-RDM (p2-RDM) method predicts that the diradical disrotatory transition state is 58.9 kcal/mol above bicyclobutane. This barrier is in agreement with previous multireference calculations as well as recent Monte Carlo and higher-order coupled cluster calculations. The p2-RDM method predicts the Nth natural-orbital occupation number of the transition state to be 0.635, revealing its diradical character. The optimized geometry from the p2-RDM method differs in important details from the complete-active-space self-consistent-field geometry used in many previous studies including the Monte Carlo calculation.

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

Perovskite Sorbent Oxygen Separation Modeling with MFiX

This document chronicles the development and implementation of computational kinetic rate models that capture absorption and desorption characteristics of the National Energy Technology Laboratory (NETL) developed perovskite, Sr 1-x Ca x FeO 3-δ . Two paths to create accurate kinetic rates were followed: (1) an isothermal rate approach where thermogravimetric (TGA) data are recast as oxygen capacities through a pseudo-second order Lagergren equation (He et al., 2009); and (2) a more traditional Arrhenius approach where experimental data are fit with a power law model to derive associate activation energies (Bulfin et al., 2020a). For reference, the mathematics and associate experimental strategies that support these derivations are included in this report. In addition, computational fluid dynamics (CFD) models were developed to utilize both kinetic rate formulations and applied to simulate oxygen uptake and release in small scale scenarios. The program Multiphase Flow with interphase eXchanges (MFiX) was used to create: (1) discrete element method (DEM) simulations of a single tube of granular perovskite experiencing isothermal O 2 -absorption; and desorption and (2) two-fluid-model (TFM) non-isothermal simulations of perovskite O 2 -absorption and desorption tubes that share a wall. Conjugate heat transfer between steel walled tubes and the perovskite bed were managed with user-defined functions. As the project moves to simulating larger scale devices that will require more robust conjugate heat transfer methods, developed kinetic rates and associate methodologies have been recast for use in the ANSYS Fluent CFD program.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Evaluation of RMI Ejecta Sourcing from Reconstructed Interfaces in 2D FLAG Flyer-Plate Simulations

This study investigates the performance of RMI ejecta sourcing from a reconstructed interface in two-dimensional FLAG ejecta flyer-plate simulations. Enabling the RMI ejecta model to operate on reconstructed interfaces has been a longstanding challenge in FLAG, limiting its applicability to experiments in which interface reconstruction is not required. Addressing this limitation is therefore an important step toward expanding the model’s predictive capability. The first phase of this effort focused on identifying a flyer-plate/target configuration that produces large ejecta clouds from two distinct shocks using the standard RMI ejecta model. A series of simulations examined the effects of flyer composition, layer thickness, and mesh resolution on shock structure and ejecta formation. VISAR traces were used to evaluate shock timing and identify configurations that produced a clear double-shock response. A copper middle-layer flyer with a doubled tin layer yielded the most favorable shock characteristics and was selected as the baseline configuration. After confirming this configuration performed as expected when running a reference FLAG version, it became the basis for three simulations: (1) a standard two-dimensional ejecta flyer-plate calculation using a reference FLAG version, (2) the same calculation with the reconstructed interface option enabled, and (3) a calculation employing ALE strategies to force the ejecta model to operate on the reconstructed interface. While the baseline simulation produced the expected ejecta response, the reconstructed interface approach did not result in measurable ejecta, indicating that this capability remains under development. These results highlight both the progress made and the remaining challenges in extending the RMI ejecta model to reconstructed interfaces.

36 MATERIALS SCIENCE↗