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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Fluid System Stability Analysis Techniques

Fluid systems, or networks, consist of multiple components that work together to achieve some desired thermofluid state. For any generic application, this fluid state can be a combination of the fluid pressure, flow rate, enthalpy, or species concentration. Fluid system components, such as pumps and valves, are often governed by nonlinear differential equations, resulting in complex component-to-component interactions. System-level fluid network stability occurs when the flow through the system can maintain a steady-state solution in the presence of small perturbations, which depends on these component interactions. System instability, however, can go undetected until issues arise during integrated system testing. This presentation explores a method for system designers to think of the fluid network as an assembly of components, each with their own thermofluid surfaces of partial stability, called nullclines. The intersections of all nullclines yields system-level solutions, called equilibrium points. When designers define operating points, they are tuning system parameters so that these equilibrium points move to the desired location in the thermofluid state plane. However, linearization theory shows us that the dynamic behavior around these points can be unstable. The local stability of these equilibrium points can be assessed analytically with eigen-analysis, or numerically by propagating state-plane samples to construct a phase portrait. Investigating a phase portrait can help designers gain a qualitative understanding of a system’s dynamic performance. This understanding can then help inform requirement definitions, component selection, and operational procedures. This presentation includes an example of the phase portrait technique on a system featuring a centrifugal pump and a back-pressure regulator (BPR). Numerical modeling of this system suggests that equilibrium points on the left-hand side of the pump curve are dynamically unstable.

Fluid System↗
NASA NTRS1972

Linear and nonlinear stability characteristics of whistlers

Linear and nonlinear propagating characteristics of right-hand polarized, slow electromagnetic, magnetoplasma waves (whistlers) are discussed in terms of stability and dispersion. An analysis of the stability of whistlers propagating at an angle to the static magnetic field is presented. A new mechanism is derived for the onset of stimulated emissions, and modulational instability for nonlinear whistlers are discussed.

Brinca, A. L.↗
NASA NTRS1994

Dynamics of Numerics and CFD

This lecture attempts to illustrate the basic ideas of how the recent advances in nonlinear dynamical systems theory (dynamics) can provide new insights into the understanding of numerical algorithms used in solving nonlinear differential equations (DEs). Examples will be given of the use of dynamics to explain unusual phenomena that occur in numerics. The inadequacy of the use of linearized analysis for the understanding of long time behavior of nonlinear problems will be illustrated, and the role of dynamics in studying the nonlinear stability, accuracy, convergence property and efficiency of using time- dependent approaches to obtaining steady-state numerical solutions in computational fluid dynamics (CFD) will briefly be explained.

Yee, H. C.↗
DOE OSTI2025

A cross-dimensional analysis of data-driven short-term load forecasting methods with large-scale smart meter data

Electricity load forecasting is essential to utility operation and power grid stability. A wide spectrum of data-driven methods, ranging from linear regression models to more recent deep learning models have been adopted to forecast electric load over the years. However, there still lacks a holistic evaluation of the applicability of conventional statistical and machine learning based algorithms with respect to different temporal and spatial scopes, computational requirements, and sensitivity of model-tuning. Enabled by a large-scale electricity load profile dataset of over 40,000 residential customers in a utility region, we conducted a cross-dimensional analysis of data-driven load forecasting methods. Three regression-based and seven deep learning algorithms with different model configurations were evaluated in terms of their overall and peak load prediction accuracy, and training burdens, across spatial aggregation levels ranging from the transformer, feeder, substation, to neighborhood. We found, first, the load forecasting accuracy is constrained by a predictability boundary, influenced by the forecasting horizon and spatial aggregation level. Specifically, RandomForest, XGBoost, TFT, TSMixer, and TiDE models achieved less than 10 % prediction error for up to 96-h ahead forecasting for district, substation, and feeder levels, while other models struggle at long-horizon predictions; Second, for winter and summer peak load dates, most models were able to predict the peak demand timing within ± 1 h, but the prediction percentage error varied by models, with TFT and TiDE models being the top performers; Third, models with similar prediction accuracy can differ in training burden by an order of magnitude. Therefore, choosing model configurations that balance prediction performance and computational resource is an important practical consideration for large-scale deployment of the machine learning based load forecasting. The outcome of this study can guide researchers and practitioners to choose the proper load forecasting algorithms based on their problem scope, required accuracy, and available resources. The predictability boundary can serve as a benchmark for electricity load forecasting problems with new algorithms and datasets.

Li, Han↗
DOE OSTI2019

Gyrokinetic simulations of m=0 mode in sheared flow Z-pinch

Axisymmetric stability properties of sheared flow Z-pinch plasmas are studied by making use of the gyrokinetic approximation in the long-wavelength limit. Numerical simulations are carried out with the high-order finite-volume code COntinuum Gyrokinetic Edge New Technology (COGENT) and are analyzed for the parameters characteristic of the FuZE experiment. Reduction of the linear growth rate with increasing shear is observed, and the results are elucidated by making use of a local dispersion relation analysis. In addition, COGENT simulations are compared with fully kinetic particle-in-cell simulations, and with an ideal magnetohydrodynamics (MHD) model. Good agreement between the gyrokinetic and fully kinetic results for the linear stability is found, with the gyrokinetic model requiring much less computational time due to its ability to step over particle gyroperiod. The ideal MHD model is found to be consistent with the kinetic models in the long-wavelength part of the spectra ($kρ_i$), while failing to adequately predict short-scale ($kρ_i$) stability. Here, $k$ is the axial wavelength vector and $ρ_i$ is the ion gyroradius.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2021

Boundary Layer Instabilities Over a Cone-Cylinder-Flare Model at Mach 6

Computations are performed to investigate the boundary-layer instabilities over a sharp cone-cylinder-flare model at zero degrees angle of attack. The model geometry and the flow conditions are selected to match the experiments conducted in the Boeing/AFOSR Mach 6 Quiet Tunnel (BAM6QT) at Purdue University. The geometry consists of a nominally sharp 5-degree half-angle cone, followed by a cylindrical segment and then a 10-degree flare. An axisymmetric separation bubble is generated as a result of the laminar shock/boundary-layer interaction in the cylinder-flare region. The comparison of the laminar flow solution and the schlieren images shows a remarkable agreement between the respective locations of both the boundarylayer edge and the reattachment shock. The predicted heat flux distribution is also in agreement with the measured values downstream of the reattachment location. The analysis of convective and global instabilities is performed for flare half angles equal to 8, 10, and 12 degrees and nosetip radii equal to 0.1, 1, and 5 mm. The linear amplification of first and second Mack mode instabilities that begin to amplify in the cone region are computed with a combination of the parabolized stability equations (PSE) and the harmonic linearized Navier-Stokes equations (HLNSE). The predicted frequency spectra of the surface pressure fluctuations associated with both planar and oblique instability waves are compared with the measured spectra at the various locations of the PCB and Kulite sensors. The comparison shows that the computational analysis captures the distinct lobes within the disturbance amplification spectra measured in the experiments, but some differences in amplification characteristics are noted at low frequencies. Overall, the oblique disturbances are found to be more amplified than the planar disturbances. To our knowledge, this represents the first successful comparison between convective instability analysis and measured surface pressure fluctuations for a hypersonic configuration with a separation bubble. Finally, the global instability analysis shows that the laminar flow becomes supercritical for flare half angles larger than 8 degrees. The unstable global mode for the experimental configuration of a 10 degrees flare and a sharp nosetip cone corresponds to a stationary three-dimensional disturbance that is concentrated in the recirculation region and achieves its maximum growth rate for an azimuthal wavenumber of 5.

Boundary layer transition↗
DOE OSTI2020

The virtual element method for linear elastodynamics models: Design, analysis, and implementation

We design the conforming virtual element method for the numerical simulation of two dimensional time-dependent elastodynamics problems. We investigate the performance of the method both theoretically and numerically. We prove the stability and the convergence of the semi-discrete approximation in the energy norm and derive optimal error estimates. We also show the convergence in the L 2 norm. The performance of the virtual element method is assessed on a set of different computational meshes, including non-convex cells up to order four in the h-refinement setting. Exponential convergence is also experimentally seen in the p-refinement setting.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2001

Integrated Modeling for the Next Generation Space Telescope "Yardstick" Concept

The so-called NASA "Yardstick" design concept for the Next Generation Space Telescope presents unique challenges for systems-level analysis. Simulations that integrate controls, optics, thermal, and structural models are required to evaluate baseline performance, study design sensitivities, and perform design optimization. An integrated modeling approach was chosen using a combination of commercial off-the-shelf and "in-house" developed codes. The resulting capability provides a foundation for linear and non-linear analysis, using both the time and frequency-domain methods. It readily allows various combinations of design parameters and environmental loads to be evaluated directly in terms of key science-related metrics, in this case the scalar RMS (root mean square) line-of-sight and RMS wavefront errors. This presentation first addresses the development of the component, or discipline, models for the Yardstick design. It will then proceed to present the integration of the component models, using linear-systems approaches, in order to support two of the most critical baseline performance analyses: jitter and thermal-elastic stability of the optical telescope assembly (OTA). The results of the jitter analysis indicate that disturbances from the reaction wheels coupled with the lightly-damped and highly-flexible structure present significant challenges to the baseline line-of-sight control architecture. Vibration isolation will be required to meet jitter error requirements. The results of the thermal-elastic analysis indicate that the mirror segment displacements due to ground-to-orbit cool-down of the telescope are within the expected capture range of the segment rigid-body control actuators. This means we will be able to align and phase the primary mirror. However, the results for the analysis of the thermal transient response following an attitude maneuver (slew) show that this telescope design is not sufficiently stable, passively, to meet the wavefront error requirements. Structural re-design is one possibility; alternatively, active thermal control of the OTA may be considered. The Yardstick integrated models were successfully used to demonstrate the feasibility of two thermal control strategies.

Mosier, Gary E.↗
NASA NTRS1976

Nonlinear Equations of Motion for Cantilever Rotor Blades in Hover with Pitch Link Flexibility, Twist, Precone, Droop, Sweep, Torque Offset, and Blade Root Offset

Nonlinear equations of motion for a cantilever rotor blade are derived for the hovering flight condition. The blade is assumed to have twist, precone, droop, sweep, torque offset and blade root offset, and the elastic axis and the axes of center of mass, tension, and aerodynamic center coincident at the quarter chord. The blade is cantilevered in bending, but has a torsional root spring to simulate pitch link flexibility. Aerodynamic forces acting on the blade are derived from strip theory based on quasi-steady two-dimensional airfoil theory. The equations are hybrid, consisting of one integro-differential equation for root torsion and three integro-partial differential equations for flatwise and chordwise bending and elastic torsion. The equations are specialized for a uniform blade and reduced to nonlinear ordinary differential equations by Galerkin's method. They are linearized for small perturbation motions about the equilibrium operating condition. Modal analysis leads to formulation of a standard eigenvalue problem where the elements of the stability matrix depend on the solution of the equilibrium equations. Two different forms of the root torsion equation are derived that yield virtually identical numerical results. This provides a reasonable check for the accuracy of the equations.

Hodges, D. H.↗
NASA NTRS1978

Radiometric correction procedure study

A comparison of MSS radiometric processing techniques identified as a preferred radiometric processing technique a procedure which equalizes the mean and standard deviation of detector-specific histograms of uncalibrated scene data. Evaluation of MSS calibration data demonstrated that the relationship between detector responses is essentially linear over the range of intensities typically observed in MSS data, and that the calibration wedge data possess a high degree of temporal stability. An analysis of the preferred radiometric processing technique showed that it could be incorporated into the MDP-MSS system without a major redesign of the system, and with minimal impact on system throughput.

Colby, C.↗
NASA NTRS1989

Stability and dynamic analysis of a slender column with curved longitudinal stiffeners

The results of a stability design study are presented for a slender column with curved longitudinal stiffeners for large space structure applications. Linear stability analyses are performed using a link-plate representation of the stiffeners to determine stiffener local buckling stresses. Results from a set of parametric analyses are used to determine an approximate explicit expression for stiffener local buckling in terms of its geometric parameters. This expression along with other equations governing column stability and mass are assembled into a determinate system describing minimum mass stiffened column design. An iterative solution is determined to solve this system and a computer program incorporating this routine is presented. Example design problems are presented which verify the solution accuracy and illustrate the implementation of the solution routine. Also, observations are made which lead to a greatly simplified first iteration design equation relating the percent increase in column mass to the percent increase in column buckling load. From this, generalizations are drawn as to the mass savings offered by the stiffened column concept. Finally, the percent increase in fundamental column vibration frequency due to the addition of deployable stiffeners is studied.

Lake, Mark S.↗
NASA NTRS1991

Nonlinear modal resonances in low-gravity slosh-spacecraft systems

Nonlinear models of low gravity slosh, when coupled to spacecraft vibrations, predict intense nonlinear eigenfrequency shifts at zero gravity. These nonlinear frequency shifts are due to internal quadratic and cubic resonances between fluid slosh modes and spacecraft vibration modes. Their existence has been verified experimentally, and they cannot be correctly modeled by approximate, uncoupled nonlinear models, such as pendulum mechanical analogs. These predictions mean that linear slosh assumptions for spacecraft vibration models can be invalid, and may lead to degraded control system stability and performance. However, a complete nonlinear modal analysis will predict the correct dynamic behavior. This paper presents the analytical basis for these results, and discusses the effect of internal resonances on the nonlinear coupled response at zero gravity.

Peterson, Lee D.↗
NASA NTRS1995

Boundary layer transition studies

A small-scale wind tunnel previously used for turbulent boundary layer experiments was modified for two sets of boundary layer transition studies. The first study concerns a laminar separation/turbulent reattachment. The pressure gradient and unit Reynolds number are the same as the fully turbulent flow of Spalart and Watmuff. Without the trip wire, a laminar layer asymptotes to a Falkner & Skan similarity solution in the FPG. Application of the APG causes the layer to separate and a highly turbulent and approximately 2D mean flow reattachment occurs downstream. In an effort to gain some physical insight into the flow processes a small impulsive disturbance was introduced at the C(sub p) minimum. The facility is totally automated and phase-averaged data are measured on a point-by-point basis using unprecedently large grids. The evolution of the disturbance has been tracked all the way into the reattachment region and beyond into the fully turbulent boundary layer. At first, the amplitude decays exponentially with streamwise distance in the APG region, where the layer remains attached, i.e. the layer is viscously stable. After separation, the rate of decay slows, and a point of minimum amplitude is reached where the contours of the wave packet exhibit dispersive characteristics. From this point, exponential growth of the amplitude of the disturbance is observed in the detached shear layer, i.e. the dominant instability mechanism is inviscid. A group of large-scale 3D vortex loops emerges in the vicinity of the reattachment. Remarkably, the second loop retains its identify far downstream in the turbulent boundary layer. The results provide a level of detail usually associated with CFD. Substantial modifications were made to the facility for the second study concerning disturbances generated by Suction Holes for laminar flow Control (LFC). The test section incorporates suction through interchangeable porous test surfaces. Detailed studies have been made using isolated holes in the impervious test plate that used to establish the Blasius base flow. The suction is perturbed harmonically and data are averaged on the basis of the phase of the disturbance, for conditions corresponding to strong suction and without suction. The technique was enhanced by using up to nine multiple probes to reduce the experimental run-time. In both cases, 3D contour surfaces in the vicinity of the hole show highly 3D TS waves which fan out in the spanwise direction forming bow-shaped waves downstream. The case without suction has proved useful for evaluating calculation methods. With suction, the perturbations on the centerline are much stronger and decay less rapidly, while the TS waves in the far field are similar to the case without suction. Downstream, the contour surfaces of the TS waves develop spanwise irregularities which eventually form into clumps. The spanwise clumping is evidence of a secondary instability that could be associated with suction vortices. Designers of porous surfaces use Goldsmith's Criterion to minimize cross-stream interactions. It is shown that partial TS wave cancellation is possible, depending on the hole spacing, disturbance frequency and free-stream velocity. New high-performance Constant Temperature Hot-Wire Anemometers were designed and built, based on a linear system theory analysis that can be extended to arbitrary order. The motivation was to achieve the highest possible frequency reponse while ensuring overall system stability. The performance is equal to or superior to commercially available instruments at about 10% of the cost. Details, such as fabrication drawings and a parts list, have been published to enable the instrument to be construced by others.

Watmuff, Jonathan H.↗
NASA NTRS2004

Proceedings from the 2nd International Symposium on Formation Flying Missions and Technologies

Topics discussed include: The Stellar Imager (SI) "Vision Mission"; First Formation Flying Demonstration Mission Including on Flight Nulling; Formation Flying X-ray Telescope in L2 Orbit; SPECS: The Kilometer-baseline Far-IR Interferometer in NASA's Space Science Roadmap Presentation; A Tight Formation for Along-track SAR Interferometry; Realization of the Solar Power Satellite using the Formation Flying Solar Reflector; SIMBOL-X : Formation Flying for High-Energy Astrophysics; High Precision Optical Metrology for DARWIN; Close Formation Flight of Micro-Satellites for SAR Interferometry; Station-Keeping Requirements for Astronomical Imaging with Constellations of Free-Flying Collectors; Closed-Loop Control of Formation Flying Satellites; Formation Control for the MAXIM Mission; Precision Formation Keeping at L2 Using the Autonomous Formation Flying Sensor; Robust Control of Multiple Spacecraft Formation Flying; Virtual Rigid Body (VRB) Satellite Formation Control: Stable Mode-Switching and Cross-Coupling; Electromagnetic Formation Flight (EMFF) System Design, Mission Capabilities, and Testbed Development; Navigation Algorithms for Formation Flying Missions; Use of Formation Flying Small Satellites Incorporating OISL's in a Tandem Cluster Mission; Semimajor Axis Estimation Strategies; Relative Attitude Determination of Earth Orbiting Formations Using GPS Receivers; Analysis of Formation Flying in Eccentric Orbits Using Linearized Equations of Relative Motion; Conservative Analytical Collision Probabilities for Orbital Formation Flying; Equations of Motion and Stability of Two Spacecraft in Formation at the Earth/Moon Triangular Libration Points; Formations Near the Libration Points: Design Strategies Using Natural and Non-Natural Ares; An Overview of the Formation and Attitude Control System for the Terrestrial Planet Finder Formation Flying Interferometer; GVE-Based Dynamics and Control for Formation Flying Spacecraft; GNC System Design for a New Concept of X-Ray Distributed Telescope; GNC System for the Deployment and Fine Control of the DARWIN Free-Flying Interferometer; Formation Algorithm and Simulation Testbed; and PLATFORM: A Formation Flying, RvD and Robotic Validation Test-bench.

Source record↗
DOE OSTI2018

A Study of the von Neumann Stability Analysis of a Semi-implicit Numerical Method Applied Over the Radial Impurity Transport Equation in Tokamak Plasma

The radial impurity transport equation for tokamak plasma is a form of diffusion–convection–reaction equation. The impurity transport equation is solved to determine the distribution of impurity (non-fuel) ion species with different ionization states perpendicular to magnetic surfaces of tokamak plasma. The equation for each charge (ionization) state Z is a non-linear, second-order in space, first-order in time, parabolic partial differential equation coupled to the previous Z − 1 and the next Z + 1 charge states of the impurity species through its reaction term. The number of differential equations to be solved simultaneously is hence determined by the number of ionization states of the impurity species studied. The solution to the set of these coupled equations can be obtained using a semi-implicit numerical method applied on it. The present study describes the application of von Neumann stability analysis over the semi-implicit numerical method applied over the radial impurity transport equation and determines a generic stability criterion for the method. The stability analysis is further illustrated using the geometry of Aditya tokamak installed at the Institute for Plasma Research Gandhinagar, India as an example. The impurity species considered is oxygen (Atomic number = 8). This leads to a set of eight coupled equations for charge states Z = 1 to 8 over which von Neumann analysis is illustrated in present study.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1993

ISAC - A tool for aeroservoelastic modeling and analysis

This paper discusses the capabilities of the Interaction of Structures, Aerodynamics, and Controls (ISAC) system of program modules. The major modeling, analysis, and data management components of ISAC are identified. Equations of motion are displayed for a Laplace-domain representation of the unsteady aerodynamic forces. Options for approximating a frequency-domain representation of unsteady aerodynamic forces with rational functions of the Laplace variable are shown. Linear time invariant state-space equations of motion that result are discussed. Model generation and analyses of stability and dynamic response characteristics are shown for an aeroelastic vehicle which illustrate some of the capabilities of ISAC as a modeling and analysis tool for aeroelastic applications.

Adams, William M., Jr.↗