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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Linear beam stability in periodic focusing systems: Krein signature and band structure

The general question of how a beam becomes unstable has been one of the fundamental research topics among beam and accelerator physicists for several decades. In this study, we revisited the general problem of linear beam stability in periodic focusing systems by applying the concepts of Krein signature and band structure. We numerically calculated the eigenvalues and other associated characteristics of one-period maps, and discussed the stability properties of single-particle motions with skew quadrupoles and envelope perturbations in high-intensity beams on an equal footing. In particular, an application of the Krein theory to envelope instability analysis was newly attempted in this study. Finally, the appearance of instabilities is interpreted as the result of the collision between eigenmodes of opposite Krein signatures and the formation of a band gap.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Computation of linear MHD instabilities with the multi-region relaxed MHD energy principle

Over the last decade, a variational principle based on a generalisation of Taylor's relaxation, referred to as multi-region relaxed magnetohydrodynamics (MRxMHDs) has been developed. The numerical solutions of the MRxMHD equilibria have been constructed using the Stepped Pressure Equilibrium Code (SPEC) (Hudson et al 2012 Phys. Plasmas 19 112502). In principle, SPEC could also be established to describe the MRxMHD stability of an equilibrium. In this work, a theoretical framework is developed to relate the second variation of the energy functional to the so-called Hessian matrix, enabling the prediction of MHD linear instabilities of cylindrical plasmas, and is implemented in SPEC. The negative and positive eigenvalues of the Hessian matrix predict the stability of an equilibrium. Verification studies of SPEC stability results with the M3D-C 1 code and the tearing mode $\Delta^{^{\prime}}$ criterion have been conducted for ideal and resistive MHD instabilities, respectively, in a pressureless cylindrical tokamak, and have shown good agreement. Overall, our stability analysis is a critical step towards understanding the MHD stability of three-dimensional MHDs where nested flux surfaces, magnetic islands and stochastic regions co-exist.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1985

The utilization of satellite data and dynamics in understanding and predicting global weather phenomena

A two layer spectral quasi-geostrophic model is used to simulate the effects of topography on the equilibria, the stability, and the long term evaluation of incipient unstable waves. The flow is forced by latitudinally dependent radiational heating. The nature of the form drag instability of high index equilibria is investigated. The proximity of the equilibrium shear to a resonant value is essential for the instability, provided the equilibrium occurs at a slightly stronger shear than resonance. The properties of the steady Hadley and Rossby required for a thermally forced rotating fluid on a sphere are further explained. An objective parameterization technique is developed for general nonlinear hydrodynamical systems. The typical structure is one in which the rates of change of the dependent variables depend on homogeneous quadratic and linear forms, as well as on inhomogeneous forcing terms. Also documented is a steady, axisymmetric model of the general circulation developed as a basis for climate stability studies. The model includes the effects of heating, rotation, and internal friction, but neglects topography. Included is further research on cloud street phenomena. Orientation angles and horizontal wavelengths of boundary layer rolls and cloud streets are determined from an analysis of a truncated spectral model of three dimensional shallow moist Boussinesq convection in a shearing environment is further explained. Relatively broadly spaced roll clouds have orientations for which the Fourier component of the roll perpendicular shear is nearly zero, but the second corresponds to narrowly spaced rolls having orientations for which the Fourier coefficients of both the perpendicular and the parallel components of the shear are nearly equal.

Shirer, H. N.↗
NASA NTRS1988

Adaptive stabilization of discrete-time systems using linear periodically time varying controllers

A direct adaptive scheme based on the use of linear time-varying periodic controllers is proposed which estimates online the periodic coefficients of the controller. It is shown that adaptive stabilization is attained for all possibly nonstably invertible plants of known order but unknown delay. Although no appeal is made to persistency of excitation arguments, a provision is needed to avoid the singularity of an estimated matrix, this property being required only for the analysis and not the control calculations.

Ortega, Romeo↗
NASA NTRS1995

Triggering of longitudinal combustion instabilities in solid rocket motors: Nonlinear combustion response

Pulsed oscillations in solid rocket motors are investigated with emphasis on nonlinear combustion response. The study employs a wave equation governing the unsteady motions in a two-phase flow, and a solution technique based on spatial- and time-averaging. A wide class of combustion response functions is studied to second-order in fluctuation amplitude to determine if, when, and how triggered instabilities arise. Conditions for triggering are derived from analysis of limit cycles, and regions of triggering are found in parametric space. Based on the behavior of model dynamical systems, introduction of linear cross-coupling and quadratic self-coupling among the acoustic modes appears to be the manner in which the nonlinear combustion response produces triggering to a stable limit cycle. Regions of initial conditions corresponding to stable pulses were found, suggesting that stability depends on initial phase angle and harmonic content, as well as the composite amplitude, of the pulse.

Wicker, J. M.↗
DOE OSTI2018

Deploy production sliding mesh capability with linear solver benchmarking.

Wind applications require the ability to simulate rotating blades. To support this use-case, a novel design-order sliding mesh algorithm has been developed and deployed. The hybrid method combines the control volume finite element methodology (CVFEM) with concepts found within a discontinuous Galerkin (DG) finite element method (FEM) to manage a sliding mesh. The method has been demonstrated to be design-order for the tested polynomial basis (P=1 and P=2) and has been deployed to provide production simulation capability for a Vestas V27 (225 kW) wind turbine. Other stationary and canonical rotating ow simulations are also presented. As the majority of wind-energy applications are driving extensive usage of hybrid meshes, a foundational study that outlines near-wall numerical behavior for a variety of element topologies is presented. Results indicate that the proposed nonlinear stabilization operator (NSO) is an effective stabilization methodology to control Gibbs phenomena at large cell Peclet numbers. The study also provides practical mesh resolution guidelines for future analysis efforts. Application-driven performance and algorithmic improvements have been carried out to increase robustness of the scheme on hybrid production wind energy meshes. Specifically, the Kokkos-based Nalu Kernel construct outlined in the FY17/Q4 ExaWind milestone has been transitioned to the hybrid mesh regime. This code base is exercised within a full V27 production run. Simulation timings for parallel search and custom ghosting are presented. As the low-Mach application space requires implicit matrix solves, the cost of matrix reinitialization has been evaluated on a variety of production meshes. Results indicate that at low element counts, i.e., fewer than 100 million elements, matrix graph initialization and preconditioner setup times are small. However, as mesh sizes increase, e.g., 500 million elements, simulation time associated with \setup-up" costs can increase to nearly 50% of overall simulation time when using the full Tpetra solver stack and nearly 35% when using a mixed Tpetra- Hypre-based solver stack. The report also highlights the project achievement of surpassing the 1 billion element mesh scale for a production V27 hybrid mesh. A detailed timing breakdown is presented that again suggests work to be done in the setup events associated with the linear system. In order to mitigate these initialization costs, several application paths have been explored, all of which are designed to reduce the frequency of matrix reinitialization. Methods such as removing Jacobian entries on the dynamic matrix columns (in concert with increased inner equation iterations), and lagging of Jacobian entries have reduced setup times at the cost of numerical stability. Artificially increasing, or bloating, the matrix stencil to ensure that full Jacobians are included is developed with results suggesting that this methodology is useful in decreasing reinitialization events without loss of matrix contributions. With the above foundational advances in computational capability, the project is well positioned to begin scientific inquiry on a variety of wind-farm physics such as turbine/turbine wake interactions.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1991

A globally convergent matrix-free algorithm for implicit time-marching schemes arising in finite element analysis in fluids

A solution procedure for solving nonlinear time-marching problems is presented. The nonsymmetric systems of equations arising from a Newton-type linearization of these time-marching problems are solved using an iterative strategy based on the generalized minimal residual (GMRES) algorithm. Matrix-free techniques leading to reduction in storage are presented. Incorporation of a linesearch algorithm in the Newton-GMRES scheme is discussed. An automatic time-increment control strategy is developed to increase the stability of the time-marching process. High-speed flow computations demonstrate the effectiveness of these algorithms.

Johan, Zdenek↗
NASA NTRS1974

Flutter of elastically supported orthotropic panels including the effects of flow angle

A theoretical panel flutter analysis and computer program was developed which is capable of analyzing an orthotropic panel with in-plane loads, at various angles of cross flow and various edge-support conditions. The resulting modal analysis uses linear piston-theory aerodynamics and includes both aerodynamic and structural damping. Calculations made for typical panels with no in-plane forces show that large reductions in dynamic presure for flutter are possible with only small changes in flow angle. The reduction in dynamic pressure with flow angle is greater for elastically supported panels than for simply supported panels. Aerodynamic damping has a significant stabilizing effect at all flow angles except zero, whereas structural damping has negligible effect.

Sawyer, J. W.↗
NASA NTRS1988

New class of control laws for robotic manipulators. I - Nonadaptive case. II - Adaptive case

A new class of exponentially stabilizing control laws for joint level control of robot arms is discussed. Closed-loop exponential stability has been demonstrated for both the set point and tracking control problems by a slight modification of the energy Lyapunov function and the use of a lemma which handles third-order terms in the Lyapunov function derivatives. In the second part, these control laws are adapted in a simple fashion to achieve asymptotically stable adaptive control. The analysis addresses the nonlinear dynamics directly without approximation, linearization, or ad hoc assumptions, and uses a parameterization based on physical (time-invariant) quantities.

Wen, John T.↗
NASA NTRS2004

Turbulent Jet Flames Into a Vitiated Coflow

Examined is the vitiated coflow flame, an experimental condition that decouples the combustion processes of flows found in practical combustors from the associated recirculating fluid mechanics. The configuration consists of a 4.57 mm diameter fuel jet into a coaxial flow of hot combustion products from a lean premixed flame. The 210 mm diameter coflow isolates the jet flame from the cool ambient, providing a hot environment similar to the operating conditions of advanced combustors; this important high temperature element is lacking in the traditional laboratory experiments of jet flames into cool (room) air. A family of flows of increasing complexity is presented: 1) nonreacting flow, 2) all hydrogen flame (fuel jet and premixed coflow), and 3) set of methane flames. This sequence of experiments provides a convenient ordering of validation data for combustion models. Laser Raman-Rayleigh-LIF diagnostics at the Turbulent Diffusion Flame laboratory of Sandia National Laboratories produced instantaneous multiscalar point measurements. These results attest to the attractive features of the vitiated coflow burner and the well-defined boundary conditions provided by the coflow. The coflow is uniform and steady, isolating the jet flame from the laboratory air for a downstream distance ranging from z/d = 50-70. The statistical results show that differential diffusion effects in this highly turbulent flow are negligible. Complementing the comprehensive set of multiscalar measurements is a parametric study of lifted methane flames that was conducted to analyze flame sensitivity to jet and coflow velocity, as well as coflow temperature. The linear relationship found between the lift-off height and the jet velocity is consistent with previous experiments. New linear sensitivities were found correlating the lift-off height to coflow velocity and temperature. A blow-off study revealed that the methane flame blows off at a common coflow temperature (1260 K), regardless of coflow or jet velocity. An explanation for this phenomenon is that entrainment of ambient air at the high lift-off heights prevents autoignition. Analysis of the results suggests that flame stabilization occurs through a combination of flame propagation, autoignition, and localized extinction processes. Proposed is an expanded view of distributed reaction combustion based on analysis of the distributions of probe volume conditions at the stabilization region of the lifted hydrogen and methane flames. Turbulent eddies the size of the flame thickness mix fuel and hot coflow across the flame front, thereby enhancing the reaction zone with autoignition of reactants at elevated temperatures; this is the reverse effect of turbulent flames in ambient air, where intense turbulence in cool mixtures result in localized extinction. Each of the three processes (i.e., flame propagation, autoignition and localized extinction) contributes to flame stabilization in varying degrees, depending on flow conditions.

Holdeman, James D.↗
NASA NTRS2015

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework [1, 2, 3]. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick [4, 5]. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, no slip boundary conditions, and the hybrid outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, thrust, ect.). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the Acoustic Velocity Potential Equation (AVPE). The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob [2, 3]. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, Sean R.↗
NASA NTRS1987

Dynamic stability and handling qualities tests on a highly augmented, statically unstable airplane

This paper describes some novel flight tests and analysis techniques in the flight dynamics and handling qualities area. These techniques were utilized during the initial flight envelope clearance of the X-29A aircraft and were largely responsible for the completion of the flight controls clearance program without any incidents or significant delays. The resulting open-loop and closed-loop frequency responses and the time history comparison using flight and linear simulation data are discussed.

Gera, Joseph↗
NASA NTRS1987

Stability and robustness of slowly time-varying linear systems

A well-known result for finite-dimensional time-varying linear systems is that if each `frozen time' is stable, then the time-varying system is stable for sufficiently slow time-variations. These results are reviewed and extended to a class of Volterra integrodifferential equations, specifically, differential equations with a convolution operator in the right-hand-side. The results are interpreted in the context of robustness of time-varying linear systems with special emphasis on analysis of gain-scheduled control systems.

Shamma, Jeff S.↗
DOE OSTI2024

Including the parallel mass flow in calculating the steady-state solutions and stability of the momentum balance equations for a quasisymmetric stellarator

The Helically Symmetric Experiment (HSX) is a quasisymmetric stellarator with minimal parallel viscous damping in a helical direction. The parallel flow (Vǁ) along the magnetic field is similarly weakly damped by viscosity. In this paper, the self-consistent steady-state parallel and poloidal momentum balance equations are used to show that a large Vǁ on the order of the ion thermal velocity can increase the ion resonant radial electric field (Er) beyond the value calculated using the typical approximation that Vǁ is zero. By altering the damping of Vǁ, either by degrading the quasisymmetry or varying the neutral density, the ion resonant Er can shift in a controllable fashion. It is shown explicitly that there exist stable and unstable steady-state solutions in the two-dimensional space of Vǁ and Er. A stability analysis of each solution is performed by calculating the eigenvalues and eigenvectors of the Jacobian. The unstable solution corresponds to a saddle point in which the eigenvalues have opposite signs. The analysis leads to the conclusion that unstable solutions occur when the derivative of the total poloidal damping with respect to Er is positive. A hysteresis in Er and Vǁ is observed when the radial current density is linearly increased to a maximum and then decreased back to zero. Jumps in the radial electric field and the parallel flow are observed as the radial current density drives the evolution from one stable point to the next. This result is similar to experimental data observed on several devices.

Physics↗
NASA NTRS1976

Applications Study of Aero-Maneuvering Orbit-to-Orbit Shuttle (AMOOS), Executive Summary

Preliminary designs and the supporting analysis for both the Aeromaneuvering Orbit-to-Orbit Shuttle (AMOOS) and the Aeromaneuvering Recovery System (AMRS) are presented. The AMOOS design is shown to yield from twice to almost three times the high energy, round-trip payloads as a purely propulsive vehicle of the same weight. Typically AMOOS can perform a crew rotation mission to equatorial geosynchronous orbit in one space shuttle launch. The flight test program analysis modeled parameters such as maximum dynamic pressure, heating rates, guidance, stability and recovery. Two model flight schedules were developed, one consisting of four flights and the other of two flights. The former is considered a very low risk, high information return program whereas the latter is a minimal cost program consistent with reasonable data returns and chance of success. The AMOOS and AMRS guidance scheme developed using linear regulator theory proved a precise and accurate guidance scheme. Both it and a classical linear system's based scheme were evaluated using 65 simulated trajectories in which the position in the entry corridor and the atmospheric density were varied randomly.

White, J.↗
NASA NTRS1989

A genuinely multi-dimensional upwind cell-vertex scheme for the Euler equations

The solution of the two-dimensional Euler equations is based on the two-dimensional linear convection equation and the Euler-equation decomposition developed by Hirsch et al. The scheme is genuinely two-dimensional. At each iteration, the data are locally decomposed into four variables, allowing convection in appropriate directions. This is done via a cell-vertex scheme with a downwind-weighted distribution step. The scheme is conservative, and third-order accurate in space. The derivation and stability analysis of the scheme for the convection equation, and the derivation of the extension to the Euler equations are given. Preconditioning techniques based on local values of the convection speeds are discussed. The scheme for the Euler equations is applied to two channel-flow problems. It is shown to converge rapidly to a solution that agrees well with that of a third-order upwind solver.

Powell, Kenneth G.↗
NASA NTRS1992

Stability of a class of interconnected evolution systems

Stability conditions for a class of interconnected systems modeled by linear abstract evolution equations and a memoryless nonlinearity are derived. These conditions are stated in terms of the passivity of each of the subsystems and can be considered as a partial generalization of the hyperstability theorem. A Liapunov function approach is used in the proof without requiring the positive definiteness of the Liapunov function. Application to the robustness analysis of the infinite-dimensional linear quadratic regulator is also discussed.

Wen, John T.↗
NASA NTRS1977

Unsteady linearized transonic flow analysis for slender bodies

An unsteady linearized formulation based on Oswatitsch-Keune's parabolic method is developed to analyze transonic flow past oscillating slender bodies. In contrast to the widely used integral transform method, it is shown that all solutions can be derived by a simpler method directly in the physical plane. By various expansion procedures, low-frequency solutions then are derived according to two clearly defined frequency ranges. Adams-Sears' iteration is employed to account for the second-order effects. Stability derivatives are compared with available theories and data. It is found that the derivatives depend more sensitively on thickness than on the reduced frequency. Finally, a critical assessment of the present method is given.

Liu, D. D.↗