Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulation System”

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

A Large Deviations Analysis of Certain Qualitative Properties of Parallel Tempering and Infinite Swapping Algorithms

Parallel tempering, or replica exchange, is a popular method for simulating complex systems. The idea is to run parallel simulations at different temperatures, and at a given swap rate exchange configurations between the parallel simulations. From the perspective of large deviations it is optimal to let the swap rate tend to infinity and it is possible to construct a corresponding simulation scheme, known as infinite swapping. In this paper we propose a novel use of large deviations for empirical measures for a more detailed analysis of the infinite swapping limit in the setting of continuous time jump Markov processes. Using the large deviations rate function and associated stochastic control problems we consider a diagnostic based on temperature assignments, which can be easily computed during a simulation. Here we show that the convergence of this diagnostic to its a priori known limit is a necessary condition for the convergence of infinite swapping. Finally, the rate function is also used to investigate the impact of asymmetries in the underlying potential landscape, and where in the state space poor sampling is most likely to occur.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Lazy Updating of hubs can enable more realistic models by speeding up stochastic simulations

To respect the nature of discrete parts in a system, stochastic simulation algorithms (SSAs) must update for each action (i) all part counts and (ii) each action's probability of occurring next and its timing. This makes it expensive to simulate biological networks with well-connected “hubs” such as ATP that affect many actions. Temperature and volume also affect many actions and may be changed significantly in small steps by the network itself during fever and cell growth, respectively. Such trends matter for evolutionary questions, as cell volume determines doubling times and fever may affect survival, both key traits for biological evolution. Yet simulations often ignore such trends and assume constant environments to avoid many costly probability updates. Such computational convenience precludes analyses of important aspects of evolution. Here we present “Lazy Updating,” an add-on for SSAs designed to reduce the cost of simulating hubs. When a hub changes, Lazy Updating postpones all probability updates for reactions depending on this hub, until a threshold is crossed. Speedup is substantial if most computing time is spent on such updates. We implemented Lazy Updating for the Sorting Direct Method and it is easily integrated into other SSAs such as Gillespie's Direct Method or the Next Reaction Method. Testing on several toy models and a cellular metabolism model showed >10× faster simulations for its use-cases—with a small loss of accuracy. Thus we see Lazy Updating as a valuable tool for some special but important simulation problems that are difficult to address efficiently otherwise.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A coupled human–natural system analysis of freshwater security under climate and population change

Limited water availability, population growth, and climate change have resulted in freshwater crises in many countries. Jordan’s situation is emblematic, compounded by conflict-induced population shocks. Integrating knowledge across hydrology, climatology, agriculture, political science, geography, and economics, we present the Jordan Water Model, a nationwide coupled human–natural-engineered systems model that is used to evaluate Jordan’s freshwater security under climate and socioeconomic changes. The complex systems model simulates the trajectory of Jordan’s water system, representing dynamic interactions between a hierarchy of actors and the natural and engineered water environment. A multiagent modeling approach enables the quantification of impacts at the level of thousands of representative agents across sectors, allowing for the evaluation of both systemwide and distributional outcomes translated into a suite of water-security metrics (vulnerability, equity, shortage duration, and economic well-being). Model results indicate severe, potentially destabilizing, declines in freshwater security. Per capita water availability decreases by approximately 50% by the end of the century. Without intervening measures, >90% of the low-income household population experiences critical insecurity by the end of the century, receiving <40 L per capita per day. Widening disparity in freshwater use, lengthening shortage durations, and declining economic welfare are prevalent across narratives. To gain a foothold on its freshwater future, Jordan must enact a sweeping portfolio of ambitious interventions that include large-scale desalinization and comprehensive water sector reform, with model results revealing exponential improvements in water security through the coordination of supply- and demand-side measures.

42 ENGINEERING↗
DOE OSTI2017

A New Modeling Technique to Analyze Safeguards Measurements in Large Systems

SafeGuards Analysis (SGA) is a computational toolbox able to simulate different safeguards scenarios across a number of different fuel cycles and at many different scales within the Matlab Simulink framework. SGA functions by simulating Material Balance Areas (MBAs) under safeguards materials control and accountability and allows the user to define the uncertainty parameters of the associated flow and inventory measurements. The simulated safeguard system uses the uncertain measurement estimates to calculate a mass-balance across the MBA. This mass balance is then evaluated by one of a number of different statistical tests to determine if a significant amount of material has been removed from the MBA. This work describes the design of SGA, the results of testing each element of the toolbox, and a number of single MBA example scenarios. In all of the test cases SGA performed as expected, and produced acceptable results from the single MBA scenario.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

RouteE: A Vehicle Energy Consumption Prediction Engine

The emergence of connected and automated vehicles and smart cities technologies create the opportunity for new mobility modes and routing decision tools, among many others. To achieve maximum mobility and minimum energy consumption, it is critical to understand the energy cost of decisions and optimize accordingly. The Route Energy prediction model (RouteE) enables accurate estimation of energy consumption for a variety of vehicle types over trips or sub-trips where detailed drive cycle data are unavailable. Applications include vehicle route selection, energy accounting and optimization in transportation simulation, and corridor energy analyses, among others. The software is a Python package that includes a variety of pre-trained models from the National Renewable Energy Laboratory (NREL). However, RouteE also enables users to train custom models using their own data sets, making it a robust and valuable tool for both fast calculations and rigorous, data-rich research efforts. The pre-trained RouteE models are established using NREL’s Future Automotive Systems Technology Simulator paired with approximately 1 million miles of drive cycle data from the Transportation Secure Data Center, resulting in energy consumption behavior estimates over a representative sample of driving conditions for the United States. Validations have been performed using on-road fuel consumption data for conventional and electrified vehicle powertrains. Transferring the results of the on-road validation to a larger set of real-world origin-destination pairs, it is estimated that implementing the present methodology in a green-routing application would accurately select the route that consumes the least fuel 90% of the time. The novel machine learning techniques used in RouteE make it a flexible and robust tool for a variety of transportation applications.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2016

Optical spectroscopy of molecular junctions: Nonequilibrium Green’s functions perspective

We consider optical spectroscopy of molecular junctions from the quantum transport perspective when radiation field is quantized and optical response of the system is simulated as photon flux. Using exact expressions for photon and electronic fluxes derived within the nonequilibrium Green function (NEGF) methodology and utilizing fourth order diagrammatic perturbation theory (PT) in molecular coupling to radiation field, we perform simulations employing realistic parameters. Results of the simulations are compared to the bare PT which is usually employed in studies on nonlinear optical spectroscopy to classify optical processes. We show that the bare PT violates conservation laws, while flux conserving NEGF formulation mixes optical processes.

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

Km‐Scale Simulations of Mesoscale Convective Systems Over South America—A Feature Tracker Intercomparison

Mesoscale convective systems (MCSs) are clusters of thunderstorms that are important in Earth's water and energy cycle. Additionally, they are responsible for extreme events such as large hail, strong winds, and extreme precipitation. Automated object-based analyses that track MCSs have become popular since they allow us to identify and follow MCSs over their entire life cycle in a Lagrangian framework. This rise in popularity was accompanied by an increasing number of MCS tracking algorithms, however, little is known about how sensitive analyses are concerning the MCS tracker formulation. Here, we assess differences between six MCS tracking algorithms on South American MCS characteristics and evaluate MCSs in kilometer-scale simulations with observational-based MCSs over 3 years. All trackers are run with a common set of MCS classification criteria to isolate tracker formulation differences. The tracker formulation substantially impacts MCS characteristics such as frequency, size, duration, and contribution to total precipitation. The evaluation of simulated MCS characteristics is less sensitive to the tracker formulation and all trackers agree that the model can capture MCS characteristics well across different South American climate zones. Dominant sources of uncertainty are the segmentation of cloud systems in space and time and the treatment of how MCSs are linked in time. Our results highlight that comparing MCS analyses that use different tracking algorithms is challenging. We provide general guidelines on how MCS characteristics compare between trackers to facilitate a more robust assessment of MCS statistics in future studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Chitah: Strong-gravitational-lens hunter in imaging surveys

Strong gravitationally lensed quasars provide powerful means to study galaxy evolution and cosmology. Current and upcoming imaging surveys will contain thousands of new lensed quasars, augmenting the existing sample by at least two orders of magnitude. To find such lens systems, we built a robot, Chitah, that hunts for lensed quasars by modeling the configuration of the multiple quasar images.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

PowerSystems.jl and PowerSimulations.jl [SWR-18-23]

Data structures in Julia to enable power systems analysis. Also, PowerSimulations.jl is a Julia package for power system modeling and simulation found at https://github.com/NREL/PowerSimulations.jl. The objectives of the package are: Provide a flexible modeling framework that can accommodate problems of different complexity and at different time-scales. Streamline the construction of large scale optimization problems to avoid repetition of work when adding/modifying model details. Exploit Julia's capabilities to improve computational performance of large scale power system simulations. The flexible modeling framework is enabled through a modular set of capabilities that enable scalable power system analysis and exploration of new analysis methods. The modularity of PowerSimulations results from the structure of the simulations enabled by the package: Simulations define a set of problems that can be solved using numerical techniques. Problems are generated by expressing model formulations against system data For example, an annual production cost modeling simulation can be created by formulating a unit commitment model against system data to assemble a set of 365 daily time-coupled scheduling problems.

Lara, Jose Daniel↗
DOE OSTI2019

Evaluating the Impact of Price-Responsive Load on Power Systems Using Integrated T&D Simulation

This paper explores the differences between simulating price-responsive load (PRL) interactions with power systems using integrated transmission and distribution (T&D) models and transmission-only (T-only) models. This analysis uses the Integrated Grid Modeling System (IGMS) software to capture 'ISO-to-appliance' simulations using a synthetic T&D model built on the PJM 5-Bus transmission network with multiple full-scale taxonomy feeders that include physics-based models of thousands of customers and PRLs. The results show important differences in the impacts of PRLs between integrated T&D and T-only models. Experiments with the synthetic integrated T&D dataset demonstrated that integrated T&D simulation revealed notably larger differences between the PRL and no-PRL cases for load, and prices compared to T-only simulation. Similarly, differences are observed between the price response of individual buildings and distribution feeders and the corresponding transmission bus in the integrated T&D simulations, which are difficult to capture in traditional T-only simulations.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Surrogate Constructed Scalable Circuits ADAPT-VQE in the Schwinger model

Inspired by recent advancements of simulating periodic systems on quantum computers, we develop a new approach, (SC)$^2$-ADAPT-VQE, to further advance the simulation of these systems. Our approach extends the scalable circuits ADAPT-VQE framework, which builds an ansatz from a pool of coordinate-invariant operators defined for arbitrarily large, though not arbitrarily small, volumes. Our method uses a classically tractable ``Surrogate Constructed'' method to remove irrelevant operators from the pool, reducing the minimum size for which the scalable circuits are defined. Bringing together the scalable circuits and the surrogate constructed approaches forms the core of the (SC)$^2$ methodology. Our approach allows for a wider set of classical computations, on small volumes, which can be used for a more robust extrapolation protocol. While developed in the context of lattice models, the surrogate construction portion is applicable to a wide variety of problems where information about the relative importance of operators in the pool is available. As an example, we use it to compute properties of the Schwinger model - quantum electrodynamics for a single, massive fermion in $1+1$ dimensions - and show that our method can be used to accurately extrapolate to the continuum limit.

Gustafson, Erik [RIACS, Mtn. View] (ORCID:00000001↗
DOE OSTI2016

TH-AB-BRA-08: Simulated Tumor Tracking in An MRI Linac for Lung Tumor Lesions Using the Monaco Treatment Planning System

Purpose: To simulate tumor tracking in an Elekta MRI-linac (MRL) and to compare this tracking method with our current ITV approach in terms of OAR sparing for lung cancer patients. Methods: Five SABR-NSCLC patients with central lung tumors were selected for reasons of potential enhancement of tumor-tissue delineation using MRI. The Monaco TPS was used to compare the current clinical ITV approach to a simulated, novel tracking method which used a 7MV MRL beam in the presence of an orthogonal 1.5 T magnetic field (4D-MRL method). In the simulated tracking scenario, achieved using the virtual couch shift (VCS), the PTV was defined using an isotropic 5mm margin applied to the GTV of each phase, as acquired from an 8-phase amplitude-binned 4DCT. These VCS plans were optimized and weighted on each phase. The dose weighting was performed using the patient-specific breathing traces. The doses were accumulated on the inhale phase. The two methods were compared by assessing the OAR DVHs. Results: The 4D-MRL method resulted in a reduced target volume (by an average of 29% over all patients). The benefits of using an MRL tracking system depended on the tumor motion amplitude and the relative OAR motion (ROM) to the target. The reduction in mean doses to parallel organs was up to 3 Gy for the heart and 2.1 Gy for the lung. The reductions in maximum doses to serial organs were up to 9.4 Gy, 5.6 Gy, and 8.7 Gy for the esophagus, spinal cord, and the trachea, respectively. Serial organs benefited from MRL tracking when the ROM was ≥ 0.3 cm despite small tumor motion amplitude in some cases. Conclusions: This work demonstrated the potential benefit for an MRL tracking system to spare OARs in SABR-NSCLC patients with central tumors. The benefits are embodied in the target volume reduction. This project was made possible with the financial support of Elekta.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

FY2020 Fourth Quarter Performance Metric: Evaluate Improvement in Simulations of Mesoscale Convective Systems from New Parameterization Developments in E3SM

Mesoscale convective systems (MCSs) consist of an assembly of cumulonimbus clouds on scales of 100 km or more and produce mesoscale circulations (Houze, 2004, 2018). As the largest form of deep convective storms, MCSs contribute to 30% – 70% of annual and warm season rainfall in the U.S. and in the global tropics (Houze 2018; Stevenson & Schumacher, 2014; Feng et al., 2019; Haberlie & Ashley, 2019). Since MCSs contribute importantly to mean and extreme precipitation in the U.S. and many other regions around the world, understanding how well they are simulated by E3SM may guide future development towards more skillful modeling of convective storms and associated hydrologic impacts. The FY2020 Second Quarter Performance Metric Report documented comparisons of MCSs in the central and eastern U.S. in a high-resolution simulation produced by E3SM v1 at 25 km resolution (Caldwell et al. 2019) with observations. MCSs in the simulation occur less frequently and produce less intense precipitation, resulting in large underestimation of MCS volumetric rain-rate compared to observations. The first and third quarter performance metric report indicated that these model biases in simulating MCSs can be attributed to model limitations in parameterizing convection, clouds, and other related processes, as well as model biases in simulating the MCS large-scale environment. In the current FY2020 Fourth Quarter Performance Metric Report, we evaluate MCSs simulated in E3SM with several new developments in convection parameterizations that are being developed for its next generation. The goal is to summarize what have been improved with the new developments and highlight what need more work in the future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Challenges & Roadmap for Beyond CMOS Computing Simulation.

Simulating HPC systems is a difficult task and the emergence of “Beyond CMOS” architectures and execution models will increase that difficulty. This document presents a “tutorial” on some of the simulation challenges faced by conventional and non-conventional architectures (Section 1) and goals and requirements for simulating Beyond CMOS systems (Section 2). These provide background for proposed short- and long-term roadmaps for simulation efforts at Sandia (Sections 3 and 4). Additionally, a brief explanation of a proof-of-concept integration of a Beyond CMOS architectural simulator is presented (Section 2.3).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Optimal Realization of Yang–Baxter Gate on Quantum Computers

Quantum computers provide a promising method to study the dynamics of many-body systems beyond classical simulation. On the other hand, the analytical methods developed and results obtained from the integrable systems provide deep insights on the many-body system. Quantum simulation of the integrable system not only provides a valid benchmark for quantum computers but is also the first step in studying integrable-breaking systems. The building block for the simulation of an integrable system is the Yang–Baxter gate. It is vital to know how to optimally realize the Yang–Baxter gates on quantum computers. Based on the geometric picture of the Yang–Baxter gates, the optimal realizations of two types of Yang–Baxter gates with a minimal number of controlled NOT (CNOT) or gates are presented. It is also shown how to systematically realize the Yang–Baxter gates via the pulse control. The different realizations on IBM quantum computers are tested and compared. It is found that the pulse realizations of the Yang–Baxter gates always have a higher gate fidelity compared to the optimal CNOT or realizations. On the basis of the above optimal realizations, the simulation of the Yang–Baxter equation on quantum computers is demonstrated. Finally, these results provide a guideline and standard for further experimental studies based on the Yang–Baxter gate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Large-scale molecular dynamics simulations of bubble collapse in water: Effects of system size, water model, and nitrogen

Molecular dynamics simulations in the microcanonical ensemble are performed to study the collapse of a bubble in liquid water using the single-site mW and the four-site TIP4P/2005 water models. To study system size effects, simulations for pure water systems are performed using periodically replicated simulation boxes with linear dimensions, L, ranging from 32 to 512 nm with the largest systems containing 8.7 × 10 6 and 4.5 × 10 9 molecules for the TIP4P/2005 and mW water models, respectively. The computationally more efficient mW water model allows us to reach converging behavior when the bubble dynamics results are plotted in reduced units, and the limiting behavior can be obtained through linear extrapolation in L –1 . Qualitative differences are observed between simulations with the mW and TIP4P/2005 water models, but they can be explained by the models’ differences in predicted viscosity and surface tension. Although bubble collapse occurs on time scales of only hundreds of picoseconds, the system sizes used here are sufficiently large to obtain bubble dynamics consistent with the Rayleigh–Plesset equation when using the models’ thermophysical properties as input. For the conditions explored here, extreme heating of the interfacial water molecules near the time of collapse is observed for the larger mW water systems (but the model underpredicts the viscosity), whereas heating is less pronounced for the TIP4P/2005 water systems because its larger viscosity contribution slows the collapse dynamics. The presence of nitrogen within the bubble only starts to affect bubble dynamics near the very end of the initial collapse, leading to an incomplete collapse and strong rebound for the mW water model. Although nitrogen is non-condensable at 300 K, it becomes highly compressed and reaches a liquid-like density near the collapse point. We find that the dissolution of nitrogen is much slower than the movement of the collapsing water front, and the re-expansion of the dense nitrogen droplet gives rise to bubble rebound. The incompatibility of the collapse and dissolution time scales should be considered for continuum-scale modeling of bubble dynamics. Finally, we also confirm that the diffusion coefficient for dissolved nitrogen is insensitive to pressure as the liquid transitions from a compressed to a stretched state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗