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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

An Inrush Current Limit Method for SiC-based Multi-level Grid-connected Converter During Low-Voltage Ride-Through

Low-voltage ride-through is one of the grid requirements for grid-connected converters, and during the transient of grid voltage change, a large inrush current could be induced due to the control delay of the converter. This inrush current is more severe in SiC-based multi-level converters due to the small filter inductance, which is selected based on harmonics consideration. In this paper, a PWM mask method is used to limit the inrush current of a cascaded H-bridge-based grid-connected converter during the transient of grid voltage change. Simulation and experiment results are provided to validate this method.

Li, Haiguo↗
DOE OSTI2022

A Hierarchical Local Electricity Market for a DER-rich Grid Edge

We report with increasing penetration of distributed energy resources (DER) in the distribution system, it is critical to design market structures that enable smooth integration of DERs. A hierarchical local electricity market (LEM) structure is proposed in this paper with a secondary market (SM) at the lower level representing secondary feeders and a primary market (PM) at the upper level, representing primary feeders, in order to effectively use DERs to increase grid efficiency and resilience. The lower level SM enforces budget, power balance, and flexibility constraints and accounts for costs related to consumers, such as their disutility, flexibility limits, and commitment reliability, while the upper level PM enforces grid physics constraints such as power balance and capacity limits, and also minimizes line losses. The hierarchical LEM is extensively evaluated using a modified IEEE-123 bus with high DER penetration, with each primary feeder consisting of at least three secondary feeders. Data from a GridLAB-D model is used to emulate realistic power injections and load profiles over the course of 24 hours. The performance of the LEM is illustrated by delineating the family of power-injection profiles across the primary and secondary feeders as well as corresponding local electricity tariffs that vary across the distribution grid. Through numerical simulations, the hierarchical LEM is shown to improve the efficiency of the market in terms of lowering overall costs, including both the distribution-level locational marginal prices (d-LMP) as well as retail tariffs paid by customers. Together, it represents an overall framework for a Distribution System Operator (DSO) who can provide the oversight for the entire LEM.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2018

New York City Impacts on a Regional Heat Wave

Abstract Extreme heat events are projected to increase in magnitude and frequency throughout this century due to increasing global temperatures, making it critically important to acquire improved understanding of their genesis and interactions with large cities. This study presents an application of the factor separation method to assess combined impacts of a synoptic scale heat wave, urban land cover, and urban energy and momentum fluxes on temperatures and winds over New York City via use of high resolution simulations (1 km grid spacing) with an urbanized WRF model. Results showed that, while the heat wave had the largest contribution to temperatures (> 8°C), urban surface factors matched it in highly urbanized areas. Surface factors matched this in highly urbanized areas during night and early morning hours, with contributions up to 5°C, when calm land breeze conditions result in a strong urban heat island. Positive interactions between all factors during morning and nighttime indicate urban heat island amplification of up to 4°C during the heat wave. Midtown Manhattan vertical cross-sections, where urban canopies are most dense, showed a change in the sign (from positive to negative) of the contribution of the urban fluxes between night and day below 500 m, possibly due to radiation blocking and increased thermal storage by buildings as well as frictional effects opposing the incoming warm air.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Resolved Gravity Waves in the Extratropical Stratosphere: Effect of Horizontal Resolution Increase from O (10) to O (1) km

Global ECMWF IFS simulations with horizontal grid spacings of 1, 4, and 9 km are used to assess gravity wave forcing (GWF) in the extratropical stratosphere. Results with important implications for GWF parameterizations at high and intermediate resolutions are presented. A doubling in the zonal-mean resolved GWF is observed when the horizontal resolution is increased from 9 to 1 km. Small-scale gravity waves with horizontal wavelengths < 100 km dominate this increase. Over most regions, excluding the polar night jet in the Antarctic spring, the total (resolved + parameterized) GWF at 9 km (4 km) is underestimated by up to 30% (15%). This implies that the parameterization of GWF is still required at 9 and 4 km horizontal resolutions. Despite the small land area in the Southern Hemisphere (SH), the resolved orographic and nonorographic GWF contribute equally to the total GWF in the SH at 1 km resolution. This is not reflected in the partitioning of the parameterized GWF, which has a significantly larger nonorographic contribution at 9 km. As a result, a zonal-mean momentum budget analysis reveals that the total GWF contributes one-third of SH springtime polar vortex deceleration at 1 km, whereas the contribution is as much as 50% at 9 km. In conclusion, this suggests that a rebalancing of the parameterized nonorographic and orographic GWF is required.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Changes in Temperature and Precipitation Extremes in Superparameterized CAM in Response to Warmer SSTs

Subdaily temperature and precipitation extremes in response to warmer SSTs are investigated on a global scale using the superparameterized (SP) Community Atmosphere Model (CAM), in which a cloud-resolving model is embedded in each CAM grid column to simulate convection explicitly. Two 10-yr simulations have been performed using present climatological sea surface temperature (SST) and perturbed SST climatology derived from the representative concentration pathway 8.5 (RCP8.5) scenario. Compared with the conventional CAM, SP-CAM simulates colder temperatures and more realistic intensity distribution of precipitation, especially for heavy precipitation. The temperature and precipitation extremes have been defined by the 99th percentile of the 3-hourly data. For temperature, the changes in the warm and cold extremes are generally consistent between CAM and SP-CAM, with larger changes in warm extremes at low latitudes and larger changes in cold extremes at mid-to-high latitudes. For precipitation, CAM predicts a uniform increase of frequency of precipitation extremes regardless of the rain rate, while SP-CAM predicts a monotonic increase of frequency with increasing rain rate and larger change of intensity for heavier precipitation. The changes in 3-hourly and daily temperature extremes are found to be similar; however, the 3-hourly precipitation extremes have a significantly larger change than daily extremes. The Clausius–Clapeyron scaling is found to be a relatively good predictor of zonally averaged changes in precipitation extremes over midlatitudes but not as good over the tropics and subtropics. The changes in precipitable water and large-scale vertical velocity are equally important to explain the changes in precipitation extremes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

On Artificial Viscosity for Shocks in Porous Materials

The treatment of downstream oscillations in shockwave simulations on discrete grids dates back to the 1940’s with the introduction of an artificial viscosity. These additional term(s) have been successful at suppressing or eliminating oscillations in smooth, well-defined equations of state. This report identifies a shortcoming of several solution methodologies and viscosity formulations that are easily observed when solid equations of state are modified with compaction models or run with a significantly lower density (i.e. ρ < ρ 0 ). The obvious solution is to increase parameter values for the artificial viscosity, but this could lead to deleterious effects in multi-material problems if the viscosity is globally set.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Evaluating the Feasibility of MCNP/FDTD EMP Calculations at Altitudes above 20km

This report investigates the use of MCNP/FDTD for 2D EMP calculations at altitudes above 20km. The objective is to determine the feasibility of performing EMP calculations for high altitude nuclear detonations using the MCNP/FDTD codes. The work in this report indicates that the issue with running such EMP calculations lies more in MCNP than in FDTD. In general, when performing EMP calculations with MCNP/FDTD, the majority of the computation time is in getting the results from MCNP. Typically, for EMP calculations with spatial domains of about 5 to 10km in both the x and z directions, the MCNP calculations take weeks to complete and the FDTD portion of the computation only takes one to two days. The problems with running MCNP on large problems are largely due to memory issues. One issue is that MCNP requires all tally grids of a simulation to reside in the memory of a single node (128GB on snow) and this, in turn, requires multiple MCNP simulations to complete large spatial/time domain problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Colorado Residential Retrofit Energy District (CoRRED) Phase I: Final Modeling Results

Electrification of buildings and transportation coupled with increased deployment of distributed energy resources (DERs) has been identified as a key step toward meeting emissions reduction goals across the U.S. Examples of pilot projects that feature innovative, 'smart' electric neighborhoods are largely focused on new construction projects, but there is a need to address the millions of existing homes so that they too may accommodate cleaner yet variable energy production in ways that benefit both the utility grid and the homes' residents. The Colorado Residential Retrofit Energy District (CoRRED) project used building and grid co-simulation tools to model an advanced energy district demonstration in an existing residential neighborhood in Denver, Colorado, and explored how existing building and utility infrastructures can be enhanced with combinations of traditional energy-efficiency retrofit measures and integration of solar panels and other DER technologies to provide better affordability and reliability. We analyzed which packages of DERs most reliably enable demand flexibility in response to a time-of-use (TOU) rate. Our results clearly indicate that incorporating DERs into efficient electrification produces much bigger utility bill savings on an annual basis than efficient electrification without DERs. While electrifying a neighborhood will increase the maximum load, batteries and solar panels can reduce load during peak hours so that the community can be a net producer during peak periods. A remaining challenge is to overcome first cost barriers to implementing energy-efficiency and DER technologies, as modest utility bill savings require long payback periods that are not practical for most homeowners.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

CRADA Number NFE-24-10110 with Qubit Engineering Inc. (CRADA Final Report)

Over the past year, the Qubit Engineering team has pushed the frontiers of power‑grid optimization, working in close collaboration with Oak Ridge National Laboratory (ORNL) and the Tennessee Valley Authority (TVA). Their progress is reflected in three newly submitted conference papers, “Unified Relational GNN Architecture for AC Optimal Power Flow Calculations in Electric Grids,” “Graph‑Based Attention Mechanisms for Solving the AC Optimal Power Flow Problem in Electrical‑Power Networks,” and “Enhanced Power‑Grid Maintenance Planning and Quantum‑Inspired Combinatorial Prospects.” These publications showcase state‑of‑the‑art graph‑neural‑network methods for AC‑OPF and novel quantum‑inspired heuristics for maintenance scheduling. Beyond the academic results, the Qubit team has converted the research into two production‑grade tools built on TVA data: Neuro‑Grid, an AI‑driven power‑flow simulator that provides instant, interactive full‑grid load‑flow visualizations, and Quanta‑Grid, a quantum‑inspired maintenance‑scheduling engine to support logistics optimization for power utilities. Together, these advances demonstrate how Qubit’s partnership with ORNL and TVA is delivering practical, physics‑grounded analytics for next‑generation grid management.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Distributed Transient Safety Verification via Robust Control Invariant Sets: A Microgrid Application

Modern safety-critical energy infrastructures are increasingly operated in a hierarchical and modular control framework which allows for limited data exchange between the modules. In this context, it is important for each module to synthesize and communicate constraints on the values of exchanged information in order to assure system-wide safety. To ensure transient safety in inverter-based microgrids, we develop a set invariance-based distributed safety verification algorithm for each inverter module. Applying Nagumo's invariance condition, we construct a robust polynomial optimization problem to jointly search for safety-admissible set of control set-points and design parameters, under allowable disturbances from neighbors. We use sum-of-squares (SOS) programming to solve the verification problem and we perform numerical simulations using grid-forming inverters to illustrate the algorithm.

Bouvier, Jean-Baptiste H.↗
DOE OSTI2024

Cholla-MHD: An Exascale-capable Magnetohydrodynamic Extension to the Cholla Astrophysical Simulation Code

Abstract We present an extension of the massively parallel, GPU native, astrophysical hydrodynamics code Cholla to magnetohydrodynamics (MHD). Cholla solves the ideal MHD equations in their Eulerian form on a static Cartesian mesh utilizing the Van Leer + constrained transport integrator, the HLLD Riemann solver, and reconstruction methods at second and third order. Cholla’s MHD module can perform ≈260 million cell updates per GPU-second on an NVIDIA A100 while using the HLLD Riemann solver and second order reconstruction. The inherently parallel nature of GPUs combined with increased memory in new hardware allows Cholla’s MHD module to perform simulations with resolutions ∼500 3 cells on a single high-end GPU (e.g., an NVIDIA A100 with 80 GB of memory). We employ GPU direct Message Passing Interface to attain excellent weak scaling on the exascale supercomputer Frontier, while using 74,088 GPUs and simulating a total grid size of over 7.2 trillion cells. A suite of test problems highlights the accuracy of Cholla’s MHD module and demonstrates that zero magnetic divergence in solutions is maintained to round off error. We also present new testing and CI tools using GoogleTest, GitHub Actions, and Jenkins that have made development more robust and accurate and ensure reliability in the future.

Astronomy & Astrophysics↗
DOE OSTI2021

Cloud droplet diffusional growth in homogeneous isotropic turbulence: bin microphysics versus Lagrangian super-droplet simulations

The increase in the spectral width of an initially monodisperse population of cloud droplets in homogeneous isotropic turbulence is investigated by applying a finite-difference fluid flow model combined with either Eulerian bin microphysics or a Lagrangian particle-based scheme. The turbulence is forced applying a variant of the so-called linear forcing method that maintains the mean turbulent kinetic energy (TKE) and the TKE partitioning between velocity components. The latter is important for maintaining the quasi-steady forcing of the supersaturation fluctuations that drive the increase in the spectral width. We apply a large computational domain (64 3 m 3 ), one of the domains considered in Thomas et al. (2020). The simulations apply 1 m grid length and are in the spirit of the implicit large eddy simulation (ILES), that is, with small-scale dissipation provided by the model numerics. This is in contrast to the scaled-up direct numerical simulation (DNS) applied in Thomas et al. (2020). Two TKE intensities and three different droplet concentrations are considered. Analytic solutions derived in Sardina et al. (2015), valid for the case when the turbulence integral timescale is much larger than the droplet phase relaxation timescale, are used to guide the comparison between the two microphysics simulation techniques. The Lagrangian approach reproduces the scalings relatively well. Representing the spectral width increase in time is more challenging for the bin microphysics because appropriately high resolution in the bin space is needed. The bin width of 0.5 µm is only sufficient for the lowest droplet concentration (26 cm -3 ). For the highest droplet concentration (650 cm -3 ), an order of magnitude smaller bin size is barely sufficient. The scalings are not expected to be valid for the lowest droplet concentration and the high-TKE case, and the two microphysics schemes represent similar departures. Finally, because the fluid flow is the same for all simulations featuring either low or high TKE, one can compare point-by-point simulation results. Such a comparison shows very close temperature and water vapor point-by-point values across the computational domain and larger differences between simulated mean droplet radii and spectral width. The latter are explained by fundamental differences in the two simulation methodologies, numerical diffusion in the Eulerian bin approach and a relatively small number of Lagrangian particles that are used in the particle-based microphysics.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1974

System support software for the Space Ultrareliable Modular Computer (SUMC)

The highly transportable programming system designed and implemented to support the development of software for the Space Ultrareliable Modular Computer (SUMC) is described. The SUMC system support software consists of program modules called processors. The initial set of processors consists of the supervisor, the general purpose assembler for SUMC instruction and microcode input, linkage editors, an instruction level simulator, a microcode grid print processor, and user oriented utility programs. A FORTRAN 4 compiler is undergoing development. The design facilitates the addition of new processors with a minimum effort and provides the user quasi host independence on the ground based operational software development computer. Additional capability is provided to accommodate variations in the SUMC architecture without consequent major modifications in the initial processors.

Hill, T. E.↗
NASA NTRS1985

A simulation of rotor-stator interaction using the Euler equations and patched grids

An unsteady Euler code to study rotor-stator interaction problem was developed. The code uses patched grids that move relative to each other to simulate the motion of the rotor airfoils with respect to the stator airfoils. The Osher integration scheme is used in conjunction with an implicit relaxation approach. The scheme is second order accurate in space and time, and is also TVD in each spatial direction. The numerical results were found to be periodic in time, thus demonstrating the capability of the integration and zonal schemes in simulating periodic time dependent flow. The pressure contours obtained are almost oscillation free because of the TVD nature of the scheme. A new procedure was developed to simulate flows about bodies that move relative to each other. This capability should prove to be very useful in the areas of rotor-stator interaction, propeller-nacelle interaction, and helicopter rotor-fuselage interaction.

Rai, M. M.↗
NASA NTRS1992

CFD applications in pump flows

The objective of the proposed paper is to develop a computational procedure that solves incompressible Navier-Stokes equations for pump flows. The solution method is based on the pseudo-compressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. The equations are solved in steadily rotating reference frames and the centrifugal force and the Coriolis force are added to the equation of motion. As a benchmark problem, the flow through the Rocketdyne inducer is numerically simulated. A coarse grid solution is obtained with a single zone by using an algebraic turbulence model. In multi-zone fine grid computation, a one-equation Baldwin-Barth turbulence model is utilized. Numerical results are compared with experimental measurements and a good agreement is found between the two. The resulting computer code is then applied to the flow analysis inside a two-stage fuel pump impeller operating at 80 percent, 100 percent, and 120 percent of design flow.

Kiris, Cetin↗
NASA NTRS1993

Particle kinetic simulation of high altitude hypervelocity flight

In this grant period, the focus has been on enhancement and application of the direct simulation Monte Carlo (DSMC) particle method for computing hypersonic flows of re-entry vehicles. Enhancement efforts dealt with modeling gas-gas interactions for thermal non-equilibrium relaxation processes and gas-surface interactions for prediction of vehicle surface temperatures. Both are important for application to problems of engineering interest. The code was employed in a parametric study to improve future applications, and in simulations of aeropass maneuvers in support of the Magellan mission. Detailed comparisons between continuum models for internal energy relaxation and DSMC models reveals that several discrepancies exist. These include definitions of relaxation parameters and the methodologies for implementing them in DSMC codes. These issues were clarified and all differences were rectified in a paper (Appendix A) submitted to Physics of Fluids A, featuring several key figures in the DSMC community as co-authors and B. Haas as first author. This material will be presented at the Fluid Dynamics meeting of the American Physical Society on November 21, 1993. The aerodynamics of space vehicles in highly rarefied flows are very sensitive to the vehicle surface temperatures. Rather than require prescribed temperature estimates for spacecraft as is typically done in DSMC methods, a new technique was developed which couples the dynamic surface heat transfer characteristics into the DSMC flow simulation code to compute surface temperatures directly. This model, when applied to thin planar bodies such as solar panels, was described in AIAA Paper No. 93-2765 (Appendix B) and was presented at the Thermophysics Conference in July 1993. The paper has been submitted to the Journal of Thermophysics and Heat Transfer. Application of the DSMC method to problems of practical interest requires a trade off between solution accuracy and computational expense and limitations. A parametric study was performed and reported in AIAA Paper No. 93-2806 (Appendix C) which assessed the accuracy penalties associated with simulations of varying grid resolution and flow domain size. The paper was also presented at the Thermophysics Conference and will be submitted to the journal shortly. Finally, the DSMC code was employed to assess the pitch, yaw, and roll aerodynamics of the Magellan spacecraft during entry into the Venus atmosphere at off-design attitudes. This work was in support of the Magellan aerobraking maneuver of May 25-Aug. 3, 1993. Furthermore, analysis of the roll characteristics of the configuration with canted solar panels was performed in support of the proposed 'Windmill' experiment. Results were reported in AIAA Paper No. 93-3676 (Appendix D) presented at the Atmospheric Flight Mechanics Conference in August 1993, and were submitted to Journal of Spacecraft and Rockets.

Haas, Brian L.↗
NASA NTRS1993

A parallel adaptive mesh refinement algorithm

Over recent years, Adaptive Mesh Refinement (AMR) algorithms which dynamically match the local resolution of the computational grid to the numerical solution being sought have emerged as powerful tools for solving problems that contain disparate length and time scales. In particular, several workers have demonstrated the effectiveness of employing an adaptive, block-structured hierarchical grid system for simulations of complex shock wave phenomena. Unfortunately, from the parallel algorithm developer's viewpoint, this class of scheme is quite involved; these schemes cannot be distilled down to a small kernel upon which various parallelizing strategies may be tested. However, because of their block-structured nature such schemes are inherently parallel, so all is not lost. In this paper we describe the method by which Quirk's AMR algorithm has been parallelized. This method is built upon just a few simple message passing routines and so it may be implemented across a broad class of MIMD machines. Moreover, the method of parallelization is such that the original serial code is left virtually intact, and so we are left with just a single product to support. The importance of this fact should not be underestimated given the size and complexity of the original algorithm.

Quirk, James J.↗
NASA NTRS2000

Cosmic Rays and Their Radiative Processes in Numerical Cosmology

A cosmological hydrodynamic code is described, which includes a routine to compute cosmic ray acceleration and transport in a simplified way. The routine was designed to follow explicitly diffusive, acceleration at shocks, and second-order Fermi acceleration and adiabatic loss in smooth flows. Synchrotron cooling of the electron population can also be followed. The updated code is intended to be used to study the properties of nonthermal synchrotron emission and inverse Compton scattering from electron cosmic rays in clusters of galaxies, in addition to the properties of thermal bremsstrahlung emission from hot gas. The results of a test simulation using a grid of 128 (exp 3) cells are presented, where cosmic rays and magnetic field have been treated passively and synchrotron cooling of cosmic ray electrons has not been included.

Ryu, Dongsu↗