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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

An integrating matrix formulation for buckling of rotating beams including the effects of concentrated masses

The integrating matrix technique of computational mechanics is extended to include the effects of concentrated masses. The stability of a flexible rotating beam with discrete masses is analyzed to determine the critical rotational speeds for buckling in the inplane and out-of-plane directions. In this problem, the beam is subjected to compressive centrifugal forces arising from steady rotation about an axis which does not pass through the clamped end of the beam. To determine the eigenvalues from which stability is assessed, the differential equations of motion are solved numerically by combining the extended integrating matrix method with an eigenanalysis. Stability boundaries for a discrete mass representation of a uniform beam are shown to asymptotically approach the stability boundaries for the corresponding continuous mass beam as the number of concentrated masses is increased. An error in the literature is also noted for the discrete mass problem concerning the behavior of the critical rotational speed for inplane buckling as the radius of rotation of the clamped end of the beam is reduced.

Lakin, William D.↗
DOE OSTI2022

High Order Strong Stability Preserving MultiDerivative Implicit and IMEX Runge--Kutta Methods with Asymptotic Preserving Properties

In this article we present a class of high order unconditionally strong stability preserving (SSP) implicit two-derivative Runge--Kutta schemes and SSP implicit-explicit (IMEX) multi-derivative Runge--Kutta schemes where the time-step restriction is independent of the stiff term. The unconditional SSP property for a method of order $p>2$ is unique among SSP methods and depends on a backward-in-time assumption on the derivative of the operator. We show that this backward derivative condition is satisfied in many relevant cases where SSP IMEX schemes are desired. We devise unconditionally SSP implicit Runge--Kutta schemes of order up to $p=4$ and IMEX Runge--Kutta schemes of order up to $p=3$. For the multiderivative IMEX schemes, we also derive and present the order conditions, which have not appeared previously. The unconditional SSP condition ensures that these methods are positivity preserving, and we present sufficient conditions under which such methods are also asymptotic preserving when applied to a range of problems, including a hyperbolic relaxation system, the Broadwell model, and the Bhatnagar--Gross--Krook kinetic equation. We present numerical results to support the theoretical results on a variety of problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

The DESC stellarator code suite. Part 1. Quick and accurate equilibria computations

Three-dimensional equilibrium codes are vital for stellarator design and operation, and high-accuracy equilibria are also necessary for stability studies. This paper details comparisons of two three-dimensional equilibrium codes: VMEC, which uses a steepest-descent algorithm to reach a minimum-energy plasma state, and DESC, which minimizes the magnetohydrodynamic (MHD) force error in real space directly. Accuracy as measured by satisfaction of MHD force balance is presented for each code, along with the computation time. It is shown that DESC is able to achieve more accurate solutions, especially near axis. The importance of higher-accuracy equilibria is shown in DESC's better agreement of stability metrics with asymptotic formulae. DESC's global Fourier–Zernike basis also yields solutions with analytic derivatives explicitly everywhere in the plasma volume, provides improved accuracy in the radial direction versus conventional finite differences and allows for exponential convergence. Further, DESC can compute a solution with the same accuracy as VMEC in order-of-magnitude less time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1973

Numerical experiments in climate stability.

Two semiempirical climate models are formulated as time-dependent problems in order to study the stability of their asymptotic steady state equilibrium solutions to perturbations in internal (in this case, initial) conditions. For fixed external conditions the time-dependent versions of these models exhibit fully transitive behavior to positive perturbations in initial conditions, slight intransivity to negative perturbations up to -18 K, and an ice-covered earth regime is obtained for extremely large negative perturbations in initial temperatures (below -18 K). The parametrization found to be most critical in these models is the albedo-temperature coupling, especially in tropical regions. The temperature albedo parametrization from an intransitive climate model used in a time-dependent form of one of the preceding models is found to be highly sensitive to negative perturbations in both internal and external conditions. Numerical experimentation with these semiempirical models shows the important role of the tropics in maintaining the equilibrium climate and suggests that the radiation balance in equatorial latitudes might have a controlling influence on the equatorward extent of the polar ice cap.

Schneider, S. H.↗
NASA NTRS1980

Strongly turbulent stabilization of electron beam-plasma interactions

The stabilization of electron beam interactions due to strongly turbulent nonlinearities is studied analytically and numerically for a wide range of plasma parameters. A fluid mode coupling code is described in which the effects of electron and ion Landau damping and linear growth due to the energetic electron beam are included in a phenomenological manner. Stabilization of the instability is found to occur when the amplitudes of the unstable modes exceed the threshold of the oscillating two-stream instability. The coordinate space structure of the turbulent spectrum which results clearly shows that soliton-like structures are formed by this process. Phenomenological models of both the initial stabilization and the asymptotic states are developed. Scaling laws between the beam-plasma growth rate and the fluctuations in the fields and plasma density are found in both cases, and shown to be in good agreement with the results of the simulation.

Freund, H. P.↗
NASA NTRS1985

Vector Lyapunov Functions for Stochastic Interconnected Systems

Theoretical paper presents set of sufficient conditions for asymptotic and exponential stability with probability 1 for class of stochastic interconnected systems. Theory applicable to complicated, large-scale mechanical or electrical systems, and, for several design problems, it reduces computational difficulty by relating stability criteria to fundamental structural features of system.

Boussalis, D.↗
NASA NTRS1984

Random equations in aerodynamics

Literature was reviewed to identify aerodynamic models which might be treated by probablistic methods. The numerical solution of some integral equations that arise in aerodynamical problems were investigated. On the basis of the numerical studies a qualitative theory of random integral equations was developed to provide information on the behavior of the solutions of these equations (in particular, boundary and asymptotic behavior, and stability) and their statistical properties without actually obtaining explicit solutions of the equations.

Bharucha-Reid, A. T.↗
NASA NTRS1999

Directional Solidification of a Binary Alloy into a Cellular Convective Flow: Localized Morphologies

A steady, two dimensional cellular convection modifies the morphological instability of a binary alloy that undergoes directional solidification. When the convection wavelength is far longer than that of the morphological cells, the behavior of the moving front is described by a slow, spatial-temporal dynamics obtained through a multiple-scale analysis. The resulting system has a "parametric-excitation" structure in space, with complex parameters characterizing the interactions between flow, solute diffusion, and rejection. The convection stabilizes two dimensional disturbances oriented with the flow, but destabilizes three dimensional disturbances in general. When the flow is weak, the morphological instability behaves incommensurably to the flow wavelength, but becomes quantized and forced to fit into the flow-box as the flow gets stronger. At large flow magnitudes the instability is localized, confined in narrow envelopes with cells traveling with the flow. In this case the solutions are discrete eigenstates in an unbounded space. Their stability boundary and asymptotics are obtained by the WKB analysis.

Chen, Y.- J.↗
DOE OSTI2021

Stability Analysis of Parameter Varying Genetic Toggle Switches Using Koopman Operators

The genetic toggle switch is a well known model in synthetic biology that represents the dynamic interactions between two genes that repress each other. The mathematical models for the genetic toggle switch that currently exist have been useful in describing circuit dynamics in rapidly dividing cells, assuming fixed or time-invariant kinetic rates. There is a growing interest in being able to model and extend synthetic biological function for growth conditions such as stationary phase or during nutrient starvation. As cells transition from one growth phase to another, kinetic rates become time-varying parameters. In this paper, we propose a novel class of parameter varying nonlinear models that can be used to describe the dynamics of genetic circuits, including the toggle switch, as they transition from different phases of growth. We show that there exists unique solutions for this class of systems, as well as for a class of systems that incorporates the microbial phenomena of quorum sensing. Further, we show that the domain of these systems, which is the positive orthant, is positively invariant. We also showcase a theoretical control strategy for these systems that would grant asymptotic monostability of a desired fixed point. We then take the general form of these systems and analyze their stability properties through the framework of time-varying Koopman operator theory. A necessary condition for asymptotic stability is also provided as well as a sufficient condition for instability. A Koopman control strategy for the system is also proposed, as well as an analogous discrete time-varying Koopman framework for applications with regularly sampled measurements.

Harrison, Jamiree (ORCID:0000000265351579)↗
NASA NTRS1996

Mixing and Transition Control Studied

Considerable progress in understanding nonlinear phenomena in both unbounded and wallbounded shear flow transition has been made through the use of a combination of high- Reynolds-number asymptotic and numerical methods. The objective of this continuing work is to fully understand the nonlinear dynamics so that ultimately (1) an effective means of mixing and transition control can be developed and (2) the source terms in the aeroacoustic noise problem can be modeled more accurately. Two important aspects of the work are that (1) the disturbances evolve from strictly linear instability waves on weakly nonparallel mean flows so that the proper upstream conditions are applied in the nonlinear or wave-interaction streamwise region and (2) the asymptotic formulations lead to parabolic problems so that the question of proper out-flow boundary conditions--still a research issue for direct numerical simulations of convectively unstable shear flows--does not arise. Composite expansion techniques are used to obtain solutions that account for both mean-flow-evolution and nonlinear effects. A previously derived theory for the amplitude evolution of a two-dimensional instability wave in an incompressible mixing layer (which is in quantitative agreement with available experimental data for the first nonlinear saturation stage for a plane-jet shear layer, a circular-jet shear layer, and a mixing layer behind a splitter plate) have been extended to include a wave-interaction stage with a three-dimensional subharmonic. The ultimate wave interaction effects can either give rise to explosive growth or an equilibrium solution, both of which are intimately associated with the nonlinear self-interaction of the three dimensional component. The extended theory is being evaluated numerically. In contrast to the mixing-layer situation, earlier comparisons of theoretical predictions based on asymptotic methods and experiments in wall-bounded shear-flow transition have been somewhat lacking in one aspect or another. The current work strongly suggests that the main weakness is the underlying asymptotic representation of the linear "part" of the problem and not the explicit modeling of the nonlinear/wave-interaction effects. Consequently, the long-wave-length/high-Reynolds-number asymptotic limit for the Blasius boundary-layer stability problem was reexamined, and a new dispersion relationship for the instability waves that is uniformly valid for both the upper- and lower branch regions to the required order of approximation was obtained. A comparison with numerical results, obtained by solving the Orr-Sommerfeld stability problem, shows that the asymptotic formula provides surprisingly good results, even for values of the frequency parameter usually encountered in experimental investigations. This is particularly evident in the dynamically important upper-branch region, where much of the nonlinear interactions in transition experiments are believed to take place. The result is important in that it can be used to greatly improve the accuracy of weakly nonlinear critical-layer-based theories, and a consistent nonlinear theory is currently under evaluation.

Source record↗
NASA NTRS1990

The inviscid stability of supersonic flow past heated or cooled axisymmetric bodies

The inviscid, linear, nonaxisymmetric, temporal stability of the boundary layer associated with the supersonic flow past axisymmetric bodies (with particular emphasis on long thin, straight circular cylinders), subject to heated or cooled wall conditions is investigated. The eigenvalue problem is computed in some detail for a particular Mach number or 3.8, revealing that the effect of curvature and the choice of wall conditions both have a significant effect on the stability of the flow. Both the asymptotic, large azimuthal wavenumber solution and the asymptotic, far downstream solution are obtained for the stability analysis and compared with numerical results. Additionally, asymptotic analyses valid for large radii of curvature with cooled/heated wall conditions, are presented. In general, important differences were found to exist between the wall temperature conditions imposed and the adiabatic wall conditions considered previously.

Shaw, Stephen J.↗
NASA NTRS1992

The inviscid stability of supersonic flow past heated or cooled axisymmetric bodies

The inviscid, linear, nonaxisymmetric, temporal stability of the boundary layer associated with the supersonic flow past axisymmetric bodies (with particular emphasis on long thin, straight circular cylinders), subject to heated or cooled wall conditions is investigated. The eigenvalue problem is computed in some detail for a particular Mach number or 3.8, revealing that the effect of curvature and the choice of wall conditions both have a significant effect on the stability of the flow. Both the asymptotic, large azimuthal wavenumber solution and the asymptotic, far downstream solution are obtained for the stability analysis and compared with numerical results. Additionally, asymptotic analyses valid for large radii of curvature with cooled/heated wall conditions are presented. In general, important differences were found to exist between the wall temperature conditions imposed and the adiabatic wall conditions considered previously.

Shaw, Stephen J.↗
NASA NTRS1971

Stability of uncertain systems

The asymptotic properties of feedback systems are discussed, containing uncertain parameters and subjected to stochastic perturbations. The approach is functional analytic in flavor and thereby avoids the use of Markov techniques and auxiliary Lyapunov functionals characteristic of the existing work in this area. The results are given for the probability distributions of the accessible signals in the system and are proved using the Prohorov theory of the convergence of measures. For general nonlinear systems, a result similar to the small loop-gain theorem of deterministic stability theory is given. Boundedness is a property of the induced distributions of the signals and not the usual notion of boundedness in norm. For the special class of feedback systems formed by the cascade of a white noise, a sector nonlinearity and convolution operator conditions are given to insure the total boundedness of the overall feedback system.

Blankenship, G. L.↗
NASA NTRS1994

The stability of two-phase flow over a swept-wing

We use numerical and asymptotic techniques to study the stability of a two-phase air/water flow above a flat porous plate. This flow is a model of the boundary layer which forms on a yawed cylinder and can be used as a useful approximation to the air flow over swept wings during heavy rainfall. We show that the interface between the water and air layers can significantly destabilize the flow, leading to traveling wave disturbances which move along the attachment line. This instability occurs for lower Reynolds numbers than in the case of the absence of a water layer. We also investigate the instability of inviscid stationary modes. We calculate the effective wavenumber and orientation of the stationary disturbance when the fluids have identical physical properties. Using perturbation methods we obtain corrections due to a small stratification in viscosity, thus quantifying the interfacial effects. Our analytical results are in agreement with the numerical solution which we obtain for arbitrary fluid properties.

Coward, Adrian↗
NASA NTRS2010

Asymptotic Linearity of Optimal Control Modification Adaptive Law with Analytical Stability Margins

Optimal control modification has been developed to improve robustness to model-reference adaptive control. For systems with linear matched uncertainty, optimal control modification adaptive law can be shown by a singular perturbation argument to possess an outer solution that exhibits a linear asymptotic property. Analytical expressions of phase and time delay margins for the outer solution can be obtained. Using the gradient projection operator, a free design parameter of the adaptive law can be selected to satisfy stability margins.

Nguyen, Nhan T.↗
NASA NTRS1995

Dynamical Approach Study of Spurious Steady-State Numerical Solutions of Nonlinear Differential Equations: Global Asymptotic Behavior of Time Discretizations - 2

The global asymptotic nonlinear behavior of 1 1 explicit and implicit time discretizations for four 2 x 2 systems of first-order autonomous nonlinear ordinary differential equations (ODES) is analyzed. The objectives are to gain a basic understanding of the difference in the dynamics of numerics between the scalars and systems of nonlinear autonomous ODEs and to set a baseline global asymptotic solution behavior of these schemes for practical computations in computational fluid dynamics. We show how 'numerical' basins of attraction can complement the bifurcation diagrams in gaining more detailed global asymptotic behavior of time discretizations for nonlinear differential equations (DEs). We show how in the presence of spurious asymptotes the basins of the true stable steady states can be segmented by the basins of the spurious stable and unstable asymptotes. One major consequence of this phenomenon which is not commonly known is that this spurious behavior can result in a dramatic distortion and, in most cases, a dramatic shrinkage and segmentation of the basin of attraction of the true solution for finite time steps. Such distortion, shrinkage and segmentation of the numerical basins of attraction will occur regardless of the stability of the spurious asymptotes, and will occur for unconditionally stable implicit linear multistep methods. In other words, for the same (common) steady-state solution the associated basin of attraction of the DE might be very different from the discretized counterparts and the numerical basin of attraction can be very different from numerical method to numerical method. The results can be used as an explanation for possible causes of error, and slow convergence and nonconvergence of steady-state numerical solutions when using the time-dependent approach for nonlinear hyperbolic or parabolic PDES.

Yee, H. C.↗