Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “grid simulator”

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 703 records · Page 39

Analyzing impact of $\mathrm{DER}$ on $\mathrm{FIDVR}$ - comparison of $\mathrm{EMT}$ simulation of a combined transmission and distribution grid with aggregated positive sequence models

The increase in percentage of distributed energy resources (DERs) brings in a degree of uncertainty in bulk power system planning and operations. Further, the presence of 1 - φ induction motors introduces the possibility of motor stalling and fault induced delayed voltage recovery (FIDVR). With newer versions of IEEE Std 1547™ there can be the opportunity to obtain dynamic voltage support from DERs. In this report, the impact of DERs providing dynamic voltage on the dynamic behavior of 1 - φ induction motors is examined in detail along with an exploration of the need for simulation platforms that go beyond positive sequence. The studies discussed here are carried out both in detailed electromagnetic time domain and in positive sequence domain for comparison.

24 POWER TRANSMISSION AND DISTRIBUTION↗

An RTDS-Based Testbed for Investigating the Impacts of Transmission-Level Disturbances on Solar PV Operation

The increasing penetration of renewable energy resources such as solar and wind via power electronic inverters is challenging grid dynamics, as well as grid planning, operation, and protection. Recently, the North American Electric Reliability Corporation (NERC) has reported a series of similar events of the unintended loss of solar generation in Southern California over a large geographic area following the transmission-level disturbances. These events highlight the importance of understanding the characteristics of the transmission-side disturbances propagating into the distribution systems and their impacts on the operation of inverter-based resources. In this paper, a real-time electromagnetic simulation testbed is constructed for real-time electromagnetic simulations to generate realistic transmission-level disturbances and investigate their impacts on the solar PV operation under different fault types and locations, solar penetration levels, and loading levels. Through the simulation analysis and grid strength assessment, it is found that the grid strength at points of integration (POIs) of solar PVs significantly affects the transient stability of solar generators. Particularly, undesirable transient stability events are more likely to occur at the weak POIs following the transmission-level disturbances. Moreover, undesirable transient stability events become severer when the transmission-level disturbance is closer to the weak POIs or the disturbances become more serious. Additionally, the impact of the transmission-level disturbances on the solar PVs at the weak POIs exacerbate with the increasing solar penetration levels and loading levels. Thus, it is important to study and develop new technologies for grid planning, operation, and protection in weak grid conditions to address the emerging issues of integrating the high penetration of solar PVs and other IBRs.

14 SOLAR ENERGY↗

Separating Physics and Dynamics Grids for Improved Computational Efficiency in Spectral Element Earth System Models

Previous studies have shown that atmospheric models with a spectral element grid can benefit from putting physics calculations on a relatively coarse finite volume grid. Here we demonstrate an alternative high-order, element-based mapping approach used to implement a quasi-equal-area, finite volume physics grid in E3SM. Unlike similar methods, the new method in E3SM requires topology data purely local to each spectral element, which trivially allows for regional mesh refinement. Simulations with physics grids defined by 2 × 2, 3 × 3, and 4 × 4 divisions of each element are shown to verify that the alternative physics grid does not qualitatively alter the model solution. The model performance is substantially affected by the reduction of physics columns when using the 2 × 2 grid, which can increase the throughput of physics calculations by roughly 60%–120% depending on whether the computational resources are configured to maximize throughput or efficiency. A pair of regionally refined cases are also shown to highlight the refinement capability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Global Gridded Crop Model Evaluation: Benchmarking, Skills, Deficiencies and Implications.

Crop models are increasingly used to simulate crop yields at the global scale, but so far there is no general framework on how to assess model performance. Here we evaluate the simulation results of 14 global gridded crop modeling groups that have contributed historic crop yield simulations for maize, wheat, rice and soybean to the Global Gridded Crop Model Intercomparison (GGCMI) of the Agricultural Model Intercomparison and Improvement Project (AgMIP). Simulation results are compared to reference data at global, national and grid cell scales and we evaluate model performance with respect to time series correlation, spatial correlation and mean bias. We find that global gridded crop models (GGCMs) show mixed skill in reproducing time series correlations or spatial patterns at the different spatial scales. Generally, maize, wheat and soybean simulations of many GGCMs are capable of reproducing larger parts of observed temporal variability (time series correlation coefficients (r) of up to 0.888 for maize, 0.673 for wheat and 0.643 for soybean at the global scale) but rice yield variability cannot be well reproduced by most models. Yield variability can be well reproduced for most major producing countries by many GGCMs and for all countries by at least some. A comparison with gridded yield data and a statistical analysis of the effects of weather variability on yield variability shows that the ensemble of GGCMs can explain more of the yield variability than an ensemble of regression models for maize and soybean, but not for wheat and rice. We identify future research needs in global gridded crop modeling and for all individual crop modeling groups. In the absence of a purely observation-based benchmark for model evaluation, we propose that the best performing crop model per crop and region establishes the benchmark for all others, and modelers are encouraged to investigate how crop model performance can be increased. We make our evaluation system accessible to all crop modelers so that other modeling groups can also test their model performance against the reference data and the GGCMI benchmark.

wheat↗

Rapid Load Transition for Integrated Solid Oxide Fuel Cell - Gas Turbine (SOFC-GT) Energy Systems: A Demonstration of the Potential for Grid Response

Rapid load transition is an essential requirement for integrated energy systems to maintain grid resilience as more renewable resources are added to the grid. Integrated solid oxide fuel cell - gas turbine (SOFC-GT) systems can provide high efficiency and low carbon emissions over a broad range of turndown. These hybrids also have the potential to enable rapid grid response. The challenge has been to demonstrate effective control strategies to manage load transitions. In the present study, a load transition of ~50% was achieved in 10 s using a novel but simple strategy. Power demand on the SOFC and the GT were ramped down concurrently. During this transition, the SOFC anode fuel was manipulated to maintain SOFC fuel utilization while the cathode inlet air flow and temperature were also manipulated to thermally protect the SOFC. This study was conducted using the Hybrid Performance (Hyper) facility at the National Energy Technology Laboratory (NETL) in a co-simulation environment with the Idaho National Laboratory (INL)'s grid-simulation. The load ramping strategy was tested using a hardware-based cyber-physical simulation methodology. The results demonstrate a high-fidelity representation of SOFC-GT hybrid dynamics and validation of the control strategy. Thermal and electrochemical transients indicated that the SOFC was well protected during rapid load turndown without violating operability constraints. This demonstration revealed the non-linear nature of tightly coupled SOFC-GT system components, especially the non-linear response of SOFC cathode air flow and inlet temperature controls. These results highlight the needs and challenges in developing adaptive automatic controls for autonomous rapid load transitions. This work demonstrates that SOFC-GT hybrids are a viable option to provide the fast-ramping characteristics essential to accommodate high levels of variable renewable power while maintaining grid resilience, reliability, and environmental performance. The results also demonstrate the utility of co-simulation in advancing the tightly-coupled integrated energy systems needed to meet goals for zero-carbon power generation.

DIRECT ENERGY CONVERSION,POWER TRANSMISSION AND DI↗

Laboratory Resources and Techniques to Evaluate Smart Home Technology: Preprint

Smart home technologies promise to make residential buildings more comfortable, less costly to operate, more energy efficient, and more integrated with the grid, but these systems are difficult to evaluate in traditional laboratory spaces. To aid in the development and evaluation of new home energy management systems (HEMS), NREL built the Systems Performance Laboratory (SPL) within the Energy Systems Integration Facility. SPL has electrical infrastructure for three residential buildings, complete with smart meters, breaker panels and a full complement of major and minor appliances in each 'home.' Cord and plug power connections mean that appliances can easily be swapped for different models. Distributed energy resources (DER), such as residential photovoltaic (PV) inverters and home batteries, are also available for integration. Sensors monitor the power consumption of all loads and controllers can simulate occupancy to drive the appliances in realistic ways. A novel hardware-in-the-loop application was developed to ensure air conditioner cycling by combining a building simulation with a physical air conditioner and a thermostat, which allows us to simulate performance in any climate with any building construction. Another application of the hardware-in-the-loop technique uses a grid simulation to model the grid impact for different control strategies related to smart home technologies. This laboratory is a powerful, flexible tool that can be employed to evaluate the effect that home energy management and other smart home technologies will have at the appliance, whole-home, and distribution feeder levels - all of which are important perspectives for different stakeholders.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

GLOBAL SIMULATIONS OF ACCRETION DISKS. I. CONVERGENCE AND COMPARISONS WITH LOCAL MODELS

Grid-based magnetohydrodynamic (MHD) simulations have proven invaluable for the study of astrophysical accretion disks. However, the fact that angular momentum transport in disks is mediated by MHD turbulence (with structure down to very small scales) raises the concern that the properties of the modeled accretion disks are affected by the finite numerical resolution of the simulation. By implementing an orbital advection algorithm into the Athena code in cylindrical geometry, we have performed a set of global (but unstratified) Newtonian disk simulations extending up to resolutions previously unattained. We study the convergence of these models as a function of spatial resolution and initial magnetic field geometry. The usual viscosity parameter ({alpha}) or the ratio of thermal-to-magnetic pressure ({beta}) is found to be a poor diagnostic of convergence, whereas the average tilt angle of the magnetic field in the (r, {phi})-plane is a very good diagnostic of convergence. We suggest that this is related to the saturation of the MHD turbulence via parasitic modes of the magnetorotational instability. Even in the case of zero-net magnetic flux, we conclude that our highest resolution simulations (with 32 zones and 64 zones per vertical scale height) have achieved convergence. Our global simulations reach resolutions comparable to those used in local, shearing-box models of MHD disk turbulence. We find that the saturation predictors derived from local simulations correspond well to the instantaneous correlations between local flux and stress found in our global simulations. However, the conservation of magnetic flux implicit in local models is not realized in our global disks. Thus, the magnetic connectivity of an accretion disk represents physics that is truly global and cannot be captured in any ab initio local model.

79 ASTRONOMY AND ASTROPHYSICS↗

Chance Constrained Optimal Power Flow: Risk-Aware Network Control under Uncertainty

When uncontrollable resources fluctuate, Optimum Power Flow (OPF), routinely used by the electric power industry to re-dispatch hourly controllable generation (coal, gas and hydro plants) over control areas of transmission networks, can result in grid instability, and, potentially, cascading outages. This risk arises be cause OPF dispatch is computed without awareness of major uncertainty, in particular fluctuations in renewable output. As a result, grid operation under the OPF with renewable variability can lead to frequent conditions where power line flow ratings are significantly exceeded. Such a condition, which is borne by simulations of real grids, would likely resulting in automatic line tripping to protect lines from thermal stress, a risky and undesirable outcome which compromises stability. Smart grid goals include a commitment to large penetration of highly fluctuating renewables, thus calling to reconsider current practices, in particular the use of standard OPF. Our Chance Constrained (CC) OPF corrects the problem and mitigates dangerous renewable fluctuations with minimal changes in the current operational procedure. Assuming availability of a reliable wind forecast parameterizing the distribution function of the uncertain generation, our CC-OPF satisfies all the constraints with high probability while simultaneously minimizing the cost of economic re-dispatch. CC-OPF allows efficient implementation, e.g. solving a typical instance over the 2746-bus Polish network in 20s on a standard laptop.

Energy Planning, Policy, & Economy(29)↗

Impact of Lidar Wind Sounding on Mesoscale Forecast

An Observing System Simulation Experiment (OSSE) was conducted to study the impact of airborne lidar wind sounding on mesoscale weather forecast. A wind retrieval scheme, which interpolates wind data from a grid data system, simulates the retrieval of wind profile from a satellite lidar system. A mesoscale forecast system based on the PSU/NCAR MM5 model is developed and incorporated the assimilation of the retrieved line-of-sight wind. To avoid the "identical twin" problem, the NCEP reanalysis data is used as our reference "nature" atmosphere. The simulated space-based lidar wind observations were retrieved by interpolating the NCEP values to the observation locations. A modified dataset obtained by smoothing the NCEP dataset was used as the initial state whose forecast was sought to be improved by assimilating the retrieved lidar observations. Forecasts using wind profiles with various lidar instrument parameters has been conducted. The results show that to significantly improve the mesoscale forecast the satellite should fly near the storm center with large scanning radius. Increasing lidar firing rate also improves the forecast. Cloud cover and lack of aerosol degrade the quality of the lidar wind data and, subsequently, the forecast.

Miller, Timothy L.↗

Geometry Modeling and Grid Generation for "Icing Effects" and "Ice Accretion" Simulations on Airfoils

There are two distinct icing-related problems for airfoils that can be simulated. One is predicting the effects of ice on the aerodynamic performance of airfoils when ice geometry is known ("icing effects" study). The other is simulating ice accretion under specified icing conditions ("ice accretion" simulation). This paper will address development of two different software packages for two-dimensional geometry preparation and grid generation for both "icing effects" and "ice accretion" studies.

Choo, Yung↗

A Sparse Grid Method for Bayesian Uncertainty Quantification with Application to Large Eddy Simulation Turbulence Models

There is wide agreement that the accuracy of turbulence models suffer from their sensitivity with respect to physical input data, the uncertainties of user-elected parameters, as well as the model inadequacy. However, the application of Bayesian inference to systematically quantify the uncertainties in parameters, by means of exploring posterior probability density functions (PPDFs), has been hindered by the prohibitively daunting computational cost associated with the large number of model executions, in addition to daunting computation time per one turbulence simulation. In this effort, we perform in this paper an adaptive hierarchical sparse grid surrogate modeling approach to Bayesian inference of large eddy simulation (LES). First, an adaptive hierarchical sparse grid surrogate for the output of forward models is constructed using a relatively small number of model executions. Using such surrogate, the likelihood function can be rapidly evaluated at any point in the parameter space without simulating the computationally expensive LES model. This method is essentially similar to those developed in Zhang et al. (Water Resour Res 49:6871–6892, 2013) for geophysical and groundwater models, but is adjusted and applied here for a much more challenging problem of uncertainty quantification of turbulence models. Through a numerical demonstration of the Smagorinsky model of two-dimensional flow around a cylinder at sub-critical Reynolds number, our approach is proven to significantly reduce the number of costly LES executions without losing much accuracy in the posterior probability estimation. Here, the model parameters are calibrated against synthetic data related to the mean flow velocity and Reynolds stresses at different locations in the flow wake. The influence of the user-elected LES parameters on the quality of output data will be discussed.

97 MATHEMATICS AND COMPUTING↗

Propulsion Simulations Using Advanced Turbulence Models with the Unstructured Grid CFD Tool, TetrUSS

A computational investigation has been completed to assess the capability of TetrUSS for exhaust nozzle flows. Three configurations were chosen for this study (1) an axisymmetric supersonic jet, (2) a transonic axisymmetric boattail with solid sting operated at different Reynolds number and Mach number, and (3) an isolated non-axisymmetric nacelle with a supersonic cruise nozzle. These configurations were chosen because existing experimental data provided a means for measuring the ability of TetrUSS for simulating complex nozzle flows. The main objective of this paper is to validate the implementation of advanced two-equation turbulence models in the unstructured-grid CFD code USM3D for propulsion flow cases. USM3D is the flow solver of the TetrUSS system. Three different turbulence models, namely, Menter Shear Stress Transport (SST), basic k epsilon, and the Spalart-Allmaras (SA) are used in the present study. The results are generally in agreement with other implementations of these models in structured-grid CFD codes. Results indicate that USM3D provides accurate simulations for complex aerodynamic configurations with propulsion integration.

Abdol-Hamid, Khaled S.↗

Effect of Terrestrial and Marine Organic Aerosol on Regional and Global Climate: Model Development, Application, and Verification with Satellite Data

In this DOE project the improvements to parameterization of marine primary organic matter (POM) emissions, hygroscopic properties of marine POM, marine isoprene derived secondary organic aerosol (SOA) emissions, surfactant effects, new cloud droplet activation parameterization have been implemented into Community Atmosphere Model (CAM 5.0), with a seven mode aerosol module from the Pacific Northwest National Laboratory (PNNL)'s Modal Aerosol Model (MAM7). The effects of marine aerosols derived from sea spray and ocean emitted biogenic volatile organic compounds (BVOCs) on microphysical properties of clouds were explored by conducting 10 year CAM5.0-MAM7 model simulations at a grid resolution 1.9° by 2.5° with 30 vertical layers. Model-predicted relationship between ocean physical and biological systems and the abundance of CCN in remote marine atmosphere was compared to data from the A-Train satellites (MODIS, CALIPSO, AMSR-E). Model simulations show that on average, primary and secondary organic aerosol emissions from the ocean can yield up to 20% increase in Cloud Condensation Nuclei (CCN) at 0.2% Supersaturation, and up to 5% increases in droplet number concentration of global maritime shallow clouds. Marine organics were treated as internally or externally mixed with sea salt. Changes associated with cloud properties reduced (absolute value) the model-predicted short wave cloud forcing from -1.35 Wm -2 to -0.25 Wm -2 . By using different emission scenarios, and droplet activation parameterizations, this study suggests that addition of marine primary aerosols and biologically generated reactive gases makes an important difference in radiative forcing assessments. All baseline and sensitivity simulations for 2001 and 2050 using global-through-urban WRF/Chem (GU-WRF) were completed. The main objective of these simulations was to evaluate the capability of GU-WRF for an accurate representation of the global atmosphere by exploring the most accurate configuration of physics options in GWRF for global scale modeling in 2001 at a horizontal grid resolution of 1° x 1°. GU-WRF model output was evaluated using observational datasets from a variety of sources including surface based observations (NCDC and BSRN), model reanalysis (NCEP/ NCAR Reanalysis and CMAP), and remotely-sensed data (TRMM) to evaluate the ability of GU-WRF to simulate atmospheric variables at the surface as well as aloft. Explicit treatment of nanoparticles produced from new particle formation in GU-WRF/Chem-MADRID was achieved by expanding particle size sections from 8 to 12 to cover particles with the size range of 1.16 nm to 11.6m. Simulations with two different nucleation parameterizations were conducted for August 2002 over a global domain at a 4º by 5º horizontal resolution. The results are evaluated against field measurement data from the 2002 Aerosol Nucleation and Real Time Characterization Experiment (ANARChE) in Atlanta, Georgia, as well as satellite and reanalysis data. We have also explored the relationship between clean marine aerosol optical properties and ocean surface wind speed using remotely sensed data from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) on board the CALIPSO satellite and the Advanced Microwave Scanning Radiometer (AMSR-E) on board the AQUA satellite. Detailed data analyses were carried out over 15 regions selected to be representative of different areas of the global ocean for the time period from June 2006 to April 2011. We show that for very low (less than 4 m s -1 ) and very high (more than 12 m s -1 ) wind speed conditions the mean CALIPSO-derived aerosol optical depth (AOD) has little dependency on the surface wind speed. For an intermediate (between 4 and 12 m s -1 ) marine AOD was linearly correlated with the surface wind speed values, with a slope of 0.0062 s m -1 . Results of our study suggest that considerable improvements to both optical properties of marine aerosols and their production mechanisms can be achieved by discriminating clean marine aerosols (or sea salt particles) from all other types of aerosols present over the ocean.

54 ENVIRONMENTAL SCIENCES↗

Simulation of wave propagation in three-dimensional random media

Quantitative error analysis for simulation of wave propagation in three dimensional random media assuming narrow angular scattering are presented for the plane wave and spherical wave geometry. This includes the errors resulting from finite grid size, finite simulation dimensions, and the separation of the two-dimensional screens along the propagation direction. Simple error scalings are determined for power-law spectra of the random refractive index of the media. The effects of a finite inner scale are also considered. The spatial spectra of the intensity errors are calculated and compared to the spatial spectra of intensity. The numerical requirements for a simulation of given accuracy are determined for realizations of the field. The numerical requirements for accurate estimation of higher moments of the field are less stringent.

Coles, William A.↗

Properties of a Variable-Delay Polarization Modulator

We investigate the polarization modulation properties of a variable-delay polarization modulator (VPM). The VPM modulates polarization via a variable separation between a polarizing grid and a parallel mirror. We find that in the limit where the wavelength is much larger than the diameter of the metal wires that comprise the grid, the phase delay derived from the geometric separation between the mirror and the grid is sufficient to characterize the device. However, outside of this range, additional parameters describing the polarizing grid geometry must be included to fully characterize the modulator response. In this paper, we report test results of a VPM at wavelengths of 350 micron and 3 mm. Electromagnetic simulations of wire grid polarizers were performed and are summarized using a simple circuit model that incorporates the loss and polarization properties of the device.

Chuss, David T.↗

Validation of a Three-Dimensional Ablation and Thermal Response Simulation Code

The 3dFIAT code simulates pyrolysis, ablation, and shape change of thermal protection materials and systems in three dimensions. The governing equations, which include energy conservation, a three-component decomposition model, and a surface energy balance, are solved with a moving grid system to simulate the shape change due to surface recession. This work is the first part of a code validation study for new capabilities that were added to 3dFIAT. These expanded capabilities include a multi-block moving grid system and an orthotropic thermal conductivity model. This paper focuses on conditions with minimal shape change in which the fluid/solid coupling is not necessary. Two groups of test cases of 3dFIAT analyses of Phenolic Impregnated Carbon Ablator in an arc-jet are presented. In the first group, axisymmetric iso-q shaped models are studied to check the accuracy of three-dimensional multi-block grid system. In the second group, similar models with various through-the-thickness conductivity directions are examined. In this group, the material thermal response is three-dimensional, because of the carbon fiber orientation. Predictions from 3dFIAT are presented and compared with arcjet test data. The 3dFIAT predictions agree very well with thermocouple data for both groups of test cases.

Chen, Yih-Kanq↗

Development of a Configurable DERMS Test System in GridAPPS-D

Distributed Energy Resources and their increased penetration within the power grid are a suitable method to provide essential reliability services. In order to do so, they must operate in aggregate. Aggregation of a wide variety of systems with varying standards and protocols is a complex task, which could benefit from rules and standards to support interoperability. To this end, the Portland State University Power Lab has developed the Energy Grid of Things project, which includes a Distributed Energy Resource Management System to provide aggregation and demonstrate the usage of an Energy Services Interface: a set of rules governing the interface between resource owners and aggregator businesses. This thesis presents a Modeling Environment that provides a simulation and test platform reflecting the complex needs of an aggregator. A grid simulator is required that is able to reflect the effects of thousands of Distributed Energy Resources operating in tandem; furthermore, these effects must be communicated to external actors such as the Grid Operator, which provides grid service requests to the aggregator and requires feedback. The Modeling Environment provides all of this functionality by leveraging the GridAPPS-D platform to develop an interoperable class-based system that provides configurable simulations with modular input and output APIs. The Modeling Environment uses abstractions to allow any feasible type of input to be reflected within the model as a Distributed Energy Resource, and serves as a proof-of-concept for the development of aggregation simulation and testing within GridAPPS-D.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗