ALBERT

Description: Solid propellant boosters for manned space flight

Description: This paper presents viewgraphs on turbulence detection and mitigation technologies in weather accident prevention. The topics include: 1) Organization; 2) Scope of Turbulence Effort; 3) Background; 4) Turbulence Detection and Mitigation Program Metrics; 5) Approach; 6) Turbulence Team Relationships; 7) WBS Structure; 8) Deliverables; 9) TDAM Changes; 10) FY-01 Results/Accomplishments; 11) Out-year Plans; and 12) Element Status.

Description: This paper presents the weather accident prevention project review during the period of June 5, through June 7, 2001. The topics include: 1) Background; 2) Guidance; 3) Plan; 4) System Elements; 5) AWIN System; 6) Market Segments; 7) Technology Development Level; 8) Aviation Safety Program Organization; 9) Partnerships; 10) NASA Facilities; 11) Timeline; 12) AWIN Research Areas; and 13) Cooperative Research with FAA. This paper is in viewgraph form.

Description: This paper presents viewgraphs on the development of the Weather Accident Prevention (WxAP) System architecture and Concept of Operation (CONOPS) activities. The topics include: 1) Background Information on System Architecture/CONOPS Activity; 2) Activity Work in Progress; and 3) Anticipated By-Products.

Description: Svante Arrhenius' research in atmospheric physics extended beyond the recent past and the near future states of the Earth, which today are at the center of sociopolitical attention. His plan encompassed all of the physical phenomena known at the time to relate to the formation and evolution of stars and planets. His two-volume textbook on cosmic physics is a comprehensive synopsis of the field. The inquiry into the possible cause of the ice ages and the theory of selective wavelength filter control led Arrhenius to consider the surface states of the other terrestrial planets, and of the ancient Earth before it had been modified by the emergence of life. The rapid escape of hydrogen and the equilibration with igneous rocks required that carbon in the early atmosphere prevailed mainly in oxidized form as carbon dioxide, together with other photoactive gases exerting a greenhouse effect orders of magnitude larger than in our present atmosphere. This effect, together with the ensuing chemical processes, would have set the conditions for life to evolve on our planet, seeded from spores spreading through an infinite Universe, and propelled, as Arrhenius thought, by stellar radiation pressure.

Description: Laboratory experiments to quantify the global production of NOx (NO + NO2) in the troposphere due to atmospheric lightning and biogenic activity in soil are presented. These laboratory experiments, as well as other studies, suggest that the global production of NOx by lightning probably ranges between 2 and 20 MT(N)y-1 of NO and is strongly dependent on the total energy deposited by lightning, a quantity not well-known. In our laboratory experiments, nitrifying micro-organisms is soil were found to be a significant source of both NO and nitrous oxide (N2O). The measured production ratio of NO to N2O averaged 2-3 for oxygen partial pressures of 0.5-10%. Extrapolating these laboratory measurements to the global scale, which is somewhat risky, suggests that nitrifying micro-organisms in soil may account for as much as 10 MT(N) y-1 of NO. Additional experiments with denitrifying micro-organisms gave an NO to N2O production ratio ranging from 2 to 4 for an oxygen partial pressure of 0.5% and a ratio of less than unity for oxygen partial pressures ranging from 1 to 20%. The production of NO and N2O, normalized with respect to micro-organism number indicates that the production of both NO and N2O by denitrifying micro-organisms is at least an order of magnitude less than production by nitrifying micro-organisms for the micro-organisms studied.

Description: A self-consistent method of determining initial conditions for the model presented by Berner, Lasaga, and Garrels (1983) (henceforth, the BLAG model) is derived, based on the assumption that the CO2 geochemical cycle was in steady state at t = -100 my (million years). This initialization procedure leads to a dissolved magnesium concentration higher than that calculated by Berner, Lasaga, and Garrels and to a low ratio of dissolved calcium to bicarbonate prior to 60 my ago. The latter prediction conflicts with the geologic record of evaporite deposits, which requires that this ratio remain greater than 0.5. The contradiction is probably caused by oversimplifications in the BLAG model, such as the neglect of the cycles of organic carbon and sulfur.

Description: Selection of high thrust launch vehicles

Description: The NASA/Goddard Space Flight Center Scanning Raman Lidar (SRL) had a highly successful deployment at the Department of Energy Cloud and Radiation Testbed (CART) Site in Billings, OK during April, 1994 for the first Intensive Operation Period (IOP) hosted there. During the IOP, the SRL operated from just after sundown to just before sunrise for all declared evenings of operation. The lidar acquired more than 123 hours of data over 15 nights with less than 1 hour of data lost due to minor system malfunction. The SRL acquired data both on the vertical and in scanning mode toward an instrumented 60 m tower during various meteorological conditions such as an intense cold frontal passage on April 15 which is the focus of this presentation.

Description: The first Atmospheric Radiation Measurement (ARM) Remote Cloud Study (RCS) Intensive Operations Period (IOP) was held during April 1994 at the Southern Great Plains (SGP) Cloud and Radiation Testbed (CART) site near Lamont, Oklahoma. This experiment was conducted to evaluate and calibrate state-of-the-art, ground based remote sensing instruments and to use the data acquired by these instruments to validate retrieval algorithms developed under the ARM program. These activities are part of an overall plan to assess general circulation model (GCM) parameterization research. Since radiation processes are one of the key areas included in this parameterization research, measurements of water vapor and aerosols are required because of the important roles these atmospheric constituents play in radiative transfer. Two instruments were deployed during this IOP to measure water vapor and aerosols and study their relationship. The NASA/Goddard Space Flight Center (GSFC) Scanning Raman Lidar (SRL) acquired water vapor and aerosol profile data during 15 nights of operations. The lidar acquired vertical profiles as well as nearly horizontal profiles directed near an instrumented 60 meter tower. Aerosol optical thickness, phase function, size distribution, and integrated water vapor were derived from measurements with a multiband automatic sun and sky scanning radiometer deployed at this site.

Description: The Airborne Visible-Infrared Imaging Spectrometer (AVIRIS) measures the total upwelling spectral radiance from 400 to 2500 nm sampled at 10 nm intervals. The instrument acquires spectral data at an altitude of 20 km above sea level, as images of 11 by up to 100 km at 17x17 meter spatial sampling. We have developed a nonlinear spectral fitting algorithm coupled with a radiative transfer code to derive the total path water vapor from the spectrum, measured for each spatial element in an AVIRIS image. The algorithm compensates for variation in the surface spectral reflectance and atmospheric aerosols. It uses water vapor absorption bands centered at 940 nm, 1040 nm, and 1380 nm. We analyze data sets with water vapor abundances ranging from 1 to 40 perceptible millimeters. In one data set, the total path water vapor varies from 7 to 21 mm over a distance of less than 10 km. We have analyzed a time series of five images acquired at 12 minute intervals; these show spatially heterogeneous changes of advocated water vapor of 25 percent over 1 hour. The algorithm determines water vapor for images with a range of ground covers, including bare rock and soil, sparse to dense vegetation, snow and ice, open water, and clouds. The precision of the water vapor determination approaches one percent. However, the precision is sensitive to the absolute abundance and the absorption strength of the atmospheric water vapor band analyzed. We have evaluated the accuracy of the algorithm by comparing several surface-based determinations of water vapor at the time of the AVIRIS data acquisition. The agreement between the AVIRIS measured water vapor and the in situ surface radiometer and surface interferometer measured water vapor is 5 to 10 percent.

Description: In the DIAL technique, the water vapor concentration profile is determined by analyzing the lidar backscatter signals for laser wavelengths tuned 'on' and 'off' a water vapor absorption line. Desired characteristics of the on-line transmitted laser beam include: pulse energy greater than or equal to 100 mJ, high-resolution tuning capability (uncertainty less than 0.25 pm), good spectral stability (jitter less than 0.5 pm about the mean), and high spectral purity (greater than 99 percent). The off-line laser is generally detuned less than 100 pm away from the water vapor line. Its spectral requirements are much less stringent. In our past research, we developed and demonstrated the airborne DIAL technique for water vapor measurements in the 720-nm spectral region using a system based on an alexandrite laser as the transmitter for the on-line wavelength and a Nd:YAG laser-pumped dye laser for the off-line wavelength. This off-line laser has been replaced by a second alexandrite laser. Diode lasers are used to injection seed both lasers for frequency and linewidth control. This eliminates the need for the two intracavity etalons utilized in our previous alexandrite laser and thereby greatly reduces the risk of optical damage. Consequently, the transmitted pulse energy can be substantially increased, resulting in greater measurement range, higher data density, and increased measurement precision. In this paper, we describe the diode injection seed source, the two alexandrite lasers, and the device used to line lock the on-line seed source to the water vapor absorption feature.

Description: The differential absorption lidar (DIAL) technique was first applied to the remote measurement of atmospheric water vapor profiles from airborne platforms in 1981. The successful interpretation of the lidar profiles relies strongly on an accurate knowledge of specific water vapor absorption line parameters: line strength, pressure broadening coefficient, pressure-induced shift coefficient and the respective temperature-dependence factors. NASA Langley Research Center has developed and is currently testing an autonomous airborne water vapor lidar system: LASE (Lidar Atmospheric Sensing Experiment). This DIAL system uses a Nd:YAG-pumped Ti:Sapphire laser seeded by a diode laser as a lidar transmitter. The tunable diode has been selected to operate in the 813-818 nm wavelength region. This 5-nm spectral interval offers a large distribution of strengths for temperature-insensitive water vapor absorption lines. In support of the LASE project, a series of spectroscopic measurements were conducted for the 16 absorption lines that have been identified for use in the LASE measurements. Prior to this work, the experimental data for this water vapor absorption band were limited - to our knowledge - to the line strengths and to the line positions.

Description: This paper addresses the design considerations and strategies for astrophysical observations as key elements of an international solar system exploration program. Emphasis is placed on the technical and programmatic challenges and opportunities associated with an evolving program of lunar-based astronomy. Both robotic and human tended facilities are discussed ranging from relatively small meter-class transit telescopes to large interferometer and filled-aperture systems.

Description: A method is presented by which measured modes and frequencies from a modal test can be used to determine the location and magnitude of damage in a space struss structure. The damage is located by computing the Euclidean distances between the measured mode shapes and the best achievable eigenvectors. The best achievable eigenvectors are the projection of the measured mode shapes onto the subspace defined by the refined analytical model of the structure and the measured frequencies. Loss of both stiffness and mass properties can be located and quantified. To examine the performance of the method when experimentally measured modes are employed, various damage detection studies using a laboratory eight-bay truss structure were conducted. The method performs well even though the measurement errors inevitably make the damage location more difficult.

Description: Following the project's first major design review, some unresolved technical issues, mainly centered on details of how to integrate Russian hardware into the U.S./international space station, remain. No 'show stoppers' were found in the review. Specific open technical issues are discussed in this article.

Description: A coupled, trajectory-based flowfield and material thermal-response analysis is presented for the European Space Agency proposed Rosetta comet nucleus sample return vehicle. The probe returns to earth along a hyperbolic trajectory with an entry velocity of 16.5 km/s and requires an ablative heat shield on the forebody. Combined radiative and convective ablating flowfield analyses were performed for the significant heating portion of the shallow ballistic entry trajectory. Both quasisteady ablation and fully transient analyses were performed for a heat shield composed of carbon-phenolic ablative material. Quasisteady analysis was performed using the two-dimensional axisymmetric codes RASLE and BLIMPK. Transient computational results were obtained from the one-dimensional ablation/conduction code CMA. Results are presented for heating, temperature, and ablation rate distributions over the probe forebody for various trajectory points. Comparison of transient and quasisteady results indicates that, for the heating pulse encountered by this probe, the quasisteady approach is conservative from the standpoint of predicted surface recession.

Description: Advanced Satellite for Cosmology and Astrophysics (ASCA) is a high-throughput X-ray astronomy observatory which is capable of simultaneous imaging and spectroscopic observations over a wide energy range 0.5-10 keV. The scientific capabilities of ASCA and some aspects related to its operation and observations are briefly described.

Description: We examine the electromagnetic (EM) bias by using retracked TOPEX altimeter data. In contrast to previous studies, we use a parameterization of the EM bias which does not make stringent assumptions about the form of the correction or its global behavior. We find that the most effective single parameter correction uses the altimeter-estimated wind speed but that other parameterizations, using a wave age related parameter of significant wave height, may also significantly reduce the repeat pass variance. The different corrections are compared, and their improvement of the TOPEX height variance is quantified.

Description: Monthly Ku band sigma(sub 0) and significant wave height (SWH) histograms from the NASA altimeter on the TOPEX/POSEIDON satellite are preseneted for January through June 1993 for three latitude bands between +/- 60 degrees. The data are compared to distributions from the Geosat mission for the same months in 1987-1989. Generally, the distributions agree quite well, although there are some seasonal/hemispherical differences. The sigma(sub 0) comparison reveals an overall bias between the two altimeters with the TOPEX sigma(sub 0) higher by about 0.7 dB, which is consistent with algorithm improvements for TOPEX. The SWH distributions show strong hemispherical/seasonal changes. The seasonal/hemispherical differences between TOPEX and Geosat are consistent for SWH and sigma(sub 0). The joint distribution of sigma(sub 0) and SWH is extremely stable friom month to month. The typical SWH is independent of sigma(sub 0) for sigma(sub 0) greater than 11.3 dB. The minimum SWH grows exponentially with wind speed. This joint distribution may be useful for understanding electromagnetic bias in altimeter measurements. Finally, altimeter data are compared to buoy values from 21 overflights of the NASA verification site near Pt. Conception, California. Wave heights agree well with an root mean square (RMS) difference of only 0.2 m. Altimeter sigma(sub 0) values are compared to buoy wind speeds. The results are consistent with the -0.7 dB sigma(sub 0) offset from the histogram comparisons.

Description: To assess the accuracy of the TOPEX altimeter data, we have reprocessed the raw altimeter waveform data using more sophisticated algorithms than those implemented in the altimeter hardware. We discuss systematic contamination of the waveform which we have observed and its effect on very long wavelength errors. We conclude that these systematic errors are responsible for a very long wavelength error whose peak-to-peak magnitude for the Ku band altimeter is of the order of 1 cm. We also examine the ability of retracked data to reduce the repeat pass variance and correct for significant wave height (SWH) and acceleration dependent errors. We find that the ground postprocessing contains SWH dependent biases which depend on the altimeter fine height correction.

Description: The routine ground processing of data from the NASA radar altimeter of TOPEX/POSEIDON includes instrument corrections for the effects of significant wave height and attitude angle changes on the altimeter's estimates of range, backscattered power, and significant wave height. This paper describes how these instrument corrections were generated and how they are applied. Detailed waveform fitting to telemetered waveform samples is use to assess the effectiveness of the corrections. There are several altimeter hardware-caused small waveform departures from the model waveforms and these departures, designated waveform 'features', are described in detailed. A consequence of the waveform features, and their positioning relationship to range rate, is that range data for ground tracks moving toward the equator may differ systematically by about a centimeter compared to range data for ground tracks moving away from the equator. The results and discussion are limited to side A of the redundant altimeter, as only side A has been operated on orbit.

Description: Results of the in-flight calibration and performance evaluation campaign for the TOPEX/POSEIDON microwave radiometer (TMR) are presented. Intercomparisons are made between TMR and various sources of ground truth, including ground-based microwave water vapor radiometers, radiosondes, global climatological models, special sensor microwave imager data over the Amazon rain forest, and models of clear, calm, subpolar ocean regions. After correction for preflight errors in the processing of thermal/vacuum data, relative channel offsets in the open ocean TMR brightness temperatures were noted at the approximately = 1 K level for the three TMR frequencies. Larger absolute offsets of 6-9 K over the rain forest indicated a approximately = 5% gain error in the three channel calibrations. This was corrected by adjusting the antenna pattern correction (APC) algorithm. AS 10% scale error in the TMR path delay estimates, relative to coincident radiosondes, was corrected in part by the APC adjustment and in part by a 5% modification to the value assumed for the 22.235 FGHz water vapor line strength in the path delay retrieval algorithm. After all in-flight corrections to the calibration, TMR global retrieval accuracy for the wet tropospheric range correction is estimated at 1.1 cm root mean square (RMS) with consistent peformance under clear, cloudy, and windy conditions.

Description: The NASA altimeter on board TOPEX/POSEIDON exploits the difference in the delays of the Ku and C band radar pulses to estimate an ionosphere correction to the range measurement. The dependence of the ionosphere correction on ocean and satellite parameters is less than 1 cm. The standard deviation of the 1-s averaged ionosphere correction depends on the height of the ocean waves and ranges from 5 to 14 mm. The accuracy of the ionosphere correction is better than 1 cm at the 1 sigma confidence level. The ionosphere correction should be averaged over 140 km (20 s) along track in order to minimize its noise without sacrificing its accuracy. Ionosphere models must achieve an independent sample spacing of 900 km or less in order to allow a single-frequency altimeter to have an ionosphere correction comparable in accuracy to that of the NASA dual-frequency altimeter.

Description: Estimates of the effectiveness of an altimetric correction, and interpretation of sea level variability as a response to atmospheric forcing, both depend upon assuming that residual errors in altimetric corrections are uncorrelated among themselves and with residual sea level, or knowing the correlations. Not surprisingly, many corrections are highly correlated since they involve atmospheric properties and the ocean surface's response to them. The full corrections (including their geographically varying time mean values), show correlations between electromagnetic bias (mostly the height of wind waves) and either atmospheric pressure or water vapor of -40%, and between atmospheric pressure and water vapor of 28%. In the more commonly used collinear differences (after removal of the geographically varying time mean), atmospheric pressure and wave height show a -30% correlation, atmospheric pressure and water vapor a -10% correlation, both pressure and water vapor a 7% correlation with residual sea level, and a bit surprisingly, ionospheric electron content and wave height a 15% correlation. Only the ocean tide is totally uncorrelated with other corrections or residual sea level. The effectiveness of three ionospheric corrections (TOPEX dual-frequency, a smoothed version of the TOPEX dual-frequency, and Doppler orbitography and radiopositioning integrated by satellite (DORIS) is also evaluated in terms of their reduction in variance of residual sea level. Smooth (90-200 km along-track) versions of the dual-frequency altimeter ionosphere perform best both globally and within 20 deg in latitude from the equator. The noise variance in the 1/s TOPEX inospheric samples is approximately (11 mm) squared, about the same as noise in the DORIS-based correction; however, the latter has its error over scales of order 10(exp 3) km. Within 20 deg of the equator, the DORIS-based correction adds (14 mm) squared to the residual sea level variance.

Description: We report a new measurement of the cosmic-ray isotopic composition of aluminum in the low-energy range form 75 to 206 MeV per nucleon.This measurement was made using the high-energy telescope of the CRS experiment on the Voyager 1 and 2 spacecraft during the time period from 1977 to 1993 with an average solar modulation level about 497 MV, roughly the same as at Earth near sunspot minimum. We obtain approximately 430 Al events of which approximately 35 are Al-26 and 395 are Al-27. The Al isotopes were separated with an average mass resolution sigma of 0.35 amu. Our interpretation of the isotopic composition of cosmic-ray aluminum is based on a standard Leaky-Box model for the interstellar propagation of cosmic-ray nuclei using the latest cross sections of the New Mexico-Saclay collaboration as well as a disk-halo diffusion model. From our observed ratio Al-26/Al-27 of 8.3 +/- 2.4 % we deduce an average interstellar density of about 0.52 (+0.26, -0.2) atoms per cu cm. This density is larger than the value of 0.28 (+0.14, -0.11) atoms per cu cm we found from an analysis of the observed abundance of the longer lived Be-10 made using data from the Voyager detectors over almost the same time interval and using essentially the same propagation program.

Description: The DET/MPS programs model and simulate the Direct Energy Transfer and Multimission Spacecraft Modular Power System in order to aid both in design and in analysis of orbital energy balance. Typically, the DET power system has the solar array directly to the spacecraft bus, and the central building block of MPS is the Standard Power Regulator Unit. DET/MPS allows a minute-by-minute simulation of the power system's performance as it responds to various orbital parameters, focusing its output on solar array output and battery characteristics. While this package is limited in terms of orbital mechanics, it is sufficient to calculate eclipse and solar array data for circular or non-circular orbits. DET/MPS can be adjusted to run one or sequential orbits up to about one week, simulated time. These programs have been used on a variety of Goddard Space Flight Center spacecraft projects. DET/MPS is written in FORTRAN 77 with some VAX-type extensions. Any FORTRAN 77 compiler that includes VAX extensions should be able to compile and run the program with little or no modifications. The compiler must at least support free-form (or tab-delineated) source format and 'do do-while end-do' control structures. DET/MPS is available for three platforms: GSC-13374, for DEC VAX series computers running VMS, is available in DEC VAX Backup format on a 9-track 1600 BPI tape (standard distribution) or TK50 tape cartridge; GSC-13443, for UNIX-based computers, is available on a .25 inch streaming magnetic tape cartridge in UNIX tar format; and GSC-13444, for Macintosh computers running AU/X with either the NKR FORTRAN or AbSoft MacFORTRAN II compilers, is available on a 3.5 inch 800K Macintosh format diskette. Source code and test data are supplied. The UNIX version of DET requires 90K of main memory for execution. DET/MPS was developed in 1990. A/UX and Macintosh are registered trademarks of Apple Computer, Inc. VMS, DEC VAX and TK50 are trademarks of Digital Equipment Corporation. UNIX is a registered trademark of AT&T Bell Laboratories.

Description: The DET/MPS programs model and simulate the Direct Energy Transfer and Multimission Spacecraft Modular Power System in order to aid both in design and in analysis of orbital energy balance. Typically, the DET power system has the solar array directly to the spacecraft bus, and the central building block of MPS is the Standard Power Regulator Unit. DET/MPS allows a minute-by-minute simulation of the power system's performance as it responds to various orbital parameters, focusing its output on solar array output and battery characteristics. While this package is limited in terms of orbital mechanics, it is sufficient to calculate eclipse and solar array data for circular or non-circular orbits. DET/MPS can be adjusted to run one or sequential orbits up to about one week, simulated time. These programs have been used on a variety of Goddard Space Flight Center spacecraft projects. DET/MPS is written in FORTRAN 77 with some VAX-type extensions. Any FORTRAN 77 compiler that includes VAX extensions should be able to compile and run the program with little or no modifications. The compiler must at least support free-form (or tab-delineated) source format and 'do do-while end-do' control structures. DET/MPS is available for three platforms: GSC-13374, for DEC VAX series computers running VMS, is available in DEC VAX Backup format on a 9-track 1600 BPI tape (standard distribution) or TK50 tape cartridge; GSC-13443, for UNIX-based computers, is available on a .25 inch streaming magnetic tape cartridge in UNIX tar format; and GSC-13444, for Macintosh computers running AU/X with either the NKR FORTRAN or AbSoft MacFORTRAN II compilers, is available on a 3.5 inch 800K Macintosh format diskette. Source code and test data are supplied. The UNIX version of DET requires 90K of main memory for execution. DET/MPS was developed in 1990. A/UX and Macintosh are registered trademarks of Apple Computer, Inc. VMS, DEC VAX and TK50 are trademarks of Digital Equipment Corporation. UNIX is a registered trademark of AT&T Bell Laboratories.

Description: The Flexible Spacecraft Dynamics and Control program (FSD) was developed to aid in the simulation of a large class of flexible and rigid spacecraft. FSD is extremely versatile and can be used in attitude dynamics and control analysis as well as in-orbit support of deployment and control of spacecraft. FSD has been used to analyze the in-orbit attitude performance and antenna deployment of the RAE and IMP class satellites, and the HAWKEYE, SCATHA, EXOS-B, and Dynamics Explorer flight programs. FSD is applicable to inertially-oriented spinning, earth oriented, or gravity gradient stabilized spacecraft. The spacecraft flexibility is treated in a continuous manner (instead of finite element) by employing a series of shape functions for the flexible elements. Torsion, bending, and three flexible modes can be simulated for every flexible element. FSD can handle up to ten tubular elements in an arbitrary orientation. FSD is appropriate for studies involving the active control of pointed instruments, with options for digital PID (proportional, integral, derivative) error feedback controllers and control actuators such as thrusters and momentum wheels. The input to FSD is in four parts: 1) Orbit Construction FSD calculates a Keplerian orbit with environmental effects such as drag, magnetic torque, solar pressure, thermal effects, and thruster adjustments; or the user can supply a GTDS format orbit tape for a particular satellite/time-span; 2) Control words - for options such as gravity gradient effects, control torques, and integration ranges; 3) Mathematical descriptions of spacecraft, appendages, and control systems- including element geometry, properties, attitudes, libration damping, tip mass inertia, thermal expansion, magnetic tracking, and gimbal simulation options; and 4) Desired state variables to output, i.e., geometries, bending moments, fast Fourier transform plots, gimbal rotation, filter vectors, etc. All FSD input is of free format, namelist construction. FSD is written in FORTRAN 77, PASCAL, and MACRO assembler for batch execution and has been implemented on a DEC VAX series computer operating under VMS. The PASCAL and MACRO routines (in addition to the FORTRAN program) are supplied as both source and object code, so the PASCAL compiler is not required for implementation. This program was last updated in 1985.

Description: The Short Orbital Flux Integration Program, SOFIP, was developed to provide a fast and accurate evaluation of the space radiation environment expected to be encountered by geocentric satellites. SOFIP uses Vette's standardized models of the terrestrial trapped particle environment in performing the evaluation. For a given trajectory, the basic analysis involves the calculation of the composite integral orbit spectrum of either protons or electrons. Additional calculations which may be performed include exposure index, peaks per orbit, percent time in electron trapping zones, differential spectrum, and solor proton fluences. The results of the SOFIP calculations should provide the analyst with a detailed picture of the radiation environment in which the vehicle will operate. This program is written in FORTRAN IV for batch execution and has been implemented on an IBM 360 computer with a central memory requirement of approximately 60K of 8 bit bytes. The SOFIP program was developed in 1979.

Description: In an effort to place payloads into orbit at the lowest possible costs, the use of air-breathing space-planes, which reduces the need to carry the propulsion system oxidizer, has been examined. As this approach would require the space-plane to fly at hypersonic speeds for periods of time much greater than that required by rockets, many factors must be considered when analyzing its benefits. The Basic Hypersonic Data and Equations spreadsheet provides data gained from three analyses of a space-plane's performance. The equations used to perform the analyses are derived from Newton's second law of physics (i.e. force equals mass times acceleration); the derivation is included. The first analysis is a parametric study of some basic factors affecting the ability of a space-plane to reach orbit. This step calculates the fraction of fuel mass to the total mass of the space-plane at takeoff. The user is able to vary the altitude, the heating value of the fuel, the orbital gravity, and orbital velocity. The second analysis calculates the thickness of a spherical fuel tank, while assuming all of the mass of the vehicle went into the tank's shell. This provides a first order analysis of how much material results from a design where the fuel represents a large portion of the total vehicle mass. In this step, the user is allowed to vary the values for gross weight, material density, and fuel density. The third analysis produces a ratio of gallons of fuel per total mass for various aircraft. It shows that the volume of fuel required by the space-plane relative to the total mass is much larger for a liquid hydrogen space-plane than any other vehicle made. This program is a spreadsheet for use on Macintosh series computers running Microsoft Excel 3.0. The standard distribution medium for this package is a 3.5 inch 800K Macintosh format diskette. Documentation is included in the price of the program. Macintosh is a registered trademark of Apple Computer, Inc. Microsoft is a registered trademark of Microsoft Corporation.

Description: The DET/MPS programs model and simulate the Direct Energy Transfer and Multimission Spacecraft Modular Power System in order to aid both in design and in analysis of orbital energy balance. Typically, the DET power system has the solar array directly to the spacecraft bus, and the central building block of MPS is the Standard Power Regulator Unit. DET/MPS allows a minute-by-minute simulation of the power system's performance as it responds to various orbital parameters, focusing its output on solar array output and battery characteristics. While this package is limited in terms of orbital mechanics, it is sufficient to calculate eclipse and solar array data for circular or non-circular orbits. DET/MPS can be adjusted to run one or sequential orbits up to about one week, simulated time. These programs have been used on a variety of Goddard Space Flight Center spacecraft projects. DET/MPS is written in FORTRAN 77 with some VAX-type extensions. Any FORTRAN 77 compiler that includes VAX extensions should be able to compile and run the program with little or no modifications. The compiler must at least support free-form (or tab-delineated) source format and 'do do-while end-do' control structures. DET/MPS is available for three platforms: GSC-13374, for DEC VAX series computers running VMS, is available in DEC VAX Backup format on a 9-track 1600 BPI tape (standard distribution) or TK50 tape cartridge; GSC-13443, for UNIX-based computers, is available on a .25 inch streaming magnetic tape cartridge in UNIX tar format; and GSC-13444, for Macintosh computers running AU/X with either the NKR FORTRAN or AbSoft MacFORTRAN II compilers, is available on a 3.5 inch 800K Macintosh format diskette. Source code and test data are supplied. The UNIX version of DET requires 90K of main memory for execution. DET/MPS was developed in 1990. A/UX and Macintosh are registered trademarks of Apple Computer, Inc. VMS, DEC VAX and TK50 are trademarks of Digital Equipment Corporation. UNIX is a registered trademark of AT&T Bell Laboratories.

Description: Anthropogenic aerosols are intricately linked to the climate system and to the hydrologic cycle. The net effect of aerosols is to cool the climate system by reflecting sunlight. Depending on their composition, aerosols can also absorb sunlight in the atmosphere, further cooling the surface but warming the atmosphere in the process. These effects of aerosols on the temperature profile, along with the role of aerosols as cloud condensation nuclei, impact the hydrologic cycle, through changes in cloud cover, cloud properties and precipitation. Unravelling these feedbacks is particularly difficult because aerosols take a multitude of shapes and forms, ranging from desert dust to urban pollution, and because aerosol concentrations vary strongly over time and space. To accurately study aerosol distribution and composition therefore requires continuous observations from satellites, networks of ground-based instruments and dedicated field experiments. Increases in aerosol concentration and changes in their composition, driven by industrialization and an expanding population, may adversely affect the Earth's climate and water supply.

Description: A laser spectrometer based on difference-frequency generation in periodically poled LiNbO3 (PPLN) has been used to quantify atmospheric formaldehyde with a detection limit of 0.32 parts per billion in a given volume (ppbV) using specifically developed data-processing techniques. With state-of-the-art fiber-coupled diode-laser pump sources at 1083 nm and 1561 nm, difference-frequency radiation has been generated in the 3.53-micrometers (2832-cm-1) spectral region. Formaldehyde in ambient air in the 1- to 10-ppb V range has been detected continuously for nine and five days at two separate field sites in the Greater Houston area operated by the Texas Natural Resource Conservation Commission (TNRCC) and the Houston Regional Monitoring Corporation (HRM). The acquired spectroscopic data are compared with results obtained by a well-established wet-chemical o-(2,3,4,5,6-pentafluorobenzyl) hydroxylamine (PFBHA) technique.

Description: Pulsed thermoelectrically cooled QC-DFB lasers operating at 15.6 micrometers were characterized for spectroscopic gas sensing applications. A new method for wavelength scanning based on repetition rate modulation was developed. A non-wavelength-selective pyroelectric detector was incorporated in the sensor configuration giving the advantage of room-temperature operation and low cost. Absorption lines of CO2 and H2O were observed in ambient air, providing information about the concentration of these species.

Description: Antenna temperatures and the corresponding geolocation data from the five sources of the Special Sensor Microwave/Imager data from the Defense Meteorological Satellite Program F11 satellite have been characterized. Data from the Fleet Numerical Meteorology and Oceanography Center (FNMOC) have been compared with data from other sources to define and document the differences resulting from different processing systems. While all sources used similar methods to calculate antenna temperatures, different calibration averaging techniques and other processing methods yielded temperature differences. Analyses of the geolocation data identified perturbations in the FNMOC and National Environmental Satellite, Data and Information Service data. The effects of the temperature differences were examined by generating rain rates using the Goddard Scattering Algorithm. Differences in the geophysical precipitation products are directly attributable to antenna temperature differences.

Description: An overview of galactic cosmic ray (GCR) interaction and transport methods, as implemented in the Langley Research Center GCR transport code, is presented. Representative results for solar minimum, exo-magnetospheric GCR dose equivalents in water are presented on a component by component basis for various thicknesses of aluminum shielding. The impact of proposed changes to the currently used quality factors on exposure estimates and shielding requirements are quantified. Using the cellular track model of Katz, estimates of relative biological effectiveness (RBE) for the mixed GCR radiation fields are also made.

Description: NASA-JPL's four-frequency telecommunication system design entails the creation and integration of a frequency-selective surface (FSS) subreflector into the high-gain antenna subsystem. The FSS design, which incorporates a periodic array of conducting elements on a kevlar/polymer composite structure, will be able to multiplex S, X, Ku, and Ka frequency-band wavelengths. Accounts are presented of the FSS's development, mechanical testing, and electrical testing.

Description: Technology is being developed to process signals from distributed sensors using distributed computations. These distributed sensors provide a new feedback capability for vibration control that has not been exploited. Additionally, the sensors proposed are of an optical and distributed nature and could be employed with known techniques of distributed optical computation (Fourier optics, etc.) to accomplish the control system functions of filtering and regulation in a distributed computer. This paper reviews a procedure for the analytic design of control systems for this application. For illustration, the procedure is applied to the problem of suppressing the vibrations of a simply supported beam. A simulator has been developed to study the effects of sensor and processing errors. An extensive study of the effects of these errors on estimation and regulation performance is presented.

Description: The Task Committee on Methods for Identification of Large Structures in Space was founded in Jul. 1984. The charter of the committee was to prepare a state-of-the-art report on methods of system identification applicable to large space structures (LSS). Funding to support preparation of the report was received in Aug. 1985 from the Air Force Rocket Propulsion Laboratory (now the Air Force Astronautics Laboratory), in the form of a contract to the ASCE. The report was completed, and published by AFRPL in Sep. 1986. The Task Committee consisted of ten members, including ASCE and AFRPL representatives. The membership represented Government, Industry, and Universities, and consisted of electrical, mechanical, and civil engineers, with backgrounds in Structural Dynamics, Optimization, and Controls. An effort was made to use consistent terminology and notation throughout the report which would be compatible with the terminology used in both the structures and controls communities.

Description: This paper describes an autonomous control concept for pointing and articulation of science instruments on the Eos (Earth observing system) NASA/NOAA platforms intended to be operational by the late 1990s. Key features of this concept include advanced control adaptation and tuning strategies which provide performance robustness over a wide range of system uncertainties and mission time criticality. System identification-control modification paradigms are synthesized to form an adaptation continuum over this extended regime of autonomous operations.

Description: Large space structures are characterized by a large number of modes, grouped frequencies, and small inherent damping. Model reduction techniques in time domain may not be effective due to small damping. The model truncation method is generally used. This method can not solve the problem of grouped frequencies, and will lose all the information about the higher order modes. A new method developed in this paper, which tries to minimize the error of interested transfer functions, makes use of all the information of the original system, and achieves improvement not only from a smaller error of transfer functions but also from better frequency distribution.

Description: A space flight experiment being developed by NASA's Office of Aeronautics and Space Technology (OAST) that uses the Space Station as a testbed to study techniques for determining the dynamic characteristics of large space structures (LSS) is described. The experiment is separate from the Space Station Program itself with research objectives outside the domain of Space Station Program objectives. A brief description of the experiment, in terms of the general objective and approach, is given along with a statement of the potential benefits to NASA and others. The bulk of material to follow deals with the experiment definition activity that is underway. The scope of an 'initial' definition study and preliminary results from supporting Space Station dynamics analyses is presented. The term initial is used to indicate that the study currently being conducted has limited objectives and is not expected to complete the required experiment definition. A follow on study is planned and is mentioned in the summary.

Description: A reanalysis of the evolution configuration module pattern for Space Station Freedom is presented. The module pattern is to consist of one habitation module, one airlock, one cupola, two Space Shuttle Orbiter pressurized docking adapters, one assured crew return vehicle, and one pressurized logistics module. Three laboratory modules are also included. A phased plan has been developed on the basis of utilization analysis being performed and derived through actual and projected user mission requirements for the growth of pressurized elements. The evolution requirement impacts to the module pattern are outlined. Application of the growth philosophy and the physical limitations resulted in a number of core element ground rules for module pattern growth. Nodes and modules are to be grown symmetrically in order to optimize poor flight control characteristics and maximize the capability for crew dual egress between elements. Core module pattern elements are not to be grown along Y-Y axis from nodes 1 or 2 due to physical clearance problems with the thermal control system radiators operational envelope.

Description: We conduct a Monte Carlo simulation of the cosmic microwave background (CMB) anisotropy in the UCSB South Pole 1991 degree-scale experiment. We examine cold dark matter cosmology with large-scale structure seeded by the Harrison-Zel'dovich hierarchy of Gaussian-distributed primordial inhomogeneities normalized to the COBE-DMR measurement of large-angle CMB anisotropy. We find it statistically implausible (in the sense of low cumulative probability F lower than 5 percent, of not measuring a cosmological delta-T/T signal) that the degree-scale cosmological CMB anisotropy predicted in such models could have escaped a detection at the level of sensitivity achieved in the South Pole 1991 experiment.

Description: Several gamma-ray bursts in the BATSE data have sufficiently long durations and complex temporal structures with pulses that appear to be spaced quasi-periodically. In order to test and quantify these periods we have applied fast Fourier transformations (FFT) to all these events. We have also performed cross spectral analyses of the FFT of the two extreme (high-low) energy bands in each case to determine the lead/lag of the pulses in different energies.

Description: A cursory examination of cosmic gamma-ray burst time profiles indicates an inhomogeneous distribution of structure. In the first approximation, there seem to be two types of profiles; smooth ones with little structure and highly variable ones with lots of structure. To put this observation to the test, we have examined the statistical nature of the profile derivative to choose which parameter might best be called the burst 'spikiness'. We have found that a good estimator is given by a count of the number of 'spikes' (defined by a specific numerical recipe) and not by the rms deviations from either a pre-burst background or any type of moving average background. The application of this parameter to 30 burst time histories shows it to be consistent over a wide range of profile types. The analysis also reveals a preferred average time between spikes of approximately 1.5 seconds.

Description: Most cosmic gamma-ray burst temporal profiles appear to be comprised of several individual pulses, many of which overlap. It is advantageous to deconvolve the temporal structures into their constituent pulses, and thereby investigate the shape, intensity and temporal distributions of the pulses as a function of energy. Such fundamental pulse descriptors would provide constraints for theoretical modeling of the burst emission process, such as indications of source size, optical depth and geometry, as a function of time. We have developed a deconvolution algorithm which treats sequences of pulse shapes that change deterministically. The algorithm, a generalization of autoregressive techniques, has been applied to a few bright bursts observed by BATSE. Results indicate that, even within short intervals, constituent pulses are not self-similarly shaped, nor do pulse shapes evolve in a simple manner throughout a burst. Hence, the direction of our future work on pulse deconvolution will focus on analysis methods which allow pulse shape to vary.

Description: Resonant Compton upscattering is commended as a mechanism that produces a hard gamma-ray spectrum while suppressing X-rays. This model, however, has severe physical and observational limitations. Effective X-ray suppression places a lower limit on the electron density; above this limit X-rays scatter multiple times, so the single-scattering approximation of this mechanism is invalid. Multiple scattering produces a spectrum that is much harder than the single-scattering spectrum. As the Thomson optical depth of a power-law electron beam approaches unity, photon spawning commences at a high rate and physically invalidates the underlying electron distribution. The Compton upscattering model is therefore only valid over a narrow range of electron densities. An observational consequence of this model is the absence of the third cyclotron resonance. Resonant scattering produces gamma-rays that propagate nearly along the magnetic field. The resonant cross section of the third harmonic, which is strongly angle dependent, falls below the Compton continuum for these gamma rays. The observation of a third cyclotron resonance in a gamma-ray burst spectrum would eliminate resonant Compton scattering as a gamma-ray burst process.

Description: Fifteen KONUS bursts and one HEAO burst showing cyclotron lines have well measured positions. If these are due to strongly magnetized neutron stars then we expect them to be associated with the galactic disk with a scale height of about 500 pc as expected for old neutron stars. The simplest way to test this hypothesis would be to test for a concentration along galactic latitude /b/ = 0. We got a marginal indication of such a concentration, and the probability that the 16 events are isotropically distributed is 0.02. We also show that this result is not in conflict with the present BATSE data set.

Description: Conclusions about the nature of gamma-ray bursts derived from the size-frequency distribution may be altered if a significant correlation exists between burst intensity and spectral shape. Moreover, if gamma-ray bursts have a cosmological origin, such a correlation may be expected to result from the expansion of the universe. We have performed a rudimentary search of the BATSE bursts for hardness/intensity correlations. The range of spectral shapes was determined for each burst by computing the ratio of the intensity in the range 100-300 keV to that in 55-300 keV. We find weak evidence for the existence of a correlation, the strongest effect being present when comparing the maximum hardness ratio for each burst with its maximum rate.

Description: We present preliminary analyses of gamma-ray burst spectra from the BATSE Spectroscopy Detectors. Our conclusions are: (1) No spectral lines have yet been detected in BATSE data from any cosmic gamma-ray burst. This is not surprising as the data for few bright bursts is available, and previous experiments saw lines in only a small fraction of the bursts. (2) Burst spectra show emission up to 20 MeV, with four of eight examined bursts having significant spectral breaks from 1 to 2 MeV. These breaks are consistent with opacity effects due to the interaction of photons with a high magnetic field. (3) Various distance independent parameters from burst spectra and time histories have no correlation with parameters related to distance. In other words, bright bursts look the same as faint bursts. This places a strong constraint on two population models of bursts. (4) The detection time of individual photons is not correlated from detector to detector, as predicted by Mitrofanov's pulsed emission model.

Description: The GRB angular/intensity distributions observed by BATSE greatly constrain models of the spatial and luminosity source distributions. Single populations of Galactic Disk and Galactic Halo sources appear to be untenable, and only a limited subset of models made up of sources in a symmetric Galactic Corona satisfy the observations. Comments are made on other classes of models.

Description: The Energetic Gamma Ray Experiment Telescope (EGRET) on the Compton Gamma Ray Observatory has observed energetic gamma ray bursts and flares. On May 3, 1991, EGRET detected a gamma ray burst both in the energy measuring NaI (Tl) scintillator and independently in the spark chamber imaging assembly. The NaI spectra were accumulated by a special BURST mode of EGRET. The spectra were measured over a range from 1 to 200 MeV, in three sequential spectra of 1,2, and 4 seconds. During the peak of the burst, six individual gamma rays were detected in the spark chamber, allowing a determination of the burst arrival direction. The intense flares of June were also detected. A solar flare on June 4 was observed to last for several minutes and for a brief time, less than a minute, had significant emission of gamma rays exceeding 150 MeV.

Description: Gamma-ray bursts are being detected with unprecedented sensitivity by the Burst and Transient Source Experiment (BATSE) on the Compton Gamma Ray Observatory since its launch in April, 1991. The experiment is detecting about one gamma-ray burst per day. A brief description is presented of the on-orbit performance of BATSE, the methods of identification of bursts, and examples of the diverse time profiles of the gamma-ray bursts observed. The most significant finding thus far is the apparent isotropy of the bursts together with the observed inhomogeneity of the sources.

Description: Using the radio through hard X-ray images of the Crab nebula to derive the spatial dependence of the electron spectrum and the magnetic field distribution from MHD flow models, we have rederived the high-energy spectrum of inverse Compton scattered gamma rays. We find agreement with the observed spectrum at TeV energies, but it is clear that the inverse Compton flux does not contribute significantly to the unpulsed nebular emission observed by COS-B from 50 to 500 MeV, which is consistent with a smooth continuation of the spectrum in hard X-rays. The emission at these energies must therefore be due to synchrotron radiation by electrons of at least PeV energies. It appears that the emission in the high-energy gamma-ray range, sensitive to the highest energy electrons in the nebula, can put interesting constraints on the acceleration mechanism.

Description: The theory of magnetospheric recirculation for cyclic electron energization and transport may apply in part to heliospheric transport of Jovian electrons if enhanced cross-IMF propagation occurs at heliospheric altitudes near and below the solar wind transition region. Low altitude, ecliptic-to-polar transport would short-circuit conventional interplanetary diffusion, facilitate rapid access to the polar heliosphere with minimal adiabatic energy losses, and provide a seed population for acceleration to 100-1000 MeV energies at the solar wind termination shock and in the heliomagnetotail.

Description: Previous applications of the Monte Carlo technique at the quasiparallel earth bow shock has motivated the extension of this technique to oblique shock geometries typical of those found in most astrophysical shock environments. In addition, such a generalization will permit the thorough examination of theoretical predictions of rapid acceleration times at quasi-perpendicular shocks. Therefore, we have embarked on the modification of our existing Monte Carlo code and in this paper outline the major technical aspects involved in developing a simulation of cosmic-ray acceleration at modified oblique nonrelativistic shocks.

Description: Since the trajectories of astrophysical charged particles are bent by magnetic fields and normally curl many times before their detection, their origin may not be inferred from their directions as is the case with photons. Fortunately, charged particles reveal their presence through interactions in many instances leading to high-energy gamma rays. Bremsstrahlung, Compton, synchrotron, and curvature radiation all generally have a monotonically decreasing energy spectra reflecting that of the parent particles, whereas nucleon-nucleon radiation has a maximum at about 70 MeV reflecting the nature of the interaction process. Gamma radiation has been seen coming from neutron stars in pulses with the same period as the radio pulsar. Solar gamma rays also have been seen, as have short bursts of gamma rays whose origin remain a mystery. Galactic diffuse gamma radiation reveals the distribution of cosmic rays in our Galaxy. Beyond our Galaxy, active galaxies are seen in gamma rays implying a huge energy in the form of cosmic rays to be present there. The Optically Violent Variable quasar 3C 279 is particularly astounding. During an active state 3C 279 was seen to be emitting approximately 10 exp 48 erg/s if its radiation is isotropic.

Description: The Burst and Transient Spectroscopy Experiment (BATSE) on the Compton Gamma-Ray Observatory has a powerful capability to provide nearly uninterrupted monitoring in the 25 keV-10 MeV range of both AGN and Galactic black hole candidates such as Cygnus X-1, using the occultation of cosmic sources by the Earth. Progress in background modeling indicates that the data accept region, or fit window tau, around the occultation step can be substantially increased over that conservatively assumed in earlier estimates of BATSE's Earth occultation sensitivity. We show samples of large-tau fits to background and source edges. As a result we expect to be able to perform long-term monitoring of Cygnus X-1 and many of the brighter AGN for the duration of the CGRO mission.

Description: Results are presented of a thermal design optimization study of the segmented GFRP primary reflector of the earth-orbiting Submillimeter Imager and Line Survey telescope. The paper examines the thermal requirements of the primary reflector and the thermal environment of the telescope and describes the thermal design of the primary reflector. Particular attention is given to the geometric math model and the thermal math model of the telescope. A summary for the steady-state thermal performance of the optimized design is presented, showing that the optimized design has reduced, by an order of magnitude, structural spatial temperature gradients, which were earlier shown to be the most significant obstacle in maintaining the required telescope figure accuracy.

Description: The research experience of SCOLE developed by the NASA Langley Research Center is reviewed. Particular attention is given to the advances made in control and modeling techniques, the experimental facilities that are now available, and the continuing needs to advance and to validate the technologies in controlling flexible structures.

Description: An analytic redundancy management approach for on-line component failure detection is described and illustrated using a simulation of the NASA Langley Spacecaft Control Laboratory Experiment (SCOLE) research facility. The SCOLE experimental apparatus simulated is a functional model of the Space Shuttle with a large, flexible, offset-feed antenna cantilevered from the payload bay. This approach uses a single, active, Kalman filter selected from a bank of filters, each element of which was previously designed to accommodate a specific failure condition. The residuals of this active filter are processed through a sequential probability ratio test (SPRT) filter to identify the failure state of the system and, hence, closing the failure accommodation loop, to select the active filter. Results are presented that illustrate the ability of the system to detect and recover from failures introduced in the angular rate and linear acceleration sensors on the SCOLE facility.

Description: One of the basic requirements in engineering analysis is the development of a mathematical model describing the system. Frequently comparisons with test data are used as a measurement of the adequacy of the model. An attempt is typically made to update or improve the model to provide a test verified analysis tool. System identification provides a systematic procedure for accomplishing this task. The terms system identification, parameter estimation, and model correlation all refer to techniques that use test information to update or verify mathematical models. The goal of system identification is to improve the correlation of model predictions with measured test data, and produce accurate, predictive models. For nonmetallic structures the modeling task is often difficult due to uncertainties in the elastic constants. A finite element model of the shell was created, which included uncertain orthotropic elastic constants. A modal survey test was then performed on the shell. The resulting modal data, along with the finite element model of the shell, were used in a Bayes estimation algorithm. This permitted the use of covariance matrices to weight the confidence in the initial parameter values as well as confidence in the measured test data. The estimation procedure also employed the concept of successive linearization to obtain an approximate solution to the original nonlinear estimation problem.

Description: Reliable structural dynamic models will be required as a basis for deriving the reduced-order plant models used in control systems for large space structures. Ground vibration testing and model verification will play an important role in the development of these models; however, fundamental differences between the space environment and earth environment, as well as variations in structural properties due to as-built conditions, will make on-orbit identification essential. The efficiency, and perhaps even the success, of on-orbit identification will depend on having a valid model of the structure. It is envisioned that the identification process will primarily involve parametric methods. Given a correct model, a variety of estimation algorithms may be used to estimate parameter values. This paper explores the effects of modeling errors and model deficiencies on parameter estimation by reviewing previous case histories. The effects depend at least to some extent on the estimation algorithm being used. Bayesian estimation was used in the case histories presented here. It is therefore conceivable that the behavior of an estimation algorithm might be useful in detecting and possibly even diagnosing deficiencies. In practice, the task is complicated by the presence of systematic errors in experimental procedures and data processing and in the use of the estimation procedures themselves.

Description: On-orbit system identification (ID) of large space systems is essential for various reasons. For example, the complex composite structure of such systems cannot be ground-tested; their structural dynamic characteristics must be known accurately in order to accomplish active control. Furthermore, such capability can be used to characterize/identify various disturbances. The identification process is consisted of four principal elements: (1) modeling, (2) the estimation algorithm, (3) input system, and (4) measurement system. These elements are highly correlated and all togerher determine the success of the identification problem. Accurate modeling of large space systems is the most important element of the identification process. Large flexible structures are non-linear and infinite dimensional systems with highly coupled parameters and low frequency packed modes. In addition, these systems are subject to stochastic and time-varying disturbances, they have structural parameters which can vary due to on-orbit assembly deployment, and operations. These systems are generally; however, represented by constant coefficient, finite order differential equations. The non-linearities, coupling and noise effects are also often neglected. Moreover, identification experiment designs which lead to highly complex optimization problems usually require the simultaneous choice of ID algorithm, sensor, and actuator type and placement. On-orbit bandwidth and power restrictions on excitation, limited data window, and restrictions on sensor/actuator type, placement and number, has led to practical questions of implementations.

Description: This paper examines the use of on-orbit identification based on Maximum Likelihood Estimation (MLE) to provide these high-order, high-accuracy control design models for large space structures (LSS's). First, it outlines a general MLE identification algorithm, together with a covariance-analysis procedure to assess algorithm performance in terms of systematic and stochastic errors. Next, it examines various simplifications appropriate for the LSS identification application. Simplified analytical performance results are presented, as are numerical results to support these analyses. Finally, a graphical interpretation of these results is given.

Description: This talk focuses on the determination of state-space models for large space systems using only the output data. The output data could be generated by the unknown or deliberate initial conditions of the space structure in question. We shall review some relevant fundamental work on the state-space modeling of sequential output data that is potentially applicable to large space structures. If formulated in terms of some generalized Markov parameters, this approach is in some sense similar to, but much simpler than, the Juang-Pappa Eigensystem Realization Algorithm (ERA) and the Ho-Kalman construction procedure.

Description: For the future space systems, on-orbit identification (ID) capability will be required to complement on-orbit control, due to the fact that the dynamics of large space structures, spacecrafts, and antennas will not be known sufficiently from ground modeling and testing. The computational requirements for ID of flexible structures such as the space station (SS) or the large deployable reflectors (LDR) are however, extensive due to the large number of modes, sensors, and actuators. For these systems the ID algorithm operations need not be computed in real-time, only in near real-time, or an appropriate mission time. Consequently the space systems will need advanced processors and efficient parallel processing algorithm design and architectures to implement the identification algorithms in near real-time. The MAX computer currently being developed may handle such computational requirements. The purpose is to specify the on-board computational requirements for dynamic and static identification for large space structures. The computational requirements for six ID algorithms are presented in the context of three examples: the JPL/AFAL ground antenna facility, the space station (SS), and the large deployable reflector (LDR).

Description: This proposal discusses a new nonlinear, nonparametric method for off-line modeling and on-line estimation of the deformation of a flexible structure undergoing rapid retargeting maneuvers. In these circumstances, the structural stiffness and damping coefficients depend on the angular acceleration omega(dot), the angular rate omega, and the square of the angular rate omega. In the single axis case, the excitation of the structure is represented by the vector u(exp T) = (omega(dot), omega(exp 2), 2(omega)), to which the structural dynamics responds as a 'bilinear' (i.e., parametrically excited) system. A similar technique for multiaxial rotations yields a bilinear model with respect to matrix valued excitations. Three methods of estimation and modeling are described in this proposal to achieve deformation state determination: (1) a method based on a feedback linearized procedure which gives an estimate of the state by means of observers installed in the deformable body; (2) off-line modeling of the deformation state of the structure by means of topological interpolators; and (3) an on-line structural state estimation method based on a combination of the two previous techniques.

Description: Viewgraphs on current analytical gaps for future needs for large space systems are presented. Topics covered include: future spacecraft; common control objectives; accuracy requirements; increasing complexity; promising approaches; and analytical gaps.

Description: We examine the COBE Differential Microwave Radiometer (DMR) data for evidence of noncosmological source contributions. The DMR maps are cross-correlated with maps of rich clusters, extragalactic IRAS sources, HEAO 1 A-2 X-ray emission, and 5 GHz radio sources. We limit the rms contributions from these sources on a 7 deg angular scale to less than 10 micro-K (95 percent confidence level) in the DMR maps, although the LMC probably contributes about 50 micro-K to a limited region of the sky. Thus, our previous interpretation that the fluctuations in the COBE DMR data are most likely due to cosmic fluctuations at the surface of last scattering remains intact. The Comptonization parameter for hot electrons traced by rich clusters is limited to delta(y) less than 2 x 10 exp -6 (95 percent confidence level) averaged over the 7 deg DMR beam.

Description: After critically reviewing observational results obtained by astronomical spacecraft in the interplanetary medium for several aspects of galactic cosmic rays (GCRs) and anomalous cosmic rays (ACRs), attention is given to spacecraft data gathered in the magnetosphere and a detailed description is given of the Anuradha cosmic-ray experiment carried by Spacelab-3. The Anuradha results discussed concern the orbit average flux and ionization state of ACRs, the origins of partially ionized galactic cosmic-ray sub-Fe and Fe ions, and the significance of enhanced abundance ratios of sub-Fe and Fe ions in GCRs inside the magnetosphere.

Description: Structural sizing and performance data are presented for two different aerobrake hexagonal heatshield panel structural concepts. One concept features a sandwich construction with an aluminum honeycomb core and thin quasi-isotropic graphite-epoxy face sheets. The other concept features a skin-rib isogrid construction with thin quasi-isotropic graphite-epoxy skins, and graphite-epoxy ribs oriented at 0, +60, and -60 deg along the panel. Linear static, linear bifurcation buckling, and nonlinear static analyses were performed to compare the structural performance of the two panel concepts and assess their feasibility for a Lunar Transfer Vehicle aerobrake application.

Description: The Indian cosmic ray experiment Anuradha, conducted onboard Spacelab 3 during April 29-May 6, 1985 was designed to obtain information on the ionization states of low-energy cosmic rays, using the geomagnetic field as a rigidity filter to place an upper limit on the ionization state of individual cosmic ray particles. This paper presents data confirming the presence of three distinct groups of energetic particles in the near-earth space: (1) low-energy (15-25 MeV/nucleon) anomalous cosmic rays that are either singly ionized or consistent with their being in singly ionized state, (2) fully ionized galactic cosmic ray ions, and (3) partially ionized iron and sub-iron group ions (which account for about 20 percent of all the iron and sub-iron group ions detected at the Spacelab 3 orbit within the magnetosphere in the energy interval 25-125 MeV/nucleon). It is argued that these partially ionized heavy ions are indeed a part of the low-energy galactic cosmic rays present in the interplanetary space.

Description: An effort is currently being carried out by the Jet Propulsion Laboratory (JPL) to study mission feasibility and to define functional requirements for various subsystems of the Space Infrared Telescope Facility (SIRTF). As a major part of this effort, structural design requirements have been derived based on the stated mission objectives. Design concerns addressed by these requirements include the limits on mass and location of the center of gravity, launch stiffness and dynamic characteristics, design loads and analysis criteria, survivability of the TITAN IV/Centaur launch environment, thermal control for maintaining a near absolute-zero operating temperature, and helium cryogen volume and storage for a five-year mission. To illustrate how the structural design requirements can be met, a point design of the SIRTF flight hardware system was developed, modeled, and analyzed. A description of the key features of this point design, along with pertinent modeling and analysis results, are discussed in this Paper.

Description: The progress in spacecraft charging is reviewed with particular attention given to the interactions of plasma and penetrating radiation with dielectrics. Topics discussed include the charging environments, elementary charging theory, the anomalies attributed to charging or discharging phenomena, and spacecraft engineering.

Description: The objectives of the Solar Probe mission and the current status of the Solar Probe thermal shield subsystem development are described. In particular, the discussion includes a brief description of the mission concepts, spacecraft configuration and shield concept, material selection criteria, and the required material testing to provide a database to support the development of the shield system.

Description: Neutron stars, relativistic and compact by nature, show great potential for the copious creation of electron-positron pairs in the magnetospheres; these rapidly cool, thermalize, and then annihilate. It is therefore expected that many neutron sources might display evidence of pair annihilation lines in the 400-500 keV range. It is shown that magnetic photon splitting, which operates effectively at these energies and in the enormous neutron star magnetic fields, can destroy an annihilation feature by absorbing line photons and reprocessing them to lower energies. In so doing, photon splitting creates a soft gamma-ray bump and a broad quasi-power-law contribution to the X-ray continuum, which is too flat to conflict with the observed X-ray paucity in gamma-ray bursts. The destruction of the line occurs in neutron stars with surface fields of 5 x 10 exp 12 G or maybe even less, depending on the size of the emission region.

Description: Arcjet thrusters are soon to be used for north/south stationkeeping on commercial communications satellites. A series of tests was performed to evaluate the possible effects of these thrusters on spacecraft charging and the degradation of thermal control material. During the tests the interaction between arcjet plumes and both charged and uncharged surfaces did not cause any significant material degradation. In addition, firing an arcjet thruster benignly reduced the potential of charged surfaces to near zero.

Description: The present multidisciplinary telescope-analysis approach, which encompasses thermal, structural, control and optical considerations, is illustrated for the case of an IR telescope in LEO; attention is given to end-to-end evaluations of the effects of mechanical disturbances and thermal gradients in measures of optical performance. Both geometric ray-tracing and surface-to-surface diffraction approximations are used in the telescope's optical model. Also noted is the role played by NASA-JPL's Integrated Modeling of Advanced Optical Systems computation tool, in view of numerical samples.

Description: In this paper a distributed parameter model for the estimation of modal characteristics of NASA Mini-Mast truss is proposed. A closed-form solution of the Timoshenko beam equation, for a uniform cantilevered beam with two concentrated masses, is derived so that the procedure and the computational effort for the estimation of modal characteristics are improved. A maximum likelihood estimator for the Timoshenko beam model is also developed. The resulting estimates from test data by using Timoshenko beam model are found to be comparable to those derived from other approaches.

Description: Parameter identification and modeling are key elements of a design and operational flight strategy for control of flexible space structures. The emphasis of the identification program is on on-orbit applications to spacecraft control. High performance robust controllers will result from advanced design synthesis techniques and on-orbit identification/system tuning. Near term goals for the program include development of an integrated on-line processing capability for flexible body parameter identification, and validation with physical structure experiments.

Description: Viewgraphs and discussion on on-orbit system parameter identification are included. Topics covered include: dynamic programming filter (DPF); cost function and estimator; frequency domain formulation structrual dynamic identification; and attributes of DPF.

Description: On-going research at The Aerospace Corporation studying the feasibility of applying adaptive control methodologies to the control of flexible space structures is described. A laboratory testbed was established to test system identification and control approaches. The laboratory set-up and controller design approach are discussed. The ARX least squares parameter estimation technique is analyzed in terms of frequency domain transfer function bias error. This analysis approach enables the determination of the effects of sampling rate, sensor type, and data prefiltering on the estimation performance. The ability to identify space structure dynamics over a range of frequencies is shown to be heavily dependent on these factors.

Description: Accurate mathematical models are clearly an important part of the design, modification, control, and damage assessment of large space structures. A critical part of model determination of any type of structure is the use of system identification (SI), a process for using measured excitation and response data to improve the form of a mathematical model or the values of the parameters in it. The reasons for using SI are manifold, but most of them involve uncertainties in either the form of or the parameters in the mathematical model. Some of the more common uncertainties are the following: the nature of the damping, the characteristics of sloppy or sticky joints, inelastic material properties, and parameters in a reduced Degree-of-Freedom (DOF) model. This paper contains a description of the SI algorithm followed by two illustrations: one in which all parameters in a shear building model are determined using pull back and quick release data, and another in which an equivalent reduced DOF model is obtained for a two story frame.

Description: This presentation will address an approach for using modal residues and roots, typically derived from vibration testing, to solve for the lumped parameters of a finite element model. The uniqueness of the approach is that the mass, stiffness, or damping matrices of a structure do not need to be known or estimated a-priori. Instead just the connectivity of the structure needs to be estimated. In addition, the presentation will demonstrate how models with orders larger than the number of modes obtained by the test program can be solved. The presentation will start by covering the basis of the approach, then illustrate a trivial application of how it works, and finally illustrate its use on a cantilevered beam using only a fraction of the mode to solve the mass, stiffness, and damping matrices.

Description: Virtually all space structure mission requirements include control of vibration, position, and possibly shape. In order to satisfy their requirements it is necessary to know the dynamic characteristics of the structure. Therefore, a means must be developed to continuously determine a dynamic model, or changes in the model, and to appropriately adapt the controls to these changes. The detection and location of physical damage would also be a very beneficial characteristic of such a procedure. A discussion outlining the methods to perform On-Orbit Model Determination is presented.

Description: For the past several years much effort has been given to the development of techniques for designing control systems for large space structures (LSS's). The main objective of these efforts has been to develop a LSS control methodology that produces designs that meet strenuous performance requirements and are robust to model inaccuracies. Unfortunately, performance and robustness are conflicting requirements. Because LSS's can not be fully tested on ground, it has become an accepted fact that the design of LSS control systems to meet performance requirements can not be completed until the LSS is placed on-orbit and tested and an accurate model is extracted from on-orbit test results. Modern MIMO sampled-data frequency response design techniques are viable candidates for designing LSS control systems. First, this paper presents techniques for performing MIMO system identification (ID) from test data. Then, techniques for improving the performance of the system ID process in the presence of noise are presented. Finally, practical utility of the system ID approaches are validated by the presentation of results obtained from application on the LSS Ground Test Facility at Marshall Space Flight Center.

Description: The Galileo scan platform is controlled in two degrees of freedom. A clock (Spin Bearing) actuator controls the relative position between the rotor and stator, and a cone actuator controls the position between the stator and the platform. Instruments on the platform are required to point to 140 micro-rad accuracy and 50 micro-rad per second stability. The system identification objectives were to identify dominant structural resonance frequencies, mode shapes, and damping ratio which exist in the transfer function between the clock actuator and the gyro sensor; and position the notch filter to limit undesirable actuator torque output to ensure stability and performance.

Description: Many modern spacecraft are complex multi-body dynamic systems where the bodies are connected by several active control systems for pointing and isolation. Mission requirements indicate that many structural modes and possibly some nonlinear effects will require characterization. Thus, even characterization at the subsystem level will become more difficult than usual. System level characterization difficulties will be compounded by the fact that only limited ground testing will be possible on the full up system, and flight testing will be restricted by an extremely limited measurements set. The object of the present discussion is the application of matrix majorant theory to the problem of assessing dynamic system performance when knowledge of the system is uncertain. We show how majorants provide an effective tool to relate required performance output to system identification test quality in terms of residual uncertainty in input-output relations, parameter values, nonlinearities, and interactions. The underlying machinery consists of the block-norm matrix which is a nonnegative matrix each of whose elements is the norm of a block of a suitably partitioned matrix. A matrix which bounds the block-norm matrix in the sense of nonnegative matrices, i.e., element by element is known as a majorant.

Description: An optimal on-orbit experiment is designed to extract the most information from an on-orbit test, subject to the constraints of the testing environment. However, simply jumping in and optimizing standard measures of information with respect to the experiment design can cause severe problems if attention is not paid to the specific needs and properties of the problem at hand. The actual criteria to be optimized depends on (among other things) the particular ID algorithm and parametrization being used. Two parametric techniques are the focus of this presentation: recursive prediction error method (RPEM) and maximum likelihood estimation (MLE).

Description: The research project which resulted in the AMI (Analytical Model Improvement) method was funded by the Langley Research Center from 1979 through 1985. The objective was to develop a method for obtaining improved dynamic models of structures by using both test data and numerical analysis. The research was successfully performed and the method was applied to a real structure having a realistic NASTRAN model with over 500 degrees of freedom. The method and its application are briefly discussed as well as other indirect benefits in related technical areas.

Description: The objective of this paper is to summarize and review several investigations on the assessment and control of structural damage in civil engineering. Specifically, the definition of structural damage is discussed. A candidate method for the evaluation of damage is then reviewed and demonstrated. Various ways of implementing passive and active control of civil engineering structures are next summarized. Finally, the possibility of applying expert systems is discussed.

Description: An eigenvector expansion method is utilized to predict eigenvalue and eigenvector derivatives due to geometric reconfiguration of a Gimbalflex fine-pointing/vibration isolation system called SAVI (Space Active Vibration Isolation). The eigenvector expansion method used is a modification of the classical method and allows for rigid body roots. Using the resulting modal derivatives, free-free nonlinear equations of motion are developed with Lagrange's Method. These equations represent a nonlinear plant model to be used in conjunction with a control system transient response simulation.

Description: Two different methods are proposed for identifying the structural properties of large orbiting space structures under ordinary service loads, and for assessing potential damage due to impact or other extreme loadings. It is shown that the behavior of a structure in a weightless environment is nonlinear due to unloaded or lightly loaded connections, an effect which not only complicates structural control, but makes the problem of system identification more difficult than for ground based systems. Both proposed methods assume that the structure is subjected to loads imposed by prescribed self stressing systems sufficient to produce repeatable internal force systems in the structure. The first method is based on statical response and requires a survey of structural displacements produced by the self stressing systems. The displacements do not have to be determined completely (i.e., in three directions at each connection), but more displacement information produces more accurate structural stiffness information. It is anticipated that displacement measurements will be taken using on-board laser measurement devices. The second technique employs dynamic stress wave measurement techniques using on-board loading devices and strain gages to track stress wave propagation through the space structure. This approach, which is in its early stages of development, relies on an analysis of transit times of impulsive stress waves and changes in transit times and wave forms due to changes in structural parameters.