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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS1980

Nonlinear grid error effects on numerical solution of partial differential equations

Finite difference solutions of nonlinear partial differential equations require discretizations and consequently grid errors are generated. These errors strongly affect stability and convergence properties of difference models. Previously such errors were analyzed by linearizing the difference equations for solutions. Properties of mappings of decadence were used to analyze nonlinear instabilities. Such an analysis is directly affected by initial/boundary conditions. An algorithm was developed, applied to nonlinear Burgers equations, and verified computationally. A preliminary test shows that Navier-Stokes equations may be treated similarly.

Dey, S. K.↗
NASA NTRS1984

Investigation of a nozzle instability on an F100 engine equipped with a digital electronic engine control

An instability in the nozzle of the F100 engine, equipped with a digital electronic engine control (DEEC), was observed during a flight evaluation on an F-15 aircraft. The instability occurred in the upper left hand corner (ULMC) of the flight envelope during augmentation. The instability was not predicted by stability analysis, closed-loop simulations of the the engine, or altitude testing of the engine. The instability caused stalls and augmentor blowouts. The nozzle instability and the altitude testing are described. Linear analysis and nonlinear digital simulation test results are presented. Software modifications on further flight test are discussed.

Burcham, F. W., Jr.↗
DOE OSTI2023

A General Autofluorescence Method to Characterize Polymerization Progress

Abstract An autofluorescence technique to characterize polymerization progress in real time/in line was developed, which functioned in the absence of typical fluorogenic groups on the monomer or polymer. The monomer dicyclopentadiene and polymer polydicyclopentadiene are hydrocarbons that lack traditional functional groups for fluorescence spectroscopy. Here, the autofluorescence of formulations containing this monomer and polymer during ruthenium‐catalyzed ring‐opening metathesis polymerization (ROMP) was harnessed for reaction monitoring. The methods fluorescence recovery after photobleaching (FRAP) and here‐developed fluorescence lifetime recovery after photobleaching (FLRAP) characterized polymerization progress in these native systems—without requiring exogenous fluorophore. (Auto)fluorescence lifetime recovery changes during polymerization correlated linearly to degree of cure, providing a quantitative link with reaction progress. These changing signals also provided relative rates of background polymerization, enabling comparison of 10 different catalyst‐inhibitor‐stabilized formulations. Multiple‐well analysis demonstrated suitability for future high‐throughput evaluation of formulations for thermosets. The central concept of the combined autofluorescence and FLRAP/FRAP method may be extendable to monitoring other polymerization reactions previously overlooked for lack of an obvious fluorescence handle.

López, Pía A.↗
DOE OSTI2025

Preserving nonlinear constraints in variational flow filtering data assimilation

Data assimilation aims to estimate the states of a dynamical system by optimally combining sparse and noisy observations of the physical system with uncertain forecasts produced by a computational model. The states of many dynamical systems of interest obey nonlinear physical constraints, and the corresponding dynamics is confined to a certain sub-manifold of the state space. Standard data assimilation techniques applied to such systems yield posterior states lying outside the manifold, violating the physical constraints. This work focuses on particle flow filters which use stochastic differential equations to evolve state samples from a prior distribution to samples from an observation-informed posterior distribution. The variational Fokker-Planck (VFP)—a generic particle flow filtering framework—is extended to incorporate non-linear, equality state constraints in the analysis. To this end, two algorithmic approaches that modify the VFP stochastic differential equation are discussed: (i) VFPSTAB, to inexactly preserve constraints with the addition of a stabilizing drift term, and (ii) VFPDAE, to exactly preserve constraints by treating the VFP dynamics as a stochastic differential-algebraic equation (SDAE). Additionally, an implicit-explicit time integrator is developed to evolve the VFPDAE dynamics. The strength of the proposed approach for constraint preservation in data assimilation is demonstrated on three test problems: the double pendulum, Korteweg-de-Vries, and the incompressible Navier-Stokes equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

High-speed X-ray imaging of the Leidenfrost collapse

Abstract The Leidenfrost layer is characterized by an insulating vapor film between a heated surface and an ambient liquid. The collapse of this film has been canonically theorized to occur from an interfacial instability between the liquid and vapor phases. The interfacial instability alone, however, is insufficient to explain the known influence of the surface on the film collapse process. In this work, we provide visual evidence for two key mechanisms governing the film collapse: the interfacial instability, and the nucleation of vapor upon multiple non-terminal liquid-solid contacts. These results were obtained by implementing high-speed X-ray imaging of the film collapse on a heated sphere submerged in liquid-water. The X-ray images were synchronized with a second high-speed visible light camera and two thermocouples to provide insight into the film formation and film collapse processes. Lastly, the dynamic film thickness was quantified by analysis of the X-ray images. This helped assess the influence of surface roughness on the disruption of the film. The results of this work encourage further investigation into non-linear stability theory to consolidate the role of the surface on the liquid-vapor interface during the film collapse process.

47 OTHER INSTRUMENTATION↗
DOE OSTI2018

Emerging investigator series: quantification of multiwall carbon nanotubes in plant tissues with spectroscopic analysis

If agricultural plants are exposed to carbon nanotubes (CNTs), they can potentially take up the CNTs from growth media and translocate them to their different tissues. In addition, agricultural application of CNTs recently attracted increasing attention, as they could promote germination, enhance crop yield, and exhibit other benefits. For evaluating the environmental effects of CNTs and optimizing their agricultural application, it is essential to quantify CNTs in plant tissues. In this study, pristine (p-) and carboxyl-functionalized (c-) multiwall CNTs (MWCNTs) were extracted from plant tissues by a sequential digestion with nitric acid (HNO 3 ) and sulfuric acid (H 2 SO 4 ). The extracted MWCNTs were stabilized with nonionic surfactant Triton X-100 and analyzed with ultraviolet-visible (UV-vis) spectroscopic analysis to measure the concentration of the MWCNTs in plant (lettuce) tissues. The MWCNT concentration was linearly correlated with the absorbance at 800 nm. The detection limit for p- and c-MWCNTs was achieved at 0.10–0.12, 0.070–0.081, 0.019–0.18 μg mg -1 for leaf, stem, and root tissues, respectively. The developed method was applied for lettuce (Lactuca sativa, cv. black seeded Simpson) hydroponically grown with 5, 10, 20 mg L -1 of p-MWCNTs and c-MWCNTs in the culture solution. We detected 0.21 ± 0.05–4.57 ± 0.39 μg mg -1 p-MWCNTs and 0.20 ± 0.17–0.75 ± 0.25 μg mg -1 c-MWCNTs in the lettuce roots, positively correlated with the dose of CNTs in solution. We have come up with a method for rapid quantification of CNTs in plant tissues using a widely-accessible technique, which can enable reliable analysis of CNTs in plant tissues and provide critical information for evaluating the environmental implications and managing agricultural application of CNTs.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1967

A modal combination computer program for dynamic analysis of structures

Computer program determines the response of a composite linear structure to sinusoidal base motion of a restrained structure or sinusoidal forces of a free structure. This program is applied to problems of testing practices and closed-loop stability of autopilot controlled space vehicles. It is written for the IBM 7094 in FORTRAN 4 language.

Bamford, R. M.↗
NASA NTRS1993

Effects of free-stream turbulence on boundary-layer transition

Traditional methods for the prediction of the transition location rest either on lump parameters like the momentum thickness or on N factors obtained from the approximate linear stability characteristics of the flow. These methods account for the effects of the disturbance environment by empirical correlation with observed transition locations. We have developed a highly efficient and accurate approach to stability analysis and transition simulation based on parabolized stability equations (PSE). Here we present a first application of this approach to realistic problems of practical interest. We study the effects of the turbulence level on transition and heat transfer in flows over flat plates, curved plates, and the stator blade of a gas-turbine model. We discuss the formulation of an input model for the PSE analysis in the light of the present knowledge of receptivity mechanisms. The input model is used to compute transition locations for various configurations and flow conditions. The results compare favorably with experimental data.

Herbert, TH.↗
NASA NTRS1952

Development of a laminar boundary layer behind a suction point

A theoretical investigation is made of the development of a laminar boundary layer behind a suction slot that is assumed to cut off part of the boundary layer without exerting any sink effect. The development, which is approximate, is based on the heat conduction equation. The heat conduction equation enters the analysis through a linearization of the Prandtl-Mises form of the boundary-layer equation.

EQUATIONS - BOUNDARY LAYER↗
NASA NTRS2012

Stability Analysis for HIFiRE Experiments

The HIFiRE-1 flight experiment provided a valuable database pertaining to boundary layer transition over a 7-degree half-angle, circular cone model from supersonic to hypersonic Mach numbers, and a range of Reynolds numbers and angles of attack. This paper reports selected findings from the ongoing computational analysis of the measured in-flight transition behavior. Transition during the ascent phase at nearly zero degree angle of attack is dominated by second mode instabilities except in the vicinity of the cone meridian where a roughness element was placed midway along the length of the cone. The growth of first mode instabilities is found to be weak at all trajectory points analyzed from the ascent phase. For times less than approximately 18.5 seconds into the flight, the peak amplification ratio for second mode disturbances is sufficiently small because of the lower Mach numbers at earlier times, so that the transition behavior inferred from the measurements is attributed to an unknown physical mechanism, potentially related to step discontinuities in surface height near the locations of a change in the surface material. Based on the time histories of temperature and/or heat flux at transducer locations within the aft portion of the cone, the onset of transition correlated with a linear N-factor, based on parabolized stability equations, of approximately 13.5. Due to the large angles of attack during the re-entry phase, crossflow instability may play a significant role in transition. Computations also indicate the presence of pronounced crossflow separation over a significant portion of the trajectory segment that is relevant to transition analysis. The transition behavior during this re-entry segment of HIFiRE-1 flight shares some common features with the predicted transition front along the elliptic cone shaped HIFiRE-5 flight article, which was designed to provide hypersonic transition data for a fully 3D geometric configuration. To compare and contrast the crossflow dominated transition over the HIFiRE-1 and HIFiRE-5 configurations, this paper also analyzes boundary layer instabilities over a subscale model of the HIFiRE-5 flight configuration that was tested in the Mach 6 quiet tunnel facility at Purdue University.

Li, Fei↗
NASA NTRS1999

Advanced Techniques for Assessment of Postural and Locomotor Ataxia, Spatial Orientation, and Gaze Stability

In addition to adapting to microgravity, major neurovestibular problems of space flight include postflight difficulties with standing, walking, turning corners, and other activities that require stable upright posture and gaze stability. These difficulties inhibit astronauts' ability to stand or escape from their vehicle during emergencies. The long-ter7n goal of the NSBRI is the development of countermeasures to ameliorate the effects of long duration space flight. These countermeasures must be tested with valid and reliable tools. This project aims to develop quantitative, parametric approaches for assessing gaze stability and spatial orientation during normal gait and when gait is perturbed. Two of this year's most important findings concern head fixation distance and ideal trajectory analysis. During a normal cycle of walking the head moves up and down linearly. A simultaneous angular pitching motion of the head keeps it aligned toward an imaginary point in space at a distance of about one meter in front of a subject and along the line of march. This distance is called the head fixation distance. Head fixation distance provides the fundamental framework necessary for understanding the functional significance of the vestibular reflexes that couple head motion to eye motion. This framework facilitates the intelligent design of counter-measures for the effects of exposure to microgravity upon the vestibular ocular reflexes. Ideal trajectory analysis is a simple candidate countermeasure based upon quantifying body sway during repeated up and down stair stepping. It provides one number that estimates the body sway deviation from an ideal sinusoidal body sway trajectory normalized on the subject's height. This concept has been developed with NSBRI funding in less than one year. These findings are explained in more detail below. Compared to assessments of the vestibuo-ocular reflex, analysis of vestibular effects on locomotor function is relatively less well developed and quantified. We are improving this situation by applying methodologies such as nonlinear orbital stability to quantify responses and by using multivariate statistical approaches to link together the responses across separate tests. In this way we can exploit the information available and increase the ability to discriminate between normal and pathological responses. Measures of stability and orientation are compared to measures such as dynamic visual acuity and with balance function tests. The responses of normal human subjects and of patients having well documented pathophysiologies are being characterized. When these studies are completed, we should have a clearer idea about normal and abnormal patterns of eye, head, and body movements during locomotion and their stability in a wide range of environments. We plan eventually to use this information to validate the efficacy of candidate neurovestibular and neuromuscular rehabilitative techniques. Some representative studies made during this year are summarized.

Wall, Conrad., III↗
NASA NTRS2004

Stochastic Stability of Sampled Data Systems with a Jump Linear Controller

In this paper an equivalence between the stochastic stability of a sampled-data system and its associated discrete-time representation is established. The sampled-data system consists of a deterministic, linear, time-invariant, continuous-time plant and a stochastic, linear, time-invariant, discrete-time, jump linear controller. The jump linear controller models computer systems and communication networks that are subject to stochastic upsets or disruptions. This sampled-data model has been used in the analysis and design of fault-tolerant systems and computer-control systems with random communication delays without taking into account the inter-sample response. This paper shows that the known equivalence between the stability of a deterministic sampled-data system and the associated discrete-time representation holds even in a stochastic framework.

Gonzalez, Oscar R.↗
NASA NTRS1979

Analysis of solar radiation pressure induced coupled librations of a gravity stabilized axisymmetric satellite

This paper presents an analysis of solar radiation pressure induced coupled librations of gravity stabilized cylindrical spacecraft with a special reference to geostationary communication satellites. The Lagrangian approach is used to obtain the corresponding equations of motion. The solar induced torques are assumed to be free of librational angles and are represented by their Fourier expansion. The response and periodic solutions are obtained through linear and nonlinear analyses, using the method of harmonic balance in the latter case. The stability conditions are obtained using Routh-Hurwitz criteria. To establish the ranges of validity the analytic response is compared with the numerical solution. Finally, values of the system parameters are suggested to make the satellite behave as desired. Among these is a possible approach to subdue the solar induced roll resonance. It is felt that the approximate analysis presented here should significantly reduce the computational efforts involved in the design and stability analysis of the systems.

Shrivastava, S. K.↗
NASA NTRS1979

A multiloop generalization of the circle criterion for stability margin analysis

In order to provide a theoretical tool suited for characterizing the stability margins of multiloop feedback systems, multiloop input-output stability results generalizing the circle stability criterion are considered. Generalized conic sectors with 'centers' and 'radii' determined by linear dynamical operators are employed to specify the stability margins as a frequency dependent convex set of modeling errors (including nonlinearities, gain variations and phase variations) which the system must be able to tolerate in each feedback loop without instability. The resulting stability criterion gives sufficient conditions for closed loop stability in the presence of frequency dependent modeling errors, even when the modeling errors occur simultaneously in all loops. The stability conditions yield an easily interpreted scalar measure of the amount by which a multiloop system exceeds, or falls short of, its stability margin specifications.

Safonov, M. G.↗
NASA NTRS1995

The initial-value problem for viscous channel flows

Plane viscous channel flows are perturbed and the ensuing initial-value problems are investigated in detail. Unlike traditional methods where traveling wave normal modes are assumed for solution, this works offers a means whereby completely arbitrary initial input can be specified without having to resort to eigenfunction expansions. The full temporal behavior, including both early time transients and the long time asymptotics, can be determined for any initial disturbance. Effects of three-dimensionality can be assessed. The bases for the analysis are: (a) linearization of the governing equations; (b) Fourier decomposition in the spanwise and streamwise directions of the flow; and (c) direct numerical integration of the resulting partial differential equations. All of the stability data that are known for such flows can be reproduced. Also, the optimal initial condition can be determined in a straight forward manner and such optimal conditions clearly reflect transient growth data that is easily determined by a rational choice of a basis for the initial conditions. Although there can be significant transient growth for subcritical values of the Reynolds number using this approach it does not appear possible that arbitrary initial conditions will lead to the exceptionally large transient amplitudes that have been determined by optimization of normal modes. The approach is general and can be applied to other classes of problems where only a finite discrete spectrum exists, such as the boundary layer for example.

Criminale, W. O.↗