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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Implicit Multirate GARK Methods

This article considers multirate generalized-structure additively partitioned Runge–Kutta methods for solving stiff systems of ordinary differential equations with multiple time scales. These methods treat different partitions of the system with different timesteps for a more targeted and efficient solution compared to monolithic single rate approaches. With implicit methods used across all partitions, methods must find a balance between stability and the cost of solving nonlinear equations for the stages. In order to characterize this important trade-off, we explore multirate coupling strategies, problems for assessing linear stability, and techniques to efficiently implement Newton iterations for stage equations. Unlike much of the existing multirate stability analysis which is limited in scope to particular methods, we present general statements on stability and describe fundamental limitations for certain types of multirate schemes. New implicit multirate methods up to fourth order are derived, and their accuracy and efficiency properties are verified with numerical tests.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1958

Estimation of Directional Stability Derivatives at Small Angles and Subsonic and Supersonic Speeds

Methods are presented for estimating the directional stability derivative increments contributed by the stabilizing surfaces of subsonic and supersonic aircraft. These methods are strictly applicable at zero angle of attack and small angles of sideslip. The procedure of totaling the incremental coefficients to obtain an estimation of the total empennage side-force and yawing-moment coefficient derivatives is also shown, together with numerical examples. A correlation is presented between estimated and experimental incremental coefficients which indicates that the methods of this report generally estimate the increment of side force gained by the addition of a panel to within +/-10 percent of the experimental value while the yawing-moment increment is generally estimated to within +/-20 percent. This is true for both subsonic and supersonic Mach numbers. An example application of the methods to one of the problems in directional stability, that of minimizing the effect of Mach number on the side-force coefficient derivative of the empennage, is discussed.

Goodwin, Frederick K.↗
NASA NTRS1989

A new uniformly valid asymptotic integration algorithm for elasto-plastic-creep and unified viscoplastic theories including continuum damage

A new scheme to integrate a system of stiff differential equations for both the elasto-plastic creep and the unified viscoplastic theories is presented. The method has high stability, allows large time increments, and is implicit and iterative. It is suitable for use with continuum damage theories. The scheme was incorporated into MARC, a commercial finite element code through a user subroutine called HYPELA. Results from numerical problems under complex loading histories are presented for both small and large scale analysis. To demonstrate the scheme's accuracy and efficiency, comparisons to a self-adaptive forward Euler method are made.

Chulya, A.↗
NASA NTRS1991

A new uniformly valid asymptotic integration algorithm for elasto-plastic creep and unified viscoplastic theories including continuum damage

A new scheme to integrate a system of stiff differential equations for both the elasto-plastic creep and the unified viscoplastic theories is presented. The method has high stability, allows large time increments, and is implicit and iterative. It is suitable for use with continuum damage theories. The scheme was incorporated into MARC, a commercial finite element code through a user subroutine called HYPELA. Results from numerical problems under complex loading histories are presented for both small and large scale analysis. To demonstrate the scheme's accuracy and efficiency, comparisons to a self-adaptive forward Euler method are made.

Chulya, Abhisak↗
DOE OSTI2024

Cover crops and poultry litter impact on soil structural stability in dryland soybean production in southeastern United States

Abstract This study explored the efficacy of soil aggregate indices in quantifying soil structural development, utilizing 5‐year field experiment data from the Southeastern United States. The experiment utilized a split‐plot design with cover crops (native vegetation as control, cereal rye (Secale cerealeL.), winter wheat (Triticum aestivum), hairy vetch (Vicia villosa), and mustard (Brassica rapa) plus cereal rye as the main factor and fertilizer source (no fertilizer as control, inorganic fertilizer with phosphorus, potassium, and elemental sulfur, and poultry litter) as the secondary factor. Aggregate size fractions were determined using the wet‐sieving method, and aggregate stability index (ASI), mean weight diameter (MWD), geometric mean diameter (GMD), and fractal dimension (FD) were calculated to assess soil structural stability. Main effects results indicated that cereal rye (55.11%) and poultry litter (50.97%) exhibited the highest ASI values. The highest MWD, GMD, and FD were observed under mustard plus cereal rye (1.187 mm), cereal rye (0.462 mm), and hairy vetch (2.573), respectively. Principal component analysis revealed that cover crops significantly improved soil aggregate structure and stability, overcoming limitations of sole fertilization practices. Regression analysis suggested that ASI, MWD, and GWD positively correlated with soil organic carbon, whereas FD negatively correlated with MWD, GMD, and ASI. Principal component analysis exhibited that FD decreased with increasing soil organic carbon, ASI, MWD, and GMD, demonstrating that lower FD values indicate enhanced soil aggregation and structure. Assessed indices, FD included, effectively gauged soil structural stability. These metrics should be prioritized in managerial decisions to support soil productivity and health in agricultural systems.

Agriculture↗
DOE OSTI2022

From Event Data to Wind Power Plant DQ Admittance and Stability Risk Assessment

This paper presents a dynamic event data-based stability risk assessment method for power grids with high penetrations of inverter-based resources (IBRs). This method relies on obtaining the IBRs' DQ admittance through dynamic event data and computing the system's eigenvalues based on the admittance models. Two critical technologies are employed in this research, including time-domain and frequency-domain data fitting and dq-frame voltage and current signal derivation. The first technology is key to obtaining the s-domain expressions from the transient response data, and the s-domain DQ admittance model from the frequency-domain measurements. The second technology is key to obtaining the dq-frame voltage and current signals from either the three-phase instantaneous measurements or the phasor measurement unit (PMU) data. The method is illustrated using data generated from a Type-4 wind power plant modeled in PSCAD. This paper demonstrates the technical feasibility of the proposed approach.

17 WIND ENERGY↗
NASA NTRS1950

Transonic Stability and Control Investigation of a 1/80-Scale Model of the Consolidated Vultee Skate 9 Seaplane, TED No. NACA DE 345: Transonic-Bump Method

An investigation of the longitudinal stability and of the all-movable horizontal tail and aileron control of a 1/80-scale reflection-plane model of the Consolidated Vultee Skate 9 seaplane has been made through a Mach number range of 0.6 to 1.16 on the transonic bump of the Langley high-speed 7- by 10-foot tunnel. At moderate lift coefficients (0.4 to 0.8) and below a Mach number of 1.0 the model was statically unstable longitudinally. The static longitudinal stability of the model at low lift coefficients increased with Mach number corresponding to a shift in aerodynamic center from 37 percent mean aerodynamic chord at a Mach number of 0.60 to 64 percent at a Mach number of 1.10. Estimates indicate that the tail deflection angle required for steady flight and for accelerated maneuvers of the Skate 9 airplane would probably not vary greatly with Mach number at sea level, but for accelerated maneuvers at altitude the tail deflection angle would probably vary erratically with Mach number. The variation of rolling-moment coefficient with aileron deflection angle was approximately linear, agreed well with theory, and held for the range of aileron deflections tested (-17.1 deg to 16.6 deg). At low lift coefficients the drag rise occurred at Mach numbers of 0.95 and 0.90 for the wing alone and the complete model, respectively. The effects of the canopy on the model were small. For the ranges investigated, angle-of-attack and Mach number changes caused no large pressure drops in the jet-engine duct.

Riebe, John M.↗
NASA NTRS2009

Bounded Linear Stability Analysis - A Time Delay Margin Estimation Approach for Adaptive Control

This paper presents a method for estimating time delay margin for model-reference adaptive control of systems with almost linear structured uncertainty. The bounded linear stability analysis method seeks to represent the conventional model-reference adaptive law by a locally bounded linear approximation within a small time window using the comparison lemma. The locally bounded linear approximation of the combined adaptive system is cast in a form of an input-time-delay differential equation over a small time window. The time delay margin of this system represents a local stability measure and is computed analytically by a matrix measure method, which provides a simple analytical technique for estimating an upper bound of time delay margin. Based on simulation results for a scalar model-reference adaptive control system, both the bounded linear stability method and the matrix measure method are seen to provide a reasonably accurate and yet not too conservative time delay margin estimation.

Nguyen, Nhan T.↗
DOE OSTI2019

Review of stability and thermal conductivity enhancements for salt hydrates

Salt hydrates can be used as phase change materials for thermal energy storage. Critical technical challenges for their widespread deployment include poor cycling stability, large supercooling, and low thermal conductivity. Here in this work, numerous enhancement techniques are reviewed. Related to stability, numerous existing methods to characterize the number of thermal cycles are summarized. Following this, 11 techniques to mitigate phase separation (plus 9 studies on macro-, micro-, and nano-encapsulation types) are reviewed. For supercooling, 38 nucleator-salt hydrate combinations to minimize subcooling are reviewed. The empirically observed effect of isomorphism in minimizing subcooling is explored and discussed, with cross-study trends depicted for 38 nucleators across 9 salt hydrates. In addition, several studies reporting combined effects on stability and supercooling are presented. Related to thermal conductivity, 32 combinations of conductivity enhancement material and salt hydrate are reviewed. For those using graphite, the dependence of conductivity enhancement on graphite mass fraction is shown across numerous studies for different graphite types. The performance and stability of calcium chloride hexahydrate, sodium sulfate decahydrate, and sodium acetate were explored and are discussed in-depth. Finally, the status of enhancement to salt hydrates is summarized, and remaining challenges are identified.

25 ENERGY STORAGE↗
NASA NTRS1975

Prediction of axial-flow instabilities in a turbojet engine by use of a multistage compressor simulation on the digital computer

A method of estimating the undistorted stall line for an axial-flow compressor by using the digital computer is presented. The method involves linearization of nonlinear dynamic equations about an operating point on a speed line, and then application of the first method of Lyapunov to determine the stability of the nonlinear system from the stability of the linear system. The method is applied to a simulation of the J85 compressor, which utilizes stage stacking and lumped volume techniques for the interstage regions to simulate steady-state and dynamic compressor performance. The stability boundary predicted by the digital simulation compares quite well with the stall line predicted by a dynamic simulation of the J85 compressor programmed on the analog computer. Since previous studies have shown that the analog-predicted stall line agrees well with the stall line of the compressor, the digital method presented is also a good means of estimating the stall line.

Daniele, C. J.↗
DOE OSTI2021

Tailoring capping-layer composition for improved stability of mixed-halide perovskites

Incorporating a low dimensional (LD) perovskite capping layer on top of a perovskite absorber, improves the stability of perovskite solar cells (PSCs). However, in the case of mixed-halide perovskites, which can undergo halide segregation into single-halide perovskites, a systematic study of the capping layer's effect on mixed-halide perovskite absorber is still lacking. This study bridges this gap by investigating how the 1D perovskite capping layers on top of MAPb(I x Br 1–x ) 3 (x = 0, 0.25, 0.5, 0.75, 1) absorbers affect the films' stability. We utilize a new method, dissimilarity matrix, to investigate the image-based stability performance of capping-absorber pair compositions across time. This method overcomes the challenge of analyzing various film colors due to bandgap difference in mixed-halide perovskites. We also discover that the intrinsic absorber stability plays an important role in the overall stability outcome, despite the capping layer's support. Within the 55 unique capping-absorber pairs, we observe a notable 1D perovskite material, 1-methoxynaphthalene-2-ethylammonium chloride (2MeO–NEA–Cl or 9-Cl), that improves the stability of MAPbI 3 and MAPb(I 0.5 Br 0.5 ) 3 by at least 8 and 1.5 times, respectively, compared to bare films under elevated humidity and temperature. Surface photovoltage results also show that the accumulation of electrostatic charges on the film surface depends on the capping layer type, which could contribute to the acceleration/deceleration of degradation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Multirate linearly-implicit GARK schemes

Many complex applications require the solution of initial-value problems where some components change fast, while others vary slowly. Multirate schemes apply different step sizes to resolve different components of the system, according to their dynamics, in order to achieve increased computational efficiency. The stiff components of the system, fast or slow, are best discretized with implicit base methods in order to ensure numerical stability. To this end, linearly implicit methods are particularly attractive as they solve only linear systems of equations at each step. This paper develops the Multirate GARK-ROS/ROW (MR-GARK-ROS/ROW) framework for linearly-implicit multirate time integration. The order conditions theory considers both exact and approximative Jacobians. The effectiveness of implicit multirate methods depends on the coupling between the slow and fast computations; an array of efficient coupling strategies and the resulting numerical schemes are analyzed. Multirate infinitesimal step linearly-implicit methods, that allow arbitrarily small micro-steps and offer extreme computational flexibility, are constructed. The new unifying framework includes existing multirate Rosenbrock(-W) methods as particular cases, and opens the possibility to develop new classes of highly effective linearly implicit multirate integrators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Wideband precision stabilization of the -18.6kV retarding voltage for the KATRIN spectrometer

The Karlsruhe Tritium Neutrino Experiment (KATRIN) measures the effective electron anti-neutrino mass with an unprecedented design sensitivity of 0.2 eV (90 % C.L.). In this experiment, the energy spectrum of beta electrons near the tritium decay endpoint is analyzed with a highly accurate spectrometer. To reach the KATRIN sensitivity target, the retarding voltage of this spectrometer must be stable to the ppm (1 × 10 -6 ) level and well known on various time scales (μs up to months), for values around -18.6 kV. A custom-designed high-voltage regulation system mitigates the impact of interference sources in the absence of a closed electric shield around the large spectrometer vessel. In this article, we describe the regulation system and its integration into the KATRIN setup. Independent monitoring methods demonstrate a stability within 2 ppm, exceeding KATRIN's specifications.

47 OTHER INSTRUMENTATION↗
NASA NTRS1989

The computation of thermo-chemical nonequilibrium hypersonic flows

Several conceptual designs for vehicles that would fly in the atmosphere at hypersonic speeds have been developed recently. For the proposed flight conditions the air in the shock layer that envelops the body is at a sufficiently high temperature to cause chemical reaction, vibrational excitation, and ionization. However, these processes occur at finite rates which, when coupled with large convection speeds, cause the gas to be removed from thermo-chemical equilibrium. This non-ideal behavior affects the aerothermal loading on the vehicle and has ramifications in its design. A numerical method to solve the equations that describe these types of flows in 2-D was developed. The state of the gas is represented with seven chemical species, a separate vibrational temperature for each diatomic species, an electron translational temperature, and a mass-average translational-rotational temperature for the heavy particles. The equations for this gas model are solved numerically in a fully coupled fashion using an implicit finite volume time-marching technique. Gauss-Seidel line-relaxation is used to reduce the cost of the solution and flux-dependent differencing is employed to maintain stability. The numerical method was tested against several experiments. The calculated bow shock wave detachment on a sphere and two cones was compared to those measured in ground testing facilities. The computed peak electron number density on a sphere-cone was compared to that measured in a flight test. In each case the results from the numerical method were in excellent agreement with experiment. The technique was used to predict the aerothermal loads on an Aeroassisted Orbital Transfer Vehicle including radiative heating. These results indicate that the current physical model of high temperature air is appropriate and that the numerical algorithm is capable of treating this class of flows.

Candler, Graham↗
DOE OSTI2024

Topological Interpretability for Deep Learning

With the growing adoption of AI-based systems across everyday life, the need to understand their decision-making mechanisms is correspondingly increasing. The level at which we can trust the statistical inferences made from AI-based decision systems is an increasing concern, especially in high-risk systems such as criminal justice or medical diagnosis, where incorrect inferences may have tragic consequences. Despite their successes in providing solutions to problems involving real-world data, deep learning (DL) models cannot quantify the certainty of their predictions. These models are frequently quite confident, even when their solutions are incorrect. This work presents a method to infer prominent features in two DL classification models trained on clinical and non-clinical text by employing techniques from topological and geometric data analysis. We create a graph of a model's feature space and cluster the inputs into the graph's vertices by the similarity of features and prediction statistics. We then extract subgraphs demonstrating high-predictive accuracy for a given label. These subgraphs contain a wealth of information about features that the DL model has recognized as relevant to its decisions. We infer these features for a given label using a distance metric between probability measures, and demonstrate the stability of our method compared to the LIME and SHAP interpretability methods. This work establishes that we may gain insights into the decision mechanism of a DL model. This method allows us to ascertain if the model is making its decisions based on information germane to the problem or identifies extraneous patterns within the data.

Spannaus, Adam↗
NASA NTRS1986

Numerical advection by conservation of second-order moments

A new, accurate, and nondiffusive method for three-dimensional advection of trace species is presented. The method preserves tracer structures by conserving the second-order moments of the spatial distribution of tracer during advection. Upstream transport and second-order tracer distribution are described, and the moments of the tracer distribution about the center of a grid box are formally defined and related to the polynomial distribution. Formulas are presented which describe how the moments of a grid box are decomposed into a unique set of moments centered about each subbox and how they are reassembled into new grid boxes. A one-dimensional example if tracer transport is given, and limits necessary to maintain positive tracer concentrations are derived. The accuracy and stability of this method are analytically examined, and numerical experiments testing the effective resolution are presented. The method is compared with other methods for numerical advection of tracers.

Prather, M. J.↗
NASA NTRS2018

Thermal Stability Analysis of Gevo Jet Fuel Using Ellipsometry

Thermal stability is an important characteristic of alternative fuels that must be evaluated before they can be used in aviation engines. This characteristic is of great importance to the effectiveness of the fuel as a coolant and to the engine's combustion performance. The thermal stability of Gevo fuel, an alcohol to jet fuel made from plant feedstock was studied. This analysis was used to make comments on the effectiveness of the current thermal stability test standard. This work was performed using a spectroscopic ellipsometer to measure the thickness of deposits left on aluminum substrates. It was observed that Gevo deposit thickness increased slowly up to 375 C and much more rapidly after that point. Similar behavior was observed in JP-8 fuel. Comparisons were also made between color standard ratings and ellipsometric thickness measurements, and it was found that in some cases, darker colors did not indicate thicker deposits. Reference tubes were used to validate the optical models used in this work, and different optical constants were found to best model the results than what are published in the ASTM D3241 test method for thermal stability.

Ellipsometry↗