Search NASA⌕ Search

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.

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

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

FLORA Equilibrium and Stability Code Archive

FLORA solves, in a 2-D domain (radial and axial dimensions) with a specified azimuthal Fourier mode, for the linearized stability of a long, thin, axisymmetric plasma equilibrium in an applied magnetic field. Before the stability equation is solved, FLORA solves a set of simple equations for pressure balance that specify the equilibrium magnetic field and plasma pressure in the long-thin limit given a simplified description of the magnetic coils. It uses an initial-value method for the linear stability problem in which an equilibrium is given an initial perturbation to its magnetic field, and the temporal behavior of the perturbation is followed. The perturbation has been Fourier expanded in the azimuthal coordinate; each azimuthal mode must be examined separately. The complex partial differential equation of motion for the perturbed radial displacement of the field lines is solved as a coupled system of two real p.d.e.'s and the solution consists of two parts, the real part and the imaginary part. The system is solved by bringing the coupling terms in each equation to the right side and using an iterative technique.

Cohen, Bruce [Lawrence Livermore National Laborato↗
DOE OSTI2026

Stabilizing Cathode–Electrolyte Interphase of Nickel-Rich Single-Crystal Cathodes for Lithium-Ion Batteries

Nickel-rich single-crystal (SC) layered oxides are promising cathode candidates for next-generation lithium-ion batteries (LIBs) owing to their high energy density and structural robustness against intergranular cracking. However, their intrinsic surface reactivity with liquid electrolytes accelerates parasitic reactions at the cathode–electrolyte interphase (CEI), leading to transition-metal dissolution, gas generation, and impedance buildup. In this work, we synthesized SC-Li x Ni 0.9 Mn 0.05 Co 0.05 O 2 (NMC9055, 1 ≤ x ≤ 1.2) using a eutectic-assisted method and investigated interface stabilization strategies. A nickel-deficient Li x Ni 0.6 Mn 0.2 Co 0.2 O 2 (NMC622, 1 ≤ x ≤ 1.2) coating was applied via evaporation-based deposition to suppress CEI degradation pathways. Structural and compositional analyses confirmed uniform shell formation and preserved particle integrity. Half-cell electrochemical testing against lithium metal revealed ∼10% higher capacity retention and improved reversibility compared with pristine SC NMC9055, particularly under high-voltage operation. In conclusion, these results highlight the critical role of controlled surface chemistry in mitigating CEI instability in nickel-rich SC cathodes, offering a pathway toward enabling durable high-energy LIBs.

Coating processes↗
DOE OSTI2022

Islet: interpolation semi-Lagrangian element-based transport

Abstract. Advection of trace species, or tracers, also called tracer transport, in models of the atmosphere and other physical domains is an important and potentially computationally expensive part of a model's dynamical core. Semi-Lagrangian (SL) advection methods are efficient because they permit a time step much larger than the advective stability limit for explicit Eulerian methods without requiring the solution of a globally coupled system of equations as implicit Eulerian methods do. Thus, to reduce the computational expense of tracer transport, dynamical cores often use SL methods to advect tracers. The class of interpolation semi-Lagrangian (ISL) methods contains potentially extremely efficient SL methods. We describe a finite-element ISL transport method that we call the interpolation semi-Lagrangian element-based transport (Islet) method, such as for use with atmosphere models discretized using the spectral element method. The Islet method uses three grids that share an element grid: a dynamics grid supporting, for example, the Gauss–Legendre–Lobatto basis of degree three; a physics parameterizations grid with a configurable number of finite-volume subcells per element; and a tracer grid supporting use of Islet bases with particular basis again configurable. This method provides extremely accurate tracer transport and excellent diagnostic values in a number of verification problems.

58 GEOSCIENCES↗
DOE OSTI2022

Process window estimation in manufacturing through Entropy-Sigma active learning

In manufacturing, there exist boundary identification problems for defining parameter spaces that meet desired thresholds on outcomes. This paper presents an Entropy-Sigma acquisition function for active learning of the process window/map in manufacturing using a Gaussian Process surrogate. Here, the method is applied to identify the stability boundary for the stability process map in machining using time-domain simulations with a periodic sampling stability metric. Results show that the proposed Entropy-Sigma method significantly outperforms Latin hypercube sampling or grid-based methods. The described method can be applied to identify the process window/map for any manufacturing application using a quantitative process outcome metric.

42 ENGINEERING↗
DOE OSTI2025

First‐Order Empirical Interpolation Method for Real‐Time Solution of Parametric Time‐Dependent Nonlinear PDEs

ABSTRACT We present a model reduction approach for the real‐time solution of time‐dependent nonlinear partial differential equations (PDEs) with parametric dependencies. A major challenge in constructing efficient and accurate reduced‐order models for nonlinear PDEs is the efficient treatment of nonlinear terms. We address this by unifying the implementation of hyperreduction methods to deal with nonlinear terms. Furthermore, we introduce a first‐order empirical interpolation method (EIM) to provide an efficient approximation of the nonlinear terms in time‐dependent PDEs. We demonstrate the effectiveness of our approach on the Allen–Cahn equation, which models phase separation, and the Buckley–Leverett equation, which describes two‐phase fluid flow in porous media. Numerical results highlight the accuracy, efficiency, and stability of the proposed method compared with both the Galerkin–Newton approach and hyper‐reduced models using the standard EIM.

Nguyen, Ngoc Cuong [Center for Computational Engin↗
DOE OSTI2020

F4‐TCNQ as an Additive to Impart Stretchable Semiconductors with High Mobility and Stability

Abstract Numerous strategies are developed to impart stretchability to polymer semiconductors. Although these methods improve the ductility, mobility, and stability of such stretchable semiconductors, they nonetheless still need further improvement. Here, it is shown that 2,3,5,6‐tetrafluoro‐7,7,8,8‐tetracyanoquinodimethane (F4‐TCNQ) is an effective molecular additive to tune the properties of a diketopyrrolopyrrole‐based (DPP‐based) semiconductor. Specifically, the addition of F4‐TCNQ is observed to improve the ductility of the semiconductor by altering the polymer’s microstructures and dynamic motions. As a p‐type dopant additive, F4‐TCNQ can also effectively enhance the mobility and stability of the semiconductor through changing the host polymer’s packing structures and charge trap passivation. Upon fabricating fully stretchable transistors with F4‐TCNQ‐DPP blended semiconductor films, it is observed that the resulting stretchable transistors possess one of the highest initial mobility of 1.03 cm 2 V −1 s −1 . The fabricated transistors also exhibit higher stability (both bias and environmental) and mobility retention under repeated strain, compared to those without F4‐TCNQ additive. These findings offer a new direction of research on stretchable semiconductors to facilitate future practical applications.

Mun, Jaewan↗
DOE OSTI2025

Spectrally Stabilized Interface Capturing Formulation and Implementation in Nek5000/NekRS

This report documents the formulation of a novel level-set method for incompressible two-phase flows in the continuous Galerkin (CG) high order spectral element framework. The overall method hinges on a novel implementation of the spectral vanishing viscosity (SVV) operator for the stabilization of linear/non-linear hyperbolic problems. The multidimensional SVV convolution kernels, which in essence, have a similar effect as a high pass filter applied to the derivatives, are formulated by exploiting the tensor product form, analogous to the construction of the usual stiffness matrix system. The resulting kernels are directionally decoupled and ensure a linear, symmetric positive definite, elliptic matrix operator. The SVV formulation is demonstrated to provide a robust stabilizing mechanism through challenging linear and non-linear hyperbolic problems, including problems pertinent to the level-set formulation. The two-phase framework conceptualized herein is based on the conservative level-set (CLS) method which represents the interface between the fluids by the 0.5 iso-contour of the smoothed Heaviside function. The CLS method is augmented with a preconditioning procedure for interface normals using the signed distance function which precludes the manifestation of spurious oscillations in the vicinty of the interface. Further, the existing mixed explicit-implicit approach for the solution of Navier-Stokes equations in Nek5000, as described in Tomboulides et al, is augmented with a pressure coefficient splitting approach for the Poisson equation, which greatly accelerated the convergence of pressure solver for two-phase systems with large density ratio. The robustness and accuracy of the overall two-phase method is demonstrated through canonical challenging problems involving high density and viscosity ratios, with and without surface tension. The two-phase formulation is wholly implemented in Nek5000 and the SVV stabilization method is implemented in NekRS, which is the essential precursor to the two-phase framework, undergoing active development.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Resilient Distributed Frequency Regulation of Renewable Generators under Communication Interruptions

Modern power systems (MPSs) face significant challenges due to the high penetration of renewable energy sources (RESs) and new types of loads such as electric vehicles (EVs). Traditional load frequency control (LFC) methods struggle with the intermittent, stochastic nature of RESs, the near-zero inertia of power-electronics-based generators, and the mobility of controllable loads and battery systems. This paper introduces a novel resilient distributed frequency regulation method to address these issues. The proposed method employs a state space model to represent the dynamic behavior of participating power sources while accounting for stochastic switching processes to model structural and parameter variations caused by disruptions such as generator connection/disconnection, communication interruptions, and physical faults. By integrating these dynamic and stochastic components, the method treats power grids as a comprehensive stochastic hybrid system. Our method enhances conventional frequency control by incorporating local stability control, neighborhood control decoupling, and coordination feedback. Theoretical analyses establish the stability, convergence, and resilience of the proposed method, and its effectiveness is validated through case studies.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2019

Extended virtual element method for the Laplace problem with singularities and discontinuities

In this paper, we propose the extended virtual element method (X-VEM) to treat singularities and crack discontinuities that arise in the Laplace problem. The virtual element method (VEM) is a stabilized Galerkin formulation on arbitrary polytopal meshes, wherein the basis functions are implicit (virtual)—they are not known explicitly nor do they need to be computed within the problem domain. In this work, suitable projection operators are used to decompose the bilinear form on each element into two parts: a consistent term that reproduces the first-order polynomial space and a correction term that ensures stability. A similar approach is pursued in the X-VEM with a few notable extensions. To capture singularities and discontinuities in the discrete space, we augment the standard virtual element space with an additional contribution that consists of the product of virtual nodal basis (partition-of-unity) functions with enrichment functions. For discontinuities, basis functions are discontinuous across the crack and for singularities a weakly singular enrichment function that satisfies the Laplace equation is chosen. For the Laplace problem with a singularity, we devise an extended projector that maps functions that lie in the extended virtual element space onto linear polynomials and the enrichment function, whereas for the discontinuous problem, the consistent element stiffness matrix has a block-structure that is readily computed. An adaptive homogeneous numerical integration method is used to accurately and efficiently (no element-partitioning is required) compute integrals with integrands that are weakly singular. Once the element projection matrix is computed, the same steps as in the standard VEM are followed to compute the element stabilization matrix. Numerical experiments are performed on quadrilateral and polygonal (convex and nonconvex elements) meshes for the problem of an $\textit{L}$-shaped domain with a corner singularity and the problem of a cracked membrane under mode III loading, and results are presented that affirm the sound accuracy and demonstrate the optimal rates of convergence in the $\textit{L}^2$ norm and energy of the proposed method.

42 ENGINEERING↗
DOE OSTI2023

CFL Optimized Forward–Backward Runge–Kutta Schemes for the Shallow-Water Equations

Abstract We present the formulation and optimization of a Runge–Kutta-type time-stepping scheme for solving the shallow-water equations, aimed at substantially increasing the effective allowable time step over that of comparable methods. This scheme, called FB-RK(3,2), uses weighted forward–backward averaging of thickness data to advance the momentum equation. The weights for this averaging are chosen with an optimization process that employs a von Neumann–type analysis, ensuring that the weights maximize the admittable Courant number. Through a simplified local truncation error analysis and numerical experiments, we show that the method is at least second-order in time for any choice of weights and exhibits low dispersion and dissipation errors for well-resolved waves. Further, we show that an optimized FB-RK(3,2) can take time steps up to 2.8 times as large as a popular three-stage, third-order strong stability-preserving Runge–Kutta method in a quasi-linear test case. In fully nonlinear shallow-water test cases relevant to oceanic and atmospheric flows, FB-RK(3,2) outperforms SSPRK3 in admittable time step by factors roughly between 1.6 and 2.2, making the scheme approximately twice as computationally efficient with little to no effect on solution quality. Significance Statement The purpose of this work is to develop and optimize time-stepping schemes for models relevant to oceanic and atmospheric flows. Specifically, for the shallow-water equations we optimize for schemes that can take time steps as large as possible while retaining solution quality. We find that our optimized schemes can take time steps between 1.6 and 2.2 times larger than schemes that cost the same number of floating point operations, translating directly to a corresponding speedup. Our ultimate goal is to use these schemes in climate-scale simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Porosity evolution during heating of copper made from powder by friction extrusion

Here, friction extrusion was used to compact and extrude solid copper rod from feedstock powders. Following extrusion, considerable porosity was observed throughout the extrudate cross section due to the entrained porosity in the feedstock material and the extrusion process. The thermal stability of the extrudate was investigated via a series of heat treatments. Porosity evolution exhibited three distinct stages—an unchanged plateau (0°C–300°C) followed by an abrupt increase (400°C–500°C) and ultimately a reduction, as the temperature increases (>500°C). The peak porosity measured was ~25%. The underlying driving force for pore evolution is described as the competition between the internal pore pressure, material strength, and sintering kinetics, as a function of temperature. The observed porosity evolution and driving force are not expected to be limited to copper. Thus, this manuscript reveals an important consideration regarding elevated-temperature thermal stability as advanced manufacturing methods involving direct powder extrusion are explored.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Measurement of Stability Diagram at IOTA at Fermilab

Nonlinear focusing elements can enhance the stability of particle beams in high-energy colliders by means of Landau Damping, through the tune spread which is introduced. We propose an experiment at Fermilab’s Integrable Optics Test Accelerator (IOTA) to investigate the influence of nonlinear focusing elements on the transverse stability of the beam. In this experiment, we employ an anti-damper, an active transverse feedback system, as a controlled mechanism to induce coherent beam instability. By utilizing the anti-damper, we can examine the impact of the nonlinear focusing element on the beam’s transverse stability. The stability diagram, a tool used to determine the system’s stability, will be measured using a recently demonstrated method at the LHC. This measurement is carried out experimentally by selecting specific threshold gains and measuring them for a range of phases. The stability diagram is represented by gei¿ on the complex plane. The experiment at IOTA adds insight towards the stability diagram measurement method by supplying a reduced machine impedance, to investigate the impedance’s effect on the stability diagram, as well as a larger range of phase measurements.

43 PARTICLE ACCELERATORS↗
DOE OSTI2019

Efficient IMEX Runge-Kutta Methods for Nonhydrostatic Dynamics

We develop new implicit-explicit Runge-Kutta (IMEX RK) methods for integrating horizontally explicit, vertically implicit (HEVI) partitionings of nonhydrostatic atmosphere models (HEVI models). These new methods, termed IMKG methods, are IMEX RK methods whose explicit part has optimal or near-optimal stability on the imaginary axis and whose implicit part is I-stable. A specialized stability region is presented for characterizing the stability of IMKG and other IMEX methods integrating HEVI models. Subsequently, we formulate two families of IMEX RK methods to enable deriving IMKG methods with a high explicit stage count for integrating HEVI models with large, stable time-steps. Here, we then derive a HEVI partitioning of the HOMME-NH nonhydrostatic model and use this model to compare the accuracy and efficiency of several IMKG methods with other IMEX RK methods from the literature.

58 GEOSCIENCES↗
DOE OSTI2018

L{sup ∞} error bound of conservative compact difference scheme for the generalized symmetric regularized long-wave (GSRLW) equations

In this paper, we design a compact finite difference scheme which preserves the original conservative properties to solve the generalized symmetric regularized long-wave equations. The existence of the difference solution is proved by the Brouwer fixed-point theorem. Applying the discrete energy method, the convergence and stability of the difference scheme is obtained, and its numerical convergence order is O(τ{sup 2}+h{sup 4}) in the L{sup ∞}-norm for u and ρ. For computing the nonlinear algebraic system generated by the compact scheme, a decoupled iterative algorithm is constructed and proved to be convergent. Numerical experiment results show that the theory is accurate and the method is efficient and reliable.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A Cyber Threat Mitigation Approach For Wide Area Control of SVCs using Stability Monitoring

We propose a stability monitoring approach for the mitigation of cyber threats directed at the wide area control (WAC) system used for coordinated control of Flexible AC Transmission Systems (FACTS) used for power oscillation damping (POD) of active power flow on inter-area tie lines. The approach involves monitoring the modes of the active power oscillation on an inter-area tie line using the Matrix Pencil (MP) method. We use the stability characteristics of the observed modes as a proxy for the presence of destabilizing cyber threats. We monitor the system modes to determine whether any destabilizing modes appear after the WAC system engages to control the POD. If the WAC signal exacerbates the POD performance, the FACTS falls back to POD using local measurements. The proposed approach does not require an expansive system-wide view of the network. We simulate replay, control integrity, and timing attacks for a test system and present results that demonstrate the performance of the SM approach for mitigation.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Large-signal Stability Analysis of Grid-forming Inverters with Equivalent-circuit Models

Here, this paper proposes an energy function-based direct method for large-signal stability assessment of grid-forming (GFM) inverters leveraging an equivalent-circuit representation of all involved control- and physical-layer dynamics. Three different primary controls, a standard inner-current outer-voltage cascaded-control architecture, output LCL filter, and reference-current saturation limiting are featured in the modeling and analysis framework. A composite energy function for the GFM inverter is obtained by summing up individual energy contributions gleaned from the circuit representation. The approach can readily be generalized to different primary controls, output-filter arrangements, and current limiters since it is based on a circuit-theoretic foundation. Numerical simulations validate the efficacy of the approach in estimating the critical clearing time following a large-signal disturbance.

24 POWER TRANSMISSION AND DISTRIBUTION↗