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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

On a Method of Study of Specific Asymptotic Stability of Solutions to a Sixth-Order Linear Integrodifferential Volterra Equation

We state sufficient conditions of the asymptotic stability on the semi-axis of solutions to a linear, homogeneous, sixth-order integrodifferential Volterra-type equation in the case where solutions of the corresponding linear, homogeneous, sixth-order differential equation are asymptotically unstable. We also present a new method and an illustrating example.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Collocation Methods for Exploring Perturbations in Linear Stability Analysis

Eigenvalue analysis is a well-established tool for stability analysis of dynamical systems. However, there are situations where eigenvalues miss some important features of physical models. For example, in models of incompressible fluid dynamics, there are examples where linear stability analysis predicts stability but transient simulations exhibit significant growth of infinitesimal perturbations. This behavior can be predicted by pseudospectral analysis. In this study, we show that an approach similar to pseudospectral analysis can be performed inexpensively using stochastic collocation methods and the results can be used to provide quantitative information about instability.In addition, we demonstrate that the results of the perturbation analysis provide insight into the behavior of unsteady flow simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Method and apparatus for implementing ultra-high stability stages with combined degrees of freedom for multiple axes of motion

A method and apparatus are provided for implementing an ultra-high stability stage with combined degrees of freedom for multiple axes of motion. The ultra-high stability stage includes a base, a driving wedge supported by the base and a following wedge supported by the driving wedge. The base and each wedge are formed of a selected stable material having predefined rigidity and low thermal expansion coefficient. Integrated air bearings, respective driving mechanics associated with each of the wedges and guiding components having selected degrees of freedom enable movement about multiple axes of motion, such as X, Y, Z translation axes, and rotation X and rotation Y axes. Another ultra-high stability stage with combined degrees of freedom for multiple axes of motion includes an intermediate wedge between the driving wedge and the following wedge to enable additional movement about a rotation X axis.

Preissner, Curt Alfred↗
DOE OSTI2023

Bound-preserving finite element approximations of the Keller–Segel equations

We report this paper aims to develop numerical approximations of the Keller–Segel equations that mimic at the discrete level the lower bounds and the energy law of the continuous problem. We solve these equations for two unknowns: the organism (or cell) density, which is a positive variable, and the chemoattractant density, which is a non-negative variable. We propose two algorithms, which combine a stabilized finite element method and a semi-implicit time integration. The stabilization consists of a nonlinear artificial diffusion that employs a graph-Laplacian operator and a shock detector that localizes local extrema. As a result, both algorithms turn out to be nonlinear and can generate cell and chemoattractant numerical densities fulfilling lower bounds. However, the first algorithm requires a suitable constraint between the space and time discrete parameters, whereas the second one does not. We design the latter to attain a discrete energy law on acute meshes. We report some numerical experiments to validate the theoretical results on blowup and nonblowup phenomena. In the blowup setting, we identify a locking phenomenon that relates the L ∞ (Ω)-norm to the L 1 (Ω)-norm limiting the growth of the singularity when supported on a macroelement.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Voltage Stability Analysis of Power Systems With Induction Motors Based on Holomorphic Embedding

Load characteristics have substantial influence on the voltage stability of power systems. The self-restorative characteristic and stalling of induction motor loads can deteriorate the voltage stability, so it is necessary to develop accurate and efficient dynamic analysis methods for voltage stability analysis of systems with induction motor loads. In this paper, a set of methods based on holomorphic embedding is proposed, which is able to solve steady states and dynamics of a power system with induction motors. The paper is the first work that applies holomorphic embedding to dynamic problems of power systems. Here, the test cases on the IEEE 14-bus system, the NPCC 140-bus system, and the Polish 2383-bus system verify both accuracy and efficiency of the proposed methods.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Deep Sigma Point Processes-Assisted Chance-Constrained Power System Transient Stability Preventive Control

Here this paper proposes a deep sigma point processes (DSPP)-assisted chance-constrained power system transient stability preventive control method to deal with uncertain renewable energy and loads-induced stability risk. The traditional transient stability-constrained preventive control is reformulated as a chance-constrained optimization problem. To deal with the computational bottleneck of the time-domain simulation-based probabilistic transient stability assessment, the DSPP is developed. DSPP is a parametric Bayesian approach that allows us to predict system transient stability with high computational efficiency while accurately quantifying the confidence intervals of the predictions that can be used to inform system instability risk. To this end, with a given preset confidence probability, we embed DSPP into the primal dual interior point method to help solve the chance-constrained preventive control problem, where the corresponding Jacobian and Hessian matrices are derived. Comparison results with other existing methods show that the proposed method can significantly speed up preventive control while maintaining high accuracy and convergence

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Physics-based stabilized finite element approximations of the Poisson–Nernst–Planck equations

We present and analyze two stabilized finite element methods for solving numerically the Poisson–Nernst–Planck equations. The stabilization we consider is carried out by using a shock detector and a discrete graph Laplacian operator for the ion equations, whereas the discrete equation for the electric potential need not be stabilized. Discrete solutions stemmed from the first algorithm preserve both maximum and minimum discrete principles. For the second algorithm, its discrete solutions are conceived so that they hold discrete principles and obey an entropy law provided that an acuteness condition is imposed for meshes. Remarkably the latter is found to be unconditionally stable. We validate our methodology through transient numerical experiments that show convergence toward steady-state solutions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

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

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.

current measurement↗
DOE OSTI2022

Liquid bicontinuous stabilized assemblies by interfacial jamming of nanoparticles

A method for preparing a stabilized assembly includes combining a first liquid phase including nanoparticles and a second, immiscible liquid phase, dissolving in the second phase a first end-functionalized polymer having a first molecular weight, and a second end-functionalized polymer having a second molecular weight, wherein the second molecular weight is greater than the first molecular weight, applying a shearing external deformation field to increase the surface area of the first phase to create a new interface, wherein the nanoparticle surfactants form a disordered, jammed assembly at the new interface, and releasing the shearing external deformation field. The polymer and the nanoparticles can interact at an interface through ligand interactions to form nanoparticle surfactants and upon releasing the external deformation field the jammed assembly at the new interface traps the first phase in a deformed state having the first liquid phase and the second liquid phase as interpenetrating domains.

Huang, Caili↗
DOE OSTI2022

Augmented reduced order models for turbulence

The authors introduce an augmented-basis method (ABM) to stabilize reduced-order models (ROMs) of turbulent incompressible flows. The method begins with standard basis functions derived from proper orthogonal decomposition (POD) of snapshot sets taken from a full-order model. These are then augmented with divergence-free projections of a subset of the nonlinear interaction terms that constitute a significant fraction of the time-derivative of the solution. The augmenting bases, which are rich in localized high wavenumber content, are better able to dissipate turbulent kinetic energy than the standard POD bases. Several examples illustrate that the ABM significantly out-performs L 2 -, H 1 - and Leray-stabilized POD ROM approaches. The ABM yields accuracy that is comparable to constraint-based stabilization approaches yet is suitable for parametric model-order reduction in which one uses the ROM to evaluate quantities of interests at parameter values that differ from those used to generate the full-order model snapshots. Several numerical experiments point to the importance of localized high wavenumber content in the generation of stable, accurate, and efficient ROMs for turbulent flows.

Kaneko, Kento↗
DOE OSTI2022

Hybrid large-eddy simulation with adaptive downstream anisotropic eddy viscosity model

Turbulence pervades most flows of engineering interest, and its prediction remains a challenge on both accuracy and cost. One promising predictive approach that reduces cost combines large-eddy simulation (LES) with simulation based on Reynolds averaged Navier–Stokes equations (RANS). This study presents a method to overcome stability and accuracy issues associated with these hybrid LES-RANS methods. The method developed involves extracting the Reynolds stresses from the upstream LES solution and then using that information to convert the downstream RANS model from a scalar eddy-viscosity model to an anisotropic nonlinear eddy-viscosity model. The method developed differs from the downstream tensorial eddy-viscosity model by being independent of the coordinate system. The method developed was evaluated by computing film cooling of a flat plate with the coolant injected through one row of circular holes. Results obtained show instabilities at the LES-to-RANS interface to be eliminated. Also, the method developed yielded solutions that compare reasonably well with those from LES, even though a significant portion of the flow is computed by the adapted anisotropic RANS model instead of LES, which significantly reduced the number of grid points and computational time needed. Since modification of the downstream RANS model is based on information extracted from the upstream LES solution, the method developed is adaptive to the problem being studied.

42 ENGINEERING↗
DOE OSTI2019

Analysis and Entropy Stability of the Line-Based Discontinuous Galerkin Method

We develop a discretely entropy-stable line-based discontinuous Galerkin method for hyperbolic conservation laws based on a flux differencing technique. By using standard entropy-stable and entropy-conservative numerical flux functions, this method guarantees that the discrete integral of the entropy is non-increasing. This nonlinear entropy stability property is important for the robustness of the method, in particular when applied to problems with discontinuous solutions or when the mesh is under-resolved. This line-based method is significantly less computationally expensive than a standard DG method. Numerical results are shown demonstrating the effectiveness of the method on a variety of test cases, including Burgers’ equation and the Euler equations, in one, two, and three spatial dimensions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Structure preserving and energy dissipative contact approaches for implicit dynamics

In this work, several structure preserving and energy dissipative contact approaches are proposed and evaluated. The time integration schemes considered are general with regard to the version of constraint type, but here the emphasis was on mortar contact. The proposed mortar contact approach conserves both linear and angular momentum for mortar contact in a novel way. The proposed time integration scheme can conserve energy or provide strict contact dissipation. In addition, the proposed scheme enforces both gap constraints and gap velocity constraints (i.e., persistency). Using a midstep time integrator often causes energy dissipation during initial impact, here this energy can be recovered in a novel way. Enforcing the gap velocity constraint mitigates the contact chatter of the contact pressure and nodes in many problems. Whereas some approaches enforce the persistency condition and gap constraints simultaneously during the solution of the equations of motion (EOM) requiring multipliers for both constraints included in the equation set, here the gap constraint is solved through the equations of motion and the persistency condition is satisfied in the time integration scheme by the velocity update after the equations of motion. It is shown that this approach is strictly dissipative in that a plastic contact condition can be achieved. Analogous to a coefficient of restitution for rigid bodies, any dissipated energy can then be returned upon release if energy conservation is desired. Structure preserving methods are good for long-time dynamics simulations and energy conserving and strictly dissipative methods can overcome stability issues associated with standard time integration algorithms such as the Newmark method.

42 ENGINEERING↗
DOE OSTI2023

Probing phase stability in CrMoNbV using cluster expansion method, CALPHAD calculations and experiments

High entropy alloys (HEA), a novel class of materials with multiple principal elements, allow us to tune the properties within a wide composition range where the solid-solution phase is stable. However, it is difficult to experimentally determine the boundaries of the single-phase region. In this paper, we calculate the phase diagram of the CrMoNbV quaternary system, as well as its constituent binary and ternary subsystems, using a combination of the cluster expansion method and CALPHAD calculations. We further verify these results by carrying out experiments at two different compositions of the full quaternary system. Furthermore, our work shows how these computational tools can enable efficient discovery and development of new HEAs.

36 MATERIALS SCIENCE↗