Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Perturbative methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,441 records · Page 80

Investigation of electrodynamic stabilization and control of long orbiting tethers

An algorithm for using electric currents to control pendular oscillations induced by various perturbing forces on the Skyhook wire is considered. Transverse and vertical forces on the tether; tether instability modes and causes during retrieval by space shuttle; simple and spherical pendulum motion and vector damping; and current generation and control are discussed. A computer program for numerical integration of the in-plane and out-of-plane displacements of the tether vs time was developed for heuristic study. Some techniques for controlling instabilities during payload retrieval and methods for employing the tether for launching satellites from the space shuttle are considered. Derivations and analyses of a general nature used in all of the areas studied are included.

Arnold, D. A.↗

Radiative transfer in spherical atmospheres

A method for defining spherical model atmospheres in radiative/convective and hydrostatic equilibrium is presented. A finite difference form is found for the transfer equation and a matrix operator is developed as the discrete space analog (in curvilinear coordinates) of a formal integral in plane geometry. Pressure is treated as a function of temperature. Flux conservation is maintained within the energy equation, although the correct luminosity transport must be assigned for any given level of the atmosphere. A perturbed integral operator is used in a complete linearization of the transfer and constraint equations. Finally, techniques for generating stable solutions in economical computer time are discussed.

Kalkofen, W.↗

A simple, accurate and efficient algorithm for unsteady transonic flow

The derivation of an algorithm for computing an airfoil's response to small unsteady perturbations is described. The analysis of the far-field boundary conditions for the flow field, and the application of the time-linearization technique to develop an accurate description of the flow field are discussed. The treatment of moving shock waves is examined. The alternating direction implicit scheme of Ballhaus and Steger (1975) is utilized for the discretization of the governing equations and boundary conditions. Numerical examples demonstrating the applicability of the time-linearization method to determine the correlation between airfoil and pressure distribution are provided. The effect of wind tunnel walls on unsteady transonic testing is studied.

Fung, K.-Y.↗

On the ability of rough surface scattering approximations to predict hydrodynamic modulation of the ocean radar cross section - A numerical study

Scattering from the ocean surface at microwave frequencies is often described by the Bragg scattering model or the two-scale model. Here, the ability of these models to predict cross section modulation is investigated by comparing Bragg and two-scale cross-sections with cross sections calculated by a numerical electromagnetic approach. The numerical method is used as the standard of comparison. For slightly rough surfaces, the Bragg model is found to be very accurate, but for much rougher surfaces, the two-scale model is superior and should be used for calculating both the absolute cross section and the change in cross section due to spatial or temporal modulation of the ocean wave height spectrum. These results are relevant to the interpretation of SAR images of ocean surface phenomena which perturb the wave height spectrum at ripple wavelengths.

Durden, Stephen L.↗

Impact of Eccentricity on East-west Stationkeeping for GPS Class of Orbits

There exists a strong relationship between eccentricity and the potential for a repeating groundtrack orbit to exhibit chaotic motion. This is true at all values of eccentricity, but, perhaps most dramatic, is that it is true even for orbits that are nearly circular. These complex motions can have a significant impact on the east-west stationkeeping process for maintaining the repeating groundtrack property of a commensurate orbit. Ely and Howell have shown that traditional stationkeeping (SK) methods are unable to maintain a repeating groundtrack in the presence of complex dynamics, such as with chaotic motion. They developed an alternate SK method that is able to maintain a repeating groundtrack for eccentric, commensurate orbits. The focus of the current study is to investigate orbits with characteristics that are similar to GPS satellites except with modestly larger eccentricities. It will be shown that at eccentricities larger than approx. .01 the chaotic regions become significant, and the need arises for a robust stationkeeping approach, such as developed in. FurtheRmore, the investigation will reveal that the influence of luni-solar perturbations contributes to the growth of eccentricity, thus increasing the probability of encountering chaotic motion during a typical satellite lifetime.

Ely, Todd A.↗

On the determination and investigation of the terrestrial ionospheric refractive indices using GEOS-3/ATS-6 satellite-to-satellite tracking data

When the radio link between two satellites (GEOS-3/ATS-6) is intercepted by the earth's ionosphere and neutral atmosphere, a change in the Doppler frequency results. Travel through the atmosphere causes the Doppler phase to be advanced in the ionosphere's portion and retarded in the neutral portion of the atmosphere. Analysis of the shortening and lengthening of the phase of the Satellite-to-Satellite Tracking (SST) data that passed within 40-700 km above the earth's surface during its ATS-6 to GEOS-3 to ATS-6 path, caused by the atmosphere, results in refractivity versus height profiles. The SST Doppler data were used directly to adjust the GEOS-3 orbit. Perturbation from the Moon, Sun and a 15th order/degree earth gravity field were included in the orbit solution. This orbit was continued through the occultation period and a model ionosphere was estimated by a least-square adjustment of the Chapman ionosphere parameters from the SST data residuals. The refractivity profile obtained by this model ionosphere was compared to a refractivity profile obtained by a direct integral inversion of the SST data residuals. Systematic differences between the 2 methods were caused by orbital errors, which propagated into the solution. The SST data yielded refractive index profiles in a novel economical manner because no additional or special on-board equipment were required.

Liu, A. S.↗

The ages of globular cluster stars - Effects of rotation on pre-main-sequence, main-sequence, and turnoff evolution

Evolutionary sequences for low-metallicity stars (Z ranging from 0.001 to 0.0001) to study the effects of internal stellar rotation on the evolutionary time scales in the pre-main sequence, the main sequence (MS), and around the MS turnoff. Although a substantial amount of angular momentum remains in the interior, rotation is only a minor perturbation on the structure and ages of globular cluster stars. Even models with large initial angular momenta have MS lifetimes that are within 1 percent of those of standard models of the same mass and composition. Therefore, rotation does not affect age estimates of globular clusters from isochrone fitting. Furthermore, the models suggest that because rotation is not likely to affect horizontal-branch (HB) morphology, it does not affect significantly age estimates from the Delta-V method. Nevertheless, the internal angular momentum in the models is consistent with observations of surface rotational velocities on the HB, which require the preservation of a large reservoir of internal angular momentum.

Deliyannis, Constantine P.↗

Waveform synthesis of surface waves in a laterally heterogeneous earth by the Gaussian beam method

The present investigation is concerned with an application of the Gaussian beam method to surface waves in the laterally heterogeneous earth. The employed method has been developed for ray tracing and synthesizing seismograms of surface waves in cases involving the laterally heterogeneous earth. The procedure is based on formulations derived by Yomogida (1985). Vertical structure of the wave field is represented by the eigenfunctions of normal mode theory, while lateral variation is expressed by the parabolic equation as in two-dimensional acoustic waves or elastic body waves. It is demonstrated that a large-amplitude change can result from a slight perturbation in the phase velocity model.

Yomogida, K.↗

Unsteady transonic aerodynamics and aeroelastic calculations at low-supersonic freestreams

A computational procedure is presented to simulate transonic unsteady flows and corresponding aeroelasticity of wings at low-supersonic freestreams. The flow is modeled by using the transonic small-perturbation theory. The structural equations of motions are modeled using modal equations of motion directly coupled with aerodynamics. Supersonic freestreams are simulated by properly accounting for the boundary conditions based on pressure waves along the flow characteristics in streamwise planes. The flow equations are solved using the time-accurate, alternating-direction implicit finite-difference scheme. The coupled aeroelastic equations of motion are solved by an integration procedure based on the time-accurate, linear-acceleration method. The flow modeling is verified by comparing calculations with experiments for both steady and unsteady flows at supersonic freestreams. The unsteady computations are made for oscillating wings. Comparisons of computed results with experiments show good agreement. Aeroelastic responses are computed for a rectangular wing at Mach numbers ranging from subtransonic to upper-transonic (supersonic) freestreams. The extension of the transonic dip into the upper transonic regime is illustrated.

Guruswamy, Guru P.↗

Azimuthal spoke propagation in Hall Effect Thrusters

Spokes are azimuthally propagating perturbations in the plasma discharge of Hall Effect Thrusters (HETs) that travel in the E × B direction and have been observed in many different systems. The propagation of azimuthal spokes are investigated in a 6 kW HET known as the H6 using ultra-fast imaging and azimuthally spaced probes. A spoke surface is a 2-D plot of azimuthal light intensity evolution over time calculated from 87,500 frames/s videos. The spoke velocity has been determined using three methods with similar results: manual fitting of diagonal lines on the spoke surface, linear cross-correlation between azimuthal locations and an approximated dispersion relation. The spoke velocity for three discharge voltages (300, 400 and 450 V) and three anode mass flow rates (14.7, 19.5 and 25.2 mg/s) yielded spoke velocities between 1500 and 2200 m/s across a range of normalized magnetic field settings. The spoke velocity was inversely dependent on magnetic field strength for low B-field settings and asymptoted at B-field higher values. The velocities and frequencies are compared to standard drifts and plasma waves such as E × B drift, electrostatic ion cyclotron, magnetosonic and various drift waves. The empirically approximated dispersion relation yielded a characteristic velocity that matched the ion acoustic speed for 5 eV electrons that exist in the near-anode and near-field plume regions of the discharge channel based on internal measurements. Thruster performance has been linked to operating mode where thrust-to-power is maximized when azimuthal spokes are present so investigating the underlying mechanism of spokes will benefit thruster operation.

Polk, James E.↗

Redundancy management of electrohydraulic servoactuators by mathematical model referencing

A description of a mathematical model reference system is presented which provides redundancy management for an electrohydraulic servoactuator. The mathematical model includes a compensation network that calculates reference parameter perturbations induced by external disturbance forces. This is accomplished by using the measured pressure differential data taken from the physical system. This technique was experimentally verified by tests performed using the H-1 engine thrust vector control system for Saturn IB. The results of these tests are included in this report. It was concluded that this technique improves the tracking accuracy of the model reference system to the extent that redundancy management of electrohydraulic servosystems may be performed using this method.

Campbell, R. A.↗

On the determination and investigation of the terrestrial ionospheric refractive indices using GEOS-3/ATS-6 satellite-to-satellite tracking data

An analysis of the shortening and lengthening of the phase of satellite-to-satellite (SST) data that passed within 40-700 km above the earth surface during its ATS-6 to GEOS-3 to ATS-6 path resulted in refractivity vs height profiles. The SST Doppler data were used directly to adjust the GEOS-3 orbit. Perturbations from the moon, sun, and a 15th-order/degree earth gravity field were included in the orbit solution. This orbit was continued through the occulation period and a model ionosphere was estimated by a least-squares adjustment of the Chapman ionosphere parameters from the SST data residuals. The refractivity profile obtained by this model ionosphere was compared to a refractivity profile obtained by a direct integral inversion of the SST data residuals. Systematic differences between the two methods were caused by orbital errors, which propagated into the solution.

Liu, A. S.↗

Direct numerical solution of the transonic perturbation integral equation for lifting and nonlifting airfoils

The linear transonic perturbation integral equation previously derived for nonlifting airfoils is formulated for lifting cases. In order to treat shock wave motions, a strained coordinate system is used in which the shock location is invariant. The tangency boundary conditions are either formulated using the thin airfoil approximation or by using the analytic continuation concept. A direct numerical solution to this equation is derived in contrast to the iterative scheme initially used, and results of both lifting and nonlifting examples indicate that the method is satisfactory.

Nixon, D.↗

Bed Rest and Orthostatic-Hypotensive Intolerance

Orthostatic tolerance may be defined as the ability of humans to maintain cerebral perfusion and consciousness upon movement from a supine or sitting position to the upright posture; for example, subjects can stand suddenly or be tilted to the head-up body position. Similar but not identical physiological responses can be induced by positive G(sub Z) (head to foot) acceleration or exposure to lower body negative pressure (LBNP). The objective is to suddenly shift blood to the lower body to determine how effectively cardiovascular and neural-hormonal compensatory responses react to maintain blood pressure. In the most precise method for measuring tolerance, individuals would be stressed until they faint (syncope). However, the potential consequences and discomforts of such a test usually prohibit such a procedure so that few investigators actually induce syncope. In a more common approach, subjects are exposed to a given level of stress, for example, head-up tilt for 15 min, and any increases in heart rate or decreases in blood pressure are interpreted as indicators of progress toward syncope. Presumably, the greater the perturbation of heart rate and blood pressure, the closer to "tolerance," i.e., point of unconsciousness. Another more appropriate approach is to induce a progressively increasing hypotensive stress until pre-determined physiological responses or pre-syncopal symptoms appear. The physiological criteria may include a sudden drop in systolic blood pressure (greater than 25 mm/min), a sudden drop in heart rate (greater than 15 beats/min), or a systolic blood pressure less than 70 mmHg. The most common pre-syncopal symptoms include lightheadedness, stomach awareness or distress, feelings of warmth, tingly skin, and light to profuse sweating. Usually a combination of physiological responses and symptoms occurs such that, on different days, the tolerance time to the same orthostatic protocol is reproducible for a given individual. The assumption is that by taking subjects to near fainting, one can determine their tolerance. This latter pre-syncopal approach is better for estimating orthostatic or hypotensive tolerance than the former measurement of heart rate and blood pressure responses to a given stress. There is considerable variability in individual responses to orthostasis. For example, some subjects are "heart-rate responders" and have a pronounced cardiovascular response similar to that when performing moderately hard aerobic exercise, whereas others may experience pre-syncopal symptoms with very little increase in heart rate. Some individuals have a slow, gradual fall in blood pressure to orthostasis, and others have little change in blood pressure until a sudden precipitous fall in pressure occurs just prior to fainting. With both tilt and LBNP tests there is a low correlation between heart-rate or blood-pressure responses to a sub-tolerance stress as a measure of pre-syncopal limited orthostatic-hypotensive tolerance.

Schneider, Suzanne M.↗

Primal and Dual Integrated Force Methods Used for Stochastic Analysis

At the NASA Glenn Research Center, the primal and dual integrated force methods are being extended for the stochastic analysis of structures. The stochastic simulation can be used to quantify the consequence of scatter in stress and displacement response because of a specified variation in input parameters such as load (mechanical, thermal, and support settling loads), material properties (strength, modulus, density, etc.), and sizing design variables (depth, thickness, etc.). All the parameters are modeled as random variables with given probability distributions, means, and covariances. The stochastic response is formulated through a quadratic perturbation theory, and it is verified through a Monte Carlo simulation.

Patnaik, Surya N.↗

Azimuthal Spoke Propagation in Hall Effect Thrusters

Spokes are azimuthally propagating perturbations in the plasma discharge of Hall Effect Thrusters (HETs) that travel in the E x B direction and have been observed in many different systems. The propagation of azimuthal spokes are investigated in a 6 kW HET known as the H6 using ultra-fast imaging and azimuthally spaced probes. A spoke surface is a 2-D plot of azimuthal light intensity evolution over time calculated from 87,500 frames/s videos. The spoke velocity has been determined using three methods with similar results: manual fitting of diagonal lines on the spoke surface, linear cross-correlation between azimuthal locations and an approximated dispersion relation. The spoke velocity for three discharge voltages (300, 400 and 450 V) and three anode mass flow rates (14.7, 19.5 and 25.2 mg/s) yielded spoke velocities between 1500 and 2200 m/s across a range of normalized magnetic field settings. The spoke velocity was inversely dependent on magnetic field strength for low B-field settings and asymptoted at B-field higher values. The velocities and frequencies are compared to standard drifts and plasma waves such as E x B drift, electrostatic ion cyclotron, magnetosonic and various drift waves. The empirically approximated dispersion relation yielded a characteristic velocity that matched the ion acoustic speed for 5 eV electrons that exist in the near-anode and near-field plume regions of the discharge channel based on internal measurements. Thruster performance has been linked to operating mode where thrust-to-power is maximized when azimuthal spokes are present so investigating the underlying mechanism of spokes will benefit thruster operation.

Cylindrical Hall Thruster (CHT)↗

Uncertainty Assessment of Space-Borne Passive Soil Moisture Retrievals

The uncertainty associated with passive soil moisture retrieval is hard to quantify, and known to be underlain by various, diverse, and complex causes. Factors affecting space-borne retrieved soil moisture estimation include: (i) the optimization or inversion method applied to the radiative transfer model (RTM), such as e.g. the Single Channel Algorithm (SCA), or the Land Parameter Retrieval Model (LPRM), (ii) the selection of the observed brightness temperatures (Tbs), e.g. polarization and incidence angle, (iii) the definition of the cost function and the impact of prior information in it, and (iv) the RTM parameterization (e.g. parameterizations officially used by the SMOS L2 and SMAP L2 retrieval products, ECMWF-based SMOS assimilation product, SMAP L4 assimilation product, and perturbations from those configurations). This study aims at disentangling the relative importance of the above-mentioned sources of uncertainty, by carrying out soil moisture retrieval experiments, using SMOS Tb observations in different settings, of which some are mentioned above. The ensemble uncertainties are evaluated at 11 reference CalVal sites, over a time period of more than 5 years. These experimental retrievals were inter-compared, and further confronted with in situ soil moisture measurements and operational SMOS L2 retrievals, using commonly used skill metrics to quantify the temporal uncertainty in the retrievals.

Quets, Jan↗

Magnetospheric electric fields deduced from drifting whistler paths.

Results of a study in which the amplitude of the E-W component E sub w of the convection electric field in the nightside magnetosphere has been inferred from the observed cross-L motions of whistler ducts within the plasmasphere, and several ducts distributed over 1 to 2 earth radii in L space and over plus or minus deg 15 around the longitude of the Eights, Antarctica, whistler station have been tracked simultaneously. The method appears capable of resolving fluctuations in E sub w with period T equal to approximately 15 min and rms amplitude as low as 0.05 mV/m. For variations with T greater than 1 hour the method has a sensitivity of the order of 0.01 mV/m. Three case studies are presented, two of which illustrate convection activity associated with relatively isolated substorms. In these two cases E sub w reversed from westward to eastward for a period following the decay of substorm bay activity. In the third case the substorm bay activity was prolonged, and E sub w remained westward and at enhanced levels until local dawn. Evidence was found that, at least in a limited longitudinal sector, perturbing substorm E sub w fields can penetrate deep within the plasmasphere. In two of the case studies comparisons of E sub w and the interplanetary magnetic-field theta component show evidence of a possible relation based on brief (less than or equal to 1 hour) southward excursions but not on long preceding southward events.

Carpenter, D. L.↗