Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation types and techniques”

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

High-Frequency Signature-Based Fault Detection for Future MV Distribution Grids: Preprint

Increased deployment of Distributed Energy Resources (DERs) and microgrids are significantly changing the operation of distribution grid. Increasing penetration of power electronics based DERs also raise several operational problems such as protection coordination issues, voltage raise, recloser synchronization, stability, and frequency control. Legacy systems are designed for unidirectional power flow. Bidirectional power flows, and varying short circuit levels are some of the major challenges in implementing these legacy protection schemes for future power grid with high penetration of power electronics based DERs. Inverter based DER’s can contribute to a maximum of 2 p.u. fault current. Interconnection standards might require even lower fault current contribution. If these protection related challenges are not addressed, it will result in limiting the penetration of DERs in distribution feeders. In this paper, the feasibility of high frequency signature generated by transients is evaluated to identify and locate faults in a distribution system. The paper will discuss the need for an accurate line modeling approach to accurately capture the high frequency signatures. These high frequency signatures are naturally occurring electromagnetic transients for events such as switching events in the network, faults events, etc. This is different than that implemented in transmission systems as the proposed approach utilizes advanced signal processing techniques to retrieve valuable information from the high frequency signatures generated during faults and does not just rely on time delay measurements between reflections of the high frequency signals. In the full paper, information of a modified IEEE distribution feeder test system to accurately model the high frequency signatures will be provided. The distribution lines will be modeled using universal line modeling (ULM) technique in an electromagnetic transients program (EMTP) modeling platform. Different test cases will be simulated and these test cases will include different fault types, multiple fault locations, capacitor bank switching, etc. The feeder system modeled for these simulations will also include power electronics based DERs. Additionally, experimental results demonstrating the high frequency signatures generated due to a transient event in a medium voltage (13.8 kV) distribution line will also be presented in the full paper. Furthermore, the analysis combined with simulation and experimental results in the full paper will help demonstrate the feasibility and efficacy of the proposed protection scheme. The developed protection scheme will not just help solve problems generated by bidirectional current flow, varying short-circuit current capacity, but also help in proliferation of DERs in distribution systems.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2025

Designing a Framework for Solving Multiobjective Simulation Optimization Problems

Multiobjective simulation optimization (MOSO) problems are optimization problems with multiple conflicting objectives, where evaluation of at least one of the objectives depends on a black-box numerical code or real-world experiment, which we refer to as a simulation. Whereas an extensive body of research is dedicated to developing new algorithms and methods for solving these and related problems, it is challenging and time-consuming to integrate these techniques into real-world production-ready solvers. This is partly because of the diversity and complexity of modern state-of-the-art MOSO algorithms and methods and partly because of the complexity and specificity of many real-world problems and their corresponding computing environments. The complexity of this problem is only compounded when introducing potentially complex and/or domain-specific surrogate-modeling techniques, problem formulations, design spaces, and data acquisition functions. Here, this paper carefully surveys the current state of the art in MOSO algorithms, techniques, and solvers, as well as problem types and computational environments where MOSO is commonly applied. We then present several key challenges in the design of a parallel multiobjective simulation optimization framework (ParMOO) and how they have been addressed. Finally, we provide two case studies demonstrating how customized ParMOO solvers can be quickly built and deployed to solve real-world MOSO problems.

engineering design optimization↗
DOE OSTI2022

Bulk properties of honeycomb lattices of superconducting microwave resonators

We have realized different honeycomb lattices for microwave photons in the 4 to 8 GHz band using superconducting spiral resonators. Each lattice comprises a few hundred sites. Two designs have been studied, one leading to two bands touching at the Dirac points and one where a gap opens at the Dirac points. Using a scanning laser technique to image the eigenmodes of this new type of photonic lattices, we are able to reconstruct their band structure. The measured bands are in excellent agreement with ab initio models that combine numerical simulations of the electromagnetic properties of the spiral resonator and analytical calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Multivariate Testing of Sampling Techniques to Address Class Imbalance in Building Use Type Classification

This study addresses the challenges inherent in building use type classification, particularly focusing on the issue of class imbalance in the training datasets for machine learning classifiers. We comprehensively analyze the efficacy of various class-balancing sampling techniques. Employing Monte Carlo simulations and Bayesian optimization, we evaluated the performance of multiple sampling methods, including Random Oversampling, Random Undersampling, SMOTE, Borderline-SMOTE, and ADASYN, across a dataset encompassing nine southeastern coastal states of the United States. Our findings reveal that simple random over- and undersampling techniques outperform more sophisticated methods. Additionally, we show inherent value in creating an imbalance in training data to effectively train a machine learning classifier for distinguishing between residential and nonresidential buildings. This study provides valuable guidance for future research on building use type classification research and lays essential groundwork for developing attribute-rich building stock datasets.

Adams, Daniel↗
DOE OSTI2021

Improving energy yield in photovoltaic modules with photonic structures

When installed outdoors, Si solar cells typically operate 20 – 30 K above ambient conditions. These elevated operating temperatures lead to diminished performance, reducing efficiency by ~0.4%/K for crystalline Si cells and decreasing module lifetime. Since much of the elevated temperature of the modules derives from parasitic absorption of sub-bandgap sunlight by the contacts, encapsulants, or other materials, reflectors that remove the sub-bandgap radiation from the module before it is absorbed promise to reduce the operating temperature. However, there are many competing strategies to reject this light: selective reflectors could be integrated into the outer module glass, at the textured Si-encapsulant interface, or on the rear. Furthermore, many of these approaches are complex and expensive to fabricate. The issue of heating can be even more critical in bifacial modules: high fractions of rear irradiance can lead to higher operating temperatures. Prior to this project, relatively little was known about the optical strategies for temperature reduction in bifacial modules. The goal of this project was to quantitatively assess optical strategies for improving the energy yield via a combination of improved anti-reflection and sub-bandgap light rejection. We studied two complementary concepts: the performance of ideal structures to determine the upper limits to performance, and low-complexity structures made from standard materials as cost-effective implementations. We also studied the performance of these reflectors in different types of modules, including Al BSF, PERC, and bifacial. Here we build from our prior accomplishments in developing mirror optimization techniques, and make use of simulation methods that accurately capture the optical properties of the modules throughout the visible and near-infrared spectrum. We connect these optical properties to energy yield under realistic, outdoor conditions. We also experimentally fabricated mirrors, incorporated them into mini-modules, and performed outdoor testing of performance at NREL. We found definitively that spectrally-selective mirrors on the outside of the module cover glass offer the best performance, both under ideal conditions and for realistic, low- complexity mirrors. Mirrors at the Si cell/encapsulant interface are limited by the multiple reflections needed to remove sub-bandgap light. Mirrors at the cell surface and mirrors on the rear also cannot remove all the sub-bandgap light, since parasitic absorption also occurs in the cover glass, encapsulant, and front contacts, depending on the module configuration. Also, cell texturing places higher demands on any cell-surface mirror because multiple interactions with the mirror are required to reject sub bandgap light which compounds losses. For bifacial modules, mirrors on the back do provide additional cooling and energy yield improvements, and the mirrors can have the same structure. In the monofacial case up to 3.3K of cooling is possible for mirrors on the front, compared to 2.2K at the cell surface and 1.2K at the rear. For bifacial modules, up to 2.4K of cooling is possible when mirrors are included at both the front and back. Our realistic designs utilize 4 – 6 layers of common photovoltaic materials such as SiO 2 , TiO 2 , and SiN x , making them more cost effective. Finally, we used our models to study thermal management in bifacial modules with different back lamination materials, finding that glass-glass laminated modules operated hotter than equivalent glass-polymer backsheet modules, due solely to the increased absorption of light incident from the rear of the module. Improving the energy yield by 2%, which is achievable with these realistic designs, would lead to LCOE reductions of $0.03/kWh by 2030. These benefits make photovoltaic panels more accessible to the public and may reduce associated maintenance costs.

14 SOLAR ENERGY↗
DOE OSTI2023

Connecting Large‐Scale Meteorological Patterns to Extratropical Cyclones in CMIP6 Climate Models Using Self‐Organizing Maps

Abstract Extratropical cyclones (ETCs) are responsible for the majority of cool‐season extreme events in the northeastern United States (NEUS), often leading to high‐impact weather conditions that can have wide‐ranging socioeconomic impacts. Evaluating the ability of climate models to adequately simulate ETC dynamics is essential for improving model performance and increasing confidence in future projections used by stakeholders and policymakers. ETCs are traditionally studied using techniques such as case studies and synoptic typing, however, these approaches can be time‐consuming, require subjective analysis, and do not necessarily identify the coincident large‐scale meteorological patterns (LSMPs). Here, we apply self‐organizing maps (SOMs) as an automated machine‐learning approach to characterize the LSMPs and associated frequency and intensity of discrete ETC events over NEUS. The dominant patterns of geopotential height variability are identified through SOM analysis of five reanalysis products during the last four decades. ETC events are tracked using TempestExtremes and are integrated with SOMs to classify the accumulated cyclone activity (ACA) associated with each pattern. We then evaluate the skill of CMIP6 historical experiments in simulating the LSMPs and ETC events identified in the SOM. Our results identify a robust bias toward more zonal patterns, with models struggling to reproduce the more amplified patterns typically associated with the highest cyclone activity. While model resolution has some impact on simulation credibility, model configuration appears to be more important in LSMP representation. The vast majority of CMIP6 models produce too few ETCs, although model errors are distributed around historical reanalyses when ACA is normalized by storm frequency.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

In situ multi-tier auto-ignition detection applied to dual-fuel combustion simulations

Here we use an anomaly detection methodology that is centered on analyzing fourth-order joint moments (co-kurtosis), particularly focusing on its application in auto-ignition of combustion problems with large numbers of species. Unsupervised anomaly detection is challenging to generalize across problem types and domains. A recent technique, centered on analyzing information in the fourth-order joint moment co-kurtosis, has shown promise, especially for high-dimensional scientific data. In this work we present developments to the co-kurtosis based anomaly detection method needed to make it effective and scalable for large-scale distributed scientific data, such as those generated by massively parallel simulations. An in situ co-kurtosis algorithm is employed as the anomaly detection method for identifying ignition kernels in simulations of turbulent combustion. Here, we extend an existing methodology which identifies regions of the domain where anomalies are present, and add another tier of anomaly detection where the individual samples contributing to the anomaly are identified. We apply this algorithm on-the-fly to a variety of turbulent reacting flow problems and compare it to the widely used (but significantly more expensive) chemical explosive mode analysis (CEMA). We demonstrate the ability of the method to detect and identify the onset of low and high temperature ignition which can be used for computational steering, as chemical and combustion anomalies occur intermittently at spatio-temporal locations unknown a priori. Finally, we apply our lightweight in situ algorithm to an exascale high-fidelity simulation with a total of 2.4 Trillion degrees of freedom, performed using an adaptive mesh refinement solver. Furthermore, through a scalability analysis, we show that the relative computational cost of this in-situ anomaly detection algorithm compared to an iteration of the reacting flow solver is negligible.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Dynamic-Phasor Model of Transformer Inrush Simulation for Unbalanced Distribution System

Abstract— Power system stability is a critical concern, with respect to the increasing distributed energy resources, microgrids, and emphasis on enhancing resiliency. For stability analysis and evaluation, modeling and simulation capabilities are important. In this paper, the authors’ previous work on modeling transformer inrush dynamics is improved with a new formulation of the key equation that converts the dynamic-phasor into real-time value. Simulation results in an electromechanical time-scale are compared with PSCAD results in an electromagnetic time-scale to demonstrate the improvement. The impact of simulation step size is discussed. A bug in PSCAD transformer setting that leads to incorrect magnetizing flux and current under Delta winding type is identified and reported. In addition, explanations and analysis are provided on two modeling techniques utilized since the previous work. This work is being incorporated into the next release of GridLAB-D. The improved transformer inrush calculation enables more accurate predictions of current amplitudes for analyzing the dynamics of distribution operation and control.

power transformers, power system simulation, satur↗
DOE OSTI2020

Structure and Dynamics of Aqueous Electrolytes Confined in 2D-TiO 2 /Ti 3 C 2 T 2 MXene Heterostructures

The synthesis of heterostructures of different two-dimensional (2D) materials offers an approach to combine advantages of different materials constituting the heterostructure and ultimately enhance their performance for applications such as electrochemical energy storage, achieving high energy, and high-power densities. Understanding the behavior of ions and solvents in confinement between these dissimilar layers is critical to understand their performance and control. Considering aqueous electrolytes, we explore the heterostructure of 2D lepidocrocite-type TiO 2 (2D-TiO 2 ) and hydroxylated or O-terminated Ti 3 C 2 MXene using ReaxFF molecular dynamics simulations and elastic/quasielastic neutron scattering techniques. Simulating a bilayer water intercalation, we find that the extent of interlayer hydration is impacted most by the surface terminations on the MXene and is marginally affected by 2D-TiO 2 . However, the introduction of 2D-TiO 2 decreases the water self-diffusion due to the notch sites (i.e., surface oxygen ridges) entrapping water molecules. Intercalating alkali cations into the heterostructures, we find that Li + is predominantly adsorbed at the 2D-TiO 2 surface instead of the MXenes with the preferential occupation of the notch sites. In contrast, Na+ forms a planar solvation with water, while K + is adsorbed both at the O-terminated MXene and 2D-TiO 2 . This behavior is altered when OH-terminated MXene is involved—the repulsion from the protons on the MXene surface forces the K + ions to be adsorbed exclusively to 2D-TiO 2 , while Na + retains some of its solvation in the water layer due to its smaller size. In OH-terminated MXenes, we see a consistent transfer of protons from the MXene surface toward 2D-TiO 2 , implying a greater capacity to store protons in the heterostructures. Of the three cations simulated, Na + hinders the proton migration the least and Li + the most because of its position near the 2D-TiO 2 surface. Furthermore, 2D-TiO 2 /MXene heterostructures are likely to exhibit a higher energy density but lower power density, especially with Na + intercalation.

25 ENERGY STORAGE↗
DOE OSTI2020

Highly-Active and Contaminant-Tolerant Cathodes for Durable Solid Oxide Fuel Cells (Final Report)

The objective of this project is to investigate the fundamental degradation methods occurring in solid oxide fuel cell (SOFC) cathodes when exposed to chromium and carbon dioxide contamination and to rationally design alkaline earth-based catalysts to increase stability and tolerance to contaminate poisoning. With a mechanistic and fundamental understanding of the degradation methods, advanced catalytic surface modifications can be applied to reduce degradation. The specific technical objectives are: (1) To identify/develop new catalysts (alkaline-earth based-) that are compatible chemically with the state-of-the-art cathode materials at high temperatures required for fabrication and with contaminates commonly encountered under operating conditions; (2) To improve the infiltration process for optimal control of the thickness, composition, and uniformity of the catalyst coatings; (3) To evaluate the electro-catalytic activity toward ORR of the chemically-stable materials when exposed to different types of contaminants using electrical conductivity relaxation measurements on bar samples and performance evaluation of catalyst-infiltrated cathodes; (4) To unravel the contamination-tolerant mechanisms of the new catalyst coatings under realistic environmental conditions (with different types of contaminants) using powerful in situ and in operando characterization techniques performed on model cells with thin-film/pattern electrodes, as guided by modeling and simulation; (5) To establish scientific basis for rational design of new catalysts of high tolerance to contaminants; (6) To validate the long term stability of modified LSCF cathodes in commercially available cells/stacks under ROC. Alkaline earth-based catalysts have been systematically explored and tested under various contamination conditions. BaCoO 3-δ (BCO) was shown to produce the best catalytic activity enhancement as well as stability in a variety of contaminating conditions, including CO 2 and chromium. The microstructural evolution was investigated with SEM, EDS, and Raman spectroscopy, showing the BCO catalyst prevents the formation of insulating SrCrO 4 by forming electrically conductive BaCrO 4 . Electrochemical relaxation measurements determined the BCO catalyst coatings increased the surface diffusion coefficient but did not significantly affect the diffusion coefficient and determined an optimum surface modification layer of 100 nm. A novel Ba 0.9 Co 0.7 Fe 0.2 Nb 0.1 O x (BCFN) catalyst was also shown to have excellent stability in chromium containing atmospheres. A novel surface sol-gel (SSG) surface modification, which offers superior thickness and compositional control, was applied to SOFCs to produce BaO catalyst coatings. SSG BaO coatings were shown to produce conformal coatings on the electrode surface, greatly increasing stability in chromium containing atmospheres. An atomic level mechanistic approach was applied to investigate the stability of PrBa 0.8 Ca 0.2 Co2O 5+δ (PBCC) with respect to a variety of common contaminants, demonstrating thermodynamically its superior stability in CO2 containing atmospheres. Finally, the best catalyst coatings demonstrated in this project, BCO and PBCC, were applied to full cells which demonstrated superior stability for over 300 hours in CO 2 and Cr atmospheres.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Highly-Active and Contaminant-Tolerant Cathodes for Durable Solid Oxide Fuel Cells

The objective of this project is to investigate the fundamental degradation methods occurring in solid oxide fuel cell (SOFC) cathodes when exposed to chromium and carbon dioxide contamination and to rationally design alkaline earth-based catalysts to increase stability and tolerance to contaminate poisoning. With a mechanistic and fundamental understanding of the degradation methods, advanced catalytic surface modifications can be applied to reduce degradation. The specific technical objectives are: (1) To identify/develop new catalysts (alkaline-earth based-) that are compatible chemically with the state-of-the-art cathode materials at high temperatures required for fabrication and with contaminates commonly encountered under operating conditions; (2) To improve the infiltration process for optimal control of the thickness, composition, and uniformity of the catalyst coatings; (3) To evaluate the electro-catalytic activity toward ORR of the chemically-stable materials when exposed to different types of contaminants using electrical conductivity relaxation measurements on bar samples and performance evaluation of catalyst-infiltrated cathodes; (4) To unravel the contamination-tolerant mechanisms of the new catalyst coatings under realistic environmental conditions (with different types of contaminants) using powerful in situ and in operando characterization techniques performed on model cells with thin-film/pattern electrodes, as guided by modeling and simulation; (5) To establish scientific basis for rational design of new catalysts of high tolerance to contaminants; (6) To validate the long term stability of modified LSCF cathodes in commercially available cells/stacks under ROC. Alkaline earth-based catalysts have been systematically explored and tested under various contamination conditions. BaCoO 3-δ (BCO) was shown to produce the best catalytic activity enhancement as well as stability in a variety of contaminating conditions, including CO 2 and chromium. The microstructural evolution was investigated with SEM, EDS, and Raman spectroscopy, showing the BCO catalyst prevents the formation of insulating SrCrO 4 by forming electrically conductive BaCrO 4 . Electrochemical relaxation measurements determined the BCO catalyst coatings increased the surface diffusion coefficient but did not significantly affect the diffusion coefficient and determined an optimum surface modification layer of 100 nm. A novel Ba 0.9 Co 0.7 Fe 0.2 Nb 0.1 O x (BCFN) catalyst was also shown to have excellent stability in chromium containing atmospheres. A novel surface sol-gel (SSG) surface modification, which offers superior thickness and compositional control, was applied to SOFCs to produce BaO catalyst coatings. SSG BaO coatings were shown to produce conformal coatings on the electrode surface, greatly4 increasing stability in chromium containing atmospheres. An atomic level mechanistic approach was applied to investigate the stability of PrBa 0.8 Ca 0.2 Co 2 O 5+δ (PBCC) with respect to a variety of common contaminants, demonstrating thermodynamically its superior stability in CO 2 containing atmospheres. Finally, the best catalyst coatings demonstrated in this project, BCO and PBCC, were applied to full cells which demonstrated superior stability for over 300 hours in CO 2 and Cr atmospheres.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2020

Highly-Active and Contaminant-Tolerant Cathodes for Durable Solid Oxide Fuel Cells

The objective of this project is to investigate the fundamental degradation methods occurring in solid oxide fuel cell (SOFC) cathodes when exposed to chromium and carbon dioxide contamination and to rationally design alkaline earth-based catalysts to increase stability and tolerance to contaminate poisoning. With a mechanistic and fundamental understanding of the degradation methods, advanced catalytic surface modifications can be applied to reduce degradation. The specific technical objectives are: (1) To identify/develop new catalysts (alkaline-earth based-) that are compatible chemically with the state-of-the-art cathode materials at high temperatures required for fabrication and with contaminates commonly encountered under operating conditions; (2) To improve the infiltration process for optimal control of the thickness, composition, and uniformity of the catalyst coatings; (3) To evaluate the electro-catalytic activity toward ORR of the chemically-stable materials when exposed to different types of contaminants using electrical conductivity relaxation measurements on bar samples and performance evaluation of catalyst-infiltrated cathodes; (4) To unravel the contamination-tolerant mechanisms of the new catalyst coatings under realistic environmental conditions (with different types of contaminants) using powerful in situ and in operando characterization techniques performed on model cells with thin-film/pattern electrodes, as guided by modeling and simulation; (5) To establish scientific basis for rational design of new catalysts of high tolerance to contaminants; (6) To validate the long term stability of modified LSCF cathodes in commercially available cells/stacks under ROC. Alkaline earth-based catalysts have been systematically explored and tested under various contamination conditions. BaCoO 3-δ (BCO) was shown to produce the best catalytic activity enhancement as well as stability in a variety of contaminating conditions, including CO 2 and chromium. The microstructural evolution was investigated with SEM, EDS, and Raman spectroscopy, showing the BCO catalyst prevents the formation of insulating SrCrO 4 by forming electrically conductive BaCrO 4 . Electrochemical relaxation measurements determined the BCO catalyst coatings increased the surface diffusion coefficient but did not significantly affect the diffusion coefficient and determined an optimum surface modification layer of 100 nm. A novel Ba 0.9 Co 0.7 Fe 0.2 Nb 0.1 O x (BCFN) catalyst was also shown to have excellent stability in chromium containing atmospheres. A novel surface sol-gel (SSG) surface modification, which offers superior thickness and compositional control, was applied to SOFCs to produce BaO catalyst coatings. SSG BaO coatings were shown to produce conformal coatings on the electrode surface, greatly 4 increasing stability in chromium containing atmospheres. An atomic level mechanistic approach was applied to investigate the stability of PrBa 0.8 Ca 0.2 Co 2 O 5+δ (PBCC) with respect to a variety of common contaminants, demonstrating thermodynamically its superior stability in CO 2 containing atmospheres. Finally, the best catalyst coatings demonstrated in this project, BCO and PBCC, were applied to full cells which demonstrated superior stability for over 300 hours in CO 2 and Cr atmospheres

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2020

Thermal Frequency Drift of 3D Printed Microwave Components

Fabrication of microwave slot array antennas and waveguide bandpass and notch filters using 3D printing has significant advantages in terms of speed and cost even for parts with high mechanical complexity. One disadvantage of Stereolithography (SLA) 3D printed, copper plated microwave components is that some SLA resins have a high Coefficient of Thermal Expansion (CTE), quoted in micrometers per meter per degree or 10 -6 per degree. Compared to typically used metals such as aluminum (CTE 24 × 10 -6 ·K -1 ) and copper (CTE 17 × 10 -6 ·K -1 ), SLA resin can have CTE above 100 × 10 -6 ·K -1 . Resonant structures experience significant frequency drift with temperature changes on the order of 10–50 °C. The issue of 3D printed microwave structures changing frequency characteristics significantly with temperature shift has not been addressed or reviewed in current literature. We measured and simulated the effect of temperature change on a slot array, cavity notch filters, and post loaded waveguide bandpass filters. We tested several types of SLA resin, different plating techniques, and also Direct Metal Laser Sintering (DMLS) and Binder Infusion metal 3D printing. Performance as a function of temperature is presented for these alternatives.

three-dimensional printing↗
DOE OSTI2023

Physics guided machine learning for multi-material decomposition of tissues from dual-energy CT scans of simulated breast models with calcifications

We introduce a physics guided data-driven method for image-based multi-material decomposition for dual-energy computed tomography (CT) scans. The method is demonstrated for CT scans of virtual human phantoms containing more than two types of tissues. The method is a physics-driven supervised learning technique. We take advantage of the mass attenuation coefficient of dense materials compared to that of muscle tissues to perform a preliminary extraction of the dense material from the images using unsupervised methods. We then perform supervised deep learning on the images processed by the extracted dense material to obtain the final multi-material tissue map. The method is demonstrated on simulated breast models with calcifications as the dense material placed amongst the muscle tissues. The physics-guided machine learning method accurately decomposes the various tissues from input images, achieving a normalized root-mean-squared error of 2.75%.

Gopalakrishnan Meena, Murali↗
DOE OSTI2022

A Fast Time-Stepping Strategy for Dynamical Systems Equipped with a Surrogate Model

Simulation of complex dynamical systems arising in many applications is computationally challenging due to their size and complexity. Model order reduction, machine learning, and other types of surrogate modeling techniques offer cheaper and simpler ways to describe the dynamics of these systems but are inexact and introduce additional approximation errors. In order to overcome the computational difficulties of the full complex models, on one hand, and the limitations of surrogate models, on the other, this work proposes a new accelerated time-stepping strategy that combines information from both. This approach is based on the multirate infinitesimal general-structure additive Runge--Kutta framework. The inexpensive surrogate model is integrated with a small time step to guide the solution trajectory, and the full model is treated with a large time step to occasionally correct for the surrogate model error and ensure convergence. Here, we provide a theoretical error analysis, and several numerical experiments, to show that this approach can be significantly more efficient than using only the full or only the surrogate model for the integration.

Surrogate models↗
DOE OSTI2025

Elastic Bayesian Model Calibration

Functional data are ubiquitous in scientific modeling. For instance, quantities of interest are modeled as functions of time, space, energy, density, etc. Uncertainty quantification methods for computer models with functional response have resulted in tools for emulation, sensitivity analysis, and calibration that are widely used. However, many of these tools do not perform well when the computer model’s parameters control both the amplitude variation of the functional output and its alignment (or phase variation). This paper introduces a framework for Bayesian model calibration when the model responses are misaligned functional data. The approach generates two types of data out of the misaligned functional responses: (1) aligned functions so that the amplitude variation is isolated and (2) warping functions that isolate the phase variation. These two types of data are created for the computer simulation data (both of which may be emulated) and the experimental data. The calibration approach uses both types so that it seeks to match both the amplitude and phase of the experimental data. The framework is careful to respect constraints that arise, especially when modeling phase variation, and is framed in a way that it can be done with readily available calibration software. In conclusion, we demonstrate the techniques on two simulated data examples and on two dynamic material science problems: a strength model calibration using flyer plate experiments and an equation of state model calibration using experiments performed on the Sandia National Laboratories’ Z-machine.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Systems and methods for direct estimation of battery parameters using only charge/discharge curves

Electrochemical models for the lithium-ion battery are useful in predicting and controlling its performance. The values of the parameters in these models are vital to their accuracy. However, not all parameters can be measured precisely, especially when destructive methods are prohibited. In some embodiments of the present disclosure, a parameter estimation approach is used to estimate the open circuit potential of the positive electrode (Up) using piecewise linear approximation together with all the other parameters of a single particle model. Up and 10 more parameters may be estimated from a single discharge curve without knowledge of the electrode chemistry using a technique such as a genetic algorithm. Different case studies were presented for estimating Up with different types of parameters of the battery model. The estimated parameters were then validated by comparing simulations at different discharge rates with experimental data.

Qi, Yanbo↗
DOE OSTI2025

On the viability of stimulated hydrogen generation from iron-rich formations

Hydrogen-based technologies present a promising solution for the global energy transition. In addition to electrolytic production, subsurface geological formations provide a potential natural source of hydrogen. Iron-rich ultramafic rocks, in particular, are favorable for hydrogen generation through natural processes such as serpentinization. Naturally occurring reactions and migration can be enhanced through various types of stimulation, including thermal, hydraulic, and chemical treatment. Through numerical simulations, we analyzed the complex interplay of factors influencing the production and migration within the subsurface, emphasizing the importance of different stimulation techniques, catalysts, and conditions. Our findings indicate that key parameters, such as damage zone permeability and width, significantly impact producible hydrogen mass. Our results indicate that a combination of large damage zone widths, high permeability, and a stimulated reaction rate of 1 × 10 -9 can yield economically viable production rates of up to 1 kg s -1 at the wellhead. Moreover, the availability of ferrous iron, rather than the serpentinization rate itself, has been identified as the primary limiting factor in achieving economically sustainable hydrogen production. In conclusion, while an unstimulated rock volume of 0.165 km 3 yields only 45t of hydrogen in two years, various stimulation techniques can increase production to 18500t.

08 - HYDROGEN↗