Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “perturbation 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 901 records · Page 50

In-Service Monitoring of Steam Pipe Systems at High Temperatures

A system and method for monitoring the properties of a fluid, such as water, in a steam pipe without mechanically penetrating the wall of the pipe. The system uses a piezoelectric transducer to launch an ultrasonic probe signal into the pipe. Reflected ultrasonic signals are captured in a transducer, which can be the same transducer that launched the probe signal. The reflected signals are subjected to data processing, which can include filtering, amplification, analog-to-digital conversion and autocorrelation analysis. A result is extracted which is indicative of a property of the fluid, such as a height of the condensed fluid, a cavitation of the condensed fluid, and a surface perturbation of the condensed fluid. The result can be recorded, displayed, and/or transmitted to another location. One embodiment of the system has been constructed and tested based on a general purpose programmable computer using instructions recorded in machine-readable non-volatile memory.

Bar-Cohen, Yoseph↗

Apparatus for and Method of Monitoring Condensed Water in Steam Pipes at High Temperature

A system and method for monitoring the properties of a fluid, such as water, in a steam pipe without mechanically penetrating the wall of the pipe. The system uses a piezoelectric transducer to launch an ultrasonic probe signal into the pipe. Reflected ultrasonic signals are captured in a transducer, which can be the same transducer that launched the probe signal. The reflected signals are subjected to data processing, which can include filtering, amplification, analog-to-digital conversion and autocorrelation analysis. A result is extracted which is indicative of a property of the fluid, such as a height of the condensed fluid, a cavitation of the condensed fluid, and a surface perturbation of the condensed fluid. The result can be recorded, displayed, and/or transmitted to another location. One embodiment of the system has been constructed and tested based on a general purpose programmable computer using instructions recorded in machine-readable non-volatile memory.

Lih, Shyh-Shiuh↗

A Super-Jupiter Orbiting A Late-Type Star: A Refined Analysis of Microlensing Event OGLE-2012-BLG-0406

We present a detailed analysis of survey and follow-up observations of microlensing event OGLE-2012-BLG-0406 based on data obtained from 10 different observatories. Intensive coverage of the light curve, especially the perturbation part, allowed us to accurately measure the parallax effect and lens orbital motion. Combining our measurement of the lens parallax with the angular Einstein radius determined from finite-source effects, we estimate the physical parameters of the lens system. We find that the event was caused by a 2.73 plus or minus 0.43 M (sub J) planet orbiting a 0.44 plus or minus 0.07 solar mass early M-type star. The distance to the lens is 4.97 plus or minus 0.29 kpc and the projected separation between the host star and its planet at the time of the event is 3.45 plus or minus 0.26 AU. We find that the additional coverage provided by follow-up observations, especially during the planetary perturbation, leads to a more accurate determination of the physical parameters of the lens.

microlensing↗

Analysis of ATS-5 data

Data from the particle detectors on ATS-5 in synchronous orbit and from the meridian line of magnetometers at a ground station are correlated for periods of well defined substorm activity. Specific substorm events are presented whose region of maximum perturbation can be accurately identified using current systems modelling techniques. It was found that the impulsive particle response at ATS-5 is correlated with changes in the auroral electrojet structure during the development of the substorm expansive phase. The response of the particle detectors appears to differ from event to event depending on the geometry of the disturbed region of the magnetosphere. It is concluded that the substorm disturbance typically begins at low latitude and propagates poleward in impulsive steps and that energetic electron enhancements are observed at ATS-5 when the poleward border of the electrojet intensifies in the latitude range of the ATS-5 field line.

Sharp, R. D.↗

Stochastic response analysis, order reduction, and output feedback controllers for flexible spacecraft

Real disturbances and real sensors have finite bandwidths. The first objective of this paper is to incorporate this finiteness in the 'open-loop modal cost analysis' as applied to a flexible spacecraft. Analysis based on residue calculus shows that among other factors, significance of a mode depends on the power spectral density of disturbances and the response spectral density of sensors at the modal frequency. The second objective of this article is to compare performances of an optimal and a suboptimal output feedback controller, the latter based on 'minimum error excitation' of Kosut. Both the performances are found to be nearly the same, leading us to favor the latter technique because it entails only linear computations. Our final objective is to detect an instability due to truncated modes by representing them as a multiplicative and an additive perturbation in a nominal transfer function. In an example problem it is found that this procedure leads to a narrow range of permissible controller gains, and that it labels a wrong mode as a cause of instability. A free beam is used to illustrate the analysis in this work.

Hablani, H. B.↗

Prospects for Eulerian CFD analysis of helicopter vortex flows

The applicability of current finite-volume CFD algorithms based on the Euler equations to the vortex flow over a helicopter in forward flight is investigated analytically. The general characteristics of the flow are reviewed; existing Euler, Navier-Stokes, perturbation, high-order, and adaptive methods are briefly characterized; and a novel Eulerian/Lagrangian approach with entropy and vorticity corrections is presented in detail. Numerical results for simple convection of a finite-core Lamb vortex moving downstream with its axis perpendicular to the flow are presented in graphs, and the possibility of extending the method to three-dimensional, viscous, and shock flows is discussed.

Drela, Mark↗

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

Approximate theory of large-amplitude wave propagation

An orbit perturbation procedure is applied to the description of monochromatic, large-amplitude, electrostatic plasma wave propagation. In the lowest-order approximation, untrapped electrons are assumed to follow constant-velocity orbits and trapped electrons are assumed to execute simple harmonic motion. The deviations of these orbits from the actual orbits are regarded as perturbations. The nonlinear damping rate and frequency shift are then obtained in terms of simple functions. The results are in good agreement with previous less approximate analyses. A significant feature of the analysis is that it treats a single wave by techniques previously applied to turbulent spectra. The analysis can consequently be extended to the case of a large-amplitude wave interacting with a lower-amplitude spectrum of waves.

Kim, H.↗

Comments on the Venus rotation pole

Possible orientations of the Venus rotation pole as a function of planetary oblateness are calculated, taking into account the variation of the orbital inclination and motion of the nodal line produced by long-term planetary perturbations and assuming the obliquity to be fully damped. If the obliquity is stabilized against solar tides by core-mantle viscous coupling, a fully damped obliquity is the expected state. An analysis of earth-based radar data from 1964 to 1977 yields a pole position which lies near damped pole positions of small oblateness, no more than about 10 to the 6th. Possible implications of this result are considered.

Ward, W. R.↗

Thermal perturbation of the Sun

Thermal perturbations of the solar convection zone can be modeled (to the first order) by perturbing the mixing length parameter alpha (equal to the ratio of the mixing length to the pressure scale height) used in the standard mixing length theory of convection. Results of such an analysis are presented and discussed in relation to recent work by others.

Twigg, L. W.↗

Control design for robust stability in linear regulators: Application to aerospace flight control

Time domain stability robustness analysis and design for linear multivariable uncertain systems with bounded uncertainties is the central theme of the research. After reviewing the recently developed upper bounds on the linear elemental (structured), time varying perturbation of an asymptotically stable linear time invariant regulator, it is shown that it is possible to further improve these bounds by employing state transformations. Then introducing a quantitative measure called the stability robustness index, a state feedback conrol design algorithm is presented for a general linear regulator problem and then specialized to the case of modal systems as well as matched systems. The extension of the algorithm to stochastic systems with Kalman filter as the state estimator is presented. Finally an algorithm for robust dynamic compensator design is presented using Parameter Optimization (PO) procedure. Applications in a aircraft control and flexible structure control are presented along with a comparison with other existing methods.

Yedavalli, R. K.↗

A controls engineering approach for analyzing airplane input-output characteristics

An engineering approach for analyzing airplane control and output characteristics is presented. State-space matrix equations describing the linear perturbation dynamics are transformed from physical coordinates into scaled coordinates. The scaling is accomplished by applying various transformations to the system to employ prior engineering knowledge of the airplane physics. Two different analysis techniques are then explained. Modal analysis techniques calculate the influence of each system input on each fundamental mode of motion and the distribution of each mode among the system outputs. The optimal steady state response technique computes the blending of steady state control inputs that optimize the steady state response of selected system outputs. Analysis of an example airplane model is presented to demonstrate the described engineering approach.

Arbuckle, P. Douglas↗

Lyapunov function-based control laws for revolute robot arms - Tracking control, robustness, and adaptive control

A new class of joint level control laws for all-revolute robot arms is introduced. The analysis is similar to a recently proposed energy-like Liapunov function approach, except that the closed-loop potential function is shaped in accordance with the underlying joint space topology. This approach gives way to a much simpler analysis and leads to a new class of control designs which guarantee both global asymptotic stability and local exponential stability. When Coulomb and viscous friction and parameter uncertainty are present as model perturbations, a sliding mode-like modification of the control law results in a robustness-enhancing outer loop. Adaptive control is formulated within the same framework. A linear-in-the-parameters formulation is adopted and globally asymptotically stable adaptive control laws are derived by simply replacing unknown model parameters by their estimates (i.e., certainty equivalence adaptation).

Wen, John T.↗

The Influence of Sound Cues on the Maintenance of Temporal Organization in the Sprague-Dawley Rat

Temporal organization is a fundamental property of living matter. From single cells to complex animals including man, most physiological systems undergo daily periodic changes in concert with environmental cues (e.g., light, temperature etc.). It is known that pulsed Environmental synchronizers, zeitgebers, (e.g. light) can modify rhythm parameters. Rhythm stability is a necessary requirement for most animal experiments. The extent to which sound can influence the circadian system of laboratory rats is poorly understood. This has implications to animal habitats in the novel environments of the Space-Laboratory or Space Station. A series of three white noise (88+/-0.82 db) zeitgeber experiments were conducted (n=6/experiment).The sound cue was introduced in the circadian free-running phase (DD-NQ) and in one additional case sound was added to the usual photoperiod (12L:12D) to determine masking effects. Circadian rhythm parameters of drinking frequency, feeding frequency, and gross locomotor activity were continuously monitored. Data analysis for these studies included macroscopic and microscopic methods. Raster plots to visually detect entrainment versus free-running period, were plotted for each animal, for all three parameters, during all sound perturbation tested. These data were processed through a series of detrending (robust locally weighted regression analyses) and complex demodulation analyses. In summary, these findings show that periodic "white" noise "influences" the rats circadian system but does not "entrain" the feeding, drinking or locomotor activity rhythms.

Winget, C. M.↗

Unsteady Extinction of Opposed Jet Ethylene/Methane HIFiRE Surrogate Fuel Mixtures vs Air

A unique idealized study of the subject fuel vs. air systems was conducted using an Oscillatory-input Opposed Jet Burner (OOJB) system and a newly refined analysis. Extensive dynamic-extinction measurements were obtained on unanchored (free-floating) laminar Counter Flow Diffusion Flames (CFDFs) at 1-atm, stabilized by steady input velocities (e.g., U(sub air)) and perturbed by superimposed in-phase sinusoidal velocity inputs at fuel and air nozzle exits. Ethylene (C2H4) and methane (CH4), and intermediate 64/36 and 15/85 molar percent mixtures were studied. The latter gaseous surrogates were chosen earlier to mimic ignition and respective steady Flame Strengths (FS = U(sub air)) of vaporized and cracked, and un-cracked, JP-7 "like" kerosene for a Hypersonic International Flight Research Experimentation (HIFiRE) scramjet. For steady idealized flameholding, the 100% C2H4 flame is respectively approx. 1.3 and approx.2.7 times stronger than a 64/36 mix and CH4; but is still 12.0 times weaker than a 100% H2-air flame. Limited Hot-Wire (HW) measurements of velocity oscillations at convergent-nozzle exits, and more extensive Probe Microphone (PM) measurements of acoustic pressures, were used to normalize Dynamic FSs, which decayed linearly with pk/pk U(sub air) (velocity magnitude, HW), and also pk/pk P (pressure magnitude, PM). Thus Dynamic Flame Weakening (DFW) is defined as % decrease in FS per Pascal of pk/pk P oscillation, namely, DFW = -100 d(U(sub air)/U(sub air),0Hz)/d(pkpk P). Key findings are: (1) Ethylene flames are uniquely strong and resilient to extinction by oscillating inflows below 150 Hz; (2) Methane flames are uniquely weak; (3) Ethylene / methane surrogate flames are disproportionately strong with respect to ethylene content; and (4) Flame weakening is consistent with limited published results on forced unsteady CFDFs. Thus from 0 to approx. 10 Hz and slightly higher, lagging diffusive responses of key species led to progressive phase lags (relative to inputs) in the oscillating flames, and caused maximum weakening. At 20 to 150 Hz, diffusion-rate-limited effects diminished, causing flames to "regain strengnth," and eventually become completely insensitive beyond 300 Hz. Detailed mechanistic understanding is needed. Overall, ethylene flames are remarkably resilient to dynamic extinction by oscillating inflows. They are the strongest, with the notable exception of H2. For HIFiRE tests, the 64%/36% surrogate disproportionally retains the high dynamic FS of ethylene, so the potential for loss of scramjet flameholding (flameout) due to low frequency oscillations is significantly mitigated.

Vaden, Sarah N.↗

Solar Radiation Pressure Binning for the Geosynchronous Orbit

Orbital maintenance parameters for individual satellites or groups of satellites have traditionally been set by examining orbital parameters alone, such as through apogee and perigee height binning; this approach ignored the other factors that governed an individual satellite's susceptibility to non-conservative forces. In the atmospheric drag regime, this problem has been addressed by the introduction of the "energy dissipation rate," a quantity that represents the amount of energy being removed from the orbit; such an approach is able to consider both atmospheric density and satellite frontal area characteristics and thus serve as a mechanism for binning satellites of similar behavior. The geo-synchronous orbit (of broader definition than the geostationary orbit -- here taken to be from 1300 to 1800 minutes in orbital period) is not affected by drag; rather, its principal non-conservative force is that of solar radiation pressure -- the momentum imparted to the satellite by solar radiometric energy. While this perturbation is solved for as part of the orbit determination update, no binning or division scheme, analogous to the drag regime, has been developed for the geo-synchronous orbit. The present analysis has begun such an effort by examining the behavior of geosynchronous rocket bodies and non-stabilized payloads as a function of solar radiation pressure susceptibility. A preliminary examination of binning techniques used in the drag regime gives initial guidance regarding the criteria for useful bin divisions. Applying these criteria to the object type, solar radiation pressure, and resultant state vector accuracy for the analyzed dataset, a single division of "large" satellites into two bins for the purposes of setting related sensor tasking and orbit determination (OD) controls is suggested. When an accompanying analysis of high area-to-mass objects is complete, a full set of binning recommendations for the geosynchronous orbit will be available.

Hejduk, M. D.↗

Contingency Trajectory Design for a Lunar Orbit Insertion Maneuver Failure by the LADEE Spacecraft

This paper presents results from a contingency trajectory analysis performed for the Lunar Atmosphere & Dust Environment Explorer (LADEE) mission in the event of a missed lunar-orbit insertion (LOI) maneuver by the LADEE spacecraft. The effects of varying solar perturbations in the vicinity of the weak stability boundary (WSB) in the Sun-Earth system on the trajectory design are analyzed and discussed. It is shown that geocentric recovery trajectory options existed for the LADEE spacecraft, depending on the spacecraft's recovery time to perform an Earth escape-prevention maneuver after the hypothetical LOI maneuver failure and subsequent path traveled through the Sun-Earth WSB. If Earth-escape occurred, a heliocentric recovery option existed, but with reduced science capacapability for the spacecraft in an eccentric, not circular near-equatorial retrograde lunar orbit.

LADEE↗

A Study of the Free Oscillations of the Earth

Published observations on the toroidal oscillations of the earth are critically reviewed. A supplementary analysis of the record obtained by the Lamont strain seismometer is presented. Eleven toroidal modes are identified, and it is concluded that the periods are known to within 1 per cent. A perturbation scheme involving the ratio of the angular velocity of the earth to the resonant frequency is used in calculating the effects due to the rotation of the earth on the resonant frequency. The free oscillations are viewed as a superposition of traveling waves. In a nonrotating system two traveling waves combine to produce a stationary standing wave. In a rotating system, the rotation distinguishes between waves that travel in the direction of rotation and those that travel in the opposite direction. Rotation removes a degeneracy and results in a splitting of a spectral peak of order l into 2 times l plus 1 peaks. The fractional displacement in frequency for the lowest-order toroidal oscillations is 1/206 and of the same order as the Q of the peak, so that splitting will probably not be observed in the toroidal oscillations. Viewed locally, rotation causes a particle to precess about a direction parallel to the axis of rotation. This precession will cause a variation of amplitude with time if the motion is recorded by an instrument with an anisotropic response function. Care is therefore needed in studying the time decay of a given spectral peak. Rotation also couples the normal coordinates so that a motion that is initially purely horizontal will develop a vertical component. It is expected that vertical seismometers should record particle motion with the toroidal frequencies. The perturbations of the toroidal oscillations due to core-mantle interaction are treated in detail. An exact expression is obtained for the rate of energy dissipated by a finitely conducting plate oscillating across a magnetic field. The energy dissipated at the core-mantle boundary due to viscous and hydromagnetic coupling is shown to be insignificant as compared with the energy dissipated within the mantle. The toroidal magnetic field leaking into the lower mantle combines with the dipole field, resulting in a stress on the mantle, tending to stiffen the lower boundary. The stress is of sufficient magnitude to produce a displacement toward higher frequency in the lower-order toroidal oscillations. Observations on the (sub 0) T (sub 2) oscillations lead to an estimate of the toroidal magnetic field in the lower mantle. A calculation of elastic energy in the low-order oscillations suggests a value of 10 (sup 18) ergs per cycles per hour for the energy density at low frequencies in the Chilean earthquake. Each mode of oscillation has a characteristic radial distribution of elastic energy associated with it. This distribution determines which parts of the earth contribute most heavily in determining a particular resonant frequency. The distribution of energy for the lower 17 modes for a homogeneous and a Gutenberg model earth is calculated. The resonant frequencies for models of the earth based on the Gutenberg and Lehmann distribution of elastic properties are presented. It is shown that the Gutenberg model earth fits the observations more closely than the Lehmann model and that a slight alteration of the Gutenberg model gives a significantly better fit to the observations. The alteration involves a lower shear-wave velocity in the lower mantle while the Gutenberg velocity distribution is maintained in the upper mantle. Various studies of the earth's oscillations coupled with surface-wave investigations substantiate Gutenberg's hypothesis of a layer of low velocity in the upper mantle. The physical conditions required for the formation of a region of low velocity are examined in detail. The results confirm Birch's earlier statement that a temperature gradient in excess of 6 degrees to 7 degrees per kilometer is needed to produce a decrease in velocity. The low-velocity layer does not require that the temperature approach or exceed the melting temperature. If tile upper mantle is homogeneous, the region of lower velocity should commence at the base of the crust and extend to 150 kilometers under the oceans and about 100 kilometers under continental regions. The distribution of thermal conductivity and radioactivity consistent with the low-velocity layer is also considered.

MacDonald, Gordon J. F.↗