Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large system simulation”

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

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

At least 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Large-scale simulations of Floquet physics on near-term quantum computers

Abstract Periodically driven quantum systems exhibit a diverse set of phenomena but are more challenging to simulate than their equilibrium counterparts. Here, we introduce the Quantum High-Frequency Floquet Simulation (QHiFFS) algorithm as a method to simulate fast-driven quantum systems on quantum hardware. Central to QHiFFS is the concept of a kick operator which transforms the system into a basis where the dynamics is governed by a time-independent effective Hamiltonian. This allows prior methods for time-independent simulation to be lifted to simulate Floquet systems. We use the periodically driven biaxial next-nearest neighbor Ising (BNNNI) model, a natural test bed for quantum frustrated magnetism and criticality, as a case study to illustrate our algorithm. We implemented a 20-qubit simulation of the driven two-dimensional BNNNI model on Quantinuum’s trapped ion quantum computer. Our error analysis shows that QHiFFS exhibits not only a cubic advantage in driving frequency ω but also a linear advantage in simulation time t compared to Trotterization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A first-principles machine-learning force field for heterogeneous ice nucleation on microcline feldspar

The formation of ice in the atmosphere affects precipitation and cloud properties, and plays a key role in the climate of our planet. Although ice can form directly from liquid water at deeply supercooled conditions, the presence of foreign particles can aid ice formation at much warmer temperatures. Over the past decade, experiments have highlighted the remarkable efficiency of feldspar minerals as ice nuclei compared to other particles present in the atmosphere. However, the exact mechanism of ice formation on feldspar surfaces has yet to be fully understood. Here, we develop a first-principles machine-learning model for the potential energy surface aimed at studying ice nucleation at microcline feldspar surfaces. The model is able to reproduce with high-fidelity the energies and forces derived from density-functional theory (DFT) based on the SCAN exchange and correlation functional. Our training set includes configurations of bulk supercooled water, hexagonal and cubic ice, microcline, and fully-hydroxylated feldspar surfaces exposed to vacuum, liquid water, and ice. In this work, we apply the machine-learning force field to study different fully-hydroxylated terminations of the (100), (010), and (001) surfaces of microcline exposed to vacuum. Our calculations suggest that terminations that do not minimize the number of broken bonds are preferred in vacuum. We also study the structure of supercooled liquid water in contact with microcline surfaces, and find that water density correlations extend up to around 10 Å from the surfaces. Finally, we show that the force field maintains a high accuracy during the simulation of ice formation at microcline surfaces, even for large systems of around 30,000 atoms. Future work will be directed towards the calculation of nucleation free energy barriers and rates using the force field developed herein, and understanding the role of different microcline surfaces on ice nucleation.

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

Designing a quantum-accurate machine-learning potential to enable large-scale simulations of deuterium under shock

Large-scale molecular dynamics of deuterium under shock can elucidate kinetic processes vital to the target design in inertial confinement fusion and high-energy-density experiments. However, modeling the complex evolution of this material from an insulating molecular state at ambient pressure to an ionized, atomic fluid under strong shock is beyond the capability of simple pair and even bond order potentials. We thus train a quantum-accurate and broadly transferable machine-learning interatomic potential for deuterium using the Chebyshev Interaction Model for Efficient Simulations framework. We show that due to an improved description of the molecular-to-atomic transition, our model is able to better reproduce the ab initio equation of state, radial distribution functions, and principal Hugoniot than bond order potentials. This represents an important step toward large-scale quantum-accurate and nonequilibrium simulations of complicated systems under dynamic changes including phase transitions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Population inversion X-ray laser oscillator

Oscillators are at the heart of optical lasers, providing stable, transform-limited pulses. Until now, laser oscillators have been available only in the infrared to visible and near-ultraviolet (UV) spectral region. In this paper, we present a study of an oscillator operating in the 5- to 12-keV photon-energy range. We show that, using the Kα 1 line of transition metal compounds as the gain medium, an X-ray free-electron laser as a periodic pump, and a Bragg crystal optical cavity, it is possible to build X-ray oscillators producing intense, fully coherent, transform-limited pulses. As an example, we consider the case of a copper nitrate gain medium generating ~ 5 × 10 10 photons per pulse with 37-fs pulse length and 48-meV spectral resolution at 8-keV photon energy. Our theoretical study and simulation of this system show that, because of the very large gain per pass, the oscillator saturates and reaches full coherence in four to six optical-cavity transits. We discuss the feasibility and design of the X-ray optical cavity and other parts of the oscillator needed for its realization, opening the way to extend X-ray–based research beyond current capabilities.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Quasi-Static Time Series Fatigue Simulation for PV Inverter Semiconductors with Long-Term Solar Profile

Power system simulations with long-term data tend to have large time steps varying from one second to a few minutes. However, for PV inverter semiconductors, the minimum thermal stresses cycle is with line frequency. This requires the time step of the fatigue simulation to be much smaller than the line period. This small time step results in poor simulation speed, especially for long-term simulations. This paper proposes a fast fatigue simulation for inverter semiconductors using the quasi-static time series (QSTS) simulation concept. The fatigue analysis typically focuses on the peak and valley values of a strain and neglect the transients from peaks to valleys. The proposed simulation utilizes this property of fatigue analysis and calculates the steady state of the semiconductor junction temperature only. The resulting time step of the fatigue simulation is 15 minutes, which is consistent with the solar dataset without losing accuracy.

fast simulation↗
DOE OSTI2023

Spatial Microsimulation and Activity Allocation in Python: An Update on the Likeness Toolkit

Understanding human security and social equity issues within human systems requires large-scale models of population dynamics that simulate high-fidelity representations of individuals and access to essential activities (work/school, social, errands, health). Likeness is a Python toolkit that provides these capabilities for Oak Ridge National Laboratory's (ORNL) UrbanPop spatial microsimulation project. In step with the initial development phase for Likeness (2021 - 2022), we built out several foundational examples of work/school and health service access. In this paper, we describe expansion and scaling of Likeness capabilities to metropolitan areas in the United States. We then provide an integrated demonstration of our methods based on a case study of Leon County, FL and perform validation exercises on 1) neighborhood demographic composition and 2) visits by demographic cohorts (gender/age) obtained from point of interest (POI) footfall data for essential services (grocery stores). Taking into account lessons learned from our case study, we scope improvements to our model as well as provide a roadmap of the anticipated Likeness development cycle into 2023 - 2024.

Tuccillo, Joe↗
DOE OSTI2023

Series FACTS Devices for Increasing Resiliency in Severe Weather Conditions

Severe weather conditions are low-probability, high-impact events that affect grid operations. The majority of power outages are caused by severe weather conditions. Grid resiliency to weather events can be enhanced by decreasing the reliance on its affected sections. One way to do this is to reduce the power flow through lines vulnerable to severe weather. If a line is disconnected, its initial power flow is distributed through the neighbor lines, which may cause congestion in the grid. FACTS devices can be used to control the power flow of lines that have a higher chance of power outages. Most previous works do not consider weather events in power flow control. In this work, a linearized optimal power flow (OPF)–based algorithm is developed to minimize the real power flow of vulnerable lines considering the thermal limits of lines to prevent infeasible solutions; the simulation is fast, making it suitable for large-scale systems. The proposed optimization problem is presented as a mixed-integer linear program (MILP), making it capable of using short-term load forecasting due to its high solution speed. The proposed optimization problem considers multiple lines with different outage probabilities and the uncertainties of the weather forecast. Moreover, it estimates the power reduction in vulnerable lines due to changes in the series FACTS devices. The performance of the proposed optimization problem is tested on IEEE 14-, 30-, and 118-bus systems for several scenarios. The results are validated with the AC power flow results from MATPOWER.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Enhancing Local Grid Resilience with Small Hydropower Hybrids: Proving the concept through demonstration, simulation, and analysis with Idaho Falls Power

Large hydropower, connected to the transmission system and typically possessing significant ability to balance grid frequency, has long been central to black start plans in regions where it is present. Small hydropower possesses most of the attributes required for black start but is often connected to distribution or sub-transmission systems and has less ability to balance grid frequency. Integrating energy storage such as batteries or ultracapacitors increases the combined asset’s ability to balance frequency. This asset integration enables a bottom-up grid restoration paradigm in which critical electric loads on the local distribution system can be powered even when the regional transmission system is down. This report presents a field demonstration conducted with Idaho Falls Power that proved this concept. The contribution of energy storage in restoring the small hydro-dominated distribution grid and its operational sensitivity across different control settings are further analyzed using high-fidelity simulations. Readers will learn about technical details on the field demonstration setup, energy storage contribution to black start, detailed simulation steps of islanded distribution grid restoration, and usage of field demonstration measurements for transient model refinements. Collectively this report points to a great opportunity for small hydropower to enhance resilience of local electric grids.

13 HYDRO ENERGY↗
DOE OSTI2021

Discrete Empirical Interpolation Method Based Dynamic Load Model Reduction

Dynamic load models add significant complexity to bulk power system time-domain simulations. The complexity is due to the large number of ordinary differential equations (ODEs) introduced by the dynamic load components such as induction motors. It is challenging to derive reduced-order models (ROMs) for dynamic loads due to the nonlinear functions in their governing equations. This paper applies the discrete empirical interpolation method enhanced proper orthogonal decomposition (DEIM-POD) to approximate the full dynamic load model with the ROM that minimizes the projection error of the nonlinear functions in dynamic load ODEs onto their dominant modes. This approach only requires evaluation of nonlinear functions at selected observation points. The observation points selected by DEIM also provide information for screening critical load buses where dynamic load model parameters contribute the most to the accuracy of ROM across multiple contingencies. The proposed approach is validated on IEEE 9-bus, WECC 179-bus and 2384-bus Polish systems.

bulk power system↗
DOE OSTI2025

AutonomieAI: An efficient and deployable vehicle energy consumption estimation toolkit

Here, this paper presents AutonomieAI, a novel toolkit designed for efficient energy estimation of vehicles across diverse trip scenarios, routes, and drive cycles, applicable to a broad range of vehicle powertrain technologies. It leverages state-of-the-art Machine Learning techniques to deliver real-time energy prediction of vehicles, enabling co-simulation with transportation level system tools and opening doors for large-scale optimization at city, network or national level. Benchmark results show that AutonomieAI achieves high accuracy, with an average percentage error below 2% for most powertrain types, and computational efficiency capable of processing over 10,000 trips per second. Applications of AutonomieAI have potential to offer the flexibility to assist in solving eco-routing problems, optimize for vehicle and powertrain selection, study charging decision behavior, and optimize for charging station placement. AutonomieAI is the result of large neural network based model architectures, trained on very large and unique high fidelity vehicle simulation data. It is lightweight, deployable, efficient and has accuracy comparable to specialized and complex physics based simulation softwares.

Autonomie↗
DOE OSTI2024

Evaluation of Autoconversion Representation in E3SMv2 Using an Ensemble of Large-Eddy Simulations of Low-Level Warm Clouds

In numerical atmospheric models that treat cloud and rain droplet populations as separate condensate categories, precipitation initiation in warm clouds is often represented by an autoconversion rate (Au), which is the rate of formation of new rain droplets through the collisions of cloud droplets. Being a function of the cloud droplet size distribution (DSD), the local Au is commonly parameterized as a function of DSD moments: cloud droplet number (n c ) and mass (q c ) concentrations. When applied in a large-scale model, the grid-mean Au must also include a correction, or enhancement factor, to account for the horizontal variability of the cloud properties across the model grid. In this study, we evaluate the Au representation in the Energy Exascale Earth System Model version 2 (E3SMv2) climate model using large-eddy simulations (LES), which explicitly resolve cloud droplet spectra, and therefore the local Au, as well as its spatial variability. The analysis of an ensemble of warm low-level cloud cases shows that the E3SMv2 formulation represents the Au reasonably well compared to the horizontally averaged explicitly computed rate from LES. The agreement, however, comes from a combination of an underestimated E3SM-tuned local Au rate and an overestimated subgrid cloud variability enhancement factor. The latter bias is traced to neglecting the horizontal variability of n c and its co-variability with q c in parameterizing the grid-mean Au.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

EMT-TS Hybrid Simulation for Large Power Grids Considering IBR-Driven Dynamics

The escalating integration of inverter-based resources (IBRs) poses new challenges to power systems by introducing fast dynamics with higher frequencies, which may need to be simulated by an electromagnetic transient (EMT) program. As an alternative to conducting EMT simulations for the entire system, which is typically time consuming, hybrid simulation between EMT and phasor-domain transient stability (TS) can greatly reduce the computational burden while preserving the detailed fast dynamics in the EMT zone. This paper establishes an EMT-TS hybrid simulation platform using open-source tools, specifically ParaEMT, GridPACK, and HELICS, which are the EMT simulator, TS simulator, and interface framework, respectively. Case studies on the 240-bus Western Electricity Coordinating Council (WECC) system demonstrate that the developed ParaEMT-HELICS-GridPACK hybrid simulator can accurately capture both slow electromechanical and fast IBR-driven dynamics with a 2.4x speedup.

electromagnetic transient↗
DOE OSTI2025

Superspin renormalization and slow relaxation in random spin systems

We develop an excited-state real-space renormalization group (RSRG-X) formalism to describe the dynamics of conserved densities in randomly interacting spin-12 systems. Our formalism is suitable for systems with U(1) and Z2 symmetries, and we apply it to chains of randomly positioned spins with dipolar XX+YY interactions, as arise in Rydberg quantum simulators and other platforms. The formalism generates a sequence of effective Hamiltonians that provide approximate descriptions for dynamics on successively smaller energy scales. These effective Hamiltonians involve “superspins”: two-level collective degrees of freedom constructed from (anti)aligned microscopic spins. Conserved densities can then be understood as relaxing via coherent collective spin flips. For the well-studied simpler case of randomly interacting nearest-neighbor XX+YY chains, the superspins reduce to single spins. Our formalism also leads to a numerical method capable of simulating the dynamics up to an otherwise inaccessible combination of large system size and late time. Focusing on disorder-averaged infinite-temperature autocorrelation functions, in particular the spin survival probability Sp¯(t), we demonstrate quantitative agreement between our algorithm and exact diagonalization (ED) at low but nonzero frequencies. Such agreement holds for chains with nearest-neighbor, next-nearest-neighbor, and long-range dipolar interactions. Our results indicate decay of Sp¯(t) slower than any power law and feature no significant deviation from the ∼1/ln2(t) asymptote expected from the infinite-randomness fixed-point of the nearest-neighbor model. We also apply the RSRG-X formalism to two-dimensional long-range systems of moderate size and find slow late-time decay of Sp¯(t).

Zhao, Yi J↗
DOE OSTI2023

Strong Dark Matter Self-interactions Diversify Halo Populations within and surrounding the Milky Way

Abstract We perform a high-resolution cosmological zoom-in simulation of a Milky Way (MW)–like system, which includes a realistic Large Magellanic Cloud analog, using a large differential elastic dark matter self-interaction cross section that reaches ≈100 cm 2 g −1 at relative velocities of ≈10 km s −1 , motivated by the diverse and orbitally dependent central densities of dwarf galaxies within and surrounding the MW. We explore the effects of dark matter self-interactions on satellite, splashback, and isolated halos through their abundance, central densities, maximum circular velocities, orbital parameters, and correlations between these variables. We use an effective constant cross section model to analytically predict the stages of our simulated halos’ gravothermal evolution, demonstrating that deviations from the collisionless R max – V max relation can be used to select deeply core-collapsed halos, where V max is a halo’s maximum circular velocity, and R max is the radius at which it occurs. We predict that a sizable fraction (≈20%) of subhalos with masses down to ≈10 8 M ⊙ is deeply core collapsed in our SIDM model. Core-collapsed systems form ≈10% of the isolated halo population down to the same mass; these isolated, core-collapsed halos would host faint dwarf field galaxies with extremely steep central density profiles. Finally, most halos with masses above ≈10 9 M ⊙ are core-forming in our simulation. Our study thus demonstrates how self-interactions diversify halo populations in an environmentally dependent fashion within and surrounding MW-mass hosts, providing a compelling avenue to address the diverse dark matter distributions of observed dwarf galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

MIRACL Co‐Simulation platform for control and operation of distributed wind in microgrid

Abstract This paper presents a co‐simulation platform (CSP) developed as a part of the “Microgrids, Infrastructure Resilience and Advanced Controls Launchpad (MIRACL)” project, hereafter called MIRACL‐CSP, to allow simulation‐based testing, demonstration, and evaluations of distributed wind under various grid operating conditions. MIRACL‐CSP provides modular interactions among the power distribution system, distributed wind, and utility decision‐making framework. A hierarchical engine for large‐scale infrastructure co‐simulations (HELICS) is used as the core engine of MIRACL‐CSP to establish time and information coordination among the MIRACL‐CSP modules. In this study, MIRACL‐CSP performance is demonstrated using the IEEE 123‐node test distribution grid modelled in GridLAB‐D and the utility decision support application modelled in Python. The functionality of MIRACL‐CSP is demonstrated through various grid operational scenarios in a microgrid and networked microgrid environment. A co‐simulation study is described that demonstrated MIRACL‐CSP capabilities for the microgrid operational scenario in the IEEE 123‐node test system. MIRACL‐CSP is a generic platform that facilitates distributed wind research for different test systems, applications, and valuation of distributed wind.

17 WIND ENERGY↗
DOE OSTI2023

Nuclear–Electronic Orbital Quantum Mechanical/Molecular Mechanical Real-Time Dynamics

Simulating the nuclear–electronic quantum dynamics of large-scale molecular systems in the condensed phase is key for studying biologically and chemically important processes such as proton transfer and proton-coupled electron transfer reactions. Herein, the real-time nuclear–electronic orbital time-dependent density functional theory (RT-NEO-TDDFT) approach is combined with a hybrid quantum mechanical/molecular mechanical (QM/MM) strategy to enable the accurate description of coupled nuclear–electronic quantum dynamics in the presence of heterogeneous environments such as solvent or proteins. The densities of the electrons and quantum protons are propagated in real time, while the other nuclei are propagated classically on the instantaneous electron–proton vibronic surface. This approach is applied to phenol bound to lysozyme, intramolecular proton transfer in malonaldehyde, and nonequilibrium excited-state intramolecular proton transfer in o-hydroxybenzaldehyde. Furthermore, these examples illustrate that the RT-NEO-TDDFT framework, coupled with an atomistic representation of the environment, allows the simulation of condensed-phase systems that exhibit significant nuclear quantum effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Towards a Techno-Economic Analysis of PCM Integrated Hybrid HVAC Systems

Thermal end uses dominate building energy consumption and are a major driver of peak demand. As heating is electrified, peak electrical power required will surge, prompting a need for innovative HVAC system designs and controls.These designs must incorporate novel technologies at the component level and new integration techniques at the system level. One such possibility involves the addition of thermal energy storage (TES) in heating and cooling equipment using a phase change material (PCM) heat exchanger. Here, thermal energy storage via phase change can be used to shift the HVAC system loads to times of lower electricity cost, reduced carbon intensity, and greater energy efficiency. Most of the current utilization of PCM in buildings involves passive components. By actively controlling when heat is stored and released from PCM, we can optimize the building HVAC system to cost-effectively meet consumer needs with the flexibility to draw on renewable energy resources when they are abundant and available. While this combination of technologies is promising in theory, simulation-based evaluation of a prototype can be difficult due to the modeling requirements at the component level and the large number of possible configurations and operating modes at the system level. To conduct this evaluation, we use the Modelica language for modeling and simulation because it enables users to represent the important physics of the problem, interchange and rearrange components in an efficient manner, and implement a range of control configurations. In this work, we considered three case studies: a portable building, a large commercial retail store, and a multifamily residential apartment unit. Each of these employs a different system design, ranging from a single package vertical unit incorporating PCM to a central plant with independent heat pump, evaporative cooling, and thermal energy storage components. This paper describes the technologies in question, presents modeling at the component and system levels, and demonstrates building energy and demand charge cost savings with local time-of-use tariffs in a hot-dry climate.

Helmns, Dre↗
DOE OSTI2021

Modeling Framework for Bulk Electric Grid Impacts from HEMP E1 and E3 Effects (Tasks 3.1 Final Report)

This report presents a framework to evaluate the impact of a high-altitude electromagnetic pulse (HEMP) event on a bulk electric power grid. This report limits itself to modeling the impact of EMP E1 and E3 components. The co-simulation of E1 and E3 is presented in detail, and the focus of the paper is on the framework rather than actual results. This approach is highly conservative as E1 and E3 are not maximized with the same event characteristics and may only slightly overlap. The actual results shown in this report are based on a synthetic grid with synthetic data and a limited exemplary EMP model. The framework presented can be leveraged and used to analyze the impact of other threat scenarios, both manmade and natural disasters. This report d escribes a Monte-Carlo based methodology to probabilistically quantify the transient response of the power grid to a HEMP event. The approach uses multiple fundamental steps to characterize the system response to HEMP events, focused on the E1 and E3 components of the event. 1) Obtain component failure data related to HEMP events testing of components and creating component failure models. Use the component failure model to create component failure conditional probability density function (PDF) that is a function of the HEMP induced terminal voltage. 2) Model HEMP scenarios and calculate the E1 coupled voltage profiles seen by all system components. Model the same HEMP scenarios and calculate the transformer reactive power consumption profiles due to E3. 3) Sample each component failure PDF to determine which grid components will fail, due to the E1 voltage spike, for each scenario. 4) Perform dynamic simulations that incorporate the predicted component failures from E1 and reactive power consumption at each transformer affected by E3. These simulations allow for secondary transients to affect the relays/protection remaining in service which can lead to cascading outages. 5) Identify the locations and amount of load lost for each scenario through grid dynamic simulation. This can be an indication of the immediate grid impacts from a HEMP event. In addition, perform more detailed analysis to determine critical nodes and system trends. 6) To help realize the longer-term impacts, a security constrained alternating current optimal power flow (ACOPF) is run to maximize critical load served. This report describes a modeling framework to assess the systemic grid impacts due to a HEMP event. This stochastic simulation framework generates a large amount of data for each Monte Carlo replication, including HEMP location and characteristics, relay and component failures, E3 GIC profiles, cascading dynamics including voltage and frequency over time, and final system state. This data can then be analyzed to identify trends, e.g., unique system behavior modes or critical components whose failure is more likely to cause serious systemic effects. The proposed analysis process is demonstrated on a representative system. In order to draw realistic conclusions of the impact of a HEMP event on the grid, a significant amount of work remains with respect to modeling the impact on various grid components.

24 POWER TRANSMISSION AND DISTRIBUTION↗