ALBERT

Description: Wind tunnel tests have been conducted on an NACA 2412 airfoil section at Reynolds number of 2.2 x 10(exp 6) and Mach number of 0.13. Detailed measurements of flow fields associated with turbulent boundary layers have been obtained at angles of attack of 12.4 degrees, 14.4 degrees, and 16.4 degrees. Pre- and post-separated velocity and pressure survey results over the airfoil and in the associated wake are presented. Extensive force, pressure, tuft survey, hot-film survey, local skin friction, and boundary layer data are also included. Pressure distributions and separation point locations show good agreement with theory for the two layer angles of attack. Boundary layer displacement thickness, momentum thickness, and shape factor agree well with theory up to the point of separation. There is considerable disparity between extent of flow reversal in the wake as measured by pressure and hot-film probes. The difference is attributed to the intermittent nature of the flow reversal.

Description: We have investigated the interaction of Io, Jupiter's innermost Galilean satellite, with the Io plasma torus. The interaction of Io with the plasma surrounding it has been a subject of interest for almost 30 years, dating from the discovery by Bigg (1964) that radio emissions from the Jovian magnetosphere are controlled by Io's position. Since that time, both ground-based and spacecraft observations have shown that Io is a unique satellite that influences the Jovian magnetosphere in important ways. In particular, material from Io is a major source of plasma for the magnetosphere, and the energy that this plasma harnesses from Jupiter's co-rotating magnetic field is an important power source for the magnetosphere. It is apparent that the local interaction of the torus plasma with Io plays a key role in the formation, composition, and energetics of the Io torus; the interaction is also highly nonlinear. We have modeled this interaction using time-dependent three-dimensional magnetohydrodynamic (MHD) simulations. During this past year, we have used NASA support to develop a new MHD code to study the interaction. As part of the Galileo spacecraft's recent successful insertion into orbit around Jupiter, the spacecraft passed within 900 km of Io's surface. Our calculations have focused on using Galileo particles and fields data to examine a question that was not resolved by the Voyager observations: Does Io have an intrinsic magnetic field? In this progress summary, we describe our efforts on this problem to date.

Description: The implementation of a two-equation k-omega turbulence model into the NPARC flow solver is described. Motivation for the selection of this model is given, major code modifications are outlined, new imputs to the code are described, and results are presented for several validation cases: an incompressible flow over a smooth flat plate, a subsonic diffuser flow, and a shock-induced separated flow. Comparison of results with the k-epsilon model indicate that the k-omega model predicts simple flows equally well whereas, for adverse pressure gradient flows, the k-omega model outperforms the other turbulence models in NPARC.

Description: This paper develops a simplified continuum (continuous with respect to time, stress, etc.) fatigue damage model for use in critical design, Life Extending Control and fault prognosis. The work is based on the local strain cyclic damage modeling method. New nonlinear explicit equation forms of cyclic damage in terms of stress amplitude are derived to facilitate the continuum modelling. Stress based continuum models are derived. Extension to plastic strain-strain rate models is also presented. Progress toward a non-zero mean stress based is presented. Also new nonlinear explicit equation forms in terms of stress amplitude are derived for this case. Application of the various models to design, control, and fault prognosis is considered.

Description: An approach for solving the compressible Euler and Navier-Stokes equations upon meshes composed of nearly arbitrary polyhedra is described. Each polyhedron is constructed from an arbitrary number of triangular and quadrilateral face elements, allowing the unified treatment of tetrahedral, prismatic, pyramidal, and hexahedral cells, as well the general cut cells produced by Cartesian mesh approaches. The basics behind the numerical approach and the resulting data structures are described. The accuracy of the mixed volume grid approach is assessed by performing a grid refinement study upon a series of hexahedral, tetrahedral, prismatic, and Cartesian meshes for an analytic inviscid problem. A series of laminar validation cases are made, comparing the results upon differing grid topologies to each other, to theory, and experimental data. A computation upon a prismatic/tetrahedral mesh is made simulating the laminar flow over a wall/cylinder combination.

Description: The eruption of Mt. Pinatubo in June 1991 produced the largest enhancement of stratospheric aerosol loading ever observed by lidar over Hampton, Virginia. Low altitude layers (less than 20 km) were the first to arrive over Hampton in early August, the result of transport associated with a tropospheric anticyclonic cell over North America. The maximum peak scattering ratio, 34 at 22.4 km, and the maximum stratospheric integrated backscatter of 0.0053 sr(exp -1), both at 694 nm, observed since the eruption were measured on February 20, 1992. After decreasing during the spring and summer of 1992, the aerosol burden increased significantly during the winter of 1992-93, evidence of a poleward winter transport from the equatorial reservoir. Over the period from February 1992 to February 1994, the stratospheric aerosol loading decreased with an average 1/e decay time of 10.1 months. The vertical distribution, intensity, and transport of Pinatubo aerosols over this site are described and compared with similar measurements after El Chichon.

Description: T-matrix computations of light scattering by polydispersions of randomly oriented nonspherical aerosols and Mie computations for equivalent spheres are compared. Findings show that even moderate nonsphericity results in suubstantial errors in the retrieved aerosol optical thickness if satellite reflectance measurements are analyzed using Mie theory. On the other hand, the use of Mie theory for nonspherical aerosols produces negligible errors in the computation of albedo and flux related quantities, provided that the aerosol size distribution and optical thickness are known beforehand. The first result can be explained by large nonspherical-spherical differences in scattering phase function, while the second result follows from small nonspherical-spherical differences in single-scattering albedo and asymmetry parameter. No cancellation of errors occurs if one consistently uses Mie theory in the retrieval algorithm and then in computing the albedo for the retrieved aerosol optical thickness.

Description: Azimuthal asymmetries in the atmospheric refractive index can lead to errors in estimated vertical and horizontal station coordinates. Daily average gradient effects can be as large as 50 mm of delay at a 7 deg elevation. To model gradients, the constrained estimation of gradient paramters was added to the standard VLBI solution procedure. Here the analysis of two sets of data is summarized: the set of all geodetic VLBI experiments from 1990-1993 and a series of 12 state-of-the-art R&D experiments run on consecutive days in January 1994. In both cases, when the gradient parameters are estimated, the overall fit of the geodetic solution is improved at greater than the 99% confidence level. Repeatabilities of baseline lengths ranging up to 11,000 km are improved by 1 to 8 mm in a root-sum-square sense. This varies from about 20% to 40% of the total baseline length scatter without gradient modeling for the 1990-1993 series and 40% to 50% for the January series. Gradients estimated independently for each day as a piecewise linear function are mostly continuous from day to day within their formal uncertainties.

Description: Particles leving the neutral sheet in the distant magnetotail at times display adiabatic trajectory sequences characterized by an inflection toward the equator and subsequent mirroring in its vicinity. We demonstrate that this low-latitude mirroring results primarily from a centrifugal deceleration due to the fast direction-changing E x B drift. This effect which we refer to as 'centrifugal trapping' appears both in guiding centre and full particle treatments. It thus does not directly relate to nonadiabatic motion. However, pitch angle scattering due to nonadiabatic neutral sheet interaction does play a role in reducing the parallel speed of the particles. We show that centrifugal trapping is an important mechanism for the confinement of the slowest (typically below the equatorial E x B drift speed) plasma sheet populations to the midplane vicinity.

Description: Estimates of the mass of dust suspended in the Martian atmosphere are derived from global and regional 9-micrometer opacity maps produced from Viking Infrared Thermal Mapper data. During the peak of the 1977b storm, a total dust mass of approximately 4.3 x 10(exp 14) g was suspended, equivalent to 4.3 x 10(exp -4) g/sq cm, or a layer 1.4 micrometers thick. During a local dust storm near Solis Planum at L(sub s) 227 deg, approximately 1.3 x 10(exp 13) g of dust were lofted, equal to about a 6-micrometer layer in that vicinity.

Description: We model an infrared outburst on Io as being due to a large, erupting lava flow which increased its area at a rate of 1.5 x 10(exp 5)/sq m and cooled from 1225 to 555 K over the 2.583-hr period of observation. The inferred effusion rate of 3 x 10(exp 5) cu m/sec for this eruption is very high, but is not unprece- dented on the Earth and is similar to the high eruption rates suggested for early lunar volcanism. Eruptions occur approxi- mately 6% of the time on Io. These eruptions provide ample resurfacing to explain Io's lack of impact craters. We suggest that the large total radiometric heat flow, 10(exp 14) W, and the size and temperature distribution of the thermal anomalies (McEwen et al. 1992; Veeder et al. 1994) can be accounted for by a series of silicate lava flows in various stages of cooling. We propose that the whole suite of Io's currently observed thermal anomalies was produced by multiple, high-eruptive-rate silicate flows within the past century.

Description: Measurements of wing buffeting, using root strain gages, were made in the NASA Langley 0.3 m cryogenic wind tunnel to refine techniques which will be used in larger cryogenic facilities such as the United States National Transonic Facility (NTF) and the European Transonic Wind Tunnel (ETW). The questions addressed included the relative importance variations in frequency parameter and Reynolds number, the choice of model material (considering both stiffness and damping) and the effects of static aeroelastic distortion. The main series of tests was made on three half models of slender 65 deg delta wings with a sharp leading edge. The three delta wings had the same planform but widely differing bending stiffnesses and frequencies (obtained by varying both the material and the thickness of the wings). It was known that the steady flow on this configuration would be insensitive to variations in Reynolds number. On this wing at vortex breakdown the spectrum of the unsteady excitation is unusual, having a sharp peak at particular frequency parameter. Additional tests were made on one unswept half-wing of aspect ratio 1.5 with an NPL 9510 aerofoil section, known to be sensitive to variations in Reynolds number at transonic speeds. The test Mach numbers were M = 0.21 and 0.35 for the delta wings and to M = 0.30 for the unswept wing. On this wing the unsteady excitation spectrum is fairly flat (as on most wings). Hence correct representation of the frequency parameter is not particularly important.

Description: Measurements of superthermal electron fluxes in the solar wind indicate that field lines within coronal mass ejections, CMEs, near and beyond 1 AU are normally connected to the Sun at both ends. However, on occasion some field lines embedded deep within CMEs appear to be connected to the Sun at only one end. Here we propose an explanation for how such field lines arise in terms of 3-dimensional reconnection close to the Sun. Such reconnection also provides a natural explanation for the flux rope topology characteristic of many CMEs as well as the coronal loops formed during long-duration, solar X-ray events. Our consideration of the field topologies resulting from 3-dimensional reconnection indicates that field lines within and near CMEs may on occasion be connected to the outer heliosphere at both ends.

Description: We investigate the orbital dynamics of small dust particles generated via the continuous micrometeoroid bombardment of the Martian moons. In addition to Mar's oblateness, we also consider the radiation pressure perturbation that is complicated by the planet's eccentric orbit and tilted rotational axis. Considering the production rates and the lifetimes of dust grains, we show that particles from Deimos with radii of about 15 micrometers are expected to dominate the population of a permanently present and tilted dust torus. This torus has an estimated peak number density of approximately equals 5 x 10(exp -12)/cu cm and an optical depth of approximately equals 4 x 10(exp -8).

Description: We examined 11 cases when the interplanetary magnetic field (IMF) was intensely northward (greater than 10 nT) for long durations of time (greater than 3 hours), to quantitatively determine an uppler limit on the efficiency of solar wind energy injection into the magnetosphere. We have specifically selected these large B(sub N) events to minimize the effects of magnetic reconnection. Many of these cases occurred during intervals of high-speed streams associated with coronal mass ejections when viscous interaction effects might be at a maximum. It is found that the typical efficiency of solar wind energy injection into the magnetosphere is 1.0 x 10(exp -3) to 4.0 x 10(exp -3), 100 to 30 times less efficient than during periods of intense southward IMFs. Other energy sinks not included in these numbers are discussed. Estimates of their magnitudes are provided.

Description: The valley network channels on the heavily cratered ancient surface of Mars suggest the presence of liquid water approximately 3.8 Gyr ago. However, the implied warm climate is difficult to explain in the context of the standard solar model, even allowing for the maximum CO2 greenhouse heating. In this paper we investigate the astronomical and planetary implications of a nonstandard solar model in which the zero-age, main-sequence Sun had a mass of 1.05 +/- 0.02 Solar Mass. The excess mass was subsequently lost in a solar wind during the first 1.2(-0.2, +0.4) Gyr of the Sun's main sequence phase. The implied mass-loss rate of 4(+3, -2) x 10(exp -11) M/yr, or about 10(exp 3) x that of the current Sun, may be detectable in several nearby young solar type stars.

Description: Hellas basin on Mars has been the site of volcanism, tectonism, and modification by fluvial, mass-wasting, and eolian processes over its more than 4-b.y. existence. Our detailed geologic mapping and related studies have resulted in the following new interpretations. The asymmetric distribution of highland massifs and other structures that define the uplifted basin rim suggest a formation of the basin by the impact of a low-angle bolide having a trajectory heading S60E. During the Late Noachian, the basin was infilled, perhaps by lava flows, that were sufficiently thick (greater than 1 km) to produce wrinkle ridges on the fill material and extensional faulting along the west rim of the basin. At about the same time, deposits buried northern Malea Planum, which are interpreted to be pyroclastic flows from Amphitrites and Peneus Paterae on the basis of their degraded morphology, topology, and the application of a previous model for pyroclastic volcanism on Mars. Peneus forms a distinctive caldera structure that indicates eruption of massive volumes of magma, whereas Amphitrites is a less distinct circular feature surrounded by a broad, low, dissected shield that suggests generally smaller volume eruptions. During the Early Hesperian, an approximately 1-to 2km-thick sequence of primarily fined-grained, eolian material was deposited on the floor of Hellas basin. Subsequently, the deposit was deeply eroded, except where armored by crater ejecta, and it retreated as much as 200-300 km along its western margin, leaving behind pedestal craters and knobby outliers of the deposit. Local debris flows within the deposit attest to concentrations of groundwater, perhaps in part brought in by outflow floods along the east rim of the basin. These floods may have deposited approximately 100-200m of sediment, subduing wrinkle ridges in the eastern part of the basin floor. During the Late Hesperian and Amazonian, eolian mantles were emplaced on the basin rim and floor and surrounding highlands. Their subsequent erosion resulted in pitted and etched plains and crater fill, irregular mesas, and pedestal craters. Local evidence occurs for the possible former presence of ground ice or ice sheets approximately 100 km across; however, we disagree with a hypothesis that suggest that the entire south rim and much of the floor of Hellas have been glaciated. Orientations of dune fields and yardangs in lower parts of Hellas basin follow directions of the strongest winds predicted by a recently published general circulation model (GCM). Transient frost and dust splotches in the region are, by contrast, related to the GCM prediction for the season in which the images they appear in were taken.

Description: Visible and near-IR refectivity, Moessbauer, and X ray diffraction data were obtained on powders of impact melt rock from the Manicouagan Impact Crater located in Quebec, Canada. The iron mineralogy is dominated by pyroxene for the least oxidized samples and by hematite for the most oxidized samples. Phyllosilicate (smectite) contents up to approximately 15 wt % were found in some heavily oxidized samples. Nanophase hematite and/or paramagnetic ferric iron is observed in all samples. No hydrous ferric oxides (e.g., goethite, lepidocrocite, and ferrihydrite) were detected, which implies the alteration occurred above 250 C. Oxidative alteration is thought to have occurred predominantly during late-stage crystallization and subsolidus cooling of the impact melt by invasion of oxidizing vapors and/or solutions while the impact melt rocks were still hot. The near-IR band minimum correlated with the extent of aleration Fe(3+)/Fe(sub tot) and ranged from approximately 1000 nm (high-Ca pyroxene) to approximately 850 nm (bulk, well-crystalline hematite) for least and most oxidized samples, respectively. Intermediate band positions (900-920 nm) are attributed to low-Ca pyroxene and/or a composite band from hematite-pyroxene assemblages. Manicouagan data are consistent with previous assignments of hematite and pyroxene to the approximately 850 and approximately 1000nm bands observed in Martian reflectivity spectra. Manicouagan data also show that possible assignments for intermediate band positions (900-920 nm) in Martian spectra are pyroxene and/or hematite-pyroxene assemblages. By analogy with impact melt sheets and in agreement with observables for Mars, oxidative alteration of Martian impact melt sheets above 250 C and subsequent erosion could produce rocks and soils with variable proportions of hematite (both bulk and nanophase), pyroxene, and phyllosilicates as iron-bearing mineralogies. If this process is dominant, these phases on Mars were formed rapidly at relativly high temperatures on a sporadic basis throughout the history of the planet. The Manicouagan samples also show that this mineralogical diversity can be accomplished at constant chemical composition, which is also indicated for Mars from the analyses of soil at the two Viking landing sites.

Description: The subject of this paper is a self-consistent, magnetohydrodynamic numerical realization for the Earth's magnetosphere which is in a quasi-steady dynamic equilibrium for a due northward interplanetary magnetic field (IMF). Although a few hours of steady northward IMF are required for this asymptotic state to be set up, it should still be of considerable theoretical interest because it constitutes a 'ground state' for the solar wind-magnetosphere interaction. Moreover, particular features of this ground state magnetosphere should be observable even under less extreme solar wind conditions. Certain characteristics of this magnetosphere, namely, NBZ Birkeland currents, four-cell ionospheric convection, a relatively weak cross-polar potential, and a prominent flow boundary layer, are widely expected. Other characteristics, such as no open tail lobes, no Earth-connected magnetic flux beyond 155 R(sub E) downstream, magnetic merging in a closed topology at the cusps, and a 'tadpole' shaped magnetospheric boundary, might not be expected. In this paper, we will present the evidence for this unusual but interesting magnetospheric equilibrium. We will also discuss our present understanding of this singular state.

Description: The stability of collisionless tearing modes is examined in the presence of curvature drift resonances and the trapped particle effects. A kinetic description for both electrons and ions is employed to investigate the stability of a two-dimensional equilibrium model. The main features of the study are to treat the ion dynamics properly by incorporating effects associated with particle trajectories in the tail fields and to include the linear coupling of trapped particle modes. Generalized dispersion relations are derived in several parameter regimes by considering two important sublayers of the reconnecting region. For a typical choice of parameters appropriate to the current sheet region, we demonstrate that localized tearing modes driven by ion curvature drift resonance effects are excited in the current sheet region with growth time of the order of a few seconds. Also, we examine nonlocal characteristics of tearing modes driven by curvature effects and show that modes growing in a fraction of a second arise when mode widths are larger than the current sheet width. Further, we show that trapped particle effects, in an interesting frequency regime, significantly enhance the growth rate of the tearing mode. The relevance of this theory for substorm onset phase and other features of the substorms is briefly discussed.

Description: The second pulsating aurora rocket (PULSAUR 2) rocket was launched into a pulsating aurora on 9 February 1994 and carried a broad range of instruments in order to perform a study of this type of aurora. The rocket measurements were complemented with a set of ground-based measurements. The particle measurements performed in the rocket are related to the ground-based optical measurements performed along the rocket trajectory. It is found that the high energy electrons are largely in phase with the measured luminosity. The EISCAT measurements carried out during the flight are reviewed. The PULSAUR 2 campaign is described. Results concerning the auroral conditions and the particle measurements are presented.

Description: The high-speed correction factor to the O(+)-O collision frequency, resulting from drift velocities between ions and neutrals, is calculated by solving the integral expression in this factor both numerically and analytically. Although the analytic solution is valid for either small or large drift velocities between ions and neutrals, for temperatures of interest and all drift velocities considered, agreement is found between analytic and detailed numerical integration results within less than 1% error. Let T(sub r) designate the average of the ion and neutral temperatures in K, and u = nu(sub d)/alpha, where nu(sub d) is the relative drift velocity in cm/s, and alpha = 4.56 x 10(exp 3) square root of T(sub r) cm/s is the thermal velocity of the O(+)-O system. Then, as u ranges from 0 to 2, the correction factor multiplying the collision frequency increases monotonically from 1 to about 1.5. An interesting result emerging from this calculation is that the correction factor for temperatures of aeronomical interest is to a good approximation independent of the temperature, depending only on the scaled velocity u.

Description: This paper describes an efficient Monte Carlo algorithm for choosing a new direction of a photon after a scattering interaction. The algorithm chooses a scattering angle by linear interpolation in a table of the inverse cumulative scattering probability. A Legendre expansion of the phase function makes it easy to apply Clenshaw's algorithm to build the interpolation table. The points in the table are close enough together that linear interpolation is accurate. With a table of 100,000 entries, we can keep the absolute and relative errors in matching the probability distribution below 10(exp -5).

Description: The relation between the solar wind input to the magetosphere, VB(sub South), and the auroral geomagnetic index AL is modeled with two linear moving-average filtering methods: linear prediction filters and a driven harmonic oscillator in the form of an electric circuit. Although the response of the three-parameter oscillator is simpler than the filter's, the methods yield similar linear timescales and values of the prediction-observation correlation and the prediction Chi(exp 2). Further the filter responses obtained by the two methods are similar in their long-term features. In these aspects the circuit model is equivalent to linear prediction filtering. This poses the question of uniqueness and proper interpretation of detailed features of the filters such as response peaks. Finally, the variation of timescales and filter responses with the AL activity level is discussed.

Description: Empirical data-based models of the magnetosphereic magnetic field have been widely used during recent years. However, the existing models (Tsyganenko, 1987, 1989a) have three serious deficiencies: (1) an unstable de facto magnetopause, (2) a crude parametrization by the K(sub p) index, and (3) inaccuracies in the equatorial magnetotail B(sub z) values. This paper describes a new approach to the problem; the essential new features are (1) a realistic shape and size of the magnetopause, based on fits to a large number of observed crossing (allowing a parametrization by the solar wind pressure), (2) fully controlled shielding of the magnetic field produced by all magnetospheric current systems, (3) new flexible representations for the tail and ring currents, and (4) a new directional criterion for fitting the model field to spacecraft data, providing improved accuracy for field line mapping. Results are presented from initial efforts to create models assembled from these modules and calibrated against spacecraft data sets.

Description: We have investigated the middle atmospheric response to the 27-day and 11-yr solar UV flux variations at low to middle latitudes using a two-dimensional photochemical model. The model reproduced most features of the observed 27-day sensitivity and phase lag of the profile ozone response in the upper stratosphere and lower mesosphere, with a maximum sensitivity of +0.51% per 1% change in 205 nm flux. The model also reproduced the observed transition to a negative phase lag above 2 mb, reflecting the increasing importance with height of the solar modulated HO(x) chemistry on the ozone response above 45 km. The model revealed the general anti-correlation of ozone and solar UV at 65-75 km, and simulated strong UV responses of water vapor and HO(x) species in the mesosphere. Consistent with previous 1D model studies, the observed upper mesospheric positive ozone response averaged over +/- 40 was simulated only when the model water vapor concentrations above 75 km were significantly reduced relative to current observations. In agreement with observations, the model computed a low to middle latitude total ozone phase lag of +3 days and a sensitivity of +0.077% per 1% change in 205 nm flux for the 27-day solar variation, and a total ozone sensitivity of +0.27% for the 11-yr solar cycle. This factor of 3 sensitivity difference is indicative of the photochemical time constant for ozone in the lower stratosphere which is comparable to the 27-day solar rotation period but is much shorter than the 11-yr solar cycle.

Description: The accuracy of Geosat satellite altimetry over the Greenland ice sheet is evaluated by comparing the measured heights to radar elevations from the airborne Greenland Aerogeophysics Project. At the center of the ice sheet where the ice surface is nearly level, surface comparisons show a fit at the 1 to 3 m level as expected, but even at moderately sloping ice regions (0.3 deg-0.6 deg), satellite altimetry mean errors in the range of 10 to 35 m are observed. These errors are found for slope-corrected and waveform-retracked data, so most previous accuracy estimates of current satellite altimetry ice sheet elevations in regions of slopping or undulating ice appear to be too optimistic.

Description: We report here on a number of examples of anomalous enhancements of eastward electric fields near sunrise in the equatorial ionospheric F-region. These examples were selected from the data base of the equatorial satellite, San Marco D (1988), which measured ionospheric electric fields during a period of solar minimum. The eastward electric fields reported correspond to vertical plasma drifts. The examples studied here are similar in signature and polarity to the pre-reversal electric field enhancements seen near sunset from ground-based radar systems. The morphology of these sunrise events, which are observed on about 14% of the morning-side satellite passes, are studied as a function of local zonal velocity, magnetic activity, geographic longitude and altitude. The nine events studied occur at locations where the zonal plasma flow is generally measured to be eastward, but reducing as a function of local time and at satellite longitudes where the magnetic declination has the opposite polarity as the declination of the sunrise terminator.

Description: A nonlinear filtering method is introduced for the study of the solar wind -- magnetosphere coupling and related to earlier linear techniques. The filters are derived from the magnetospheric state, a representation of the magnetospheric conditions in terms of a few global variables, here the auroral electrojet indices. The filters also couple to the input, a representation of the solar wind variables, here the rectified electric field. Filter-based iterative prediction of the indices has been obtained for up to 20 hours. The prediction is stable with respect to perturbations in the initial magnetospheric state; these decrease exponentially at the rate of 30/min. The performance of the method is examined for a wide range of parameters and is superior to that of other linear and nonlinear techniques. In the magnetospheric state representation the coupling is modeled as a small number of nonlinear equations under a time-dependent input.

Description: Magnetic reconnection between the interplanetary magnetic field (IMF) and the geomagnetic field is thought to play a major role in the transfer of solar wind momentum and energy to the magnetosphere. As the angle between the IMF and the geomagnetic field is changed at the bow of the magnetosphere, the topological record of the location of the reconnection region should be recorded in the magnetosheath and on the magnetopause along the flanks of the tail, because the super fast flow freezes strong magnetic gradients formed in the bow reconnection regions into the plasma downstream. In this report, we present results from a three-dimensional, magnetohydrodynamic (MHD), global numerical simulation code for the location of the separatrix between unconnected IMF magnetosheath field lines and reconnected field lines which penetrate the magnetopause and connect to the polar ionosphere. The angle between the IMF direction and the line where the separatrix crosses the magnetopause is shown to be a sensitive function of the IMF clock angle. We also explain how this behavior can be used to derive an approximate relation for the dependence of the cross-polar voltage on the IMF clock angle. We conclude with a note of caution concerning the importance of physical boundary conditions in magnetoplasma simulations.

Description: Many problems in geophysical and astrophysical convection systems are characterized by fast rotation and spherical shell geometry. The combined effects of Coriolis forces and spherical shell geometry produce a unique spatial symmetry for the convection pattern in a rapidly rotating spherical shell. In this paper, we first discuss the general spatial symmetries for rotating spherical shell convection. A special model, a spherical shell heated from below, is then used to illustrate how and when the spatial symmetries are broken. Symmetry breaking occurs via a sequence of spatial transitions from the primary conducting state to the complex multiple-layered columnar structure. It is argued that, because of the dominant effects of rotation, the sequence of spatial transitions identified from this particular model is likely to be generally valid. Applications of the spatial symmetry breaking to planetary convection problems are also discussed.

Description: Initial observations of a newly documented type of optical emission above thunderstorms are reported. 'Blue jets', or narrowly collimated beams of blue light that appear to propagate upwards from the tops of thunderstorms, were recorded on B/W and color video cameras for the first time during the Sprites94 aircraft campaign, June-July, 1994. The jets appear to propagate upward at speeds of about 100 km/s and reach terminal altitudes of 40-50 m. Fifty six examples were recorded during a 22 minute interval during a storm over Arkansas. We examine some possible mechanisms, but have no satisfactory theory of this phenomenon.

Description: It is shown from flux transfer event (FTE) occurrence statistics, observed as a function of MLT by the ISEE satellites, that recent 2-dimensional analytic theories of the effects of pulsed Petschek reconnection predict FTEs to contribute between 50 and 200 kV to the total reconnection voltage when the magnetosheath field points southward. The upper limit (200 kV) allows the possibility that FTEs provide all the antisunward transport of open field lines into the tail lobe. This range is compared with the voltages associated with series of FTEs signatures, as inferred from ground-based observations, which are in the range 10-60 kV. We conclude that the contribution could sometimes be made by a series of single, large events; however, the voltage is often likely to be contributed by several FTEs at different MLT.

Description: New investigations of the core ion motion within high-latitude topside ionosphere and near-Earth magnetosphere, using data from the Dynamics Explorer (DE) retarding ion mass spectrometer (RIMS), reveal the existence of significant regions of downward moving O(+). The occurences of downgoing versus upgoing O(+) are not clearly separable in terms of either polar zenith angle or Kp but are well distinquished by the direction of the z component of the interplanetary magnetic field (IMF). On the average, down flow dominates when IMFB2 less than O, while upflow dominates for IMFB2 greater than 0. Combining cross-field convection velocities derived from an empirical convection electric field model with the observed parallel velocities yields a two dimensional ion velocity field. This velocity field is consistent with a senario which has O(+) of cusp/cleft and auroral zone origin concvecting into the polar cap and, because of the dominance of gravitational energy over the upward kinetic energy, falling back into the inonsphere. This provides additional confirmation of the results of studies of the cleft ion fountain. Estimates of the flux of O(+) in the upflowing and downflowing regions for Lambda greater than or equal to 60 deg give a total upflow of approximately 6 X 10(exp 25) ions/sec for IMFB2 greater than 0 and total upflow and downflow of approximately 4 X 10 (exp 25) ions/sec and 1 x 10(exp 25) ions/sec, respectively, for IMFB less than 0. In all cases the magnitude of the dayside outflow is consistent with previous work on upwelling ions. While the magnitudes vary for high and low Kp the ratios of upward to downward flow are roughly the same at approximately 1.7. The downflowing O(+) shows a correlation with the magnitude of the outflow of light ions in the same region but the cause and effect of this relationship is not distinquishable. hable.

Description: We have developed the numerical algorithm for the computation of transient viscoelastic responses in the time domain for a radially stratified Earth model. Stratifications in both the elastic parameters and the viscosity profile have been considered. The particular viscosity profile employed has a viscosity maximum with a constrast of O(100) in the mid lower mantle. The distribution of relaxation times reveals the presence of a continuous spectrum situated between O(100) and O(exp 4) years. The principal mode is embedded within this continuous spectrum. From this initial-value approach we have found that for the low degree harmonics the non-modal contributions are comparable to the modal contributions. For this viscosity model the differences between the time-domain and normal-mode results are found to decrease strongly with increasing angular order. These calculations also show that a time-dependent effective relaxation time can be defined, which can be bounded by the relaxation times of the principal modes.

Description: The effects on double-probe electric field measurements induced by electron density and temperature gradients are investigated. We show that on some occasions such gradients may lead to marked spurious electric fields if the probes are assumed to lie at the same probe potential with repect to the plasma. The use of a proper bias current will decrease the magnitude of such an error. When the probes are near the plasma potential, the magnitude of these error signals, delta Epsilon, can vary as delta Epsilon approx. T(sub e)(Delta n(sub e)/n(sub e)) + 0.5 Delta T(sub e), where T(sub e) is the electron temperature, Delta n(sub e)/ n(sub e) the relative electron density variation between the two sensors, and Delta T(sub e) the electron temperature difference between the two sensors. This not only implies that the error signals will increase linearly with the density variations but also that such signatures grow with Delta T(sub e) i.e., such effects are 10 times larger in a 10-eV plasma than in a 1-eV plasma. This type of error is independent of the probe separation distance provided the gradient scale length is much larger than the distance. The largest errors occur when the probes are near to the plasma potential. During the crossing of a small structure (e.g, a double layer) the error signal appears as a bipolar signature. Our analysis shows that errors in double-probe measurements caused by plasma gradients are not significant at large scale (much greater than 1 km) plasma boundaries, and may only be important in cases where small-scale (less than 1 km), internal gradient structures exist. Bias currents tailored for each plasma parameter regime (i.e., variable bias current) would improve the double-probe response to gradient effects considerably.

Description: The evolution of the volcanic debris plume originating from the June 1991 eruption of Mt. Pinatubo has been monitored since its genesis using a ground-based backscatter lidar facility sited at the Jet Propulsion Laboratory (JPL). Both absolute and relative pre- and post-Pinatubo backscatter observations are in accord with Mie scattering projections based on measured aerosol particle size distributions reported in the literature. The post-Pinatubo column-integrated backscatter coefficient peaked approximately 400 days after the eruption, and the observed upper boundary of the aerosol column subsided at a rate of approximately 200 m/mon.

Description: Using Global Positioning System (GPS) receivers, we reoccupied several leveling benchmarks on the Kenai Peninsula of Alaska which had been surveyed by conventional leveling immediately following the March 27, 1964, Prince William Sound earthquake (M(sub w) = 9.3). By combining the two sets of measurements with a new, high-resolution model of the geoid in the region, we were able to determine the cumulative 1993-1964 postseismic vertical displacement. We find uplift at all of our benchmarks, relative to Seward, Alaska, a point that is stable according to tide gauge data. The maximum uplift of about 1 m occurs near the middle of the peninsula. The region of maximum uplift appears to be shifted northwest relative to the point of maximum coseismic subsidence. If we use tide gauge data at Nikishka and Seward to constrain the vertical motion, then the observed uplift has a trenchward tilt (down to the southeast) as well as an arching component. To explain the observations, we use creep-at-depth models. Most acceptable models require a fault slip of about 2.75 m, although this result is not unique. If the slip has been continuous since the 1964 earthquake, then the average slip rate is nearly 100 mm/yr, twice the plate convergence rate. Comparing the net uplift achieved in 29 years with that observed over 11 years in an adjacent region southeast of Anchorage, Alaska, we conclude that the rate of uplift is decreasing. A further decrease in the uplift rate is expected as the 29-year averaged displacement rate is about twice the plate convergence rate and therefore cannot be sustained over the entire earthquake cycle.

Description: Narrow enhancements of electron precipitation, with energy and flux well above typical values, have been observed with Dynamics Explorer 2 (DE 2) in the cusp/cleft region. The electron flux in the energy range 0.2-1 keV was 2 orders of magnitude higher in these structures than in the magnetosheath and were seen in approximaetly 80% of DE 2 cusp crossings at ionospheric altitudes. Typically, there was more than one electron structure in each cusp crossing. The position of these structures showed a systematic variation: for poleward ion dispersion (energy decreases with increasing latitude), electron structures were seen more often on the equatorial boundary of the cusp, while for equatorward ion dispersion (energy decreases with decreasing latitude), electron structures were more often seen on the poleward boundary. This suggests that the electron structures are associated with newly reconnected field lines. The electron spectra suggest that field-aligned acceleration processes could produce the electron structures, first near the boundary of the cusp/cleft during the reconnection of field lines and then in the cusp/cleft during the motion of reconnected flux tubes through the polar ionosphere.

Description: A permanent Global Positioning System (GPS) receiver at Casa Diablo Hot Springs, Long Valley Caldera, California was installed in January, 1993, and has operated almost continuously since then. The data have been transmitted daily to the Jet Propulsion Laboratory (JPL) for routine analysis with data from the Fiducial Laboratories for an International Natural sciences Network (FLINN) by the JPL FLINN analysis center. Results from these analyses have been used to interpret the on going deformation at Long Valley, with data excluded from periods when the antenna was covered under 2.5 meters of snow and from some periods when Anti Spoofing was enforced on the GPS signal. The remaining time series suggests that uplift of the resurgent dome of Long Valley Caldera during 1993 has been 2.5 +/- 1.1 cm/yr and horizontal motion has been 3.0 +/- 0.7 cm/yr at S53W in a no-net-rotation global reference frame, or 1.5 +/- 0.7 cm/yr at S14W relative to the Sierra Nevada block. These rates are consistent with uplift predicted from frequent horizontal strain measurements. Spectral analysis of the observations suggests that tidal forcing of the magma chamber is not a source of the variability in the 3 dimensional station location. These results suggest that remotely operated, continuously recording GPS receivers could prove to be a reliable tool for volcanic monitoring throughout the world.

Description: A comparison is made of the vertical distribution of high-level cloud tops derived from the Stratospheric Aerosol and Gas Experiment II (SAGE II) occultation measurements and from the International Satellite Cloud Climatology Project (ISCCP) for all Julys and Januarys in 1985 to 1990. The results suggest that ISCCP overestimates the pressure of high-level clouds by up to 50-150 mbar, particularly at low latitudes. This is caused by the frequent presence of clouds with diffuse tops (greater than 50% time when cloudy events are observed). The averaged vertical extent of the diffuse top is about 1.5 km. At midlatitudes where the SAGE II and ISCCP cloud top pressure agree best, clouds with distinct tops reach a maximum relative proportion of the total level cloud amount (about 30-40%), and diffuse-topped clouds are reduced to their minimum (30-40%). The ISCCP-defined cloud top pressure should be regarded not as the material physical height of the clouds but as the level which emits the same infrared radiance as observed. SAGE II and ISCCP cloud top pressures agree for clouds with distinct tops. There is also an indication that the cloud top pressures of optically thin clouds not overlying thicker clouds are poorly estimated by ISCCP at middle latitudes. The average vertical extent of these thin clouds is about 2.5 km.

Description: Global high-level clouds identified in Stratospheric Aerosol and Gas Experiment II (SAGE II) occultation measurements for January and July in the period 1985 to 1990 are compared with near-nadir-looking observations from the International Satellite Cloud Climatology Project (ISCCP). Global and zonal mean high-level cloud amounts from the two data sets agree very well, if clouds with layer extinction coefficients of less than 0.008/km at 1.02 micrometers wavelength are removed from the SAGE II results and all detected clouds are interpreted to have an average horizontal size of about 75 km along the 200 km transimission path length of the SAGE II observations. The SAGE II results are much more sensitive to variations of assumed cloud size than to variations of detection threshold. The geographical distribution of cloud fractions shows good agreement, but systematic regional differences also indicate that the average cloud size varies somewhat among different climate regimes. The more sensitive SAGE II results show that about one third of all high-level clouds are missed by ISCCP but that these clouds have very low optical thicknesses (less than 0.1 at 0.6 micrometers wavelength). SAGE II sampling error in monthly zonal cloud fraction is shown to produce no bias, to be less than the intraseasonal natural variability, but to be comparable with the natural variability at longer time scales.

Description: Enhancements in the fluxes of relativistic electrons trapped within the Earth's magnetosphere have been measured by the high-energy particle spectrometer, part of the particle environment monitor on the upper atmosphere research satellite (UARS). The largest increase in the electron fluxes with energies greater than 1 MeV observed on UARS from October 1991 through July 1994 was in early May 1992. The fluxes of trapped electrons in the drift loss cone and locally precipitating electrons showed differing buildup and decay rates as a function of invariant latitude. Increases of more than 2 orders of magnitude were observed in drift loss cone fluxes at magnetic latitudes of 40 deg-66 deg and in precipitating fluxes from 48 deg to 66 deg. The energy flux contained in the most intense local precipitation observed was approximately 0.1 erg/sq cm/s, entering the atmosphere and creating up to 1000 ion pairs/cu cm/s at 55-km altitude. The daily averaged energy flux from directly precipitating electrons with energies greater than 1 MeV deposited greater than 10(exp 20) erg/d worldwide into the atmosphere for the period May 12-21, 1992, producing greater than 10(exp 31) odd nitrogen molecules below 60-km altitude.

Description: Low temperature Zn diffusion in GaSb, where the minimum temperature was 450 C, was studied. The pseudo-closed box (PCB) method was used for Zn diffusion into GaAs, AlGaAs, InP, InGaAs and InGaAsP. The PCB method avoids the inconvenience of sealed ampoules and proved to be simple and reproducible. The special design of the boat for Zn diffusion ensured the uniformality of Zn vapor pressure across the wafer surface, and thus the uniformity of the p-GaSb layer depth. The p-GaSb layers were studied using Raman scattering spectroscopy and the x-ray rocking curve method. As for the postdiffusion processing, an anodic oxidation was used for a precise thinning of the diffused GaSb layers. The results show the applicability of the PCB method for the large-scale production of the GaSb structures for solar cells.

Description: Indium phosphide is considered to be a strong contender for many photovoltaic space applications because of its radiation resistance and its potential for high efficiency. An overview of recent progress is presented, and possible future research directions for indium phosphide space solar cells are discussed. The topics considered include radiation damage studies and space flight experiments.

Description: Understanding solar cell response to pulsed laser outputs is important for the evaluation of power beaming applications. The time response of high efficiency GaAs and silicon solar cells to a 25 nS monochromatic pulse input is described. The PC-1D computer code is used to analyze the cell current during and after the pulse for various conditions.

Description: Four-terminal mechanically stacked solar cells were developed for advanced space arrays with line-focus reflective concentrators. The top cells are based on AlGaAs/GaAs multilayer heterostructures prepared by low temperature liquid phase epitaxy. The bottom cells are based on heteroepitaxial InP/InGaAs liquid phase epitaxy or on homo-junction GaSb, Zn-diffused structures. The sum of the highest reached efficiencies of the top and bottom cells is 29.4 percent. The best four-terminal tandems have an efficiency of 27 to 28 percent. Solar cells were irradiated with 1 MeV electrons and their performances were determined as a function of fluence up to 10(exp 16) cm(exp-2). It was shown that the radiation resistance of developed tandem cells is similar to the most radiative stable AlGaAs/GaAs cells with a thin p-GaAs photoactive layer.

Description: An investigation of solar cells based on AlGaAs/GaAs heterostructures with an internal Bragg reflector as the back-surface reflector is presented. The Bragg reflector is grown by low pressure metalorganic chemical vapor deposition on n-GaAs substrates in a horizontal resistively heated reactor. The Bragg reflector with its maximum reflectance centered at a wavelength of 860 nm consists of 12 pairs of AlAs/GaAs layers. The resulting Bragg reflector has a thickness of 0.072 micrometers for AlAs and 0.059 micrometers for GaAs. The multi-layered quasi-dielectric stack selectively reflects weakly absorbed photons with energies near to the GaAs band gap for a second pass through the photoactive region, thus increasing the photocurrent. The use of the Bragg reflector allows the external quantum efficiency to be increased in the long wavelength of the spectrum. The use of the Bragg reflector and an antireflective coating and prismatic cover allowed an efficiency of 23.4 percent to be obtained.

Description: The high efficiency solar cells based on multilayer AlGaAs/GaAs heterostructures, prepared by low temperature liquid phase epitaxy (LPE), were developed and tested. An investigation of the low temperature LPE process for the crystallization of AlGaAs heterostructures of as high as 24.0 to 24.7 percent under AMO conditions at concentration ratios of 20 to 100x, were reached. Developed solar cells show substantial radiation resistance to the damage induced by 3.75 MeV electrons.

Description: A large set of bow shock crossings (i.e., 1392) observed by 17 spacecraft has been used to explore the three-dimensional shape and location of the Earth's bow shock and its dependence on solar wind and interplanetary magnetic field (IMF) conditions. This study investigates deviations from gas dynamic flow models associated with the magnetic terms in the magnetohydrodynamic (MHD) equations. Empirical models predicting the statistical position and shape of the bow shock for arbitrary values of the solar wind pressure, IMF, and Alfvenic Mach number (M(sub A)) have been derived. The resulting data set has been used to fit three-dimensional bow shock surfaces and to explore the variations in these surfaces with sonic (M(sub S)), Alfvenic (M(sub A)) and magnetosonic (M(sub MS)) Mach numbers. Analysis reveals that among the three Mach numbers, M(sub A) provides the best ordering of the least square bow shock curves. The subsolar shock is observed to move Earthward while the flanks flare outward in response to decreasing M(sub A); the net change represents a 6-10% effect. Variations due to changes in the IMF orientation were investigated by rotating the crossings into geocentric interplanetary medium coordinates. Past studies have suggested that the north-south extent of the bow shock surface exceeds the east-west dimension due to asymmetries in the fast mode Mach cone. This study confirms such a north-south versus east-west asymmetry and quantifies its variation with M(sub S), M(sub A), M(sub MS), and IMF orientation. A 2-7% effect is measured, with the asymmetry being more pronounced at low Mach numbers. Combining the bow shock models with the magnetopause model of Roelof and Sibeck (1993), variations in the magnetosheath thickness at different local times are explored. The ratio of the bow shock size to the magnetopause size at the subpolar point is found to be 1.46; at dawn and dusk, the ratios are found to be 1.89 and 1.93, respectively. The subsolar magnetosheath thickness is used to derive the polytropic index gamma according to the empirical relation of Spreiter et al. (1966). The resulting gamma = 2.3 suggests the empirical formula is inadequate to describe the MHD interaction between the solar wind and the magnetosphere.

Description: Radiation resulting from interaction between the effluent cloud of a space shuttle thruster and the ambient atmosphere was observed with a spectograph aboard the shutttle. The spectral measurements were made between 400 and 800 nm with a resolutoion of 3 nm. The primary emissions are identified as NO2, HNO, O(1)D, and O(1)S. These are the first observations od O(1)S emission in the shuttle plume. These data are compared with the previous measurements, and possible excitation mechanisms are discussed. The results are also compared with a Monte Carlo simulation of thruster plume-atmosphere interaction radiation.

Description: We report multi-instrument observations during an isolated substorm on 17 October 1989. The European Incoherent Scatter (EISCAT) radar operated in the SP-UK-POLI mode measuring ionospheric convection at latitudes 71 deg Lambda - 78 deg Lambda. Sub-Auroral Magnetometer Network (SAMNET) and the EISCAT Magnetometer Cross provide information on the timing of substorm expansion phase onset and subsequent intensifications, as well as the location of the field aligned and ionospheric currents associated with the substorm current wedge. Interplanetary Monitoring Platform-8 (IMP-8) magnetic field data are also included. Evidence of a substorm growth phase is provided by the equatorward motion of a flow reversal boundary across the EISCAT radar field of view at 2130 MLT, following a southward turning of the interplanetary magnetic field (IMF). We infer that the polar cap expanded as a result of the addition of open magnetic flux in the tail lobes during this interval. The flow reversal boundary, which is a lower limit to the polar cap boundary, reached an invariant latitude equatorward of 71 deg Lambda by the time of the expansion phase onset. We conclude that the substorm onset region in the ionosphere, defined by the westward electrojet, mapped to a part of the tail radially earthward of the boundary between open and closed magnetic flux, the distant neutral line. Thus the substorm was not initiated at the distant neutral line, although there is evidence that it remained active during the expansion phase.

Description: The impact of the extraterrestrial object that formed the Chicxulub crater in the northwestern Yucatan peninsula of Mexico is the leading suspect for the extinction of the dinosaurs. This article reports on a Planetary Society expedition to Albion Island in the Rio Hondo region of Belize to investigate evidence supporting the impact theory.

Description: Direct current (DC) electric field and ion density measurements near density depletion regions (that is, equatorial plasma bubbles) are used to estimate the vertical neutral wind speed. The measured zonal electric field in a series of density depletions crossed by the San Marco D satellite at 01.47-01.52 Universal Time (UT) on 25 October 1988, can be explained if a downward neutral wind of 15-30 m/s exists. Simultaneously, the F-region plasma was moving downward at a speed of 30-50 m/s. These events appear in the local time sector of 23.00-23.15 in which strong downward neutral winds may occur. Indeed, airglow measurements suggest that downward neutral velocities of 25-50 m/s are possible at time near midnight in the equatorial F-region.

Description: An overview of the observations of backstreaming electrons in the foreshock and the mechanisms that have been proposed to explain their properties will be presented. A primary characteristic of observed foreshock electrons is that their velocity distributions are spatially structured in a systematic way depending on distance from the magnetic field line which is tangent to the shock. There are two interrelated aspects to explaining the structure of velocity distributions in the foreshock, one involving the acceleration mechanism and the other, propagation from the source to the observing point. First, the source distribution of electrons energized by the shock must be determined along the shock surface. Proposed acceleration mechanisms include magnetic mirroring of incoming solar wind particles and mechanisms involving transmission of particles through the shock. Secondly, the kinematics of observable electrons streaming away from a curved shock with an initial parallel velocity and a downstream perpendicular velocity component due to the motional electric field must be determined. This is the context in which the observations and their explanations will be reviewed.

Description: Based on the conservation of chemical elements in chemical reactions, a rule is proved that the number of boundary conditions given by densities and/or nonzero velocities should not be less than the number of chemical elements in the system, and the boundary conditions for species given by densities and velocities should include all elements in the system. Applications of this rule to Mars are considered. It is shown that the problem of the CO2-H2O chemistry in the lower and middle atmosphere of Mars, say, in the range of 0-80 km does not have a unique solution, if only CO2 and H2O densities are given at the lower boundary, and the remaining boundary conditions are fluxes. Two examples of models of this type are discussed. Two models of the photochemistry of the Martian atmosphere, with and without nitrogen chemistry, are considered. The oxygen nonthermal escape ratio of 1.2 x 10(exp 8)/cu cm/s is given at 240 km and is balanced with the total hydrogen escape rate within an uncertainty of 1% for both models. Both models fit the measured O2 and CO mixing ratios, the O3 abundance, and the O2 1.27-micrometer dayglow almost within the uncertainties of the measured values, though the model without nitrogen chemistry fits better. The importance of nitrogen chemistry in the lower and middle atmosphere of Mars depends on a fine balance between production of NO and N in the upper atmosphere which is not known within the required accuracy.

Description: Goethite-bearing samples with values of Mn(s) (Mn/(Mn+Fe) mole fraction) up to 0.206 were synthesized by precipitation from alkaline solution. Samples with Mn(s) less than or equal 0.061 were single-phase Mn-goethites: samples with higher Mn(s) values contained another Mn-bearing phase (probably jacobsite). Mn-hematites were prepared by dehydroxylation of corresponding Mn-goethites at 500 C. Orthorhombic a and b unit cell dimensions of Mn-goethites changed in a linear manner with Mn(s), but not at rates predicted by the Vegrad law. Hexagonal unit cell dimensions of Mn-hematites did not vary with Mn(s). Moessbauer parameters isomer shift (IS), quadrupole splitting (QS), and hyperfine field (B(sub hf)) were measured at 293 and 15 K. For all single-phase Mn-goethites and Mn-hematites (Mn(s) less than or equal 0.061), magnetic splitting was observed at both temperatures. At 293 K, small but systematic decreases in B(sub hf) were observed with increasing Mn substitution; IS and QS were not dependent on Mn(s). Mn substitution strongly lowered the Morin transition temperature of hematite. At 15 K, the Morin transition was not present for Mn(s) greater than 0.020(4). The saturation magnetization of Mn-goethites and Mn-hematites (Mn(s) less than or equal 0.061) was the expected zero (within error) for antiferromagnetic goethite and for hematites obtained from dehydroxylation of goethites. Mn-geothites with Mn(s) greater than 0.061 were magnetic because of the presence of strongly magnetic jacobsite. For reflectivity spectra, bands resulting from MN(3+) were centered near 454 and 596 nm for Mn-goethites and near 545 and 700 nm for Mn-hematites. There is evidence for a approximately 700 nm band in spectral data for Martian bright regions, but association of it with Mn(3+) is not a unique interpretation. Comparison of laboratory and Martian spectral data implies that Mn(s) less than 0.032 for the Mn(3+) content of Martian hematites.

Description: Ferric-iron-bearing materials play an important role in the interpretation of visible to near-IR Mars spectra, and they may play a similarly important role in the analysis of new mid-IR spacecraft spectral observations to be obtained over the next decade. We review exisiting data on mid-IR transmission spectra of ferric oxides/oxyhydroxides and present new transmission spectra for ferric-bearing materials spanning a wide range of mineralogy and crystallinity. These materials include 11 samples of well-crystallized ferric oxides (hematite, maghemite, and magnetite) and ferric oxyhydroxides (goethite, lepidocrocite). We also report the first transmission spectra for purely nanophase ferric oxide samples that have been shown to exhibit spectral similarities to Mars in the visible to near-IR and we compare these data to previous and new transmission spectra of terrestial palagonites. Most of these samples show numerous, diagnostic absorption features in the mid-IR due to Fe(3+) - 0(2-) vibrational transitions, structural and/or bound OH, and/or silicates. These data indicate that high spatial resolution, moderate spectral resolution mid-IR ground-based and spacecraft observations of Mars may be able to detect and uniquely discriminate among different ferric-iron-bearing phases on the Martian surface or in the airborne dust.

Description: This paper focuses on the comparison of cloud amounts derived from an atmospheric general circulation model (AGCM), Satellite-observed clouds, and Ground-based cloud observations. Unlike Earth Radiation Budget Experiment (ERBE)-type comparisons it does not mix potential errors in the cloud amount with those in the radiation code embedded in the model. Long term cloud climatologies were used to compare global cloud amounts and regional seasonal cycles. The AGCM successfully reproduced the signatures of the warm pool and North Pacific seasonal cycle cloudiness but failed in the low stratus region off the coast of South America, a known problem for AGCMs. The data sets also reproduced the anomaly signature associated with El Nino in the warm pool region, but the model amounts were lower. Global results had a similar success rate, with the model generally producing lower total cloud compared to the satellite and in situ measurements. To compare cloud vertical distributions the cloud height may need to be validated using the corresponding radiation fields. Unfortunately there were also some large discrepancies between the two observed cloud data sets. While tremendously improved over the last decade the character of the observed cloud data sets, must be substantially enhanced before they will be useful in validating AGCMs by any but the crudest levels of comparison.

Description: The Photogrammetric Appendage Structural Dynamics Experiment (PASDE) is a Hitchhiker payload scheduled to fly as part of the International Space Station (ISS) Phase-1 flight program to the Russian Space Station Mir. The objective of the first flight of PASDE on STS-74 is to obtain video images of the Mir Kvant-2 solar array response to various structural dynamic excitation events. This experiment will demonstrate the use of photogrammetric techniques for on-orbit structural dynamics measurements. Photogrammetric measurements will provide a low cost alternative to appendage mounted accelerometers to the ISS program. The PASDE experiment hardware consists of three instruments each containing two video cameras, two video tape recorders, a modified video signal time inserter, and associated avionics boxes. The instruments were designed and built at the NASA Langley Research Center, and are integrated into standard Hitchhiker canisters at the NASA Goddard Space Flight Center. The Hitchhiker canisters are then installed into the Space Shuttle cargo bay in locations selected to achieve good video coverage and photogrammetric geometry. The measurement resolution of the instruments is expected to be on the order of 0.25 cm (0.1 in.).

Description: The Thermal Energy Storage-1 (TES-1) is a flight experiment that flew on the Space Shuttle Columbia (STS-62), in March 1994, as part of the OAST-2 mission. TES-1 is the first experiment in a four experiment suite designed to provide data for understanding the long duration microgravity behavior of thermal energy storage fluoride salts that undergo repeated melting and freezing. Such data have never been obtained before and have direct application for the development of space-based solar dynamic (SD) power systems. These power systems will store solar energy in a thermal energy salt such as lithium fluoride or calcium fluoride. The stored energy is extracted during the shade portion of the orbit. This enables the solar dynamic power system to provide constant electrical power over the entire orbit. Analytical computer codes have been developed for predicting performance of a spaced-based solar dynamic power system. Experimental verification of the analytical predictions is needed prior to using the analytical results for future space power design applications. The four TES flight experiments will be used to obtain the needed experimental data. This paper will focus on the flight results from the first experiment, TES-1, in comparison to the predicted results from the Thermal Energy Storage Simulation (TESSIM) analytical computer code. The TES-1 conceptual development, hardware design, final development, and system verification testing were accomplished at the NASA lewis Research Center (LeRC). TES-1 was developed under the In-Space Technology Experiment Program (IN-STEP), which sponsors NASA, industry, and university flight experiments designed to enable and enhance space flight technology. The IN-STEP Program is sponsored by the Office of Space Access and Technology (OSAT).

Description: Movement of water as vapor in the atmosphere is a fundamental process in the Earth's hydrological cycle. Investigations of spatial and time scales of water vapor transport in the atmosphere are important areas of research. Water vapor transmits energy as a function of its abundance across the spectrum. This is shown in the 400- to 2500-nm spectral region where the transmission of the terrestrial atmosphere has been modeled using the MODTRAN radiative transfer code for a range of water vapor abundances. Based on these model results, spectra measured by the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) have been used to investigate the movement of water vapor at 20-m spatial resolution over an 11-by-30-km area at approximately 15-minute time intervals (1.25 hours total). AVIRIS measures the upwelling spectral radiance from 400 to 2500 nm at 10-nm spectral intervals and collects images of 11-by-up-to-1000 km at 20-m spatial resolution. Data are collected at a rate of 1 km of flight line per 4.5 seconds. A set of five AVIRIS flight lines was acquired in rapid succession over Rogers Dry Lake, CA on May 18, 1993 at 18:59, 19:13, 19:29, 19:47, and 19:59 UTC. Rogers Dry lake is located 2 hours north of Los Angeles, CA at 34.84 degrees north latitude and 117.83 degrees west longitude in the Mojave Desert.

Description: The Voyager 1 spacecraft flew by Jupiter on March 5, 1979. Spacecraft navigation was performed with radio tracking data from NASA's Deep Space Network. In the years since then, there has been a great deal of progress in the definition of celestial reference frames and in determining the orbit and orientation of the Earth. Using these improvements, the radio metric range and Doppler data acquired from the Voyager 1 spacecraft near its encounter with Jupiter have been reanalyzed to determine the plane-of-sky position of Jupiter with much greater accuracy than was possible at the time of the encounter. The position of Jupiter at the time of encounter has been determined with an accuracy of 40 nrad in right ascension and 140 nrad in declination with respect to the celestial reference frame defined by the International Earth Rotation Service. This position estimate has been done to improve the ephemeris of Jupiter prior to the upcoming encounter of the Galileo spacecraft with Jupiter.

Description: The Boeing Company, under contract to the Electric Power Research Institute (EPRI), has completed a test program on the Mod-2 wind turbines at Goodnoe Hills, Washington. The objectives were to update fatigue load spectra, discern site and machine differences, measure vortex generator effects, and to evaluate rotational sampling techniques. This paper shows the test setup and loads instrumentation, loads data comparisons and test/analysis correlations. Test data are correlated with DYLOSAT predictions using both the NASA interim turbulence model and rotationally sampled winds as inputs. The latter is demonstrated to have the potential to improve the test/analysis correlations. The paper concludes with an assessment of the importance of vortex generators, site dependence, and machine differences on fatigue loads. The adequacy of prediction techniques used are evaluated and recommendations are made for improvements to the methodology.

Description: Analytical investigation of dynamic stall on HAWT (horizontal-axis wind turbines) rotor loads was conducted. Dynamic stall was modeled using the Gormont approach on the MOD-2 rotor, treating the blade as a rigid body teetering about a fixed axis. Blade flapwise bending moments at station 370 were determined with and without dynamic stall for spatial variations in local wind speed due to wind shear and yaw. The predicted mean flapwise bending moments were found to be in good agreement with test results. Results obtained with and without dynamic stall showed no significant difference for the mean flapwise bending moment. The cyclic bending moments calculated with and without dynamic stall effects were substantially the same. None of the calculated cyclic loads reached the level of the cyclic loads measured on the MOD-2 using the Boeing five-minute-average technique.

Description: Empirical equations are presented with which to model rotationally-sampled (R-S) turbulence for input to structural-dynamic computer codes and the calculation of wind turbine fatigue loads. These equations are derived from R-S turbulence data which were measured at the vertical-plane array in Clayton, New Mexico. For validation, the equations are applied to the calculation of cyclic flapwise blade loads for the NASA/DOE Mod-2 2.5-MW experimental HAWT's (horizontal-axis wind turbines), and the results compared to measured cyclic loads. Good correlation is achieved, indicating that the R-S turbulence model developed in this study contains the characteristics of the wind which produce many of the fatigue loads sustained by wind turbines. Empirical factors are included which permit the prediction of load levels at specified percentiles of occurrence, which is required for the generation of fatigue load spectra and the prediction of the fatigue lifetime of structures.

Description: This paper describes the engineering thought process behind the failure analysis, redesign, and rework of the flight hardware for the Brilliant Eyes Thermal Storage Unit (BETSU) experiment. This experiment was designed to study the zero-g performance of 2-methylpentane as a suitable phase change material. This hydrocarbon served as the cryogenic storage medium for the BETSU experiment which was flown 04 Mar 94 on board Shuttle STS-62. Ground testing had indicated satisfactory performance of the BETSU at the 120 Kelvin design temperature. However, questions remained as to the micro-gravity performance of this unit; potential deviations in ground (1 g) versus space flight (0 g) performance, and how the unit would operate in a realistic space environment undergoing cyclical operation. The preparations and rework performed on the BETSU unit, which failed initial flight qualification, give insight and lessons learned to successfully develop and qualify a space flight experiment.

Description: A computational procedure is presented for evaluating the sensitivity coefficients of the dynamic axisymmetric response of viscoplastic shells of revolution. The analytical formulation is based on Reissner's large deformation shell theory with the effects of transverse shear deformation, rotatory inertia and moments turning around the normal to the middle surface included. The material model is chosen to be isothermal viscoplasticity, and an associated flow rule is used with a von Mises effective stress. A mixed formulation is used with the fundamental unknowns consisting of six stress resultants, three generalized displacements and three velocity components. Spatial discretization is performed using finite elements, with discontinuous stress resultants across element interfaces. The temporal integration is performed by using an explicit central difference scheme (leap-frog method) with an implicit constitutive update. The sensitivity coefficients are evaluated using a direct differentiation approach. Numerical results are presented for a spherical cap subjected to step loading, and a circular plate subjected to impulsive loading. The sensitivity coefficients are generated by evaluating the derivatives of the response quantities with respect to thickness, mass density, Young's modulus, and two of the material parameters characterizing the viscoplastic response. Time histories of the response and sensitivity coefficients are presented, along with spatial distributions of these quantities at selected times.

Description: The Solar Array Module Plasma Interactions Experiment (SAMPIE) is a flight experiment that flew on the Space Shuttle Columbia (STS-62) in March 1994, as part of the OAST-2 mission. The overall objective of SAMPIE was to determine the adverse environmental interactions within the space plasma of low earth orbit (LEO) on modern solar cells and space power system materials which are artificially biased to high positive and negative direct current (DC) voltages. The two environmental interactions of interest included high voltage arcing from the samples to the space plasma and parasitic current losses. High voltage arcing can cause physical damage to power system materials and shorten expected hardware life. parasitic current losses can reduce power system efficiency because electric currents generated in a power system drain into the surrounding plasma via parasitic resistance. The flight electronics included two programmable high voltage DC power supplies to bias the experiment samples, instruments to measure the surrounding plasma environment in the STS cargo bay, and the on-board data acquisition system (DAS). The DAS provided in-flight experiment control, data storage, and communications through the Goddard Space Flight Center (GSFC) Hitchhiker flight avionics to the GSFC Payload Operations Control Center (POCC). The DAS and the SAMPIE POCC computer systems were designed for telescience operations; this paper will focus on the experiences of the SAMPIE team regarding telescience development and operations from the GSFC POCC during STS-62. The SAMPIE conceptual development, hardware design, and system verification testing were accomplished at the NASA Lewis Research Center (LeRC). SAMPIE was developed under the In-Space Technology Experiment Program (IN-STEP), which sponsors NASA, industry, and university flight experiments designed to enable and enhance space flight technology. The IN-STEP Program is sponsored by the Office of Space Access and Technology (OSAT).

Description: A coordinated effort has been underway over the past four years to elevate unstructured-grid methodology to a mature level. The goal of this endeavor is to provide a validated capability to non-expert users for performing rapid aerodynamic analysis and design of complex configurations. The Euler component of the system is well developed, and is impacting a broad spectrum of engineering needs with capabilities such as rapid grid generation and inviscid flow analysis, inverse design, interactive boundary layers, and propulsion effects. Progress is also being made in the more tenuous Navier-Stokes component of the system. A robust grid generator is under development for constructing quality thin-layer tetrahedral grids, along with a companion Navier-Stokes flow solver. This paper presents an overview of this effort, along with a perspective on the present and future status of the methodology.

Description: Hypervelocity impact experiments were performed to further test the survivability of carbonaceous impactors and to determine potential products that may have been synthesized during impact. Diamonds were launched by the Ames two-stage light gas gun into Al plate at velocities of 2.75 and 3.1 km sec(exp -1). FESEM imagery confirms that diamond fragments survived in both experiments. Earlier experiments found that diamonds were destroyed on impact above 4.3 km sec(exp -1). Thus, the upper stability limit for diamond on impact into Al, as determined from our experimental conditions, is between 3.1 and 4.3 km sec(exp -1). Particles of the carbonaceous chondrite Nogoya were also launched into Al at a velocity of 6.2 km sec (exp -1). Laser desorption (L (exp 2) MS) analyses of the impactor residues indicate that the lowest and highest mass polycyclic aromatic hydrocarbons (PAH's) were largely destroyed on impact; those of intermediate mass (202-220 amu) remained at the same level or increased in abundance. In addition, alkyl-substituted homologs of the most abundant pre-impacted PAH's were synthesized during impact. These results suggest that an unknown fraction of some organic compounds can survive low to moderate impact velocities and that synthesized products can be expected to form up to velocities of, at least, 6.5 km sec(exp -1). We also present examples of craters formed by a unique microparticle accelerator that could launch micron-sized particles of almost any coherent material at velocities up to approximately 15 km sec(exp -1). Many of the experiments have a direct bearing on the interpretation of LDEF craters.

Description: Interpretation of the wealth of impact data available from the Long Duration Exposure Facility, in terms of the absolute and relative populations of space debris and natural micrometeoroids, requires three dimensional models of the distribution of impact directions, velocities and masses of such particles, as well as understanding of the impact processes. Although the stabilized orbit of LDEF provides limited directional information, it is possible to determine more accurate impact directions from detailed crater morphology. The applicability of this technique has already been demonstrated but the relationship between crater shape and impactor direction and velocity has not been derived in detail. We present the results of impact experiments and simulations: (1) impacts at micron dimensions using the Unit's 2MV Van de Graaff accelerator; (2) impacts at mm dimensions using a Light Gas Gun; and (3) computer simulations using AUTODYN-3D from which an empirical relationship between crater shape and impactor velocity, direction and particle properties we aim to derive. Such a relationship can be applied to any surface exposed to space debris or micrometeoroid particles for which a detailed pointing history is available.

Description: A 5.2 mm crater in Al-metal represents the largest found on LDEF. We have examined this crater by field emission scanning electron microscopy (FESEM), energy dispersive spectroscopy (EDS) and time-of-flight/secondary ion mass spectroscopy (TOF-SIMS) in order to determine if there is any evidence of impactor residue. Droplet and dome-shaped columns, along with flow features, are evidence of melting. EDS from the crater cavity and rim show Mg, C, O and variable amounts of Si, in addition to Al. No evidence for a chondritic impactor was found, and it hypothesized that the crater may be the result of impact with space debris.

Description: On July 4,1997, the Mars Pathfinder spacecraft lands on Mars and starts conducting technological and scientific experiments. One experiment, the Alpha-Proton-X-ray Spectrometer, uses a sensor head placed against rocks and soil to determine their composition. To guarantee proper placement, a deployment mechanism mounted on the Mars Rover aligns the sensor head to within 20 deg of the rock and soil surfaces. In carrying out its task, the mechanism mimics the action of a human hand and arm. Consisting of a flexible wrist, a parallel link arm, a brush dc motor actuator, and a revolutionary non-pyrotechnic fail-safe release device, the mechanism correctly positions the sensor head on rocks as high as 0.29 m and on targets whose surfaces are tilted as much as 45 deg from the nominal orientation of the sensor head face. The mechanism weighs less than 0.5 kg, can withstand 100 g's, and requires less than 2.8 N x m of actuation torque. The fail-safe coupler utilizes Cerrobend, a metal alloy that melts at 60 C, to fuse the actuator and the rest of the mechanism together. A film heater wrapped around the coupler melts the metal, and Negator springs drive the mechanism into its stowed position. The fail-safe actuates using 6.75 Watts for 5 minutes in the event of an actuator failure.

Description: The finite deformation field of a plane strain Mode 1 crack in a hyperelastic and incompressible material was examined under the assumptions of small scale nonlinearity. Finite element analyses were performed for two different material laws, a Neo-Hookean material and a third order invariant of a Rivlin material. The numerical results for both materials were compared to the appropriate theoretical asymptotic solution. A local cavitation locus surrounding the crack tip was identified for the Neo-Hookean material. For the third order invariant Rivlin material, maximum values of the dominant stress component were found close to the surface of the crack, above and below the deformed crack tip.

Description: AstroPower is developing InGaAsSb thermophotovoltaic (TPV) devices. This photovoltaic cell is a two-layer epitaxial InGaAsSb structure formed by liquid-phase epitaxy on a GaSb substrate. The (direct) bandgap of the In(1 - x)Ga(x)As(1 -y)Sb(y) alloy is 0.50 to 0.55 eV, depending on its exact alloy composition (x, y); and is closely lattice-matched to the GaSb substrate The use of the quaternary alloy, as opposed to a ternary alloy - such as, for example, InGaAs/InP - permits low bandgap devices optimized for 1000 to 1500 C thermal sources with, with at the time, near-exact lattice matching to the GaSb substrate. Lattice-matching is important since even a small degree of lattice mismatch degrades device performance and reliability and increases processing complexity. For bandgaps of 0.52 eV,Fo internal quantum efficiencies as high as 95% have been measured at a wavelength of 2 microns. At 1 micron wavelengths, internal quantum efficiencies of 55% have been observed. The open-circuit voltage at currents of 0.3 A/sq cm is 0.220 volts and 0.260 V for current densities of 2 A/sq cm. Fill factors of 56% have also been measured. These preliminary results lead to the conclusion that the GaSb-based quaternary compounds provide a viable and high performance energy conversion solution for thermophotovoltaic systems operating with 1000 to 1500 C source temperatures.

Description: InP p(+)nn(+) MOCVD solar cells were irradiated with 0.2 MeV and 10 MeV protons to a fluence of 10(exp 13)/sq cm. The degradation of power output, IV behavior, carrier concentration and defect concentration were observed at intermediate points throughout the irradiations. The 0.2 MeV proton irradiated solar cells suffered much greater and more rapid degradation in power output than those irradiated with 10 meV protons. The efficiency losses were accompanied by larger increases in the recombination currents in the 0.2 MeV proton irradiated solar cells. The low energy proton irradiations also had a larger impact on the series resistance of the solar cells. Despite the radiation induced damage, the carrier concentration in the base of the solar cells showed no reduction after 10 MeV or 0.2 MeV proton irradiations and even increased during irradiation with 0.2 MeV protons. In a DLTS study of the irradiated samples, the minority carrier defects H4 and H5 at E(v) + 0.33 and E(v) + 0.52 eV and the majority carrier defects E7 and E10 at E(c)- 0.39 and E(c)-0.74 eV, were observed. The defect introduction rates for the 0.2 MeV proton irradiations were about 20 times higher than for the 10 MeV proton irradiations. The defect E10, observed here after irradiation, has been shown to act as a donor in irradiated n-type InP and may be responsible for obscuring carrier removal. The results of this study are consistent with the much greater damage produced by low energy protons whose limited range causes them to stop in the active region of the solar cell.

Description: We have carried out three (piggyback) radon-related projects aboard the KAO. The first, which was limited to upper tropospheric measurements while in level flight, revealed the systematic occurrence of unexpectedly high radon concentrations in this region of the atmosphere. The second project was an instrument development project, which led to the installation of an automatic radon measurement system aboard the NASA ER-2 High Altitude Research Aircraft. In the third, we installed a new system capable of collecting samples during the normal climb and descent of the KAO. The results obtained in these projects have resulted in significant contributions to our knowledge of atmospheric transport processes, and are currently playing a key role in the validation of global circulation and transport models.

Description: A data base comprised of all available ionosphere satellite ion composition measurements - the Goddard Comprehensive Ionosphere Data Base (GCID) - has been set upon optical disks for convenient merging and accessing of data from different satellites. This data has recentely been expanded to include all accessible satellite electron density and plasma temperature measurements. This paper demonstrates with a couple examples, the potential of GCID as an alternative to empirical models for undertanding ionsophere physics and chemistry. Through binning techniques GCID provdes a tool that complements and improves on the International Reference Ionosphere model (IRI) in delineating the topology of high latitude ion composition. GCID provides at a glance a measure of the local statistical variabilty of the ion compostions that is not in empirical models such as IRI. Furthermusing the data to determine the spatial and geophysical parameter range over which the minor ion species are approximately in chemical equilibrium, the number of data points available for empirical ion composition models can be increased by using the statiscally more reliable neutral and electron empirical models to derive ions composition. Currently available empirical models are not capable of adequately defining the comples high latitude distribution - a measure of the variablity is needed.

Description: High voltage power systems in space will interact with the space plasma in a variety of ways. One of these, snapover, is characterized by sudden enlargement of the current collection area across normally insulating surfaces generating enhanced electron current collection. Power drain on solar array power systems results from this enhanced current collection. Optical observations of the snapover phenomena in the laboratory indicates a functional relation between glow area and bia potential as a consequence of the fold/cusp bifurcation in chaos theory. Successful characterizations of snapover as a chaotic phenomena may provide a means of snapover prevention and control through chaotic synchronization.

Description: The PASP-Plus (Photovoltaic Array Space Power Plus Diagnostics) program is a photovoltaic experiment which is flying on the Air Force satellite APEX (Advanced Photovoltaic and Electronics Experiments). The satellite was launched on August 3, 1994 with a Pegasus low-cost launch vehicle. There are two other small experiments on APEX, however PASP-Plus is the largest, uses the most power, and accounts for the largest portion of the data requirements. The satellite is in an elliptical orbit with an apogee of 2552 km and a perigee of 363 km. The inclination is 70 degrees. The PASP-Plus experiment consists of twelve photovoltaic panels containing a total of sixteen separate cell modules. Two of the modules are concentrator modules, while the rest are planar. There are several different solar cell types flying on PASP-Plus including silicon, GaAs on germanium substrates, InP, amorphous silicon, and three multi-bandgap cells. The purpose of this paper is to present some of the data from the first year of the PASP-Plus flight. Cell performance and module thermal performance will be discussed as well as other relevant data.

Description: The Astro Edge solar array is a new and innovative reflective low concentrator power generating system which has been selected for the CTA Incorporate/Lockheed Martin Clark spacecraft under the NASA Small Spacecraft Technology Initiative (SSTI) program. In support of this program, Astro Aerospace Corporation has produced one qualification and two flight solar array wings to support a July 1996 launch. The Astro Edge solar array was selected as a new technology to benefit future NASA, military and commercial missions by providing high specific power, high deployed stiffness, low stowed volume, low risk, and cost reduction features which meet the agency's 'better, faster, cheaper' goals. This novel array accounts for five of the thirty-six advanced technologies which the Clark spacecraft will demonstrate. A brief SSTI Astro Edge solar array program overview is presented. Completed qualification and acceptance testing is discussed. Finally, the major discriminators which make the Astro Edge solar array 'better, faster, cheaper' technology are provided.

Description: Gallium phosphide (GaP) energy converters may be successfully deployed to provide new mission capabilities for spacecraft. Betavoltaic power supplies based on the conversion of tritium beta decay to electricity using GaP energy converters can supply long term low-level power with high reliability. High temperature solar cells, also based on GaP, can be used in inward-bound missions greatly reducing the need for thermal dissipation. Results are presented for GaP direct conversion devices powered by Ni-63 and compared to the conversion of light emitted by tritiarated phosphors. Leakage currents as low as 1.2 x 10(exp -17) A/sq cm have been measured and the temperature dependence of the reverse saturation current is found to have ideal behavior. Temperature dependent IV, QE, R(sub sh), and V(sub oc) results are also presented. These data are used to predict the high-temperature solar cell and betacell performance of GaP devices and suggest appropriate applications for the deployment of this technology.

Description: Literature on solar array angle of incidence corrections was found to be sparse and contained no tabular data for support. This lack along with recent data on 27 GaAs/Ge 4 cm by 4 cm cells initiated the analysis presented in this paper. The literature cites seven possible contributors to angle of incidence effects: cosine, optical front surface, edge, shadowing, UV degradation, particulate soiling, and background color. Only the first three are covered in this paper due to lack of sufficient data. The cosine correction is commonly used but is not sufficient when the incident angle is large. Fresnel reflection calculations require knowledge of the index of refraction of the coverglass front surface. The absolute index of refraction for the coverglass front surface was not known nor was it measured due to lack of funds. However, a value for the index of refraction was obtained by examining how the prediction errors varied with different assumed indices and selecting the best fit to the set of measured values. Corrections using front surface Fresnel reflection along with the cosine correction give very good predictive results when compared to measured data, except there is a definite trend away from predicted values at the larger incident angles. This trend could be related to edge effects and is illustrated by a use of a box plot of the errors and by plotting the deviation of the mean against incidence angle. The trend is for larger deviations at larger incidence angles and there may be a fourth order effect involved in the trend. A chi-squared test was used to determine if the measurement errors were normally distributed. At 10 degrees the chi-squared test failed, probably due to the very small numbers involved or a bias from the measurement procedure. All other angles showed a good fit to the normal distribution with increasing goodness-of-fit as the angles increased which reinforces the very small numbers hypothesis. The contributed data only went to 65 degrees from normal which prevented any firm conclusions about extreme angle effects although a trend in the right direction was seen. Measurement errors were estimated and found to be consistent with the conclusions that were drawn. A controlled experiment using coverglasses and cells from the same lots and extending to larger incidence angles would probably lead to further insight into the subject area.

Description: Past NASA missions to Mars, Jupiter, and the outer planets were powered by radioisotope thermal generators (RTG's). Although these devices proved to be reliable, their high cost and highly toxic radioactive heat source has made them far less desirable for future planetary missions. This has resulted in a renewed search for alternate energy sources, some of them being photovoltaic (PV) and thermophotovoltaic (TPV). Both of these alternate energy sources convert light/thermal energy directly into electricity. In order to create a viable PV and TPV data base for planetary mission planners and cell designers, we have compiled low temperature low intensity (LILT) I-V data on single junction and multi-junction high efficiency solar cells. The cells tested here represent the latest photovoltaic technology. Using this LILT data to calculate dI(sub SC)/dT, dV(sub OC)/dT, dFF/dT, and also as a function of intensity, an accurate prediction of cell performance under the AMO spectrum can be determined. When combined with QUantum efficiency at Low Temperature (QULT) data, one can further enhance the data by adding spectral variations to the measurements. This paper presents an overview of LILT measurements and is only intended to be used as a guideline for material selection and performance predictions. As single junction and multi-junction cell technologies emerge, new test data must be collected. Cell materials included are Si, GaAs/Ge, GainP/GaAs/Ge, InP, InGaAs/InP, InP/InGaAs/InP, and GainP. Temperatures range as low as -175 C and intensities range from 1 sun to .02 suns.

Description: When determining the best solar cell technology for a particular space flight mission, accurate prediction of solar cell performance in a space radiation environment is essential. The current methodology used to make such predictions requires extensive experimental data measured under both electron and proton irradiation. Due to the rising cost of accelerators and irradiation facilities, such extensive data sets are expensive to obtain. Moreover, with the rapid development of novel cell designs, the necessary data are often not available. Therefore, a method for predicting cell degradation based on limited data is needed. Such a method has been developed at the Naval Research Laboratory based on damage correlation using 'displacement damage dose' which is the product of the non-ionizing energy loss (NIEL) and the particle fluence. Displacement damage dose is a direct analog of the ionization dose used to correlate the effects of ionizing radiations. In this method, the performance of a solar cell in a complex radiation environment can be predicted from data on a single proton energy and two electron energies, or one proton energy, one electron energy, and Co(exp 60) gammas. This method has been used to accurately predict the extensive data set measured by Anspaugh on GaAs/Ge solar cells under a wide range of electron and proton energies. In this paper, the method is applied to InP solar cells using data measured under 1 MeV electron and 3 MeV proton irradiations, and the calculations are shown to agree well with the measured data. In addition to providing accurate damage predictions, this method also provides a basis for quantitative comparisons of the performance of different cell technologies. The performance of the present InP cells is compared to that published for GaAs/Ge cells. The results show InP to be inherently more resistant to displacement energy deposition than GaAs/Ge.

Description: Advanced materials play an important role in electrochemical energy devices such as batteries, fuel cells, and electrochemical capacitors. They are being used as both electrodes and electrolytes. Sol-gel processing is a versatile solution technique used in fabrication of ceramic materials with tailored stoichiometry, microstructure, and properties. The application of sol-gel processing in the fabrication of advanced electrochemical energy materials will be presented. The potentials of sol-gel derived materials for electrochemical energy applications will be discussed along with some examples of successful applications. Sol-gel derived metal oxide electrode materials such as V2O5 cathodes have been demonstrated in solid-slate thin film batteries; solid electrolytes materials such as beta-alumina for advanced secondary batteries had been prepared by the sol-gel technique long time ago; and high surface area transition metal compounds for capacitive energy storage applications can also be synthesized with this method.

Description: Grid generation plays an integral part in the solution of computational fluid dynamics problems for aerodynamics applications. A major difficulty with standard structured grid generation, which produces quadrilateral (or hexahedral) elements with implicit connectivity, has been the requirement for a great deal of human intervention in developing grids around complex configurations. This has led to investigations into unstructured grids with explicit connectivities, which are primarily composed of triangular (or tetrahedral) elements, although other subdivisions of convex cells may be used. The existence of large gradients in the solution of aerodynamic problems may be exploited to reduce the computational effort by using high aspect ratio elements in high gradient regions. However, the heuristic approaches currently in use do not adequately address this need for high aspect ratio unstructured grids. High aspect ratio triangulations very often produce the large angles that are to be avoided. Point generation techniques based on contour or front generation are judged to be the most promising in terms of being able to handle complicated multiple body objects, with this technique lending itself well to adaptivity. The eventual goal encompasses several phases: first, a partitioning phase, in which the Voronoi diagram of a set of points and line segments (the input set) will be generated to partition the input domain; second, a contour generation phase in which body-conforming contours are used to subdivide the partition further as well as introduce the foundation for aspect ratio control, and; third, a Steiner triangulation phase in which points are added to the partition to enable triangulation while controlling angle bounds and aspect ratio. This provides a combination of the advancing front/contour techniques and refinement. By using a front, aspect ratio can be better controlled. By using refinement, bounds on angles can be maintained, while attempting to minimize the number of Steiner points.

Description: Approximately 20 sq m of protective thermal blankets, largely composed of Teflon, were retrieved from the Long Duration Exposure Facility (LDEF) after the spacecraft had spent approximately 5.7 years in space. Examination of these blankets revealed that they contained thousands of hypervelocity impact features ranging from micron-sized craters to penetration holes several millimeters in diameter. We conducted impact experiments in an effort to reproduce such features and to -- hopefully -- understand the relationships between projectile size and the resulting crater or penetration-hole diameter over a wide range of impact velocity. Such relationships are needed to derive the size- and mass-frequency distribution and flux of natural and man-made particles in low-Earth orbit. Powder propellant and light-gas guns were used to launch soda-lime glass spheres of 3.175 mm (1/8 inch) nominal diameter (Dp) into pure Teflon FEP targets at velocities ranging from 1 to 7 km/s. Target thickness (T) was varied over more than three orders of magnitude from infinite halfspace targets (Dp/T less than 0.1) to very thin films (Dp/T greater than 100). Cratering and penetration of massive Teflon targets is dominated by brittle failure and the development of extensive spall zones at the target's front and, if penetrated, the target's rear side. Mass removal by spallation at the back side of Teflon targets may be so severe that the absolute penetration-hole diameter (Dh) can become larger than that of a standard crater (Dc) at relative target thicknesses of Dp/T = 0.6-0.9. The crater diameter is infinite halfspace Teflon targets increases -- at otherwise constant impact conditions -- with encounter velocity by a factor of V0.44. In contrast, the penetration-hole size is very thin foils (Dp/T greater than 50) is essentially unaffected by impact velocity. Penetrations at target thicknesses intermediate to these extremes will scale with variable exponents of V. Our experimental matrix is sufficiently systematic and complete, up to 7 km/s, to make reasonable recommendations for the velocity-scaling of Teflon craters and penetrations. We specifically suggest that cratering behavior and associated equations dominate all impacts in which the shock-pulse duration of the projectile (tp) is shorter than that of the target (tt). We also demonstrate that each penetration hole from space-retrieved surfaces may be assigned a unique projectile size, provided an impact velocity is known or assumed. This calibration seems superior to the traditional ballistic-limit approach.

Description: A computational procedure is presented for evaluating the sensitivity coefficients of the dynamic frictional contact/impact response of axisymmetric composite structures. The structures are assumed to consist of an arbitrary number of perfectly bonded homogeneous anisotropic layers. The material of each layer is assumed to be hyperelastic, and the effect of geometric nonlinearity is included. The sensitivity coefficients measure the sensitivity of the response to variations in different material, lamination and geometric parameters of the structure. A displacement finite element model is used for the discretization. The normal contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of the nodal displacements, and the Lagrange multipliers associated with the contact conditions. The Lagrange multipliers are allowed to be discontinuous at interelement boundaries. Tangential contact conditions are incorporated by using a penalty method in conjunction with the classical Coulomb's friction model. Temporal integration is performed by using Newmark method. The Newton-Raphson iterative scheme is used for the solution of the resulting nonlinear algebraic equations, and for the determination of the contact region, contact conditions (sliding or sticking), and the contact pressures. The sensitivity coefficients are evaluated by using a direct differentiation approach. Numerical results are presented from the frictional contact/impact response of a composite spherical cap impacting on a rigid plate.

Description: A computational strategy is presented for postbuckling and nonlinear static analyses of large complex structures on distributed-memory parallel computers. The strategy is designed for message-passing parallel computer systems. The key elements of the proposed strategy are: (1) a multiple-parameter reduced basis technique; (2) a nested dissection (or multilevel substructuring) ordering scheme; (3) parallel assembly of global matrices; and (4) a parallel sparse equation solver. The effectiveness of the strategy is assessed by performing thermomechanical postbuckling analyses of stiffened composite panels with cutouts, and nonlinear large-deflection analyses of High Speed Civil Transport models on three distributed-memory computers. The numerical studies presented demonstrate the advantages of nested dissection-based solvers over traditional skyline-based solvers on distributed-memory machines.

Description: A computational strategy is presented for the nonlinear static and postbuckling analyses of large complex structures on massively parallel computers. The strategy is designed for distributed-memory, message-passing parallel computer systems. The key elements of the proposed strategy are: (1) a multiple-parameter reduced basis technique; (2) a nested dissection (or multilevel substructuring) ordering scheme; (3) parallel assembly of global matrices; and (4) a parallel sparse equation solver. The effectiveness of the strategy is assessed by applying it to thermo-mechanical postbuckling analyses of stiffened composite panels with cutouts, and nonlinear large-deflection analyses of HSCT models on Intel Paragon XP/S computers. The numerical studies presented demonstrate the advantages of nested dissection-based solvers over traditional skyline-based solvers on distributed memory machines.

Description: A reduced basis technique and a computational procedure are presented for generating the nonlinear vibrational response, and evaluating the first-order sensitivity coefficients of thin-walled composite frames. The sensitivity coefficients are the derivatives of the nonlinear frequency with respect to the material and lamination parameters of the frame. A mixed formulation is used with the fundamental unknowns consisting of both the generalized displacements and stress resultants in the frame. The flanges and webs of the frames are modeled by using geometrically nonlinear two-dimensional shell and plate finite elements. The computational procedure can be conveniently divided into three distinct steps. The first step involves the generation of various-order perturbation vectors, and their derivatives with respect to the material and lamination parameters of the frame, using the Linstedt-Poincare perturbation technique. The second step consists of using the perturbation vectors as basis vectors, computing the amplitudes of these vectors and the nonlinear frequency of vibration, via a direct variational procedure. The third step consists of using the perturbation vectors, and their derivatives, as basis vectors and computing the sensitivity coefficients of the nonlinear frequency via a second application of the direct variational procedure. Numerical results are presented for semicircular thin-walled frames with I and J section, showing the convergence of the nonlinear frequency and the sensitivity coefficients obtained by both the reduced-basis and perturbation techniques.

Description: An international conference on high-latitude ionospheric modeling produced 27 papers in the areas of ionospheric mapping, electron density and distribution, ion density and distribution, ionospheric storems, ionospheric composition, and ionospheric sounding techniques. Upgrades to the International Reference Ionosphere (IRI) model were proposed in several papers.

Description: Two symposia related to the exploration of the Venus and Mars atmospheres produced papers which described the results of spacecraft-based and ground-based observations of these two planets. The topics addressed included solar wind interactions, molecular ions, planetary magnetospheres, planetary ionospheres, planetary thermospheres, gravity waves, the homopause, radiative transfer, radiative balance, atmospheric water, planetary lightning, very low frequency bursts, cosmic-ray ionization, planetary loss, and future missions to study Mars aeronomy, energetic particles near Mars, radar exploration of Mars, and the Mars Environmental Survey (MESUR) mission.

Description: The experimental performance of a multijunction monolithic In(0.53)Ga(0. 47)As power converter under blackbody irradiation is reported. Eight InGaAs PN junctions grown epitaxially on a semi-insulating wafer were monolithically integrated in series to boost the approximately 0.4 V photovoltage per typical InGaAs junction to over 3 volts for the 1 sq cm chip. This chip was originally designed and characterized for free-space 1.3 micron laser power beaming. This is the first report of such a multijunction TPV. This is not a traditional tandem cell in which the junctions are stacked vertically. The junctions are each about 1mm long by 1 cm wide and are laterally connected across the 1 sq cm device area. This multijunction design has the potential for lower I(sup 2)R power loss since the smaller PN junction area limits the current to one-eighth that of the equivalent surface area. In essence, the current is traded for voltage to avoid the I(sup 2)R loss, analogous to the way power utilities avoid I(sup 2)P loss in high-tension power lines, by transforming the high current, low voltage generated at a power plant into a high voltage at a low current before transmitting the power over great distances.

Description: Progress is reported with respect to the development of ultra-lightweight, high performance, thin, light trapped GaAs solar cells for advanced space power systems. Conversion efficiencies of over 17.7% have been demonstrated for a 3 micron thick, 1 sq cm silicone bonded solar cell. This results in a specific power of over 1020 W/kg. Device parameters were 1.011 V open circuit voltage, 80% fill factor, and a short-circuit current density of 29.5 mA/sq cm . In addition to silicone bonding, the use of electrostatic bonding to attach the coverglass support to the front surface enables an ultra-thin, all back contact design that survives processing temperatures greater than 750 C. This also results in a 10% reduction of the cell weight for a potential specific power of 1270 W/kg. All back contact, ultra-thin, electrostatically bonded GaAs solar cell prototypes have been completed demonstrating an open circuit voltage of 1 volt for a cell base thickness of 1 micron with a 0.5 micron emitter. This technology will result in a revolutionary improvement in survivability, performance, and manufacturability of lightweight GaAs solar cell products for future Earth-orbiting science and space exploration missions. The thin, electrostatically bonded, all back contact GaAs device technology has multiple uses for specialty high performance solar cells and other optoelectronic devices.

Description: Multi-junction solar cells are attractive for space applications because they can be designed to convert a larger fraction of AMO into electrical power at a lower cost than single-junction cells. The performance of multi-junction cells is much more sensitive to the spectral irradiance of the illuminating source than single-junction cells. The design of high efficiency multi-junction cells for space applications requires matching the optoelectronic properties of the junctions to AMO spectral irradiance. Unlike single-junction cells, it is not possible to carry out quantum efficiency measurements using only a monochromatic probe beam and determining the cell short-circuit current assuming linearity of the quantum efficiency. Additionally, current-voltage characteristics can not be calculated from measurements under non-AMO light sources using spectral-correction methods. There are reports in the literature on characterizing the performance of multi junction cells by measuring and convoluting the quantum efficiency of each junction with the spectral irradiance; the technique is of limited value for the characterization of cell performance under AMO power-generating conditions. We report the results of research to develop instrumentation and techniques for characterizing multi junction solar cells for space . An integrated system is described which consists of a standard lamp, spectral radiometer, dual-source solar simulator, and personal computer based current-voltage and quantum efficiency equipment. The spectral radiometer is calibrated regularly using the tungsten-halogen standard lamp which has a calibration based on NIST scales. The solar simulator produces the light bias beam for current-voltage and cell quantum efficiency measurements. The calibrated spectral radiometer is used to 'fit' the spectral irradiance of the dual-source solar simulator to WRL AMO data. The quantum efficiency apparatus includes a monochromatic probe beam for measuring the absolute cell quantum efficiency at various voltage biases, including the voltage bias corresponding to the maximum-power point under AMO light bias. The details of the procedures to 'fit' the spectral irradiance to AMO will be discussed. An assessment of the role of the accuracy of the 'fit' of the spectral irradiance and probe beam intensity on measured cell characteristics will be presented. quantum efficiencies were measured with both spectral light bias and AMO light bias; the measurements show striking differences. Spectral irradiances were convoluted with cell quantum efficiencies to calculate cell currents as function of voltage. The calculated currents compare with measured currents at the 1% level. Measurements on a variety of multi-junction cells will be presented. The dependence of defects in junctions on cell quantum efficiencies measured under light and voltage bias conditions will be presented. Comments will be made on issues related to standards for calibration, and limitations of the instrumentation and techniques. Expeditious development of multi-junction solar cell technology for space presents challenges for cell characterization in the laboratory.

Description: It is well known that the behavior of III-V compound based solar cells is largely controlled by their surface, since the majority of light generated carriers (63% for GaAs and 79% for InP) are created within 0.2 mu m of the surface of the illuminated cell. Consequently, the always observed high surface recombination velocity (SRV) on these cells is a serious limiting factor for their high efficiency performance, especially for those with p-n junction made by either thermal diffusion or ion implantation. A good surface passivation layer, ideally a grown oxide as opposed to a deposited one, will cause a significant reduction in the SRV without adding interface problems, thus improving the performance of III-V compound based solar cells. Another significant benefit to the overall performance of the solar cells can be achieved by a substantial reduction of their large surface optical reflection by the use of a well designed antireflection (AR) coating. In this paper, we demonstrate the effectiveness of using a chemically grown thermally and chemically stable oxide, not only for surface passivation but also as an integral part of a 3-layer AR coating for thermally diffused p+n InP solar cells. A phosphorus-rich interfacial oxide, In(PO3)3, is grown at the surface of the p+ emitter using an etchant based on HNO3, o-H3PO4 and H2O2. This oxide has the unique properties of passivating the surface as well as serving as an efficient antireflective layer yielding a measured record high AMO open-circuit voltage of 890.3 mV on a thermally diffused InP(Cd,S) solar cell. Unlike conventional single layer AR coatings such as ZnS, Sb2O3, SiO or double layer AR coatings such as ZnS/MgF2 deposited by e-beam or resistive evaporation, this oxide preserves the stoichiometry of the InP surface. We show that it is possible to design a three-layer AR coating for a thermally diffused InP solar cell using the In(PO3)3 grown oxide as the first layer and Al2O3 and MgF2 as the second and third layers respectively, so as to yield an overall theoretical reflectance of less than 2%. Since chemical oxides are readily grown on III-V semiconductors materials, the technique of using the grown oxide layer to both passivate the surface as well as serve as the first of a multilayer AR coating should work well for all III-V compound-based solar cells.

Description: High-efficiency, heteroepitaxial (HE) InP solar cells, grown on GaAs, Si or Ge substrates, are desirable for their mechanically strong, light-weight and radiation-hard properties. However, dislocations, caused by lattice mismatch, currently limit the performance of the HE cells. This occurs through shunting paths across the active photovoltaic junction and by the formation of deep levels. In previous work we have demonstrated that plasma hydrogenation is an effective and stable means to passivate the electrical activity of dislocations in specially designed HE InP test structures. In this work, we present the first report of successful hydrogen passivation in actual InP cell structures grown on GaAs substrates by metalorganic chemical vapor deposition (MOCVD). We have found that a 2 hour exposure to a 13.56 MHz hydrogen plasma at 275 C reduces the deep level concentration in HE n+n InP cell structures from as-grown values of approximately 10(exp 15)/cm(exp -3), down to 1-2 x 10(exp 13)/cm(exp -3). The deep levels in the p-type base region of the cell structure match those of our earlier p-type test structures, which were attributed to dislocations or related point defect complexes. All dopants were successfully reactivated by a 400 C, 5 minute anneal with no detectable activation of deep levels. I-V analysis indicated a subsequent approximately 10 fold decrease in reverse leakage current at -1 volt reverse bias, and no change in the forward biased series resistance of the cell structure which indicates complete reactivation of the n+ emitter. Furthermore, electrochemical C-V profiling indicates greatly enhanced passivation depth, and hence hydrogen diffusion, for heteroepitaxial structures when compared with identically processed homoepitaxial n+p InP structures. An analysis of hydrogen diffusion in dislocated InP will be discussed, along with comparisons of passivation effectiveness for n+p versus p+n heteroepitaxial cell configurations. Preliminary hydrogen-passivated HE InP cell results will also be presented.