Search NASA⌕ Search

SEARCH · Search NASA

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 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Granger Causality for prediction in Dynamic Mode Decomposition: Application to power systems

Here, the dynamic mode decomposition (DMD) technique extracts the dominant modes characterizing the innate dynamical behavior of the system within the measurement data. For appropriate identification of dominant modes from the measurement data, the DMD algorithm necessitates ensuring the quality of the input measurement data sequences. On that account, for validating the usability of the dataset for the DMD algorithm, the paper proposed two conditions: Persistence of excitation (PE) and the Granger Causality Test (GCT). The virtual data sequences are designed with the hankel matrix representation such that the dimensions of the subspace spanning the essential system modes are increased with the addition of new state variables. The PE condition provides the lower bound for the trajectory length, and the GCT provides the order of the model. Satisfying the PE condition enables estimating an approximate linear model, but the predictability with the identified model is only assured with the temporal causation among data searched with GCT. The proposed methodology is validated with the application for coherency identification (CI) in a multi-machine power system (MMPS), an essential phenomenon in transient stability analysis. The significance of PE condition and GCT is demonstrated through various case studies implemented on 22 bus six generator system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2018

Fault Detection for Dynamical Systems using Differential Geometric and Concurrent Learning Approach

This paper presents a fault identification and estimation approach that brings together the differential geometric concept of observability codistribution with data driven concurrent learning. In order to identify faults in presence of unknown disturbances, we utilize the differential geometric approach to design a coordinate transformation, to find a subspace in which the effect of disturbances and system faults can be segregated. We then use concurrent learning to estimate magnitude of the constant fault. We illustrate the approach to fault isolation for a linear helicopter dynamics and a spherical pendulum dynamics. We use Lyapunov stability analysis to show that the fault estimate by concurrent learning converges to the actual fault value, and then illustrate the design of a recovery controller.

Fault detection and diagnosis↗
DOE OSTI2014

The iterative thermal emission method: A more implicit modification of IMC

For over 40 years, the Implicit Monte Carlo (IMC) method has been used to solve challenging problems in thermal radiative transfer. These problems typically contain regions that are optically thick and diffusive, as a consequence of the high degree of “pseudo-scattering” introduced to model the absorption and reemission of photons from a tightly-coupled, radiating material. IMC has several well-known features that could be improved: a) it can be prohibitively computationally expensive, b) it introduces statistical noise into the material and radiation temperatures, which may be problematic in multiphysics simulations, and c) under certain conditions, solutions can be nonphysical, in that they violate a maximum principle, where IMC-calculated temperatures can be greater than the maximum temperature used to drive the problem. We have developed a variant of IMC called iterative thermal emission IMC, which is designed to have a reduced parameter space in which the maximum principle is violated. ITE IMC is a more implicit version of IMC in that it uses the information obtained from a series of IMC photon histories to improve the estimate for the end of time step material temperature during a time step. A better estimate of the end of time step material temperature allows for a more implicit estimate of other temperature-dependent quantities: opacity, heat capacity, Fleck factor (probability that a photon absorbed during a time step is not reemitted) and the Planckian emission source. We have verified the ITE IMC method against 0-D and 1-D analytic solutions and problems from the literature. These results are compared with traditional IMC. We perform an infinite medium stability analysis of ITE IMC and show that it is slightly more numerically stable than traditional IMC. We find that significantly larger time steps can be used with ITE IMC without violating the maximum principle, especially in problems with non-linear material properties. The ITE IMC method does however yield solutions with larger variance because each sub-step uses a different Fleck factor (even at equilibrium)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Effect of inlet throttling on thermohydraulic instability in a large scale water-based RCCS: An experimental study

The objective of the present experimental study is to investigate the effect of inlet throttling on the thermohydraulic stability of a large scale water-based Reactor Cavity Cooling System (RCCS). The test was performed using the water-based Natural convection Shutdown heat removal Test Facility (NSTF) at Argonne, which represented a ½ axial scale and 12.5° sector slice of the full scale Framatome 625 MW t SC-HTGR RCCS concept. A two-phase steady state was first established through direct condensate refill, followed by increased inlet throttling over 10 stages, corresponding to a loss coefficient K over the range of 0.05–653. With the inlet throttling gradually increased, the system experienced a unique transition process between stabilization and destabilization. Through a stability analysis, three instability mechanisms were identified in the present test, including a compound mechanism due to both natural circulation oscillations (NCOs) and density wave oscillations (DWOs), Type-II DWOs, and geysering.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Broadband impedance modeling and single bunch instabilities estimations of the advanced light source upgrade project

We report the Advanced Light Source Upgrade (ALS-U) is a 4th generation diffraction-limited soft x-ray radiation source, consisting of a new accumulator ring (AR) and a new storage ring (SR). In both rings coupling-impedance driven instabilities need careful evaluation to ensure meeting the machine’s high-performance goals. This paper presents the workflow followed in building the impedance models and the beam-stability analysis based on those models. We follow the commonly accepted approach of separating the resistive-wall and the geometric parts of the impedance; the former is obtained by analytical formulas, the latter by numerical electro-magnetic codes (primarily CST Studio software) with perfectly-conducting boundary conditions. Impedance budgets are established and pseudo-Green functions calculated to be used in beam dynamics studies. We also present various ways to cross-check simulation results for reliable impedance modeling. Finally, the crucial single-bunch instability current thresholds for various operation modes are determined and discussed.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2021

Migration of inclusions in a matrix due to a spatially varying interface energy

An interfacial energy can be a function of a bulk field such as temperature or electric field. Here we find that in a system with a gradient in temperature or electric potential, the resulting variation in interfacial energy can induce a particle to migrate by either surface or bulk diffusion. For a circular particle under a constant unidirectional gradient in the bulk field, the field dependence of the interfacial energy induces a circular particle to move as a circle at a constant velocity along the direction from higher to lower interfacial energy. A linear stability analysis of this steady state migration suggests that perturbations will damp out as time evolves, and thus under these conditions a migrating circular particle is morphologically stable. Other spatial distributions of interface energy can lead to the distortion of an initially circular shaped particle during migration. A phase field model is developed that captures these distortions and verifies the theoretical results mentioned above.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Simulating ice-shelf extent using damage mechanics

Abstract Inaccurate representations of iceberg calving from ice shelves are a large source of uncertainty in mass-loss projections from the Antarctic ice sheet. Here, we address this limitation by implementing and testing a continuum damage-mechanics model in a continental scale ice-sheet model. The damage-mechanics formulation, based on a linear stability analysis and subsequent long-wavelength approximation of crevasses that evolve in a viscous medium, links damage evolution to climate forcing and the large-scale stresses within an ice shelf. We incorporate this model into the BISICLES ice-sheet model and test it by applying it to idealized (1) ice tongues, for which we present analytical solutions and (2) buttressed ice-shelf geometries. Our simulations show that the model reproduces the large disparity in lengths of ice shelves with geometries and melt rates broadly similar to those of four Antarctic ice shelves: Erebus Glacier Tongue (length ~ 13 km), the unembayed portion of Drygalski Ice Tongue (~ 65 km), the Amery Ice Shelf (~ 350 km) and the Ross Ice Shelf (~ 500 km). These results demonstrate that our simple continuum model holds promise for constraining realistic ice-shelf extents in large-scale ice-sheet models in a computationally tractable manner.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

ARC physics basis–magnetohydrodynamics

ARC is designed to produce 400 ⁢MW of net electricity and prove the commercial feasibility of a fusion power plant. In order to achieve this goal ARC has to operate with optimal core performance in a stationary scenario that minimises wear on the first wall and divertor. This requires avoiding or mitigating magnetohydrodynamic (MHD) instabilities which have the potential to not only degrade the plasma core but also lead to deleterious transient heat loads on plasma facing components. Therefore, this work aims at characterising the MHD stability of the high performance ARC scenario and inform the design of error field correction coils. Firstly, simulations of vertical displacement events show that an in-vessel coil is not needed and instead the poloidal shaping coils can be used to control vertical stability. These simulations also inform the demands on the corresponding coil power supplies. Stability analysis of the ideal kink mode with or without a conducting wall and kinetic effects suggests that the ARC baseline scenario operates deeply in the stable region. Using RDCON, tearing modes at the 𝑚/𝑛 =2/1 and 3/2 surfaces (with poloidal mode number 𝑚, and toroidal mode number 𝑛) are shown to be linearly stable, and including thermal transport effects in the rational surfaces lead to further stabilisation. However, other transient plasma instabilities can seed neoclassical tearing modes (NTMs). The marginally stable width of NTMs in ARC strongly depends on the internal inductance and can fall below 0.1% of the normalised poloidal flux. Furthermore, an empirical cross-machine model of the 𝑛 =1 error field leading to a disruption predicts a critical error field larger than SPARC but smaller than ITER. Three-dimensional coils can be designed with the Generalised Purturbed Equilbium Code based on a simple model that calculates the maximum correctable error field that is limited by the neoclassical toroidal viscosity torque. Broad scans of different coil geometries identify a set of 2 rows of off-midplane coils to be a suitable solution. It is also determined that such a set of three-dimensional coils is capable of correcting 𝑛 =2 error fields to some degree and creating strong enough 𝑛 =2 or 𝑛 =3 edge resonant perturbation fields for the suppression of edge-localised modes at reasonable coil currents. The final design of the first ARC will be further informed by results from SPARC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Complex Orbitals, Multiple Local Minima, and Symmetry Breaking in Perdew–Zunger Self-Interaction Corrected Density Functional Theory Calculations

Implentation of seminumerical stability analysis for calculations using the Perdew-Zunger self-interaction correction is described. It is shown that real-valued solutions of the Perdew-Zunger equations for gas phase atoms are unstable with respect to imaginary orbital rotations, confirming that a proper implementation of the correction requires complex-valued orbitals. The orbital density dependence of the self-interaction corrected functional is found to lead to multiple local minima in the case of the acrylic acid, H 6 , and benzene molecules. Finally, in the case of benzene, symmetry breaking that results in incorrect ground state geometry is found to occur, erroneously leading to alternating bond lengths in the molecule.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Simultaneous Decomposition and Dewetting of Nanoscale Alloys: A Comparison of Experiment and Theory

In this work, we consider the coupled process of phase separation and dewetting of metal alloys of nanoscale thickness deposited on solid substrates. The experiments involve applying nanosecond laser pulses that melt the Ag40Ni60 alloy films in two setups: either on thin supporting membranes or on bulk substrates. These two setups allow for extracting both temporal and spatial scales on which the considered processes occur. The theoretical model involves a longwave version of the Cahn–Hilliard formulation used to describe spinodal decomposition, coupled with an asymptotically consistent longwave-based description of dewetting that occurs due to destabilizing interactions between the alloy and the substrate, modeled using the disjoining pressure approach. Careful modeling, combined with linear stability analysis and fully nonlinear simulations, leads to results consistent with the experiments. In particular, we find that the two instability mechanisms occur concurrently, with the phase separation occurring faster and on shorter temporal scales. The modeling results show a strong influence of the temperature dependence of relevant material properties, implying that such a dependence is crucial for the understanding of the experimental findings. The agreement between theory and experiment suggests the utility of the proposed theoretical approach in helping to develop further experiments directed toward formation of metallic alloy nanoparticles of desired properties.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Morphological instability of aqueous dissolution of silicate glasses and minerals [Nonlinear dynamics of aqueous dissolution of silicate glasses and minerals: Morphological instability]

Understanding of aqueous dissolution of silicate glasses and minerals is of great importance to both Earth science and materials science. Silicate dissolution exhibits complex temporal evolution and spatial pattern formations. Recently, we showed how observed complexity could emerge from a simple self-organizational mechanism: dissolution of the silica framework in a material could be catalyzed by the cations released from the reaction itself. This mechanism enables us to systematically predict many key features of a silicate dissolution process including the occurrence of a sharp corrosion front (vs. a leached surface layer), oscillatory dissolution and multiple stages of the alteration process (e.g., an alteration rate resumption at a late stage of glass dissolution). Here, through a linear stability analysis, we show that this same mechanism can also lead to morphological instability of an alteration front, which, in combination with oscillatory dissolution, can potentially lead to a whole suite of patterning phenomena, as observed on archaeological glass samples, including wavy dissolution fronts, growth rings, incoherent bandings of alteration products, and corrosion pitting. Here, the result thus further demonstrates the importance of the proposed self-accelerating mechanism in silicate material degradation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Substituted SrFeO 3 as robust oxygen sorbents for thermochemical air separation: correlating redox performance with compositional and structural properties

Thermochemical air separation via cyclic redox reactions of oxide-based oxygen sorbents has the potential to achieve high energy efficiency. Although a number of promising sorbents have been investigated, further improvements in sorbent performance through a fundamental understanding of the structure–performance relationships are highly desirable. In this study, we systematically investigated the effects of A and B site dopants on the oxygen uptake/release properties (i.e., vacancy formation energy, reduction enthalpy, oxygen release temperature, and oxygen capacity) of the SrFeO 3 family of perovskites as oxygen sorbents. A monotonic correlation between DFT calculated oxygen vacancy formation energy and oxygen release temperature demonstrates the effectiveness of DFT for guiding sorbent selection. Combining vacancy formation energy with stability analysis, dopants such as Ba and Mn have been identified for tuning the redox property of SrFeO 3 sorbents, and increasing the oxygen capacity for temperature and pressure swings when compared to undoped SrFeO 3 . The Mn doped sample proved to be highly stable, with less than a 3% decrease in capacity over 1000 cycles. Although the dynamic nature of the redox process makes it difficult to use a single vacancy formation energy as the descriptor, a systematic approach was developed to correlate the oxygen storage capacities with the sorbents’ compositional properties and vacancy formation energies. The combination of DFT calculations with experimental studies from this study provides a potentially effective strategy for developing improved sorbents for thermochemical air separation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Bandgap engineering of SrZrS 3 chalcogenide perovskite via substitutional doping for photovoltaic applications: a first-principles DFT study

Strontium zirconium sulfide (SrZrS 3 ) has garnered significant attention for photovoltaic (PV) applications due to its excellent optoelectronic properties, high chemical and moisture stability, and non-toxicity. However, the bandgaps of both the α- and β-phases lie outside the optimum ranges for both single-junction solar cells (SJSCs) and tandem solar cells (TSCs), thereby limiting their applications in PV technologies. In this study, we employed hybrid density functional theory to engineer the band gaps of α- and β-SrZrS 3 through substitutional doping. We considered three dopants (Hf, Sn, and Ti) at the Zr site of SrZrS 3 with varying doping concentrations and investigated their effects on the structural, electronic, and optical properties of the materials. We found that Sn and Ti doping effectively lowers the band gaps of both α- and β-SrZrS 3 , whereas Hf doping increases them. For x values up to 0.25, the band gaps of the α-SrZr 1−x Sn x S 3 and α-SrZr 1−x Ti x S 3 are within the optimum range for SJSCs, and those of β-SrZr 1−x Sn x S 3 and β-SrZr 1−x Ti x S 3 lie within the optimum range for Si/perovskite as well as perovskite/perovskite TSCs. The three dopants exhibited significant effects on the optical properties of both α- and β-SrZrS 3 , including the absorption coefficient, energy-loss functions, reflectivity, and refractivity spectra. Thermodynamic stability analysis revealed that for both phases, SrZr 1−x Hf x S 3 can be synthesized via exothermic processes, whereas the formation of SrZr 1−x Ti x S 3 and SrZr1−xSn x S 3 is endothermic and hence, not thermodynamically favorable. Further analysis showed that SrZr 1−x Ti x S 3 in both α- and β-phases are stable under thermodynamic equilibrium conditions, whereas SrZr 1−x Sn x S 3 is prone to dissociation into ternary phases (SrZrS 3 and SrSnS 3 ), especially at higher doping concentrations. These results show that Ti doping is effective in tuning the band gaps of α- and β-SrZrS 3 toward the optimal values for PV applications.

14 SOLAR ENERGY↗
DOE OSTI2025

Analytical Identification Method of Generalized Short‐Circuit Ratio Using Phasor Measurement Units

This paper introduces a novel analytical approach for the identification of the admittance matrix and the generalized short-circuit ratio (gSCR) in power systems integrated with renewable energy sources. The proposed method leverages voltage and current measurements from phasor measurement units (PMUs) to construct a least squares objective function, which is then solved using matrix calculus and partial derivatives. Unlike conventional optimization algorithms, this approach provides an analytical solution that substantially reduces data requirements, enabling the efficient and accurate identification of the gSCR with smaller datasets. Additionally, its fixed computational complexity allows for real-time updates as new data are collected, ensuring continuous refinement of the system of equations and enabling rapid, precise gSCR calculations. The method also exhibits strong robustness against measurement noise, making it well-suited for practical applications in dynamic power systems. The combination of reduced data requirements, real-time adaptability, noise robustness and fixed computational load establishes this method as a highly efficient and reliable tool for real-time power system stability analysis. Case studies on an EPRI 36-bus system demonstrate the method's effectiveness, highlighting its accuracy in closely matching true gSCR values, even under diverse disturbances and noisy conditions.

Han, Zelei [Hohai University, Nanjing (China)] (OR↗
DOE OSTI2023

Optimal Dispatch for the US Virgin Islands to Increase Renewable Rates in Saint Croix

The US Virgin Islands, like many other small island nations, face significant challenges related to its energy supply. Their dependency on imported fossil fuels not only drives up energy costs, but also contributes to different environmental problems. Nevertheless, their generally sunny and windy climate makes this archipelago an ideal location for renewable energy development. Increasing the use of sources such as wind, and solar can help reduce their energy dependency, lower their levelized cost of electricity, and improve the health and wellbeing of local communities. In this paper, an optimal dispatcher is employed to analyze the implications on energy security of employing different technology options according to the plans of the local system operator. In this context, short term stability is ensured by inserting operational constraints in the energy management system based on stability analysis developed in PSCAD. The results highlight the usefulness of integrating PV and battery system, while further analysis are required to justify integrating wind power.

energy islands↗
DOE OSTI2014

Overstability of acoustic waves in strongly magnetized anisotropic magnetohydrodynamic shear flows

We present a linear stability analysis of the perturbation modes in anisotropic magnetohydrodynamic (MHD) flows with velocity shear and strong magnetic field. Collisionless or weakly collisional plasma is described within the 16-momentum MHD fluid closure model that takes into account not only the effect of pressure anisotropy but also the effect of anisotropic heat fluxes. In this model, the low frequency acoustic wave is revealed into a standard acoustic mode and higher frequency fast thermo-acoustic and lower frequency slow thermo-acoustic waves. It is shown that thermo-acoustic waves become unstable and grow exponentially when the heat flux parameter exceeds some critical value. It seems that velocity shear makes thermo-acoustic waves overstable even at subcritical heat flux parameters. Thus, when the effect of heat fluxes is not profound acoustic waves will grow due to the velocity shear, while at supercritical heat fluxes the flow reveals compressible thermal instability. Anisotropic thermal instability should be also important in astrophysical environments, where it will limit the maximal value of magnetic field that a low density ionized anisotropic flow can sustain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗