Search NASA⌕ Search

SEARCH · Search NASA

Results for “Asymptotic stability”

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 289 records · Page 16

Co-Spin With Symmetry Axis Stabilization, and De-Spin for Asteroid Capture

Consideration is given to attitude control associated with capturing a free-flying asteroid using an axisymmetric spacecraft. Asymptotically stable controllers are designed to align the spacecraft axis of symmetry with a line of descent that is fixed in the asteroid, and to eliminate all relative angular velocity before capture takes place. An analytical expression is presented for the torque required to maintain alignment of the axes of symmetry of the spacecraft and an axisymmetric asteroid. After the asteroid is securely captured, the angular velocity of the rigid composite body relative to an inertial frame is arrested; we present a controller that is asymptotically stable and stays within specified thrust limits.

Shen, Haijun↗

Omega-Limit Sets and Input-To-State Stability in Power Grids with Switching Equilibria

This paper studies a power transmission system with both conventional generators (CGs) and distributed energy assets (DEAs) providing frequency control. We consider an operating condition with demand aggregating two dynamic components: one that switches between different values on a finite set, and one that varies smoothly over time. Such dynamic operating conditions may result from protection scheme activations, external cyber-attacks, or due to the integration of dynamic loads, such as data centers. Mathematically, the dynamics of the resulting system are captured by a system that switches between a finite number of vector fields - or modes -, with each mode having a distinct equilibrium point induced by the demand aggregation. To analyze the stability properties of the resulting switching system, we leverage tools from hybrid dynamic inclusions and the concept of ..omega.. -limit sets from sets. Specifically, we characterize a compact set that is semi-globally practically asymptotically stable under the assumption that the switching frequency and load variation rate are sufficiently slow. For arbitrarily fast variations of the load, we use a level-set argument with multiple Lyapunov functions to establish input-to-state stability of a larger set and with respect to the rate of change of the loads. The theoretical results are illustrated via numerical simulations on the IEEE 39-bus test system.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Recursive identification and tracking of parameters for linear and nonlinear multivariable systems

The problem of identifying constant and variable parameters in multi-input, multi-output, linear and nonlinear systems is considered, using the maximum likelihood approach. An iterative algorithm, leading to recursive identification and tracking of the unknown parameters and the noise covariance matrix, is developed. Agile tracking, and accurate and unbiased identified parameters are obtained. Necessary conditions for a globally, asymptotically stable identification process are provided; the conditions proved to be useful and efficient. Among different cases studied, the stability derivatives of an aircraft were identified and some of the results are shown as examples.

Sidar, M.↗

Recursive identification and tracking of parameters for linear and non-linear multivariable systems

The problem of identifying constant and variable parameters in multi-input, multi-output, linear and nonlinear systems is considered, using the maximum likelihood approach. An iterative algorithm, leading to recursive identification and tracking of the unknown parameters and the noise covariance matrix, is developed. Agile tracking and accurate and unbiased identified parameters are obtained. Necessary conditions for a globally asymptotically stable identification process are provided; the conditions proved to be useful and efficient. Among different cases studied, the stability derivatives of an aircraft were identified and some of the results are shown as examples.

Sidar, M.↗

Numerical simulation of receptivity and transition in a boundary layer on a flat plate with a suction hole

A Navier-Stokes solver has been modified to allow for unsteady boundary conditions to be imposed at the inflow boundary. Because the streamwise derivatives are discretized with three-point compact differences, there is no need for special treatment of the points adjacent to the inflow boundary, provided that all necessary derivatives are properly specified at the boundary points. The unsteady boundary conditions can be obtained from a previous computation, from a solution of the parabolized stability equations (PSE), or from an analytical exact or approximate solution of the Navier-Stokes equations, e.g. an asymptotic expansion solution to model a sound wave interacting with the boundary layer.

Fasel, Hermann F.↗

Interface Morphology During Crystal Growth: Effects of Anisotropy and Fluid Flow

The effect of a parallel shear flow and anisotropic interface kinetics on the onset of instability during growth from a supersaturated solution is analyzed. The model used for anisotropy is based on the microscopic picture of step motion. A shear flow (linear Couette flow or asymptotic suction profile) parallel to the crystal solution interface in the same direction as the step motion decreases interface stability. A shear flow counter to the step motion enhances stability and for sufficiently large shear rates the interface is absolutely morphologically stable. For large wave numbers, the perturbed flow field can be neglected and a simple analytic approximation for the stability-instability demarcation is found.

Coriell, S. R.↗

An Eigenvalue Analysis of finite-difference approximations for hyperbolic IBVPs

The eigenvalue spectrum associated with a linear finite-difference approximation plays a crucial role in the stability analysis and in the actual computational performance of the discrete approximation. The eigenvalue spectrum associated with the Lax-Wendroff scheme applied to a model hyperbolic equation was investigated. For an initial-boundary-value problem (IBVP) on a finite domain, the eigenvalue or normal mode analysis is analytically intractable. A study of auxiliary problems (Dirichlet and quarter-plane) leads to asymptotic estimates of the eigenvalue spectrum and to an identification of individual modes as either benign or unstable. The asymptotic analysis establishes an intuitive as well as quantitative connection between the algebraic tests in the theory of Gustafsson, Kreiss, and Sundstrom and Lax-Richtmyer L(sub 2) stability on a finite domain.

Warming, Robert F.↗

An eigenvalue analysis of finite-difference approximations for hyperbolic IBVPs

The eigenvalue spectrum associated with a linear finite-difference approximation plays a crucial role in the stability analysis and in the actual computational performance of the discrete approximation. The eigenvalue spectrum associated with the Lax-Wendroff scheme applied to a model hyperbolic equation was investigated. For an initial-boundary-value problem (IBVP) on a finite domain, the eigenvalue or normal mode analysis is analytically intractable. A study of auxiliary problems (Dirichlet and quarter-plane) leads to asymptotic estimates of the eigenvalue spectrum and to an identification of individual modes as either benign or unstable. The asymptotic analysis establishes an intuitive as well as quantitative connection between the algebraic tests in the theory of Gustafsson, Kreiss, and Sundstrom and Lax-Richtmyer L (sub 2) stability on a finite domain.

Warming, Robert F.↗

A precision test of averaging in AdS/CFT

We reconsider the role of wormholes in the AdS/CFT correspondence. We focus on Euclidean wormholes that connect two asymptotically AdS or hyperbolic regions with $\mathbb{S}$ 1 × $\mathbb{S}$ d-1 boundary. There is no solution to Einstein’s equations of this sort, as the wormholes possess a modulus that runs to infinity. To find on-shell wormholes we must stabilize this modulus, which we can do by fixing the total energy on the two boundaries. Such a wormhole gives the saddle point approximation to a non-standard problem in quantum gravity, where we fix two asymptotic boundaries and constrain the common energy. Crucially the dual quantity does not factorize even when the bulk is dual to a single CFT, on account of the fixed energy constraint. From this quantity we extract a smeared version of the microcanonical spectral form factor. For a chaotic theory this quantity is self-averaging, i.e. well-approximated by averaging over energy windows, or over coupling constants. We go on to give a precision test involving the microcanonical spectral form factor where the two replicas have slightly different coupling constants. In chaotic theories this form factor is known to smoothly decay at a rate universally predicted in terms of one replica physics, provided that there is an average either over a window or over couplings. We compute the expected decay rate for holographic theories, and the form factor from a wormhole, and the two exactly agree for a wide range of two-derivative effective field theories in AdS. This gives a precision test of averaging in AdS/CFT. Our results interpret a number of confusing facts about wormholes and factorization in AdS and suggest that we should regard gravitational effective field theory as a mesoscopic description, analogous to semiclassical mesoscopic descriptions of quantum chaotic systems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Stability of an oscillated fluid with a uniform density gradient

Instabilities in a fluid with a constant density gradient that is subject to arbitrarily oriented oscillatory accelerations are considered. With the Boussinesq approximation and for the case of an unbounded fluid, transformation to Lagrangian coordinates allows the reduction of the problem to an ordinary differential equation for each three-dimensional wavenumber. The problem has three parameters: the nondimensional amplitude R of the base-state oscillation, the nondimensional level of background steady acceleration, which for some cases can be represented in terms of a local (in time) Richardson number Ri, and the Prandtl number Pr. Some general bounds on stability are derived. For Pr = 1 closed-form solutions are found for impulse (delta function) accelerations and a general asymptotic solution is constructed for large R and general imposed accelerations. The asymptotic solution takes advantage of the fact that at large R wave growth is concentrated at 'zero points'. These are times when the effective vertical wavenumber passes through zero. Kelvin-Helmholtz instabilities are found to dominate at low R, while Rayleigh-Taylor instabilities dominate at high R. At high R, the uniform shear of the Kelvin-Helmholtz case tends to distort and weaken instability waves. With unsteady flows, Ri = 1/4 is no longer an instability limit. Significant instabilities have been found for sinusoidal forcing for Ri up to 0.6.

Jacqmin, David↗

Convection in the Physical Vapor Transport Process: Thermosolutal - Part 2

We consider the effect of an inert gas on the diffusive-convective physical vapor transport process. We investigate the case when the temperature gradient is stabilizing and the concentration gradient is destabilizing for a wide parametric range. When an inert gas is present, the thermal and solutal convection oppose each other. The solutal field is destabilizing while the thermal field and the advective-diffusive flux stabilize the flow field. When the pressure of the inert component is increased, the stabilizing effect of the advective-diffusive flux is decreased. The intensity of convection as well as the oscillatory transient time increases. Below, the critical Rayleigh number, the nonlinear dynamics of the flow field show an oscillatory approach to steady state. For parametric values in the neighborhood of the critical Rayleigh number, the flow field undergoes a chaotic transient which settles to a periodic state. The asymptotic state of the flow field shows that growth and amalgamation of cells yields an overturning motion which results in an asymmetric cellular structure. The low gravity environment yields the stabilizing advective-diffusive flow which results in uniform temperature and concentration gradients near the crystal interface.

Duval, Walter M. B.↗

On regulators with a prescribed degree of stability

Several important aspects of the Regulator with a Prescribed Degree of Stability (RPDS) methodology and its applications are considered. The solution of the time varying RPDS problem as well as the characterization of RPDS closed loop eigenstructure properties are obtained. Based on the asymptotic behavior of RPDS root loci, a one step algorithm for designing Regulators with Prescribed Damping Ratio (RPDR) is developed. The robustness properties of RPDS are characterized in terms of the properties of the return difference and the inverse return difference matrices for the RPDS state feedback loop. This class of regulators is found to possess excellent multiloop margins with respect to stability and degree of stability properties. The ability of RPDS design to tolerate changing operating conditions and unmodelled dynamics are illustrated with a multiterminal dc/ac power system example. The output feedback realization of RPDS requires the use of Linear Quadratic Gaussian (LQG) methodology.

Ng, P. T. P.↗

On the nonlinear interfacial instability of rotating core-annular flow

The interfacial stability of rotating core-annular flows is investigated. The linear and nonlinear effects are considered for the case when the annular region is very thin. Both asymptotic and numerical methods are used to solve the flow in the core and film regions which are coupled by a difference in viscosity and density. The long-term behavior of the fluid-fluid interface is determined by deriving its nonlinear evolution in the form of a modified Kuramoto-Sivashinsky equation. We obtain a generalization of this equation to three dimensions. The flows considered are applicable to a wide array of physical problems where liquid films are used to lubricate higher or lower viscosity core fluids, for which a concentric arrangement is desired. Linearized solutions show that the effects of density and viscosity stratification are crucial to the stability of the interface. Rotation generally destabilizes non-axisymmetric disturbances to the interface, whereas the centripetal forces tend to stabilize flows in which the film contains the heavier fluid. Nonlinear affects allow finite amplitude helically travelling waves to exist when the fluids have different viscosities.

Coward, Aidrian V.↗

Performance of a low-power subsonic-arc-attachment arcjet thruster

A subsonic-arc-attachment thruster design was scaled from a 30 kW 1960's vintage thruster to operate at nominally 3 kW. Performance measurements were obtained over a 1-4 kW power range using hydrogen as the propellant. Several modes of operation were identified and were characterized by varying degrees of voltage instability. A stability map was developed showing that the voltage oscillations were brought upon by elevated current or propellant levels. At a given specific energy level the specific impulse increased asymptotically with increased flow rates. Comparisons of performance were made between radial and tangential propellant injection. When the vortex flow was eliminated using radial injection, the operating voltages were lower at a given current, and the specific impulse and efficiency decreased. Tests were also conducted to determine the effects of background pressure on operation. Performance data were obtained at pressures of 0.047 Pa and 18 Pa. For a given specific energy level, the performance increased with a decrease in facility background pressure. Lowering the background pressure also caused a dramatic change in the voltage-current characteristic and the voltage stability, a phenomenon not previously reported with conventional supersonic-arc-attachment thrusters.

Sankovic, John M.↗

Performance of a low-power subsonic-arc-attachment arcjet thruster

A subsonic-arc-attachment thruster design was scaled from a 30 kW 1960's vintage thruster to operate at nominally 3 kW. Performance measurements were obtained over a 1-4 kW power range using hydrogen as the propellant. Several modes of operation were identified and were characterized by varying degrees of voltage instability. A stability map was developed showing that the voltage oscillations were brought upon by elevated current or propellant levels. At a given specific energy level the specific impulse increased asymptotically with increased flow rates. Comparisons of performance were made between radial and tangential propellant injection. When the vortex flow was eliminated using radial injection, the operating voltages were lower at a given current, and the specific impulse and efficiency decreased. Tests were also conducted to determine the effects of background pressure on operation, and performance data were obtained at pressures of 0.047 Pa and 18 Pa. For a given specific energy level, the performance increased with a decrease in facility background pressure. Lowering the background pressure also caused a dramatic change in the voltage-current characteristic and the voltage stability, a phenomenon not previously reported with conventional supersonic-arc-attachment thrusters.

Sankovic, John M.↗

Asymptotic vacuum solution at tokamak X-point tip

In the H-mode regime of diverted tokamaks, the presence of strong pressure gradients in the pedestal gives rise to a sizable bootstrap current, together with the Ohmic and Pfirsch–Schlueter currents, close to the separatrix. For such equilibria, the presence of finite current density close to the separatrix requires the reexamination of equilibrium properties. It is almost universally assumed that the two branches of the separatrix (the stable and unstable manifolds) are straight as they cross at the X-point. However, the opposite angles of the plasma-filled segment and vacuum one cannot be equal if the current density does not vanish at the separatrix on the plasma side. We solve this difficulty by chipping off a thin layer of plasma edge so that the sharp corner of the plasma-filled segment becomes a hyperbola. Using the conformal transformation, we found that in the assumption of a hyperbolic boundary, the X point moves beyond the plasma boundary to fall in the vacuum region. An acute angle of the plasma-filled segment leads to an obtuse opposite angle of vacuum segment and vice versa. In the case of an acute angle of the plasma-filled segment, the new X point shifts inside the X point formed by the asymptotes of a hyperbolic boundary; in the case of an obtuse angle of the plasma-filled segment, the new X point shifts outside the X point formed by the asymptotes of a hyperbolic plasma boundary. Furthermore, the results are important for understanding the X point features, which affect the tokamak edge stability and transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonlinear stability of Taylor's vortex array

It is proved that the two-dimensional Taylor vortex array, which is an exact unsteady solution of the Navier-Stokes equation, is globally and asymptotically stable in the mean with respect to three-dimensional periodic disturbances. A time-dependent bound on the decay rate of the kinetic energy of disturbances is obtained.

Lin, S. P.↗

A performance comparison of integration algorithms in simulating flexible structures

Asymptotic formulas for the characteristic root errors as well as transfer function gain and phase errors are presented for a number of traditional and new integration methods. Normalized stability regions in the lambda h plane are compared for the various methods. In particular, it is shown that a modified form of Euler integration with root matching is an especially efficient method for simulating lightly-damped structural modes. The method has been used successfully for structural bending modes in the real-time simulation of missiles. Performance of this algorithm is compared with other special algorithms, including the state-transition method. A predictor-corrector version of the modified Euler algorithm permits it to be extended to the simulation of nonlinear models of the type likely to be obtained when using the discretized structure approach. Performance of the different integration methods is also compared for integration step sizes larger than those for which the asymptotic formulas are valid. It is concluded that many traditional integration methods, such as RD-4, are not competitive in the simulation of lightly damped structures.

Howe, R. M.↗