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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A Simple Comparison of Biochemical Systems Theory and Metabolic Control Analysis

This paper explores some basic concepts of Biochemical Systems Theory (BST) and Metabolic Control Analysis (MCA), two frameworks developed to understand the behavior of biochemical networks. Initially introduced by Savageau, BST focuses on system stability and employs power laws in modeling biochemical systems. On the other hand, MCA, pioneered by authors such as Kacser and Burns and Heinrich and Rapoport, emphasizes linearization of the governing equations and describes relationships (known as theorems) between different measures. Despite apparent differences, both frameworks are shown to be equivalent in many respects. Through a simple example of a linear chain, the paper demonstrates how BST and MCA yield identical results when analyzing steady-state behavior and logarithmic gains within biochemical pathways. This comparative analysis highlights the interchangeability of concepts such as kinetic orders, elasticities and other logarithmic gains.

FOS: Biological sciences↗
DOE OSTI2021

Modal dynamics and flutter analysis of floating offshore vertical axis wind turbines

Floating platforms provide vertical axis wind turbines (VAWTs) significant advantages over offshore horizontal axis wind turbines, which has led to increased research interest in offshore VAWT technology. Here, this paper examines dynamic stability of floating VAWTs and presents structural dynamic and aero-elastic flutter analysis of utility-scale 5 MW floating VAWTs by coupling a rotor finite element model with a linearized model for the floating system. The coupled floating VAWT model is utilized to study the effects of tower height, number of blades, blade tapering scheme, and impacts of the floating platform on modal properties and flutter. The study provides important new results and guidance to floating VAWT designers: (1) tower height increase leads to lower tower frequencies more prone to resonance, (2) by moving from land-based foundations to a floating TLP foundation, tower-mode frequency is increased significantly (by 49%–68%) which mitigates resonance and increases the flutter RPM (154% increase), (3) increasing the number of blades reduces the tower mode frequencies for both land-based and floating configuration, and (4) blade chord tapering can have a significant impact on increasing the tower mode flutter RPM. Further, analysis of low-frequency rigid body modes shows a safe, higher flutter RPM relative to the operating RPM.

17 WIND ENERGY↗
NASA NTRS1974

Steady-state decoupling and design of linear multivariable systems

A constructive criterion for decoupling the steady states of a linear time-invariant multivariable system is presented. This criterion consists of a set of inequalities which, when satisfied, will cause the steady states of a system to be decoupled. Stability analysis and a new design technique for such systems are given. A new and simple connection between single-loop and multivariable cases is found. These results are then applied to the compensation design for NASA STOL C-8A aircraft. Both steady-state decoupling and stability are justified through computer simulations.

Thaler, G. J.↗
NASA NTRS1976

Dynamics of flexible spinning satellites with radial wire antennas

A dynamic analysis is presented for a spin-stabilized spacecraft employing four radial wire antennas with tip masses. The satellite is modelled as a 14-degree-of-freedom system, and the linearized equations of motion are found. The lowest-order vibrational modes and natural frequencies of the gyroscopically coupled system are then determined. Because the satellite spin rate is decreased by antenna deployment, a spin-up maneuver is needed. The response of the time-varying mode equations during spin-up is found for the planar modes in terms of Bessel functions and a Struve function of order -0.25. The particular solution is expressed in various forms, including an infinite series of Bessel functions and a particularly useful asymptotic expansion. An error formula for the latter is derived, showing that it gives good accuracy. A simple approximation to the complementary function is obtained using the Wenzel-Kramers-Brillouin method, and the phase error in the approximation is shown to be small.

Longman, R. W.↗
NASA NTRS1976

Digital controllers for VTOL aircraft

Using linear-optimal estimation and control techniques, digital-adaptive control laws have been designed for a tandem-rotor helicopter which is equipped for fully automatic flight in terminal area operations. Two distinct discrete-time control laws are designed to interface with velocity-command and attitude-command guidance logic, and each incorporates proportional-integral compensation for non-zero-set-point regulation, as well as reduced-order Kalman filters for sensor blending and noise rejection. Adaptation to flight condition is achieved with a novel gain-scheduling method based on correlation and regression analysis. The linear-optimal design approach is found to be a valuable tool in the development of practical multivariable control laws for vehicles which evidence significant coupling and insufficient natural stability.

Stengel, R. F.↗
NASA NTRS1987

Effect of roughness on the stability of boundary layers

An analysis is conducted on the effect of imperfections consisting of humps and dips on the stability of incompressible flows over flat plates. The mean flow is calculated using interacting boundary layers. Linear quasiparallel spatial stability is used to calculate the growth rates and mode shapes of two-dimensional disturbances. Then, the amplification factor is computed. A search for the most dangerous frequency is conducted based on an amplification factor of 9 in the shortest distance. Correlations are made with the transition experiment of Walker and Greening using the e sup 9 method.

Nayfeh, Ali H.↗
NASA NTRS1992

Hypersonic vehicle control law development using H infinity and mu-synthesis

Applicability and effectiveness of robust control techniques to a single-stage-to-orbit (SSTO) airbreathing hypersonic vehicle on an ascent accelerating path and their effectiveness are explored in this paper. An SSTO control system design problem, requiring high accuracy tracking of velocity and altitude commands while limiting angle of attack oscillations, minimizing control power usage and stabilizing the vehicle all in the presence of atmospheric turbulence and uncertainty in the system, was formulated to compare results of the control designs using H infinity and mu-synthesis procedures. The math model, an integrated flight/propulsion dynamic model of a conical accelerator class vehicle, was linearized as the vehicle accelerated through Mach 8. Controller analysis was conducted using the singular value technique and the mu-analysis approach. Analysis results were obtained in both the frequency and the time domains. The results clearly demonstrate the inherent advantages of the structured singular value framework for this class of problems. Since payload performance margins are so critical for the SSTO mission, it is crucial that adequate stability margins be provided without sacrificing any payload mass.

Gregory, Irene M.↗
NASA NTRS1992

Hypersonic vehicle control law development using H infinity and mu-synthesis

Applicability and effectivness of robust control techniques to a single-stage-to-orbit (SSTO) airbreathing hypersonic vehicle on an ascent accelerating path and their effectivness are explored in this paper. An SSTO control system design problem, requiring high accuracy tracking of velocity and altitude commands while limiting an angle of attack oscillations, minimizing control power usage and stabilizing the vehicle all in the presence of atmospheric turbulence and uncertainty in the system, was formulated to compare results of the control designs using H infinity and mu-synthesis procedures. The math model, an integrated flight/propulsion dynamic model of a conical accelerator class vehicle, was linearized as the vehicle accelerated through Mach 8. Controller analysis was conducted using the singular value technique and the mu-analysis approach. Analysis results were obtained in both the frequency and the time domains. The results clearly demonstrate the inherent advantages of the structured singular value framework for this class of problems. Since payload performance margins are so critical for the SSTO mission, it is crucial that adequate stability margins be provided without sacrificing any payload mass.

Gregory, Irene M.↗
NASA NTRS1996

Development of Three-Dimensional Flow Code Package to Predict Performance and Stability of Aircraft with Leading Edge Ice Contamination

In the work reported herein, a simplified, uncoupled, zonal procedure is utilized to assess the capability of numerically simulating icing effects on a Boeing 727-200 aircraft. The computational approach combines potential flow plus boundary layer simulations by VSAERO for the un-iced aircraft forces and moments with Navier-Stokes simulations by NPARC for the incremental forces and moments due to iced components. These are compared with wind tunnel force and moment data, supplied by the Boeing Company, examining longitudinal flight characteristics. Grid refinement improved the local flow features over previously reported work with no appreciable difference in the incremental ice effect. The computed lift curve slope with and without empennage ice matches the experimental value to within 1%, and the zero lift angle agrees to within 0.2 of a degree. The computed slope of the un-iced and iced aircraft longitudinal stability curve is within about 2% of the test data. This work demonstrates the feasibility of a zonal method for the icing analysis of complete aircraft or isolated components within the linear angle of attack range. In fact, this zonal technique has allowed for the viscous analysis of a complete aircraft with ice which is currently not otherwise considered tractable.

Aircraft icing↗
DOE OSTI2018

Diagnostic and Power Feed Upgrades to the MAIZE Facility

The Michigan Accelerator for Inductive Z-pinch experiments (MAIZE) is a one stage, 1 MA, Linear Transformer Driver (LTD) at the University of Michigan. The MAIZE facility has been the site of various experimental studies, including wire array z-pinches for x-ray source development and cylindrical foil loads for implosion stability analysis. Here, in order to better understand and investigate these and other experiments, a 4- frame extreme ultraviolet (EUV) camera was implemented, the optical stability of a 12-frame laser shadowgraphy system was improved, B-dot probes were upgraded, and a Rogowski coil was constructed. New conical load hardware was also developed to reduce the inductance of the transmission line and thus improve the current delivered to the loads.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1990

Bearingless rotor aeromechanical stability measurements and correlations using nonlinear aerodynamics

The aeromechanical stability of a 1/8th Froude scaled bearingless rotor model was investigated experimentally in a wind tunnel. Both shaft-fixed and shaft-free conditions were examined to study the aeroelastic stability of a bearingless rotor without the incorporation of auxiliary dampers. This wind tunnel investigation generated a set of stability data for four different advance ratios, and a wide range of collective pitch settings. Theoretical analysis was performed using the newly developed University of Maryland Advanced Rotorcraft Code (UMARC). For analysis, the blade is modeled as an elastic beam undergoing flap bending, lag bending, elastic twist, and axial deformation. Blade response is calculated using a finite element method in time. Nonlinear aerodynamic effects are included by using a semiempirical stall modeling. The linearized periodic rotor perturbation equations in the nonrotating frame are solved for stability roots using Floquet transition matrix theory, as well as constant coefficient approximation. The predicted results are compared with the experimental data.

Wang, James M.↗
DOE OSTI2018

RF transient analysis and stabilization of the phase and energy of the proposed PIP-II LINAC

Here, this paper describes a recent effort to develop and benchmark a simulation tool for the analysis of RF transients and their compensation in an H- linear accelerator. Existing tools in this area either focus on electron LINACs or lack fundamental details about the LLRF system that are necessary to provide realistic performance estimates. In our paper we begin with a discussion of our computational models followed by benchmarking with existing beam-dynamics codes and measured data. We then analyze the effect of RF transients and their compensation in the PIP-II LINAC, followed by an analysis of calibration errors and how a Newton’s Method based feedback scheme can be used to regulate the beam energy to within the specified limits.

43 PARTICLE ACCELERATORS↗
NASA NTRS1997

Tailored Excitation for Multivariable Stability-Margin Measurement Applied to the X-31A Nonlinear Simulation

Safety and productivity of the initial flight test phase of a new vehicle have been enhanced by developing the ability to measure the stability margins of the combined control system and vehicle in flight. One shortcoming of performing this analysis is the long duration of the excitation signal required to provide results over a wide frequency range. For flight regimes such as high angle of attack or hypersonic flight, the ability to maintain flight condition for this time duration is difficult. Significantly reducing the required duration of the excitation input is possible by tailoring the input to excite only the frequency range where the lowest stability margin is expected. For a multiple-input/multiple-output system, the inputs can be simultaneously applied to the control effectors by creating each excitation input with a unique set of frequency components. Chirp-Z transformation algorithms can be used to match the analysis of the results to the specific frequencies used in the excitation input. This report discusses the application of a tailored excitation input to a high-fidelity X-31A linear model and nonlinear simulation. Depending on the frequency range, the results indicate the potential to significantly reduce the time required for stability measurement.

Bosworth, John T.↗
DOE OSTI2025

Enhanced MPM framework with multipatch isogeometric analysis for geotechnical applications

Achieving stable stress solutions at large strains using the Material Point Method (MPM) is challenging due to the accumulation of errors associated with geometry discretization, cell-crossing noise, and volumetric locking. Several simplified attempts exist in the literature to mitigate these errors, including higher-order frameworks. However, the stability of the MPM solution in such frameworks has been limited to simple geometries and the single-phase formulation (i.e., neglecting pore fluid). Although never explored, multipatch isogeometric analysis offers desirable qualities to simulate complex geometries while mitigating errors in the MPM. The degree of required high-order spatial integration has also never been investigated to infer a minimum limit for the stability of the stress solution in MPM. This paper presents a general-purpose numerical framework for simulating stable stresses in porous media, capturing both near incompressibility and multiphase interactions. First, the numerical framework is presented considering Non-Uniform Rational B-splines (NURBS) to perform isogeometric analysis (IGA) in MPM. Additionally, a volumetric strain smoothing algorithm is used to alleviate errors associated with volumetric locking. Second, the manifestation of cell-crossing errors is assessed via a series of problems with orders ranging from linear to cubic interpolation functions. Third, the use of NURBS is investigated and verified for problems with circular geometries. Finally, multipatch analysis is deployed to simulate plane strain and 3D penetration in soils, considering nearly incompressible elastoplastic (total stress) analysis and fully-coupled hydro-mechanical (effective stress) analysis. The stability of the solution is also analyzed for different constitutive models. From the results, it can be concluded that the framework using cubic interpolation functions with strain smoothing is the most convenient, presenting stable stress solutions for a broad range of multiphase geotechnical applications.

58 GEOSCIENCES↗
DOE OSTI2021

Implicit highly-coupled single-ion Hall-MHD formulation for hybrid particle-in-cell codes

The rudiments of a particle-based single-fluid two-temperature magnetohydrodynamic (MHD) algorithm have been outlined in Thoma et al. (2013). The extension of this algorithm to include the effect of Hall physics is described in this paper. An implicit leapfrog version of the algorithm, which allows timesteps large compared to the resistive decay time and other relevant timescales, has recently been added to a hybrid particle-in-cell code. In standard MHD the Hall term in the generalized Ohm’s law can often be neglected when the Hall parameter is small. This term must, however, be retained in regimes where it is non-negligible. The retention of displacement current in Maxwell’s equations avoids the numerical difficulties associated with the whistler mode, which are encountered in standard explicit Hall-MHD codes, and allows the algorithm to be incorporated into hybrid particle-in-cell codes, for which particles may migrate from a kinetic to fluid to MHD description based upon local ambient plasma conditions. A highly-coupled implicit Hall-MHD formalism is presented, in which displacement current can either be retained or neglected. Even when displacement current is neglected, the highly-coupled implicit formalism avoids the restrictive timesteps for the whistler mode in explicit Hall-MHD codes. A comparison of numerical and analytic dispersion analysis demonstrates the feasibility of this approach and establishes relevant constraints to assure numerical stability. The implementation of the algorithm is described, and test simulation results in 1D and 2D in both linear and nonlinear regimes are presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

A Crystalline Polyimide Porous Organic Framework for Selective Adsorption of Acetylene over Ethylene

The separation of acetylene from ethylene is a crucial process in the petrochemical industry, as even small acetylene impurities can lead to premature termination of ethylene polymerization. Herein, we present the synthesis of a robust, crystalline naphthalene diimide porous aromatic framework via imidization of linear naphthalene-1,4,5,8-tetracarboxylic dianhydride and triangular tris(4-aminophenyl)amine. The resulting material, PAF-110, exhibits impressive thermal and long-term structural stability, as indicated by thermogravimetric analysis and powder X-ray diffraction characterization. Gas adsorption characterization reveals that PAF-110 has a capacity for acetylene that is more than twice its ethylene capacity at 273 K and 1 bar, and it exhibits a moderate acetylene selectivity of 3.9 at 298 K and 1 bar. Complementary computational investigation of each guest binding in PAF-110 suggests that this affinity and selectivity for acetylene arises from its stronger electrostatic interaction with the carbonyl oxygen atoms of the framework. As a result, to the best of our knowledge, PAF-110 is the first crystalline porous organic material to exhibit selective adsorption of acetylene over ethylene, and its properties may provide insight into the further optimized design of porous organic materials for this key gas separation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

MHD stability trends and improved performance of LHD inward-shifted configurations: The role of the neutral beam current drive and thermal plasma density

The aim of the present study is to analyze the effect of the neutral beam current drive (NBCD), thermal plasma density, and NBI operational regime on the stability of pressure gradient-driven modes (PGDM) and Alfvén eigenmodes (AE) in LHD inward-shifted configurations. The stabilization of n/m = 1/2 PGDM (n toroidal mode and m poloidal mode) is observed in the discharge 167 800 during the co-NBCD phase. The iota profile evolution measured by motional stark effect diagnostic may indicate the iota profile up-shift caused by the co-NBCD can induce a non-resonant transition of the rational surface 1/2 before the mode stabilization. The evolution of the iota profile and continuum gaps in the discharge 167 805 during the ctr-NBCD phase leads to the stabilization of the AE, caused by the narrowing of the continuum gap as the iota profile down-shift. Opposite stability trends are identified for PGDM and AE stability with respect to the thermal plasma density. A larger thermal plasma density (larger thermal β) further enhances PGDM although the continuum gaps are narrower leading to configurations with stable AEs. The linear stability of AEs is analyzed using the gyro-fluid FAR3d code to reproduce the AE stability trends observed in the experiments with respect to the NBCD and thermal plasma density. The analysis of hypothetical scenarios dedicated to study different NBI operational regimes with respect to EP energy, and β and radial density profiles indicate off-axis NBI operation shows a higher EP β threshold to destabilize AEs compared to on-axis configuration. This is explained by the presence of a TAE gap in the inner plasma region, easily destabilized by an on-axis NBI injection. The control of the NBCD and thermal plasma in the discharge 167 800 shows a transitory stabilization of PGDM and AEs, as well as an improved discharge performance identified by an increment of the neutron fluxes.

Varela, Jacobo↗