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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

All sources NASA NTRS DOE OSTI

NASA NTRS1975

A theoretical evaluation of rigid baffles in suppression of combustion instability

An analytical technique for the prediction of the effects of rigid baffles on the stability of liquid propellant combustors is presented. This analysis employs both two and three dimensional combustor models characterized by concentrated combustion sources at the chamber injector and a constant Mach number nozzle. An eigenfunction-matching method is used to solve the linearized partial differential equations describing the unsteady flow field for both models. Boundary layer corrections to this unsteady flow are used in a mechanical energy dissipation model to evaluate viscous and turbulence effects within the flow. An integral stability relationship is then employed to predict the decay rate of the oscillations. Results of this analysis agree qualitatively with experimental observations and show that sufficient dissipation exists to indicate that the proper mechanism of baffle damping is a fluid dynamic loss.

Combustion instability↗

NASA NTRS2013

Estimation of Aerodynamic Stability Derivatives for Space Launch System and Impact on Stability Margins

This paper describes the techniques involved in determining the aerodynamic stability derivatives for the frequency domain analysis of the Space Launch System (SLS) vehicle. Generally for launch vehicles, determination of the derivatives is fairly straightforward since the aerodynamic data is usually linear through a moderate range of angle of attack. However, if the wind tunnel data lacks proper corrections then nonlinearities and asymmetric behavior may appear in the aerodynamic database coefficients. In this case, computing the derivatives becomes a non-trivial task. Errors in computing the nominal derivatives could lead to improper interpretation regarding the natural stability of the system and tuning of the controller parameters, which would impact both stability and performance. The aerodynamic derivatives are also provided at off nominal operating conditions used for dispersed frequency domain Monte Carlo analysis. Finally, results are shown to illustrate that the effects of aerodynamic cross axis coupling can be neglected for the SLS configuration studied

Pei, Jing↗

NASA NTRS1988

Stabilization and robustness of non-linear unity-feedback system - Factorization approach

The paper is a self-contained discussion of a right factorization approach in the stability analysis of the nonlinear continuous-time or discrete-time, time-invariant or time-varying, well-posed unity-feedback system S1(P, C). It is shown that a well-posed stable feedback system S1(P, C) implies that P and C have right factorizations. In the case where C is stable, P has a normalized right-coprime factorization. The factorization approach is used in stabilization and simultaneous stabilization results.

Desoer, C. A.↗

DOE OSTI2019

Modulational instability and soliton trains in a model for two-mode fiber ring lasers

A model of orthogonally polarized two-field fiber ring laser with a linear gain is considered, with emphasis on the continuous-wave stability and the existence of soliton trains. The continuous-wave stability analysis is carried out within the framework of the modulational-instability approach, the variations of the gain spectrum with the modulation frequency and characteristic parameters of the model give rise to a rich variety of stability features including single-band and multiband stability regions. Seeking for pulse structures of the model, the two coupled cubic complex Ginzburg–Landau equations describing individual mode propagations are transformed into a set of coupled, first-order nonlinear ordinary-differential equations for the amplitudes and phases of the two modes. Numerical simulations of the last set of coupled equations indicate that in the anomalous dispersion regime, envelopes of the two fields are periodic trains of pulses the amplitudes of which are affected by the linear gain.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Study of Alfvén eigenmode stability in Quasi-Poloidal Stellarator (QPS) plasma using a Landau closure model

The aim of this study is to analyze the linear stability of Alfvén eigenmodes (AE) in the QPS device heated by a tangential neutral beam injector (NBI). The analysis is performed using the gyro-fluid code FAR3d, that solves the reduced MHD equations for the thermal plasma coupled with moments of the kinetic equation for the energetic particles (EP). The AE stability is calculated in several operational regimes of the tangential NBI: EP β between 0.001 and 0.1, EP energy between 12 and 180 keV and different radial locations of the beam. The analysis is performed for vacuum and finite β equilibria as well as QPS configurations with two and three periods. The EP β threshold in the vacuum case is 0.001 and the AE frequency is lower as the energy of the EP population decreases. Toroidal Alfvén eigenmodes with f = 80–120 kHz and elliptical AE between f = 120–350 kHz are triggered between the middle-outer plasma region ($r/a \gt 0.5$). The AE stability improves in the simulations with finite β equilibria and three period configurations with respect to the vacuum case with two periods because the continuum gaps are slender, leading to a higher threshold of the EP β, above 0.03 for the AEs triggered by the helical mode families. Helical effects are not strong enough to destabilize Helical Alfvén eigenmodes, the AEs with the largest growth rates are triggered by the n = 1 and n = 2 toroidal families.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1971

Concerning the origin of novae and U Geminorum stars

The dynamical stability of the contact component in a semi-detached binary system is investigated by a linear time dependent analysis. The boundary condition imposed on the photosphere by the Roche lobe is assumed to be one of constant pressure over the whole surface of the lobe, which is treated as being spherically symmetric. The eigenvalues of the linear equation of motion are evaluated for sequences of stellar envelope models under the constraint of this boundary condition. An extensive region in the Hertzprung-Russell diagram is found to be dynamically unstable. This region coincides almost precisely with that predicted by a previous quasi-static study. Time scales for the growth of the instability are also evaluated. Color observations of several novae and U Geminorum stars are discussed in an attempt to provide semi-empirical tests of the proposed theoretical model of the outbursts. In particular the track in the H-R diagram of the U Geminorum star, EM Cygni, is shown to be consistent with an outburst which is purely a temperature change at constant radius, as is demanded by the theory.

Bath, G. T.↗

NASA NTRS1974

Stability region maximization by decomposition-aggregation method

This work is to improve the estimates of the stability regions by formulating and resolving a proper maximization problem. The solution of the problem provides the best estimate of the maximal value of the structural parameter and at the same time yields the optimum comparison system, which can be used to determine the degree of stability of the Skylab. The analysis procedure is completely computerized, resulting in a flexible and powerful tool for stability considerations of large-scale linear as well as nonlinear systems.

Siljak, D. D.↗

NASA NTRS1974

On the use of first order rotor dynamics in multiblade coordinates

This paper is directed to the question of how to represent most efficiently rotor/body coupling in a linear flight dynamics analysis. Rigid body pitch, roll and vertical motions are considered for the rotor/body coupling studies. Flapping stability limits, eigenvalues, transient responses to control step inputs, to step gusts and to random gusts are determined for a hypothetical hingeless compound helicopter operating up to .8 advance ratio. Data are obtained for the basic helicopter and for the craft with two simple control feedback systems. While complete periodic system modeling is necessary for determining flapping stability limits and vibrations, constant system modeling using first order dynamics in each of the multiblade rotor coordinates was found to be adequate for rotor-craft stability and response computations.

Hohenemser, K. H.↗

NASA NTRS1977

TRIM-STAB-Aerospace vehicle trim and stability

Program calculates stability and control characteristics of aerospace vehicles, which are usally calculated by either constant coefficient time-slice linear analysis or by nonlinear simulation.

Barnes, R. E.↗

NASA NTRS1999

N-S Simulations of Crow-Type Instabilities in Vortex Wake

Three-dimensional Navier-Stokes simulations of the Crow instability of wake vortices are conducted using large-eddy simulations without background turbulence. Sinusoidal displacement has been specified as the initial perturbation for the vortex system. The results have shown that the minimum Crow instability wavelength is about one vortex spacing shorter than predicted by Crow's linear stability theory. The planar- standing-wave-angle value and the amplitude indifference behavior agree with Crow's analysis. Simulations with periodicity in the axial direction have indicated minor influence of axial flow on the Crow instability.

Zheng, Z. C.↗

DOE OSTI2023

Mathematical solutions in internal dose assessment: A comparison of Python-based differential equation solvers in biokinetic modeling

Abstract In biokinetic modeling systems employed for radiation protection, biological retention and excretion have been modeled as a series of discretized compartments representing the organs and tissues of the human body. Fractional retention and excretion in these organ and tissue systems have been mathematically governed by a series of coupled first-order ordinary differential equations (ODEs). The coupled ODE systems comprising the biokinetic models are usually stiff due to the severe difference between rapid and slow transfers between compartments. In this study, the capabilities of solving a complex coupled system of ODEs for biokinetic modeling were evaluated by comparing different Python programming language solvers and solving methods with the motivation of establishing a framework that enables multi-level analysis. The stability of the solvers was analyzed to select the best performers for solving the biokinetic problems. A Python-based linear algebraic method was also explored to examine how the numerical methods deviated from an analytical or semi-analytical method. Results demonstrated that customized implicit methods resulted in an enhanced stable solution for the inhaled 60 Co (Type M) and 131 I (Type F) exposure scenarios for the inhalation pathway of the International Commission on Radiological Protection (ICRP) Publication 130 Human Respiratory Tract Model (HRTM). The customized implementation of the Python-based implicit solvers resulted in approximately consistent solutions with the Python-based matrix exponential method ( expm ). The differences generally observed between the implicit solvers and expm are attributable to numerical precision and the order of numerical approximation of the numerical solvers. This study provides the first analysis of a list of Python ODE solvers and methods by comparing their usage for solving biokinetic models using the ICRP Publication 130 HRTM and provides a framework for the selection of the most appropriate ODE solvers and methods in Python language to implement for modeling the distribution of internal radioactivity.

61 RADIATION PROTECTION AND DOSIMETRY↗

NASA NTRS1990

The effect of steady aerodynamic loading on the flutter stability of turbomachinery blading

An aeroelastic analysis is presented which accounts for the effect of steady aerodynamic loading on the aeroelastic stability of a cascade of compressor blades. The aeroelastic model is a two degree of freedom model having bending and torsional displacements. A linearized unsteady potential flow theory is used to determine the unsteady aerodynamic response coefficients for the aeroelastic analysis. The steady aerodynamic loading was caused by the addition of airfoil thickness and camber and steady flow incidence. The importance of steady loading on the airfoil unsteady pressure distribution is demonstrated. Additionally, the effect of steady loading on the tuned flutter behavior and flutter boundaries indicates that neglecting either airfoil thickness, camber or incidence could result in nonconservative estimates of flutter behavior.

Smith, Todd E.↗

NASA NTRS1991

The effect of steady aerodynamic loading on the flutter stability of turbomachinery blading

Smith, Todd E.↗

DOE OSTI2024

Safety First: Stability and Dissipation of Line-tied Force-free Flux Tubes in Magnetized Coronae

Abstract Magnetized plasma columns and extended magnetic structures with both footpoints anchored to a surface layer are an important building block of astrophysical dissipation models. Current loops shining in X-rays during the growth of plasma instabilities are observed in the corona of the Sun and are expected to exist in highly magnetized neutron star magnetospheres and accretion disk coronae. For varying twist and system sizes, we investigate the stability of line-tied force-free flux tubes and the dissipation of twist energy during instabilities using linear analysis and time-dependent force-free electrodynamics simulations. Kink modes ( m = 1) and efficient magnetic energy dissipation develop for plasma safety factors q ≲ 1, where q is the inverse of the number of magnetic field line windings per column length. Higher-order fluting modes ( m > 1) can distort equilibrium flux tubes for q > 1 but induce significantly less dissipation. In our analysis, the characteristic pitch μ ˜ 0 of flux-tube field lines determines the growth rate ( ∝ μ ˜ 0 3 ) and minimum wavelength of the kink instability ( ∝ μ ˜ 0 − 1 ). We use these scalings to determine a minimum flux tube length for the growth of the kink instability for any given μ ˜ 0 . By drawing analogies to idealized magnetar magnetospheres with varying regimes of boundary shearing rates, we discuss the expected impact of the pitch-dependent growth rates for magnetospheric dissipation in magnetar conditions.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1972

Research on combustion instability and application to solid propellant rocket motors. II.

Review of the current state of analyses of combustion instability in solid-propellant rocket motors, citing appropriate measurements and observations. The work discussed has become increasingly important, both for the interpretation of laboratory data and for predicting the transient behavior of disturbances in full-scale motors. Two central questions are considered - namely, linear stability and nonlinear behavior. Several classes of problems are discussed as special cases of a general approach to the analysis of combustion instability. Application to motors, and particularly the limitations presently understood, are stressed.

Culick, F. E. C.↗

NASA NTRS1981

Comparison of closed loop model with flight test results

An analytic technique capable of predicting the landing characteristics of proposed aircraft configurations in the early stages of design was developed. In this analysis, a linear pilot-aircraft closed loop model was evaluated using experimental data generated with the NT-33 variable stability in-flight simulator. The pilot dynamics are modeled as inner and outer servo loop closures around aircraft pitch attitude, and altitude rate-of-change respectively. The landing flare maneuver is of particular interest as recent experience with military and other highly augmented vehicles shows this task to be relatively demanding, and potentially a critical design point. A unique feature of the pilot model is the incorporation of an internal model of the pilot's desired flight path for the flare maneuver.

George, F. L.↗

NASA NTRS1991

Effect of periodic accelerations on interface stability in a multilayered fluid configuration

The increasing number of research opportunities in a microgravity environment will benefit not only fundamental studies in fluid dynamics, but also technological applications such as those involving materials processing. In particular, fluid configurations that involve fluid-fluid interfaces would occur in a variety of experimental investigations. This work investigates the stability of a configuration involving fluid-fluid interfaces in the presence of a time-dependent (periodic) forcing. The fluid configuration is multilayered and infinite in extent. The analysis is linear and inviscid, and the acceleration vector is oriented perpendicular to each interface. A Floquent analysis is employed, and the resulting algebraic eigensystem is truncated. Nondimensional parameters appear in the algebraic system. A numerical study is performed to elucidate the regions of instability and the effect of parameter variation on the fluid configuration stability.

Lyell, M. J.↗

NASA NTRS1974

Maximum likelihood identification of aircraft stability and control derivatives

Application of a generalized identification method to flight test data analysis. The method is based on the maximum likelihood (ML) criterion and includes output error and equation error methods as special cases. Both the linear and nonlinear models with and without process noise are considered. The flight test data from lateral maneuvers of HL-10 and M2/F3 lifting bodies are processed to determine the lateral stability and control derivatives, instrumentation accuracies, and biases. A comparison is made between the results of the output error method and the ML method for M2/F3 data containing gusts. It is shown that better fits to time histories are obtained by using the ML method. The nonlinear model considered corresponds to the longitudinal equations of the X-22 VTOL aircraft. The data are obtained from a computer simulation and contain both process and measurement noise. The applicability of the ML method to nonlinear models with both process and measurement noise is demonstrated.

Mehra, R. K.↗