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,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Stability of spherical converging shock wave

Based on Guderley's self-similar solution, stability of spherical converging shock wave is studied. A rigorous linear perturbation theory is developed, in which the growth rate of perturbation is given as a function of the spherical harmonic number ℓ and the specific heats ratio γ. Numerical calculation reveals the existence of a γ-dependent cut-off mode number ℓ{sub c}, such that all the eigenmode perturbations for ℓ > ℓ{sub c} are smeared out as the shock wave converges at the center. The analysis is applied to partially spherical geometries to give significant implication for different ignition schemes of inertial confinement fusion. Two-dimensional hydrodynamic simulations are performed to verify the theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Implementation of a partitioned algorithm for simulation of large CSI problems

The implementation of a partitioned numerical algorithm for determining the dynamic response of coupled structure/controller/estimator finite-dimensional systems is reviewed. The partitioned approach leads to a set of coupled first and second-order linear differential equations which are numerically integrated with extrapolation and implicit step methods. The present software implementation, ACSIS, utilizes parallel processing techniques at various levels to optimize performance on a shared-memory concurrent/vector processing system. A general procedure for the design of controller and filter gains is also implemented, which utilizes the vibration characteristics of the structure to be solved. Also presented are: example problems; a user's guide to the software; the procedures and algorithm scripts; a stability analysis for the algorithm; and the source code for the parallel implementation.

Alvin, Kenneth F.↗
NASA NTRS1984

Orbital stability in combined uniform axial and three-dimensional wiggler magnetic fields for free-electron lasers

Zachary Phys. Rev. A 29 (6), 3224 (1984) recently analyzed the instability of relativistic-electron helical trajectories in combined uniform axial and helical wiggler magnetic fields when the radial variation of the wiggler field is taken into account. It is shown here that the type 2 instability comprised of secular terms growing linearly in time, identified by Zachary and earlier by Diament Phys. Rev. A 23 (5), 2537 (1981), is an artifact of simple perturbation theory. A multiple-time-scale perturbation analysis reveals a nonsecular evolution on a slower time scale which accommodates an arbitrary initial perturbation. It is shown that, in the absence of exponential instability, the electron seeks a modified helical orbit more appropriate to its perturbed state and oscillates stably about it. Thus, the perturbed motion is oscillatory but nonsecular, and hence the helical orbits are stable.

Johnston, S.↗
DOE OSTI2019

Analysis of the Detonation Wave Structure in a Linearized Rotating Detonation Engine

The rotating detonation engine is increasingly favored as a viable pressure gain combustion technology for both propulsion and power generation applications. Practical designs involve the discrete injection of fuel and air, which then partially mix to produce the reactive mixture that is processed by a continuously moving detonation wave within the detonation chamber. One of the fundamental challenges in the successful operation of rotating detonation engines is the design of a robust fuel-oxidizer injection system. To minimize pressure losses while ensuring efficient mixing that promotes stable detonation, the injection has to satisfy many different constraints. In this work, the impact of discrete injection on wave propagation is studied in an effort to understand the impact of fuel stratification on detonation stability. To this end, the structure of the detonation wave within a linearized rotating detonation engine with fully premixed hydrogen-air and hydrogen-oxygen injection at various equivalence ratios is analyzed. The direct numerical simulation of the linearized model detonation engine is compared with experimental results, and wave behavior is correlated to the wave structure and flow turbulence. As a result, the flow properties and chemical composition of the gases across the detonation wave are studied to examine the characteristics of the reaction zone.

42 ENGINEERING↗
DOE OSTI2016

Finite time step and spatial grid effects in δf simulation of warm plasmas

This paper introduces a technique for analyzing time integration methods used with the particle weight equations in δf method particle-in-cell (PIC) schemes. The analysis applies to the simulation of warm, uniform, periodic or infinite plasmas in the linear regime and considers the collective behavior similar to the analysis performed by Langdon for full-f PIC schemes [1,2]. We perform both a time integration analysis and spatial grid analysis for a kinetic ion, adiabatic electron model of ion acoustic waves. An implicit time integration scheme is studied in detail for δf simulations using our weight equation analysis and for full-f simulations using the method of Langdon. It is found that the δf method exhibits a CFL-like stability condition for low temperature ions, which is independent of the parameter characterizing the implicitness of the scheme. The accuracy of the real frequency and damping rate due to the discrete time and spatial schemes is also derived using a perturbative method. The theoretical analysis of numerical error presented here may be useful for the verification of simulations and for providing intuition for the design of new implicit time integration schemes for the δf method, as well as understanding differences between δf and full-f approaches to plasma simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Fragmentation of faceted crystalline wires

The reduction of surface energy is a well-known driving force for the fragmentation of nanowires. Fragmentation can limit the service life of those small low-dimensional building blocks or can be controlled to induce beneficial shape changes. For isotropic surface energy, fragmentation is analogous to the classic Rayleigh-Plateau instability of liquid jets. However, commonly synthesized crystalline nanowires have strongly anisotropic surface energies and exhibit facets coinciding with cusps in the Wulff plot. Depending on growth orientation, different fragmentation behaviors have been seen in nanowires. Using phase-field simulations, we show that fragmentation of faceted nanowires with cubic crystal symmetry may occur by a finite-amplitude nonlinear instability, as opposed to a Rayleigh-Plateau-like linear instability, depending on nanowire growth orientation. Here, we carry out a weakly nonlinear analysis based on sharp interface theory to characterize the faceted nanowire shape corresponding to a surface-energy saddle point. The analysis predicts that the minimum amplitude of a periodic shape perturbation to trigger fragmentation increases with cusp strength but decreases inversely proportionally to the perturbation wavelength for long wavelengths, in good quantitative agreement with phase-field predictions. The results provide the theoretical foundation to predict nanowire stability as a function of length and surface energy in diverse applications.

36 MATERIALS SCIENCE↗
NASA NTRS1992

Adaptive control of Space Station with control moment gyros

An adaptive approach to Space Station attitude control is investigated. The main components of the controller are the parameter identification scheme, the control gain calculation, and the control law. The control law is a full-state feedback space station baseline control law. The control gain calculation is based on linear-quadratic regulator theory with eigenvalues placement in a vertical strip. The parameter identification scheme is a recursive extended Kalman filter that estimates the inertias and also provides an estimate of the unmodeled disturbances due to the aerodynamic torques and to the nonlinear effects. An analysis of the inertia estimation problem suggests that it is possible to estimate Space Station inertias accurately during nominal control moment gyro operations. The closed-loop adaptive control law is shown to be capable of stabilizing the Space Station after large inertia changes. Results are presented for the pitch axis.

Bishop, Robert H.↗
NASA NTRS1995

Sensitivity analysis for aeroacoustic and aeroelastic design of turbomachinery blades

A new method for computing the effect that small changes in the airfoil shape and cascade geometry have on the aeroacoustic and aeroelastic behavior of turbomachinery cascades is presented. The nonlinear unsteady flow is assumed to be composed of a nonlinear steady flow plus a small perturbation unsteady flow that is harmonic in time. First, the full potential equation is used to describe the behavior of the nonlinear mean (steady) flow through a two-dimensional cascade. The small disturbance unsteady flow through the cascade is described by the linearized Euler equations. Using rapid distortion theory, the unsteady velocity is split into a rotational part that contains the vorticity and an irrotational part described by a scalar potential. The unsteady vorticity transport is described analytically in terms of the drift and stream functions computed from the steady flow. Hence, the solution of the linearized Euler equations may be reduced to a single inhomogeneous equation for the unsteady potential. The steady flow and small disturbance unsteady flow equations are discretized using bilinear quadrilateral isoparametric finite elements. The nonlinear mean flow solution and streamline computational grid are computed simultaneously using Newton iteration. At each step of the Newton iteration, LU decomposition is used to solve the resulting set of linear equations. The unsteady flow problem is linear, and is also solved using LU decomposition. Next, a sensitivity analysis is performed to determine the effect small changes in cascade and airfoil geometry have on the mean and unsteady flow fields. The sensitivity analysis makes use of the nominal steady and unsteady flow LU decompositions so that no additional matrices need to be factored. Hence, the present method is computationally very efficient. To demonstrate how the sensitivity analysis may be used to redesign cascades, a compressor is redesigned for improved aeroelastic stability and two different fan exit guide vanes are redesigned for reduced downstream radiated noise. In addition, a framework detailing how the two-dimensional version of the method may be used to redesign three-dimensional geometries is presented.

Lorence, Christopher B.↗
NASA NTRS1998

Supersonic Coaxial Jets: Noise Predictions and Measurements

The noise from perfectly expanded coaxial jets was measured in an anechoic chamber for different operating conditions with the same total thrust, mass flow, and exit area. The shape of the measured noise spectrum at different angles to the jet axis was found to agree with spectral shapes for single, axisymmetric jets. Based on these spectra, the sound was characterized as being generated by large turbulent structures or fine-scale turbulence. Modeling the large scale structures as instability waves, a stability analysis was conducted for the coaxial jets to identify the growing and decaying instability waves in each shear layer and predict their noise radiation pattern outside the jet. When compared to measured directivity, the analysis identified the region downstream of the outer potential core, where the two shear layers were merging, as the source of the peak radiated noise where instability waves, with their origin in the inner shear layer, reach their maximum amplitude. Numerical computations were also performed using a linearized Euler equation solver. Those results were compared to both the results from the instability wave analysis and to measured data.

Dahl, Milo D.↗
DOE OSTI2021

Operating Stresses and Their Effects on Degradation of LSM-Based SOFC Cathodes

The performance of solid oxides fuel cells (SOFCs) with four different Mn excess of lanthanum strontium manganite (LSM) -based cathodes were examined under different temperatures (1,000 °C, 900 °C), current densities (0, 380, and 760 mA cm -2 ), and cathode atmospheres (p O 2 )=0.1,0.15,0.21) for durations ranging from 58 to 1,008 h. Each yttria-stabilized zirconia (YSZ) electrolyte-supported “button” cell had a porous Ni/YSZ composite anode and a porous LSM/YSZ composite cathode. The cells’ output voltage versus time were recorded and electrochemical impedance spectroscopy (EIS) and linear sweep voltammetry (LSV) measurements were performed every 24 hours. The total area specific resistance (ASR) was calculated from these measurements. The values of ASR from both EIS and LSV were comparable between each other but lower than that from durability test. Distribution of relaxation times (DRT) analysis was performed to investigate the electrochemical processes and their corresponding relaxation frequencies as well as their attributions to the total ASR. The series resistance and parallel resistance that obtained from the equivalent circuit fit of the Nyquist plot were higher at low temperature and low p O 2 regardless of LSM compositions. The area of the peaks derived from DRT analysis increased with time during individual tests. The microstructures of the tested cells were examined using scanning electron microscopy and energy-dispersive x-ray spectroscopy (SEM/EDS). Manganese oxide particles were observed near the cathode-electrolyte interface after prolonged (1,008 h) testing in air of a cell with LSM of 11% manganese excess, and after a short test (58 h) under low oxygen (p O 2 = 0.10) in a cell with LSM of 5% manganese excess. However, the role of MnOx in cell degradation is still unknown.

08 HYDROGEN↗
DOE OSTI2016

One-dimensional mercury(II) halide coordination polymers of 3,6-bis(2-pyridyl)-1,2,4,5-tetrazine ligand: Synthesis, crystal structure, spectroscopic and DFT studies

Three new coordination polymers, [Hg(μ-bptz)X{sub 2}]{sub n} (X=Cl (1), Br (2)) and [Hg{sub 2}(μ-bptz)(μ-I){sub 2}I{sub 2}]{sub n} (3) (bptz=3,6-bis(2-pyridyl)-1,2,4,5-tetrazine) were synthesized. X-ray structural analysis indicated that compounds 1 and 2 are composed of one-dimensional (1D) linear chains while the compound 3 has 1D stair-stepped structure. The electronic band structure along with density of states (DOS) calculated by the DFT method indicates that compound 1 and 2 are direct band gap semiconductors; however, compound 3 is an indirect semiconductor. The linear optical properties of the compounds are also calculated by DFT method. According to the DFT calculations, the observed emission band of the compounds in solid state is due to electron transfer from an excited bptz-π* state (CBs) to the top of VBs. {sup 1}H NMR spectra of the compounds indicate that, in solution phase, the compounds don’t decompose completely. Thermal stability of the compounds is studied using TG, DTA methods. - Graphical abstract: Synthesis, crystal structure and emission spectra of [Hg(μ-bptz)X{sub 2}]{sub n} (X=Cl and Br) and [Hg{sub 2}(μ-bptz)(μ-I){sub 2}I{sub 2}]{sub n} are presented. The electronic band structure and linear optical properties of the compounds are calculated by the DFT method. - Highlights: • Three 1D Hg(II) halide coordination polymers with bptz ligand have been prepared. • The structures of the compounds are determined by single crystal XRD. • DFT calculations show that [Hg(μ-bptz)X{sub 2}]{sub n} (X=Cl and Br) have a direct band gap. • DFT calculations show that [Hg{sub 2}(μ-bptz)(μ-I){sub 2}I{sub 2}]{sub n} has an indirect band gap. • The compounds show an intraligand electron transfer emission band in solid state.

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

Effect of resonant magnetic perturbations on microturbulence in DIII-D pedestal

Vacuum resonant magnetic perturbations (RMP) applied to otherwise axisymmetric tokamak plasmas produce in general a combination of non-resonant effects that preserve closed flux surfaces (kink response) and resonant effects that introduce magnetic islands and/or stochasticity (tearing response). The effect of the plasma kink response on the linear stability and nonlinear transport of edge turbulence is studied using the gyrokinetic toroidal code GTC for a DIII-D plasma with applied n = 2 vacuum RMP. GTC simulations use the 3D equilibrium of DIII-D discharge 158103 (Nazikian et al 2015 Phys. Rev. Lett. 114 105002), which is provided by nonlinear ideal MHD VMEC equilibrium solver in order to include the effect of the plasma kink response to the external field but to exclude island formation at rational surfaces. Analysis using the GTC simulation results reveal no increase of growth rates for the electrostatic drift wave instability and for the electromagnetic kinetic-ballooning mode in the presence of the plasma kink response to the RMP. Moreover, nonlinear electrostatic simulations show that the effect of the 3D equilibrium on zonal flow damping is very weak and found to be insufficient to modify turbulent transport in the electrostatic turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1987

Investigation of temperature and concentration oscillations in the directional solidification of Pb-Sn-Te

Directional solidification of the pseudobinary compound semiconductor material Pb sub 1-x Sn sub x Te by the Bridgman crystal growth process will be studied. Natural convection in the molten sample will be visualized with a novel electrochemical cell technique that employs the solid electrolyte material yttria-stabilized zirconia. Mass transfer by both diffusion and convection will be measured by detecting the motion of oxygen tracer in the liquid. Additional applications for electrochemical cells in semiconductor crystal growth are suggested. Unsteady convection in the melt will also be detected by the appearance of temperature oscillations. The purpose of this study is to experimentally characterize the overstable conditions for a Pb sub 1-x Sn sub x Te melt in the vertical Bridgman crystal growth technique and use a linear analysis to predict the onset of convection for this system.

Anderson, T. J.↗
NASA NTRS1988

Guaranteed properties for nonlinear gain scheduled control systems

An analysis is presented of two types of nonlinear gain-scheduled control systems: scheduling on a reference trajectory and scheduling on the plant output. Conditions are given which guarantee stability, robustness, and performance properties of the global gain-scheduled designs. These conditions confirm and formalize popular notions regarding gain-scheduled designs; e.g., that the scheduling variable should vary slowly and capture the plant's nonlinearities. These results extend the work of Shamma and Athans (1987), which addressed the case of linear plants whose dynamics depends on exogenous parameters.

Shamma, Jeff S.↗
NASA NTRS1957

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↗
DOE OSTI2016

Transition to turbulence and noise radiation in heated coaxial jet flows

Laminar-turbulent transition and noise radiation of a parametrized set of subsonic coaxial jet flows with a hot primary (core) stream are investigated numerically by Large-Eddy Simulation (LES) and direct noise computation. This study extends our previous research on local linear stability of heated coaxial jet flows by analyzing the nonlinear evolution of initially laminar flows disturbed by a superposition of small-amplitude unstable eigenmodes. First, a baseline configuration is studied to shed light on the flow dynamics of coaxial jet flows. Subsequently, LESs are performed for a range of Mach and Reynolds numbers to systematically analyze the influences of the temperature and the velocity ratios between the primary and the secondary (bypass) stream. The results provide a basis for a detailed analysis of fundamental flow-acoustic phenomena in the considered heated coaxial jet flows. Increasing the primary-jet temperature leads to an increase of fluctuation levels and to an amplification of far-field noise, especially at low frequencies. Strong mixing between the cold bypass stream and the hot primary stream as well as the intermittent character of the flow field at the end of the potential core lead to a pronounced noise radiation at an aft angle of approximately 35{sup ∘}. The velocity ratio strongly affects the shear-layer development and therefore also the noise generation mechanisms. Increasing the secondary-stream velocity amplifies the dominance of outer shear-layer perturbations while the disturbance growth rates in the inner shear layer decrease. Already for r{sub mic} > 40R{sub 1}, where r{sub mic} is the distance from the end of the potential core and R{sub 1} is the core-jet radius, a perfect 1/r{sub mic} decay of the sound pressure amplitudes is observed. The potential-core length increases for higher secondary-stream velocities which leads to a shift of the center of the dominant acoustic radiation in the downstream direction.

42 ENGINEERING↗
NASA NTRS1978

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.↗
NASA NTRS1991

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.↗