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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Identifying Potential Sub-Synchronous Oscillations Using Impedance Scan Approach: Preprint

This paper presents an impedance scan study on the real-world power system in Australia that has observed 17 to 20 Hz intermittent subsynchronous oscillations. Through the impedance scan of each of the IBRs, both individually and collectively, potential resonance modes are identified. The impedance scans were carried out using electromagnetic transient PSCAD models of the network. The network comprises of site-specific, black-boxed models of IBRs supplied by the generators. The impedance scan approach was divided in three major steps: (1) the first step is to identify IBRs where impedance analysis needs to be performed based on the magnitude of oscillations observed at their points of interconnection (POIs); (2) the second step focuses on performing impedance scans at selected IBRs in single-machine infinite-bus (SMIB) configuration to identify internal resonance modes of an IBR and to evaluate if any of these modes become unstable under certain grid conditions; (3) the third step performs impedance scans at a few selected IBRs while connected to the wider network model to obtain the impedance response of both the IBR and the grid - this step evaluates control interactions among IBRs. The impedance scan study found that a few IBRs have an underdamped resonance mode at around 17 Hz, which becomes unstable under a certain operating condition. Another interesting finding was that certain IBRs increase the effective grid impedance seen by another IBR in proximity for a particular operating condition. Under this operating condition, the resulting subsynchronous oscillation mode is more pronounced due to a combination of control interactions among IBRs through the transmission network under study and a resonance mode inside the IBRs.

control interactions↗
NASA NTRS2011

The Grid Density Dependence of the Unsteady Pressures of the J-2X Turbines

The J-2X engine was originally designed for the upper stage of the cancelled Crew Launch Vehicle. Although the Crew Launch Vehicle was cancelled the J-2X engine, which is currently undergoing hot-fire testing, may be used on future programs. The J-2X engine is a direct descendent of the J-2 engine which powered the upper stage during the Apollo program. Many changes including a thrust increase from 230K to 294K lbf have been implemented in this engine. As part of the design requirements, the turbine blades must meet minimum high cycle fatigue factors of safety for various vibrational modes that have resonant frequencies in the engine's operating range. The unsteady blade loading is calculated directly from CFD simulations. A grid density study was performed to understand the sensitivity of the spatial loading and the magnitude of the on blade loading due to changes in grid density. Given that the unsteady blade loading has a first order effect on the high cycle fatigue factors of safety, it is important to understand the level of convergence when applying the unsteady loads. The convergence of the unsteady pressures of several grid densities will be presented for various frequencies in the engine's operating range.

Schmauch, Preston B.↗
DOE OSTI2018

Two-step perfectly matched layer for arbitrary-order pseudo-spectral analytical time-domain methods

Numerical simulation of an electrodynamic system in empty space requires the implementation of open boundary conditions (BC) to terminate the solution of Maxwell’s equations on the boundaries of the computational domain. The Perfectly Matched Layer (PML) has become the method of choice for open BC with wave equations, as it is straightforward and relatively easy to implement, and offers very efficient and user-adjustable absorption rates. PMLs are most often employed with the Finite-Difference Time-Domain (FDTD) algorithm, which in its most common implementation offers second-order accuracy in space and time on Cartesian grids. Yet, simulations (including some class of electromagnetic Particle-In-Cell simulations) that require higher precision may resort to higher-order Maxwell solvers employing extended finite-difference stencils, or even to pseudo-spectral Maxwell solvers, for which a general, versatile and efficient formulation of the PML has been missing so far. Here we propose a novel “two-step” formulation of the PML that is simple, very versatile and can be used as is with any Maxwell solver. In particular, it is applicable to a large class of Maxwell solvers including the arbitrary-order Pseudo-Spectral Analytical Time-Domain (PSATD) solver, which offers arbitrarily low numerical dispersion when increasing solver order and becomes dispersion-free at infinite order. Analysis and numerical simulations demonstrate that the new formulation is as efficient as the standard PML formulation, both for the FDTD and the PSATD implementations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

The Impacts of Rotational Mixing on the Precipitation Simulated by a Convection Permitting Model

With increased availability of computational resources, regional and global scale convection-permitting model (CPM, Δx ~ 1–10 km) simulations are becoming more common. CPMs have improved accuracy in their representation of deep convection and mesoscale convective systems (MCSs) compared to coarser resolution models. However, CPMs still exhibit convective cloud and precipitation biases relative to observations, notably a lesser frequency of light precipitation rates and greater frequency of heavy precipitation rates. In this work we hypothesize that these CPM biases are related to under-resolved mixing between convective updrafts and their surrounding environment. To test this hypothesis, we introduce a parameterization to the Weather Research and Forecasting model (WRF) that adds a small angular rotation of the grid-scale flow about the axis perpendicular to the plane of convective drafts. This rotated flow is then allowed to alter advection of moisture and hydrometeors. The effects of such mixing on precipitation characteristics are evaluated in month-long 4-km grid spacing simulations over the Amazon. The enhanced mixing transports moisture and condensate from convective cores to other areas including downdrafts. This increases the frequency of low-precipitable water and light precipitation. It also decreases the frequency of intense precipitation from isolated deep convection and MCSs, increases cloud top temperatures, reduces radar echo-top heights, and increases overall precipitation by altering the relationship of precipitation with precipitable water, in better agreement with observations. The results suggest when optimized using multiple observations, such an approach may provide a path toward more accurate representation of convection and precipitation statistics in convection-permitting simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Conservative velocity mappings for discontinuous Galerkin kinetics

Continuum computational kinetic plasma models evolve the distribution function of a plasma species f s on a phase-space grid over time. In many problems of interest the distribution function has limited extent in velocity space; hence, using a uniform, highly refined mesh would be costly and slow. Nonuniform velocity grids can reduce the computational cost by placing more degrees of freedom where f s is appreciable and fewer where it is not. In this work we introduce a first-of-its kind discontinuous Galerkin approach to nonuniform velocity-space discretization using mapped velocity coordinates. This new method is presented in the context of a gyrokinetic model used to study magnetized plasmas. We create discretizations of collisionless and collisional terms using mappings in a way that exactly conserves particles and energy. Numerical tests of such properties are presented, and we show that this new discretization can reproduce earlier gyrokinetic simulations using grids with up to 6–60 times fewer cells and 22X-60X speed-ups depending on dimensionality, geometry and plasma parameters.

Discontinuous Galerkin↗
DOE OSTI2019

Influence of turbulence and thermophysical fluid properties on cavitation erosion predictions in channel flow geometries

Cavitation and cavitation-induced erosion have been observed in fuel injectors in regions of high acceleration and low pressure. Although these phenomena can have a large influence on the performance and lifetime of injector hardware, questions still remain on how these physics should be accurately and efficiently represented within a computational fluid dynamics model. While several studies have focused on the validation of cavitation predictions within canonical and realistic injector geometries, it is not well documented what influence the numerical and physical parameters selected to represent turbulence and phase change will have on the predictions for cavitation erosion propensity and severity. In this work, a range of numerical and physical parameters are evaluated within the mixture modeling approach in CONVERGE to understand their influence on predictions of cavitation, condensation and erosion. Particular attention is paid to grid resolution, turbulence model and near-wall treatment, fuel surrogate properties, and non-condensable gas content. Assessment of cavitation predictions are conducted through comparison of measured and predicted mass flow rates and cavitation probability distributions for flow through a channel with a sharp inlet. Predictions for hydrodynamic impact loading and cavitation erosion are compared with the experimentally measured incubation period and critical site for erosion. Based on these findings, recommendations are provided for modeling turbulent cavitating flows, using the single fluid mixture modeling approach, to improve predictions for cavitation-induced erosion. In particular, to capture the fluid dynamic phenomena characterizing cavitation cloud formation, development and shedding, a Large Eddy simulation with grid resolution as fine as 2.50 μm is recommended. The assumed concentration of non-condensable gas content is observed to have a strong influence on the predicted cavitation erosion severity, which motivates the need for dissolved gas concentration measurements for future cavitation erosion experimental studies. Using the best practices established in this work, good agreement is achieved between the measured and predicted cavitation parameters, as well as the critical site for cavitation-induced erosion.

Magnotti, Gina Maureen↗
DOE OSTI2017

Investigating the dependence of SCM simulated precipitation and clouds on the spatial scale of large-scale forcing at SGP [Investigating the scale dependence of SCM simulated precipitation and cloud by using gridded forcing data at SGP]

Large-scale forcing data, such as vertical velocity and advective tendencies, are required to drive single-column models (SCMs), cloud-resolving models, and large-eddy simulations. Previous studies suggest that some errors of these model simulations could be attributed to the lack of spatial variability in the specified domain-mean large-scale forcing. This study investigates the spatial variability of the forcing and explores its impact on SCM simulated precipitation and clouds. A gridded large-scale forcing data during the March 2000 Cloud Intensive Operational Period at the Atmospheric Radiation Measurement program's Southern Great Plains site is used for analysis and to drive the single-column version of the Community Atmospheric Model Version 5 (SCAM5). When the gridded forcing data show large spatial variability, such as during a frontal passage, SCAM5 with the domain-mean forcing is not able to capture the convective systems that are partly located in the domain or that only occupy part of the domain. This problem has been largely reduced by using the gridded forcing data, which allows running SCAM5 in each subcolumn and then averaging the results within the domain. This is because the subcolumns have a better chance to capture the timing of the frontal propagation and the small-scale systems. As a result, other potential uses of the gridded forcing data, such as understanding and testing scale-aware parameterizations, are also discussed.

58 GEOSCIENCES↗
DOE OSTI2026

Temporal Convolutional Network Using Empirical Mode Decomposition to Detect Faults in Grid Connected Systems

Grid-connected power electronic systems require timely and reliable fault detection to prevent equipment damage and reduce downtime. This paper presents a forecasting-based anomaly detection pipeline that decomposes voltage and current measurements into intrinsic mode functions (IMFs) using empirical mode decomposition (EMD), then trains a causal temporal convolutional network (TCN) on normal-operation IMF data to predict short-horizon future dynamics. Deviations between forecasts and observations are summarized as reliability-weighted residual scores and thresholded per sensor using robust statistics with temporal persistence constraints to suppress false positives. To reduce runtime, EMD is performed on downsampled signals for detection, while raw-rate EMD is applied only within a short region of interest for high-frequency interpretability near detected events. Results on a simulated grid-connected converter system demonstrate that IMF-domain forecasting improves anomaly separability relative to raw-signal forecasting and provides interpretable evidence of faults across decomposition channels.

Sutton, Elizabeth [ORNL] (ORCID:0009000078885935)↗
DOE OSTI2023

Input files and WRF run directory for a LASSO-CACTI simulation

Tar file containing files necessary for reproducing a given Weather Research and Forecasting (WRF) simulation from the Large-Eddy Simulation (LES) Atmospheric Radiation Measurement (ARM) Symbiotic Simulation and Observation (LASSO) deep-convection scenario for the Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign. The LASSO-CACTI simulations span grid spacings from 7.5 km to 100 m for convection near the Sierras de Córdoba mountain range, roughly centered on the ARM Mobile Facility. More information can be found at https://www.arm.gov/capabilities/modeling/lasso. Types of files in this tar include the initial and boundary conditions, WRF namelist, and other files present in a WRF run directory.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Inverter Model Validation and Calibration Using Phasor Measurement Unit Data

As the penetration of inverter-based renewable energy resources increases in the power grid, especially at the distribution and microgrid levels, the need to accurately represent them in planning studies increases as well. However, due to the lack of well-established standard procedures, and vendor reluctance towards the detailed sharing of proprietary models, automated dynamic model validation and parameter calibration tools for inverter based resources (IBRs) remain scarce. This work presents a model validation and parameter calibration platform for representing IBRs with generic phasor-domain models. Phasor measurements of power system events are used for continuous validation using the data playback method, and model parameters are re-calibrated if a significant mismatch between measurements and model response is observed. Unique features of the proposed platform include- (a) an iterative Bayesian optimization approach towards parameter calibration to address a possible mismatch between the structures of generic models implemented in simulation softwares and actual commercial inverters, (b) error metrics designed to account for a possible mismatch between the time resolution of simulation and measurements, and (c) analysis of the measurement-simulation mismatch to provide guidance to engineering personnel regarding model shortcomings. The performance of the platform has been illustrated using both simulated data and field measurements to validate/calibrate inverter models in GridLAB-D.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Hydrodynamic simulations of white dwarf–white dwarf mergers and the origin of R Coronae Borealis stars

We study the properties of double white dwarf (DWD) mergers by performing hydrodynamic simulations using the new and improved adaptive mesh refinement code OCTO-TIGER. We follow the orbital evolution of DWD systems of mass ratio $q=0.7$ for tens of orbits until and after the merger to investigate them as a possible origin for R Coronae Borealis (RCB) type stars. We reproduce previous results, finding that during the merger, the helium WD donor star is tidally disrupted within 20–80 min since the beginning of the simulation onto the accretor carbon–oxygen WD, creating a high temperature shell around the accretor. We investigate the possible helium burning in this shell and the merged object’s general structure. Specifically, we are interested in the amount of oxygen-16 dredged-up from the accretor to the hot shell and the amount of oxygen-18 produced. This is critical as the discovery of very low oxygen-16 to oxygen-18 ratios in RCB stars pointed out the merger scenario as a favourable explanation for their origin. A small amount of hydrogen in the donor may help keep the oxygen-16 to oxygen-18 ratios within observational bounds, even if moderate dredge-up from the accretor occurs. In addition, we perform a resolution study to reconcile the difference found in the amount of oxygen-16 dredge-up between smoothed-particle hydrodynamics and grid-based simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Advanced grid-forming (GFM) inverter controls, modeling and system impact study for inverter dominated grids

Our main goals in this project are to understand: 1) the stability of the power system with very high penetration of PV or other inverter-based resources (IBR), 2) the interactions among IBRs and traditional generators, and 3) the impact of the GFM inverters on bulk power system. To achieve these goals, we took a holistic research and development approach. First, an impedance-based large system modeling and stability analysis tool was developed to aid the large IBR system small-signal stability analysis. Second, consensus-based control and transient overload ride through were developed to control and coordinate multiple GFM inverters in a large system to maintain transient stability during big grid events. Third, high fidelity electromagnetic transient (EMT) IBR models were developed and integrated into a large reduced WECC model to study the system impact from IBRs. Lastly, extensive validation tests of the GFM control which was developed in this project and implemented in two 1MW commercial solar inverters were performed using different validation platforms, e.g., control hardware-in-the-loop (CHIL) platform, power hardware- in-the-loop (PHIL) platform, and a real PV plant.

14 SOLAR ENERGY↗
DOE OSTI2023

Realistic large eddy and dispersion simulation experiments during project sagebrush phase 1

Typical dispersion simulations use mesoscale models at grid spacing greater than 1 km to simulate pollutant transport. Advances in computing power have allowed Large Eddy Simulations (LES), (grid spacing less than 100 m) to be run with boundary conditions in the grayscale regime (grid size of 1 km–100 m). We used the WRF mesoscale/LES model and the HYSPLIT particle dispersion model to simulate SF 6 transport from Project Sagebrush Phase 1, using a grayscale-aware scheme with 4 nests of 3 km, 1 km, 200 m, 67 m. The WRF/LES simulation was in satisfactory agreement with meteorological observations but the HYSPLIT SF6 simulations indicated insufficient sub-grid turbulence kinetic energy, especially in the lowest 100 m. This resulted in over predictions of the surface SF6 concentration and insufficient vertical transport. Best agreement with the Sagebrush observations was obtained with the WRF/LES large-scale and sub-grid turbulence, enhanced with a HYSPLIT parameterized turbulence profile.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Subsetted model output for process-tendency-related variables from a LASSO-CACTI WRF simulation

Process-tendency-related model output, mostly for the Thompson microphysics scheme (when used), from a Weather Research and Forecasting (WRF) simulation from the Large-Eddy Simulation (LES) Atmospheric Radiation Measurement (ARM) Symbiotic Simulation and Observation (LASSO) deep-convection scenario for the Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign. The LASSO-CACTI simulations span grid spacings from 7.5 km to 100 m for convection near the Sierras de Córdoba mountain range, roughly centered on the ARM Mobile Facility. More information can be found at https://www.arm.gov/capabilities/modeling/lasso. This version of the output is a collection of variables subsetted from the raw WRF output to reduce the file overhead for users not needing the full raw dataset. The subset files are in netCDF format with the height-based variables on the raw, model eta levels.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Subsetted model output for tracer-related variables from a LASSO-CACTI WRF simulation

Inert-tracer-related model output from a Weather Research and Forecasting (WRF) simulation from the Large-Eddy Simulation (LES) Atmospheric Radiation Measurement (ARM) Symbiotic Simulation and Observation (LASSO) deep-convection scenario for the Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign. The LASSO-CACTI simulations span grid spacings from 7.5 km to 100 m for convection near the Sierras de Córdoba mountain range, roughly centered on the ARM Mobile Facility. More information can be found at https://www.arm.gov/capabilities/modeling/lasso. This version of the output is a collection of variables subsetted from the raw WRF output to reduce the file overhead for users not needing the full raw dataset. The subset files are in netCDF format with the height-based variables on the raw, model eta levels.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

On the Influence of Swell Propagation Angle on Surface Drag

Abstract A combination of turbulence-resolving large-eddy simulations and observations are used to examine the influence of swell amplitude and swell propagation angle on surface drag. Based on the analysis a new surface roughness parameterization with nonequilibrium wave effects is proposed. The surface roughness accounts for swell amplitude and wavelength and its relative motion with respect to the mean wind direction. The proposed parameterization is tested in uncoupled three-dimensional Weather and Research Forecasting (WRF) simulations at grid sizes near 1 km where we explore potential implications of our modifications for two-way coupled atmosphere–wave models. Wind–wave misalignment likely explains the large scatter in observed nondimensional surface roughness under swell-dominated conditions. Andreas et al.’s relationship between friction velocity and the 10-m wind speed under predicts the increased drag produced by misaligned winds and waves. Incorporating wave-state (speed and direction) influences in parameterizations improves predictive skill. In a broad sense, these results suggest that one needs information on winds and wave state to upscale buoy measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

The Building ENergy Demand (BEND) Model: A Comprehensive Composite Building Energy Modeling Framework

The Building ENergy Demand (BEND) model was developed at Pacific Northwest National Laboratory (PNNL) to simulate climate-dependent hourly building energy demand with the ability to disaggregate down to or aggregate up to any geographic area of interest, including counties, states, electric utility service territories, and census regions. The BEND modeling framework is designed to study the effects of climate change, economics, technology advancement, and various mitigation policies on aggregate building energy demand. To make the process computationally feasible, BEND reduces granularity in building characteristics and weather forcing while maintaining enough detail to minimize systematic differences (biases) between BEND results and data from the Energy Information Administration (EIA). Based on the Commercial Building Energy Consumption Survey (CBECS) and Residential Energy Consumption Survey (RECS) databases, BEND establishes a set of representative buildings with predefined sets of building types, sizes, and vintages. These representative buildings are defined for each of the four U.S. census regions and the five climate zones used by the EIA for building surveys (hereafter “EIA climate zones”). BEND’s representative buildings are based on a statistical analysis of the CBECS and RECS databases that infers detailed parameters such as construction type, thermal properties, internal loads, equipment and mechanical systems, and operation schedules from the high-level descriptors in those databases (building type, vintage, size, and location). The representative buildings are simplified and consolidated into a set of prototype buildings for the purpose of simulating building energy demand. Combining the prototype buildings with weather files from selected locations and years, BEND creates a set of simulation-ready building models to be run using the EnergyPlus simulation engine (EnergyPlus). The base spatial unit of each BEND simulation is a subregion that is the intersection of an International Energy Conservation Code climate zone, a census region, and an EIA climate zone; analysis for each spatial unit is influenced by the weather file BEND selects to represent the climate-similar area. The major simulation outputs of interest to BEND are hourly energy consumption values for each simulated year, possibly for multiple fuel types. Based on building-stock weights derived from the CBECS/RECS data, the hourly output streams are aggregated to a single, 8760-hour load profile for each of the major building categories—residential and commercial—with the possibility of segregating by fuel types. To correct for biases inherent in the use of a limited suite of representative buildings, BEND is calibrated against observational targets of total annual building energy consumption and hourly load profiles (including peak demand) obtained for the decision-relevant geographic scales of interest (e.g., total electric load data for individual balancing authorities available from EIA). Based on available data such as population distribution, BEND is able to aggregate both up and down scales, allowing the study of building energy demand over a range of spatial scales (e.g., states or counties within the Western Electricity Coordinating Council, Eastern Interconnection, or the entire U.S.) and/or subsets of building types (e.g., office buildings with particular size ranges).

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1993

Moving body overset grid methods for complete aircraft tiltrotor simulations

A hypothetical, but realistic, set of flight conditions for the V-22 aircraft is established to facilitate rigorous testing of a new domain connectivity algorithm, and to carry out an overset grid proof-of-concept tiltrotor simulation. Relative motion and interference effects between the V-22 airframe and rotor-blades are directly simulated within the context of an unsteady, thin-layer Navier-Stokes computation. The domain connectivity algorithm is verified to perform at rates equal to or greater than those realized previously for store-separation-like applications. The feasibility of carrying out unsteady Navier-Stokes analyses of rotorcraft problems is demonstrated.

Meakin, Robert L.↗