Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Grid disturbance”

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 OSTI2026

Model Data Archive Associated with Manuscript "Fire-altered Carbon Pools Create Disturbance Memory in Stream Dissolved Organic Carbon"

This data package supports the publication “Fire-altered Carbon Pools Create Disturbance Memory in Stream Dissolved Organic Carbon” by Li et al. (2026). The package contains processed model inputs, configuration files, restart files, simulation outputs, scripts, and visualization products used to evaluate post-fire dissolved organic carbon (DOC) dynamics in the Naches River Watershed, Washington, USA, following the 2021 Schneider Springs Fire. The modeling workflow couples ELM-BGC, the biogeochemistry-enabled Energy Exascale Earth System Model Land Model; ATS, the Advanced Terrestrial Simulator for integrated surface-subsurface hydrology; and PFLOTRAN, a reactive transport model for multicomponent aqueous geochemistry. Together, these models simulate how wildfire-induced changes in vegetation, litter, coarse woody debris, and soil organic matter influence DOC production, transport, and reaction from burned hillslopes to stream networks. The archive includes preprocessed meteorological, geospatial, hydrologic, and biogeochemical forcing data; ELM-BGC-derived DOC source terms; ATS mesh files; PFLOTRAN reactive-transport inputs; model configuration files; spin-up and transient restart files; watershed-scale diagnostic outputs; stream concentration time series; and figures or visualization files used to inspect and reproduce key results. File types include Hierarchical Data Format 5 (HDF5) files for gridded forcing and model-coupling data, model input and configuration files for ELM-BGC, ATS, and PFLOTRAN, restart and simulation-output files generated by the modeling workflow, tabular or time-series diagnostic outputs, scripts for post-processing and figure generation, and image or visualization products associated with the manuscript. Use of the package depends on the intended task. Re-running the simulations requires the relevant modeling software, including ELM-BGC, ATS, and PFLOTRAN as ATS's geochemical engine. Inspecting outputs and reproducing figures requires Python with scientific plotting libraries such as Matplotlib, and three-dimensional model outputs may be viewed with ParaView. Geographic information system files or maps may be inspected with ArcGIS Pro or comparable GIS software. The data package is intended to enable traceability, reuse, and partial reproduction of the coupled land-to-watershed hydro-biogeochemical modeling workflow used to test how wildfire disturbance affects terrestrial carbon pools and downstream DOC dynamics.

ATS↗
NASA NTRS1993

Analysis of rocket beacon transmissions for computerized reconstruction of ionospheric densities

Three methods are described to obtain ionospheric electron densities from transionospheric, rocket-beacon TEC data. First, when the line-of-sight from a ground receiver to the rocket beacon is tangent to the flight trajectory, the electron concentration can be obtained by differentiating the TEC with respect to the distance to the rocket. A similar method may be used to obtain the electron-density profile if the layer is horizontally stratified. Second, TEC data obtained during chemical release experiments may be interpreted with the aid of physical models of the disturbed ionosphere to yield spatial maps of the modified regions. Third, computerized tomography (CT) can be used to analyze TEC data obtained along a chain of ground-based receivers aligned along the plane of the rocket trajectory. CT analysis of TEC data is used to reconstruct a 2D image of a simulated equatorial plume. TEC data is computed for a linear chain of nine receivers with adjacent spacings of either 100 or 200 km. The simulation data are analyzed to provide an F region reconstruction on a grid with 15 x 15 km pixels. Ionospheric rocket tomography may also be applied to rocket-assisted measurements of amplitude and phase scintillations and airglow intensities.

Bernhardt, P. A.↗
DOE OSTI2024

Modeling of Induction Motors and Variable Speed Drives for Multi-Domain System Simulations Using Modelica and the OpenIPSL Library

This paper introduces an innovative method for characterizing, implementing, and validating both three-phase and single-phase induction motor models, accompanied by a variable speed drive model. The primary goal is to investigate interactions between the electrical power grid and other dynamic domains (e.g., thermofluidic) that impact motor/load drive behavior. Our approach involves establishing a mechanical interface based on a physically meaningful equation linking motor torque/speed to the electrical model in the phasor domain. This allows seamless integration of diverse domain subsystems into a unified multi-domain model using Modelica v4.0.0 and the OpenIPSL library v3.0.1, overcoming co-simulation limitations. The proposed model, which requires only one Modelica-compliant tool for simulation, introduces additional dynamics through the mechanical interface, enabling explicit simulation of load disturbances based on constitutive physics. This deepens our understanding of dynamic interactions between the electrical power domain and other subsystems connected through the motor. We detail the modeled components using mathematical equations and textual descriptions, emphasizing the Modelica modeling approach. Simulation examples validate the implementation, demonstrating the multi-domain modeling capabilities of the newly developed components.

Computer Science↗
NASA NTRS2013

A Major Solar Eruptive Event in July 2012: Defining Extreme Space Weather Scenarios

A key goal for space weather studies is to define severe and extreme conditions that might plausibly afflict human technology. On 23 July 2012, solar active region 1520 (approx. 141 deg W heliographic longitude) gave rise to a powerful coronal mass ejection (CME) with an initial speed that was determined to be 2500+/-500 km/s. The eruption was directed away from Earth toward 125 deg W longitude. STEREO-A sensors detected the CME arrival only about 19 h later and made in situ measurements of the solar wind and interplanetary magnetic field. In this paper, we address the question of what would have happened if this powerful interplanetary event had been Earthward directed. Using a well-proven geomagnetic storm forecast model, we find that the 23-24 July event would certainly have produced a geomagnetic storm that was comparable to the largest events of the twentieth century (Dst approx. - 500 nT). Using plausible assumptions about seasonal and time-of-day orientation of the Earth's magnetic dipole, the most extreme modeled value of storm-time disturbance would have been Dst= - 1182 nT. This is considerably larger than estimates for the famous Carrington storm of 1859. This finding has far reaching implications because it demonstrates that extreme space weather conditions such as those during March of 1989 or September of 1859 can happen even during a modest solar activity cycle such as the one presently underway. We argue that this extreme event should immediately be employed by the space weather community to model severe space weather effects on technological systems such as the electric power grid.

solar longitude↗
DOE OSTI2026

Post-Event Fault Identification with Machine Learning for Protection System Validation

Power system protection devices have transitioned over the past few decades from mechanical to analog devices, then to solid state and finally digital. Relays and their associated critical network of equipment have significantly increased in complexity. Even internally, relays have gained significant intricacy, with relatively simple overcurrent or differential functions now being assisted by a myriad of other functions. This is necessary as the grid becomes more complex, but it brings increased difficulty in monitoring and upkeep. Misoperation caused by accidental improper relay settings or deliberate malicious actions is a constant challenge faced by all utilities. These improper settings can be difficult to identify and may require exhaustive post-mortem analysis, typically after a major outage event has already occurred. A mechanism is needed for monitoring the behavior of protection systems to validate that their performance falls within expectations. Relays that fail to isolate a fault or trip when there is no system disturbance can be flagged for settings review in situations where this behavior may not have been noticed due to manual restoration or backup protection operations. This work presents a concept for a machine learning (ML) system capable of validating the performance of protection systems by identifying fault events and characterizing protection system responses based solely on available current and voltage measurements. As a first step in its development, an experimental dataset is generated, and a random forest model is implemented with high accuracy in distinguishing four power system scenarios.

24 - POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Electrical Systems of Pumped Storage Hydropower Plants: Electrical Generation, Machines, Power Electronics, and Power Systems

Adjustable-speed (AS) pumped storage hydropower (PSH) technology has the potential to become a large, consistent contributor to grid stability enabling higher penetrations of wind and solar energy on the future U.S. power system. AS-PSH has high-value characteristics, such as fast response to provide ancillary services to the grid, because it is a power converter interface with the grid (like battery storage), but at the same time it has the energy content large enough to supply both short-term (seconds-to-minutes) and long-term (minutes-to-hours) of energy needs, like more conventional power plants. However, designs must be optimized to lower the capital expenditure (CapEx) and to provide a high-quality grid interface capability (e.g. power quality, ancillary service provider, fault-tolerant or fault ride-through capability), which is a primary factor in the acceptance of AS-PSH into a utility’s generation mix. This CapEx will be greatly affected by the cost savings associated with the civil structure, turbine design, power electronics, control systems, or unique generator designs. A holistic design must be considered to get a full picture of the benefits of the technology proposed. The AS-PSH can be controlled to reduce the impact of transient disturbance on a power system and at the same time can be controlled to minimized subsequent component fatigue and potential oscillation modes within the plant, with the overall impact in reducing the operational expenditures (OpEx).Generating clean power to meet standards such as Institute of Electrical and Electronics Engineers (IEEE) 519 and International Electrotechnical Commission (IEC) 1000-3-2 will be a continuing challenge. For many technology developers, however, improved AS-PSH technologies will become a key component of generator-storage systems of the future given the prospects of increased performance and decreasing costs, and the ever-increasing penetration of renewables (e.g wind power and solar power).

13 HYDRO ENERGY↗
NASA NTRS2019

Importance of Tree- and Species-Level Interactions With Wildfire, Climate, and Soils in Interior Alaska: Implications for Forest Change Under A Warming Climate

The boreal zone of Alaska is dominated by interactions between disturbances, vegetation, and soils. These interactions are likely to change in the future through increasing permafrost thaw, more frequent and intense wildfires, and vegetation change from drought and competition. We utilize an individual tree-based vegetation model, the University of Virginia Forest Model Enhanced (UVAFME), to estimate current and future forest conditions across sites within interior Alaska. We updated UVAFME for application within interior Alaska, including improved simulation of permafrost dynamics, litter decay, nutrient dynamics, fire mortality, and post-fire regrowth. Following these updates, UVAFME output on species-specific biomass and stem density was comparable to inventory measurements at various forest types within interior Alaska. We then simulated forest response to climate change at specific inventory locations and across the Tanana Valley River Basin on a 2 × 2 km2 grid. We derived projected temperature and precipitation from a five-model average taken from the CMIP5 archive under the RCP 4.5 and 8.5 scenarios. Results suggest that climate change and the concomitant impacts on wildfire and permafrost dynamics will result in overall decreases in biomass (particularly for spruce (Picea spp.)) within the interior Tanana Valley, despite increases in quaking aspen (Populus tremuloides) biomass, and a resulting shift towards higher deciduous fraction. Simulation results also predict increases in biomass at cold, wet locations and at high elevations, and decreases in biomass in dry locations, under both moderate (RCP 4.5) and extreme (RCP 8.5) climate change scenarios. These simulations demonstrate that a highly detailed, species interactive model can be used across a large region within Alaska to investigate interactions between vegetation, climate, wildfire, and permafrost. The vegetation changes predicted here have the capacity to feed back to broader scale climate-forest interactions in the North American boreal forest, a region which contributes significantly to the global carbon and energy budgets.

forest modeling↗
DOE OSTI2025

InterGraph-CPS: A Graph-Theoretic Approach to Characterize Cross-Domain Cyber-Physical Interdependencies and Uncertainties in Electric Grid Systems for Improved Decision-Making in Operation and Response

Critical infrastructure systems such as the electric grid are increasingly cyber-physical; yet, despite the cyber-physical characteristics of critical infrastructure systems, the physical process system and communication/control network system are traditionally analyzed in siloes. As these systems become more cyber-physical, it is crucial that models and methods are available to assess the cyber physical system (CPS) interdependencies, characteristics, and event propagation for improved planning, operation, and response. Thus, we proposed an integrated structural and temporal CPS interdependency analysis framework, InterGraph-CPS, that provides insight into the CPS function during normal operation as well as disturbances. This integrated structural and temporal interdependency framework is uniquely designed for assessing CPSs by account for the challenges of analyzing cyber and physical data streams together due to data availability, data type, and time scale differences. By leveraging both structural (e.g., graph analysis) and temporal (e.g., data analytics) techniques, different CPS behaviors and configurations can be accounted for.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

A Cooperative Game Theory-based Approach to Under-frequency Load Shedding Control

This paper proposes a cooperative game theory-based under-frequency load shedding (UFLS) approach for frequency stability and control in power systems. UFLS is a crucial factor for frequency stability and control especially in power grids with high penetration of renewable energy sources and restructured power systems. Conventional UFLS methods, most of which are off-line, usually shed fixed amounts of predetermined loads based on a predetermined schedule which can lead to over or under curtailment of load. This paper presents a co-operative game theory-based two-stage strategy to effectively and precisely determine locations and amounts of loads to be shed for UFLS control. In the first stage, the total amount of loads to be shed, also referred to as deficit in generation or the disturbance power, is computed using the initial rate of change of frequency (ROCOF) referred to the equivalent inertial center. In the second stage, the Shapley value, one of the solution concepts of cooperative game theory, is used to determine load shedding amounts and locations. The proposed method is implemented on the reduced 9-bus 3-machine Western Electricity Coordinating Council (WECC) system and simulated on Real-time Digital Simulators (RTDS). The results show that the proposed UFLS approach can effectively return the system to normal state after disturbances.

42 ENGINEERING↗
NASA NTRS2001

A Cloud-Resolving Simulation of Hurricane Bob (1991): Storm Structure and Eyewall Buoyancy

A numerical simulation of Hurricane Bob (1991) is conducted using the Penn State University-National Center for Atmospheric Research mesoscale model MM5 with a horizontal grid spacing of 1.3 Km on the finest nested mesh The model produces a realistic hurricane that intensifies slowly during the period of fine-scale simulation. Time-averaged results reveal the effects of storm motion. vertical shear, beta gyres and deformation forcing on the structure of radial inflow, vertical motion, and precipitation. Instantaneous model fields show that radial inflow in the eyewall is very intense near the surface but transitions to strong low-level outflow near the top of the boundary layer. The low-level structure is modulated by a wavenumber 2 disturbance that rotates around the eyewall at half the speed of the maximum tangential winds and is consistent with a vortex Rossby edge wave. The statistical distribution of vertical velocity in the eyewall indicates that the eyewall is composed of a small number of intense updrafts that account for the majority of the upward mass flux rather than a more gradual and symmetric eyewall circulation, consistent with the concept of hot towers. Tongues of high equivalent potential temperature, Theta(sub e), are seen along the inner edge of the eyewall updraft and within the low-level outflow. This air originates from outside of the eyewall with the highest theta(sub e) air coming from the layer closest to the surface after penetrating closest to the center. Occasionally, high Theta(sub e), air within the eye is drawn into the eyewall updrafts. The high Theta(sub e), air rising within the eyewall is shown to be associated with positive eyewall buoyancy with sufficient convective available potential energy along its path to produce relatively strong convective updrafts. Although the requirements for conditional symmetric instability are met within the eyewall and the air parcel trajectories follow slanted paths, the radial displacement of air parcels in the low-level outflow moves the air parcel sufficiently far away from the upper- warm core that the air becomes unstable to vertical displacements. Hence, convective instability rather than symmetric instability accounts for the stronger updrafts in the eyewall.

Braun, Scott A.↗
DOE OSTI2024

Passive Islanding Detection of Inverter-Based Resources in a Noisy Environment

Islanding occurs when a load is energized solely by local generators and can result in frequency and voltage instability, changes in current, and poor power quality. Poor power quality can interrupt industrial operations, damage sensitive electrical equipment, and induce outages upon the resynchronization of the island with the grid. This study proposes an islanding detection method employing a Duffing oscillator to analyze voltage fluctuations at the point of common coupling (PCC) under a high-noise environment. Unlike existing methods, which overlook the noise effect, this paper mitigates noise impact on islanding detection. Power system noise in PCC measurements arises from switching transients, harmonics, grounding issues, voltage sags and swells, electromagnetic interference, and power quality issues that affect islanding detection. Transient events like lightning-induced traveling waves to the PCC can also introduce noise levels exceeding the voltage amplitude by more than seven times, thus disturbing conventional detection techniques. The noise interferes with measurements and increases the nondetection zone (NDZ), causing failed or delayed islanding detection. The Duffing oscillator nonlinear dynamics enable detection capabilities at a high noise level. The proposed method is designed to detect the PCC voltage fluctuations based on the IEEE standard 1547 through the Duffing oscillator. For the voltages beyond the threshold, the Duffing oscillator phase trajectory changes from periodic to chaotic mode and sends an islanded operation command to the inverter. The proposed islanding detection method distinguishes switching transients and faults from an islanded operation. Experimental validation of the method is conducted using a 3.6 kW PV setup.

Energy & Fuels↗
NASA NTRS2000

Single Airfoil Gust Response Problem

An unsteady aerodynamic code, called GUST3D, has been developed to solve the single, non-constant coefficient, inhomogeneous convective wave equation for flows with periodic vortical disturbance. The code uses a frequency-domain approach with second-order central differences and a Sommerfeld radiation condition in the far field. It has been extensively validated on model problems with analytical solutions. GUST3D requires as input certain mean flow quantities which are calculated separately by a potential flow solver. This solver calculates the mean flow using a Gothert's Rule approximation. On the airfoil surface, it uses the solution calculated by the potential code FL036. A figure shows the mean pressure along the airfoil surface. To calculate the unsteady pressure, GUST3D was run oil systematically refined grids to obtain a converged solution at each frequency. It was found that 24 points per wavelength was sufficient for convergence. The location of the outer grid boundary was also varied to check for sensitivity to the far-field boundary condition.

Scott, James R.↗
NASA NTRS1984

An efficient coordinate transformation technique for unsteady, transonic aerodynamic analysis of low aspect-ratio wings

An efficient coordinate transformation technique is presented for constructing grids for unsteady, transonic aerodynamic computations for delta-type wings. The original shearing transformation yielded computations that were numerically unstable and this paper discusses the sources of those instabilities. The new shearing transformation yields computations that are stable, fast, and accurate. Comparisons of those two methods are shown for the flow over the F5 wing that demonstrate the new stability. Also, comparisons are made with experimental data that demonstrate the accuracy of the new method. The computations were made by using a time-accurate, finite-difference, alternating-direction-implicit (ADI) algorithm for the transonic small-disturbance potential equation.

Guruswamy, G. P.↗
DOE OSTI2024

Adaptive Online Model Update Algorithm for Predictive Control in Networked Systems

In this article, we introduce an adaptive on-line model update algorithm designed for predictive control applications in networked systems, particularly focusing on power distribution systems. Unlike traditional methods that depend on historical data for offline model identification, our approach utilizes real-time data for continuous model updates. This method integrates seamlessly with existing online control and optimization algorithms and provides timely updates in response to real-time changes. This methodology offers significant advantages, including a reduction in the communication network bandwidth requirements by minimizing the data exchanged at each iteration and enabling the model to adapt after disturbances. Furthermore, our algorithm is tailored for non-linear convex models, enhancing its applicability to practical scenarios. The efficacy of the proposed method is validated through a numerical study, demonstrating improved control performance using a synthetic IEEE test case.

data-driven model predictive control↗
NASA NTRS2004

Single Airfoil Gust Response Problem: Category 3, Problem 1

An unsteady aerodynamic code, called GUST3D (ref. 3), has been developed to solve equation (8) for flows with periodic vortical disturbances. The code uses a frequency-domain approach with second-order central differences and a pressure radiation condition in the far field. GUST3D requires as input certain mean flow quantities which are calculated separately by a potential flow solver. The solver calculates the mean ow using a Gothert's Rule approximation (ref. 3). On the airfoil surface, it uses the solution calculated by the potential code FLO36 (ref. 4). Figures 1-2 show the mean pressure along the airfoil surface for the two airfoil geometries. In Figures 3 - 8, we present the RMS pressure on the airfoil surface. Each figure shows three GUST3D solutions (calculated on grids with different far-field boundary locations). Three solutions are shown to provide some indication of the numerical uncertainty in the results. Figures 9 - 13 present the acoustic intensity. We again show three solutions per case. Note that no results are presented for the k1 = k2 = 2.0 loaded airfoil case, as an acceptable solution could not be obtained. A few comments need to be made about the results shown. First, since the last Workshop, the GUST3D code has been substantially upgraded. This includes implementing a more accurate far-field boundary condition (ref. 5) and developing improved gridding capabilities. This is the reason for any differences that may exist between the present results and results from the last Workshop. Second, the intensity results on the circle R = 4C were obtained using a Kirchoff method (ref. 6). The Kirchoff surface was the circle R = 2C. Finally, the GUST3D code is most accurate for low reduced frequencies. A new domain decomposition approach (ref. 7) has been developed to improve accuracy. Both the single domain and domain decomposition approaches were used in generating the present results.

Scott, James R.↗
NASA NTRS1992

Validation of three-dimensional incompressible spatial direct numerical simulation code: A comparison with linear stability and parabolic stability equation theories for boundary-layer transition on a flat plate

Spatially evolving instabilities in a boundary layer on a flat plate are computed by direct numerical simulation (DNS) of the incompressible Navier-Stokes equations. In a truncated physical domain, a nonstaggered mesh is used for the grid. A Chebyshev-collocation method is used normal to the wall; finite difference and compact difference methods are used in the streamwise direction; and a Fourier series is used in the spanwise direction. For time stepping, implicit Crank-Nicolson and explicit Runge-Kutta schemes are used to the time-splitting method. The influence-matrix technique is used to solve the pressure equation. At the outflow boundary, the buffer-domain technique is used to prevent convective wave reflection or upstream propagation of information from the boundary. Results of the DNS are compared with those from both linear stability theory (LST) and parabolized stability equation (PSE) theory. Computed disturbance amplitudes and phases are in very good agreement with those of LST (for small inflow disturbance amplitudes). A measure of the sensitivity of the inflow condition is demonstrated with both LST and PSE theory used to approximate inflows. Although the DNS numerics are very different than those of PSE theory, the results are in good agreement. A small discrepancy in the results that does occur is likely a result of the variation in PSE boundary condition treatment in the far field. Finally, a small-amplitude wave triad is forced at the inflow, and simulation results are compared with those of LST. Again, very good agreement is found between DNS and LST results for the 3-D simulations, the implication being that the disturbance amplitudes are sufficiently small that nonlinear interactions are negligible.

Joslin, Ronald D.↗
NASA NTRS1999

Effects of G-Jitter on Directional Solidification of a Binary Alloy

A study of directional solidification of a weak binary alloy (specifically, Bi - 1 at% Sn) based on the fixed grid single domain approach is being undertaken. The enthalpy method is used to solve for the temperature field over the computational domain including both the solid and liquid phases; latent heat evolution is treated with the aid of an effective specific heat coefficient. A source term accounting for the release of solute into the liquid during solidification has been incorporated into the solute transport equation. The vorticity-stream function formulation is used to describe thermo-solutal convection in the liquid region. In this paper we present a numerical simulation of g-jitter. A background gravity of 1 microgram has been assumed, and new results for the effects of periodic disturbances over a range of amplitudes and frequencies on solute field and segregation have been presented.

Benjapiyaporn, C.↗
NASA NTRS1989

Local gravity disturbance estimation from multiple-high-single-low satellite-to-satellite tracking

The idea of satellite-to-satellite tracking in the high-low mode has received renewed attention in light of the uncertain future of NASA's proposed low-low mission, Geopotential Research Mission (GRM). The principal disadvantage with a high-low system is the increased time interval required to obtain global coverage since the intersatellite visibility is often obscured by Earth. The U.S. Air Force has begun to investigate high-low satellite-to-satellite tracking between the Global Positioning System (GPS) of satellites (high component) and NASA's Space Transportation System (STS), the shuttle (low component). Because the GPS satellites form, or will form, a constellation enabling continuous three-dimensional tracking of a low-altitude orbiter, there will be no data gaps due to lack of intervisibility. Furthermore, all three components of the gravitation vector are estimable at altitude, a given grid of which gives a stronger estimate of gravity on Earth's surface than a similar grid of line-of-sight gravitation components. The proposed Air Force mission is STAGE (Shuttle-GPS Tracking for Anomalous Gravitation Estimation) and is designed for local gravity field determinations since the shuttle will likely not achieve polar orbits. The motivation for STAGE was the feasibility to obtain reasonable accuracies with absolutely minimal cost. Instead of simulating drag-free orbits, STAGE uses direct measurements of the nongravitational forces obtained by an inertial package onboard the shuttle. The sort of accuracies that would be achievable from STAGE vis-a-vis other satellite tracking missions such as GRM and European Space Agency's POPSAT-GRM are analyzed.

Jekeli, Christopher↗