SEARCH · Search NASA

Results for “stabilization method”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Voltage stability smart meter for analyzing voltage data and controlling an electrical power source and/or an electric appliance

Systems and methods for voltage stability monitoring and active/reactive power support are disclosed herein. In some embodiments, a smart electric meter of an end user in a grid power system can measure the voltage supplied to the end user via the grid power system, and can analyze the voltage data to detect critical voltage characteristics. The critical voltage characteristics may indicate that a voltage collapse event is likely. The smart electric meter can further estimate a voltage stability margin based on the voltage data. If necessary, the smart electric meter can control an electrical power source and/or an electric appliance positioned at or near the end user to increase the voltage stability margin.

Min, Liang↗

DOE OSTI2023

The Legendre Polynomial Axial Expansion Method

This work presents a new formulation of the axial expansion transport method explicitly using Legendre polynomials for arbitrarily high-order expansions. This new formulation also features an alternative method of axial leakage calculation to allow for nonextruded flat source region meshes. This alternative axial leakage is introduced alongside a balance equation requirement to ensure that neutron balance is preserved in the coarse mesh for a given axial leakage formulation, which allows for effective coarse mesh finite difference acceleration. A matrix exponential table method is derived to allow for fast computations of arbitrarily high-order matrix exponentials for this work and precludes the need for further research into matrix exponential calculations for this method. Numerical results are presented that demonstrate the stability of the axial expansion method in systems with voidlike regions, showcase the speedup from matrix exponential tables, and investigate the axial convergence of the method in terms of both expansion order and mesh size.

Herring, Nicholas↗

DOE OSTI2022

High Order Strong Stability Preserving MultiDerivative Implicit and IMEX Runge--Kutta Methods with Asymptotic Preserving Properties

In this article we present a class of high order unconditionally strong stability preserving (SSP) implicit two-derivative Runge--Kutta schemes and SSP implicit-explicit (IMEX) multi-derivative Runge--Kutta schemes where the time-step restriction is independent of the stiff term. The unconditional SSP property for a method of order $p>2$ is unique among SSP methods and depends on a backward-in-time assumption on the derivative of the operator. We show that this backward derivative condition is satisfied in many relevant cases where SSP IMEX schemes are desired. We devise unconditionally SSP implicit Runge--Kutta schemes of order up to $p=4$ and IMEX Runge--Kutta schemes of order up to $p=3$. For the multiderivative IMEX schemes, we also derive and present the order conditions, which have not appeared previously. The unconditional SSP condition ensures that these methods are positivity preserving, and we present sufficient conditions under which such methods are also asymptotic preserving when applied to a range of problems, including a hyperbolic relaxation system, the Broadwell model, and the Bhatnagar--Gross--Krook kinetic equation. We present numerical results to support the theoretical results on a variety of problems.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Method of making a distributed optical fiber sensor having enhanced Rayleigh scattering and enhanced temperature stability, and monitoring systems employing same

A method of making an optical fiber sensor device for distributed sensing includes generating a laser beam comprising a plurality of ultrafast pulses, and focusing the laser beam into a core of an optical fiber to form a nanograting structure within the core, wherein the nanograting structure includes a plurality of spaced nanograting elements each extending substantially parallel to a longitudinal axis of optical fiber. Also, an optical fiber sensor device for distributed sensing includes an optical fiber having a longitudinal axis, a core, and a nanograting structure within the core, wherein the nanograting structure includes a plurality of spaced nanograting elements each extending substantially parallel to the longitudinal axis of the optical fiber. Also, a distributed sensing method and system and an energy production system that employs such an optical fiber sensor device.

42 ENGINEERING↗

DOE OSTI2022

Method of making a distributed optical fiber sensor having enhanced Rayleigh scattering and enhanced temperature stability, and monitoring systems employing same

Chen, Peng Kevin↗

DOE OSTI2018

Methods of computing steady-state voltage stability margins of power systems

In steady-state voltage stability analysis, as load increases toward a maximum, conventional Newton-Raphson power flow Jacobian matrix becomes increasingly ill-conditioned so power flow fails to converge before reaching maximum loading. A method to directly eliminate this singularity reformulates the power flow problem by introducing an AQ bus with specified bus angle and reactive power consumption of a load bus. For steady-state voltage stability analysis, the angle separation between the swing bus and AQ bus can be varied to control power transfer to the load, rather than specifying the load power itself. For an AQ bus, the power flow formulation is only made up of a reactive power equation, thus reducing the size of the Jacobian matrix by one. This reduced Jacobian matrix is nonsingular at the critical voltage point, eliminating a major difficulty in voltage stability analysis for power system operations.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2023

A Stability Analysis Tool for Bulk Power Systems Using Black-Box Models of Inverter-Based Resources

This paper presents a small-signal stability analysis tool for large-scale power systems with high penetration of inverter-based resources (IBRs). Firstly, a network transfer function matrix (NTFM), which represents the information of the system topology, transmission lines, loads, IBRs locations, etc., is derived to model the entire power system network. Secondly, small-signal perturbation method is applied to obtain the sequence impedance/admittance models of the block-box IBRs considering the frequency cross-coupling effects. With the obtained NTFM as well as IBR models, a multi-input, multi-output (MIMO) feedback system is constructed, and the generalized Nyquist criterion (GNC)-based stability method is employed to analyze the stability of the entire power system. Furthermore, based on the developed stability analysis method, sensitivity analysis is conducted on an unstable case to identify which parameter has a high impact on the system stability. As a result, different test cases based on a modified IEEE 14- bus system as well as a reduced 240-bus WECC system are studied to verify the proposed stability analysis tool.

42 ENGINEERING↗

DOE OSTI2018

Numerical Method for Fractional Advection-Diffusion Equation with Heredity

We propose a method of construction of difference schemes for fractional partial differential equations with delay in time. For the fractional equation with two-sided diffusion, fractional transfer in time, and a functional aftereffect, we construct an implicit difference scheme. We use the shifted Grünwald–Letnikov formulas for the approximation of fractional derivatives with respect to spatial variables and the L1-algorithm for the approximation of fractional derivatives in time. Also we use piecewise constant interpolation and extrapolation by extending the discrete prehistory of the model in time. The algorithm is a fractional analog of a purely implicit method; on each time step, it is reduced to the solution of linear algebraic systems. We prove the stability of the method and find its order of convergence.

97 MATHEMATICS AND COMPUTING↗

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 OSTI2018

An improved three-point bending test method for the investigation of nanosecond laser dicing of ultrathin Si dies with Cu stabilization layer

Highlights: • A selective chemical etching method has been developed for removal of metallization layers on ultrathin Si dies for three-point bending test. • The complex elastoplastic stress distribution in the die is simplified to an elastic stress distribution. • The die sidewall fracture strength can be determined using the standard three-point bending tensile stress equation. • The characteristic fracture strengths of ultrathin Si dies were correlated to the laser-induced microstructures of the die sidewall. - Abstract: The sidewall fracture strength of ultrathin Si dies with frontside and backside metallization layers cannot be determined directly from three-point bending (3PB) test. The complex 3PB elastoplastic stress distribution makes it difficult to determine the die sidewall fracture strength by analytical method or finite element analysis due to various theoretical and practical constraints. In this work, a selective chemical etching method has been developed for complete removal of the frontside and backside metallization layers on the ultrathin die without affecting the critical microstructure and strength of the die sidewall induced by laser dicing. By removing all the metallization layers on the die, the complex 3PB elastoplastic stress distribution is simplified to an elastic stress distribution. This enables a straightforward calculation of the die sidewall fracture strength using the standard 3PB tensile stress equation for a homogeneous linear elastic material with a rectangular cross section. Two chemical etchants were evaluated and found to be suitable but FeCl{sub 3} is preferred because of the much faster etching time. With this improved 3PB test method, the effect of nanosecond laser dicing on ultrathin Si dies with backside Cu layer was investigated. The die frontside and backside characteristic fracture strengths of ultrathin Si dies with varying backside Cu thickness were determined and correlated to TEM analyses of the compositions and microstructures at the fracture origins.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Automatically adaptive stabilized finite elements and continuation analysis for compaction banding in geomaterials

Under compressive creep, viscoplastic solids experiencing internal mass transfer processes can accommodate singular cnoidal wave solutions as material instabilities at the stationary wave limit. These instabilities appear when the loading rate is significantly faster than the material's capacity to diffusive internal perturbations, leading to localized failure features (e.g., cracks and compaction bands). These cnoidal waves, generally found in fluids, have strong nonlinearities that produce periodic patterns. Due to the singular nature of the solutions, the applicability of the theory is currently limited. Additionally, practical simulation tools require proper regularization to overcome the challenges that singularity induces. We focus on the numerical treatment of the governing equation using a nonlinear approach building on a recent adaptive stabilized finite element method. This automatic refinement method provides an error estimate that drives mesh adaptivity, a crucial feature for the problem at hand. We compare the performance of this adaptive strategy against analytical and standard finite element solutions. We then investigate the sensitivity of the diffusivity ratio, the parameter controlling the process, and identify multiple possible solutions with several stress peaks. We also show the evolution of the spacing between peaks for all solutions as a function of that parameter.

42 ENGINEERING↗

DOE OSTI2025

An automated integrated web-based smart tool for open stope design

The Stability Graph is a widely used tool for the design of open stopes in underground mining. Many users of the Stability Graph still apply this design method manually. Although the manual approach has benefits, using multiple graphs and stability number computation charts for each stope surface is time-consuming, even for the experienced mining engineer. Current practice in the use of the method also limits data sharing. This paper presents a StopeSoft web-based tool for open stope stability prediction that is developed on the basis of the Stability Graph method and is available at openstope.com. StopeSoft incorporates flexibility in terms of Stability Graph options and incorporates additional critical factors often overlooked. As a web-based tool, StopeSoft encourages and makes data sharing possible globally, focused on expanding the database and improving the current limitations of the Stability Graph to provide practical, reliable solutions for mining engineers, consultants, and academics. The StopeSoft automated process facilitates the process of open stope stability prediction, saving time and minimizing potential human errors. Statistical treatment of the data accounts for the variability of input parameters to emphasize the probabilistic nature of the Stability Graph method. The probabilistic interpretation of the stability states of stope surfaces eliminates the false feeling of absolute stope performance based on its location on the Stability Graph , as implied by the deterministic approach.

58 GEOSCIENCES↗

DOE OSTI2024

Impedance Scan of Inverter-Based Resources and Diesel Generator for Stability Analysis: Preprint

Impedance-based methods are widely used for power system stability analysis with inverter-based resources (IBRs), e.g., assessing dynamic interactions between the power grid and an IBR, control interactions between multiple IBRs, and the sub-synchronous oscillation and damping phenomenon. Since it is difficult to get a numerical model 100% matching with the hardware IBR, using the hardware inverter directly to obtain its output impedance has become a prominent approach nowadays. Therefore, this article presents the impedance scan using hardware IBRs, and also a hardware diesel generator as it still stays with the grid before the grid completely goes to renewable. The devices under test (DuTs) for the impedance scan includes two 3-..phi.., 480 V, 60 Hz commercial grid-forming IBRs (one of 250 kVA and another of 125 kVA rating) in series with ..delta..-Y transformers, one 3-..phi.., 480 V, 60 Hz commercial grid-following IBR (of 125 kVA rating), and a 3-..phi.., 480 V, 60 Hz commercial diesel generator (of 187.5 kVA rating). Using voltage signals perturbed with sub-, inter-, and higher harmonic components, and measuring the current response, the positive-sequence impedances are computed via an offline- based post-analysis. Moreover, best-fit transfer functions are estimated that closely resemble the measured data points of the positive-sequence impedances. Based on the observations from various outcomes of the hardware experiments, this article also provides some fundamental insights on the equivalent positive- sequence impedance of a combination of multiple hardware components by comparing the estimated and the empirically computed impedances. A comparative insight on the damping capability of the DuTs using the positive-sequence impedances of the hardware is also discussed.

grid following inverter↗

DOE OSTI2022

Using distributed power electronics-based devices to improve the voltage and frequency stability of distribution systems

Methods include, in response to a line frequency variation of a power grid, adjusting a voltage setpoint of a voltage regulator coupled to the power grid at a grid edge to maintain a voltage at the grid edge, wherein the adjusting the regulated voltage setpoint is configured to reduce the line frequency variation to stabilize the line frequency of the power grid. Apparatus include a voltage regulator configured to couple to a power grid at a grid edge and to maintain a voltage at the grid edge, wherein the voltage regulator is further configured to adjust a voltage setpoint of a voltage regulator in response to a line frequency variation of the power grid to reduce the line frequency variation and stabilize the line frequency of the power grid.

Du, Wei↗

DOE OSTI2023

Role of spin-orbit coupling on crystal-field splitting and phase-stability of rare-earth based layered intermetallic

Layered and incommensurate heterostructures have attracted much attention for the occurrence of superconductivity and charge density waves with the possibility of intercalating foreign atoms. However, lanthanide-based heterostructures where (Eu x Sm 1-x )S and TaS 2 are alternatively stacked have been scarcely investigated. In this work, we performed phase stability, bonding behavior and electronic-structure analysis of (Eu x Sm 1-x )TaS 3 using first-principles density-functional theory methods. Our phase stability analysis suggests 50 at.% solubility of Eu in (Eu x Sm 1-x )TaS 3 compared to Eu solubility in all proportions in cubic SmS. The instability of Eu beyond 50 at.% in (Eu x Sm 1-x )TaS 3 was attributed to higher density of Eu-4f states at the Fermi-level. Based on band position calculated from spin-orbit coupling effect, we constructed a qualitative schematic of possible crystal-field analysis. The local change in bond-length and bond-angle around Eu-S in (Eu x Sm 1-x )TaS 3 correlate well with our crystal-field analysis. We believe that quantum mechanical insights provided in this work will be useful to understand other complex heterostructures.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Numerical multi-scaling method to solve the linear stochastic partial differential equations

In this paper, we propose a multi-scaling method to solve a class of the stochastic partial differential equations (SPDEs) of Itô-type. By this method, we approximate the solution of the SPDEs with strongly elliptic operator on a finite domain. The method is based on B-spline wavelet approximations; some of these functions are reshaped to satisfy on boundary conditions exactly. We prove the consistency and stability of the method. Illustrative examples are included to demonstrate the validity and applicability of our proposed technique.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2018

Cost-effective Manufacturing and Morphological Stabilization of Nanostructured Cathodes for Commercial SOFCs

The overall goals of the program were to develop methods for preparing and stabilizing SOFC electrodes by infiltration that can be manufactured on a practical scale, then transfer those methods to cells in stacks. Cathodes prepared by infiltration of nanoparticles into a porous scaffold of the electrolyte have demonstrated outstanding electrochemical performance and mechanical strength but are difficult to manufacture due to the large number of infiltration steps required to produce a conductive composite and show poor stability due to sintering of the nanoparticles. At Penn, the main goal was to develop a fabrication method that uses only a small number of infiltration steps. Infiltration with molten salts and electrodeposition were pursued but deemed impractical. A method in which an electronically conductive scaffold consisted of a composite of a mixed conductor, LSF (La 0.9 Sr 0.1 FeO 3 ), with yttria-stabilized zirconia (YSZ) was found to be promising. While some cation mixing between LSF and YSZ was observed during calcination, there was no evidence for new phases in XRD and incorporation of Zr into the perovskite could be minimized by reducing the Sr content of the LSF. To stabilize nanoparticles infiltrated into the LSF-YSZ composite scaffold, the USC team examined the effects of Atomic Layer Deposition with ZrO 2 and infiltrated nanoparticles. Significantly enhanced stability was observed with the infiltrated nanoparticles. The FuelCell Energy team prepared 4”x4” cells in which an LSF-YSZ composite scaffold was cofired with the YSZ electrolyte and infiltrated with LSCF nanoparticles. Apparently because the LSF-YSZ scaffold layer was very thin in this case, the cation mixing that occurred during calcination resulted in a poor scaffold-electrolyte interface which caused reduced performance.

36 MATERIALS SCIENCE↗

DOE OSTI2018

An exponential time-integrator scheme for steady and unsteady inviscid flows

In this report an exponential time-integrator scheme of second-order accuracy based on the predictor–corrector methodology, denoted PCEXP, is developed to solve multi-dimensional nonlinear partial differential equations pertaining to fluid dynamics. The effective and efficient implementation of PCEXP is realized by means of the Krylov method. The linear stability and truncation error are analyzed through a one-dimensional model equation. The proposed PCEXP scheme is applied to the Euler equations discretized with a discontinuous Galerkin method in both two and three dimensions. The effectiveness and efficiency of the PCEXP scheme are demonstrated for both steady and unsteady inviscid flows. The accuracy and efficiency of the PCEXP scheme are verified and validated through comparisons with the explicit third-order total variation diminishing Runge–Kutta scheme (TVDRK3), the implicit backward Euler (BE) and the implicit second-order backward difference formula (BDF2). For unsteady flows, the PCEXP scheme generates a temporal error much smaller than the BDF2 scheme does, while maintaining the expected acceleration at the same time. Moreover, the PCEXP scheme is also shown to achieve the computational efficiency comparable to the implicit schemes for steady flows.

97 MATHEMATICS AND COMPUTING↗