Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

The stability of motion of satellites with flexible appendages

The mathematical formulation associated with the problem of stability of motion of a satellite consisting of a main rigid body and three (or less) pairs of flexible rods is presented. The rods are capable of flexure in two orthogonal directions. Whereas the rotational motion of the body is described by generalized coordinates depending on time alone, the elastic displacements of the rods depend both on spatial position and time. Assuming no external torques, there exist motion integrals in the form of momentum integrals. These integrals can be regarded as constraint equations relating the system velocities, and used to reduce the number of variables describing the motion. The stability analysis has been carried out by means of an extension of the Liapunov direct method. Since the elastic vibrations result in energy dissipation, it is shown that the equilibrium position is asymptotically stable if the Hamiltonian is positive definite and unstable if it can take negative values in the neighborhood of the equilibrium. Determining the sign definiteness of the Hamiltonian is complicated by the fact that it contains spatial derivatives of the elastic displacements. Two methods are presented to cope with this problem. The first, the standard modal analysis in conjunction with series truncation, develops criteria in terms of infinite series associated with the natural modes and frequencies of the elastic rods. The second, the method of integral coordinates, yields closed-form stability criteria involving the system parameters, such as the body moments of inertia, the length and mass distribution of the rods, the lowest natural frequencies of the rods, and the satellite spin velocity.

Meirovitch, L.↗
DOE OSTI2021

Mitigating Injection‐Induced Seismicity Along Basement Faults by Extraction: Application to 2016–2018 Pohang Earthquakes

Abstract The risk of induced earthquake is considered as the major obstacle for sustainable subsurface energy activities. As a mitigation approach, mechanistic risk assessment can be essential to minimize the potential of seismic hazard. Coupled hydro‐mechanical modeling approach is integrated with geomechanical analysis to better assess the occurrence of seismicity induced by injection and/or extraction, which will provide a basis to update operational conditions (e.g., number and location of mitigating wells, rate and duration of operations) for mitigating earthquake nucleation. The Coulomb stability analysis and seismicity rate estimates from this generic study indicate that simultaneous extraction can relieve injection‐induced pore‐pressure buildup, but also transfer and accumulate elastic energy along the faults. Field‐scale simulation with variation in well design for the Pohang EGS site, South Korea, show that multiple wells aligned across a preexisting fault form gradients of pressure and stress fields normal to the fault, which can enhance the direct impact of diffusion process and poroelastic stressing on the fault stability perturbation. Our mechanistic study suggests multiphysics‐oriented mitigation strategies that consider (1) possible mechanical deformation, caused by hydro‐mechanical coupling effects, and (2) proper configuration of well operations based on site‐specific geological characteristics.

Chang, K. W.↗
NASA NTRS1982

Stability of mixed time integration schemes for transient thermal analysis

A current research topic in coupled-field problems is the development of effective transient algorithms that permit different time integration methods with different time steps to be used simultaneously in various regions of the problems. The implicit-explicit approach seems to be very successful in structural, fluid, and fluid-structure problems. This paper summarizes this research direction. A family of mixed time integration schemes, with the capabilities mentioned above, is also introduced for transient thermal analysis. A stability analysis and the computer implementation of this technique are also presented. In particular, it is shown that the mixed time implicit-explicit methods provide a natural framework for the further development of efficient, clean, modularized computer codes.

Liu, W. K.↗
NASA NTRS1992

The effect of gravity modulation on thermosolutal convection in an infinite layer of fluid

The effect of time-periodic vertical gravity modulation on the onset of thermosolutal convection in an infinite horizontal layer with stress-free boundaries is investigated using Floquet theory for the linear stability analysis. Situations for which the fluid layer is stably stratified in either the fingering or diffusive regimes of double-diffusive convection are considered. Results are presented both with and without steady background acceleration. Modulation may stabilize an unstable base solution or destabilize a stable base solution. In addition to synchronous and subharmonic response to the modulation frequency, instability in the double diffusive system can occur via a complex conjugate mode. In the diffusive regime, where oscillatory onset occurs in the unmodulated system, regions of resonant instability occur and exhibit strong coupling with the unmodulated oscillatory frequency. The response to modulation of the fundamental instability of the unmodulated system is described both analytically and numerically; in the double-diffusive system this mode persists under subcritical conditions as a high-frequency lobe.

Saunders, B. V.↗
NASA NTRS1992

Granular instability in fluidized beds: A small-over-large instability

For different particle diameters arranged in a gas fluidized bed, a new gravitational instability is predicted. The first order stability analysis yields conditions for growth of the most unstable wavenumber both as a function of particle diameter and dimensionless Bond number. Gravity acts to distribute the greater fluid mass associated with smaller particles, while surface tension stabilizes the entire bed against accelerations. The phenomenon, called the 'Small-over-Large' (SL) instability, evolves somewhat akin to a Rayleigh-Taylor mechanism and in addition to fundamental interest, carries implications for predicting solute redistribution and gravitational alignments in particle beds and chemical separation columns.

Noever, David A.↗
NASA NTRSDate not supplied

Impedance Measurements of Motor Drives and Supplies in NASA NEAT Facility

Many electrified aircraft configurations consist of a number of power components operating on the same shared electrical bus. Understanding the impedance performance of electrical loads and sources is key to understanding, and designing for, acceptable overall vehicle power quality and power system stability. This paper discusses impedance measurements made on the NASA Electric Aircraft Testbed (NEAT) Electrical Power System (EPS) at NASA Glenn Research Center. First, an overall discussion of the NEAT facility configuration during testing is presented; then details of the impedance measurement approach and the various testing configurations are discussed. Next, input impedance measurements at NEAT, including sources (DC supplies) and loads (multiple motor drives, electric machines, and resistive load banks, with long interconnecting cable leads), are performed under a number of conditions, and results, including load model comparisons and stability analysis, are presented and discussed.

electrified aircraft propulsion↗
DOE OSTI2026

Validation of the GFS model for gyrokinetic stability of NSTX pedestal data

This study presents a large database validation of the gyro fluid system (GFS) model for linear gyrokinetic stability for high-mode (H-mode) edge transport barrier conditions in the national spherical torus experiment (NSTX) tokamak. The database of linear stability calculations with the CGYRO gyrokinetic code was produced using plasma profile measurements from NSTX discharges to identify kinetic ballooning modes (KBM), trapped electron modes (TEM), and micro-tearing modes (MTM) that limit the pressure profile gradient in the H-mode barrier. A novel Bayesian optimization approach determines optimal resolution parameters for GFS specifically for spherical tokamak pedestal conditions. Our results demonstrate that GFS, with optimized resolution, can achieve accurate linear stability analysis in NSTX pedestal conditions for reduced resolution compared to CGYRO. GFS can accurately find the KBM, TEM, and MTM instability branches. Parametric analysis reveals that GFS accuracy in this extreme pedestal parameter range is degraded for low magnetic shear and near the separatrix conditions. These findings establish GFS as a fast linear eigenmode solver for spherical tokamak pedestal gyrokinetic stability and demonstrate a systematic methodology for determining the optimum resolution settings.

Yang, Minglei [Oak Ridge National Laboratory (ORNL↗
NASA NTRS1990

On the stability of swept wing laminar boundary layers including curvature effects

A linear stability analysis that encompasses curvature effects has been conducted in wind tunnel experiments on a swept NACA 64(2)-A015 wing, and published transition-onset results have been correlated with computed N-factor values. A strong stabilizing influence is noted upon the growth of the crossflow disturbance, when the flow is accelerated in regions of high body curvature. The maximum amplified crossflow disturbances were in all cases travelling waves; when TS waves reached their maximum, the N-factors at transition lay in the 9.9-13.8 range. Stabilization due to curvature effects was less pronounced in cases where acceleration occurred over a large portion of chord.

Viken, J.↗
NASA NTRS1983

Mixed time integration schemes for transient conduction forced-convection analysis

A partition procedure for forced-convection conduction transient problems is presented. Mixed time partitions are defined wherein coupled conduction force-matrix equations are discretized using an implicit integration method, followed by derivation of a mixed time integration technique. Explicit-implicit and explicit-explicit partitions are performed for a stability analysis for transient conditions, e.g., those found in an actively air-cooled engine and airframe structure.

Liu, W. K.↗
NASA NTRS1995

Hydrodynamic Stability of Multicomponent Droplet Gasification in Reduced Gravity

This investigation addresses the problem of hydrodynamic stability of a two-component droplet undergoing spherically-symmetrical gasification. The droplet components are assumed to have characteristic liquid species diffusion times that are large relative to characteristic droplet surface regression times. The problem is formulated as a linear stability analysis, with a goal of predicting when spherically-symmetric droplet gasification can be expected to be hydrodynamically unstable from surface-tension gradients acting along the surface of a droplet which result from perturbations. It is found that for the conditions assumed in this paper (quasisteady gas phase, no initial droplet temperature gradients, diffusion-dominated gasification), surface tension gradients do not play a role in the stability characteristics. In addition, all perturbations are predicted to decay such that droplets were hydrodynamically stable. Conditions are identified, however, that deserve more analysis as they may lead to hydrodynamic instabilities driven by capillary effects.

Aharon, I.↗
NASA NTRS2020

Streak Instability Analysis for BOLT Configuration

The BOLT configuration designed for an upcoming flight experiment exemplifies the complexities of hypersonic boundary layer transition in the presence of multiple and potentially interacting instability mechanisms. A numerical study of the boundary layer instabilities over this configuration has been performed with the goal of applying state-of-the-art tools in multidimensional stability analysis that complement the other ongoing research related to the BOLT flight campaign. Specifically, the paper outlines modal instability characteristics of the boundary layer streaks adjacent to the minor-axis symmetry plane of a 1/3 rd scale BOLT configuration that was tested in the low-disturbance hypersonic wind tunnel at Purdue University and both quiet and conventional tunnels at Texas A&M University (TAMU). The basic state computations agree well with the heat flux measurements in the Purdue quiet tunnel as well as with the off-body velocity field described in a previous DNS study at the University of Minnesota. The encouraging agreement between the predicted frequency spectra of the streak instability in the region of increased boundary layer thickness and the surface pressure spectra measured during the TAMU experiments suggests that these streak instabilities may play an important role during the transition process near the symmetry plane. The modest value of the predicted peak N-factor (N max ≈ 5) for the streak instabilities provides a potential explanation for why transition is not observed under quiet flow conditions, but may occur in a noisy wind tunnel. Similar instability analyses are performed for the flight configuration at two selected flow conditions from the ascent portion of the planned flight trajectory. The predicted values of the peak N-factor suggest that transition onset in the vicinity of the center plane is possible at the ascent time of t = 26.98 seconds (N max ≈ 12.5 at X/L = 0.77), but the streak instabilities weaken sufficiently at the later time of t = 28.88 seconds (N max ≈ 11.5 at X/L ≈ 1.00) and are less likely to cause transition at that time.

Laminar turbulent transition↗
NASA NTRS1984

Theoretical analyses of baroclinic flows

Completed and ongoing research activities are discussed briefly, including a three-dimensional, linear stability analysis of the baroclinic Hadley cell and a numerical model of the baroclinic flow between two rotating concentric spheres. This model simulates axisymmetric flow in the Atmospheric General Circulation Experiment configuration. A computer code designed to solve the strongly nonlinear stability problem for the Eady basic state is mentioned.

Antar, B. N.↗
NASA NTRS1987

A cyclotron-maser instability associated with a nongyrotropic distribution

A stability analysis for the cyclotron-maser instability in the presence of a nongyrotropic electron distribution is presented. The model configuration describes a uniformly magnetized cold ambient plasma that contains a relatively diffuse suprathermal electron component coherently bunched in gyrophase. The stability of perturbations propagating parallel to the ambient magnetic field is considered, and substantial growth rates are found to occur. The results are contrasted with those found for a comparable gyrotropic loss-cone distribution, and it is found that the nongyrotropic instability is characterized by substantially higher growth rates.

Freund, H. P.↗
NASA NTRS1990

Adaptive quadrilateral and triangular finite-element scheme for compressible flows

The development of an adaptive mesh refinement procedure for analyzing high-speed compressible flows using the finite-element method is described. This new adaptation procedure, which uses both quadrilateral and triangular elements, was implemented with two explicit finite-element algorithms - the two-step Taylor-Galerkin and the multistep Galerkin-Runge-Kutta schemes. A von Neumann stability analysis and a rotating 'cosine hill'problem demonstrate the instability of the Taylor-Galerkin scheme when coupled with the adaptation procedure. For the same adaptive refinement scheme, the Galerkin-Runge-Kutta procedure yields stable solutions within its explicit stability limit. The utility of this new adaptation procedure for the prediction of compressible flow features is illustrated for inviscid problems involving strong shock interactions at hypersonic speeds.

Ramakrishnan, R.↗
NASA NTRS1991

Robust control of systems with real parameter uncertainty and unmodelled dynamics

During this research period we have made significant progress in the four proposed areas: (1) design of robust controllers via H infinity optimization; (2) design of robust controllers via mixed H2/H infinity optimization; (3) M-delta structure and robust stability analysis for structured uncertainties; and (4) a study on controllability and observability of perturbed plant. It is well known now that the two-Riccati-equation solution to the H infinity control problem can be used to characterize all possible stabilizing optimal or suboptimal H infinity controllers if the optimal H infinity norm or gamma, an upper bound of a suboptimal H infinity norm, is given. In this research, we discovered some useful properties of these H infinity Riccati solutions. Among them, the most prominent one is that the spectral radius of the product of these two Riccati solutions is a continuous, nonincreasing, convex function of gamma in the domain of interest. Based on these properties, quadratically convergent algorithms are developed to compute the optimal H infinity norm. We also set up a detailed procedure for applying the H infinity theory to robust control systems design. The desire to design controllers with H infinity robustness but H(exp 2) performance has recently resulted in mixed H(exp 2) and H infinity control problem formulation. The mixed H(exp 2)/H infinity problem have drawn the attention of many investigators. However, solution is only available for special cases of this problem. We formulated a relatively realistic control problem with H(exp 2) performance index and H infinity robustness constraint into a more general mixed H(exp 2)/H infinity problem. No optimal solution yet is available for this more general mixed H(exp 2)/H infinity problem. Although the optimal solution for this mixed H(exp 2)/H infinity control has not yet been found, we proposed a design approach which can be used through proper choice of the available design parameters to influence both robustness and performance. For a large class of linear time-invariant systems with real parametric perturbations, the coefficient vector of the characteristic polynomial is a multilinear function of the real parameter vector. Based on this multilinear mapping relationship together with the recent developments for polytopic polynomials and parameter domain partition technique, we proposed an iterative algorithm for coupling the real structured singular value.

Chang, Bor-Chin↗