Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “grid operators”

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 1,279 records · Page 71

Momentary Cessation: Improving Dynamic Performance and Modeling of Utility-Scale Inverter Based Resources During Grid Disturbances

Sandia National Laboratories worked with NERC staff to provide stakeholder guidance in responding to a May 2018 NERC alert regarding dynamic performance and modeling issues for utility-scale inverter-based resources. The NERC alert resulted from event analyses for grid disturbances that occurred in southern California in August 2016 and October 2017. Those disturbances revealed the use of momentary cessation of transmission connected inverter-based generation- a short time period when they ceased to inject current into the grid, counter to desired transmission operation. The event analyses concluded that, in many cases, the Western Interconnection system models used to determine planning and operating criteria do not reflect the actual behavior of solar plants, resulting in overly optimistic planning assessments and substandard operational responses. This technical report summarizes the gaps between the models and actual performance that were observed at those times, and the guidance that Sandia and NERC provided to owners of solar PV power plants, transmission planners, transmission operators and planning/reliability coordinators to modify existing models to reflect that actual performance

24 POWER TRANSMISSION AND DISTRIBUTION↗

Big PanDa Workflow Management on Titan for High Energy and Nuclear Physics and for Future Extreme Scale Scientific Application

Over a three year period, from 2016-2019, this project demonstrated the scientific benefits of integrating the Titan supercomputer at Oak Ridge Leadership Computing Facility into traditional high throughput grid based distributed computing systems managed by PanDA, the workflow management system used for the execution of all distributed computing applications by the ATLAS experiment at the Large Hadron Collider. PanDA manages millions of batch jobs daily at hundreds of clusters worldwide on request by thousands of physicist users, and processes more than an exabyte of data annually using grid middleware. High levels of operational use of Titan was sustained by PanDA in order to meet the physics goals of ATLAS. The success of this project led to the use of other supercomputers worldwide by ATLAS, and to the adoption of PanDA by other experiments and other scientists. Multiple innovative operational and computer science research goals were achieved supporting the use of supercomputers for scientific domains with large scale distributed data and distributed processing needs.

97 MATHEMATICS AND COMPUTING↗

Ensemble Manufacturing Techniques for Steam Turbine Components Across Length Scales

Faster design to manufacturing requirements were sought for steam turbine components to meet the changing demands of today’s power grid of improved efficiency through operating temperature increases and enhanced operational flexibility from baseload to cyclic operations. Advances in multiple advanced manufacturing (AM) process enabled redesign of turbine components for extreme environments with the potential to reduce cost. AM is of particular interest to improve component functionality, higher temperature capability, and superior durability in large scale steam turbine applications. AM methods have an immense potential to open-up the design space by working directly with the 3D model to produce near-net shapes, enable fast design-manufacturing iterations, thereby significantly reducing product cost and lead-time up to 25 % from current baseline. However, these benefits cannot be fully realized due to the potential for unknown AM processing defects and their resultant effect upon component performance in service. Siemens is partnering with Oak Ridge National Laboratory (ORNL), Electric Power Research Institute (EPRI), and Connecticut Center for Advanced Technology (CCAT) to advance the knowledge of complex process-material interactions for desired microstructures and properties that are closely interlinked to component geometries across different length scales. The proposed program utilized an ensemble of multidisciplinary technologies to accelerate the development of materials, high-throughput experiments for their qualification, and design flexibility/topology optimization for repair/redesign of components to address critical failure mechanisms for improved performance and increased reliability of existing power plant components. The proposed activities, if successfully demonstrated for identified components, will enable paradigm shift in customized manufacturing and accelerated qualification/certification towards increased steam turbine component durability, increased turbine efficiency, and reduced CO 2 emissions in load-following environments compared to today’s technology. Technology maturation is built into the project as successful research will include customized process-component down-selection enabling AM methodologies to be incorporated directly into the existing supply chain. The specific activities of the proposed effort are: 1. Topology optimization of down-selected steam turbine parts that are amenable to additive and hybrid manufacturing for cost/performance improvement. 2. Process-structure-property relationships for five AM processes for steam turbine materials of interest to compare with conventional materials. 3. Perform part/assembly build process using advanced additive/hybrid machine tools followed by quality inspection of the built components for insight into qualification for production scale-up. Steam turbine rig testing of printed components under targeted, well monitored and characterized environmental conditions of for performance comparison of baseline and redesigned components.

36 MATERIALS SCIENCE↗

Alternating Direction Decomposition with Strong Bounding and Convexification (ADDSBC) for Solving Security Constrained AC Unit Commitment Problems

This project aims to develop efficient and robust computational methods for solving the security-constrained unit commitment and alternating current optimal power flow problem (SC-UC-ACOPF). The SC-UC-ACOPF problem is at the center of the short-term operation of the U.S. Power Grid. It is solved every week, every day, and every 10 minutes to plan for the optimal action of electricity generation and consumption by minimizing the generation cost and maintaining power system reliability against potential disruptions of equipment failures. In mathematical terms, SC-UC-ACOPF is a challenging large-scale mixed-integer nonlinear optimization model. This means that the decisions involve both discrete variables, e.g. the turning on and off of generators and switching of transmission lines and transformers, and continuous decisions, e.g. the amount of energy generated by each generator and the power flows in the power grid. The physics of the power flow is described by nonlinear equations involving real and reactive power and bus voltages. Another key feature is the large number of contingencies, i.e. the system needs to stay reliable in face of failure of any one equipment, such as transmission lines and generators. The U.S. power grids are extremely complicated and large scale with more than 5,000 generators, 50,000 buses, and 100,000 high-voltage transmission lines, making the SC-UC-ACOPF a very large-scale computation challenge. The research developed in this project aims to solve the SC-UC-ACOPF problems in the three timescales, i.e. weekly, daily, and every 10-min. The proposed computational methods are built on a principled algorithmic approach of decomposition and penalization. More specifically, the algorithm develops spatial and temporal decomposition by exploiting the strong temporal coupling and weak spatial coupling of the UC problem and the complementary feature, i.e. weak temporal coupling and strong spatial coupling of the ACOPF problem. The algorithm also leverages recent progresses in strong convex relaxation of ACOPF. A unique feature of the proposed approach is that it generates a valid, global upper bound on the optimal maximum profit. In this way, a global optimality gap is available to measure the quality of the solution. To further speed up computation, the research team has developed a plethora of effective heuristics to strengthen the iterative penalty-based decomposition framework. For instance, a heuristic is developed to construct inner approximations of the time coupling constraints within the time decoupled problems. Contingencies are pre-screened and low-rank matrix computation is exploited to find the almost unique solution to each contingency. A novel heuristic for line switching is proposed and tested with positive impacts on instances where line switching is beneficial. Taking a systematic approach and carefully handling every detail of the problem pays off. The TIM-GO’s performance throughout the trials and the final event was stellar. TIM-GO garnered the second highest total prize money and is ranked in the top three positions across all categories of comparison.

97 MATHEMATICS AND COMPUTING↗

Photovoltaic Cable Jackets: A Comparison of Representative Products Using Combined-Accelerated Stress Testing [Slides]

Photovoltaic (PV) cables facilitate the distribution of electricity collected from modules to the energy grid. Durable cabling enables continuous operation of PV installations, whereas cables with a lifetime less than the modules must be replaced - reducing electricity generation and adding to the operating expense. This study primarily focusses on the aging of the key cable types using the combined-accelerated stress testing (C-AST) protocol. Representative cables for utility, building, and control/auxiliary applications were examined. Cable jacket materials examined include: polyolefin, polyethylene, polyamide, poly(vinyl chloride), chlorinated polyethylene, thermoplastic elastomer, and ethylene propylene diene monomer rubber. Specimen characterizations applied include: optical microscopy, mechanical profilometry, instrumented indentation, scanning electron microscopy (SEM), energy-dispersive X-ray spectroscopy (EDS), Fourier-transform infrared spectroscopy (FTIR), thermogravimetric analysis (TGA), and differential scanning calorimetry (DSC). A variety of performance and durability characteristics were observed, depending on the base material, polymer formulation, and jacket color. The results from C-AST are analyzed and discussed relative to a recent industry survey on electronic balance of system components in addition to a recent study where similar cables were aged using steady state ultraviolet weathering (International Electrotechnical Commission Technical Specification 62788-7-2). Recommendations are made for the screening, industry qualification, and service life prediction of PV cable jackets.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Making Wireless Power Matter (CRADA Final Report)

The Resonant Link project team investigated the integration of their novel, self-resonant magnetic structure for electromagnetic coils into inductive wireless power transfer systems. Relative to conventional wireless power coils, their multi-layer self-resonant structure (the “MSRS”) reduces losses by up to 6x and lowers to BOM cost by up to 10x. These power-dense, cost-effective coils are a promising path towards the first value-add wireless power systems in a variety of applications, from phones to medical devices to kitchen appliances to mobile robots to electric vehicles. The researchers explored development of the power electronics for wireless power transfer systems using the multi-layer self-resonant structure (MSRS) concept that efficiently transmits power over air gaps for in-home applications, including through-wall power transfer, appliances built around the Ki and Qi standards, and other potential use cases This project specifically focused on the power electronics design, development, and testing for these applications. Resonant Link’s proprietary wireless power technology holds the promise of matching the performance of wired chargers in size and efficiency, especially when operated from a DC distribution grid.

42 ENGINEERING↗

Energy-Flow-Environment Linkage Map

Understanding how flexibility in environmental requirements can facilitate co-optimization of hydropower production outcomes and environmental outcomes is critical for future grid decision-making and operations as renewable energy resources increase. The environmental and power system outcomes connectivity linkage maps presented here provide a framework for conceptually and quantitatively linking power system outcomes to environmental outcomes through hydropower flow decisions. The Executive Summary Map serves as a starting for exploring the links between hydropower system performance outcomes and environmental outcomes. The centralized topic is “Flow from Hydropower System,” and connects the hydropower operations through “Flow through turbines” and “Non-turbine flows”, and to environmental outcomes through Reservoir elevation” and “Flow downstream of the hydropower system”. To the left of these central topics, “Hydro-mechanical operations” are linked through “Hydro-electrical operations” to “Hydropower performance outcomes” (“Reliability”, “Resilience”, “Revenue”, “Emissions”). On the right, environmental outcomes are grouped together by their physical location: “Upstream Outcomes” (“Upstream geomorphology”, “Upstream recreation”, “Upstream habitat”, “Upstream biota and biodiversity”, “Upstream water quality and greenhouse gas”), Outcomes relevant to both “Upstream/downstream or dam interface” (“Navigation”, “Dam safety and maintenance", “Human health”, “Water supply”, “Flood control”, ”Fish passage”), and “Downstream outcomes“ (e.g., “Downstream geomorphology”, “Downstream recreations”, “Downstream habitat”, “Downstream biota and biodiversity”, “Downstream water quality and greenhouse gas”). Each of these subtopics (e.g., “Hydro-mechanical operations”, “Hydro-electrical operations”, “Upstream geomorphology”, “Upstream recreation”) is further explored through their corresponding submaps. The “Read Me” file provides more detail information on map navigation. The “Models and tools database” file provides detailed information of models and tools presented in the maps.

13 HYDRO ENERGY↗

The BANDIT computer program for the reduction of matrix bandwidth for NASTRAN

A computer program, BANDIT, is described which was developed for use as a preprocessor to NASTRAN to automatically resequence the grid point numbers for reduced matrix bandwidth. The BANDIT program accepts a standard NASTRAN data deck as input, resequences the grid point numbers for reduced bandwidth, and generates a set of SEQGP cards for insertion into the NASTRAN bulk data deck. BANDIT uses the Cuthill-McKee strategy for resequencing grid points. The program is operational on the CDC 6000 series of computers, the IBM 360, and the Univac 1108.

Everstine, G. C.↗

Spatial land-use inventory, modeling, and projection/Denver metropolitan area, with inputs from existing maps, airphotos, and LANDSAT imagery

A landscape model was constructed with 34 land-use, physiographic, socioeconomic, and transportation maps. A simple Markov land-use trend model was constructed from observed rates of change and nonchange from photointerpreted 1963 and 1970 airphotos. Seven multivariate land-use projection models predicting 1970 spatial land-use changes achieved accuracies from 42 to 57 percent. A final modeling strategy was designed, which combines both Markov trend and multivariate spatial projection processes. Landsat-1 image preprocessing included geometric rectification/resampling, spectral-band, and band/insolation ratioing operations. A new, systematic grid-sampled point training-set approach proved to be useful when tested on the four orginal MSS bands, ten image bands and ratios, and all 48 image and map variables (less land use). Ten variable accuracy was raised over 15 percentage points from 38.4 to 53.9 percent, with the use of the 31 ancillary variables. A land-use classification map was produced with an optimal ten-channel subset of four image bands and six ancillary map variables. Point-by-point verification of 331,776 points against a 1972/1973 U.S. Geological Survey (UGSG) land-use map prepared with airphotos and the same classification scheme showed average first-, second-, and third-order accuracies of 76.3, 58.4, and 33.0 percent, respectively.

Tom, C.↗

Performance characteristics of ring-cusp thrusters with xenon propellant

The performance characteristics and operating envelope of several 30-cm ring-cusp ion thrusters with xenon propellant were investigated. Results indicate a strong performance dependence on the discharge chamber boundary magnetic fields and resultant distribution of electron currents. Significant improvements in discharge performance over J-series divergent-field thrusters were achieved for large throttling ranges, which translate into reduced cathode emission currents and reduced power dissipation which should be of significant benefit for operation at thruster power levels in excess of 10 kW. Mass spectrometry of the ion beam was documented for both the ring-cusp and J-series thrusters with xenon propellant for determination of overall thruster efficiency, and lifetime. Based on the lower centerline values of doubly charged ions in the ion beam and the lower operating discharge voltage, the screen grid erosion rate of the ring-cusp thruster is expected to be lower than the divergent-field J-series thruster by a factor of 2.

Patterson, M. J.↗

A linearized Euler analysis of unsteady flows in turbomachinery

A method for calculating unsteady flows in cascades is presented. The model, which is based on the linearized unsteady Euler equations, accounts for blade loading shock motion, wake motion, and blade geometry. The mean flow through the cascade is determined by solving the full nonlinear Euler equations. Assuming the unsteadiness in the flow is small, then the Euler equations are linearized about the mean flow to obtain a set of linear variable coefficient equations which describe the small amplitude, harmonic motion of the flow. These equations are discretized on a computational grid via a finite volume operator and solved directly subject to an appropriate set of linearized boundary conditions. The steady flow, which is calculated prior to the unsteady flow, is found via a Newton iteration procedure. An important feature of the analysis is the use of shock fitting to model steady and unsteady shocks. Use of the Euler equations with the unsteady Rankine-Hugoniot shock jump conditions correctly models the generation of steady and unsteady entropy and vorticity at shocks. In particular, the low frequency shock displacement is correctly predicted. Results of this method are presented for a variety of test cases. Predicted unsteady transonic flows in channels are compared to full nonlinear Euler solutions obtained using time-accurate, time-marching methods. The agreement between the two methods is excellent for small to moderate levels of flow unsteadiness. The method is also used to predict unsteady flows in cascades due to blade motion (flutter problem) and incoming disturbances (gust response problem).

Hall, Kenneth C.↗

Performance characteristics of ring-cusp thrusters with xenon propellant

The performance characteristics and operating envelope of several 30-cm ring-cusp ion thrusters with xenon propellant were investigated. Results indicate a strong performance dependence on the discharge chamber boundary magnetic fields and resultant distribution of electron currents. Significant improvements in discharge performance over J-series divergent-field thrusters were achieved for large throttling ranges, which translate into reduced cathode emission currents and reduced power dissipation which should be of significant benefit for operation at thruster power levels in excess of 10 kW. Mass spectrometer of the ion beam was documented for both the ring-cusp and J-series thrusters with xenon propellant for determination of overall thruster efficiency, and lifetime. Based on the lower centerline values of doubly charged ions in the ion beam and the lower operating discharge voltage, the screen grid erosion rate of the ring-cusp thruster is expected to be lower than the divergent-field J-series thruster by a factor of 2.

Patterson, M. J.↗

Calculation of unsteady flows in turbomachinery using the linearized Euler equations

A method for calculating unsteady flows in cascades is presented. The model, which is based on the linearized unsteady Euler equations, accounts for blade loading shock motion, wake motion, and blade geometry. The mean flow through the cascade is determined by solving the full nonlinear Euler equations. Assuming the unsteadiness in the flow is small, then the Euler equations are linearized about the mean flow to obtain a set of linear variable coefficient equations which describe the small amplitude, harmonic motion of the flow. These equations are discretized on a computational grid via a finite volume operator and solved directly subject to an appropriate set of linearized boundary conditions. The steady flow, which is calculated prior to the unsteady flow, is found via a Newton iteration procedure. An important feature of the analysis is the use of shock fitting to model steady and unsteady shocks. Use of the Euler equations with the unsteady Rankine-Hugoniot shock jump conditions correctly models the generation of steady and unsteady entropy and vorticity at shocks. In particular, the low frequency shock displacement is correctly predicted. Results of this method are presented for a variety of test cases. Predicted unsteady transonic flows in channels are compared to full nonlinear Euler solutions obtained using time-accurate, time-marching methods. The agreement between the two methods is excellent for small to moderate levels of flow unsteadiness. The method is also used to predict unsteady flows in cascades due to blade motion (flutter problem) and incoming disturbances (gust response problem).

Hall, Kenneth C.↗

The large discretization step method for time-dependent partial differential equations

A new method for the acceleration of linear and nonlinear time dependent calculations is presented. It is based on the Large Discretization Step (LDS) approximation, defined in this work, which employs an extended system of low accuracy schemes to approximate a high accuracy discrete approximation to a time dependent differential operator. Error bounds on such approximations are derived. These approximations are efficiently implemented in the LDS methods for linear and nonlinear hyperbolic equations, presented here. In these algorithms the high and low accuracy schemes are interpreted as the same discretization of a time dependent operator on fine and coarse grids, respectively. Thus, a system of correction terms and corresponding equations are derived and solved on the coarse grid to yield the fine grid accuracy. These terms are initialized by visiting the fine grid once in many coarse grid time steps. The resulting methods are very general, simple to implement and may be used to accelerate many existing time marching schemes.

Haras, Zigo↗

High-order Two-way Artificial Boundary Conditions for Nonlinear Wave Propagation with Backscattering

When solving linear scattering problems, one typically first solves for the impinging wave in the absence of obstacles. Then, by linear superposition, the original problem is reduced to one that involves only the scattered waves driven by the values of the impinging field at the surface of the obstacles. In addition, when the original domain is unbounded, special artificial boundary conditions (ABCs) that would guarantee the reflectionless propagation of waves have to be set at the outer boundary of the finite computational domain. The situation becomes conceptually different when the propagation equation is nonlinear. In this case the impinging and scattered waves can no longer be separated, and the problem has to be solved in its entirety. In particular, the boundary on which the incoming field values are prescribed, should transmit the given incoming waves in one direction and simultaneously be transparent to all the outgoing waves that travel in the opposite direction. We call this type of boundary conditions two-way ABCs. In the paper, we construct the two-way ABCs for the nonlinear Helmholtz equation that models the laser beam propagation in a medium with nonlinear index of refraction. In this case, the forward propagation is accompanied by backscattering, i.e., generation of waves in the direction opposite to that of the incoming signal. Our two-way ABCs generate no reflection of the backscattered waves and at the same time impose the correct values of the incoming wave. The ABCs are obtained for a fourth-order accurate discretization to the Helmholtz operator; the fourth-order grid convergence is corroborated experimentally by solving linear model problems. We also present solutions in the nonlinear case using the two-way ABC which, unlike the traditional Dirichlet boundary condition, allows for direct calculation of the magnitude of backscattering.

Fibich, Gadi↗

ICEG2D (v2.0) - An Integrated Software Package for Automated Prediction of Flow Fields for Single-Element Airfoils With Ice Accretion

An integrated geometry/grid/simulation software package, ICEG2D, is being developed to automate computational fluid dynamics (CFD) simulations for single- and multi-element airfoils with ice accretions. The current version, ICEG213 (v2.0), was designed to automatically perform four primary functions: (1) generate a grid-ready surface definition based on the geometrical characteristics of the iced airfoil surface, (2) generate high-quality structured and generalized grids starting from a defined surface definition, (3) generate the input and restart files needed to run the structured grid CFD solver NPARC or the generalized grid CFD solver HYBFL2D, and (4) using the flow solutions, generate solution-adaptive grids. ICEG2D (v2.0) can be operated in either a batch mode using a script file or in an interactive mode by entering directives from a command line within a Unix shell. This report summarizes activities completed in the first two years of a three-year research and development program to address automation issues related to CFD simulations for airfoils with ice accretions. As well as describing the technology employed in the software, this document serves as a users manual providing installation and operating instructions. An evaluation of the software is also presented.

Thompson David S.↗

Translation Optics for 30 cm Ion Engine Thrust Vector Control

Data were obtained from a 30 cm xenon ion thruster in which the accelerator grid was translated in the radial plane. The thruster was operated at three different throttle power levels, and the accelerator grid was incrementally translated in the X, Y, and azimuthal directions. Plume data was obtained downstream from the thruster using a Faraday probe mounted to a positioning system. Successive probe sweeps revealed variations in the plume direction. Thruster perveance, electron backstreaming limit, accelerator current, and plume deflection angle were taken at each power level, and for each accelerator grid position. Results showed that the thruster plume could easily be deflected up to six degrees without a prohibitive increase in accelerator impingement current. Results were similar in both X and Y direction.

Haag, Thomas↗

Study of Near-Surface Models in Large-Eddy Simulations of a Neutrally Stratified Atmospheric Boundary Layer

Large-eddy simulation (LES) is a widely used technique in armospheric modeling research. In LES, large, unsteady, three dimensional structures are resolved and small structures that are not resolved on the computational grid are modeled. A filtering operation is applied to distinguish between resolved and unresolved scales. We present two near-surface models that have found use in atmospheric modeling. We also suggest a simpler eddy viscosity model that adopts Prandtl's mixing length model (Prandtl 1925) in the vicinity of the surface and blends with the dynamic Smagotinsky model (Germano et al, 1991) away from the surface. We evaluate the performance of these surface models by simulating a neutraly stratified atmospheric boundary layer.

Senocak, I.↗