Search NASA⌕ Search

SEARCH · Search NASA

Results for “prediction correction”

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 271 records · Page 15

Voyager Photometry of Saturn's A Ring

Saturn's A Ring samples a wide range of dynamical environments, from the relatively unperturbed, optically thick inner region to the outer part of the ring, which contains numerous density waves. We analyze Voyager images of the A Ring to determine how the reflectivity of different radial regions varies with lighting and viewing geometry. We model our data with a classical radiative transfer code that includes the illumination of the rings by the Sun and Saturn. The particles in the inner and mid-A Ring have Bond albedos near 0.5 and are more backscattering than satellites of comparable albedo. The region outside the Encke Gap becomes progressively less backscattering with increasing radius. Particle properties change abruptly outside the Keeler Gap; particles here have an albedo near 0.6 and a Lambert-like phase function. In contrast with previous suggestions, the abundance of free, submicrometer "dust" is small throughout the entire A Ring; this conclusion holds even in the outermost A Ring, which is strongly perturbed by density waves. Models derived from low-phase data, assuming only macroscopic particles, correctly predict the highphase reflectivity of the outer A Ring and individual strong density waves in the mid-A Ring. However, the inner and mid-A Ring are typically darker at high phase by a factor of two than our models predict. This discrepancy may be due to the reduced multiple scattering from a layer in which the particles are more closely packed. We have also studied the quadrupole azimuthal brightness asymmetry of the A Ring. The asymmetry has a full amplitude of 35% in the mid-A Ring in low-phase Voyager 2 images. We present results on its behavior and possible implications for the structure of the rings. Finally, we compare our results with studies using other data sets to synthesize our current understanding of the nature of the A Ring.

Dones, Luke↗

Evaluation of a vortex-based subgrid stress model using DNS databases

The performance of a SubGrid Stress (SGS) model for Large-Eddy Simulation (LES) developed by Misra k Pullin (1996) is studied for forced and decaying isotropic turbulence on a 32(exp 3) grid. The physical viability of the model assumptions are tested using DNS databases. The results from LES of forced turbulence at Taylor Reynolds number R(sub (lambda)) approximately equals 90 are compared with filtered DNS fields. Probability density functions (pdfs) of the subgrid energy transfer, total dissipation, and the stretch of the subgrid vorticity by the resolved velocity-gradient tensor show reasonable agreement with the DNS data. The model is also tested in LES of decaying isotropic turbulence where it correctly predicts the decay rate and energy spectra measured by Comte-Bellot & Corrsin (1971).

Misra, Ashish↗

Hubble Space Telescope Observations of the Luminous IRAS Source FSC 10214+4724: A Gravitationally Lensed Infrared Quasar

With a redshift of 2.3, the IRAS source FSC 10214+4724 is apparently one of the most luminous objects known in the universe. We present an image of FSC 10214+4724 at 0.8 pm obtained with the Hubble Space Telescope (HST) WFPC2 Planetary Camera. The source appears as an unresolved (less then 0.06) arc 0.7 long, with significant substructure along its length. The center of curvature of the arc is located near an elliptical galaxy 1.18 to the north. An unresolved component 100 times fainter than the arc is clearly detected on the opposite side of this galaxy. The most straightforward interpretation is that FSC 10214+4724 is gravitationally lensed by the foreground elliptical galaxy, with the faint component a counter-image of the IRAS source. The brightness of the arc in the HST image is then magnified by approx. 100, and the intrinsic source diameter is approx. 0.0l (80 pc) at 0.25 microns rest wavelength. The bolometric luminosity is probably amplified by a smaller factor (approx. 30) as a result of the larger extent expected for the source in the far-infrared. A detailed lensing model is presented that reproduces the observed morphology and relative flux of the arc and counterimage and correctly predicts the position angle of the lensing galaxy. The model also predicts reasonable values for the velocity dispersion, mass, and mass-to-light ratio of the lensing galaxy for a wide range of galaxy redshifts. A redshift for the lensing galaxy of -0.9 is consistent with the measured surface brightness profile from the image, as well as with the galaxy's spectral energy distribution. The background lensed source has an intrinsic luminosity approx. 2 x 10(exp 13) L(solar mass) and remains a highly luminous quasar with an extremely large ratio of infrared to optical/ultraviolet luminosity.

Eisenhardt, Peter R.↗

On the Relationship Between Transit Velocity of Interplanetary Shocks and Solar Active Processes

Recently, it was reported that preferential relationships exist between the transit velocity V(sub T) of earthward-directed interplanetary shocks and solar active processes, in particular, eruptive filaments outside active regions (the size of the erupting filament L(sub f)) and solar flares (the value of the X-ray characteristic J). Unfortunately, statistical testing of the proposed associations was not accomplished, nor was the 'geo-effectiveness' of the events adequately described. Reported here are the results of a re-examination of the 21 eruptive filaments (SSC-EF events) and 26 X-ray flares (SSC-F events) that have been associated with storm sudden commencements (SSCs) at Earth. Simple statistical testing refutes the claim that a preferential relationship exists between V(sub T) and L(sub F), while it supports the claim that one exists between V(sub T) and J. More importantly, the inferred relationship between V(sub T) and J is found to be more complicated than previously thought. In particular, it now appears that SSC-F events may be separable into two groups, based on the value of J: a low-J group (J less than 56), in which V(sub T) varies directly with J, and a high-J group (J greater than 56), in which V(sub T) varies inversely with J. As a whole, high-J events are associated with shocks of higher average transit velocity than those of low-J events, and SSC-F events with shocks of higher average transit velocity than those of SSC-EF events. Further, high-J events tend to be of greater X-ray class ( greater than M3), longer duration (greater then 80 min), and are more likely to be associated with type II/IV radio emission (9 of 12) than low-J events. They also tend to occur in magnetically complex (gamma/delta configuration) active regions (10 of 12) that are large in area extent (area greater than 445 millionths of a solar hemisphere) on the day of flaring (9 of 12). Of the 9 solar proton events that affected the Earth's environment that were found to be associated with SSC-F events, six were high-J events. Concerning 'geo-effectiveness', there appears to be no preferential relationship between the value of the J-parameter and the most negative value of the Dst geomagnetic index Dst(min) following the SSC, which is found to usually occur at 6-14 h after SSC onset (18 of 26) and which ranged in value from -1 to -249 (having a median value of about -75). Of the 26 SSC-F events, only 14 can be associated with a Dst(min) less than or equal to -75, and of these only 7 were high-J events. Of the 14 storm-related events (i.e. Dst(min) less than or equal to -75), three have previously been identified as being either 'magnetic clouds' or 'bidirectional flows', both manifestations of earthward-directed coronal mass ejections (CMEs). Superposed epoch analyses of selected solar wind parameters and Dst during the interval of storm-related SSC-F events demonstrate that geoeffective SSC-F events tend to be associated with solar wind flows that are faster, greater in magnetic field strength, and have a rotating field which has a strong southward component shortly after SSC onset, in comparison to SSC-F events that do not have Dst(min) less than or equal to 75. Therefore, it is inferred that geoeffective SSC-F events are probably fast earthward-directed CMEs. Although no single parameter is found that can serve as a predictor of high-skill level for determining the geoeffectiveness of an SSC-F event prior to its occurrence at Earth, one finds that knowledge of the flare's hemispheric location and appearance or lack of appearance of a two-ribbon structure is sufficient to correctly predict the geoeffectiveness of 20 out of 25 of the SSC-F events (80%). Surprisingly, the association or lack of association of metric type II/IV radio emission as a characteristic for determining the geoeffectiveness of the SSC-F events proved unfruitful, as did, to a lesser extent, the duration of the X-ray emission.

Wilson, Robert M.↗

Structure-Property Relationships of Bismaleimides

The purpose of this research was to control and systematically vary the network topology of bismaleimides through cure temperature and chemistry (addition of various coreactants) and subsequently attempt to determine structure-mechanical property relationships. Characterization of the bismaleimide structures by dielectric, rheological, and thermal analyses, and density measurements was subsequently correlated with mechanical properties such as modulus, yield strength, fracture energy, and stress relaxation. The model material used in this investigation was 4,4'-BismaleiMidodIphenyl methane (BMI). BMI was coreacted with either 4,4'-Methylene Dianiline (MDA), o,o'-diallyl bisphenol A (DABA) from Ciba Geigy, or Diamino Diphenyl Sulfone (DDS). Three cure paths were employed: a low- temperature cure of 140 C where chain extension should predominate, a high-temperature cure of 220 C where both chain extension and crosslinking should occur simultaneously, and a low-temperature (140 C) cure followed immediately by a high-temperature (220 C) cure where the chain extension reaction or amine addition precedes BMI homopolymerization or crosslinking. Samples of cured and postcured PMR-15 were also tested to determine the effects of postcuring on the mechanical properties. The low-temperature cure condition of BMI/MDA exhibited the highest modulus values for a given mole fraction of BMI with the modulus decreasing with decreasing concentration of BMI. The higher elastic modulus is the result of steric hindrance by unreacted BMI molecules in the glassy state. The moduli values for the high- and low/high-temperature cure conditions of BMI/MDA decreased as the amount of diamine increased. All the moduli values mimic the yield strength and density trends. For the high-temperature cure condition, the room- temperature modulus remained constant with decreasing mole fraction of BMT for the BMI/DABA and BMI/DDS systems. Postcuring PMR-15 increases the modulus over that of the cured material even though density values of cured and postcured PMR were essentially the same. Preliminary results of a continuous and intermittent stress relaxation experiment for BMI:MDA in a 2:1 molar ratio indicate that crosslinking is occurring when the sample is in the undeformed state. Computer simulation of properties such as density, glass transition temperature, and modulus for the low- temperature cure conditions of BMI/MDA and BMI/DABA were completed. The computer modeling was used to help further understand and confirm the structure characterization results. The simulations correctly predicted the trends of these properties versus mole fraction BMI and were extended to other BMI/diamine systems.

Tenteris-Noebe, Anita D.↗

Validation of Vortex-Lattice Method for Loads on Wings in Lift-Generated Wakes

A study is described that evaluates the accuracy of vortex-lattice methods when they are used to compute the loads induced on aircraft as they encounter lift-generated wakes. The evaluation is accomplished by the use of measurements made in the 80 by 120 ft Wind Tunnel of the lift, rolling moment, and downwash in the wake of three configurations of a model of a subsonic transport aircraft. The downwash measurements are used as input for a vortex-lattice code in order to compute the lift and rolling moment induced on wings that have a span of 0.186, 0.510, or 1.022 times the span of the wake-generating model. Comparison of the computed results with the measured lift and rolling-moment distributions the vortex-lattice method is very reliable as long as the span of the encountering or following wing is less than about 0.2 of the generator span. As the span of the following wing increases above 0.2, the vortex-lattice method continues to correctly predict the trends and nature of the induced loads, but it overpredicts the magnitude of the loads by increasing amounts.

Rossow, Vernon J.↗

Time-Resolved O3 Chemical Chain Reaction Kinetics Via High-Resolution IR Laser Absorption Methods

Excimer laser photolysis in combination with time-resolved IR laser absorption detection of OH radicals has been used to study O3/OH(v = 0)/HO2 chain reaction kinetics at 298 K, (i.e.,(k(sub 1) is OH + 03 yields H02 + 02 and (k(sub 2) is H02 + 03 yields OH + 202). From time-resolved detection of OH radicals with high-resolution near IR laser absorption methods, the chain induction kinetics have been measured at up to an order of magnitude higher ozone concentrations ([03] less than or equal to 10(exp 17) molecules/cu cm) than accessible in previous studies. This greater dynamic range permits the full evolution of the chain induction, propagation, and termination process to be temporally isolated and measured in real time. An exact solution for time-dependent OH evolution under pseudo- first-order chain reaction conditions is presented, which correctly predicts new kinetic signatures not included in previous OH + 03 kinetic analyses. Specifically, the solutions predict an initial exponential loss (chain "induction") of the OH radical to a steady-state level ([OH](sub ss)), with this fast initial decay determined by the sum of both chain rate constants, k(sub ind) = k(sub 1) + k(sub 2). By monitoring the chain induction feature, this sum of the rate constants is determined to be k(sub ind) = 8.4(8) x 10(exp -14) cu cm/molecule/s for room temperature reagents. This is significantly higher than the values currently recommended for use in atmospheric models, but in excellent agreement with previous results from Ravishankara et al.

Kulcke, Axel↗

Long-Wavelength Rupturing Instability in Surface-Tension-Driven Benard Convection

A liquid layer with a free upper surface and heated from below is subject to thermocapillary-induced convective instabilities. We use very thin liquid layers (0.01 cm) to significantly reduce buoyancy effects and simulate Marangoni convection in microgravity. We observe thermocapillary-driven convection in two qualitatively different modes, short-wavelength Benard hexagonal convection cells and a long-wavelength interfacial rupturing mode. We focus on the long-wavelength mode and present experimental observations and theoretical analyses of the long-wavelength instability. Depending on the depths and thermal conductivities of the liquid and the gas above it, the interface can rupture downwards and form a dry spot or rupture upwards and form a high spot. Linear stability theory gives good agreement to the experimental measurements of onset as long as sidewall effects are taken into account. Nonlinear theory correctly predicts the subcritical nature of the bifurcation and the selection between the dry spot and high spots.

Swift, J. B.↗

Experimental and Computational Boundary-Layer Studies in a Supersonic Two-Dimensional Nozzle

Calculations have been carried out on the adiabatic laminar boundary layer developing on the surface of a two-dimensional supersonic nozzle, consisting of contoured nozzle blocks and flat sidewalls. Two- and three-dimensional Navier-Stokes codes, as well as two-dimensional boundary-layer codes have been employed. Thee codes have been adapted to the characteristics of a specific wind tunnel nozzle, so that their numerical results could be directly compared with experimental data obtained in the same nozzle. Such comparisons have been made for the boundary-layer growth on the contoured nozzle, and for the boundary-layer growth, surface streamlines and surface shear on the sidewalls. The three-dimensional Navier-Stokes code was found to be the only one to correctly predict the mean boundary-layer flow on both the sidewalls and the contoured nozzle. Theory and experiment both indicate that the sidewall flow is highly three-dimensional, with non-uniform shear, comer vortices and a boundary layer strongly distorted by cross flows induced by lateral pressure gradients.

Brogan, Torence P.↗

Experimental and Computational Boundary-Layer Studies In a Supersonic Two-Dimensional Nozzle

Experiments and computations were carried out on the adiabatic laminar boundary layer developing along the surfaces of a two-dimensional supersonic nozzle, consisting of upper and lower contoured nozzle blocks and flat sidewalls. Two- and three-dimensional Navier-Stokes codes, as well as two-dimensional boundary-layer codes were employed. These codes were adapted to the characteristics of a specific wind tunnel nozzle, so that their numerical results could be directly compared with experimental data obtained in the same nozzle. Such comparisons were made for the boundary-layer growth on the nozzle contoured surfaces, and for the boundary-layer growth, surface streamlines and surface shear on the sidewalls. The three-dimensional Navier-Stokes code was found to be the only one to correctly predict the mean laminar boundary-layer flow on both the sidewalls and the contoured surfaces. Theory and experiment both indicated that the sidewall flow is highly three-dimensional, with non-uniform shear, corner vortices and a boundary layer strongly distorted by cross flows induced by lateral pressure gradients.

Demetriades, Anthony↗

Hybrid Verification of an Interface for an Automatic Landing

Modern commercial aircraft have extensive automation which helps the pilot by performing computations, obtaining data, and completing procedural tasks. The pilot display must contain enough information so that the pilot can correctly predict the aircraft's behavior, while not overloading the pilot with unnecessary information. Human-automation interaction is currently evaluated through extensive simulation. In this paper, using both hybrid and discrete-event system techniques, we show how one could mathematically verify that an interface contains enough information for the pilot to safely and unambiguously complete a desired maneuver. We first develop a nonlinear, hybrid model for the longitudinal dynamics of a large civil jet aircraft in an autoland/go-around maneuver. We find the largest controlled subset of the aircraft's flight envelope for which we can guarantee both safe landing and safe go-around. We abstract a discrete procedural model using this result, and verify a discrete formulation of the pilot display against it. An interface which fails this verification could result in nondeterministic or unpredictable behavior from the pilot's point of view.

Oishi, Meeko↗

Unfolding of Proteins: Thermal and Mechanical Unfolding

We have employed a Hamiltonian model based on a self-consistent Gaussian appoximation to examine the unfolding process of proteins in external - both mechanical and thermal - force elds. The motivation was to investigate the unfolding pathways of proteins by including only the essence of the important interactions of the native-state topology. Furthermore, if such a model can indeed correctly predict the physics of protein unfolding, it can complement more computationally expensive simulations and theoretical work. The self-consistent Gaussian approximation by Micheletti et al. has been incorporated in our model to make the model mathematically tractable by signi cantly reducing the computational cost. All thermodynamic properties and pair contact probabilities are calculated by simply evaluating the values of a series of Incomplete Gamma functions in an iterative manner. We have compared our results to previous molecular dynamics simulation and experimental data for the mechanical unfolding of the giant muscle protein Titin (1TIT). Our model, especially in light of its simplicity and excellent agreement with experiment and simulation, demonstrates the basic physical elements necessary to capture the mechanism of protein unfolding in an external force field.

Hur, Joe S.↗

Combustion Of Metals In Reduced Gravity And Extraterrestrial Environments

The recent focus of this research project has been to model the combustion of isolated metal droplets and, in particular, to couple the existing theories and formulations of phenomena such as condensation, reaction kinetics, radiation, and surface reactions to formulate a more complete combustion model. A fully transient, one-dimensional (spherical symmetry) numerical model that uses detailed chemical kinetics, multi-component molecular transport mechanisms, condensation kinetics, and gas phase radiation heat transfer was developed. A coagulation model was used to simulate the particulate formation of MgO. The model was used to simulate the combustion of an Mg droplet in pure O2 and CO2. Methanol droplet combustion is considered as a test case for the solution method for both quasi-steady and fully transient simulations. Although some important processes unique to methanol combustion, such as water absorption at the surface, are not included in the model, the results are in sufficient agreement with the published data. Since the major part of the heat released in combustion of Mg, and in combustion of metals in general, is due to the condensation of the metal oxide, it is very important to capture the condensation processes correctly. Using the modified nucleation theory, an Arrhenius type rate expression is derived to calculate the condensation rate of MgO. This expression can be easily included in the CHEMKIN reaction mechanism format. Although very little property data is available for MgO, the condensation rate expression derived using the existing data is able to capture the condensation of MgO. An appropriate choice of the reference temperature to calculate the rate coefficients allows the model to correctly predict the subsequent heat release and hence the flame temperature.

Abbud-Madrid, A.↗

Measurements of the secondary particle energy spectra in the Space Shuttle

Measurements of the energy spectra of secondary particles produced by galactic cosmic rays and trapped protons due to the nuclear interactions of these particles with the Shuttle shielding provide a powerful tool for validating radiation transport codes. A code validated in this way can be used to better estimate the dose and dose equivalent to body organs, measurements that cannot be made directly. The principal cause of single event upsets in electronic devices in the region of the South Atlantic Anomaly is secondary particles, and even in the region of galactic cosmic radiation a significant fraction is produced by secondary particles. In this paper, we describe the first direct measurements of the energy spectra of secondary protons, deuterons, tritons, 3He and 4He produced by galactic cosmic rays inside the Space Shuttle using a charged particle spectrometer. A comparison of these spectra with radiation transport code HZETRN showed reasonably good agreement for secondary protons. However, the code seriously underestimated the flux of all other light ions. The code has been modified to include pick-up and knock-on processes. The modified code leads to good agreement for deuterons and 3He but not for other light ions. This revised code leads to about 10% higher dose equivalent than the original code under moderate shielding, if we assume that higher charge ion fluxes are correctly predicted by the model.

short duration↗

Numerical Modeling of Conjugate Heat Transfer in Fluid Network

Fluid network modeling with conjugate heat transfer has many applications in Aerospace engineering. In modeling unsteady flow with heat transfer, it is important to know the variation of wall temperature in time and space to calculate heat transfer between solid to fluid. Since wall temperature is a function of flow, a coupled analysis of temperature of solid and fluid is necessary. In cryogenic applications, modeling of conjugate heat transfer is of great importance to correctly predict boil-off rate in propellant tanks and chill down of transfer lines. In TFAWS 2003, the present author delivered a paper to describe a general-purpose computer program, GFSSP (Generalized Fluid System Simulation Program). GFSSP calculates flow distribution in complex flow circuit for compressible/incompressible, with or without heat transfer or phase change in all real fluids or mixtures. The flow circuit constitutes of fluid nodes and branches. The mass, energy and specie conservation equations are solved at the nodes where as momentum conservation equations are solved at the branches. The proposed paper describes the extension of GFSSP to model conjugate heat transfer. The network also includes solid nodes and conductors in addition to fluid nodes and branches. The energy conservation equations for solid nodes solves to determine the temperatures of the solid nodes simultaneously with all conservation equations governing fluid flow. The numerical scheme accounts for conduction, convection and radiation heat transfer. The paper will also describe the applications of the code to predict chill down of cryogenic transfer line and boil-off rate of cryogenic propellant storage tank.

Majumdar, Alok↗

Surviving atmospheric spacecraft breakup

Spacecraft travel higher and faster than aircraft, making breakup potentially less survivable. As with aircraft breakup, the dissipation of lethal forces via spacecraft breakup around an organism is likely to greatly increase the odds of survival. By employing a knowledge of space and aviation physiology, comparative physiology, and search-and-rescue techniques, we were able to correctly predict and execute the recovery of live animals following the breakup of the space shuttle Columbia. In this study, we make what is, to our knowledge, the first report of an animal, Caenorhabditis elegans, surviving the atmospheric breakup of the spacecraft that was supporting it and discuss both the lethal events these animals had to escape and the implications for search and rescue following spacecraft breakup.

Disasters↗

Aeroelastic Calculations of Quiet High- Speed Fan Performed

An advanced high-speed fan was recently designed under a cooperative effort between the NASA Glenn Research Center and Honeywell Engines & Systems. The principal design goals were to improve performance and to reduce fan noise at takeoff. Scale models of the Quiet High-Speed Fan were tested for operability, performance, and acoustics. During testing, the fan showed significantly improved noise characteristics, but a self-excited aeroelastic vibration known as flutter was encountered in the operating range. Flutter calculations were carried out for the Quiet High-Speed Fan using a three-dimensional, unsteady aerodynamic, Reynolds-averaged Navier-Stokes turbomachinery code named "TURBO." The TURBO code can accurately model the viscous flow effects that can play an important role in various aeroelastic problems such as flutter with flow separation, flutter at high loading conditions near the stall line (stall flutter), and flutter in the presence of shock and boundary-layer interaction. Initially, calculations were performed with no blade vibrations. These calculations were at a constant rotational speed and a varying mass flow rate. The mass flow rate was varied by changing the backpressure at the exit boundary of the computational domain. These initial steady calculations were followed by aeroelastic calculations in which the blades were prescribed to vibrate harmonically in a natural mode, at a natural frequency, and with a fixed interblade phase angle between adjacent blades. The AE-prep preprocessor was used to interpolate the in-vacuum mode shapes from the structural dynamics mesh onto the computational fluid dynamics mesh and to smoothly propagate the grid deformations from the blade surface to the interior points of the grid. The aeroelastic calculations provided the unsteady aerodynamic forces on the blade surface due to blade vibrations. These forces were vector multiplied with the structural dynamic mode shape to calculate the work done on the blade during one vibration period, then this result was converted to an aerodynamic damping. Flutter occurs when the aerodynamic damping becomes negative, if structural damping is ignored. The results of these aeroelastic calculations are summarized in this plot of aerodynamic damping versus mass flow rate at a constant rotational speed. As the backpressure is increased, the mass flow rate through the fan decreases and the fan operating point moves towards the stall line. The aeroelastic calculations showed that the aerodynamic damping decreases as the stall line is approached, as observed during testing. In addition, the aeroelastic calculations with the TURBO code correctly predicted the aeroelastic parameters: the most unstable vibration mode and interblade phase angle, as observed during testing. The Quiet High-Speed Fan demonstrated significant noise reductions during testing, but flutter imposed limits on its operating range. The accurate calculation of the aeroelastic characteristics using the TURBO code is a significant step toward eliminating flutter from the operating range and toward realizing the benefits of reduced fan noise. The aeroelastic calculations described here were performed under a grant by University of Toledo researchers in collaboration with Glenn's researchers.

Bakhle, Milind A.↗

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