Search NASA⌕ Search

SEARCH · Search NASA

Results for “Grid 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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Photovoltaic and Behind-the-Meter Battery Storage: Advanced Smart Inverter Controls and Field Demonstration

Electric utilities have little visibility of the electrical distribution system, and consequently, limited diagnostic capabilities. The distribution grid was designed for a unidirectional power flow, where energy is supplied by few large centralized power plants; however, this is set to change to meet California’s aggressive decarbonization goals. The large-scale deployment of distributed renewable generation, such as photovoltaics (PV), can have many negative effects on an unprepared grid. As illustrated in the California “duck curve,” PV generation modifies load profiles during the day, which causes steep ramping in the evening. This effect indicates the need to revise the electrical grid’s design and operation. This project sought to (a) create a centralized resource in California to test and validate distribution technology controls with industry standard PV, storage, and high-fidelity sensors, (b) support strategic and operational decisions for new grid architectures by providing simulation models, (c) promote new ways to control clusters of PV smart inverters, in accordance with California Rule 21, and (d) innovate, develop, and field test a predictive controller to maximize economic profit for the customer while supporting the grid. The controller was built using the state-of-the-art model predictive control methodology to optimally control behind-the-meter PV and battery storage. In consideration of the duck curve, the controller optimally controls the battery dispatch by charging during excess generation periods and discharging during the critical afternoon demand hours. The controller was evaluated in annual simulations and revealed the potential cost-effectiveness of behind-the-meter battery storage. The simulations showed that the annual electricity bill could be reduced by as much as 35 percent, with a payback period of the investment in battery storage in about 6 years – significantly shorter than the manufacturer’s 10-year warranty. All developed simulation models, the grid event library, and the MPC controller are open-source and available online.

Model Predictive control↗
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 OSTI2017

Assessment of Vertical Mesh Refinement in Concurrently Nested Large-Eddy Simulations Using the Weather Research and Forecasting Model

To facilitate multiscale simulation using the Weather Research and Forecasting Model, vertical mesh refinement for one-way concurrent nested simulation was recently introduced. Grid refinement in the vertical dimension removes the requirement of different grid aspect ratios on the bounding versus the nested domain, such that results from refinement are in the horizontal directions only, and thereby can also reduce numerical errors on the bounding domain arising from large aspect ratios in the presence of complex terrain. Herein, the impacts of vertical grid refinement on the evolving downstream flow in nested large-eddy simulations are evaluated in relation to other model configuration choices, including turbulence subfilter-scale (SFS) stress models, mesh configuration, and alternative methods for calculating several near-surface flow parameters. Although vertical nesting requires coarsening of the vertical grid on the bounding domain, leading to a smaller range of resolved turbulence scales in the nest’s lateral boundary conditions, parameter values within the nested domains are generally only minimally impacted, relative to nesting using the same vertical grid on each domain. Two dynamic SFS models examined herein generally improved the simulated mean wind speed, turbulence kinetic energy, stresses and spectra, on both domains, and accelerated equilibration rates within nested domains, relative to two constant coefficient models. A new method of extrapolating horizontal velocity components to near-surface locations at nested domain lateral boundaries, and a correction to the calculation of deformation elements near the surface, are each shown to slightly alter the mean parameter values, yet only minimally impact equilibration rates within the nested domain.

14 SOLAR ENERGY↗
DOE OSTI2012

CENTRAL REGIONS OF BARRED GALAXIES: TWO-DIMENSIONAL NON-SELF-GRAVITATING HYDRODYNAMIC SIMULATIONS

The inner regions of barred galaxies contain substructures such as off-axis shocks, nuclear rings, and nuclear spirals. These substructures may affect star formation, and control the activity of a central black hole (BH) by determining the mass inflow rate. We investigate the formation and properties of such substructures using high-resolution, grid-based hydrodynamic simulations. The gaseous medium is assumed to be infinitesimally thin, isothermal, and non-self-gravitating. The stars and dark matter are represented by a static gravitational potential with four components: a stellar disk, a bulge, a central BH, and a bar. To investigate various galactic environments, we vary the gas sound speed, c{sub s} , as well as the mass of the central BH, M{sub BH}. Once the flow has reached a quasi-steady state, off-axis shocks tend to move closer to the bar major axis as c{sub s} increases. Nuclear rings shrink in size with increasing c{sub s} , but are independent of M{sub BH}, suggesting that the ring position is not determined by the Lindblad resonances. Rings in low-c{sub s} models are narrow since they are occupied largely by gas on x{sub 2}-orbits and well decoupled from nuclear spirals, while they become broad because of large thermal perturbations in high-c{sub s} models. Nuclear spirals persist only when either c{sub s} is small or M{sub BH} is large; they would otherwise be destroyed completely by the ring material on eccentric orbits. The shape and strength of nuclear spirals depend sensitively on c{sub s} and M{sub BH} such that they are leading if both c{sub s} and M{sub BH} are small, weak trailing if c{sub s} is small and M{sub BH} is large, and strong trailing if both c{sub s} and M{sub BH} are large. While the mass inflow rate toward the nucleus is quite small in low-c{sub s} models because of the presence of a narrow nuclear ring, it becomes larger than 0.01 M{sub Sun} yr{sup -1} when c{sub s} is large, providing a potential explanation of nuclear activity in Seyfert galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Phase transitions of macromolecular microsphere composite hydrogels based on the stochastic Cahn–Hilliard equation

We use the stochastic Cahn–Hilliard equation to simulate the phase transitions of the macromolecular microsphere composite (MMC) hydrogels under a random disturbance. Based on the Flory–Huggins lattice model and the Boltzmann entropy theorem, we develop a reticular free energy suit for the network structure of MMC hydrogels. Taking the random factor into account, with the time-dependent Ginzburg-Landau (TDGL) mesoscopic simulation method, we set up a stochastic Cahn–Hilliard equation, designated herein as the MMC-TDGL equation. The stochastic term in the equation is constructed appropriately to satisfy the fluctuation-dissipation theorem and is discretized on a spatial grid for the simulation. A semi-implicit difference scheme is adopted to numerically solve the MMC-TDGL equation. Some numerical experiments are performed with different parameters. The results are consistent with the physical phenomenon, which verifies the good simulation of the stochastic term.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Incorporation of aerosol into the COSPv2 satellite lidar simulator for climate model evaluation

Abstract. Atmospheric aerosol has substantial impacts on climate, air quality and biogeochemical cycles, and its concentrations are highly variable in space and time. A key variability to evaluate within models that simulate aerosol is the vertical distribution, which influences atmospheric heating profiles and aerosol–cloud interactions, to help constrain aerosol residence time and to better represent the magnitude of simulated impacts. To ensure a consistent comparison between modeled and observed vertical distribution of aerosol, we implemented an aerosol lidar simulator within the Cloud Feedback Model Intercomparison Project (CFMIP) Observation Simulator Package version 2 (COSPv2). We assessed the attenuated total backscattered (ATB) signal and the backscatter ratios (SRs) at 532 nm in the U.S. Department of Energy's Energy Exascale Earth System Model version 1 (E3SMv1). The simulator performs the computations at the same vertical resolution as the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP), making use of aerosol optics from the E3SMv1 model as inputs and assuming that aerosol is uniformly distributed horizontally within each model grid box. The simulator applies a cloud masking and an aerosol detection threshold to obtain the ATB and SR profiles that would be observed above clouds by CALIOP with its aerosol detection capability. Our analysis shows that the aerosol distribution simulated at a seasonal timescale is generally in good agreement with observations. Over the Southern Ocean, however, the model does not produce the SR maximum as observed in the real world. Comparison between clear-sky and all-sky SRs shows little differences, indicating that the cloud screening by potentially incorrect model clouds does not affect the mean aerosol signal averaged over a season. This indicates that the differences between observed and simulated SR values are due not to sampling errors, but to deficiencies in the representation of aerosol in models. Finally, we highlight the need for future applications of lidar observations at multiple wavelengths to provide insights into aerosol properties and distribution and their representation in Earth system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Methods for estimating X-ray machine output through measurement and simulation

We report ball grid arrays are increasingly being applied in the electronics industry and may require X-ray inspection to ensure the integrity and correct placement of solder pins. However, as the architecture of integrated circuits continues to narrow while simultaneously growing more complex, the risk of electronic failure due to radiation damage increases. While medical X-ray devices have been held to high standards and are repeatedly shown to be well characterized, devices used for electronic inspection are often lacking detailed characterization. This study presents unique methods to solve for important properties in X-ray inspection devices such as source to object distance and energy spectrum. This information can then be applied to Monte Carlo models to achieve better overall dose estimates to electronics, which will lead to superior manufactured products. Since X-ray devices can vary greatly in source characteristics, this work investigates spectral measurement and Monte Carlo representation of three X-ray devices. For a Philips SRO 33 100 medical diagnostic device, the spectral output followed expected trends given by the prediction software SpekCalc and Spektr. For the Dage XD7500NT, direct measurement showed a spectral artifact that through the use of Gafchromic films, was shown to be a contributing effect in the dose output. For the Rad Source RS1800, a high powered irradiation device, direct spectral measurement was not achieved. However, a Monte Carlo model using an assumed spectra was found to match ion chamber measurements to a high degree.

47 OTHER INSTRUMENTATION↗
DOE OSTI2012

Calculating frequency at loads in simulations of electro-mechanical transients

This paper introduces a new method for calculating frequency at an electrical load in simulations of electro-mechanical transients. The method is designed for simulation studies that require accurate models of sensors such as phasor measurement units and F-Net devices, which measure frequency at locations away from generating plants. These sensors are poised to become critical components in the control systems of electrical power grids, and therefore simulation tools that incorporate accurate models of these devices are essential. The method proposed here corrects two drawbacks of using numerically computed phase angle derivatives to approximate frequency. First, it eliminates spurious spikes in frequencies calculated at loads. Second, it eliminates instabilities induced by the simulator in studies of frequency responsive loads. The proposed method is derived from a simplified model of the generators and loads in an electrical system, but in the final analysis does not depend critically on these simplifications and is therefore applicable to more sophisticated models. The method is demonstrated with the simplified model applied to the IEEE 14 and 300 bus systems.

42 ENGINEERING↗
DOE OSTI2018

On the metal/ZnO contacts in a sliding-bending piezoelectric nanogenerator

We carry out molecular dynamics (MD) simulations and analyze the piezoelectric and mechanical properties of a ZnO nanowire (NW) under mechanical bending in a sliding-bending piezoelectric nanogenerator. The computational model contains a vertically aligned ZnO NW along the [0001] direction and a Pt (111) metal tip. MD simulations are performed to mimic the interfacial sliding dynamics of an atomic force microscope (AFM) metal tip sliding over a ZnO NW. Ionic charges for Zn and O in the ZnO NW are used to calculate the piezoelectric potential distributions in a grid-point matrix. Simulation results reveal the very detailed distributions of the piezoelectric potential within the ZnO NW, which are different from the continuum field theory predictions. The very similar variations of the lateral friction and contact forces versus the sliding distance of the AFM metal tip over the ZnO NW during repeated forward-backward sliding suggest that the metal/ZnO contact is robust. Moreover, the substantial contact area at the metal/ZnO interface suggests that a Schottky contact should be established at the interface, which needs further investigations.

Chemistry↗
DOE OSTI2014

On the diversity of magnetic interactions in close-in star-planet systems

Magnetic interactions between close-in planets and their host star can play an important role in the secular orbital evolution of the planets, as well as the rotational evolution of their host. As long as the planet orbits inside the Alfvén surface of the stellar wind, the magnetic interaction between the star and the planet can modify the wind properties and also lead to direct angular momentum transfers between the two. We model these star-planet interactions using compressible magnetohydrodynamic (MHD) simulations, and quantify the angular momentum transfers between the star, the planet, and the stellar wind. We study the cases of magnetized and non-magnetized planets and vary the orbital radius inside the Alfvén surface of the stellar wind. Based on a grid of numerical simulations, we propose general scaling laws for the modification of the stellar wind torque, for the torque between the star and the planet, and for the planet migration associated with the star-planet magnetic interactions. We show that when the coronal magnetic field is large enough and the star is rotating sufficiently slowly, the effect of the magnetic star-planet interaction is comparable to tidal effects and can lead to a rapid orbital decay.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Amplitude flux, probability flux, and gauge invariance in the finite volume scheme for the Schrödinger equation

The time-dependent Schrödinger equation can be put in a probability conserving, gauge invariant form, on arbitrary structured grids via finite volume discretization. The gauge terms in the discrete system cancel with a portion of the amplitude flux to produce abbreviated flux functions. The resulting time translation operator is strictly unitary, and is compatible with an efficient operator splitting scheme that allows for multi-dimensional simulation with complex grid geometries. Moreover, the abbreviated amplitude flux is necessary to the construction of a conservative probability current. This construction turns out to be important when computing Bohmian trajectories in multi-dimensions. Bohmian trajectories are useful in the interpretation of quantum mechanical phenomena such as tunneling ionization, and provide a bridge between quantum and classical regimes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Second release of the CoRe database of binary neutron star merger waveforms

Abstract We present the second data release of gravitational waveforms from binary neutron star (BNS) merger simulations performed by the Computational Relativity ( CoRe ) collaboration. The current database consists of 254 different BNS configurations and a total of 590 individual numerical-relativity simulations using various grid resolutions. The released waveform data contain the strain and the Weyl curvature multipoles up to ℓ = m = 4 . They span a significant portion of the mass, mass-ratio, spin and eccentricity parameter space and include targeted configurations to the events GW170817 and GW190425. CoRe simulations are performed with 18 different equations of state, seven of which are finite temperature models, and three of which account for non-hadronic degrees of freedom. About half of the released data are computed with high-order hydrodynamics schemes for tens of orbits to merger; the other half is computed with advanced microphysics. We showcase a standard waveform error analysis and discuss the accuracy of the database in terms of faithfulness. We present ready-to-use fitting formulas for equation of state-insensitive relations at merger (e.g. merger frequency), luminosity peak, and post-merger spectrum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Battery Performance, Thermal, and Life Modeling for Southern California Edison (Cooperative Research and Development Final Report)

The objective of this project is to develop battery performance models, thermal models and life models for SCE to use in simulation and analysis of grid services involving energy storage. NREL will fit the models to three separate chemistries using test data provided by SCE. The dataset for each chemistry will include electrical, thermal and aging response to different temperatures and cycling conditions, measured under a variety of cell-level and module-level experiments described in SCE test protocol documents. Once tuned to a test dataset, the battery system model software will provide predictions of battery energy and power loss, and thus cycle and calendar life, for any energy storage grid service power profile of interest to SCE. The automated battery life modeling and simulation tool can also be applied to any other grid battery systems as long as necessary input data is available to accelerate battery model development and battery lifetime analysis.

25 ENERGY STORAGE↗
DOE OSTI2018

Advanced Inverter Voltage Controls: Simulation and Field Pilot Findings

This report describes work performed by the Hawaiian Electric Companies and the National Renewable Energy Laboratory (NREL) to model and simulate advanced inverter grid-support functions and to validate and expand on those simulations through a field pilot study. This work builds on earlier research, referred to as the Voltage Regulation Operational Strategies study. The objective of both is to investigate functionalities available in most photovoltaic (PV) inverters to modulate active and reactive power autonomously based on local voltage measurements for the purpose of mitigating off-nominal grid voltage conditions. Specifically of interest are volt/volt-ampere reactive (VAR) control and volt/watt control, the effect of those functions on quasi-steady-state feeder voltages, and the impact of the functions on PV energy production. Because volt/VAR and volt/watt control autonomously adjust inverter output based on local conditions without requiring communication with any other devices, they are good candidates to increase PV hosting capacity when it is limited by voltage constraints in a feeder with very large numbers of PV systems.

14 SOLAR ENERGY↗
DOE OSTI2026

High fidelity actuator line data from 9 turbine wind farm simulations using ExaWind

This data was generated with the ExaWind code suite (https://github.com/Exawind) to investigate the performance of different Active Wake Mixing turbine control in a wind farm situated in a stable atmospheric boundary layer. All cases correspond to a 3x3 wind farm in a 10km x 10km domain using a total mesh size that varied between 1.6 X 10^9 to 1.85 X 10^9 grid cells. The simulations were run across 1800-2000 GPUs on Frontier. The case description and data generation process is fully documented in Yalla, G. R., Brown, K., Cheung, L., Houck, D., deVelder, N., and Balaji, J. (2025). "Estimating annual energy production of wake mixing control strategies including comparisons to wake steering." Wind Energy Sciences (https://doi.org/10.5194/wes-2025-250).

17 WIND ENERGY↗