Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear stability analysis”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

A multiloop generalization of the circle criterion for stability margin analysis

In order to provide a theoretical tool suited for characterizing the stability margins of multiloop feedback systems, multiloop input-output stability results generalizing the circle stability criterion are considered. Generalized conic sectors with 'centers' and 'radii' determined by linear dynamical operators are employed to specify the stability margins as a frequency dependent convex set of modeling errors (including nonlinearities, gain variations and phase variations) which the system must be able to tolerate in each feedback loop without instability. The resulting stability criterion gives sufficient conditions for closed loop stability in the presence of frequency dependent modeling errors, even when the modeling errors occur simultaneously in all loops. The stability conditions yield an easily interpreted scalar measure of the amount by which a multiloop system exceeds, or falls short of, its stability margin specifications.

Safonov, M. G.↗
NASA NTRS1995

The initial-value problem for viscous channel flows

Plane viscous channel flows are perturbed and the ensuing initial-value problems are investigated in detail. Unlike traditional methods where traveling wave normal modes are assumed for solution, this works offers a means whereby completely arbitrary initial input can be specified without having to resort to eigenfunction expansions. The full temporal behavior, including both early time transients and the long time asymptotics, can be determined for any initial disturbance. Effects of three-dimensionality can be assessed. The bases for the analysis are: (a) linearization of the governing equations; (b) Fourier decomposition in the spanwise and streamwise directions of the flow; and (c) direct numerical integration of the resulting partial differential equations. All of the stability data that are known for such flows can be reproduced. Also, the optimal initial condition can be determined in a straight forward manner and such optimal conditions clearly reflect transient growth data that is easily determined by a rational choice of a basis for the initial conditions. Although there can be significant transient growth for subcritical values of the Reynolds number using this approach it does not appear possible that arbitrary initial conditions will lead to the exceptionally large transient amplitudes that have been determined by optimization of normal modes. The approach is general and can be applied to other classes of problems where only a finite discrete spectrum exists, such as the boundary layer for example.

Criminale, W. O.↗
NASA NTRS2010

Development and Breakdown of Goertler Vortices in High Speed Boundary Layers

The nonlinear development of G rtler instability over a concave surface gives rise to a highly distorted stationary flow in the boundary layer that has strong velocity gradients in both spanwise and wall-normal directions. This distorted flow is susceptible to strong, high frequency secondary instability that leads to the onset of transition. For high Mach number flows, the boundary layer is also subject to the second mode instability. The nonlinear development of G rtler vortices and the ensuing growth and breakdown of secondary instability, the G rtler vortex interactions with second mode instabilities as well as oblique second mode interactions are examined in the context of both internal and external hypersonic configurations using nonlinear parabolized stability equations, 2-D eigenvalue analysis and direct numerical simulation. For G rtler vortex development inside the Purdue Mach 6 Ludwieg tube wind tunnel, multiple families of unstable secondary eigenmodes are identified and their linear and nonlinear evolution is examined. The computation of secondary instability is continued past the onset of transition to elucidate the physical mechanisms underlying the laminar breakdown process. Nonlinear breakdown scenarios associated with transition over a Mach 6 compression cone configuration are also explored.

Li, Fei↗
NASA NTRSDate not supplied

A Gaussian Process Enhancement to Linear Parameter Varying Models

Simulation and analysis for modern engineering systems now routinely requires the merging of multiple disciplines, physical-domains, time-scales, and data sets — all at ever increasing levels. These capabilities are especially needed in the domain of Advanced Air Mobility, where rapidly emerging vehicle designs are significantly more complex, while having to be both cost-effective and safe. To meet these engineering challenges, machine learning methods are an attractive option for merging models and data across multiple areas while providing uncertainty quantification and maintaining computational efficiency. This paper examines the use of Gaussian process machine learning to generalize and enhance the commonly used class of quasi-Linear Parameter Varying models for fast full-envelope simulation while also supporting control system design and analysis with model uncertainty. Gaussian process machine learning is selected because it: can fuse multiple data sets, enables an easy trade-off between data fitting and smoothing, provides model uncertainty quantification, scales well with increasing complexity, and does not generally require starting from a large training data set. To demonstrate the benefits of the approach, a robust stability analysis with Gaussian process uncertainty is shown for a NASA reference design of an electric quad-rotor air-taxi concept vehicle with motor parameter uncertainty.

Gaussian Process↗
NASA NTRS1992

Nonlinear stability and control study of highly maneuverable high performance aircraft, phase 2

This research should lead to the development of new nonlinear methodologies for the adaptive control and stability analysis of high angle-of-attack aircraft such as the F18 (HARV). The emphasis has been on nonlinear adaptive control, but associated model development, system identification, stability analysis and simulation is performed in some detail as well. Various models under investigation for different purposes are summarized in tabular form. Models and simulation for the longitudinal dynamics have been developed for all types except the nonlinear ordinary differential equation model. Briefly, studies completed indicate that nonlinear adaptive control can outperform linear adaptive control for rapid maneuvers with large changes in alpha. The transient responses are compared where the desired alpha varies from 5 degrees to 60 degrees to 30 degrees and back to 5 degrees in all about 16 sec. Here, the horizontal stabilator is the only control used with an assumed first-order linear actuator with a 1/30 sec time constant.

Mohler, R. R.↗
DOE OSTI2023

Oxide derived Cu nanofibril assembly for enhanced nonenzymatic glucose sensing

In this report, we have prepared copper nanofibril assembly via thermal oxidation followed by electrochemical reduction processes, exhibiting superior glucose detection ability. The morphological analysis evidenced the formation of rough islands or nanofibril structure on the Cu surface depending on initial thermal oxidation temperature. The glucose detection performance was investigated by performing cyclic voltammetry and chronoamperometry with varying glucose concentration. A high sensitivity (4131.57 μA mM – 1 cm – 2 ), low detection limit (1.41 μM), wider linear range (0–3.9 mM) and long–term stability (30 days) have been recorded for electrode thermally oxidized at 400 °C followed by electrochemical reduction process (rCu_400). The sensitivity is almost three times higher in comparison to the planner Cu surface. This significantly enhanced glucose sensing ability rCu_400 has been attributed to the nanofibril morphology, origination of Cu (111) facet and formation of a stable oxide layer evidenced by scanning electron microscopy, X-ray diffraction and energy dispersive spectroscopy analysis. Despite higher sensitivity, rCu_400 electrode does not show any response to the chloride ion, dopamine, ascorbic acid and uric acid. Finally, these results indicate that the Cu nanofibril structure prepared via simple oxidation/reduction process can be an excellent candidate to be used as an electrode for glucose sensing application.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Stability of hairy black holes in shift-symmetric scalar-tensor theories via the effective field theory approach

Shift-symmetric Horndeski theories admit an interesting class of Schwarzschild-de Sitter black hole solutions exhibiting time-dependent scalar hair. The properties of these solutions may be studied via a bottom-up effective field theory (EFT) based on the background symmetries. This is in part possible by making use of a convenient coordinate choice — Lemaître-type coordinates — in which the profile of the Horndeski scalar field is linear in the relevant time coordinate. We construct this EFT, and use it to understand the stability of hairy black holes in shift-symmetric Horndeski theories, providing a set of constraints that the otherwise-free functions appearing in the Horndeski Lagrangian must satisfy in order to admit stable black hole solutions. The EFT is analyzed in the decoupling limit to understand potential sources of instability. Further, we also perform a complete analysis of the EFT with odd-parity linear perturbations around general spherically symmetric space-time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A General Autofluorescence Method to Characterize Polymerization Progress

An autofluorescence technique to characterize polymerization progress in real time/in line was developed, which functioned in the absence of typical fluorogenic groups on the monomer or polymer. The monomer dicyclopentadiene and polymer polydicyclopentadiene are hydrocarbons that lack traditional functional groups for fluorescence spectroscopy. Here, the autofluorescence of formulations containing this monomer and polymer during ruthenium-catalyzed ring-opening metathesis polymerization (ROMP) was harnessed for reaction monitoring. The methods fluorescence recovery after photobleaching (FRAP) and here-developed fluorescence lifetime recovery after photobleaching (FLRAP) characterized polymerization progress in these native systems—without requiring exogenous fluorophore. (Auto)fluorescence lifetime recovery changes during polymerization correlated linearly to degree of cure, providing a quantitative link with reaction progress. These changing signals also provided relative rates of background polymerization, enabling comparison of 10 different catalyst–inhibitor-stabilized formulations. Multiple-well analysis demonstrated suitability for future high-throughput evaluation of formulations for thermosets. Furthermore, the central concept of the combined autofluorescence and FLRAP/FRAP method may be extendable to monitoring other polymerization reactions previously overlooked for lack of an obvious fluorescence handle.

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

W7-X and the sawtooth instability: towards realistic simulations of current-driven magnetic reconnection

Magnetic reconnection in W7-X is studied by means of global numerical simulations in a series of models of increasing complexity. The magnetic geometry ranges from that of a cylinder to the full three-dimensional field of W7-X, and the equations solved range from ideal magnetohydrodynamics (MHD) to gyrokinetics. We simulate plasmas from the first operation phase with electron cyclotron current drive (ECCD). These are characterized by an equilibrium magnetic field featuring an ECCD-distorted 'humped' profile of the rotational transform ι, with ι = 1 in two radial locations. Such plasmas generally show sawtooth activity, hence motivating the present study. We pay particular attention to the role of equilibrium current density gradients in the destabilization of reconnecting modes. When the equilibrium temperature and density gradients are artificially suppressed (to eliminate the pressure gradient drive), the perturbed electrostatic potential is radially localized between the locations at which ι = 1. This is shown with a purely collisionless gyrokinetic model, in cylindrical geometry. In the real toroidal geometry of W7-X, for a non-ideal MHD model including a uniform resistivity, electron inertia and (numerical) viscosity, the same qualitative behaviour is observed. In particular, even if a resonant (m, n) = (1, -1) perturbation is initialized, the most unstable mode is the (m, n) = (-4, 4), where m and n are the poloidal and toroidal mode numbers, respectively. Other modes are destabilized due to geometric coupling. The growth rate of this instability scales as η 1/3 , where η is the plasma resistivity, thus suggesting that ECCD drives ideal MHD stable W7-X plasmas towards non-ideal marginality. An ideal magnetohydrodynamic analysis confirms the result. A fluid-kinetic hybrid version of the EUTERPE code shows that gyrokinetic ions have a stabilizing effect on these modes. For W7-X relevant collisionalities, the growth rate scales linearly with the electron skin depth, d e . Implications of our results for sawtoothing W7-X operation are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2014

Stability Assessment and Tuning of an Adaptively Augmented Classical Controller for Launch Vehicle Flight Control

Recently, a robust and practical adaptive control scheme for launch vehicles [ [1] has been introduced. It augments a classical controller with a real-time loop-gain adaptation, and it is therefore called Adaptive Augmentation Control (AAC). The loop-gain will be increased from the nominal design when the tracking error between the (filtered) output and the (filtered) command trajectory is large; whereas it will be decreased when excitation of flex or sloshing modes are detected. There is a need to determine the range and rate of the loop-gain adaptation in order to retain (exponential) stability, which is critical in vehicle operation, and to develop some theoretically based heuristic tuning methods for the adaptive law gain parameters. The classical launch vehicle flight controller design technics are based on gain-scheduling, whereby the launch vehicle dynamics model is linearized at selected operating points along the nominal tracking command trajectory, and Linear Time-Invariant (LTI) controller design techniques are employed to ensure asymptotic stability of the tracking error dynamics, typically by meeting some prescribed Gain Margin (GM) and Phase Margin (PM) specifications. The controller gains at the design points are then scheduled, tuned and sometimes interpolated to achieve good performance and stability robustness under external disturbances (e.g. winds) and structural perturbations (e.g. vehicle modeling errors). While the GM does give a bound for loop-gain variation without losing stability, it is for constant dispersions of the loop-gain because the GM is based on frequency-domain analysis, which is applicable only for LTI systems. The real-time adaptive loop-gain variation of the AAC effectively renders the closed-loop system a time-varying system, for which it is well-known that the LTI system stability criterion is neither necessary nor sufficient when applying to a Linear Time-Varying (LTV) system in a frozen-time fashion. Therefore, a generalized stability metric for time-varying loop=gain perturbations is needed for the AAC.

VanZwieten, Tannen↗
DOE OSTI2018

Theoretical analysis and numerical simulation for a hyperbolic equation with Dirichlet and acoustic boundary conditions

This paper is concerned with a theoretical and numerical study for the initial-boundary value problem for a linear hyperbolic equation with variable coefficient and acoustic boundary conditions. On the theoretical results, we prove the existence and uniqueness of global solutions, and the uniform stability of the total energy. Numerical simulations using the finite element method associated with the finite difference method are employed, for one-dimensional and two-dimensional cases, to validate the theoretical results. In addition, numerically the uniform decay rate for energy and the order of convergence of the approximate solution are also shown.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

The Integrated Reference Region Analysis for Parallel DFIGs’ Interfacing Inductors

Although the traditional design of doubly-fed induction generators (DFIGs)’ interfacing inductors consider the peak ripple of the grid-side converter (GSC)’s output current, it does not consider the inductors’ impact on DFIGs’ smallsignal stability. Located in series between the GSC and the stator, the inappropriate selection of the interfacing inductor can easily result to system instability. Therefore, this paper first proposes the integrated reference region analysis for small wind farm's interfacing inductors to improve the traditional design method. Firstly, a new detailed parallel DFIGs’ smallsignal model that considers output currents’ coupling is built in d-q coordinate system with state-space approach. The model focuses on representing the operating states of parallel DFIGs in wind farm. Secondly, the linearized state-space matrix of parallel DFIGs is decomposed into nominal-value matrix and location matrix. Considering the traditional design requirements, the integrated reference region for interfacing inductor is proposed through spectral radius and bialternate matrix sum (BMS). Furthermore, it can provide better guidance for parameter selecting and stabilization method researches. Finally, the simulation and experimental results show that the proposed integrated reference region for parallel DFIGs’ interfacing inductors is accurate and instructive.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Plasticizer-Free Gradient-Crosslinked Polyurethane Electrolyte for Room-Temperature Solid-State Lithium Batteries

Polymer-based solid-state electrolytes are promising for next-generation lithium metal batteries, yet their limited ionic conductivity and mechanical stability at ambient conditions remain key challenges. Herein, we report a novel gradient crosslinked polymer electrolyte (PU/PUA/PU) synthesized via a sequential in situ UV-curing process that integrates a mechanically robust polyurethane acrylate (PUA) core with soft linear polyurethane (PU) interfaces. This all-solid membrane operates without any liquid plasticizer and the interfacial PU layers ensure low interfacial resistance and intimate electrode contact, while the PUA core provides enhanced dimensional stability and dendrite suppression. As a result, the gradient electrolyte delivers an impressive ∼2.6 × 10 −4 S cm −1 ionic conductivity at 25 °C (two orders of magnitude higher than conventional PEO) and remains electrochemically stable >5 V (vs Li + /Li). Structural analysis confirms the formation of a well-defined crosslinked network with suppressed crystallinity and expanded interchain spacing. Electrochemical impedance spectroscopy (EIS) and linear sweep voltammetry (LSV) further validate these properties. When applied in a Li||NMC811 cell, the PU/PUA/PU electrolyte delivers stable cycling with high capacity at both 60 °C and 25 °C, demonstrating its potential for room-temperature solid-state battery applications without reliance on plasticizers or heating. In conclusion, this gradient design offers a practical path toward ambient-condition operation of solid-state lithium batteries, providing a paradigm for overcoming the traditional trade-off between ionic conductivity and mechanical robustness.

25 ENERGY STORAGE↗
NASA NTRS1994

Hypersonic vehicle model and control law development using H(infinity) and micron synthesis

The control system design for a Single Stage To Orbit (SSTO) air breathing vehicle will be central to a successful mission because a precise ascent trajectory will preserve narrow payload margins. The air breathing propulsion system requires the vehicle to fly roughly halfway around the Earth through atmospheric turbulence. The turbulence, the high sensitivity of the propulsion system to inlet flow conditions, the relatively large uncertainty of the parameters characterizing the vehicle, and continuous acceleration make the problem especially challenging. Adequate stability margins must be provided without sacrificing payload mass since payload margins are critical. Therefore, a multivariable control theory capable of explicitly including both uncertainty and performance is needed. The H(infinity) controller in general provides good robustness but can result in conservative solutions for practical problems involving structured uncertainty. Structured singular value mu framework for analysis and synthesis is potentially much less conservative and hence more appropriate for problems with tight margins. An SSTO control system requires: highly accurate tracking of velocity and altitude commands while limiting angle-of-attack oscillations, minimized control power usage, and a stabilized vehicle when atmospheric turbulence and system uncertainty are present. The controller designs using H(infinity) and mu-synthesis procedures were compared. An integrated flight/propulsion dynamic mathematical model of a conical accelerator vehicle was linearized as the vehicle accelerated through Mach 8. Vehicle acceleration through the selected flight condition gives rise to parametric variation that was modeled as a structured uncertainty. The mu-analysis approach was used in the frequency domain to conduct controller analysis and was confirmed by time history plots. Results demonstrate the inherent advantages of the mu framework for this class of problems.

Gregory, Irene M.↗
DOE OSTI2020

Distributed non-negative matrix factorization with determination of the number of latent features

The holistic analysis and understanding of the latent (that is, not directly observable) variables and patterns buried in large datasets is crucial for data-driven science, decision making and emergency response. Such exploratory analyses require devising unsupervised learning methods for data mining and extraction of the latent features, and non-negative matrix factorization (NMF) is one of the prominent such methods. NMF is based on compute-intense non-convex constrained minimization, which, for large datasets requires fast and distributed algorithms. However, current parallel implementations of NMF fail to estimate the number of latent features. In practice, identifying these features is both difficult and significant for pattern recognition and latent feature analysis, especially for large dense matrices. Here, we introduce a distributed NMF algorithm coupled with distributed custom clustering followed by a stability analysis on dense data, which we call DnMFk, to determine the number of latent variables. The results on synthetic data and the classical Swimmer data set demonstrate the accuracy of model determination while scaling nearly linearly across multiple processors for large data. Further, we employ DnMFk to determine the number of hidden features from a terabyte matrix.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1973

Thermomechanical behavior of a polynorbornadiene

A high-temperature linear amorphous hydrocarbon polymer synthesized cationically from 2,5-norbornadiene was studied thermomechanically with respect to physical transitions and stability in nitrogen. The glass transition temperature was determined to be 320 C (at less than 1 cps), which is probably the highest known for a linear hydrocarbon addition polymer. The thermomechanical technique of torsional braid analysis, together with thermogravimetric analysis, differential thermal analysis, infrared analysis, and solubility studies, was used to investigate the sequential events of the glass transition and degradation. The polymer is of particular interest since it is a high-temperature plastic which in the bulk form would probably need to be processed at high speeds in the vicinity of the glass transition temperature in an inert atmosphere. The presence of tertiary hydrogen atoms should render it amenable to degradation by the earth's environment.

Roller, M. B.↗
NASA NTRS1989

On the role of ozone in the stability of Rossby normal modes

The role of ozone in the linear stability of Rossby normal modes is examined in a continuously stratified, extratropical baroclinic atmosphere. The flow is described by coupled equations for the quasi-geostraphic potential vorticity and ozone volume mixing ratio. A perturbation analysis is carried out under the assumption of weak diabatic heating, which is generated by Newtonian cooling and dynamics-ozone interaction. An expression for the propagation and growth characteristics is obtained analytically in terms of the vertically averaged wave activity, which depends explicitly on the wave spatial structure, photochemistry, and basic state distributions of wind, temperature, and ozone mixing ratio. Calculations show that stationary internal modes, whose amplitudes are largest in the stratosphre, are destabilized by dynamics-ozone interaction and Newtonian cooling, with e-folding times on the order of 20-40 days.

Nathan, Terrence R.↗
NASA NTRS1997

Protein Crystal Movements and Fluid Flows During Microgravity Growth

The growth of protein crystals suitable for X-ray crystal structure analysis is an important topic. The methods of protein crystal growth are under increasing study whereby different methods are being compared via diagnostic monitoring including Charge Coupled Device (CCD) video and interferometry. The quality (perfection) of protein crystals is now being evaluated by mosaicity analysis (rocking curves) and X-ray topographic images as well as the diffraction resolution limit and overall data quality. Choice of a liquid-liquid linear crystal growth geometry and microgravity can yield a spatial stability of growing crystals and fluid, as seen in protein crystallization experiments on the unmanned platform EURICA. A review is given here of existing results and experience over several microgravity missions. The results include CCD video as well as interferometry during the mission, followed, on return to earth, by rocking curve experiments and full X-ray data collection on LMS and earth control lysozyme crystals. Diffraction data recorded from LMS and ground control apocrustacyanin C(sub 1) crystals are also described.

Boggon, Titus J.↗