Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stability analysis”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Small-Signal Stability Analysis of Low-Inertia Power Grids with Inverter-Based Resources and Synchronous Condensers: Preprint

With the potential environmental impacts of conventional fossil fuels and the technology advances of grid interactive power electronics, inverter-based resources (IBRs) are playing a crucial role in modern power grids. However, due to the lack of physical inertia, IBRs are vulnerable to system disturbances. Further, conventional grid-following (GFL) inverters cannot actively establish the frequency and voltage, which could not be a appropriate solution for 100% inverter-based power grids. Grid-forming (GFM) control, as an emerging inverter control function, has drawn great attention in recent years, while it is still under development and pre-mature for large-scale field deployment. A potential solution for gradually transitioning into a 100% inverter-based power grid is to combine both IBRs and synchronous machines into a mixed generation portfolio. More importantly, synchronous machines can be operated as synchronous condensers (SCs) until they are retired. It is also worth mentioning that the stability analysis is of great importance to ensure the smooth transition from existing power grids to 100% inverter-based system, which is the main focus of this work with detailed modeling and modal analysis.

100% renewable↗
DOE OSTI2024

Holistic Small-Signal Stability Analysis for Large-Scale Inverter-Intensive Power Systems with Coupled and Full-Order Dynamics from Control Systems and Power Networks

The increasing penetration of inverter-based resources (IBRs) into the existing power systems introduces tremendous benefits for enhanced sustainability but also poses inevitable challenges in terms of insufficient inertia, potential instability, and complex network dynamics, among others. However, the additional coupling introduced by the interactions among gridfollowing (GFL) and grid-forming (GFM) IBRs and the other components (i.e., synchronous generators [SGs], loads, and network, etc.) has not been clearly explored. A holistic, scalable, and quantitative stability analysis framework with the control systems and power networks is still missing. Here, in this paper, to fill in the technical gaps, a holistic small-signal model of the entire system with both rotating generation units and IBRs is established. An extended power flow model with operation dynamics from both generator control schemes and power networks is proposed to provide the varying steady-state operating points for small-signal modeling. The proposed method is compared with MATLAB solvers, and the results show that the proposed approach has a minimum calculation time, which can be less than 12 seconds for a large-scale power system with up to 2,000 buses. Furthermore, a quantitative method is developed to identify the impacts of IBRs on system performance with emphases on the potential stability issues with GFL IBRs, additional benefits of employing GFM IBRs, the feasibility of replacing SGs with GFM IBRs, and the impact of penetration level of different kinds of generation units. Finally, a field island power system is used to verify the proposed approach, and hardware-in-the-loop (HIL) tests are provided to further demonstrate the effectiveness of the proposed analysis.

14 SOLAR ENERGY↗
DOE OSTI2022

The role of plasma instabilities in relativistic radiation-mediated shocks: stability analysis and particle-in-cell simulations

ABSTRACT Relativistic radiation-mediated shocks are likely formed in prodigious cosmic explosions. The structure and emission of such shocks are regulated by copious production of electron–positron pairs inside the shock-transition layer. It has been pointed out recently that substantial abundance of positrons inside the shock leads to a velocity separation of the different plasma constituents, which is expected to induce a rapid growth of plasma instabilities. In this paper, we study the hierarchy of plasma microinstabilities growing in an electron-ion plasma loaded with pairs and subject to a radiation force. Linear stability analysis indicates that such a system is unstable to the growth of various plasma modes which ultimately become dominated by a current filamentation instability driven by the relative drift between the ions and the pairs. These results are validated by particle-in-cell simulations that further probe the non-linear regime of the instabilities, and the pair-ion coupling in the microturbulent electromagnetic field. Based on this analysis, we derive a reduced-transport equation for the particles via pitch-angle scattering in the microturbulence and demonstrate that it can couple the different species and lead to non-adiabatic compression via a Joule-like heating. The heating of the pairs and, conceivably, the formation of non-thermal distributions, arising from the microturbulence, can affect the observed shock-breakout signal in ways unaccounted for by current single-fluid models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

FLAW STABILITY ANALYSIS OF SURFACE CRACKS IN DOE STANDARD CANISTERS UNDER OPERATION LOADS AND WELDING RESIDUAL STRESSES

There are over 3000 commercial spent nuclear fuel (SNF) storage canisters that are made of stainless steel and planned for multi-purpose functions of SNF handling: interim storage, transportation, and ultimate disposal at a future SNF disposition location. To date, many multi-purpose canisters (MPC) are located in coastal regions for long-term storage. The canisters are fabricated by welding, and post-welding heat treatment is not required for relieving welding residual stresses (WRS). As a result, these canisters may be susceptible to chloride induced stress corrosion cracking (CISCC) due to the chloride-bearing marine salts. The flaw stability analysis of MPC canisters was performed at SRNL (i.e., PVP2016-4935 and PVP2016-63887) for through-wall flaws and for surface flaws respectively, where WRS was considered, and the failure assessment diagram (FAD)-based fracture mechanics method codified by API 579-1/ASME FFS-1-2007 Edition was adopted.

DOE standard canister↗
DOE OSTI2021

A hierarchical gray-box dynamic modeling methodology for direct-expansion cooling systems to support control stability analysis; Méthodologie de modélisation dynamique hiérarchique de type boîte grise pour les systèmes de refroidissement à détente directe afin de soutenir l'analyse de stabilité de la commande

In this paper, a gray-box dynamic modeling approach for direct-expansion cooling systems is presented. The overall approach incorporates a multi-stage training procedure that consists of 1) identification of component sub-models from quasi-steady-state performance data, 2) system model integration with estimation of refrigerant charge and 3) fine tuning of thermal capacitances of the evaporator and condenser to capture the system dynamic responses. Compared to traditional physics-based models, the proposed modeling approach has advantages including reduced engineering efforts in the model development phase, improved computational efficiency and enhanced prediction accuracy. The modeling method was validated using a 3-ton variable-speed heat pump and proved to be capable of accurately predicting the system transient behaviors over a wide range of operating conditions. The established dynamic model was then applied for control stability analysis, with a specific goal of determining a proper control execution time step. The case study results showed that the stable control execution time step could change significantly, from 3 sec to 19 sec, as the operating conditions and control settings vary, and a proper selection of the execution time step is critical to ensure stable and reliable operations.

42 ENGINEERING↗
DOE OSTI2025

Stability analysis of Alfvén eigenmodes excited by ion cyclotron resonance heating on EAST

Alfvén wave instabilities driven by energetic particles are common in fusion devices. Understanding their behavior is essential for the good confinement of fusion plasma. For this purpose, a series of experiments are performed on EAST to investigate the excitation of toroidal Alfvén eigenmodes (TAEs). Experimentally, it was found that AEs with frequencies around 80, 134, 157 kHz are excited by ion cyclotron resonance heating (ICRH) hydrogen minority heating scenario. Moreover, the excitation of AEs is independent of the wall-coating materials, but determined with the generation and confinement of the fast ions. Statistical analysis suggests that the TAEs can only be excited when the ICRH power is larger than 2 MW and the H 98 factor is larger than 1.15. In line with the experiments, a set of simulations using TORIC, ASCOT, and NOVA-C is performed. The simulation results show that with the experimental density profile and safety factor, the measured TAE with frequency of 134 kHz is well reproduced by taking into account the plasma toroidal rotation frequency. However, the appearance of the 80 kHz mode cannot be captured in the simulations. The ICRH-generated fast hydrogen ions mainly have banana orbits and that their perpendicular velocity is much larger than that of the parallel ions. The stability analysis with NOVA-C suggests that the radial gradient of the fast H ion distribution is the only driving source of TAEs. For TAEs locate in the center of Alfvén continuum, the D ion Landau damping is the dominant damping term; For most other TAEs which have intersections with the continuum boundary, the continuum damping becomes the dominant damping mechanism.

Alfvén eigenmodes↗
DOE OSTI2020

Stability Analysis of Interface Conditions for Ocean-Atmosphere Coupling

In this paper, we analyze the stability of different coupling strategies for multidomain PDEs that arise in general circulation models used in climate simulations. We focus on fully coupled ocean–atmosphere models that are needed to represent and understand the complicated interactions of these two systems, becoming increasingly important in climate change assessment in recent years. Numerical stability issues typically arise because of different time-stepping strategies applied to the coupled PDE system. In particular, the contributing factors include using large time steps, lack of accurate interface flux, and single-iteration coupling. We investigate the stability of the coupled ocean–atmosphere models for various interface conditions such as the Dirichlet–Neumann condition and the bulk interface condition, which is unique to climate modeling. By analyzing a simplified model, we demonstrate here how the parameterization of the bulk condition and other numerical and physical parameters affect the coupling stability and establish stability conditions for different coupling strategies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Nonideal stability analysis of differentially rotating plasmas with global curvature effects

The linear stability of global nonaxisymmetric modes in differentially rotating, magnetized, nonideal plasma is critical to classifying turbulence and transport phenomena. We investigate the competition between the local magneto-rotational instability (MRI) and the magneto-curvature instability (MCI)—a distinct nonaxisymmetric low-frequency curvature-driven global branch that appears alongside MRI. Here, to accomplish this, we developed a nonideal global spectral method, which is validated against NIMROD code simulations. This spectral approach allows for the direct derivation of an extended effective potential formalism and a resistive Alfvénic resonance condition, providing a framework for direct analysis of energy contributions and confinement mechanisms. Our study reveals that the global, low-frequency MCI persists at low magnetic Reynolds numbers (Rm), whereas the localized, high-frequency MRI is stabilized by diffusive broadening of its structure around its Alfvénic resonances. Consequently, we identify the global MCI branch as the primary onset mechanism for nonaxisymmetric magnetohydrodynamic instability in systems with finite curvature, e.g., astrophysical rotators. We establish distinct parameter regimes for mode dominance: MCI prevails in geometrically moderate-thickness disks with intermediate curvature and radial gaps, while MRI dominates in thin, low-curvature disks with large radial gaps. Mode competition is also highly sensitive to the flow profile, particularly vorticity and its gradient, with nonuniform shear profiles exhibiting more robust instability due to flow curvature (i.e., the second derivative of the flow profile) and shear contributions. A key outcome is the development of spectral diagrams derived from the global spectral method. These diagrams comprehensively map dominant instabilities and their characteristics, offering a predictive tool for critical onset parameters (i.e., flow curvature, magnetic field, and Rm) and facilitating the interpretation of experimental and simulation results. Notably, these diagrams demonstrate that the global MCI is generally the sole unstable mode at the initial onset of nonaxisymmetric instability.

Haywood, Alexander [Princeton Univ., NJ (United St↗
DOE OSTI2023

Hierarchical Control and Stability Analysis for a Nonisolated Grid-Tied DC Energy Router Integrating Energy Storage and Partial Distributed Generation

This article proposes a nonisolated dc converter-based energy router (dc-ER) and its operating strategy. Here, the intent is to integrate energy storage (ES), distributed generation (DG), local load, and dc power grid in an autonomous and more efficient way. The ES/DG power ports are coupled with each other in partial-type connections, which obtains higher voltage supply gain and direct input–output power transmission. High-voltage supply gain allows a wide solar power tracking range and the possibility of optimal battery charging/discharging. Direct input–output power transmission improves the energy conversion efficiency. In this article, the operating modes of dc-ER are first analyzed, followed by the optimal hardware design counting the DG current ripple minimization and the limitation caused by the power flow direction. Second, the mathematical model of dc-ER is deduced. The hierarchical control is then proposed to manage port energy in a flexible manner. The stability is analyzed by using impedance modeling. Finally, experimental and simulation results verify the superiorities of the proposed topology in terms of flexible operating mode transition and high-power conversion efficiency.

25 ENERGY STORAGE↗
DOE OSTI2025

Stability analysis of the Eulerian–Lagrangian finite volume methods for nonlinear hyperbolic equations in one space dimension

In this paper, we construct a novel Eulerian–Lagrangian finite volume (ELFV) method for nonlinear scalar hyperbolic equations in one space dimension. It is well known that the exact solutions to such problems may contain shocks though the initial conditions are smooth, and direct numerical methods may suffer from restricted time step sizes. To relieve the restriction, we propose an ELFV method, where the space-time domain was separated by the partition lines originated from the cell interfaces whose slopes are obtained following the Rakine–Hugoniot junmp condition. Unfortunately, to avoid the intersection of the partition lines, the time step sizes are still limited. To fix this gap, we detect effective troubled cells (ETCs) and carefully design the influence region of each ETC, within which the partitioned space-time regions are merged together to form a new one. Then with the new partition of the space-time domain, we theoretically prove that the proposed first-order scheme with Euler forward time discretization is total-variation-diminishing and maximum-principle-preserving with at least twice larger time step constraints than the classical first order Eulerian method for Burgers’ equation. Numerical experiments verify the optimality of the designed time step sizes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Comparative Small-Signal Stability Analysis of Grid-Forming and Grid-Following Inverters in Low-Inertia Power Systems: Preprint

Inverter-based resources (IBRs) are gradually replacing conventional fossil energy in modern power systems, which could make the system susceptible due to insufficient inertia, especially when moving towards a 100% inverter-based power grid. Grid-forming (GFM) and grid-following (GFL) controls are two common inverter control algorithms. However, the corresponding stability issues along with a comparative study on low-inertia power systems should still be investigated. This paper focuses on small-signal stability and conducts a comprehensive analysis of GFM and GFL controls, with emphases on the interactions between synchronous generators and IBRs. The detailed modeling of GFM and GFL inverters, and synchronous generators (SGs) are established. The coupling mechanism of IBRs and SGs is then conducted through participation factor analysis. The impacts of critical parameters on system stability are finally identified. The theoretical analyses are validated through time-domain simulation in MATLAB/Simulink.

grid-following↗
DOE OSTI2022

Explainable Bayesian Neural Network for Probabilistic Transient Stability Analysis Considering Wind Energy

While several data-driven models have been developed for transient stability assessment, how to consider the uncertainties from load and renewable generations and provide interpretation of data-driven assessment results are still open. This paper proposes an explainable Bayesian Neural Network (BNN) for probabilistic transient stability assessment (TSA). By extracting the uncertainties from loads and wind farms, the BNN model can make a reliable prediction and quantify the prediction uncertainties. We also develop the Gradient Shap algorithm to make the global and local explanations for the probabilistic TSA model, a significant advantage over existing black-box data-driven methods. Numerical results on the modified IEEE 39-bus system show that the proposed method outperforms the existing methods in terms of prediction accuracy and uncertainty quantification capabilities. The explainability of the proposed method allows system operators to design preventive controls for enhancing system stability.

Bayesian Neural Network↗
DOE OSTI2024

Partitioning of Large-Scale Power Electronics-Based Power Systems for Small-Signal Stability Analysis

The nodal admittance matrix (NAM)-based approach is suitable for analyzing the small-signal stability of large-scale power electronics-based power systems (PEPSs) as it preserves the system structure by utilizing the admittance matrix. Previously, NAM-based area partition has been proposed, which divides the system into various subareas and interconnections for easier analysis of the low-dimension matrix compared to the entire system-based high-dimension matrix. However, no partition algorithm has been presented for the NAM-based area partition method. This paper focuses on implementing the spectral partitioning algorithm for partitioning large-scale PEPSs into a low-dimension matrix to reduce the computation complexity of the analysis. These spectral components facilitate data transformation into a new space, enabling the application of traditional clustering methods like k-means. To evaluate the performance of the partitioning method, the subareas and interconnections obtained from the spectral clustering algorithm are incorporated into the NAM-based area partition method for a large system with 140 buses. The computational times of the original method, where the NAM-based criterion is directly applied to the entire system, are compared with those of the NAM-based partition method in MATLAB. PSCAD simulations of the whole system and the obtained subareas are conducted to validate the effectiveness of the proposed algorithm.

Nupur, Nupur↗
DOE OSTI2021

Linear stability analysis of two fluid columns of different densities and viscosities in a gravity field

The linear stability of a vertical interface separating two miscible fluid columns of different densities and viscosities under the influence of gravity is investigated. This flow possesses a time-dependent reference state (each column accelerates at different rates owing to their different densities) and the interface thickness grows as the square root of time (by diffusion). Numerical integration of the linear initial-value problem is carried out and discussed in detail as a function of vertical and spanwise wavenumbers and the flow parameters. Adjoint-based optimization is performed in order to determine initial conditions that lead to maximum growth of disturbances in finite time. Results indicate that the rate of growth of the perturbation energy at small wavenumbers (less affected by viscosity initially) is dominated by two-dimensional modes (no spanwise variation). Substantial transient growth is observed at higher wave modes initially, followed by asymptotic decay of the perturbations at large time. Sensitivity of perturbation growth with respect to initial time, density and viscosity ratios is investigated. This work is complementary to previous inviscid analysis of this configuration, which showed that the interface was unconditionally unstable at all wave modes, even in the presence of surface tension, and that instability grew as the exponential of time squared.

Mechanics↗
DOE OSTI2021

Flame stability analysis of flame spray pyrolysis by artificial intelligence

Flame spray pyrolysis (FSP) is a process used to synthesize nanoparticles through the combustion of an atomized precursor solution; this process has applications in catalysts, battery materials, and pigments. Current limitations revolve around understanding how to consistently achieve a stable flame and the reliable production of nanoparticles. Machine learning and artificial intelligence algorithms that detect unstable flame conditions in real time may be a means of streamlining the synthesis process and improving FSP efficiency. In this study, the FSP flame stability is first quantified by analyzing the brightness of the flame's anchor point. This analysis is then used to label data for both unsupervised and supervised machine learning approaches. The unsupervised learning approach allows for autonomous labeling and classification of new data by representing data in a reduced dimensional space and identifying combinations of features that most effectively cluster it. The supervised learning approach, on the other hand, requires human labeling of training and test data but is able to classify multiple objects of interest (such as the burner and pilot flames) within the video feed. The accuracy of each of these techniques is compared against the evaluations of human experts. Both the unsupervised and supervised approaches can track and classify FSP flame conditions in real time to alert users of unstable flame conditions. This research has the potential to autonomously track and manage flame spray pyrolysis as well as other flame technologies by monitoring and classifying the flame stability.

42 ENGINEERING↗
DOE OSTI2021

Impedance Methods for Analyzing Stability Impacts of Inverter-Based Resources: Stability Analysis Tools for Modern Power Systems

Power systems around the world are undergoing a major transformation because of the increasing shares of renewable energy, growing deployment of energy storage systems, proliferation of distributed energy resources, electrification of other sectors, and so on. In the United States, wind and solar provided almost 10% of electricity in 2019. The U.S. Energy Information Administration, in its 2020 Annual Energy Outlook, forecasted that the share of electricity from renewables will reach 38% by 2050, of which more than 80% will come from wind and solar. Wind and solar, along with battery energy storage systems, interface with the grid using power electronic inverters; hence, they are collectively referred to as inverter-based resources (IBRs). The increasing annual share of electricity from IBRs in a power system means that during more times of the year, the system will operate at a much higher concentration of IBRs. Figure 1 presents the hourly share of wind and solar generation in the Electric Reliability Council of Texas (ERCOT) system in Texas in 2019. While the annual wind share was at 20%, the instantaneous percentage share was much higher. Moments of high shares of IBRs (>50%) will continue to grow as more IBRs will be deployed in a power system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Quantitative Small-Signal Stability Analysis of Grid-Forming Inverters in Large-Scale Power Systems

In this paper, the challenges induced by the integration of inverter-based resources (IBRs) into existing power systems are discussed. Although the grid-forming (GFM) control has been used to enhance system performance in multiple sections of electric power systems (e.g., transmission and distribution systems, microgrids, among others), the coupling mechanism between GFM inverter and various types of other components in modern power systems should be clearly explored. Existing analyses mainly focus on qualitative investigations that verify the effectiveness of deploying GFM inverters for stability enhancement. A holistic system-level quantitative analysis is still missing to verify the detailed benefits of GFM inverters and their interactions with other generation units (e.g., synchronous generators [SGs], conventional grid-following [GFL] inverters, synchronous condensers, among others). To fill in the technical gaps, a holistic small-signal model of the entire system with both rotating generation units and IBRs is established along with the customized power flow analysis for varying steady-state operating points. Further, a quantitative method is proposed to identify the impacts of GFM inverters on system performance with an emphasis on a quantifiable stability margin and a complete modal analysis at the system level. Finally, a realistic test system is used to verify the proposed approach.

grid-following↗