Search NASA⌕ Search

SEARCH · Search NASA

Results for “MD simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Molecular dynamics simulations of radiation response of $\mathrm{LiAlO_2}$ and $\mathrm{LiAl_5O_8}$

Here, we report the findings of a computational investigation of defect production and migration in two lithium aluminate ceramics: LiAlO 2 and LiAl 5 O 8 . Although the first ceramic, LiAlO 2 has been widely investigated computationally, especially via molecular dynamics (MD), there have been no MD studies done on LiAl 5 O 8 up to date, due to the lack of interatomic potentials for this system. This work implemented an existing set of Li-Al-O coulombic and Buckingham potential parameters for the LiAl 5 O 8 system from the literature, followed by validating it based on experimental values of density, crystallinity, bond lengths and melting point. The Li + diffusion in the two ceramics was examined and was found to be one order of magnitude slower in LiAl 5 O 8 (D 0 = 3.17 x 10 -11 m 2 /s) as compared to that in LiAlO 2 (D 0 = 4.02 x 10 -10 m 2 /s) at 600 K. The lithium vacancy migration barrier of 4.15 eV in LiAl 5 O 8 was more than three times that in LiAlO 2 (1.31 eV) possibly due to stronger ordering in the Li-O 6 octahedra of LiAl 5 O 8 than in Li-O 4 tetrahedra of LiAlO 2 . The Li displacement threshold energy (E d ) in LiAl 5 O 8 was found to have a median value of 68 eV which is much higher than that in LiAlO 2 (40 eV). The greater E d for Li signifies a lower likelihood for defect formation and clustering and thus lower tendency for amorphization in LiAl 5 O 8 . The presented results show the difference in defect dynamics in the two ceramics that may help us understand the susceptibility of the two ceramics to amorphization caused by irradiation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Insights into Sorption and Molecular Transport of Aqueous Glucose into Zeolite Nanopores

Liquid-phase heterogeneous catalysis using zeolites is important for biomass conversion to fuels and chemicals. There is a substantial body of work on gas-phase sorption in zeolites with different topologies; however, studies investigating the diffusion of complex molecules in liquid medium into zeolitic nanopores are scarce. Here, we present a molecular dynamics study to understand the sorption and diffusion of aqueous β-d-glucose into β-zeolite silicate at T = 395 K and P = 1 bar. Through 2-μs-long molecular dynamics trajectories, we reveal the role of the solvent, the kinetics of the pore filling, and the effect of the water model on these properties. We find that the glucose and water loading is a function of the initial glucose concentration. Although the glucose concentration increases monotonically with the initial glucose concentration, the water loading exhibits a nonmonotonic behavior. At the highest initial concentration (~20 wt %), we find that the equilibrium loading of glucose is approximately five molecules per unit cell and displays a weak dependence on the water model. Glucose molecules follow a single-file diffusion in the nanopores due to confinement. Here, the dynamics of glucose and water molecules slows significantly at the interface. The average residence time for glucose molecules is an order of magnitude larger than that in the bulk solution, while it is about twice as large for the water molecules. Our simulations reveal critical molecular details of the glucose molecule’s local environment inside the zeolite pore relevant to catalytic conversion of biomass to valuable chemicals.

09 BIOMASS FUELS↗
DOE OSTI2020

Performance and Feature Improvements in Parareal-based Power System Dynamic Simulation

In recent years, a novel Parareal-based approach has been developed for fast transient simulations of large power system interconnections. Parareal belongs to the class of Parallel-in-time algorithms for solution of systems of differential-algebraic equations in parallel over an interval of time. The selection of a reasonably fast and accurate coarse solution is crucial to improve the performance of Parareal algorithm. Semi-analytical solution methods are one promising approach to achieve this goal. They have been investigated, and some preliminary results are presented here. In addition, Parareal-based simulator has been expanded to enable co-simulation with OpenDSS, a widely used open-source distribution system simulator. Preserving the parallel nature of the Parareal approach and taking advantage of the parallel capabilities of the latest versions of OpenDSS, each distribution system can be solved in their entirety on different processors in parallel within the main Parareal simulator. This paper also presents the structure of the transmission and distribution co-simulation and some results with different dynamic models of inverter-based resources in the distribution systems.

Park, Byungkwon↗
DOE OSTI2022

Multi-level impacts of climate change and supply disruption events on a potato supply chain: An agent-based modeling approach

Context: The world is experiencing frequent extreme weather events like droughts, snowstorms, and shifting of seasons due to climate change. Increased frequency and severity of these extreme weather events threaten food security because agriculture depends on climate conditions. Impacts of climate change on the agricultural system not only occur at the grower’s level, but also cascade to other levels along the supply chain. Objective: In this study, we aim to quantify a wide range of economic impacts of different extreme climate events on different stages of a food supply chain. Methods: We chose the potato supply chain in Idaho as a case study. We developed a multi-echelon supply chain simulation model using an agent-based modeling (ABM) approach with five types of agents—farmers, shippers, processors, retailers, and logistics companies. In addition to the business-as-usual (BAU) scenario, we designed two climate-related disruption events—drought and snowstorm. We quantified the heterogeneous impacts at different stages of the supply chain for both fresh and processed potatoes using key performance indicators (KPIs) including revenues, prices, lead times, traded quantities, and food waste quantity. Results and Conclusion: The impacts of the disruption events are different on different agents in the supply chain for different product categories. The price hike of fresh potatoes is far higher than processed potatoes during disruption events. This price hike makes consumers switch to processed potatoes, which require more fresh potatoes as an input that further reinforces the price hike. However, because processed potatoes have an elastic demand, once their prices go up due to higher input cost, their demand drops. Non-contracted farmers gain additional revenues from the disruption events, whereas contracted farmers incur a loss due to lock-in price and lower than usual harvest. Significance: The methodology developed in this study could be applied to other food and agricultural supply chains for understanding the vulnerabilities at agent levels due to climate change disruption events. Furthermore, the findings would help develop mitigation strategies or policies to improve the well-being of the overall supply chain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Inferring assembly-curving trends of bacterial micro-compartment shell hexamers from crystal structure arrangements

Bacterial microcompartments (BMC) are complex macromolecular assemblies that participate in varied chemical processes in about one fourth of bacterial species. BMC-encapsulated enzymatic activities are segregated from other cell contents by means of semipermeable shells, justifying why BMC are viewed as prototype nano-reactors for biotechnological applications. Herein, we undertook a comparative study of bending propensities of BMC hexamers (BMC-H), the most abundant shell constituents. Published data show that some BMC-H, like β -carboxysomal CcmK, tend to assemble flat whereas other BMC-H often build curved objects. Inspection of available crystal structures presenting BMC-H in tiled arrangements permitted us to identify two major assembly modes with a striking connection with experimental trends. All-atom molecular dynamics (MD) supported that BMC-H bending is triggered robustly only from the arrangement adopted in crystals by BMC-H that experimentally form curved objects, leading to very similar arrangements to those found in structures of recomposed BMC shells. Simulations on triplets of planar-behaving hexamers, which were previously reconfigured to comply with such organization, confirmed that bending propensity is mostly defined by the precise lateral positioning of hexamers, rather than by BMC-H identity. Finally, an interfacial lysine was pinpointed as the most decisive residue in controlling PduA spontaneous curvature. Globally, results presented herein should contribute to improve our understanding of the variable mechanisms of biogenesis characterized for BMC, and of possible strategies to regulate BMC size and shape.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Quantum Algorithm for Linear Non-unitary Dynamics with Near-Optimal Dependence on All Parameters

We introduce a family of identities that express general linear non-unitary evolution operators as a linear combination of unitary evolution operators, each solving a Hamiltonian simulation problem. This formulation can exponentially enhance the accuracy of the recently introduced linear combination of Hamiltonian simulation (LCHS) method [An, Liu, and Lin, Physical Review Letters, 2023]. For the first time, this approach enables quantum algorithms to solve linear differential equations with both optimal state preparation cost and near-optimal scaling in matrix queries on all parameters.

Applied Dynamical Systems↗
DOE OSTI2023

A Scalable PDC Placement Technique for Fast and Resilient Monitoring of Large Power Grids

The wide-area measurement system (WAMS) is a key enabler of real-time monitoring of power grids. The essential goals of WAMS design are fast and resilient data transfer from phasor measurement units (PMU) to phasor data concentrators (PDC). We propose a scalable two-stage PDC placement technique for minimizing the end-to-end delay while maintaining resiliency. In the prescreening stage, the plausible candidates of PDC configurations are identified based on a graph theory-based multi-median function (MMF). Here, in this article, a computationally efficient meta-heuristic algorithm is used to address scalability. In the candidate selection stage, two different algorithms, namely, Suurballe's and Dijkstra's, are employed to identify the best of those plausible PDC configurations as the final design. This technique not only minimizes the hop paths between PMUs and PDCs, but also ensures network resiliency against single PMU, PDC, or communication link failure by incorporating the roles of PMUs in power grid observability into routing policy. Simulation results on the IEEE 57-bus test power system and the 2000-bus test power system demonstrate the effectiveness and scalability of the proposed technique.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Self-consistent modeling of tokamak edge plasma transport with lithium sources

Magnetic confinement fusion devices require effective heat and particle exhaust solutions on the divertor plates to operate sustainably, especially under reactor-relevant conditions. Liquid lithium divertors have been proposed to address two major challenges: control of excessive heat flux to plasma-facing components through vapor shielding and minimization of core plasma contamination from impurities. The National Spherical Torus Experiment-Upgrade (NSTX-U) will explore lithium as a divertor material due to its potential to meet both objectives. We present a self-consistent coupling framework between the plasma boundary transport code UEDGE and the lithium wall transport code Wall–Li to evaluate the feasibility and operational limits of lithium-based divertors. The model aims to optimize lithium sourcing levels to prevent core plasma contamination via fuel dilution while ensuring divertor protection through vapor shielding. This integrated framework, applicable to any tokamak with lithium sources, dynamically adjusts lithium sourcing based on local plasma conditions and surface temperature. The coupled model is tested using NSTX-like geometry and plasma conditions to assess its performance and reliability. Wall–Li calculates lithium fluxes from plasma-facing components, incorporating physical sputtering, thermally enhanced sputtering, and evaporation driven by surface temperature and ion flux. These fluxes are reintroduced into UEDGE as neutral lithium atoms, enabling simulation of their transport and distribution within the plasma. UEDGE computes plasma and neutral transport, surface heat flux, and iteratively feeds this information back to Wall–Li. A small time step is employed to ensure numerical stability and convergence, enabling accurate simulations over typical tokamak discharge durations. This integrated modeling approach provides a robust tool for identifying operational regimes that balance effective lithium sourcing with minimal core plasma contamination, offering critical insights for optimizing lithium-based divertor systems in current and future fusion devices.

Magnetic confinement fusion↗
DOE OSTI2025

Exponentially Reduced Circuit Depths Using Trotter Error Mitigation

Product formulas are a popular class of digital quantum simulation algorithms due to their conceptual simplicity, low overhead, and performance, which often exceeds theoretical expectations. Recently, Richardson extrapolation and polynomial interpolation have been proposed to mitigate the Trotter error incurred by the use of these formulas. This work provides a rigorous, general analysis of these techniques for computing time-evolved observables, simplifying the interpolation algorithm in the process, and shows that extrapolation generically improves the performance of product formulas for this task. We demonstrate that, to achieve error 𝜖 in a simulation of time 𝑇 using a 𝑝 ⁢th-order product formula with extrapolation, circuit depths of 𝑂⁡(𝑇 1+1/𝑝 ⁢polylog (1/𝜖)) are sufficient—an exponential improvement in the precision over product formulas alone. Furthermore, we prove that these algorithms achieve commutator scaling, and improve the 𝑇 complexity for the interpolation algorithm. By relaxing the requirement of performing exact Chebyshev interpolation, our simplified algorithm eliminates the need for fractional implementations of Trotter steps, reducing computational overhead. Finally, we show these techniques can be combined with the classical shadows method to estimate many time-evolved local observables. Taken together, our findings provide the strongest evidence yet for the utility of Trotter error-mitigation techniques in algorithmic applications.

quantum algorithms & computation↗
DOE OSTI2021

Calibration of the Diffusivity Predictions of Centipede Using Approximate Bayesian Computation and Applications in Nyx (Engineering Scale) and Xolotl-MARMOT (Meso-Scale) Simulations

Fission gas evolution and release in UO 2 nuclear fuel are important fuel performance metrics and occur in several distinct stages: 1) nucleation, growth and resolution of intra-granular bubbles, 2) diffusion to grain boundaries and 3) nucleation and growth of bubbles at grain boundaries, which eventually form a connected network (percolation) enabling release of gas from grain boundaries through connections to triple junctions, grain edges or free surfaces. The NE-SciDAC project is developing several computational tools to model this problem, which are connected in a hierarchical multi-scale framework. The information transfer in the multi-scale framework is a critical step that, in addition to best-estimates, should include uncertainty quantification. Despite taking a first-principles multi-scale approach, there is a need to perform parameter calibration to ensure consistency with available experimental data. In the present study, uncertainty quantification (UQ) and parameter calibration is demonstrated for one of the lower length scale codes in the multi-scale framework (Centipede) and then the results, including instances of the propagated uncertainties, are used in other codes within the framework, specifically Nyx and Xolotl-MARMOT. We calibrated the model parameters in Centipede, a computer code used to predict diffusivities of uranium (U) and xenon (Xe) in the context of the simulation of fission gas in uranium oxide (UO 2 ) nuclear fuel. The Centipede code depends on 183 parameters, all of which are subject to uncertainty. The three data sets used in our calibration effort are taken from the literature. This data is available as a set of measurements, including measurement errors. Our goal is to calibrate a statistical model that predicts both the value of the measurement and the uncertainty associated with the measurement. We perform a Bayesian calibration of the model parameters using a dedicated approximate Bayesian computation (ABC) likelihood function. To avoid excessive computational costs, we replace the expensive Centipede simulation code by a higher-order surrogate model, constructed using only the 9 most important parameters. These important parameters are identified by a preliminary global sensitivity analysis (GSA) study. Among the important parameters are T0 (the temperature at which UO 2 is perfectly stoichiometric) and Hf_pO2 (the temperature dependence of the oxygen (O) partial pressure) that should be considered as operating conditions to be estimated along with the other parameters. We consider two different cases: one where we define one set of these operating conditions for all data sets, and one where we define distinct operating condition parameters for each data set. The Xe diffusivities predicted by the latter case show distinct features that could not be observed in the former. Next, we use the diffusivity predictions by Centipede as input to Nyx, a reduced order fuel performance code focused on gas behavior alone, in order to estimate quantities associated with inter-granular bubble formation at conditions specified by the experiments. Finally, the diffusivities obtained from the calibrated Centipede runs were used in coupled Xolotl-MARMOT simulations of intra- and inter-granular gas evolution. The results are compared to simulations using the baseline diffusivities from Turnbull et al.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Self Configuring Digital Twin for Optimizing E-Waste Recycling

Electronic waste recycling industry needs a decision support tool for optimizing their processes to become cost competitive. We developed a software called CMAT, Comprehensive Manufacturing Assessment Tool. The aim of CMAT is to provide the e-waste recycling companies with a fully customizable decision support framework that analyzes the optimal supply chain configurations. The software optimizes the logistics operations, helps to identify the best recycling process configuration, and generates valuable insights regarding the economic performance of different categories of e-waste. The ultimate purpose of the tool is to assist the users developing a digital twin of their processes and to provide insights on questions pertinent to the e-waste recycling industry including how to increase efficiency and reduce costs, energy consumption, and greenhouse gas emissions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Fine-Tuning Microporosity of Crystalline Vanadomolybdate Frameworks for Selective Adsorptive Separation of Kr from Xe

Selective adsorptive capture and separation of chemically inert krypton (Kr) and xenon (Xe) noble gases with very low ppmv concentrations in air and industrial off-gases constitute an important technological challenge. Here, using a synergistic combination of experiment and theory, the microporous crystalline vanadomolybdates (MoVO x ) as highly selective Kr sorbents are studied in detail. By varying the Mo/V ratios, we show for the first time that their one-dimensional (1D) pores can be fine-tuned for the size-selective adsorption of Kr over the larger Xe with selectivities reaching >100. Using extensive electronic structure calculations and grand canonical Monte Carlo simulations, the competition between Kr uptake with CO 2 and N 2 was also investigated. As most materials reported so far are selective toward the larger, more polarizable Xe than Kr, this work constitutes an important step toward robust Kr-selective sorbent materials. Furthermore this work highlights the potential use of porous crystalline transition metal oxides as energy-efficient and selective noble gas capture sorbents for industrial applications.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Reactive flash sintering of TiZrN and TiAlN ternary metal nitrides

Here, this study demonstrated the reactive flash sintering (RFS) for two powder mixtures: TiN-ZrN (both conducting) and AlN-TiN (TiN conducting but AlN insulating), targeting ternary metal nitrides (TMN) of Ti 0.5 Zr 0.5 N and Ti 0.5 Al 0.5 N, respectively. A constant volage and pressure (e.g., 8 V DC, ∼15 MPa) at room temperature triggered the flash (current density up to 27 A/mm 2 ) without pre-heating, and the entire RFS process finished in a few minutes. For TiZrN, the flash was instantaneous whereas for TiAlN, there was a long incubation before the flash followed by a quick and dramatic flash. Both conventional ex situ XRD and in situ synchrotron study had been carried out. They showed a uniform Ti 0.57 Zr 0.43 N solution formed in RFS and persisted upon cooling, while (Ti, Al) N solid solution formed at high temperature was not stable and likely went through a very quick phase separation in the cooling process. The final products from RFS had been characterized using SEM/EDS for microstructure. Both TiZrN and TiAlN were dense. Distribution of Ti, Zr, and N was uniform for TiZrN; for TiAlN, Ti and N distribution was uniform, while association of Al with oxygen was observed. TGA-DSC revealed the onset oxidation temperature for TiZrN was comparable to TiN and ZrN, while it was higher by ∼200 °C for TiAlN, likely due to the formation of Al 2 O 3 . In terms of mechanical properties such as hardness or fracture toughness, forming a single-phase solid solution (like TiZrN) does not offer obvious benefits. while large grain size from RFS seemed to be unfavorable. Future optimization of RFS condition and in-depth study by both experiments and simulation are needed to fully understand the composition-processing-structure-property relationships for such TMN from the reactive flash sintering process.

Flash Sintering↗
DOE OSTI2025

Housing-Performance Atlas of Baltimore Row Homes: Archetype-Based Multi-Hazard Baseline of Energy, Heat, Survivability, and Durability

Baltimore’s historic row-home neighborhoods face escalating risks to energy, heat, and durability under intensifying climate stress. This study develops a Housing-Performance Atlas that quantifies multi-hazard performance for eight representative archetypes using DesignBuilder/EnergyPlus Version 7.3.1.003, under Baltimore TMY3 boundary conditions. Performance is evaluated across the following four adaptation domains: energy use intensity, passive survivability during 72 h outage events, roof overheating exposure (>150 °F exceedance hours), and material service life derived from ISO 15686 and synthesized into Lean and Full Deficit Indices for comparative resilience ranking. Results show that EUI ranged from 46.7 to 67.6 kBtu ft −2 ·yr −1 , survivability from 0 to 23 h, and roof temperatures exceeded 150 °F for 150–210 h, shortening roof service life by up to 10 years. Composite Lean and Full Deficit Indices ranged 7.8–92.4, ranking Model 5 (end-unit, flat roof, two-story with basement) as the most resilient configuration and Model 8 (end-unit, pitched roof, three-story above-grade) as the least resilient due to compounded overheating and energy losses. Heat-related domains accounted for nearly 70% of overall resilience deficits, confirming thermal safety and roof reflectivity as retrofit priorities. The Housing-Performance Atlas establishes a reproducible diagnostic framework linking simulation, service life, and resilience metrics to guide cost-effective, climate-responsive retrofits in Baltimore’s aging urban housing stock.

Housing-Performance atlas↗
DOE OSTI2021

Natural Charge-Transfer Analysis: Eliminating Spurious Charge-Transfer States in Time-Dependent Density Functional Theory via Diabatization, with Application to Projection-Based Embedding

For many types of vertical excitation energies, linear-response time-dependent density functional theory (LR-TDDFT) offers a useful degree of accuracy combined with unrivaled computational efficiency, although charge-transfer excitation energies are often systematically and dramatically underestimated, especially for large systems and those that contain explicit solvent. As a result, low energy electronic spectra of solution-phase chromophores often contain tens to hundreds of spurious charge-transfer states, making LR-TDDFT needlessly expensive in bulk solution. More nefariously, intensity borrowing by the low-energy charge-transfer states can affect intensities of the valence excitations even if those excitation energies are accurate. At higher excitation energies, it is difficult to distinguish spurious CT states from genuine charge-transfer-to-solvent (CTTS) excitations. In this work, we introduce an automated diabatization scheme that enables fast and effective screening of the CTTS acceptor space in bulk solution. Our procedure introduces the concept of “natural charge-transfer orbitals”, which provide a means to isolate characteristic pairs of orbitals that are most likely to participate in a CTTS excitation. The projection of these orbitals onto solvent-centered virtual orbitals provides a criterion for defining the most important solvent molecules in a given excitation. We apply this method to analyze an ab initio molecular dynamics (MD) trajectory of I-(aq) and report the lowest-energy CTTS band in the absorption spectrum. Our results are in excellent agreement with experimental measurements for bulk I-(aq), and only one-third of the water molecules in the I-(H2O)96 simulation cell need to be described with LR-TDDFT in order to obtain excitation energies that are converged to < 0.1 eV. The tools introduced herein will improve the accuracy, efficiency, and usability of LR-TDDFT in solution-phase environments

Carter-Fenk, Kevin D.↗
DOE OSTI2025

Covariance-Free Bifidelity Control Variates Importance Sampling for Rare Event Reliability Analysis

Multifidelity modeling has been steadily gaining attention as a tool to address the problem of exorbitant model evaluation costs that makes the estimation of failure probabilities a significant computational challenge for complex real-world problems, particularly when failure is a rare event. To implement multifidelity modeling, estimators that efficiently combine information from multiple models/sources are necessary. In past works, the variance reduction techniques of control variates (CV) and importance sampling (IS) have been leveraged for this task. In this paper, we present the CVIS framework—a creative take on a coupled CV and IS estimator for bifidelity reliability analysis. The framework addresses some of the practical challenges of the CV method by using an estimator for the control variate mean and sidestepping the need to estimate the covariance between the original estimator and the control variate through a clever choice for the tuning constant. Furthermore, the task of selecting an efficient IS distribution is also considered, with a view towards maximally leveraging the bifidelity structure and maintaining expressivity. Additionally, a diagnostic is provided that indicates both the efficiency of the algorithm as well as the relative predictive quality of the models utilized. Finally, the behavior and performance of the framework is explored through analytical and numerical examples.

Markov chain Monte Carlo↗
DOE OSTI2026

Techno-Economic Analysis of Recycling Strategies for Catalyst and Acid During Catalytic Graphitization

With the aim of meeting the urgent demand for active anode materials (AAM) in energy storage systems, bio-based graphite (biographite) emerges as an affordable solution to de-risk the turbulent supply chain of critical minerals. Anode grade biographite requires high crystallinity and purity, which can be achieved by catalytic graphitization with iron, followed by acid washing. Therefore, a well-conceived process integration that recycles catalyst can be the starting point to commercialization. This study evaluates closed-loop catalyst recovery, and byproducts valorization scenarios through a technoeconomic framework to help understand the scale-up potential of biographite. For the acid washing, three reactors in series meet the required biographite purity at 99.95%. Iron and acid recovery can reduce material consumption and waste generation by ~95%, albeit at the expense of ~80% increase in capital costs. Recovery scenarios present similar capital and operational expenses, yielding minimum selling prices (MSP) near $6 kg-1 of biographite. Monte Carlo methodology reveals that feedstock price accounts for ~60% of MSP variance, followed by plant capacity ~20%. The likelihood of reaching a competitive profit margin of 30% in the U.S. AAM market sits at 85% average for recovery scenarios, and 103% when iron oxide is sold as byproduct. Additionally, an IRR >= 15% can be achieved for half of Monte Carlo simulations, representing promising early-stage results. Biographite production offers a strategic pathway to stabilize the anode market beyond China by integrating established technologies for a scalable, economically viable, and sustainable process. The role of catalyst recovery and byproducts utilization is critical for advancing the biomaterials industry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Entropy Stable Conservative Flux Form Neural Networks

We propose an entropy-stable conservative flux form neural network (CFN) to predict the dynamics of unknown governing conservation laws. The design of the network is based on the entropy-stable, second-order, and non-oscillatory Kurganov-Tadmor (KT) scheme. The proposed entropy-stable CFN, hereafter referred to as ESCFN, uses slope limiting as a denoising mechanism, ensuring accurate predictions in both noisy and sparse observation environments, as well as in both smooth and discontinuous regions. Importantly, our method is designed to predict long term dynamics of the unknown conservation law exclusively from a short temporal window of observed data, that is, without oracle knowledge of the PDE or later-time solution profiles. Numerical experiments demonstrate that the ESCFN achieves both stability and conservation while maintaining accuracy over extended time domains, and successfully predicts shock propagation speeds in long-term simulations. Furthermore, it is also robust to both noisy and sparse data environments.

Hyperbolic conservation laws↗