Engineering Papers⌕ Search

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.

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

At least 1,045 records · Page 58

Determination of Longitudinal Stability and Control Characteristics from Free-Flight Model Tests with Results at Transonic Speeds for Three Airplane Configurations

A test technique and data analysis method has been developed for determining the longitudinal aerodynamic characteristics from free-flight tests of rocket-propelled models. The technique makes use of accelerometers and an angle-of-attack indicator to permit instantaneous measurements of lift, drag, and pitching moments. The data, obtained during transient oscillations resulting from control-surface disturbances, are analyzed by essentially nonlinear direct methods (such as cross plots of the variation of lift coefficient with angle of attack) and by linear indirect methods by using the equations of motion for a transient oscillation. The analysis procedure has been set forth in some detail and the feasibility of the method has been demonstrated by data measured through the transonic speed range on several airplane configurations. It was shown that the flight conditions and dynamic similitude factors for the tests described were reasonably close to typical full-scale airplane conditions.

Gillis, Clarence L↗

Aeroelastic stability and response of horizontal axis wind turbine blades

The coupled flap-lag-torsion equations of motion of an isolated horizontal axis wind turbine blade are formulated. Quasi-steady blade-element strip theory was applied to derive the aerodynamic operator which includes boundary layer type gradient winds. The final equations which have periodic coefficients were solved in order to obtain the aeroelastic response and stability of large horizontal axis wind turbine blade. A new method of generating an appropriate time-dependent equilibrium position (required for the stability analysis) has been implemented. Representative steady-state responses and stability boundaries, applicable mainly to an existing blade design (NASA/-ERDA MOD-0), are presented. The results indicate that the MOD-0 configuration is a basically stable design and that blade stability is not sensitive to offsets between blade elastic axis and aerodynamic center. Blade stability appears to be sensitive to precone. The tower shadow (or wake) has a considerable effect on the flap response but leaves blade stability unchanged. Finally, it was found that non linear terms in the equations of motion can significantly affect the linearized stability boundaries, however, these terms have a negligible effect on blade response at operating conditions.

Kottapalli, S. B. R.↗

A modal aeroelastic analysis scheme for turbomachinery blading

An aeroelastic analysis is developed which has general application to all types of axial-flow turbomachinery blades. The approach is based on linear modal analysis, where the blade's dynamic response is represented as a linear combination of contributions from each of its in-vacuum free vibrational modes. A compressible linearized unsteady potential theory is used to model the flow over the oscillating blades. The two-dimensional unsteady flow is evaluated along several stacked axisymmetric strips along the span of the airfoil. The unsteady pressures at the blade surface are integrated to result in the generalized force acting on the blade due to simple harmonic motions. The unsteady aerodynamic forces are coupled to the blade normal modes in the frequency domain using modal analysis. An iterative eigenvalue problem is solved to determine the stability of the blade when the unsteady aerodynamic forces are included in the analysis. The approach is demonstrated by applying it to a high-energy subsonic turbine blade from a rocket engine turbopump power turbine. The results indicate that this turbine could undergo flutter in an edgewise mode of vibration.

Smith, Todd E.↗

Computation of Steady and Unsteady Laminar Flames: Theory

In this paper we describe the numerical analysis underlying our efforts to develop an accurate and reliable code for simulating flame propagation using complex physical and chemical models. We discuss our spatial and temporal discretization schemes, which in our current implementations range in order from two to six. In space we use staggered meshes to define discrete divergence and gradient operators, allowing us to approximate complex diffusion operators while maintaining ellipticity. Our temporal discretization is based on the use of preconditioning to produce a highly efficient linearly implicit method with good stability properties. High order for time accurate simulations is obtained through the use of extrapolation or deferred correction procedures. We also discuss our techniques for computing stationary flames. The primary issue here is the automatic generation of initial approximations for the application of Newton's method. We use a novel time-stepping procedure, which allows the dynamic updating of the flame speed and forces the flame front towards a specified location. Numerical experiments are presented, primarily for the stationary flame problem. These illustrate the reliability of our techniques, and the dependence of the results on various code parameters.

Hagstrom, Thomas↗

Quantifying Fractional Ground Cover on the Climate Sensitive High Plains Using AVIRIS and Landsat TM Data

The High Plains is an economically important and climatologically sensitive region of the United States and Canada. The High Plains contain 100,000 sq km of Holocene sand dunes and sand sheets that are currently stabilized by natural vegetation. Droughts and the larger threat of global warming are climate phenomena that could cause depletion of natural vegetation and make this region susceptible to sand dune reactivation. This thesis is part of a larger study that is assessing the effect of climate variability on the natural vegetation that covers the High Plains using Landsat 5 and Landsat 7 data. The question this thesis addresses is how can fractional vegetation cover be mapped with the Landsat instruments using linear spectral mixture analysis and to what accuracy. The method discussed in this thesis made use of a high spatial and spectral resolution sensor called AVIRIS (Airborne Visible and Infrared Imaging Spectrometer) and field measurements to test vegetation mapping in three Landsat 7 sub-scenes. Near-simultaneous AVIRIS images near Ft. Morgan, Colorado and near Logan, New Mexico were acquired on July 10, 1999 and September 30, 1999, respectively. The AVIRIS flights preceded Landsat 7 overpasses by approximately one hour. These data provided the opportunity to test spectral mixture algorithms with AVIRIS and to use these data to constrain the multispectral mixed pixels of Landsat 7. The comparisons of mixture analysis between the two instruments showed that AVIRIS endmembers can be used to unmix Landsat 7 data with good estimates of soil cover, and reasonable estimates of non-photosynthetic vegetation and green vegetation. Landsat 7 derived image endmembers correlate with AVIRIS fractions, but the error is relatively large and does not give a precise estimate of cover.

Warner, Amanda Susan↗

Supersonic Traveling Crossflow Wave Characteristics in Ground and Flight Tests

This paper continues analysis of data acquired in a series of supersonic crossflow transition experiments on a 65° swept- wing model tested in both wind tunnel and flight. Flow visualization data are summarized to provide context for the swept-wing instability mechanisms of interest. Unsteady pressure measurements obtained near the wing model leading edge (LE) are studied using 2nd- and 3rd-order statistical methods to extract traveling crossflow wave characteristics from these data. Comparisons between tunnel and flight data with boundary-layer (BL) Linear Stability Theory (LST) predictions provide insights into BL transition phenomena similarities and differences observed in tests using the same swept-wing model configurations. Comparisons of unsteady pressure data results for wing LE configurations with and without distributed-roughness-element (DRE) patterns installed, suggest traveling crossflow disturbance growth is reduced when this type of flow control is applied. A brief stud y of different background surface roughness levels on the wing LE suggests the potential for another flow control approach to promote laminar flow on swept wings using streamwise-biased surface finishes.

Lewis R Owens↗

Analysis of Aircraft, Radiosonde and Radar Observations in Cirrus Clouds Observed During FIRE II: The Interactions Between Environmental Structure, Turbulence and Cloud Microphysical Properties

Ways to determine the turbulence intensity and the horizontal variability in cirrus clouds have been investigated using FIRE-II aircraft, radiosonde and radar data. Higher turbulence intensities were found within some, but not all, of the neutrally stratified layers. It was also demonstrated that the stability of cirrus layers with high extinction values decrease in time, possibly as a result of radiative destabilization. However, these features could not be directly related to each other in any simple manner. A simple linear relationship was observed between the amount of horizontal variability in the ice water content and its average value. This was also true for the extinction and ice crystal number concentrations. A relationship was also suggested between the variability in cloud depth and the environmental stability across the depth of the cloud layer, which requires further investigation.

Smith, Samantha A.↗

Model Validation for Propulsion - On the TFNS and LES Subgrid Models for a Bluff Body Stabilized Flame

With advances in computational power and availability of distributed computers, the use of even the most complex of turbulent chemical interaction models in combustors and coupled analysis of combustors and turbines is now possible and more and more affordable for realistic geometries. Recent more stringent emission standards have enticed the development of more fuel-efficient and low-emission combustion system for aircraft gas turbine applications. It is known that the NOx emissions tend to increase dramatically with increasing flame temperature. It is well known that the major difficulty, when modeling the turbulence-chemistry interaction, lies in the high non-linearity of the reaction rate expressed in terms of the temperature and species mass fractions. The transport filtered density function (FDF) model and the linear eddy model (LEM), which both use local instantaneous values of the temperature and mass fractions, have been shown to often provide more accurate results of turbulent combustion. In the present, the time-filtered Navier-Stokes (TFNS) approach capable of capturing unsteady flow structures important for turbulent mixing in the combustion chamber and two different subgrid models, LEM-like and EUPDF-like, capable of emulating the major processes occurring in the turbulence-chemistry interaction will be used to perform reacting flow simulations of a selected test case. The selected test case from the Volvo Validation Rig was documented by Sjunnesson.

bluff body↗

Steady vibrations of wing of circular plan form

This paper treats the problem of determining the lift, moment, and induced drag of a thin wing of circular plan form in uniform incompressible flow on the basis of linearized theory. As contrasted to a similar paper by Kinner, in which the acceleration potential method was used, the present paper utilizes the concept of the velocity potential. Calculations of the lift and moment are presented for several deformed shapes. It is shown that considerable deviations exist between the strip theory analysis and the more exact theory. The lift, moment, and induced drag are also determined for a harmonically oscillatory circular plan form wing. As contrasted to a similar paper by Schade, in which the acceleration potential method was used, the present paper utilizes the concept of the velocity potential. Expressions for lift, moment, and induced drag are given and finally specialized to the case of a slowly oscillating circular wing.

FLUTTER - WINGS↗

On a variational formulation of the weakly nonlinear magnetic Rayleigh–Taylor instability

The magnetic-Rayleigh–Taylor (MRT) instability is a ubiquitous phenomenon that occurs in magnetically-driven Z-pinch implosions. It is important to understand this instability since it can decrease the performance of such implosions. In this work, I present a theoretical model for the weakly nonlinear MRT instability. I obtain such a model by asymptotically expanding an action principle, whose Lagrangian leads to the fully nonlinear MRT equations. After introducing a suitable choice of coordinates, I show that the theory can be cast as a Hamiltonian system, whose Hamiltonian is calculated up to the sixth order in a perturbation parameter. The resulting theory captures the harmonic generation of MRT modes. It is shown that the amplitude at which the linear magnetic-Rayleigh–Taylor instability exponential growth saturates depends on the stabilization effect of the magnetic-field tension. Overall, the theory provides an intuitive interpretation of the weakly nonlinear MRT instability and provides a systematic approach for studying this instability in more complex settings.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Global linear and nonlinear gyrokinetic simulations of tearing modes

Abstract To better understand multi-scale interactions between global tearing modes and microturbulence in the Madison Symmetric Torus (MST) reversed-field pinch (RFP), the global gyrokinetic code Gene is modified to describe global tearing mode instability via a shifted Maxwellian (SM) distribution consistent with experimental equilibria. The implementation of the SM is tested and benchmarked by comparisons with different codes and models. Good agreement is obtained in code-code and code-theory comparisons. Linear stability of tearing modes of a non-reversed MST discharge is studied. A collisionality scan is performed to the lowest order unstable modes ( n = 5, n = 6) and shown to behave consistently with theoretical scaling. The nonlinear evolution is simulated, and saturation is found to arise from mode coupling and transfer of energy from the most unstable tearing mode to small-scale stable modes mediated by the m = 2 tearing mode. The work described herein lays the foundation for nonlinear simulation and analysis of the interaction of tearing modes and gyroradius-scale instabilities in RFP plasmas.

Physics↗

Handbook of structural stability part III : buckling of curved plates and shells

Available theories and test data on buckling of curved plates and shells are reviewed. For torsion and external-pressure loadings, the test data are correlated in terms of linear buckling theories for both the elastic and inelastic ranges. The cases which exhibit a marked disagreement between linear theory and test data include those of curved plates and cylinders under axial compression, cylinders under bending, and spherical plates under external pressure. These cases have been analyzed by a unified semiempirical approach for both the elastic and inelastic ranges which is satisfactory for analysis and design purposes. The effects of internal pressure on buckling of elements under uniaxial loads are discussed and data on various combined loadings are presented in interaction form. (author)

Gerard, George↗

Spectral Analysis of Regular Material Point Method and its Application to Study High Pressure Reverse Osmosis Membrane Compaction and Embossing

Material Point Method (MPM) is gaining widespread interest in applied continuum mechanics. The fact that all the continuum properties are stored on the particles (or material points) and the governing equations are solved on these material points makes MPM extremely suited to problems involving severe material deformations, such as crack propagation, soil movement, and fluid flows. Despite its popularity, only a few studies have focused on the numerical properties of MPM. This presentation introduces a global spectral analysis of the regular material point method. Contrary to previous studies, the analysis focuses on the numerical properties of the method in the spectral space. The amplification factor is derived as a function of the non- dimensional wave numbers. It provides insights into the stability and dissipative properties of the method for various CFL and Fourier numbers. The effect of the grid shape functions, number of particles per cell and their locations inside the grid cell are also analyzed. The EXAGOOP MPM solver (https://github.com/NREL/Exagoop.git) is developed at the National Renewable Energy Laboratory as a part of the NAWI UHPRO project and is based on the AMReX framework. A single-level, uniform cartesian grid is used as the background mesh, while the particle class in AMReX is used to manage the material point operations. Linear hat and B-splines are used as grid shape functions, while the time integration is performed using explicit Euler time integration. EXAGOOP is both CPU and GPU compatible and has been demonstrated to work well on multiple compute architectures. The performance of EXAGOOP on various computing architectures is presented along with its application to study compaction and embossing of high-pressure reverse osmosis membranes. The MPM solution accurately reproduces the membrane deformation. The deformed pore size and structure simulated using MPM also agree well with experimental SEM images.

material point method↗

Attitude control of a spinning flexible spacecraft

The dynamics of rotational motion of a spinning orbiting spacecraft consisting of two rigid bodies connected by a flexible joint and arbitrary number of flexible appendages (two of which are flexible massless booms having masses on their tips) is analyzed. Active attitude control is provided by momentum exchange devices (e.g. control moment gyroscopes) or a mass expulsion system. The linearized equations of motion describing the vehicle are presented, and a large scale digital simulation that has been developed at the Marshall Space Flight Center is presented. A simplified model of the geometrically complex vehicle is selected to make it analytically tractable. The simplified model consists of a single rigid core body with two attached flexible massless booms having tip masses. The states of the vehicle are defined as small perturbations about its steady-state spin. An analysis is performed to determine the domain of stability.

Seltzer, S. M.↗

Instabilities in line-driven stellar winds. I - Dependence on perturbation wavelength

An analytical study is presented of the stability of absorption line-driven flows, such as found in stellar winds, in the presence of small-amplitude disturbances. A generalized calculation of the perturbed direct extinction force is performed and the evolution of the perturbation into the nonlinear regime and the dynamical results of perturbation in the scattered radiation field are examined. An expression is derived for the wavenumber variation of the perturbed line force from the milieu of nonoverlapping lines which have a power-law distribution in opacity. A linear dispersion analysis is carried out to model the growth and propagation of radiative-acoustic waves in absorption line-driven flows, which are found unstable to perturbations that may be smaller than the Sobolev length. No damping mechanism was found that would eliminate the absorption line-driven flows.

Owocki, S. P.↗

Vibration analysis of rotor systems using reduced subsystem models

A general impedance method using reduced submodels has been developed for the linear dynamic analysis of rotor systems. Formulated in terms of either modal or physical coordinates of the subsystems, the method enables imbalance responses at specific locations of the rotor systems to be efficiently determined from a small number of 'master' degrees of freedom. To demonstrate the capability of this impedance approach, the Space Shuttle Main Engine high-pressure oxygen turbopump has been investigated to determine the bearing loads due to imbalance. Based on the same formulation, an eigenvalue analysis has been performed to study the system stability. A small 5-DOF model has been utilized to illustrate the application of the method to eigenvalue analysis. Because of its inherent characteristics of allowing formulation of reduced submodels, the impedance method can significantly increase the computational speed.

Fan, Uei-Jiun↗

Analysis and Tests of Reinforced Carbon-Epoxy/Foam-Core Sandwich Panels with Cutouts

The results of a study of a low-cost structurally efficient minimum-gage shear-panel design that can be used in light helicopters are presented. The shear-panel design is based on an integrally stiffened syntactic-foam stabilized-skin with an all-bias-ply tape construction for stabilized-skin concept with an all-bias-ply tape construction for the skins. This sandwich concept is an economical way to increase the panel bending stiffness weight penalty. The panels considered in the study were designed to be buckling resistant up to 100 lbs/in. of shear load and to have an ultimate strength of 300 lbs/in. The panel concept uses unidirectional carbon-epoxy tape on a syntactic adhesive as a stiffener that is co-cured with the skin and is an effective concept for improving panel buckling strength. The panel concept also uses pultruded carbon-epoxy rods embedded in a syntactic adhesive and over-wrapped with a bias-ply carbon-epoxy tape to form a reinforcing beam which is an effective method for redistributing load around rectangular cutout. The buckling strength of the reinforced panels is 83 to 90 percent of the predicted buckling strength based on a linear buckling analysis. The maximum experimental deflection exceeds the maximum deflection predicted by a nonlinear analysis by approximately one panel thickness. The failure strength of the reinforced panels was two and a half to seven times of the buckling strength. This efficient shear-panel design concept exceeds the required ultimate strength requirement of 300 lbs/in by more than 100 percent.

Baker, Donald J.↗

Digital Real-Time Simulation and Power Quality Analysis of a Hydrogen-Generating Nuclear-Renewable Integrated Energy System

This paper investigates the challenges and solutions associated with integrating a hydrogen-generating nuclear-renewable integrated energy system (NR-IES) under a transactive energy framework. The proposed system directs excess nuclear power to hydrogen production during periods of low grid demand while utilizing renewables to maintain grid stability. Using digital real-time simulation (DRTS) in the Typhoon HIL 404 model, the dynamic interactions between nuclear power plants, electrolyzers, and power grids are analyzed to mitigate issues such as harmonic distortion, power quality degradation, and low power factor caused by large non-linear loads. A three-phase power conversion system is modeled using the Typhoon HIL 404 model and includes a generator, a variable load, an electrolyzer, and power filters. Active harmonic filters (AHFs) and hybrid active power filters (HAPFs) are implemented to address harmonic mitigation and reactive power compensation. The results reveal that the HAPF topology effectively balances cost efficiency and performance and significantly reduces active filter current requirements compared to AHF-only systems. During maximum electrolyzer operation at 4 MW, the grid frequency dropped below 59.3 Hz without filtering; however, the implementation of power filters successfully restored the frequency to 59.9 Hz, demonstrating its effectiveness in maintaining grid stability. Future work will focus on integrating a deep reinforcement learning (DRL) framework with real-time simulation and optimizing real-time power dispatch, thus enabling a scalable, efficient NR-IES for sustainable energy markets.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗