Search NASA⌕ Search

SEARCH · Search NASA

Results for “computational frameworks”

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 361 records · Page 20

Scalable Control Co-design for Resilient-by-Design Cyber Physical Systems

Critical infrastructure networks, such as power and transportation networks, are often modelled as cyber-physical systems. With ever increasing complexity of these systems, there is a need for newer and more relevant metrics and design tools that will co-optimize the physical system components and control policies to guarantee resilience against cyber and natural threats. To this end, a simulation-based control co-design computational framework that will concurrently determine the system and control parameters of a cyber-physical system to meet pre-specified resilience, operational and economic objectives has been developed. The capabilities of the developed co-design engine are demonstrated by designing the physical components and control parameters of a microgrid system that will meet its resiliency objectives when subjected to various cyber and physical threats.

42 ENGINEERING↗

Land Model Testbed: Accelerating Development, Benchmarking and Analysis of Land Surface Models

A Land Model Testbed (LMT), designed to provide a computational framework for systematically assessing model fidelity and supporting rapid development of complex multiscale models, offers a general-purpose workflow for conducting large ensemble simulations of multiple land surface models, post-processing large volumes of model output, and evaluating model results. It leverages existing tools for launching model simulations and the International Land Model Benchmarking (ILAMB) package for assessing model fidelity through comparison with best-available observational datasets. Increased complexity and proliferation of uncertain parameters in process representations in land surface models has driven the need for frequent and intensive testing and evaluating of models to quantify uncertainties and optimize parameters such that results are consistent with observations. The LMT described here meets these needs by providing tools to run thousands of ensemble simulations simultaneously and post-process their output files, by automating execution of an enhanced version of ILAMB with site-specific benchmarks and multivariate functional relationships, and by offering ensemble diagnostics and a customizable dashboard for displaying model performance metrics and associated graphics. We envision the LMT capabilities will serve as a foundational computational resource for a proposed user facility focused on terrestrial multiscale model--data integration.

Sreepathi, Sarat↗

Vibrational Signatures of Electronic Properties in Renewable-Energy Catalysis

The objective of this research program is to discover and develop new approaches for ab initio computational simulations of molecular motion. Specifically, this program aims to decipher the connection between the underlying electronic structure and the accordant molecular vibrations in reactive ions and radicals for renewable-energy purposes. Recent experimental progress in ion sources and optical spectroscopies has unearthed considerable new inner-sphere detail for these complexes, but the connection between these spectral signatures and mechanistic information often remains elusive, to the continued frustration of experimentalists. For this purpose, new anharmonic vibrational frequency methods, along with a publicly deployed software package, will be developed. Working closely with committed experimental collaborators, this conceptual and computational framework will be used to explain the results of new spectroscopy experiments, focusing specifically on the inner-shell mechanisms of renewable-energy catalysis. The oxidation half of catalytic water-splitting chemistry will be a central focus, along with fundamental studies of the manner in which strong ions and radicals activate solvent as a chemical species. The resulting products of the research program will include openly available software and algorithms for the ab initio simulation of challenging vibrational spectra, as well as critical mechanistic insight into energy-focused catalytic processes that are opaque to other existing analytical techniques.

42 ENGINEERING↗

Mechanochemical spinodal decomposition: a phenomenological theory of phase transformations in multi-component, crystalline solids

Here, we present a phenomenological treatment of diffusion-driven martensitic phase transformations in multi-component crystalline solids that arise from non-convex free energies in mechanical and chemical variables. The treatment describes diffusional phase transformations that are accompanied by symmetry-breaking structural changes of the crystal unit cell and reveals the importance of a mechanochemical spinodal, defined as the region in strain-composition space, where the free-energy density function is non-convex. The approach is relevant to phase transformations wherein the structural order parameters can be expressed as linear combinations of strains relative to a high-symmetry reference crystal. The governing equations describing mechanochemical spinodal decomposition are variationally derived from a free-energy density function that accounts for interfacial energy via gradients of the rapidly varying strain and composition fields. A robust computational framework for treating the coupled, higher-order diffusion and nonlinear strain gradient elasticity problems is presented. Because the local strains in an inhomogeneous, transforming microstructure can be finite, the elasticity problem must account for geometric nonlinearity. An evaluation of available experimental phase diagrams and first-principles free energies suggests that mechanochemical spinodal decomposition should occur in metal hydrides such as ZrH 2-2c . The rich physics that ensues is explored in several numerical examples in two and three dimensions, and the relevance of the mechanism is discussed in the context of important electrode materials for Li-ion batteries and high-temperature ceramics.

36 MATERIALS SCIENCE↗

A cross-platform execution engine for the quantum intermediate representation

Hybrid languages like the quantum intermediate representation (QIR) are essential for programming systems that mix quantum and conventional computing models, while execution of these programs is often deferred to a system-specific implementation. Here, we develop the QIR Execution Engine (QIR-EE) for parsing, interpreting, and executing QIR across multiple hardware platforms. QIR-EE uses LLVM to execute hybrid instructions specifying quantum programs and, by design, presents extension points that support customized runtime and hardware environments. We demonstrate an implementation that uses the XACC quantum hardware-accelerator library to dispatch prototypical quantum programs on different commercial quantum platforms and numerical simulators, and we validate execution of QIR-EE on IonQ, Quantinuum, and IBM hardware. Our results highlight the efficiency of hybrid executable architectures for handling mixed instructions, managing mixed data, and integrating with quantum computing frameworks to realize cross-platform execution.

LLVM↗

Impact of Feed Injection Strategies on Fluidization Dynamics for Biomass Thermochemical Conversion

To better understand some of the key parameters that control biomass conversion processes in dense granular beds, an efficient computational framework for large-scale simulations of dense, reactive particulate flows using a Lagrange-Euler approach has been developed. This framework is applied here to the investigation of feed injection in a hot fluidized bed reactor, and how it may impact the biomass conversion dynamics. A simple, pseudo-two dimensional configuration is adopted to facilitate the parametric study. Chemical processes are modeled using global kinetics that accurately reproduce particle mass loss and gas release. A posteriori analysis of particle heating rate, mixing and segregation, along with products distribution and residence time inside the reactor is performed for different injection strategies, and compared to non-reactive cases. Results highlight some non-trivial coupling between chemistry and flow dynamics.

09 BIOMASS FUELS↗

Building a Prospective LCA Framework to Analyze Emerging Technologies in a Dynamic System Context

NREL's Lifecycle Analysis Integration into Opensource Numerical models (LiAISON) framework computes temporally explicit life cycle impacts and resource uses for specific technologies (foreground) in a dynamic system context (background). LiAISON computes results for a series of environmental mid-points enabling an analysis of prospective tradeoffs of emerging technologies toward 2100. This prospective feature is of critical importance when analyzing present-day emerging technologies whose large-scale impacts during deployment phases will occur in different, future system contexts. LiAISON systematically accounts for dynamic system changes by applying an integrated background of the future energy-economy-land-climate system, generated by exogenous integrated assessment models (IAMs). Using IAM scenarios, LiAISON generates a time-series of life cycle inventory (LCI) databases, which are then used to calculate the impacts per functional unit per time step. This expands current practice of using static, future system assumptions, e.g., a specific grid-mix each year. Further, IAM scenarios are provided in a standardized format of shared-socioeconomic pathways (SSP) and representative concentration pathways (RCP) combinations. These are coherent, regularly, published, and peer-reviewed scenario combinations that establish a reproducible and standardized societal and climate mitigation futures context. They are comparable across IAMs and expand the system boundary of the traditional LCA by including dimensions such as societal and behavioral changes. We apply the framework to assess two emerging Power-to-Hydrogen processes, high temperature electrolysis using solid oxide fuel cell (HT-SOE) and polymer electrolyte membrane electrolysis (PEME). We compare the technologies to a baseline Hydrogen production process via steam methane reforming. Despite the decarbonized electricity systems' beneficial effects on the PtH2 processes' carbon intensities, we find environmental tradeoffs, which require technology improvements via learning-by-doing to be alleviated. Future work via ongoing collaborations will focus on linking the framework to other energy-economy-land-climate models and open-source life cycle inventory databases.

ENERGY PLANNING, POLICY, AND ECONOMY↗

Mixed Arlequin method for multiscale poromechanics problems: Mixed Arlequin method for multiscale poromechanics problems

An Arlequin poromechanics model is introduced to simulate the hydro-mechanical coupling effects of fluid-infiltrated porous media across different spatial scales within a concurrent computational framework. A two-field poromechanics problem is first recast as the twofold saddle point of an incremental energy functional. We then introduce Lagrange multipliers and compatibility energy functionals to enforce the weak compatibility of hydro-mechanical responses in the overlapped domain. Here, to examine the numerical stability of this hydro-mechanical Arlequin model, we derive a necessary condition for stability, the twofold inf–sup condition for multi-field problems, and establish a modified inf–sup test formulated in the product space of the solution field. We verify the implementation of the Arlequin poromechanics model through benchmark problems covering the entire range of drainage conditions. Finally, through these numerical examples, we demonstrate the performance, robustness, and numerical stability of the Arlequin poromechanics model.

42 ENGINEERING↗

Representing ductile damage with the dual domain material point method

In this study, we incorporate a ductile damage material model into a computational framework based on the Dual Domain Material Point (DDMP) method. As an example, simulations of a flyer plate experiment involving ductile void growth and material failure are performed. The results are compared with experiments performed on high purity tantalum. We also compare the numerical results obtained from the DDMP method with those obtained from the traditional Material Point Method (MPM). Effects of an overstress model, artificial viscosity, and physical viscosity are investigated. Our results show that a physical bulk viscosity and overstress model are important in this impact and failure problem, while physical shear viscosity and artificial shock viscosity have negligible effects. A simple numerical procedure with guaranteed convergence is introduced to solve for the equilibrium plastic state from the ductile damage model.

36 MATERIALS SCIENCE↗

HPC-Driven Modeling with ML-Based Surrogates for Magnon-Photon Dynamics in Hybrid Quantum System

Here, we introduce a hybrid computational framework that merges HPC-based numerical solvers with physics-informed ML surrogates for efficient modeling of magnon-photon interactions. By running short-duration, high-fidelity Maxwell-LLG simulations and feeding their results into an ML model, we substantially cut simulation time while achieving accurate predictions across larger spatiotemporal domains.

Accuracy↗

Consistent Feature Extraction From Vector Fields: Combinatorial Representations and Analysis Under Local Reference Frames

This dissertation presents research on addressing some of the contemporary challenges in the analysis of vector fields—an important type of scientific data useful for representing a multitude of physical phenomena, such as wind flow and ocean currents. In particular, new theories and computational frameworks to enable consistent feature extraction from vector fields are presented. One of the most fundamental challenges in the analysis of vector fields is that their features are defined with respect to reference frames. Unfortunately, there is no single “correct” reference frame for analysis, and an unsuitable frame may cause features of interest to remain undetected, thus creating serious physical consequences. This work develops new reference frames that enable extraction of localized features that other techniques and frames fail to detect. As a result, these reference frames objectify the notion of “correctness” of features for certain goals by revealing the phenomena of importance from the underlying data. An important consequence of using these local frames is that the analysis of unsteady (time-varying) vector fields can be reduced to the analysis of sequences of steady (timeindependent) vector fields, which can be performed using simpler and scalable techniques that allow better data management by accessing the data on a per-time-step basis. Nevertheless, the state-of-the-art analysis of steady vector fields is not robust, as most techniques are numerical in nature. The residing numerical errors can violate consistency with the underlying theory by breaching important fundamental laws, which may lead to serious physical consequences. This dissertation considers consistency as the most fundamental characteristic of computational analysis that must always be preserved, and presents a new discrete theory that uses combinatorial representations and algorithms to provide consistency guarantees during vector field analysis along with the uncertainty visualization of unavoidable discretization errors. Together, the two main contributions of this dissertation address two important concerns regarding feature extraction from scientific data: correctness and precision. The work presented here also opens new avenues for further research by exploring more-general reference frames and more-sophisticated domain discretizations.

97 MATHEMATICS AND COMPUTING↗

Validating quantum-classical programming models with tensor network simulations

The exploration of hybrid quantum-classical algorithms and programming models on noisy near-term quantum hardware has begun. As hybrid programs scale towards classical intractability, validation and benchmarking are critical to understanding the utility of the hybrid computational model. In this paper, we demonstrate a newly developed quantum circuit simulator based on tensor network theory that enables intermediate-scale verification and validation of hybrid quantum-classical computing frameworks and programming models. We present our tensor-network quantum virtual machine (TNQVM) simulator which stores a multi-qubit wavefunction in a compressed (factorized) form as a matrix product state, thus enabling single-node simulations of larger qubit registers, as compared to brute-force statevector simulators. Our simulator is designed to be extensible in both the tensor network form and the classical hardware used to run the simulation (multicore, GPU, distributed). The extensibility of the TNQVM simulator with respect to the simulation hardware type is achieved via a pluggable interface for different numerical backends (e.g., ITensor and ExaTENSOR numerical libraries). We demonstrate the utility of our TNQVM quantum circuit simulator through the verification of randomized quantum circuits and the variational quantum eigensolver algorithm, both expressed within the eXtreme-scale ACCelerator (XACC) programming model.

97 MATHEMATICS AND COMPUTING↗

GEOS Code Development Road Map

GEOS is a massively parallel computational framework designed to enable HPC-based simulations of subsurface reservoir stimulation activities with the goal of optimizing current operations and evaluating innovative stimulation methods. GEOS will enable coupling of different solvers associated with the various physical processes occurring during reservoir stimulation in unique and sophisticated ways, adapted to various geologic settings, materials and stimulation methods. The overall architecture of the framework includes consistent data structures and will allow incorporation of additional physical and materials models as demanded by future applications. Along with predicting the initiation, propagation and reactivation of fractures, GEOS will also generate a seismic source term that can be linked with seismic wave propagation codes to generate synthetic microseismicity at surface and downhole arrays. Similarly, the output from GEOS can be linked with existing fluid/thermal transport codes. GEOS can also be linked with existing, non-intrusive uncertainty quantification schemes to constrain uncertainty in its predictions and sensitivity to the various parameters describing the reservoir and stimulation operations. We anticipate that an implicit-explicit 3D version of GEOS, including a preliminary seismic source model, will be available for parametric testing and validation against experimental and field data by Oct. 1, 2013.

02 PETROLEUM↗

A Monte Carlo model for 3D grain evolution during welding

Welding is one of the most wide-spread processes used in metal joining. However, there are currently no open-source software implementations for the simulation of microstructural evolution during a weld pass. Here we describe a Potts Monte Carlo based model implemented in the SPPARKS kinetic Monte Carlo computational framework. The model simulates melting, solidification and solid-state microstructural evolution of material in the fusion and heat-affected zones of a weld. The model does not simulate thermal behavior, but rather utilizes user input parameters to specify weld pool and heat-affect zone properties. Weld pool shapes are specified by Bezier curves, which allow for the specification of a wide range of pool shapes. Pool shapes can range from narrow and deep to wide and shallow representing different fluid flow conditions within the pool. Surrounding temperature gradients are calculated with the aide of a closest point projection algorithm. Furthermore, the model also allows simulation of pulsed power welding through time-dependent variation of the weld pool size. Example simulation results and comparisons with laboratory weld observations demonstrate microstructural variation with weld speed, pool shape, and pulsed-power.

36 MATERIALS SCIENCE↗

Prediction of Probabilistic Shock Initiation Thresholds of Energetic Materials Through Evolution of Thermal-Mechanical Dissipation and Reactive Heating

The ignition threshold of an energetic material (EM) quantifies the macroscopic conditions for the onset of self-sustaining chemical reactions. The threshold is an important theoretical and practical measure of material attributes that relate to safety and reliability. Historically, the thresholds are measured experimentally. In this work, we present a new Lagrangian computational framework for establishing the probabilistic ignition thresholds of heterogeneous EM out of the evolutions of coupled mechanical-thermal-chemical processes using mesoscale simulations. Furthermore, the simulations explicitly account for microstructural heterogeneities, constituent properties, and interfacial processes and capture processes responsible for the development of material damage and the formation of hotspots in which chemical reactions initiate. The specific mechanisms tracked include viscoelasticity, viscoplasticity, fracture, post-fracture contact, frictional heating, heat conduction, reactive chemical heating, gaseous product generation, and convective heat transfer. To determine the ignition threshold, the minimum macroscopic loading required to achieve self-sustaining chemical reactions with a rate of reactive heat generation exceeding the rate of heat loss due to conduction and other dissipative mechanisms is determined. Probabilistic quantification of the processes and the thresholds are obtained via the use of statistically equivalent microstructure sample sets (SEMSS). The predictions are in agreement with available experimental data.

36 MATERIALS SCIENCE↗

Neural Network Enhanced RKPM for Electrochemical-Mechanical Coupled Damage Modeling of Energy Storage Materials

Energy storage materials undergo significant charge cycling, which makes understanding their reliability and durability fundamental in predicting performance and service life. Strong electrochemical-mechanical coupling and highly anisotropic material properties contribute to the formation and propagation of micro-cracking, largely along material interfaces and grain boundaries. For Li-ion batteries, for example, lithium moving between electrodes during charging and discharging process causes expansion and contraction of grains, and the strongly anisotropic and nonlinearly [Li]-dependent grain material properties can cause grains to expand into and contract away from each other, leading to chemo-mechanical cracking. In the first part of this work, a RKPM based computational framework for solving the coupled solid-phase lithium conservation with Fickian diffusion and the lithium concentration dependent anisotropic mechanical problem subjected to a highly nonlinear Butler-Volmer boundary condition is introduced. The choice of RKPM completeness conditions for lithium concentration and mechanical deformation fields, and the variational consistency condition for the domain integration of the coupled problem is first determined. In the second part of this work, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1] is leveraged to accurately capture damage and crack propagation throughout the material, by learning the location, orientation, and sharpness of discontinuity while allowing for a coarser nodal distribution than that is necessary for capturing sharp solution transitions using traditional mesh-based methods. NN-RKPM is used to inform how crack opening and closure in turn affect the coupled chemical equations and material microstructure.

damage modeling↗

Improving Uintah's Scalability Through the Use of Portable Kokkos-Based Data Parallel Tasks

The University of Utah’s Carbon Capture Multidisciplinary Simulation Center (CCMSC) is using the Uintah Computational Framework to predict performance of a 1000 MWe ultra-supercritical clean coal boiler. The center aims to utilize the Intel Xeon Phi-based DOE systems, Theta and Aurora, through the Aurora Early Science Program by using the Kokkos C++ library to enable node-level performance portability. This paper describes infrastructure advancements and portability improvements made possible by the integration of Kokkos within Uintah. This integration marks a step towards consolidating Uintah’s MPI+PThreads and MPI+CUDA hybrid parallelism approaches into a single MPI+Kokkos approach. Scalability results are presented that compare serial and data parallel task execution models for a challenging radiative heat transfer calculation, central to the center’s predictive boiler simulations. Here, these results demonstrate both good strong-scaling characteristics to 256 Knights Landing (KNL) processors on the NSF Stampede system, and show the KNL-based calculation to compete with prior GPU-based results for the same calculation.

97 MATHEMATICS AND COMPUTING↗

Multi-modal Energy-optimal Trip Scheduling in Real-time (METS-R) for Transportation Hubs (Final Report)

This report summarizes the work performed under the award number EE0008524. The project develops the Multi-modal Energy-optimal Trip Scheduling in Real-time (METS-R) platform as the next-generation transportation solution based on autonomous electric vehicles (AEV) serving passenger trips from and to urban transportation hubs, to substantially reduce transportation energy consumption. Extensive data collection and analyses were first conducted to understand the demand patterns and energy consumption of hub-based on-road trips. Then, a data-driven framework that consists of an analytical module and a simulation module was proposed. For the analytical module, five planning + operation tools were developed to support the planning and energy-efficient operations of urban AEV services: the charging station planning that robotically allocates charging supplies based on the stationary charging demand distribution; the transit planning and demand adaptive scheduling model that efficiently generates\ candidate transit routes from hubs to other places and dynamically adjusts the transit time table to fit the current demand; the online energy-efficient routing that learns the energy-optimal paths from observations of link-level energy consumption in real-time; the hub-based ridesharing that matches trip requests together with account for the uncertainty of future trip demand and vehicle supply; and finally, the integrated demand prediction and anomaly detection pipeline that leverages the flight/train time table and support other planning/operation tools. To demonstrate the performance of these tools, a scalable high-performance agent-based simulator was built. We divided the urban space into multiple service zones where each zone was considered as an agent for passenger generation and vehicle charging. Two types of AEV agents were coded to model two types of mobility services: AEV taxi and AEV transit. For the AEV taxi, the team implemented the functions of pickup/drop-off passengers, energy-efficient routing, ridesharing, fleet rebalancing, and recharging. For the AEV bus, the team implemented the functions of demand-adaptive route scheduling, passenger boarding, and recharging. A high-performance computing framework was introduced to receive various profiling information (such as link energy updates, vehicle speed) from the simulator instances and communicate the operational commands back to the instances. The numerical experiments show that each of the proposed operational algorithms can reduce energy consumption and improve system efficiency. Furthermore, there exists the need to collectively consider multiple planning + operational strategies as multiple strategies can influence each other in terms of performance impacts. Recommendations for future work related to AEV planning and simulation are discussed.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗