SEARCH · Engineering Papers
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.
This presentation provides a high-level overview of NREL's smart home hardware-in-the-loop testing. It was presented at the Fourth International Workshop on Grid Simulator Testing of Energy Systems and Wind Turbine Powertrains, held April 25-26, 2017, hosted by NREL and Clemson University at the Energy Systems Integration Facility in Golden, Colorado.
Building on models of electrical power systems, and on powerful mathematical techniques including optimization, model predictive control, and simluation, this project investigated important issues related to the stable operation of power grids. A topic of particular focus was cascading failures of the power grid: simulation, quantification, mitigation, and control. We also analyzed the vulnerability of networks to component failures, and the design of networks that are responsive to and robust to such failures. Numerous other related topics were investigated, including energy hubs and cascading stall of induction machines
Computational simulations using structured overset grids with the Launch Ascent and Vehicle Aerodynamics (LAVA) solver framework are presented for predicting oblique shock/plume interaction effects to near-field sonic boom signatures. Standard second-order accurate as well as higher-resolution numerical discretizations are utilized and compared in the study. The numerical approach is compared with supersonic wind-tunnel data for three cases. The cases include an empty wind-tunnel at the operating conditions, an isolated shockgenerating diamond wedge within the tunnel, and a nozzle with diamond wedge configuration at five different nozzle pressure ratios. Solution sensitivity to numerical discretization is analyzed. Favorable comparisons between the computational results and experimental data of near-field pressure signatures are obtained. A simple prediction method for plume induced shock deflection is developed and results are compared with the CFD data.
Power grid modernization and wide-spread distributed energy resource (DER) integration are driving an explosion of interest in advanced algorithms and new optimization schemes. However, there has historically been a lack of open-access realistically sized and detailed datasets-particularly for the distribution system-and simulation tools able to model these systems. This talk will present recent advances with very large-scale (>>1M nodes) distribution data sets and the use of co-simulation to tractably model large integrated transmission-distribution grids and their interactions with new controls/optimization.
We describe the load flow formulation and the solution algorithms available in the Electric power Grid Simulator (EGSim) software toolkit. EGSim contains tools aimed at simulating static load flow solutions for electric power grids. It parses power grid models described in IEEE Common Data Format, and generates solutions for the bus voltages and voltage angles, and real and reactive power values through the transmission lines. The software, written in C++, implements both Gauss-Seidel and Newton solution methods. Example results for the 118 bus models and 300 bus models are also presented.
Electricity grid operators around the world face a dual challenge; withstanding increasingly severe weather and the longer term impacts of climate change, while simultaneously decarbonizing. Extreme weather events such as heat waves and droughts are rising in both severity and frequency, which is threatening the reliability of electricity grids through increased demand, generation capacity losses, and equipment failures. Consequently, incorporating hydrometeorological stressors into computational power systems analysis is becoming an even more critical tool in long term planning and short term operations. However, there is a general lack of open-source customizable grid simulation software capable of exhaustively stress testing the grid under hydrometeorological uncertainty, and/or examining potential risk mitigation pathways. A related, persistent challenge for power system modelers is striking an appropriate balance between model fidelity (e.g. spatial scale and time resolution) and computational tractability (wall clock run-time). In this study, we are proposing a solution to this problem with open-source software that allows users to seamlessly customize the scale and track the accuracy of grid operations models. Our approach allows users to search over numerous model parameters (network topology, mathematical formulation, economic hurdle rates, and transmission line scaling) to identify model instantiations that accommodate experimental design. Further, we use this approach to demonstrate the importance of including extreme weather events in model validation and model selection. Focusing on the occurrence of heatwaves and droughts in the U.S. Western Interconnection, we examine role of extreme events in balancing tradeoffs between model fidelity and run-time at the model design stage.
The objective of this thesis is to further develop and test a stochastic model of turbulent combustion in recirculating flows. There is a requirement to increase the accuracy of multi-dimensional combustion predictions. As turbulence affects reaction rates, this interaction must be more accurately evaluated. In this work a more physically correct way of handling the interaction of turbulence on combustion is further developed and tested. As turbulence involves randomness, stochastic modeling is used. Averaged values such as temperature and species concentration are found by integrating the probability density function (pdf) over the range of the scalar. The model in this work does not assume the pdf type, but solves for the evolution of the pdf using the Monte Carlo solution technique. The model is further developed by including a more robust reaction solver, by using accurate thermodynamics and by more accurate transport elements. The stochastic method is used with Semi-Implicit Method for Pressure-Linked Equations. The SIMPLE method is used to solve for velocity, pressure, turbulent kinetic energy and dissipation. The pdf solver solves for temperature and species concentration. Thus, the method is partially familiar to combustor engineers. The method is compared to benchmark experimental data and baseline calculations. The baseline method was tested on isothermal flows, evaporating sprays and combusting sprays. Pdf and baseline predictions were performed for three diffusion flames and one premixed flame. The pdf method predicted lower combustion rates than the baseline method in agreement with the data, except for the premixed flame. The baseline and stochastic predictions bounded the experimental data for the premixed flame. The use of a continuous mixing model or relax to mean mixing model had little effect on the prediction of average temperature. Two grids were used in a hydrogen diffusion flame simulation. Grid density did not effect the predictions except for peak temperature and tangential velocity. The hybrid pdf method did take longer and required more memory, but has a theoretical basis to extend to many reaction steps which cannot be said of current turbulent combustion models.
This paper provides experimental results of a utility-scale (500 kW) solar PV inverter during simulated grid disturbances. This device was tested using simulations of system voltage during three different disturbances; simulated distribution and transmission faults, and voltage records that were collected during the Blue Cut Fire Event. The methods used for testing the three different scenarios of grid disturbances are detailed and the test setup is described. Additionally, the test results for each of the grid events are provided and an analysis of the test results is given. The test results show that the tested PV inverter's response to grid disturbances depends in part on the fault-induced phase shift and harmonic content of the system voltage. These factors can impact the inverters ability to accurately calculate system frequency and can cause erroneous tripping on miscalculated frequency. The information provided in this paper can inform utilities, manufacturers, developers, and system planners of the potential performance of solar PV inverters during grid events.
For many atmospheric monitoring applications, networks of measurement sites—such as the radionuclide stations of the International Monitoring System—can be sparse. With measurement locations potentially hundreds to thousands of kilometers from a release it is important to quantify the effects of physical processes on transport and dispersion of plumes between source and measurement locations. This study addresses the effects of background sources and topography resolution near the release location of radionuclides. We use the Weather Research and Forecasting (WRF) model with inline chemistry to investigate (1) how an additional, time-varying source of 133 Xe, such as an operational medical isotope production facility, contributes to activity concentration measurements at monitoring sites, and (2) how complex topography influences on atmospheric conditions near emission sources impact plume concentrations at varying distances from the source. Two 133 Xe emission sources, including (1) a high flux rate of short duration representing an explosive event, and (2) a variable and continuous background source, are simulated. The continuous background source contributes significantly to total 133 Xe concentrations at several monitoring stations. Further, a WRF simulation at 9 km horizontal resolution is compared with a nested grid simulation, where the innermost domain has a resolution of 1 km. Increased topographic resolution leads to an improved representation of plume responses to local winds, with topographic influences greatest at locations closest to the sources. Differences between the two domain resolutions decrease at greater distances from the sources, as plumes have time to spread and mix and are influenced by synoptic scale circulation patterns that are represented similarly in both simulations.
This report summarizes modeling efforts for the simulation of wave climate along the U.S. West Coast using an unstructured, nested-grid modeling approach that incorporates a global-regional nested grid using WAVEWATCH III® and the high-resolution UnSWAN (Unstructured-grid Simulating Waves Nearshore) model. Wave resource and long-term wave bulk parameters were simulated and subsequently validated with wave buoy data within the model domain. Overall, model results match observations at buoys well. The model hindcast was able to reproduce the seasonal variation of the sea state with large waves that occur in the winter and early spring months when wind forcing is strong and the calm sea state during the summer when wind forcing is weak. The nested-grid modeling framework employed in this study provides a powerful and efficient modeling approach to accurately simulate wave climate at regional and long-term temporal scales with sufficiently fine resolutions in the nearshore region. Results can be used to inform the deployment of wave energy converters and other nearshore instrumentation/devices in nearshore regions, assist in prioritizing hotspots for near-term market opportunities and development, and provide wave resource assessment data and information to inform resource estimates for marine and hydrokinetic energy technologies.
This dataset provides 0.125° gridded simulations of soil moisture and crop grain carbon for the Contiguous United States (CONUS), generated using the Community Land Model version 5 (CLM5) with the biogeochemistry module enabled. The data covers a historical baseline (1980–2015) and mid-century future projections (2020–2055). Future projections are organized into two sets of scenarios to distinguish the impacts of different drivers: (1) Atmospheric Only (ATM-only): These scenarios apply future atmospheric forcings while holding land use and land cover change (LULCC) at historical baseline levels. The atmospheric forcings represent moderately versus severely hotter/drier atmospheric conditions (dynamically downscaled perturbed thermodynamics simulations based on CMIP6 SSP245 and SSP585 warming signals), each with cooler versus hotter Earth System Model temperature sensitivity instantiations. These scenarios are identified in the file names as atm45cooler, atm45hotter, atm85cooler, and atm85hotter. (2) Coupled Atmospheric and Land-Use (LAND+ATM): These scenarios apply future atmospheric forcing together with future LULCC by pairing atmospheric pathways with lower versus higher population/economic growth scenarios representing Shared Socioeconomic Pathways 3 and 5 (SSP3 and SSP5). These scenarios are identified in the file names as atm45cooler_ssp3, atm45hotter_ssp3, atm45cooler_ssp5, atm45hotter_ssp5, atm85cooler_ssp3, atm85hotter_ssp3, atm85cooler_ssp5, and atm85hotter_ssp5. Please refer to the README file for detailed information on file structure, variables, units, and data formats.
Thermochemical nonequilibrium flow simulation capabilities have been previously implemented, verified, and validated for central processing unit (CPU) systems in NASA’s unstructured-grid computational fluid dynamics solver FUN3D. Many exascale-class high-performance computing systems will rely on graphics processing unit (GPU) architectures for high throughput and energy efficiency; thus, CPU-based scientific computing software unable to effectively utilize these systems must be updated. In this work, we present a CUDA C++ implementation of FUN3D’s thermochemical nonequilibrium flow simulation capabilities targeting NVIDIA Tesla GPUs. An overview of the porting and optimization strategy is described and performance comparisons with other recent architectures are presented. Scaling to thousands of GPUs is demonstrated, yielding computational performance equivalent to that of several million CPU cores. The implementation enables efficient, high-fidelity, scale-resolving simulations of thermochemical nonequilibrium flows for many applications including atmospheric entry, hypersonics, and combustion.
The devastating impact of extreme weather-related events is increasingly evident on power grids, especially on distribution grids. The severity of their potential impact calls for 1) developing a suitable resilience assessment framework to capture the system performance and 2) assessing relevant mitigative strategies to lessen the impact of such events. This paper proposes a framework to identify grid vulnerabilities using the energy-at-risk concept to select, disconnect and isolate grid portions due to a resilience event. The proposed framework mainly consists of two steps; i) processing the utility's available infrastructure, i.e., a network model, possible switching combinations, and outage information for those combinations, as a graph-based database, and ii) implementing a novel optimal switching algorithm leveraging database and grid simulated metrics. These switching actions are generated to implement load curtailment in a rolling manner during anticipated grid scarcity conditions. In this study, a test case is created using two taxonomy feeders and is simulated against an extreme temperature event, e.g., a long, relatively cold, and prolonged freeze peak, thereby creating stress on the grid. It is demonstrated that the proposed framework allows utilities to predict the energy-at-risk during such resilience events and design suitable outage management strategies.
The time-averaged information from a 10-million-cell multiphase computational fluid dynamic simulation of a coal gasifier was used as the inlet boundary condition for fine grid simulations of the region near the coal inlet. The simulation results show a distinct coal jet region, characterized by high void fraction and low solids temperature. The jets from right and left inlets have different behaviors caused by the asymmetry in the inlet to the gasifier. The medium grid resolution with a high order discretization scheme offers better jet resolution at lower computational cost than a fine grid with a low order discretization scheme.
Observation-based modeling case studies of continental boundary layer clouds have been developed to study cloudy boundary layers, aerosol influences upon them, and their representation in cloud- and global-scale models. Three 60-hour case study periods span the temporal evolution of cumulus, stratiform, and drizzling boundary layer cloud systems, representing mixed and transitional states rather than idealized or canonical cases. Based on in-situ measurements from the RACORO field campaign and remote-sensing observations, the cases are designed with a modular configuration to simplify use in large-eddy simulations (LES) and single-column models. Aircraft measurements of aerosol number size distribution are fit to lognormal functions for concise representation in models. Values of the aerosol hygroscopicity parameter, κ, are derived from observations to be ~0.10, which are lower than the 0.3 typical over continents and suggestive of a large aerosol organic fraction. Ensemble large-scale forcing datasets are derived from the ARM variational analysis, ECMWF forecasts, and a multi-scale data assimilation system. The forcings are assessed through comparison of measured bulk atmospheric and cloud properties to those computed in 'trial' large-eddy simulations, where more efficient run times are enabled through modest reductions in grid resolution and domain size compared to the full-sized LES grid. Simulations capture many of the general features observed, but the state-of-the-art forcings were limited at representing details of cloud onset, and tight gradients and high-resolution transients of importance. Methods for improving the initial conditions and forcings are discussed. The cases developed are available to the general modeling community for studying continental boundary clouds.
Data-driven models of power system inverter-based resources are desired to run simulations faster than with detailed electromagnetic transient models, to hide proprietary design details, to support control system design applications, and to aggregate the effects of distributed energy resources. This paper applies a customized Hammerstein Wiener framework to train block diagram models from thousands of electromagnetic transient simulations or experimental test records. The block diagram models integrate with larger grid simulations as voltagecontrolled current sources or current-controlled voltage sources for several simulators. Guidelines for block architecture and training are presented. Three-phase balanced, three-phase unbalanced, and single-phase examples all achieve an acceptable root mean square error of no more than 0.05 per-unit.
This research project will estimate the value to the United States electric grid of deploying hydrogen technology (such as electrolyzers and hydrogen-fueled generation) under projected conditions of high renewable penetration. The analysis will advance the state of the art in systems-level cost-benefit analysis of hydrogen technology for the electric grid by incorporating production cost modeling results in the analysis. Large-scale grid simulation tools will be used to evaluate total system production cost and grid operation when hydrogen technology is deployed for applications such as energy storage and demand response. Scenarios will include one or more future grid mixes in the Western Interconnect (WI) with a high proportion of intermittent renewables. Electric Power Research Institute (EPRI) will work with four utility companies to refine scenarios. Results of the analysis will include comparing the net cost of hydrogen to other technologies for long duration storage, and of power-to-gas (P2G) scenarios including merchant hydrogen sale and hydrogen-fueled generation.
Observation-based modeling case studies of continental boundary layer clouds have been developed to study cloudy boundary layers, aerosol influences upon them, and their representation in cloud- and global-scale models. Three 60 h case study periods span the temporal evolution of cumulus, stratiform, and drizzling boundary layer cloud systems, representing mixed and transitional states rather than idealized or canonical cases. Based on in situ measurements from the Routine AAF (Atmospheric Radiation Measurement (ARM) Aerial Facility) CLOWD (Clouds with Low Optical Water Depth) Optical Radiative Observations (RACORO) field campaign and remote sensing observations, the cases are designed with a modular configuration to simplify use in large-eddy simulations (LES) and single-column models. Aircraft measurements of aerosol number size distribution are fit to lognormal functions for concise representation in models. Values of the aerosol hygroscopicity parameter, kappa, are derived from observations to be approximately 0.10, which are lower than the 0.3 typical over continents and suggestive of a large aerosol organic fraction. Ensemble large-scale forcing data sets are derived from the ARM variational analysis, European Centre for Medium-Range Weather Forecasts, and a multiscale data assimilation system. The forcings are assessed through comparison of measured bulk atmospheric and cloud properties to those computed in "trial" large-eddy simulations, where more efficient run times are enabled through modest reductions in grid resolution and domain size compared to the full-sized LES grid. Simulations capture many of the general features observed, but the state-of-the-art forcings were limited at representing details of cloud onset, and tight gradients and high-resolution transients of importance. Methods for improving the initial conditions and forcings are discussed. The cases developed are available to the general modeling community for studying continental boundary clouds.