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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Visualizing Stereodynamics in Cold Collisions through Shape Resonance Wavefunctions

Shape resonances in cold collisions are often strongly affected by stereodynamics. These resonances can sometimes be assigned by single-channel quantization along the scattering coordinate. However, sensitive steric control of collision implies strong anisotropy of the interaction potential energy surface, which usually leads to coupling among multiple scattering channels in the strongly interacting region. Hence, the resonances might be the result of quantization in a multidimensional space. Using the cold rotationally inelastic collision between para-H 2 (v 1 = 1, j 1 = 2) and HF (v 2 = 0, j 2 = 0) as an example, we analyze four low-lying shape resonances via diagonalization of the full-dimensional Hamiltonian with a stabilization method. While some resonances can indeed be assigned with a single partial wave, others apparently involve more than one scattering channel. Furthermore, a new model based on these resonance wavefunctions is developed to better understand the stereodynamics through shape resonances in cold collisions.

Collisions↗
DOE OSTI2019

Remedial Treatment of Corroded Iron Objects by Environmental Aeromonas Isolates

Using bacteria to transform reactive corrosion products into stable compounds represents an alternative to traditional methodsployed in iron conservation. Two environmental Aeromonas strains (CA23 and CU5) were used to transform ferric iron corrosion products (goethite and lepidocrocite) into stable ferrous iron-bearing minerals (vivianite and siderite). A genomic and transcriptomic approach was used to analyze the metabolic traits of these strains and to evaluate their pathogenic potential. Although genes involved in solid-phase iron reduction were identified, key genes present in other environmental iron-reducing species are missing from the genome of CU5. Several pathogenicity factors were identified in the genomes of both strains, but none of these was expressed under iron reduction conditions. Additional in vivo tests showed hemolytic and cytotoxic activities for strain CA23 but not for strain CU5. Both strains were easily inactivated using ethanol and heat. Nonetheless, given a lesser potential for a pathogenic lifestyle, CU5 is the most promising candidate for the development of a bio-based iron conservation method stabilizing iron corrosion. Based on all the results, a prototype treatment was established using archaeological items. On those, the conversion of reactive corrosion products and the formation of a homogenous layer of biogenic iron minerals were achieved. This study shows how naturally occurring microorganisms and their metabolic capabilities can be used to develop bio-inspired solutions to the problem of metal corrosion.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Camellia: A Rapid Development Framework for Finite Element Solvers

We report the discontinuous Petrov–Galerkin (DPG) methodology of Demkowicz and Gopalakrishnan guarantees the optimality of the finite element solution in a user-controllable energy norm, and provides several features supporting adaptive schemes. The approach provides stability automatically; there is no need for carefully derived numerical fluxes (as in DG schemes) or for mesh-dependent stabilization terms (as in stabilized methods). In this paper, we focus on features of Camellia that facilitate implementation of new DPG formulations; chief among these is a rich set of features in support of symbolic manipulation, which allow, e.g., bilinear formulations in the code to appear much as they would on paper. Many of these features are general in the sense that they can also be used in the implementation of other finite element formulations. In fact, because DPG’s requirements are essentially a superset of those of other finite element methods, Camellia provides built-in support for most common methods. We believe, however, that the combination of an essentially “hands-free” finite element methodology as found in DPG with the rapid development features of Camellia are particularly winsome, so we focus on use cases in this class. In addition to the symbolic manipulation features mentioned above, Camellia offers support for one-irregular adaptive meshes in 1D, 2D, 3D, and space-time. It provides a geometric multigrid preconditioner particularly suited for DPG problems, and supports distributed parallel execution using MPI. For its load balancing and distributed data structures, Camellia relies on packages from the Trilinos project, which simplifies interfacing with other computational science packages. Camellia also allows loading of standard mesh formats through an interface with the MOAB package. Camellia includes support for static condensation to eliminate element-interior degrees of freedom locally, usually resulting in substantial reduction of the cost of the global problem. Finally, we include a discussion of the variational formulations built into Camellia, with references to those formulations in the literature, as well as an MPI performance study.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Structure preserving transport stabilized compatible finite element methods for magnetohydrodynamics

Here, we present compatible finite element space discretizations for the ideal compressible magnetohydrodynamic equations. The magnetic field is considered both in div- and curl-conforming spaces, leading to a strongly or weakly preserved zero-divergence condition, respectively. The equations are discretized in space such that transfers between the kinetic, internal, and magnetic energies are consistent, leading to a preserved total energy. We also discuss further adjustments to the discretization required to additionally achieve magnetic helicity preservation. Finally, we describe new transport stabilization methods for the magnetic field equation which maintain the zero-divergence and energy conservation properties, including one method which also preserves magnetic helicity. The methods' preservation and improved stability properties are confirmed numerically using a steady state and a magnetic dynamo test case.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Downwinding for preserving strong stability in explicit integrating factor Runge-Kutta methods

Here strong stability preserving (SSP) Runge-Kutta methods are desirable when evolving in time problems that have discontinuities or sharp gradients and require nonlinear non-inner-product stability properties to be satisfied. Unlike the case for L 2 linear stability, implicit methods do not significantly alleviate the time-step restriction when the SSP property is needed. For this reason, when handling problems with a linear component that is stiff and a nonlinear component that is not, SSP integrating factor Runge-Kutta methods may offer an attractive alternative to traditional time-stepping methods. The strong stability properties of integrating factor Runge-Kutta methods where the transformed problem is evolved with an explicit SSP Runge-Kutta method with non-decreasing abscissas was recently established. However, these methods typically have smaller SSP coefficients (and therefore a smaller allowable time-step) than the optimal SSP Runge-Kutta methods, which often have some decreasing abscissas. In this work, we consider the use of downwinded spatial operators to preserve the strong stability properties of integrating factor Runge-Kutta methods where the Runge-Kutta method has some decreasing abscissas. We present the SSP theory for this approach and present numerical evidence to show that such an approach is feasible and performs as expected. However, we also show that in some cases the integrating factor approach with explicit SSP Runge-Kutta methods with non-decreasing abscissas performs nearly as well, if not better, than with explicit SSP Runge-Kutta methods with downwinding. In conclusion, while the downwinding approach can be rigorously shown to guarantee the SSP property for a larger timestep, in practice using the integrating factor approach by including downwinding as needed with optimal explicit SSP Runge-Kutta methods does not necessarily provide significant benefit over using explicit SSP Runge-Kutta methods with non-decreasing abscissas.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Materials for stabilizing semiconductors and methods of making the same

The present disclosure relates to a device that includes a first layer having an active material and a stabilizing material, where the active material includes a semiconductor, the stabilizing material includes at least one of an oligomer, an elastomer, a polymer, and/or a resin, and the stabilizing material provides to the device an improved performance metric compared to a device constructed of the first layer but constructed of only the active material (i.e., in the absence of the stabilizing material).

Zhu, Kai↗
DOE OSTI2020

sta-imex (scripts to investigate stability of IMEX methods)

Set of Matlab scripts to investigate properties of IMEX time stepping methods. Scripts include discretization of a system of linearized equations. One script plots numerical dispersion against analytical dispersion of the system. Another set of scripts uses the system to break it into stiff and non stiff part to investigate stability of Implicit-Explicit Runge-Kutta methods. Scripts are based on Thuburn, J. and Woollings, T. J., Vertical Discretizations for Compressible Euler Equation Atmospheric Models Giving Optimal Representation of Normal Modes , 2005 and Lock, S.-J. and Wood, N. and Weller, H., Numerical analyses of Runge-Kutta implicit-explicit schemes for horizontally explicit, vertically implicit solutions of atmospheric models, 2014 expanding ideas from these two papers. Scripts are used to develop IMEX schemes. SAND2020-12594 M Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525.

Guba, Oksana↗
DOE OSTI2020

Affinity propagation clustering of full-field, high-spatial-dimensional measurements for robust output-only modal identification: A proof-of-concept study

Determination of the model order is a challenging problem in system identification, especially in output-only or operational modal identification where some modes are weakly excited. While existing methods such as the stabilization diagram method (spectral information) are effective, they do not scale to high-dimensional data, which is usually needed for high-fidelity characterization of structural dynamics and has been made available in the emerging full-field measurement techniques using optical methods such as photogrammetry and laser vibrometers. In this proof-of-concept study we present a new non-parametric, data-driven approach for robust output-only identification of high-spatial-dimensional modal parameters of basic structures by efficiently processing and interactively exploiting the full-field measurement (i.e., very dense spatial measurement points). Specifically, we first over-estimate the system model once, producing a pool of candidate modes associated with their modal frequencies and full-field, high-spatial-dimensional mode shapes. This is accomplished by a data-driven method termed affinity propagation clustering (APC), where the active clusters, which are the active modes in our formulations, emerge from the "message-passing" procedure and does not require a pre-determination of the cluster number (mode or model order). Next, rather than using the spectral information to distinguish the physical and spurious modes in the stabilization diagram method, we exploit and visualize the spatial, full-field mode shape associated with each candidate mode to do so. We conduct extensive experiments on basic structural models with comparisons to a few existing methods. The results indicate that the new method is computationally efficient for identifying high-spatial-dimensional modal parameters, and robust to identify weak modes by exploiting the full-field measurement. We also discuss its applicability and limitations for structures with complex geometry (shapes).

42 ENGINEERING↗
DOE OSTI2022

Embedded pairs for optimal explicit strong stability preserving Runge–Kutta methods

We construct a family of embedded pairs for optimal explicit strong stability preserving Runge–Kutta methods of order 2 ≤ p ≤ 4 to be used to obtain numerical solution of spatially discretized hyperbolic PDEs. In this construction, the goals include non-defective property, large stability region, and small error values as defined in Dekker and Verwer (1984) and Kennedy et al. (2000). The new family of embedded pairs offer the ability for strong stability preserving (SSP) methods to adapt by varying the step-size. Through several numerical experiments, we assess the overall effectiveness in terms of work versus precision while also taking into consideration accuracy and stability.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Adaptive Model Reduction for Parareal in Time Method for Transient Stability Simulations

Real time or faster than real time simulation can enable system operators to foresee the effect of crucial contingencies on the power system dynamics and take timely actions to prevent system instability. Parareal in time method uses concurrent computations on different segments of the time domain of interest to speed up the dynamic simulations. This paper describes the application of an adaptive nonlinear model reduction method in improving computational speed of the Parareal solver. The proposed method adaptively switches between a hybrid system with selective linearization and a completely linear system based on the size of a disturbance. The functions in the hybrid system are linearized based on the electrical distance between specific generators and the area where disturbances originated. The proposed method is tested on the 327-machine 2383-bus Polish system.

Osipov, Denis↗
DOE OSTI2024

One-step sputtering of MoSSe metastable phase as thin film and predicted thermodynamic stability by computational methods

Abstract We present the fabrication of a MoS 2−x Se x thin film from a co-sputtering process using MoS 2 and MoSe 2 commercial targets with 99.9% purity. The sputtering of the MoS 2 and MoSe 2 was carried out using a straight and low-cost magnetron radio frequency sputtering recipe to achieve a MoS 2−x Se x phase with x = 1 and sharp interface formation as confirmed by Raman spectroscopy, time-of-flight secondary ion mass spectroscopy, and cross-sectional scanning electron microscopy. The sulfur and selenium atoms prefer to distribute randomly at the octahedral geometry of molybdenum inside the MoS 2−x Se x thin film, indicated by a blue shift in the A 1g and E 1 g vibrational modes at 355 cm −1 and 255 cm −1 , respectively. This work is complemented by computing the thermodynamic stability of a MoS 2−x Se x phase whereby density functional theory up to a maximum selenium concentration of 33.33 at.% in both a Janus-like and random distribution. Although the Janus-like and the random structures are in the same metastable state, the Janus-like structure is hindered by an energy barrier below selenium concentrations of 8 at.%. This research highlights the potential of transition metal dichalcogenides in mixed phases and the need for further exploration employing low-energy, large-scale methods to improve the materials’ fabrication and target latent applications of such structures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Impedance Methods for Analyzing Stability Impacts of Inverter-Based Resources: Stability Analysis Tools for Modern Power Systems

Power systems around the world are undergoing a major transformation because of the increasing shares of renewable energy, growing deployment of energy storage systems, proliferation of distributed energy resources, electrification of other sectors, and so on. In the United States, wind and solar provided almost 10% of electricity in 2019. The U.S. Energy Information Administration, in its 2020 Annual Energy Outlook, forecasted that the share of electricity from renewables will reach 38% by 2050, of which more than 80% will come from wind and solar. Wind and solar, along with battery energy storage systems, interface with the grid using power electronic inverters; hence, they are collectively referred to as inverter-based resources (IBRs). The increasing annual share of electricity from IBRs in a power system means that during more times of the year, the system will operate at a much higher concentration of IBRs. Figure 1 presents the hourly share of wind and solar generation in the Electric Reliability Council of Texas (ERCOT) system in Texas in 2019. While the annual wind share was at 20%, the instantaneous percentage share was much higher. Moments of high shares of IBRs (>50%) will continue to grow as more IBRs will be deployed in a power system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2019

Systems and methods for improved stability of power systems

The embodiments described herein provide for a system including a processor. The processor is configured to select at least one grid system contingency from a plurality of grid system contingencies. The processor is further configured to derive one or more eigen-sensitivity values based on the at least on grid system contingency. The processor is also configured to derive one or more control actions at least partially based on the eigen-sensitivity values. The processor is additionally configured to apply the one or more control actions for generation re-dispatch of a grid system.

Baone, Chaitanya Ashok↗
DOE OSTI2025

A review of high order strong stability preserving two-derivative explicit, implicit, and IMEX methods

High order strong stability preserving time discretizations ensure the nonlinear non-inner-product strong stability properties of spatial discretizations suited for the stable simulation of hyperbolic PDEs in a wide variety of application areas including fluid dynamics, magnetohydrodynamics, semiconductor devices, electromagnetics, and astrophysics. Over the past decade multiderivative time-stepping have been increasingly used for the time-evolution hyperbolic PDEs, so that the strong stability properties of these methods have become important. In this work we review sufficient conditions for a two-derivative multistage method to preserve the strong stability properties of spatial discretizations in a forward Euler and different conditions on the second derivative. In particular we present the strong stability preserving theory for explicit and implicit two-derivative Runge–Kutta schemes, including a special condition on the second derivative under which these implicit methods may be unconditionally strong stability preserving. This special condition is natural for the stiff component of wide range of plasma physics problems, and can be useful in the context of strong stability preserving implicit-explicit multi-derivative Runge–Kutta schemes, where the time-step restriction is then independent of the stiff term. Lastly, we present the strong stability preserving theory for implicit-explicit multi-derivative general linear methods, and some novel second and third order methods where the time-step restriction is independent of the stiff term.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Method of thermally-stabilizing an oxygen transport membrane-based reforming system

A method of operating an oxygen transport membrane based reforming system employing one or more packs of thermally coupled panels of reformer tubes and oxygen transport membrane ("OTM") reactors close to thermo-neutral point is provided. The method produces syngas by converting a hydrocarbon-containing feed, such as natural gas in the reformer tubes of a pack by endothermic steam reforming reactions. The heat required for endothermic reforming reactions is provided by exothermic oxidizing reactions occurring inside the OTM reactors of the pack. At a thermo-neutral point the heat released by exothermic reactions matches the heat required to support endothermic reactions and heat losses in the pack. The method modulates the flow rate of hydrocarbon-containing feed and/or steam-to-carbon ratio of the combined feed to the pack to maintain the surface temperature of oxygen transport membrane reactors below a target maximum temperature. The syngas product oxygen-to-carbon ratio is maintained within a desired target range such that the OTM based reforming system is operated close to thermo-neutral point.

42 ENGINEERING↗
DOE OSTI2024

Numerical integration in the virtual element method with the scaled boundary cubature scheme

Abstract The virtual element method (VEM) is a stabilized Galerkin method on meshes that consist of arbitrary (convex and nonconvex) polygonal and polyhedral elements. A crucial ingredient in the implementation of low‐ and high‐order VEM is the numerical integration of monomials and nonpolynomial functions over such elements. In this article, we apply the recently proposed scaled boundary cubature (SBC) scheme to compute the weak form integrals in various virtual element formulations over polygonal and polyhedral meshes. In doing so, we demonstrate the flexibility of the approach and the accuracy that it delivers on a broad suite of boundary‐value problems in 2D and 3D over polytopes with affine faces as well as on elements with curved boundaries. In addition, the use of the SBC scheme is exemplified in an enriched Poisson formulation of the VEM in which weakly singular functions are required to be integrated. This study establishes the SBC method as a simple, accurate and efficient integration scheme for use in the VEM.

Chin, Eric B.↗
DOE OSTI2020

Analysis of the SBP-SAT Stabilization for Finite Element Methods Part I: Linear Problems

In the hyperbolic community, discontinuous Galerkin (DG) approaches are mainly applied when finite element methods are considered. As the name suggested, the DG framework allows a discontinuity at the element interfaces, which seems for many researchers a favorable property in case of hyperbolic balance laws. On the contrary, continuous Galerkin methods appear to be unsuitable for hyperbolic problems and there exists still the perception that continuous Galerkin methods are notoriously unstable. To remedy this issue, stabilization terms are usually added and various formulations can be found in the literature. However, this perception is not true and the stabilization terms are unnecessary, in general. In this paper, we deal with this problem, but present a different approach. We use the boundary conditions to stabilize the scheme following a procedure that are frequently used in the finite difference community. Here, the main idea is to impose the boundary conditions weakly and specific boundary operators are constructed such that they guarantee stability. This approach has already been used in the discontinuous Galerkin framework, but here we apply it with a continuous Galerkin scheme. No internal dissipation is needed even if unstructured grids are used. Further, we point out that we do not need exact integration, it suffices if the quadrature rule and the norm in the differential operator are the same, such that the summation-by-parts property is fulfilled meaning that a discrete Gauss Theorem is valid. This contradicts the perception in the hyperbolic community that stability issues for pure Galerkin scheme exist. In numerical simulations, we verify our theoretical analysis.

97 MATHEMATICS AND COMPUTING↗