SEARCH · Search NASA

Results for “Advanced Modeling and 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.

SourceResearch connectionFrom yearThrough year

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

At least 1,477 records · Page 82

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Computational Methods in Heterogeneous Catalysis

The unprecedented ability of computations to probe atomic-level details of catalytic systems holds immense promise for the fundamentals-based bottom-up design of novel heterogeneous catalysts, which are at the heart of the chemical and energy sectors of industry. Here, we critically analyze recent advances in computational heterogeneous catalysis. First, we will survey the progress in electronic structure methods and atomistic catalyst models employed, which have enabled the catalysis community to build increasingly intricate, realistic, and accurate models of the active sites of supported transition-metal catalysts. We then review developments in microkinetic modeling, specifically mean-field microkinetic models and kinetic Monte Carlo simulations, which bridge the gap between nanoscale computational insights and macroscale experimental kinetics data with increasing fidelity. Here, we finally review the advancements in theoretical methods for accelerating catalyst design and discovery. Throughout the review, we provide ample examples of applications, discuss remaining challenges, and provide our outlook for the near future.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Development of a high charge 10 GeV laser electron accelerator

Low-density meter-scale plasma waveguides produced in meter-scale supersonic gas jets have paved the way for recent demonstrations of all-optical multi-gigaelectronvolt laser wakefield acceleration (LWFA). This paper reviews recent advances by the University of Maryland, which have enabled these results, focusing on the development of elongated supersonic gas jets up to ∼1 m in length, experimental and simulation studies of plasma waveguide formation, and a new three-stage model for relativistic pulse propagation dynamics in these waveguides. We also present results from recent LWFA experiments conducted at the Laboratory for Advanced Lasers and Extreme Photonics at Colorado State University demonstrating high charge, low divergence electron bunches to ∼10 GeV, with laser-to-electron beam efficiency of at least ∼30%.

Rockafellow, E. [Univ. of Maryland, College Park, ↗

DOE OSTI2019

Efficient Simulation and Abuse Modeling of Mechanical-Electrochemical-Thermal Phenomena in Li-Ion Batteries

This poster covered simultaneous coupling of electrochemical-thermal models with mechanical deformation in lithium ion batteries. Efficiency and stability of mechanical models was significantly enhanced by implementing electrochemical models into LS-DYNA using User-Defined Elements. Six case studies were built and licensed out to participants from Industry for initial testing and their feedback is being incorporated into these tools. Dynamic response of the cells was incorporated by measuring mechanical response of components at strain rates as high as 250 /s. Temperature range for property measurements was expanded (as high as 200 degrees C) to account for property changes at high temperatures experienced by cell components under battery abuse. Multi-cell validation has been expanded to include four different sets of experimental data, with support from various partners. Complex failure modes and fracture response are currently being investigated. These are still very challenging, given the limited amount of prior work available in the literature.

25 ENERGY STORAGE↗

DOE OSTI2019

Using Computational Fluid Dynamics of Wind Simulations Coupled With Weather Data to Calculate Dynamic Line Ratings

Dynamic line rating is a technology that allows for the rating of electrical conductors to be calculated based on local weather conditions rather than using ”worst-case” assumptions of weather conditions. The static ratings are typically based on weather conditions that are conservative in nature, and thus dynamic line rating may provide room to increase the current capacity. The local weather conditions may estimate conditions more accurately near the transmission lines accounting for terrain such as nearby hills. One way to estimate the local weather conditions is through the use of computational fluid dynamics. The results of a wind field simulation can be coupled with sparsely located weather stations to provide an accurate assessment of the wind speeds along the path of the conductor. Here, the wind field simulations are calculated using a steady state Reynolds-averaged Navier-Stokes model. These results along with local solar radiance measurements can be used to provide estimates of the steady state ampacity with the standard IEEE equations. An advanced test case for this process involves performing these simulations in a complex terrain, Hells Canyon, the deepest river gorge in North America. The results show that by using weather-based sensors without considering localized wind conditions, the available ampacity may be over-predicted significantly.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2024

2019 Budget Request for the DOE Computational Science Graduate Fellowship (CSGF) Grant

The Department of Energy Computational Science Graduate Fellowship (DOE CSGF) is necessary to meet the continual challenging national workforce needs that arise as computational science and engineering problems continue to grow in scope and complexity. Computational science and engineering (CSE) is a multidisciplinary approach that uses scientific computing to solve practical problems methods and to supply technical tools across the scientific discovery spectrum. In particular, the DOE CSGF emphasizes high-performance computing (HPC) that enables CSE that advances science and engineering in directions important to the DOE and the economy in general. Over the past half-century, HPC has been an essential tool for DOE’s success. During this period, important missions, such as nuclear stockpile stewardship, have turned to HPC as an essential technology. Entire science disciplines, such as biology and cosmology, have been transformed through the augmentation of scientific observation via HPC. At government laboratories and in industry, DOE CSGF alumni are helping push traditional HPC boundaries while contributing to discoveries in high-energy physics, renewable energy, fusion-reactor design, additive manufacturing, nanomaterials for next-generation batteries and transistors, and turbine and advanced nuclear reactor modeling. In addition, HPC is used to address national health needs that will eventually point to cures both by helping cancer researchers manage and analyze huge troves of data, by simulating biological mechanisms, and by accelerating drug development — including continuing to rise to the challenge of pandemic-related research. A 2023 report from the ASCAC Subcommittee on American Competitiveness and Innovation to the ASCR office, “Can the United States Maintain Its Leadership in High-Performance Computing?” says of the Program, “The CSGF program provides a barometer for disciplines that will be of interest to future DOE computing.” An explosion in scientific and technological data has driven the need for increasingly sophisticated HPC to transform those data into scientific understanding. With access to more and more data and the proliferation of HPC, Machine Learning and Artificial Intelligence are experiencing a renaissance, complementing the now well-established use of computational simulation. Indeed, in its September 2020 subcommittee report on “AI/ML, Data Intensive Science and High-Performance Computing”, the DOE Advanced Scientific Computing Advisory Committee (ASCAC) explicitly called for a fellowship program to train computational and data scientists to tackle exascale and data-intensive computing challenges. This collaboration of empirical and theory-based modeling will increasingly inform federal policymakers whose decisions affect American society and future generations, and it requires highly skilled and intellectually agile computational scientists who can support the fast-moving DOE National Laboratory research environment. In fact, the DOE CSGF program has explicitly and consistently addressed this need.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Using Homogenized Models to Explore the Effect of Fracture Densities on Weathering

Despite its importance, only a few researchers have incorporated the effects of fracturing into models of reactive transport for rock weathering. Here we explore 2D simulations that describe weathering under conditions of diffusive and advective transport within heterogeneous media consisting of rocky blocks and fractures. In our simulations, the Darcy velocities vary in space and time and depend on weathering processes within the rock matrix. We explore simulations with saturated and unsaturated flow for weathering bedrock that consists of blocks separated by inert or weathered material. The simulations show that a simplified homogenized model can approximate exact solutions for some of the simulated columns and hills and can allow exploration of coupling between flow and reaction in fractured rock. These hillslope simulations document that, even in the presence of 2D water flow, i) an increase in fracture density results in faster weathering advance rates; and ii) the water table locates deeper for a rock system that is weathered and fractured rather than weathered and unfractured. Some of these patterns have also been observed for natural systems. But these simulations also highlight how simplified models that do not use appropriate averaging of heterogeneities can be inaccurate in predicting weathering rate for natural systems. For example, if water flows both vertically and laterally through the vadose zone of a hill, then a prediction of the depth of regolith that is based on modeling strictly unidirectional downward infiltration will be unrealistically large. Likewise, if the fracture density observed near the land surface is used in a model to predict depth of weathering for a system where the fracture density decreases downward, the model will overestimate regolith depth. Learning how to develop accurately homogenized models could thus enable better conceptual models and predictions of weathering advance in natural systems.

58 GEOSCIENCES↗

DOE OSTI2025

AI‐Driven Defect Engineering for Advanced Thermoelectric Materials

Thermoelectric materials offer a promising pathway to directly convert waste heat to electricity. However, achieving high performance remains challenging due to intrinsic trade-offs between electrical conductivity, the Seebeck coefficient, and thermal conductivity, which are further complicated by the presence of defects. This review explores how artificial intelligence (AI) and machine learning (ML) are transforming thermoelectric materials design. Advanced ML approaches including deep neural networks, graph-based models, and transformer architectures, integrated with high-throughput simulations and growing databases, effectively capture structure-property relationships in a complex multiscale defect space and overcome the “curse of dimensionality”. This review discusses AI-enhanced defect engineering strategies such as composition optimization, entropy and dislocation engineering, and grain boundary design, along with emerging inverse design techniques for generating materials with targeted properties. Finally, it outlines future opportunities in novel physics mechanisms and sustainability, highlighting the critical role of AI in accelerating the discovery of thermoelectric materials.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Optimization of direct air capture processes using reactive transport models of adsorption-desorption cycles

In this study, we develop and implement a reactive transport model in COMSOL Multiphysics® to address the challenges of direct air carbon capture. The model is validated against experimental data and used to simulate the cyclic steady state of the adsorption-desorption process. The optimization of this model is achieved through advanced trust-region methods integrated with Gaussian Processes. Key decision variables, including adsorption and desorption times, desorption temperature and pressure, input velocity, bed porosity, column length, and radius were optimized to minimize the capture cost. After optimization, a sensitivity analysis revealed the complex interplay between the decision variables and their effect on the specific energy and cost of removing the CO 2 . We optimized the capture cost while taking into account the trade-off between energy consumption and productivity. The resulting minimum capture cost was determined to be 265.2 $/t-CO 2 , which aligns with expected values reported in the literature. Numerical results suggest the effectiveness of the optimization strategies applied, and underscore the importance of simultaneous decision variable selection in improving the performance in direct air capture processes. We also extend the modeling approach to a 2D axisymmetric model to better visualize CO₂ uptake and temperature profiles, revealing significant radial gradients during the regeneration step. As a main drawback, this enhanced model comes with a computational cost approximately 40 times higher than that of the 1D model.

Adsorption-desorption process↗

DOE OSTI2022

Social vulnerability and power loss mitigation: A case study of Puerto Rico

The increasing occurrence of extreme weather events urges us to reevaluate the resiliency and vulnerability aspects of our most critical infrastructures — such as power grids — as their failures result in both economic loss and severe human hardship. Seen through the lens of alleviating human suffering, it is crucial to be able to identify critical system components of the infrastructure for targeted hardening given resource constraints. This effort is of particular importance in islanded areas such as Puerto Rico where hurricanes are frequent and resources are limited, and where the spatially diverse effects of power loss on human suffering are all the more severe. Recent studies on evaluating infrastructure networks during extreme weather events have taken a simulation based approach that incorporates a variety of component models, such as weather realizations, topological network models, fragility models, and power flow models to estimate expected loss of service. Here, in this work, we expand such a Component Based Event Simulation (CBES) methodology proposed in the literature and integrate it with a social vulnerability modeling component. This paradigm-advancing approach of synthesizing the cutting edge capability of power network modeling and the social impacts of the power transmission network failure is demonstrated for the island of Puerto Rico. Our work exemplifies the efficacy of this integrated modeling framework in developing a decision metric for targeted transmission line hardening.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2019

Three-dimensional flow model development for thermal mixing and stratification modeling in reactor system transients analyses

Mixing, thermal-stratification, and mass transport phenomena in large pools or enclosures play major roles for the safety of reactor systems. Depending on the fidelity requirement and computational resources, various modeling methods, from the 0-D perfect mixing model to 3-D Computational Fluid Dynamics (CFD) models, are available. Each is associated with its own advantages and shortcomings. It is very desirable to develop an advanced and efficient thermal mixing and stratification modeling capability embedded in a modern system analysis code to improve the accuracy of reactor safety analyses and to reduce modeling uncertainties. An advanced system analysis tool, SAM, is being developed at Argonne National Laboratory for advanced non-LWR reactor safety analysis. While SAM is being developed as a system-level modeling and simulation tool, a reduced-order three-dimensional module is under development to model the mull-dimensional flow and thermal mixing and stratification in large enclosures of reactor systems. This paper provides an overview of the three-dimensional finite element flow model in SAM, including the governing equations, stabilization scheme, and solution methods. Additionally, several verification and validation tests are presented, including lid-driven cavity flow, natural convection inside a cavity, laminar flow in a channel of parallel plates. Here, based on the comparisons with the analytical solutions and experimental results, it is demonstrated that the developed 3-D fluid model can perform very well for a wide range of flow problems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Simulating Pyrocumulonimbus Clouds Using a Multiscale Wildfire Simulation Framework

Pyrocumulonimbus (pyroCb) clouds, driven by extreme fires under favorable meteorological conditions, can inject smoke into the stratosphere at magnitudes comparable to those of moderate volcanic eruptions, potentially altering the global radiative balance and atmospheric composition. However, simulating pyroCb is particularly challenging in Earth system models. Using the Energy Exascale Earth System Model (E3SM), we developed a novel global multiscale framework to model pyroCb events in California, which includes a high‐resolution fire radiative power time series, a one‐dimensional plume‐rise parameterization, a fire‐induced vertical water vapor transport scheme, and a surface wildfire sensible heat flux representation. Our simulation successfully reproduces many pyroCb features, including cloud height, spatiotemporal evolution, and convective intensity in comparison with satellite and ground‐based observations. Sensitivity experiments show that realistic pyroCb simulation depends on vertical water vapor transport. These advances provide a basis for future exploration of pyroCb impacts at regional and global scales within climate models.

E3SM↗

DOE OSTI2026

Finite-element-based simulations of electrodes for CO 2 cascade reduction reactions

The multielectron reduction of CO 2 to liquid fuels could be a path to scalable energy storage, but reaching this goal requires major advances in catalysis and systems engineering. Cascade catalysis, which couples sequential reactions without isolating intermediates, has emerged as a promising route to enhance selectivity and efficiency in CO 2 reduction (CO 2 R). In this review, we examine how finite-element-based simulations of continuum model [finite element method (FEM)] approaches are being used to analyze and guide CO 2 R cascade systems. We first outline the fundamentals of cascade catalysis and recent advances in catalytic materials (metallic, molecular, and hybrid architectures). We then focus on FEM developments at the electrode and device scales, emphasizing how these models capture transport phenomena, local microenvironments, and geometry-dependent effects. To clarify design principles, we present case studies of cascade electrodes organized in systems without and with integrated semiconductors. We further emphasize the integration of FEM with multiscale frameworks (density functional theory, molecular dynamics, kinetic Monte Carlo) and its role in bridging atomic-level insights with device-level performance. Finally, we identify current limitations and future prospects, including improved boundary conditions, coupling with operando experiments, and machine learning-accelerated model development. Together, these insights provide design principles for next-generation CO 2 R cascade systems for efficient solar fuel production.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Mechanical cosmology: Simulating scalar fluctuations in expanding universes using synthetic mechanical lattices

Inspired by recent advances in observational astrophysics and continued explorations in the field of analog gravity, we discuss the prospect of simulating models of cosmology within the context of synthetic mechanical lattice experiments. We focus on the physics of expanding universe scenarios described by the Friedmann-Lemaître-Robertson-Walker (FLRW) metric. Specifically, quantizing scalar fluctuations in a background FLRW spacetime leads to a quadratic bosonic Hamiltonian with temporally varying pair production terms. Here we present a mapping that provides a one-to-one correspondence between these classes of cosmology models and feedback-coupled mechanical oscillators. As proof of principle, we then perform experiments on a synthetic mechanical lattice composed of such oscillators. We simulate two different FLRW expansion scenarios with universes dominated by vacuum energy and matter and discuss our experimental results. Published by the American Physical Society 2025

Rhyno, Brendan (ORCID:0000000235266998)↗

DOE OSTI2025

Real-Time Lifetime Prediction of Semiconductor Devices Using Hardware-in-the-Loop

This paper presents a unique approach to enable real-time lifespan prediction of semiconductor power modules using a Hardware-in-the-Loop (HIL) system. By integrating the module's overall loss characteristics-specifically switching and conduction losses-with a thermoelectric model of the thermal management system, this research demonstrates that the model can dynamically estimates the junction temperature profile of the semiconductor devices in response to a changing torque demand profile for the motor drive system. This capability enables continuous monitoring of the module's operational time and cumulative stress induced on the devices to compute accumulated remaining lifetime or time-to-failure (TTF). This study provides an architectural framework for the HIL system with high-fidelity component models of multiple physical domains, allowing simulation of dynamic behaviors of a closely-coupled motor drive system. The advanced real-time computation and measurement functionalities of the HIL system allow for both dynamic lifetime calculations based on simulated data and aggregate lifetime predictions utilizing historical data. Moreover, this paper details an algorithm that not only computes cumulative damage but also synthesizes these data into a comprehensive aggregated lifetime metric. This methodology can enhance the maintenance scheduling strategies and operational reliability of semiconductor devices in critical applications, ultimately extending their service life while optimizing performance.

hardware-in-the-loop (HIL)↗

DOE OSTI2018

Recent Advancements in the Numerical Simulation of Surface Irradiance for Solar Energy Applications

This paper briefly reviews the National Renewable Energy Laboratory's recent efforts to develop all-sky solar irradiance models for solar energy applications. The Fast All-sky Radiation Model for Solar applications (FARMS) uses the simulation of clear-sky transmittance and reflectance and a parameterization of cloud transmittance and reflectance to rapidly compute broadband irradiances on horizontal surfaces. The accuracy of FARMS is comparable to that of two-stream approximation, but it is approximately 1,000 times faster. A FARMS for Narrowband Irradiance over Tilted surfaces (FARMS-NIT) has been developed to compute spectral irradiances on photovoltaic (PV) panels in 2,002 wavelength bands. FARMS-NIT has been extended to bifacial PV panels by accounting for solar radiation reaching the backside of a PV panel.

14 SOLAR ENERGY↗

DOE OSTI2018

Interaction Between Proppant Compaction and Single-/Multiphase Flows in a Hydraulic Fracture

Understanding proppant-transport and -deposition patterns in a hydraulic fracture is vital for effective and economical production of petroleum hydrocarbons. In this research, a numerical-modeling approach, combining the discrete-element method (DEM) with single-/multiphase lattice Boltzmann (LB) simulation, was adopted to advance understanding of the interaction between reservoir depletion, proppant-particle compaction, and single-/multiphase flows in a hydraulic fracture. DEM was used to simulate effective-stress increase and the resultant proppant-particle movement and rearrangement during the process of reservoir depletion caused by hydrocarbon production. Simulated pore structure of the proppant pack was extracted and used as boundary conditions in the LB simulation to calculate the time-dependent permeability of the proppant pack. We first validated the DEM-LB numerical work flow, and the simulated proppant-pack permeabilities as functions of effective stress were in good agreement with laboratory measurements. Furthermore, three proppant packs with the same average particle diameter but different diameter distributions were generated to study the role of proppant-size heterogeneity (variation in particle diameter). Specifically, we used the coefficient of variation (COV) of diameter, defined as the ratio of standard deviation of diameter to mean diameter, to characterize the heterogeneity of particle size. We obtained proppant-pack porosity, permeability, and fracture-width reduction (closure distance) as functions of effective stress. Under the same effective stress, a proppant pack with a higher diameter COV had lower porosity and permeability and larger fracture-width reduction. This was because the high diameter COV gave rise to a wider diameter distribution of proppant particles; smaller particles were compressed into the pore space between larger particles with increasing stress, leading to larger closure distance and lower porosity and permeability. With multiphase LB simulation, relative permeability curves were obtained, which are critical for larger-scale reservoir simulations under various effective stresses. The relative permeability of the oil phase was more sensitive to changes in geometry and stress, compared with the wetting phase. This was because the oil phase occupied larger pores; compaction of the proppant pack affected the structure of the pores, because the pores were farther from the grain contacts and thus more susceptible to collapse. It is also interesting to note that when effective stress increased continuously, oil relative permeability increased first and then decreased. This nonmonotonic behavior was the result of the nonmonotonic development of pore structure and oil connectivity under increasing stress.

Fan, Ming↗

DOE OSTI2017

(U) Ristra Next Generation Code Report

LANL’s Weapons Physics management (ADX) and ASC program office have defined a strategy for exascale-class application codes that follows two supportive, and mutually risk-mitigating paths: evolution for established codes (with a strong pedigree within the user community) based upon existing programming paradigms (MPI+X); and Ristra (formerly known as NGC), a high-risk/high-reward push for a next-generation multi-physics, multi-scale simulation toolkit based on emerging advanced programming systems (with an initial focus on data-flow task-based models exemplified by Legion [5]). Development along these paths is supported by the ATDM, IC, and CSSE elements of the ASC program, with the resulting codes forming a common ecosystem, and with algorithm and code exchange between them anticipated. Furthermore, solution of some of the more challenging problems of the future will require a federation of codes working together, using established-pedigree codes in partnership with new capabilities as they come on line. The role of Ristra as the high-risk/high-reward path for LANL’s codes is fully consistent with its role in the Advanced Technology Development and Mitigation (ATDM) sub-program of ASC (see Appendix C), in particular its emphasis on evolving ASC capabilities through novel programming models and data management technologies.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2019

Detailed Modeling of Physical Processes in Electron Sources for Accelerator Applications

At present, electron sources are essential in a wide range of applications – from common technical use to exploring the nature of matter. Depending on the application requirements, different methods and materials are used to generate electrons. State of-the-art accelerator applications set a number of often-contradicting requirements for electron sources (e.g., high quantum efficiency vs. high spin-polarization, high current density vs. long lifetime, etc). Development of advanced electron sources includes modeling and design of cathodes, material growth, fabrication of cathodes, and cathode testing. The detailed simulation and modeling of physical processes are required in order to shed light on the exact mechanisms of electron emission and to develop new-generation electron sources with optimized efficiency. The purpose of the present work is to study physical processes in prominent electron sources and to develop scientific tools that could be used to predict electron emission from novel nano-structured materials. In particular, the Monte Carlo model has been implemented and used to study photoemission from Negative Electron Affinity (NEA) GaAs widely used to provide spin-polarized electron beams for many accelerator facilities, including Thomas Jefferson National Accelerator Facility. According to our best knowledge, this is the first to date scientific tool to simulate the complete from-excitation-to-emission-into-vacuum photoemission of spin-polarized electrons from NEA GaAs based on the details of material band structure. Moreover, an effective image-processing algorithm has been developed and used to study field emission properties of Nitrogen-incorporated UltraNanoCrystalline Diamond ((N)UNCD) films. The algorithm processes large sets of image micrographs representing field emission sites formed on the anode plate and therefore provides significant information about the field emission area and the current density. The field emission model from convenient semiconductors has been modified and used to study the influence of the space-charge region and the electric-field-dependent electron mobility on the field emission characteristics of (N)UNCD field emitters.

Chubenko, Oksana↗