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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
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 OSTI2018

Solving numerically nonlinear systems of balance laws by multivariate sigmoidal functions approximation

A numerical method for solving systems of nonlinear one-dimensional balance laws, based on multivariate sigmoidal functions approximation, is developed. Constructive approximation theorems are first established in both, uniform and L{sup p} norms. A priori as well as a posteriori error estimates are derived for the numerical solutions, when various kinds of sigmoidal functions, such as unit step, logistic and Gompertz functions, are chosen. The residual of the numerical method is also estimated. Numerical examples are given to test the performance of the algorithm, a comparison with the Godunov method is made concerning accuracy and computational cost. Finally, the numerical stability of the method is analyzed.

97 MATHEMATICS AND COMPUTING↗
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↗
DOE OSTI2021

Stabilized copper selenide thermoelectric materials and methods of fabrication thereof

A thermoelectric composition is provided that includes a nanocomposite comprising a copper selenide (Cu 2 Se) matrix having a plurality of nanoinclusions comprising copper metal selenide (CuMSe 2 ) distributed therein. M may be selected from the group consisting of: indium (In), aluminum (Al), gallium (Ga), antimony (Sb), bismuth (Bi), and combinations thereof. The thermoelectric composition has an average figure of merit (ZT) of greater than or equal to about 1.5 at a temperature of less than or equal to about 850K (about 577° C.). Methods of making such a thermoelectric nanocomposite material by a sequential solid-state transformation of a CuSe 2 precursor are also provided.

Poudeu-Poudeu, Pierre Ferdinand↗
DOE OSTI2022

Stabilized electrodes for ion batteries and methods of making the same

The present disclosure relates to a silicon nanoparticle that includes a core of silicon having an outer surface that includes *SiH3-xRx, where *Si is a silicon atom on the outer surface of the core, the first layer covers at least a portion of the outer surface, R is a ligand that includes at least one of —O—R′, —C—R′, —N—R′, —Si—R′, and/or —S—R′, and R′ is a functional group that includes at least one of carbon, oxygen, nitrogen, hydrogen, sulfur, and/or phosphorus.

Neale, Nathan Richard↗
DOE OSTI2022

Stabilized copper selenide thermoelectric materials and methods of fabrication thereof

A thermoelectric composition is provided that includes a nanocomposite comprising a copper selenide (Cu 2 Se) matrix having a plurality of nanoinclusions comprising copper metal selenide (CuMSe 2 ) distributed therein. M may be selected from the group consisting of: indium (In), aluminum (Al), gallium (Ga), antimony (Sb), bismuth (Bi), and combinations thereof. The thermoelectric composition has an average figure of merit (ZT) of greater than or equal to about 1.5 at a temperature of less than or equal to about 850K (about 577° C.). Methods of making such a thermoelectric nanocomposite material by a sequential solid-state transformation of a CuSe2 precursor are also provided.

Poudeu-Poudeu, Pierre Ferdinand↗
DOE OSTI2019

Enhancing Power System Transient Stability by Virtual Synchronous Generator Control Using Wide-Area Measurements

Virtual synchronous generator control provides a potential solution to connect power electronics interfaced sources to the power system. With its flexible control, the emulated inertia may be adaptively changed to enhance system small-signal stability. However, previously proposed methods may have adverse impact to system transient stability. This paper proposes to improve the methods by using wide-area measurements, namely center of inertia frequency. Experimental results are presented to both demonstrate the limitations of previous control schemes, and the effectiveness of the proposed one. The practical implementation of this method is discussed.

Ma, Yiwei↗