Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stability properties”

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 235 records · Page 13

Nonlinear development and secondary instability of Gortler vortices in hypersonic flows

In a hypersonic boundary layer over a wall of variable curvature, the region most susceptible to Goertler vortices is the temperature adjustment layer over which the basic state temperature decreases monotonically to its free stream value. Except for a special wall curvature distribution, the evolution of Goertler vortices trapped in the temperature adjustment layer will in general be strongly affected by the boundary layer growth through the O(M sup 3/2) curvature of the basic state, where M is the free stream Mach number. Only when the local wavenumber becomes as large as of order M sup 3/8, do nonparallel effects become negligible in the determination of stability properties. In the latter case, Goertler vortices will be trapped in a thin layer of O(epsilon sup 1/2) thickness which is embedded in the temperature adjustment layer; here epsilon is the inverse of the local wavenumber. A weakly nonlinear theory is presented in which the initial nonlinear development of Goertler vortices in the neighborhood of the neutral position is studied and two coupled evolution equations are derived. From these, it can be determined whether the vortices are decaying or growing depending on the sign of a constant which is related to wall curvature and the basic state temperature.

Fu, Yibin B.↗

Collisionless drift-tearing modes in the magnetopause

The linear stability properties of collisionless drift-tearing modes are analyzed in a modified Harris equilibrium model of the magnetopause. Particular attention is paid to the relevance of the parametric behavior of growth rates to the 'magnetic percolation' theory of flux transfer event formation (Galeev et al., 1986). Numerical methods are used to solve the drift-tearing eigenmode equations and the results are compared with those previously obtained by analytical methods. The analytical results are found to correctly model important parametric dependencies but to typically overestimate the rate of growth. The eigenmode equations are numerically difficult, and an integration scheme utilizing Ricatti transforms is developed to affect their solution.

Gladd, N. T.↗

The inviscid stability of supersonic flow past a sharp cone

The effects of lateral curvature on the development of supersonic laminar inviscid boundary-layer flow on a sharp cone with adiabatic wall conditions are investigated analytically, with a focus on the linear temporal inviscid stability properties. The derivation of the governing equations and of a 'triply generalized' inflexion condition is outlined, and numerical results for freestream Mach number 3.8 are presented in extensive graphs and characterized in detail. A third instability mode related to the viscous mode observed by Duck and Hall (1990) using triple-deck theory is detected and shown to be more unstable and to have larger growth rates than the second mode in some cases. It is found that the 'sonic' neutral mode is affected by the lateral curvature and becomes a supersonic neutral mode.

Duck, Peter W.↗

Application of Krylov exponential propagation to fluid dynamics equations

This paper presents an application of matrix exponentiation via Krylov subspace projection, to the solution of fluid dynamics problems. The main idea is to approximate the operation exp(A)v by means of a projection-like process onto a Krylov subspace. This results in a computation of an exponential matrix vector product similar to the one above but of a much smaller size. Time integration schemes can then be devised to exploit this basic computational kernel. The motivation of this approach is to provide time-integration schemes that are essentially of an explicit nature but which have good stability properties.

Saad, Y.↗

Thermal phases of the interstellar medium in galaxies

This review deals with the theory of multiphase media in astrophysical systems. The basic reasons for the existence of multiple thermal phases, and the fundamental connection between multiphase media and thermal instability are discussed. Important examples of multiphase media are described, and the interactions among phases, i.e., mass exchange driven by thermal conduction and hydrodynamic ablation are examined. Mass exchange may compete with radiative heating and cooling for control of the thermal state of the hot phase, and may alter the thermal stability properties of the system.

Begelman, Mitchell C.↗

Application of Krylov exponential propagation to fluid dynamics equations

An application of matrix exponentiation via Krylov subspace projection to the solution of fluid dynamics problems is presented. The main idea is to approximate the operation exp(A)v by means of a projection-like process onto a krylov subspace. This results in a computation of an exponential matrix vector product similar to the one above but of a much smaller size. Time integration schemes can then be devised to exploit this basic computational kernel. The motivation of this approach is to provide time-integration schemes that are essentially of an explicit nature but which have good stability properties.

Saad, Youcef↗

Ion Beam Textured and Coated Surfaces Experiment (IBEX)

Ion beam textured and commercial materials suitable for use in space power systems were flown in low Earth orbit on the Long Duration Exposure Facility (LDEF) for 5.8 years. Because of their location on LDEF (98 deg from the ram direction), the 36 materials were primarily exposed to vacuum ultraviolet radiation, thermal cycling, the vacuum of space, the micrometeoroid environment, and grazing incidence atomic oxygen. Measurements of solar absorptance and thermal emittance (pre- and post-flight) showed no changes for almost all of the materials, except for the S-13G and Kapton and coated Kapton samples. The optical property stability of ion beam textured surfaces and most other surfaces indicates that they are functionally durable to the synergistic rigors of the space environment.

Mirtich, Michael J.↗

Singularity perturbed zero dynamics of nonlinear systems

Stability properties of zero dynamics are among the crucial input-output properties of both linear and nonlinear systems. Unstable, or 'nonminimum phase', zero dynamics are a major obstacle to input-output linearization and high-gain designs. An analysis of the effects of regular perturbations in system equations on zero dynamics shows that whenever a perturbation decreases the system's relative degree, it manifests itself as a singular perturbation of zero dynamics. Conditions are given under which the zero dynamics evolve in two timescales characteristic of a standard singular perturbation form that allows a separate analysis of slow and fast parts of the zero dynamics.

Isidori, A.↗

Implicit treatment of diffusion terms in lower-upper algorithms

A method is presented which allows diffusion terms to be treated implicitly in the lower-upper (LU) algorithm (which is a commonly used method for solving 'compressible' Euler and Navier-Stokes equations) so that the algorithm's good stability properties will not be impaired. The new method generalizes the concept of LU factorization from that associated with the sign of eigenvalues to that associated with backward- and forward-difference operators without regard to eigenvalues. The method is verified in a turbulent boundary layer study.

Shih, T. I.-P.↗

An implicit method for the nonlinear modelling and simulation of piezoceramic actuators displaying hysteresis

Experiments have shown that piezoceramic materials display a nonlinear relationship between the applied electric field and the actuation strain. This relationship also displays a substantial hysteresis upon reversal of the applied field. In this paper, piezoceramic actuator models are incorporated into the structural equations of motion to arrive at a set of nonlinear actively controlled structural equations of motion. A new implicit algorithm for determining the time history of the actively controlled structure is presented. The algorithm employs the trapezoidal rule for stepping the equations forward in time. The algorithm is compared to an explicit algorithm and is shown to provide greater numerical accuracy. In addition, the numerical stability and convergence characteristics are presented via example. The new algorithm displays convergence and stability properties that are comparable to the standard trapezoidal rule used for dynamic systems without hysteresis.

Leigh, Timothy D.↗

Hypersonic stability and transition

An evaluation is conducted of the development status of growing normal modes in the computation of flow fields over hypersonic vehicles, giving attention to the effects of surface cooling, pressure gradients, and suction on the growth rates of the normal modes. Also noted are the effects of bluntness on the stability properties of hypersonic vehicles, and the efficacy of the e exp N method of transition estimation at supersonic and hypersonic speeds. A survey is made of the precautions that must be taken to ensure the effectiveness of experimental and computational studies of hypersonic flow phenomena.

Reshotko, Eli↗

Ion mixing in the plasma sheet boundary layer by drift instabilities

The linear stability properties of collisionless drift instabilities are analyzed in a Harris equilibrium model of the plasma sheet boundary layer (PSBL). The strearmng ions with drift-type instabilities driven in the PSBL are considered. The fluid approximation leads to growth but predicts that the mode width approaches the gyroradius of the energetic ions. Thus an integral equation theory for the modes is developed taking into account that in the PSBL the curvature drift is weak compared with the grad-B drift. The exact wave particle resonance is kept in the nonlocal response functions. Plasma density, temperature, and magnetic gradient drift motions are taken into account. The drift modes produce an anomalous cross-field momentum transport mixing the PSBL ions on the time scale of tens of seconds. A nonlinear simulation is performed which shows the coalescence of the small scale, fast growing modes into large-scale vortices. The relation between these collective modes and plasma sheet transport phenomena is discussed including the comparison with the competing plasma mixing from single-particle stochasticity.

Horton, W.↗

Effect of nonzero surface admittance on receptivity and stability of compressible boundary layer

The effect of small-amplitude short-scale variations in surface admittance on the acoustic receptivity and stability of two-dimensional compressible boundary layers is examined. In the linearized limit, the two problems are shown to be related both physically and mathematically. This connection between the two problems is used, in conjunction with some previously reported receptivity results, to infer the modification of stability properties due to surface permeability. Numerical calculations are carried out for a self-similar flat-plate boundary layer at subsonic and low supersonic speeds. Variations in mean suction velocity at the perforated admittance surface can also induce receptivity to an acoustic wave. For a subsonic boundary layer, the dependence of admittance-induced receptivity on the acoustic-wave orientation is significantly different from that of the receptivity produced via mean suction variation. The admittance-induced receptivity is generally independent of the angle of acoustic incidence, except in a relatively narrow range of upstream-traveling waves for which the receptivity becomes weaker. However, this range of angles is precisely that for which the suction-induced receptivity tends to be large. At supersonic Mach numbers, the admittance-induced receptivity to slow acoustic models is relatively weaker than that in the case of the fast acoustic modes. We also find that purely real values for the surface admittance tend to have a destabilizing effect on the evolution of an instability wave over a slightly permeable surface. The limits on the validity of the linearized approximation are also assessed in one specific case.

Choudhari, Meelan↗

Equilibrium, stability, and orbital evolution of close binary systems

We present a new analytic study of the equilibrium and stability properties of close binary systems containing polytropic components. Our method is based on the use of ellipsoidal trial functions in an energy variational principle. We consider both synchronized and nonsynchronized systems, constructing the compressible generalizations of the classical Darwin and Darwin-Riemann configurations. Our method can be applied to a wide variety of binary models where the stellar masses, radii, spins, entropies, and polytropic indices are all allowed to vary over wide ranges and independently for each component. We find that both secular and dynamical instabilities can develop before a Roche limit or contact is reached along a sequence of models with decreasing binary separation. High incompressibility always makes a given binary system more susceptible to these instabilities, but the dependence on the mass ratio is more complicated. As simple applications, we construct models of double degenerate systems and of low-mass main-sequence star binaries. We also discuss the orbital evoltuion of close binary systems under the combined influence of fluid viscosity and secular angular momentum losses from processes like gravitational radiation. We show that the existence of global fluid instabilities can have a profound effect on the terminal evolution of coalescing binaries. The validity of our analytic solutions is examined by means of detailed comparisons with the results of recent numerical fluid calculations in three dimensions.

Lai, Dong↗

Hydrodynamics of binary coalescence. 2: Polytropes with gamma = 5/3

We present a new numerical study of the equilibrium and stability properties of close binary systems. We use the smoothed-particle hydrodynamics (SPH) technique both to construct accurate equilibrium configurations in three dimensions and to follow their hydrodynamic evolution. We adopt a simple polytropic equation of state p = K(sub rho)(exp gamma) with gamma = 5/3 and K = constant within each star, applicable to low-mass degenerate dwarfs as well as low-mass main-sequence stars. For degenerate configurations, we set the two polytropic constants equal, K = K prime, independent of the mass ratio. For main-sequence stars, we adjust K and K prime so as to obtain a simple mass-radius relation of the form R/R prime = M/M prime, where R prime and M prime are the radius and mass of the secondary. Along a sequence of binary equilibrium configurations for two identical stars, we demonstrate the existence of both secular and dynamical instabilities, confirming directly the results of recent analytic work. We use the SPH method to calculate the nonlinear development of the dynamical instability and to determine the final fate of the system. We find that the two stars merge together into a single, rapidly rotating object in just a few orbital periods. Equilibrium sequences are also constructed for systems containing two nonidentical stars. These sequences terminate at a Roche limit, which we can determine very accurately using SPH. For two low-mass main-sequence stars with mass ratio q approximately less than 0.4 we find that the (synchronized) Roche limit configuration is secularly unstable. We discuss the implications of our results for the evolution of double white-dwarf systems and W Ursae Majoris binaries.

Rasio, Frederic A.↗

Building a better leapfrog

In stellar dynamical computer simulations, as well as other types of simulations using particles, time step size is often held constant in order to guarantee a high degree of energy conservation. In many applications, allowing the time step size to change in time can offer a great saving in computational cost, but variable-size time steps usually imply a substantial degradation in energy conservation. We present a meta-algorithm' for choosing time steps in such a way as to guarantee time symmetry in any integration scheme, thus allowing vastly improved energy conservation for orbital calculations with variable time steps. We apply the algorithm to the familiar leapfrog scheme, and generalize to higher order integration schemes, showing how the stability properties of the fixed-step leapfrog scheme can be extended to higher order, variable-step integrators such as the Hermite method. We illustrate the remarkable properties of these time-symmetric integrators for the case of a highly eccentric elliptical Kepler orbit and discuss applications to more complex problems.

Hut, Piet↗

A new criterion for bar-forming instability in rapidly rotating gaseous and stellar systems. 2: Nonaxisymmetric form

We have previously introduced the parameter alpha as an indicator of stability to m = 2 nonaxisymmetric modes in rotating, self-gravitating, axisymmetric, gaseous (alpha less than or approximately equal to 0.34) and stellar (alpha less than or approximately equal to 0.25) systems. This parameter can be written as alpha = (ft/2)(exp 1/2), where t is defined as T/(absolute value of W), T is the total rotational kinetic energy, W is the total gravitational potential energy, and f is a function characteristic of the topology/connectedness and the geometric shape of a system. In this paper, we extend the stability criterion to nonaxisymmetric equilibrium systems by determining empirically the appropriate form of the function f for ellipsoids and elliptical disks and cylinders. We test the validity of this extension of the stability indicator alpha by considering its predictions for previously published, gaseous and stellar, nonaxisymmetric models. The formulation and critical values account accurately for the stability properties of m =2 modes in gaseous Riemann S-type ellipsoids (including the Jacobi and Dedekind ellipsoids) and elliptical Riemann disks as well as in stellar elliptical Freeman disks and cylinders: all these systems are dynamically stable except the stellar elliptical Freeman disks that exhibit a relatively small region of m = 2 dynamical instability. A partial disagreement in the case of stellar Freeman ellipsoids in maximum rotation may be due to the fact that the region of instability has not been previously determined with sufficient accuracy.

Christodoulou, Dimitris M.↗

Active Control of Wind Tunnel Noise

The need for an adaptive active control system was realized, since a wind tunnel is subjected to variations in air velocity, temperature, air turbulence, and some other factors such as nonlinearity. Among many adaptive algorithms, the Least Mean Squares (LMS) algorithm, which is the simplest one, has been used in an Active Noise Control (ANC) system by some researchers. However, Eriksson's results, Eriksson (1985), showed instability in the ANC system with an ER filter for random noise input. The Restricted Least Squares (RLS) algorithm, although computationally more complex than the LMS algorithm, has better convergence and stability properties. The ANC system in the present work was simulated by using an FIR filter with an RLS algorithm for different inputs and for a number of plant models. Simulation results for the ANC system with acoustic feedback showed better robustness when used with the RLS algorithm than with the LMS algorithm for all types of inputs. Overall attenuation in the frequency domain was better in the case of the RLS adaptive algorithm. Simulation results with a more realistic plant model and an RLS adaptive algorithm showed a slower convergence rate than the case with an acoustic plant as a delay plant. However, the attenuation properties were satisfactory for the simulated system with the modified plant. The effect of filter length on the rate of convergence and attenuation was studied. It was found that the rate of convergence decreases with increase in filter length, whereas the attenuation increases with increase in filter length. The final design of the ANC system was simulated and found to have a reasonable convergence rate and good attenuation properties for an input containing discrete frequencies and random noise.

Hollis, Patrick↗