Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transient stability assessment”

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 37 records · Page 2

A NASTRAN DMAP alter for linear buckling analysis under dynamic loading

A modification to the NASTRAN solution sequence for transient analysis with direct time integration (COSMIC NASTRAN rigid format 9) was developed and incorporated into a DMAP alter. This DMAP alter calculates the buckling stability of a dynamically loaded structure, and is used to predict the onset of structural buckling under stress-wave loading conditions. The modified solution sequence incorporates the linear buckling analysis capability (rigid format 5) of NASTRAN into the existing Transient solution rigid format in such a way as to provide a time dependent eigensolution which is used to assess the buckling stability of the structure as it responds to the impulsive load. As a demonstration of the validity of this modified solution procedure, the dynamic buckling of a prismatic bar subjected to an impulsive longitudinal compression is analyzed and compared to the known theoretical solution. In addition, a dynamic buckling analysis is performed for the analytically less tractable problem of the localized dynamic buckling of an initially flawed composite laminate under transverse impact loading. The addition of this DMAP alter to the transient solution sequence in NASTRAN facilitates the computational prediction of both the time at which the onset of dynamic buckling occurs in an impulsively loaded structure, and the dynamic buckling mode shapes of that structure.

Aiello, Robert A.↗

Multi-Timescale Integrated Dynamics and Scheduling for Solar (MIDAS-Solar) (Final Technical Report)

Solar photovoltaic (PV) installations have experienced unprecedented growth in the United States. PV will become not only an energy producer but also a necessary provider of ancillary services at multiple timescales. Conventional methods to simulate power systems operations - such as long-term production simulation (which typically considers schedules from hours to minutes by using an optimization framework) and short-term transient studies (which simulate dynamics from seconds to sub-seconds using state variables and differential equations) - are not sufficient for studying the multiple-timescale variation of solar generation and its impact on system reliability. Long-term system economics and short-term system dynamics are highly coupled, particularly when the penetration level of renewable generation is extremely high, because the uncertainty and variability of solar generation will impact both power system steady-state and dynamic performance. This project helps meet and exceed the U.S. Department of Energy Office of Energy Efficiency and Renewable Energy Solar Energy Technologies Office goal of systems integration by directly addressing this stability and reliability challenge for power grid planning and operation. We have developed a temporally comprehensive, closed-loop simulation model, named Multi-timescale Integrated Dynamics and Scheduling (MIDAS), that seamlessly simulates power system operations from economic scheduling (day-ahead to hours) to dynamic response analysis (seconds to sub-seconds). For schedules with very high levels of inverter-based resources (IBRs), up to and including 100%, the stability of grid controls has been evaluated through electromagnetic transient (EMT) simulations and power-hardware-in-the-loop (PHIL) simulations of key transient events at key schedule points. Specifically, MIDAS provides: 1) a closed-loop simulation framework for simulating timescales from economic scheduling to dynamic stability analysis; 2) machine learning-based stability assessment; 3) EMT modeling and analysis for large-scale power systems; 4) MIDAS PHIL test bed. We worked with Hawaii Electric Companies to apply the MIDAS study framework to a Maui grid study. The entire island's transmission system was modeled in detail - from a yearly scheduling model, to a second-level frequency dynamic model, down to a sub-second-scale EMT model to address critical stability issues. The project demonstrated how MIDAS can help system planners and operators assess system reliability and stability while the power grid is marching toward a high-renewable, high-IBR future. In this Maui grid study, we found that 100% instantaneous IBR operation is achievable in EMT simulation and PHIL testing, and grid planners and operators might need new analysis/simulation tools to assess grid reliability and stability in the scheduling stage. MIDAS will bring Maui and other systems closer to 100% clean and stable energy futures. (In this study, we examined transient stability. Other topics necessary for 100% IBR operation, such as protection and resource adequacy, were not examined.)

100% Renewables↗

Critical Clearing Time for Droop-Controlled Grid-Forming Inverters With Circular Current Limiting and Virtual Impedance Current Limiting

The critical clearing time is a crucial parameter for assessing the large-signal stability of grid-forming (GFM) inverters. Variations in current limiters can significantly impact an inverter's transient behavior during large disturbances. Here, this letter introduces a method based on the power-angle relationship to evaluate the fault recovery of a droop-controlled GFM inverter under a three-phase bolted fault scenario, considering both the circular current limiter and the virtual impedance method. High-fidelity simulations validate the accuracy of the proposed formulas in estimating the critical clearing time of the GFM inverter.

42 ENGINEERING↗

Fuel Flow Modeling under Thermal Fluidics Considerations within Operational Domains of Molten Salt Reactors

Molten Salt Reactors (MSRs) offer significant versatility supported by their design flexibility and safety advantages due to high temperatures, design customization capabilities, strong negative reactivity feedback effects, and salt chemistry characteristics. Different control schedules and different mass flow rates in the MSR primary systems impact neutronics and thermal hydraulics coupling. The result is two-fold. First, it leads to adaptability in various applications when the customization of a MSR unit is desired. Second, it calls for the quantification of dynamics and safety characteristics to assess the operational domain and make sure reactor stability and inherent safety characteristics are maintained. This paper provides an in-depth transient analysis of the operational implications of MSRs having heat-generating fuel and heat-transporting salt mixed forming liquid fuel salts contained within primary systems. A coupling model was developed to simulate fuel flow in molten salt reactors. A zero-dimensional reactor kinetics model was used with a one-dimensional heat transfer model in the reactor core. The temperature reactivity feedbacks of the fuel salt resulting from the thermal-hydraulics model were used in the reactor kinetics model to complete the coupled code. The model was applied to simulate the Molten Salt Breeder Reactor (MSBR) response at steady state and transients. The developed model was shown to be a suitable tool for the dynamic analysis of molten salt reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Transient Stability Study of IEEE 2800-2022 Standard Negative-Sequence-Current Compliant GFM IBRs

The primary function of grid-forming (GFM) inverter-based resources (IBRs) is to maintain voltage and fre- quency/phase angle stability; however, during unbalanced grid faults, IEEE Std. 2800-2022 requires IBRs to inject both positive- and negative-sequence reactive current. Enforcing compliance with this standard could compromise the stability of GFM IBRs under such fault conditions. This paper presents a preliminary transient stability study focused on the power angle stability of GFM IBRs operating in compliance with IEEE Std. 2800 during unbalanced faults. The goal is to assess potential stability issues associated with the standard's current injection requirements under fault ride-through conditions. The results show that injecting positive-sequence current with d-axis priority during unbalanced faults better preserves the power angle stability than the q-axis priority while still enabling negative-sequence current injection. These findings suggest that, to ensure stability, GFM IBRs might need to partially comply with IEEE Std. 2800 by prioritizing positive-sequence current injection along the d-axis.

14 SOLAR ENERGY↗

Use of satellite data and modeling to assess the influence of stratospheric processes on the troposphere

Current research focused on the examination of the effects of ozone heating and Newtonian cooling on the linear stability of transient planetary waves. A finite difference numerical model was developed in which the basic state wind, temperature, and ozone distributions were constructed using ground-based and satellite data. Preliminary results indicate that ozone dynamics interaction may play a more important role in the stability and maintenance of planetary waves than previously thought, particularly in summer when the mean solar zenith angle is smallest and thus ozone heating largest. Examination of the finite amplitude interactions among radiation, ozone, and dynamics was carried out. Self-consistent, coupled equations governing the weakly nonlinear interactions are currently being analyzed to provide a better understanding of wave-mean flow interactions and ozone transport in a continuously stratified model of the troposphere-stratosphere coupled system.

Nathan, Terrence↗

High-fidelity multiphysics load following and accidental transient modeling of microreactors using NEAMS tools: Application of NEAMS codes to perform multiphysics modeling analyses of micro-reactor concepts

The feasibility of modeling microreactors using high-fidelity models with the Nuclear Energy Advanced Modeling and Simulation (NEAMS) tools is investigated in this report. Three overarching questions guided this research: can NEAMS tools readily be applied for high-fidelity multiphysics modeling of different types of transients in microreactor designs; how accurate are the results obtained; and are improvements needed in accuracy or user experience of NEAMS tools, especially considering newly developed capabilities? This work builds upon FY-2022 work, and two microreactor concepts considering heat pipe (HP-MR) and gas-cooled (GC-MR) technologies were further analyzed using high-fidelity multiphysics simulations. The NEAMS tools considered and coupled within the MultiApp environment are Griffin for neutronics, BISON for thermo-mechanics, Sockeye for heat pipe modeling (in HP-MR), SAM for 1D Fluid – 3D solid modeling of coolant channels and system modeling of balance of plant components (in GC-MR), and the SWIFT code for hydrogen redistribution in hydride moderator. The Heat Pipe MicroReactor (HP-MR) concept was further analyzed in FY-2023 to demonstrate the stochastic TRISO failure modeling capability in BISON to check operational limits of the TRISO fuel. A new full-core Gas-Cooled MicroReactor (GC-MR) model was developed based on the initial assembly-model used in Y-2022 and used for steady-state and accidental depressurization transient simulations. Accuracy of the simulations performed was assessed through 1) verification analyses completed on the different physics with code-to-code comparison, and 2) validation of the multiphysics simulations based on modeling of the Kilopower Reactor Using Stirling Technology (KRUSTY) experiment. In FY-2023, the mesh and model of KRUSTY was updated to closely match publicly available data, and the neutronic model was verified and validated against experimental control rod worth measurements. The multiphysics model of KRUSTY was developed and used for steady-state analysis and for modeling reactivity insertion transient. The calculated power increase and stabilization agrees well with experimental data following adjustment in fuel thermal expansion coefficient. As an important component of this project, the ANL team gathered experience with a wide range of NEAMS tools: the MOOSE Mesh System, Griffin, BISON, SWIFT, Sockeye, SAM, Workbench, and the MOOSE MultiApp System, and provided assessment of new capabilities. Noteworthy are the user assessment of the “vapor-only” flow model in Sockeye and development of a multiphysics startup transient in HP-MR unit cell for use as tutorial in Sockeye. The full-core GC-MR model was used for assessment of SAM for balance of plant modeling and for demonstrating the SWIFT code capability for hydrogen redistribution modeling in multiphysics transient analyses. In this process, several bugs/issues were identified and reported to developers. Finally, the assembly GC-MR model developed in FY-2022 coupling Griffin, BISON and SAM through flow blockage and rod ejection transients was published to the National Reactor Innovation Center (NRIC) Virtual Test Bed (VTB). The Heat Pipe MicroReactor (HP-MR) concept high-fidelity multiphysics coupling of Griffin/BISON/Sockeye in load-following and heat pipe failure transients was also published on the VTB. Those submissions are enabling thorough review of these models as well as wide distribution to industry, regulator, and university users. In this analysis, several new research questions were uncovered, and follow-up analyses are recommended to further improve some models, consider additional transients, and continue development of VTB models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

ML-Based Power System Stability Assessment Considering Network Topology Changes: WECC 20,000+ Bus System Case Study

Modern power grids are fast-changing and thus require real-time monitoring and online stability assessment. With the rapid development of machine learning (ML) techniques, using data-driven models to provide fast and accurate estimations of power system stability marginal information, such as frequency nadir for frequency stability and critical clearing time (CCT) for transient stability, have become possible. However, despite the numerous research on ML-based methods for frequency nadir and CCT prediction, there is limited work on the impact of different network topology changes. Furthermore, most previous studies only focused on small or synthetic systems, and there is a lack of research on actual large power system models. In this paper, the above issues are addressed by studying the actual U.S. Western Electricity Coordinating Council (WECC) system model with more than 20,000 buses. Massive simulations are conducted in PowerWorld Simulator to study the impact of various topology change scenarios on both frequency stability and transient stability. System operating information is extracted from the success dispatch cases of various network topologies to generate a comprehensive dataset for ML-based models. Two ML methods, random forest (RF) and multilayer perceptron (MLP) neural network, are trained and tested for both frequency nadir prediction and CCT prediction. Test results have proven the models are capable of online stability assessment for large power networks such as the WECC system with sufficient accuracy.

critical clearing time↗

Oscillation Risks of Grid-Following and Grid-Forming Inverter-Based Resources in Series-Compensated Networks

Here, this paper investigates the dynamic behavior of a grid-connected inverter-based resource (IBR) when connected radially to a series compensated line. Potential interactions between the series compensation and the IBR have been identified for both types: grid-following (GFL) or grid-forming (GFM). The study begins with electromagnetic transient (EMT) simulations to demonstrate stability issues. Subsequently, nonlinear analytical models are formulated in the dq frame, validated against the EMT simulation, and leveraged to assess eigenvalues and participation factors. Influencing factors of the dominant oscillation modes have been identified. The analysis results show that series compensation may make a mode associated with the synchronization unit unstable. Furthermore, customized feedback systems are built for the synchronizing loop. Series compensation can increase the sensitivity of the voltage phase angle towards the synchronizing angle, and introduce phase lag in the real power response towards the synchronizing angle. These factors may cause interactions with the phase-locked loop in GFL-IBR systems and with power-based synchronization in GFM-IBR systems, potentially leading to instability.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Analysis of Grid-Forming Inverter Controls for Grid-Connected and Islanded Microgrid Integration

Autonomous grid-forming (GFM) inverter testbeds with scalable platforms have attracted interest recently. In this study, a self-synchronized universal droop controller (SUDC) was adopted, tested, and scaled in a small network and a test feeder using a real-time simulation tool to operate microgrids without synchronous generators. We presented a novel GFM inverter control adoption to better understand the dynamic behavior of the inverters and their scalability, which can impact the distribution system (DS). This paper provides a steady-state and transient analysis of the GFM power inverter controller via simulation to better understand voltage and frequency stabilization and ensure that the critical electric loads are not affected during a prolonged power outage. The controllers of the GFM inverter are simulated in HYPERSIM to examine voltage and frequency fluctuations. This analysis includes assessing the black start capability for photovoltaic microgrids, both grid-connected and islanded, during transient fault conditions. The high photovoltaic PV penetration levels open exciting opportunities and challenges for the DS. The GFM inverter control demonstrated appropriate response times for synchronization, connection, and disconnection to the grid. The DS has become more resilient and independent of fossil fuels by increasing the penetration of inverter-based distributed energy resources (DERs).

Ward, Laura (ORCID:0000000345501581)↗

Vibration attenuation using magnetic suspension isolation

The paper derives the vibration transfer function for a magnetic suspension in generic form, with specific reference to the Annular Suspension and Pointing System for Shuttle payloads. It is shown that a magnetically suspended pointer can exhibit considerably improved dynamic stability over conventional means, particularly at high frequencies and for transient disturbances. The amount of improvement to be anticipated is quantified to allow first-cut assessments for specific applications.

Hamilton, B. J.↗

Compressor Stall Recovery Through Tip Injection Assessed

Aerodynamic stability is a fundamental limit in the compressor design process. The development of robust techniques for increasing stability has several benefits: enabling higher loading and fewer blades, increasing safety throughout a mission, increasing tolerance to stage mismatch during part-speed operation and speed transients, and providing an opportunity to match stages at the compressor maximum efficiency point, thus reducing fuel burn. Mass injection upstream of the tip of a high-speed axial compressor rotor is a stability enhancement approach known to be effective in suppressing stall in tip-critical rotors if the injection is activated before stall occurs. This approach to stall suppression requires that a reliable stall warning system be available. Tests have recently been performed to assess whether steady injection can also be used to recover from fully developed stall. If mass injection is effective in recovering from stall quickly enough to avoid structural damage or loss of engine power, then a stall warning system may not be required. The stall recovery tests were performed on a transonic compressor rotor at its design tip speed of 1475 ft/sec using four injectors evenly spaced around the compressor case upstream of the rotor. The injectors were connected to an external air source. In an actual engine application, the injected air would be supplied with compressor bleed air. The injectors were isolated from the air source by a fast-acting butterfly valve. With the injectors turned off, the compressor was throttled into stall. Air injection was then activated with no change in throttle setting by opening the butterfly valve. The compressor recovered from stall at a fixed throttle setting with the aid of tip injection. The unsteady operating characteristic of the rotor was measured during these tests using high-response pressure sensors located upstream and downstream of the rotor. The figure shows the results, where the unsteady pressure and mass flow are superimposed on the steady operating characteristic. The total injected mass flow was equal to 1.3 percent of the compressor flow. The solid line with no solid squares on it denotes the operating point during the beginning of throttle closure and the initial drop into stall. The gray traces denote the operating point during an additional throttle closure that occurred over the next 1200 rotor revolutions (4 sec). The dashed line denotes the recovery from stall that occurred during 90 rotor revolutions (0.3 sec) after the injectors were activated with no change in throttle setting. Tip injection not only recovers the compressor from stall, but also restores the compressor to its pre-stall level of pressure rise. In contrast, standard stall recovery schemes such as compressor bleed, stator vane actuation, or engine throttle modulation result in a loss of pressure rise across the compressor, which results in a loss of engine power.

Suder, Ken L.↗

Transonic stability and control of aircraft using CFD methods

Implementation of a capability to calculate longitudinal short-period response in the CAP-TSD (Computational Aeroelasticity Program - Transonic Small Disturbance) finite-difference code is described. The code, developed recently at the NASA Langley Research Center, is capable of solving steady and unsteady flows about complete aircraft configurations and is used primarily for aeroelastic calculations in the critical transonic speed range. The longitudinal short-period equations of motion in state-space form have been coupled to the time-accurate lift and moment calculated by the program. Transient responses to an elevator pulse for free-flying aircraft demonstrate the new capability. A trim routine is also added to the code to obtain trim automatically during steady-state flow field convergence. Stability and control derivatives are estimated from the calculated transient response by a maximum likelihood estimation program. Results for a fighter configuration and a general aviation configuration are presented to assess the capability.

Vinh, Lam-Son↗

The initial-value problem for viscous channel flows

Plane viscous channel flows are perturbed and the ensuing initial-value problems are investigated in detail. Unlike traditional methods where traveling wave normal modes are assumed for solution, this works offers a means whereby completely arbitrary initial input can be specified without having to resort to eigenfunction expansions. The full temporal behavior, including both early time transients and the long time asymptotics, can be determined for any initial disturbance. Effects of three-dimensionality can be assessed. The bases for the analysis are: (a) linearization of the governing equations; (b) Fourier decomposition in the spanwise and streamwise directions of the flow; and (c) direct numerical integration of the resulting partial differential equations. All of the stability data that are known for such flows can be reproduced. Also, the optimal initial condition can be determined in a straight forward manner and such optimal conditions clearly reflect transient growth data that is easily determined by a rational choice of a basis for the initial conditions. Although there can be significant transient growth for subcritical values of the Reynolds number using this approach it does not appear possible that arbitrary initial conditions will lead to the exceptionally large transient amplitudes that have been determined by optimization of normal modes. The approach is general and can be applied to other classes of problems where only a finite discrete spectrum exists, such as the boundary layer for example.

Criminale, W. O.↗

Modeling Framework for Data Center

This chapter highlights the critical need for advanced modeling of data centers due to their rapidly increasing energy consumption and impact on grid reliability. Driven by the demand for AI applications, data centers are projected to consume a significant portion of US energy by 2028, putting stress on an already challenged power grid. The chapter emphasizes the importance of "fast" time-scale models to understand the dynamic interactions between data centers and the grid, especially given the rapid power fluctuations of AI workloads. It outlines a modeling framework that includes both offline and real-time EMT domain simulations, detailing the necessary representations for various components like utility interfaces, transformers, IT loads, UPS, cooling loads, Battery Energy Storage Systems (BESS), generators, protection systems, and higher-level control systems. While standard simulation tools like PSCAD offer basic models, custom development is often required to accurately capture the unique and fast-changing behaviors of modern data centers. The chapter also discusses key metrics and test cases for validating these models, focusing on transient load responses, protection relay coordination, and demand flexibility. Finally, it addresses the challenges of modeling large-scale data centers, such as computational complexity and the trade-off between model fidelity and practicality, suggesting hybrid modeling approaches as a solution. The overarching goal is to create a robust framework that helps assess data center impacts on grid stability, identify vulnerabilities, and inform the development of standards for reliable integration of these large loads into the bulk power system.

25 ENERGY STORAGE↗

Design practices for controlling spacecraft charging interactions

A design guidelines handbook prepared to provide criteria for assessing and minimizing spacecraft charging interactions is described. An evaluation philosophy of analyzing specific satellite designs in a substorm environment specification with NASCAP is proposed. Criteria for possible discharges are given and a technique for computing the discharge transients is outlined. The charging of a three axis stabilized satellite is examined to illustrate the philosophy. Possible discharge locations are found and transients computed. The effect of changing selected surface coatings is evaluated and found to substantially reduce charging levels.

Stevens, N. J.↗

Smooth Pursuit of Flicker-Defined Motion

We examined the pursuit response to stimuli defined by space-variant flicker of a dense random dot carrier pattern. On each frame, every element of the pattern could change polarity, with a probability given by a two-dimensional Gaussian distribution. A normal distribution produces a circular region of twinkle, while inverting the distribution results in a spot of static texture in a twinkling surround. In this latter case, the carrier texture could be stationary, or could move with the twinkle modulator,thereby producing first-order motion in the region of the spot. While the twinkle-defined spot produces a strong sensation of motion, the complementary stimulus defined by the absence of twinkle does not;when viewed peripherally, it appears to move in steps even when the generating distribution moves smoothly. We examined pursuit responses to these stimuli using two techniques: 1) the eye movement correlogram, obtained by cross-correlating eye velocity with the velocity of a randomly-moving stimulus; and 2) delayed visual feedback, where transient stabilization of a target can produce spontaneous oscillations of the eye, with a period empirically observed to vary linearly with the applied delay. Both techniques provide an estimate of the internal processing time, which can be as short as100 milliseconds for a first-order target. Assessed by the correlogram method, the response to flicker defined motion is delayed by more than 100 milliseconds, and significantly weaker (especially in the vertical dimension). When initially presented in the delayed feedback condition, purely saccadic oscillation is observed. One subject eventually developed smooth oscillations (albeit with significant saccadic intrusions), showing a period-versus-delay slope similar to that observed for first-order targets. This result is somewhat surprising, given that we interpret the slope of the period-versus-delayfunction as reflecting the balance between position- and velocity-sensitive inputs to pursuit.

smooth pursuit↗

Effects of buoyancy on gas jet diffusion flames - Experiment and theory

Theoretical and experimental research on the effects of buoyancy on gas-jet diffusion flames is described. Part of this research involves an assessment of existing data obtained under reduced-gravity conditions. The results show that uncertainties in the current understanding of flame structure exist and further research is required before reliable predictions of ignition, stabilization, and propagation of flames under microgravity conditions can be made. Steady-state and transient theories have been developed and used in the analysis of existing drop-tower data and new data obtained from a stationary experiment involving inverted flames. The result of this research has led to the definition of a microgravity experiment to be performed in space.

Edelman, R. B.↗