Search NASA⌕ Search

SEARCH · Search NASA

Results for “mean square radius”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1974

A numerical study of electromagnetic scattering from ocean-like surfaces

The integral equations describing electromagnetic scattering from cylindrical, perfectly conducting surfaces are formulated and numerical results are presented. The results are compared with those obtained using approximate methods such as physical optics, geometrical optics, and perturbation theory. The integral equation solutions show that the surface radius of curvature must be greater than 2.5 wavelengths for either the physical optics or geometrical optics to give satisfactory results. It has also been shown that perturbation theory agrees with the exact fields as long as the root-mean-square surface roughness is less than one tenth of a wavelength.

Lentz, R. R.↗
NASA NTRS2002

Equal-Curvature X-Ray Telescopes

We introduce a new type of x-ray telescope design; an Equal-Curvature telescope. We simply add a second order axial sag to the base grazing incidence cone-cone telescope. The radius of curvature of the sag terms is the same on the primary surface and on the secondary surface. The design is optimized so that the on-axis image spot at the focal plane is minimized. The on-axis RMS (root mean square) spot diameter of two studied telescopes is less than 0.2 arc-seconds. The off-axis performance is comparable to equivalent Wolter type 1 telescopes.

Saha, Timo T.↗
NASA NTRS1988

Spatial Distribution of Lightning Strikes to Ground During Small Thunderstorms in Florida

The spatial patterns of the strike points produced by cloud-to-ground lightning under three small thunderstorms have been analyzed to determine the area flash density as a function of radius from the storm center, the distribution of the nearest-neighbor distances, and the distribution of the horizontal distances between successive flashes. The storm average flash densities range from about 0.8 to 1.6 Fl/km squared, and the average lightning fluxes range from 0.03 to 0.05 Fl/km squared/min. The mean nearest-neighbor distances are about 0.7 km and smaller, but are still in good agreement with a theory that assumes an infinite and uniform flash density. The mean distance between successive flashes ranges from 3.2 to 4.2 km, and a sizable fraction of this variation could be due to channel geometry.

Krider, E. Philip↗
NASA NTRS1993

Three-dimensional Structure of Straight and Curved Plane Wakes

Although the plane wake is marked by the formation of strong spanwise vortices, the initially two-dimensional Karman-like vortices soon develop a three-dimensional structure in the form of secondary streamwise vortices. So far, this streamwise vortex structure has been studied mostly through flow visualization and at relatively low Reynolds numbers. The primary objective of the present program was to investigate the origin and evolution of the three-dimensional structure of straight and curved plane wakes at relatively high Reynolds numbers (Re(sub b) = 28,000) through detailed measurements of the mean and turbulent properties at several streamwise locations. The experiments were conducted in three phases. In the first phase, the development of a straight plane wake was investigated. In the second phase, the effects of imposed streamwise curvature on the wake development were examined. The streamwise curvature was of constant radius and very mild in terms of the curvature ratio (b/(square root of R) is less than 2 percent). In both the first and second phases, the role of initial conditions was examined in wakes generated from both untripped (laminar) and tripped (turbulent) initial boundany layers. In the third phase, the effects of injecting streamwise vorticity and the effects of increased Reynolds number on the tripped wake structure and development were investigated.

Weygandt, James H.↗
NASA NTRS1983

On obtaining the forward phase functions of Saturn ring features from radio occultation observations

It is noted that the near-forward scattering functions of particles in Saturn ring features are related to 3.6 cm radio occultation power spectra by a Fredholm integral equation of the first kind. The equation reduces to an algebraic system of equation whose solution by usual inversion techniques (that is, least mean squares) is ruled out by the near singularity of the forward transformation matrix. A combination of constrained linear inversion and a filtering algorithm based on eigenvector decomposition of the matrix reduces the instabilities; this yields derived phase functions valid over the range of zero to about 12 mrad. It is noted that these functions represent the collective forward diffraction lobe of particles greater than about 1 m in radius. Since multiple scattering of the signal is a significant effect, the measured phase functions must be adjusted to obtain the singly scattered component. This single-scattering correction is examined for two physical models, namely the monolayer and the classical discrete random slab; in addition, the fraction of opacity in submeter particles for each model for particular ring features is estimated.

Zebker, H. A.↗
NASA NTRS1991

Magellan - Electrical and physical properties of Venus' surface

Magellan probes Venus'surface by 12.6-cm-wavelength vertical and oblique radar scattering and measures microwave thermal emission. Emissivity and root-mean-square slope maps between 330 deg and 30 deg E and 90 deg N and 80 deg S are dissimilar, although some local features are exceptions. Inferred surface emissivities typically are 0.85, but vary from 0.35 at Maxwell to 0.95 northeast of Gula Mons and other locations. Lowest emissivities appear in topographically high areas; this relation suggests that a phase change or differences in chemical weathering occur at about 6055-kilometer radius. Initial results indicate that there are significant variations in the surface scattering function.

Tyler, G. Leonard↗
NASA NTRS2011

Influence of Convection and Aerosol Pollution on Ice Cloud Particle Effective Radius

Satellite observations show that ice cloud effective radius (r(sub e)) increases with ice water content (IWC) but decreases with aerosol optical thickness (AOT). Using least-squares fitting to the observed data, we obtain an analytical formula to describe the variations of r(sub e) with IWC and AOT for several regions with distinct characteristics of r(sub e) -IWC-AOT relationships. As IWC directly relates to convective strength and AOT represents aerosol loading, our empirical formula provides a means to quantify the relative roles of dynamics and aerosols in controlling r(sub e) in different geographical regions, and to establish a framework for parameterization of aerosol effects on r(sub e) in climate models.

Jiang, J. H.↗
NASA NTRS1999

Buoyancy Driven Shear Flows of Bubble Suspensions

In this work the gas volume fraction and the root-mean-squared fluid velocity are measured in buoyancy driven shear flows of bubble suspensions in a tall, inclined, rectangular channel. The experiments are performed under conditions where We << 1a nd Re >> 1, for which comparisons are made with kinetic theory and numerical simulations. Here Re = gamma(a(exp 2)/nu is the Reynolds number and We = rho(gamma(exp 2))a(exp 3)/sigma is the Weber number; gamma is the shear rate, a is the bubble radius, nu is the kinematic viscosity of the liquid, rho is the density of the liquid, and sigma is the surface tension of the gas/liquid interface. Kang et al. calculated the bubble phase pressure and velocity variance of sheared bubble suspensions under conditions where the bubbles are spherical and the liquid phase velocity field can be approximated using potential flow theory, i.e. We= 0 and Re >> 1. Such conditions can be achieved in an experiment using gas bubbles, with a radius of O(0.5mm), in water. The theory requires that there be no average relative motion of the gas and liquid phases, hence the motivation for an experimental program in microgravity. The necessity of performing preliminary, Earth based experiments, however, requires performing experiments where the gas phase rises in the liquid, which significantly complicates the comparison of experiments with theory. Rather than comparing experimental results with theory for a uniform, homogeneous shear flow, experiments can be compared directly with solutions of the averaged equations of motion for bubble suspensions. This requires accounting for the significant lift force acting on the gas phase when the bubbles rise parallel to the average velocity of the sheared suspension. Shear flows can be produced in which the bubble phase pressure gradient, arising from shear induced collisions amongst the bubbles, balances a body force (centrifugal or gravitational) on the gas phase. A steady, non-uniform gas volume fraction can be measured, from which the bubble phase pressure gradient can be obtained and compared to theory and numerical simulations. The presence of bounding walls further complicates the experiments, since the detailed interactions of the bubbles with bounding walls is not well understood, especially in the presence of gravity, where the momentum and energy exchange depends on the inclination of the wall.

Koch, D. L.↗
NASA NTRS1979

A model of the Venus atmosphere from radio, radar, and occultation observations

A model is presented of the atmosphere and surface of Venus which best fits in the least-squares sense the available radio-brightness, radar cross-section, radio interferometric, and Mariner 5 and 10 radio occultation observations. The determinations of the radius of the planet obtained by others from radar time-delay measurements are included in the data set. The values of the adjusted parameters are: molar fraction of CO2 = 95 plus or minus 3%; fraction of combined nitrogen and argon = 5 plus or minus 3%; total atmospheric opacity at a wavelength of 1 cm = 19.4 plus or minus 1.3; mean radius of the surface = 6050.7 plus or minus 0.8 km; mean dielectric constant of the surface = 4.1 plus or minus 0.2, and percentage of total opacity due to chemical species other than CO2 = 45 plus or minus 12. The model temperature and pressure at the mean surface are 755 K and 91.4 atm, respectively.

Muhleman, D. O.↗
NASA NTRS1999

Dynamics of Single Chains of Suspended Ferrofluid Particles

We present an experimental study of the dynamics of isolated chains made of super-paramagnetic particles under the influence of a magnetic field. The motivation of this work is to understand if the chain fluctuations exist and, if it does, how does the fluctuation affect chain aggregation. We find that single chains strongly fluctuate and that the characteristic frequency of their fluctuations is inversely proportional to the magnetic field strength. The higher the field the lower the characteristic frequency of the chain fluctuations. In the high magnetic field limit, chains behave like rigid rods without any internal motions. In this work, we used ferrofluid particles suspended in water. These particles do not have any intrinsic magnetization. Once a magnetic field is applied, a dipole moment is induced in each particle, proportional to the magnetic field. A dipolar magnetic interaction then occurs between particles. If dipole-dipole magnetic energy is higher than the thermal energy, the result is a structure change inside the dipolar fluid. The ratio of these two energies is expressed by a coupling constant lambda as: lambda = (pi(a(exp 3))(chi(exp 2))(mu(sub 0))(H(sub 0))(exp 2))/18kT Where a is the particle radius, mu(sub 0) is the vacuum magnetic permeability, H(sub 0) the applied magnetic field, k the Boltzmann constant and T the absolute temperature. If lambda > 1, magnetic particles form chains along the field direction. The lateral coalescence of several chains may form bigger aggregates especially if the particle volume fraction is high. While many studies and applications deal with the rheological properties and the structural changes of these dipolar fluids, this work focuses on the understanding of the chain dynamics. In order to probe the chain dynamics, we used dynamic light scattering (DLS) in self-beating mode as our experimental technique. The experimental geometry is such that the scattering plane is perpendicular to the magnetic field. Therefore, only motions in this plane are probed. A very dilute sample of a ferrofluid emulsion with a particle volume fraction of 10(exp -5) is used in this experiment. We chose such a low volume fraction to avoid multiple light scattering as well as lateral chain-chain aggregation. DLS measures the dynamic structure factor S(q,t) of the sample (q is the scattering wave vector, t is the time). In the absence of the magnetic field, identical particles of ferrofluid droplets are randomly distributed and S(q,t) reduces to exp(-q(exp 2)2D(sub 0)t). D(sub 0)=(kT/(6(pi)(eta)(a)) is the diffusion coefficient of Brownian particles (where Xi = (6(pi)(eta)(a)) is the Stokes frictional coefficient of a spherical particle in a fluid of viscosity eta). If interactions or polydispersity can not be ignored, an effective diffusion coefficient is introduced. Formally, D(sub eff) is defined as: D(sub eff) = - q(exp -2) partial derivative of (ln(S(q,t)) with respect to time, as t goes to 0. D(sub eff) reduces to D(sub 0) if no interactions and only a few particles size are present. Therefore, we can use DLS to measure particle size. The particle radius was found to be a=0.23 mu m with 7% of polydispersity. In this case, if we vary the scattering angle theta (and so q) we do not have any change in the measured diffusion coefficient: it is q-independent. When a magnetic field is applied, particles aggregate into chains if lambda > 1. We first studied the kinetics of the chain formation when lambda = 406. At a fixed scattering angle, we measured diffusion coefficient D(sub eff) as a function of time. Experimentally, we find that D(sub eff) decreases monotonously with time. Physically, this means that chains are becoming longer and longer. Since we are only sensitive to motions in the scattering plane and since chains have their main axis perpendicular to this plane, the measured diffusion coefficient is the trans-verse diffusion coefficient. We can relate D(sub eff) to the mean number of particles per chain N(t) at a given time and to the diffusion coefficient of an isolated particle D(sub 0) as D(sub eff)=f(N(t))D(sub 0). Since f(N) is known from other recent work, N can be expressed as a function of the time. We found a square root dependency: N(t) proportional to the square root of t. As expected for very low volume fraction, this behavior is characteristic of a diffusion-limited aggregation as suggested by several authors and by our previous work. In this study, we focus on the dependence of the effective diffusion coefficient on the scattering angle and the magnetic field strength. After the magnetic field is applied (lambda = 406) for a long time, typically 6 hours, kinetics of chain formation becomes very slow. Chain size does not vary much over the next hour period. Thus, we can perform different interesting experiments. First, at a fixed magnetic field, we measure the effective diffusion coefficient as a function of the scattering angle (from 5 to 130 deg). Our results show that the measured diffusion coefficient increases linearly with the scattering angle: D(sub eff) proportional to q. If we do the same experiment for different lambda values, D(sub eff) depends on lambda as D(sub eff) proportional to lambda(exp -1/2). We also find for different lambda values that the same asymptotic D(sub eff) value is obtained when q approaches zero. The angle dependency of D(sub eff) suggests that an additional motion exists besides chain drifting. Chain size is constant during experiment, which was verified by measuring the same diffusion coefficient at the beginning and at the end of the angle switching. If chains are rigid, D(sub eff) is independent of q. Therefore, we found that D(sub eff) not only measures the motion of the entire chain but also its internal fluctuations. These internal motions are the fluctuations of the particles in the chain. To understand the q dependency of D(sub eff), let us look at the probing length used. In our study, the characteristic length scale probed is l=2pi/q which is in the range of 0.9<l/a<20. When l is much larger than the particle radius (l/a going to infinity) we are mainly sensitive to the center of mass diffusion of the chain. D(sub eff) is then the diffusion of the entire chain and depends only on N but not on q and lambda. The value D(sub eff) allows us to obtain the number of particles per chain N. In the opposite limit (i.e. l/a < 1), we are sensitive to motions on the size of individual particles. The main contributions to the measured diffusion coefficient come from internal motions of the chains (i.e. particles' fluctuations inside the chain). We investigated also the effect of the chain size on D(sub eff). For the same value of the magnetic interaction lambda and different chain sizes, we found that the slope of D(sub eff) versus q decreases when N increases as a power law D(sub eff) proportional to N(exp -0.7). This is an important behavior because this means that the longer the chain the slower its internal fluctuations besides the whole chain motion.

Cutillas, S.↗
NASA NTRS1989

Ginzburg-Landau theory for the solid-liquid interface of bcc elements. II - Application to the classical one-component plasma, the Wigner crystal, and He-4

The previously developed Ginzburg-Landau theory for calculating the crystal-melt interfacial tension of bcc elements to treat the classical one-component plasma (OCP), the charged fermion system, and the Bose crystal. For the OCP, a direct application of the theory of Shih et al. (1987) yields for the surface tension 0.0012(Z-squared e-squared/a-cubed), where Ze is the ionic charge and a is the radius of the ionic sphere. Bose crystal-melt interface is treated by a quantum extension of the classical density-functional theory, using the Feynman formalism to estimate the relevant correlation functions. The theory is applied to the metastable He-4 solid-superfluid interface at T = 0, with a resulting surface tension of 0.085 erg/sq cm, in reasonable agreement with the value extrapolated from the measured surface tension of the bcc solid in the range 1.46-1.76 K. These results suggest that the density-functional approach is a satisfactory mean-field theory for estimating the equilibrium properties of liquid-solid interfaces, given knowledge of the uniform phases.

Zeng, X. C.↗
NASA NTRS1994

Mean sea surface and geoid gradient comparisons with TOPEX altimeter data

Cycles 4 to 54 of TOPEX data have been analyzed through comparisons with the mean sea surface given on the disturbed geophysical data record (GDR). Two inverted barometer correction procedures were considered for the data reduction. One used a constant atmospheric pressure for all data while the one adopted for use, for most computations, introduced a cycle average pressure. The maximum difference between the two estimates was 3.0 cm with a clear annual signal. With the modified correction the TOPEX sea surface was compared to The Ohio State University (OSU) mean sea surface, given on the GDR, to estimate three translations ( delta x = -2.3 cm; delta y = 25.0 cm; delta z = -0.3 cm) and a bias (43.3 cm) between the two surfaces. The only significant translation is delta y which indicates the reference frame of the TOPEX system differs from that used in the OSU mean sea surface system. The bias between the TOPEX mean sea surface and the OSU mean sea surface was used to estimate an equatorial radius of 6,378,136.55 m based on an 18-cm biased estimate of the TOPEX altimeter. Examination of the average difference, by cycle, between the TOPEX sea surface and the OSU mean sea surface suggested a bias change of 3.1 +/- 2.2 mm/yr with a positive sign indicating the average ocean surface is rising or the altimeter measured distance is decreasing. Models were implemented that solved directly for a bias, bias rate annual/semiannual, and tide correction terms. The computations indicated that a simultaneous solution for this bias, bias rate, and annual/semiannual terms gave the most accurate results. Nonsimultaneous solutions led to slightly different bias rate values. The root mean square difference between the TOPEX sea surface and OSU sea surface, after translation and bias correction, was +/- 17 cm for a typical cycle. Some locations were indentified where the difference could reach 2.3 cm and were repeated over several cycles indicating errors in the mean sea surface. Most of the large differences occur in regions lacking altimeter data prior to the TOPEX/POSEIDON mission and/or areas of significant bathymetric signature. Geoid gradients are needed for the reduction of individual track data to a reference track. The accuracy of the determination of such gradients was determined through the comparison of predicted along-track gradients to the observed gradients. Among four mean sea surfaces tested the best agreement was found with the OSU mean sea surface placed on the TOPEX geophysical data record.

Rapp, Richard H.↗
NASA NTRS1982

Equivalent source modeling of the main field using MAGSAT data

An iterative least squares estimation algorithm with the capability for including a priori statistical information was implemented to recover multiple magnetic dipole models of the Earth's main magnetic field. The dipoles are fixed to a specified radius at or below the core-mantle boundary and centered on equal area blocks. The algorithm can solve for dipole magnitudes only (fixed orientations), or allow full freedom of orientation and solve for vector components. External field parameters and observatory anomaly biases can also be estimated simultaneously. Time dependence is modeled using first time derivatives for dipole vector components. Single-epoch and time dependent dipole models are derived using MAGSAT and observatory annul means data. Equivalent spherical harmonic representation are computed in closed form from the dipole models and compared with truncated spherical harmonic models estimated in the standard way from the same data sets. In particular, a 21 deg spatial resolution model based on 93 dipoles was computed based on observatory annual means data and a selected MAGSAT data set and was compared with candidate IGRF 1975 models and their 1980 secular variation.

Mayhew, M. A.↗
NASA NTRS2002

XTE J1908+094

Goddard Space Flight Center reports the serendipitous discovery of a new x-ray transient, XTE J1908+094, in RXTE (Rossi X-ray Timing Explorer) PCA (Proportional Counter Array) observations of the soft-gamma-ray repeater SGR 1900+14, triggered following the burst activity on Feb. 17-18 (GCN 1253). These observations failed to detect the 5.2-s SGR pulsations, pointing towards a possible new source as the origin of the high x-ray flux. An RXTE PCA scan of the region around SGR 1900+14 on Feb. 21 was consistent with emission only from known sources (and no new sources). However, the scans required SGR 1900+14 to be 20 times brighter than its quiescent flux level (GCN 1256). A Director's Discretionary Time Chandra observation on Mar. 11 showed that the SGR was quiescent and did not reveal any new source within the Chandra ACIS (Advanced CCD (charge coupled device) Imaging Spectrometer) field-of-view. A subsequent RXTE PCA scan on Mar. 17, taken in combination with the first scan, required that a new source be included in the fit. The best-fit position is R.A. 19h 08m 50s, Decl. = +9 22 deg .5 (equinox J2000.0; estimated 2 deg systematic error radius), or approximately 24 deg away from the SGR source. The source spectrum (2-30 kev) can be best fit with a power-law function including photoelectric absorption (column density N_h = 2.3 x 10(exp 22), photon index = 1.55). Iron line emission is present, but may be due to the Galactic ridge. Between Feb. 19 and Mar. 17, the source flux (2-10 keV) has risen from 26 to 64 mCrab. The power spectrum is flat between 1 mHz and 0.1 Hz, falling approximately as 1/f**0.5 up to 1 Hz. At 1 Hz is seen a broad quasiperiodic oscillation peak and a break to a 1/f**2 power law, which continues to 4 Hz. The fractional rms (root mean square) amplitude from 1 mHz to 4 Hz is 43 percent. No coherent pulsations are seen between 0.001 and 1024 Hz. The authors conclude that XTE J1908+094 is a new blackhole candidate.

Woods, P. M.↗
NASA NTRS1996

Supersonic Channel Concept for Enhancement of Lift/Drag Ratio

The concept studied during the summer NASA/ASEE Fellowship provides a means of lowering drag and a means for directional control of supersonic and hypersonic vehicles. Low drag and efficient directional control are essential for the success of aircraft, atmospheric entry vehicles, missiles, and other vehicles in supersonic and hypersonic flight. Drag reduction can result in increased vehicle range, increased speed, improved fuel efficiency, increased lift/drag ratio, and increased climb rate. For high supersonic and hypersonic vehicles heat transfer considerations dictate the design of the nose and leading edge. The heat transfer to such vehicles is most severe at stagnation points which occur on the leading edges and nose of the vehicle. Theoretical formulations, experimental data, and semi-empirical formulas all agree in the fact that stagnation point heat transfer is inversely proportional to the square root of the nose or leading edge radius. Thus, the noses and leading edges of supersonic and hypersonic vehicles are typically blunted so that the heat transfer and structural loads will be manageable. However, much of the wave drag experienced by these vehicles is due to nose blunting.

Ruffin, Stephen M.↗
NASA NTRS2002

Discovery of a 450 Hz Quasi-Periodic Oscillation from the Microquasar GRO J1655-40 with the Rossi X-ray Timing Explorer

We report the discovery with the Proportional Counter Array on board the Rossi X-Ray Timing Explorer of a 450 Hz quasi-periodic oscillation (QPO) in the hard X-ray flux from the Galactic microquasar GRO J1655-40. This is the highest frequency QPO modulation seen to date from a black hole. The QPO is detected only in the hard X-ray band above approx. 13 keV. It is both strong and narrow, with a typical rms (root mean square) amplitude of 4.5% in the 13-27 keV range and a width of approx. 40 Hz (FWHM). For two observations in which we detect the 450 Hz QPO, a previously known approx. 300 Hz QPO is also observed in the 2-13 keV band. We show that these two QPOs sometimes appear simultaneously, thus demonstrating the first detection of a pair of high-frequency QPOs in a black hole system. Prior to this, pairs of high-frequency QPOs have been detected only in neutron star systems. GRO J1655-40 is one of only a handful of black hole systems with a good dynamical mass constraint. For a nonrotating black hole with mass between 5.5 and 7.9 solar masses, the innermost stable circular orbit (ISCO) ranges from 45 to 70 km. For any mass in this range the radius at which the orbital frequency reaches 450 Hz is less than the ISCO radius, indicating that, if the modulation is caused by Kepler motion, the black hole must have appreciable spin. If the QPO frequency is set by the orbital frequency of matter at the ISCO, then for this mass range the dimensionless angular momentum lies in the range 0.15 < j < 0.5. Moreover, if the modulation is caused by oscillation modes in the disk or Lense-Thirring precession, then this would also require a rapidly rotating hole. We briefly discuss the implications of our findings for models of X-ray variability in black holes and neutron stars.

Strohmayer, Tod E.↗
NASA NTRS2002

Observer-Based Magnetic Bearing Controller Developed for Aerospace Flywheels

A prototype of a versatile, observer-based magnetic bearing controller for aerospace flywheels was successfully developed and demonstrated on a magnetic bearing test rig (see the photograph) and an actual flywheel module. The objective of this development included a fast, yet low risk, control development process, and a robust, high-performance controller for a large variety of flywheels. This required a good system model, an efficient development procedure, and a model-based controller that addressed the key problems associated with flywheel and bearing imbalance, sensor error, and vibration. The model used in this control development and tuning procedure included the flexible rotor dynamics and motor-induced vibrations. Such a model was essential for low-risk scheduling of speed-dependent control parameters and for reliable evaluation of novel control strategies. The successfully tested control prototype utilized an extended Kalman filter to estimate the true rotor principal-axis motion from the raw sensor position feedback. For control refinement, the extended Kalman filter also estimated and eliminated the combined effects of mass-imbalance and sensor runouts from the input data. A key advantage of the design based on the extended Kalman filter is its ability to accurately estimate both the rotor's principal-axis position and gyroscopic rates with the least amount of phase lag. This is important for control parameter scheduling to dampen the gyroscopic motions. Because of large uncertainties in the magnetic bearing and imbalance characteristics, this state-estimation scheme alone is insufficient for containing the rotor motion within the desired 1-mil excursion radius. A nonlinear gain adjustment based on an estimation of the principal-axis orbit size was needed to provide a coarse (nonoptimal), but robust, control of the orbit growth. Control current minimization was achieved with a (steepest gradient) search of synchronous errors in the principal-axis position input data. Actual flywheel tests of this observer-based controller (developed entirely in-house) at the NASA Glenn Research Center showed that the model correctly predicted the rotor orbit growth as a function of rotational speed, and it demonstrated the capability of gain adjustments to arrest this growth. Data from these tests on an actual flywheel module spun to 26,000 rpm proved that the controller was able to contain the shaft motion to within much less than 0.5 mils of radial excursion with axis currents less than 300 mA in root-mean-square estimate. The test speed range was limited because of thermal expansion concerns for this particular flywheel unit, not because of any deficiency in the controller. Simulations for this unit indicated that the controller should be robust up to its top operating speed of 60,000 rpm. Aside from these important achievements, and most significantly, it took less than 1 week to adapt this controller from the simple test rig to the actual flywheel and to demonstrate full five-axis levitation and control. This demonstration showed that both the controller and the model-based development and tuning framework are easily adaptable to a wide range of rotors and bearing configurations and, hence, are capable of reducing design risks and costs for many future flywheel technology developments.

Le, Dzu K.↗
NASA NTRS2011

Influence of Vibrotactile Feedback on Controlling Tilt Motion After Spaceflight

We hypothesize that adaptive changes in how inertial cues from the vestibular system are integrated with other sensory information leads to perceptual disturbances and impaired manual control following transitions between gravity environments. The primary goals of this ongoing post-flight investigation are to quantify decrements in manual control of tilt motion following short-duration spaceflight and to evaluate vibrotactile feedback of tilt as a sensorimotor countermeasure. METHODS. Data is currently being collected on 9 astronaut subjects during 3 preflight sessions and during the first 8 days after Shuttle landings. Variable radius centrifugation (216 deg/s, <20 cm radius) in a darkened room is utilized to elicit otolith reflexes in the lateral plane without concordant canal or visual cues. A Tilt-Translation Sled (TTS) is capable of synchronizing pitch tilt with fore-aft translation to align the resultant gravitoinertial vector with the longitudinal body axis, thereby eliciting canal reflexes without concordant otolith or visual cues. A simple 4 tactor system was implemented to provide feedback when tilt position exceeded predetermined levels in either device. Closed-loop nulling tasks are performed during random tilt steps or sum-of-sines (TTS only) with and without vibrotactile feedback of chair position. RESULTS. On landing day the manual control performance without vibrotactile feedback was reduced by >30% based on the gain or the amount of tilt disturbance successfully nulled. Manual control performance tended to return to baseline levels within 1-2 days following landing. Root-mean-square position error and tilt velocity were significantly reduced with vibrotactile feedback. CONCLUSIONS. These preliminary results are consistent with our hypothesis that adaptive changes in vestibular processing corresponds to reduced manual control performance following G-transitions. A simple vibrotactile prosthesis improves the ability to null out tilt motion within a limited range of motion disturbances.

Wood, S. J.↗