ALBERT

Description: During the period December 23,1997 and December August 31,2004, we accomplished the development of 2 CFD codes for DNS/LES/RANS simulation of turbine cascade flows, namely LESTool and UNCLE. LESTool is a structured code making use of 5th order upwind differencing scheme and UNCLE is a second-order-accuracy unstructured code. LESTool has both Dynamic SGS and Spalart's DES models and UNCLE makes use of URANS and DES models. The current report provides a description of methodologies used in the codes.

Description: Bangui anomaly is the name given to one of the Earth s largest crustal magnetic anomalies and the largest over the African continent. It covers two-thirds of the Central African Republic and therefore the name derives from the capitol city-Bangui that is also near the center of this feature. From surface magnetic survey data Godivier and Le Donche (1962) were the first to describe this anomaly. Subsequently high-altitude world magnetic surveying by the U.S. Naval Oceanographic Office (Project Magnet) recorded a greater than 1000 nT dipolar, peak-to-trough anomaly with the major portion being negative (figure 1). Satellite observations (Cosmos 49) were first reported in 1964, these revealed a 40nT anomaly at 350 km altitude. Subsequently the higher altitude (417-499km) POGO (Polar Orbiting Geomagnetic Observatory) satellite data recorded peak-to-trough anomalies of 20 nT these data were added to Cosmos 49 measurements by Regan et al. (1975) for a regional satellite altitude map. In October 1979, with the launch of Magsat, a satellite designed to measure crustal magnetic anomalies, a more uniform satellite altitude magnetic map was obtained. These data, computed at 375 km altitude recorded a -22 nT anomaly (figure 2). This elliptically shaped anomaly is approximately 760 by 1000 km and is centered at 6%, 18%. The Bangui anomaly is composed of three segments; there are two positive anomalies lobes north and south of a large central negative field. This displays the classic pattern of a magnetic anomalous body being magnetized by induction in a zero inclination field. This is not surprising since the magnetic equator passes near the center of this body.

Description: The Hubble Space Telescope (HST) was launched April 24, 1990, and was deployed April 25 into low Earth orbit (LEO). It was soon discovered that the metal poles holding the solar arrays were expanding and contracting as the telescope orbited the Earth passing between the sunlight and the Earth s shadow. The expansion and contraction, although very small, was enough to cause the telescope to shake because of thermal-induced jitters, a detrimental effect when trying to take pictures millions of miles away. Therefore, the European Space Agency (ESA, the provider of the solar arrays) built new solar arrays (SA-11) that contained bi-stem thermal shields which insulated the solar array metal poles. These thermal shields were made of 2 mil thick aluminized-Teflon fluorinated ethylene propylene (FEP) rings fused together into a circular bellows shape. The new solar arrays were put on the HST during an extravehicular activity (EVA), also called an astronaut space walk, during the first servicing mission (SM1) in December 1993. An on-orbit photograph of the HST with the SA-11, and a close up of the bellows-like structure of the thermal shields is provided in Figure 1.

Description: Ancient geologic and hydrologic phenomena on Mars observed through the magnetic data provide windows to the ancient past through the younger Argyre and Hellas impacts, the northern plains basement and the rock materials that mantle the basement, and the Tharsis and Elysium magmatic complexes (recently referred to as superplumes). These signatures, coupled with highly degraded macrostructures (tectonic features that energetic planet during its embryonic development (0.5 Ga or so of activity) with an active dynamo and magnetosphere. One such window into the ancient past occurs northwest of the Hellas impact basin in Arabia Tern. Arabia Terra is one of the few water-rich equatorial regions of Mars, as indicated I through impact crater and elemental information. This region records many unique characteristics, including predominately Noachian materials, a highland-lowland boundary region that is distinct from other boundary regions, the presence of very few macrostructures when compared to the rest of the cratered highlands, the largest region of fretted terrain on Mars, outflow channels such as Mamers Valles that do not have obvious origins, and distinct albedo, thermal inertia, gravity, magnetic, and elemental signatures.

Description: Pavillion Lake is 5.7km long and an average of 0.8 km in width, and is located in Marble Canyon in the interior of British Columbia, Canada. It is a slightly alkaline, freshwater lake with a maximum-recorded depth of 65m. The basin walls of Pavilion Lake are lined with microbialite structures that are oriented perpendicularly to the shoreline, and which are found from depths of 5 meters to the bottom of the photic zone (light levels 1% of ambient; approximately 30m depth). These structures are speculated to have begun formation nearly 11,000 years ago, after the glacial retreat of the Cordilleran Ice Sheet. They are likely a distinctive assemblage of freshwater calcite microbialites, which display micromorphologies possibly related to the ancient Epiphyton and Girvanella classes of calcareous organosedimentary structures.

Description: The east rim of the Hellas basin and the surrounding highlands comprise a geologically significant region for evaluating volatile abundance, volatile distribution and cycling, and potential changes in Martian environmental conditions. This region of the Martian surface exhibits landforms shaped by a diversity of geologic processes and has a well-preserved geologic record, with exposures of Noachian, Hesperian, and Amazonian units, as well as spans a wide range in both latitude and elevation due to the magnitude of Hellas basin. In addition, geologically contemporaneous volcanism and volatile-driven activity in the circum-Hellas highlands provide important ingredients for creating habitats for potential Martian life.

Description: A preliminary study of the impact of the north-central Pacific circulation in the subtropical stratosphere on ozone variability locally observed by lidar is presented.

Description: Tvashtar Catena (63 N, 120 W) is one of the most interesting features on Io. This chain of large paterae (caldera-like depressions) has exhibited highly variable volcanic activity in a series of observations. Tvashtar is the type example of a style of volcanism seen only at high latitudes, with short-lived Pele-type plumes and short-lived by intense thermal events. Evidence for a hot spot at Tvashtar was first detected in an eclipse observation in April 1997 (orbit G7) by the Solid State Imager (SSI) on the Galileo Spacecraft. Tvashtar was originally targeted for observation at higher resolution in the close flyby in November 1999 (I25) because of its interesting large-scale topography. There are relatively few but generally larger paterae at high latitudes on Io. I25 images revealed a 25 km long, 1-2 km high lava curtain via a pattern of saturation and bleeding in the CCD image, which requires very high temperatures.

Description: Hot spots are manifestations of Io s mechanism of internal heating and heat transfer. Therefore, the global distribution of hot spots and their power output has important implications for how Io is losing heat. The end of the Galileo mission is an opportune time to revisit studies of the distribution of hot spots on Io, and to investigate the distribution of their power output.

Description: Ionian paterae are a class of volcanic feature that are characterized by irregular craters with steep walls, flat floors, and arcuate margins that may or may not exhibit nesting. Loki (310 W, 12 N) is Io's largest patera at approx.200 km in diameter (Figure 1), and may account for 15% of Io's total heat flow. Earth-based infrared data, as well as information collected using the Galileo Near-Infrared Mapping Spectrometer (NIMS) and the Photopolarimeter Radiometer (PPR) have been used to interpret Loki s eruption style. Debate continues over whether Loki s occasional (periodic or not) temperature increases are due to an overturning lava lake within the patera, or to an eruption of surface flows on the patera floor. Interpretation of model results and comparisons with active terrestrial lava lakes suggest that Loki behaves quite differently from active lava lakes on Earth, and that surface flows (rather than an overturning lava lake) are a more likely explanation of Loki's thermal brightening.

Description: Now that the Galileo spacecraft s tour of the Jupiter system is over, we seek to integrate all available datasets in the hopes of understanding Io as completely as possible. We have compiled information about the morphologies and locations of paterae (volcano-tectonic depressions), mountains, and hotspots on Io in a single database. It is our hope that an analysis of the spatial and temporal relationships between these features will provide more indications of the nature of the crust of Io and the mechanisms leading to these features formation. Since Io s tidal heat escapes through its crust, more knowledge about the crust will lead to an understanding of internal processes, such as magma generation and delivery to the surface, and magnitude and orientation of internal stresses.

Description: On 25 August 1992, the Microwave Limb Sounder (MLS) on the Upper Atmosphere Research Satellite observed a significant enhancement in the abundance of lower stratospheric methyl cyanide (CH3CN) at 100??hPa (~16??km altitude) in a small region off the east coast of Florida.

Description: We summarize 24 years (1978??2) of ice export estimates and examine, over a 9-year record, the associated variability in the time-varying upward-looking sonar (ULS) thickness distributions of the Fram Strait.

Description: A three-dimensional (3-D) Global Assimilative Ionospheric Model (GAIM) is currently being developed by a joint University of Southern California and Jet Propulsion Laboratory (JPL) team. To estimate the electron density on a global grid, GAIM uses a first-principles ionospheric physics model and the Kalman filter as one of its possible estimation techniques.

Description: A variance analysis technique is developed here to extract gravity wave (GW) induced temperature fluctuations from NOAA AMSU-A (Advanced Microwave Sounding Unit-A) radiance measurements. By carefully removing the instrument/measurement noise, the algorithm can produce reliable GW variances with the minimum detectable value as small as 0.1 K2. Preliminary analyses with AMSU-A data show GW variance maps in the stratosphere have very similar distributions to those found with the UARS MLS (Upper Atmosphere Research Satellite Microwave Limb Sounder). However, the AMSU-A offers better horizontal and temporal resolution for observing regional GW variability, such as activity over sub-Antarctic islands.

Description: The interplanetary shock/electric field event of 5-6 November 2001 is analyzed using ACE interplanetary data. The consequential ionospheric effects are studied using GPS receiver data from the CHAMP and SAC-C satellites and altimeter data from the TOPEX/ Poseidon satellite. Data from ~100 ground-based GPS receivers as well as Brazilian Digisonde and Pacific sector magnetometer data are also used. The dawn-to-dusk interplanetary electric field was initially ~33 mV/m just after the forward shock (IMF BZ = -48 nT) and later reached a peak value of ~54 mV/m 1 hour and 40 min later (BZ = -78 nT). The electric field was ~45 mV/m (BZ = -65 nT) 2 hours after the shock. This electric field generated a magnetic storm of intensity DST = -275 nT. The dayside satellite GPS receiver data plus ground-based GPS data indicate that the entire equatorial and midlatitude (up to +/-50(deg) magnetic latitude (MLAT)) dayside ionosphere was uplifted, significantly increasing the electron content (and densities) at altitudes greater than 430 km (CHAMP orbital altitude). This uplift peaked ~2 1/2 hours after the shock passage. The effect of the uplift on the ionospheric total electron content (TEC) lasted for 4 to 5 hours. Our hypothesis is that the interplanetary electric field ''promptly penetrated'' to the ionosphere, and the dayside plasma was convected (by E x B) to higher altitudes. Plasma upward transport/convergence led to a ~55-60% increase in equatorial ionospheric TEC to values above ~430 km (at 1930 LT). This transport/convergence plus photoionization of atmospheric neutrals at lower altitudes caused a 21% TEC increase in equatorial ionospheric TEC at ~1400 LT (from ground-based measurements). During the intense electric field interval, there was a sharp plasma ''shoulder'' detected at midlatitudes by the GPS receiver and altimeter satellites. This shoulder moves equatorward from -54(deg) to -37(deg) MLAT during the development of the main phase of the magnetic storm. We presume this to be an ionospheric signature of the plasmapause and its motion. The total TEC increase of this shoulder is ~80%. Part of this increase may be due to a "superfountain effect." The dayside ionospheric TEC above ~430 km decreased to values ~45% lower than quiet day values 7 to 9 hours after the beginning of the electric field event. The total equatorial ionospheric TEC decrease was ~16%. This decrease occurred both at midlatitudes and at the equator. We presume that thermospheric winds and neutral composition changes produced by the storm-time Joule heating, disturbance dynamo electric fields, and electric fields at auroral and subauroral latitudes are responsible for these decreases.

Description: Studies of Martian surface geomorphology and detection of near-surface water ice by the Mars Odyssey gamma ray spectrometer suggest that Mars may have had a water-rich past. While 2 to 5 wt.% of carbonate has been detected in the Martian dust [1,2], no spectral evidence for significant deposits of carbonates or sulfates has been found to date. Most investigations into Mars aqueous mineralogy have been global in scope with only a few regional studies (e.g., [3]). We are searching for localized deposits in putative lacustrine basins utilizing a basin flow model to identify basins with large drainage areas. Such basins are more likely to accumulate high concentrations of aqueous minerals than deep basins which drain only small regions.

Description: Sulfates are likely to be present on Mars as indicated by the sulfur abundances measured at the Viking and Pathfinder landing sites (approx. 5-10% by weight SO3) [1-3] and because of Mars strongly oxidizing environment. Telescopic observations of Mars tentatively identified weak sulfate bands in near infrared [4] and thermal infrared [5] data. The currently orbiting midinfrared instruments (TES, THEMIS) and the Mini-TES on the Mars Exploration Rover landers may enable a positive identification [6] and determination of the chemistry of the sulfates. Critically important to the identification of these minerals is the presence of their spectra in a spectral library. There exist approximately 370 sulfate-mineral species [7]. Sulfate minerals occur in volcanic, hydrothermal, evaporitic, and chemical-weathering environments.

Description: The global dichotomy divides the northern lowlands from the southern highlands, except where interrupted by relatively young volcanic provinces and impact basins. An elevation change of 2-4 km is typical across the dichotomy, and more than 6 km locally, over distances of several 100s km to as much as 1300 km [1,2]. A variety of exogenic and endogenic formation models have been proposed. Distinguishing between these models would help constrain the overall thermal evolution of the planet, possibly timing of core formation, and the associated mantle heat flux over time. A first step is to determine whether or not gravitational relaxation plays a role in modifying the boundary. Nimmo and Stevenson [3] examined 10 profiles across the dichotomy and used models of gravitational relaxation to conclude the relaxation has not occurred. In this study we begin by considering the geologic history in detail as inputs for modeling [4].

Description: The origin of the Martian crustal dichotomy remains a puzzle that when solved can provide an insight to the geological and geophysical evolution of Mars. In this study we model crustal relaxation in order to better constrain the original topographic shape, rheology, and temperature of the Martian crust. Our approach is to model the detailed geologic history of the Ismenius region of Mars, including slope, strain, and timing of faulting [1]. This region may contain the best preserved section of the dichotomy boundary as it is relatively unaffected by large impacts and erosion. So far the only study Martian crustal relaxation [2] suggests that the original topographic shape of the dichotomy is preserved. However, in this area strain from faulting implies at least some relaxation [1].

Description: A computational heat transfer design methodology was developed to study the dual-engine linear aerospike plume-induced base-heating environment during one power-pack out, in ascent flight. It includes a three-dimensional, finite volume, viscous, chemically reacting, and pressure-based computational fluid dynamics formulation, a special base-bleed boundary condition, and a three-dimensional, finite volume, and spectral-line-based weighted-sum-of-gray-gases absorption computational radiation heat transfer formulation. A separate radiation model was used for diagnostic purposes. The computational methodology was systematically benchmarked. In this study, near-base radiative heat fluxes were computed, and they compared well with those measured during static linear aerospike engine tests. The base-heating environment of 18 trajectory points selected from three power-pack out scenarios was computed. The computed asymmetric base-heating physics were analyzed. The power-pack out condition has the most impact on convective base heating when it happens early in flight. The source of its impact comes from the asymmetric and reduced base bleed.

Description: Numerous studies have documented the effect of El Nino-Southern Oscillation (ENSO) on rainfall in many regions of the globe. The question of whether ENSO is the single most important factor in interannual rainfall variability has received less attention, mostly because the kind of data that would be required to make such an assessment were simply not available. Until 1979 the evidence linking El Nino with changes in rainfall around the world came from rain gauges measuring precipitation over land masses and a handful of islands. From 1980 until the launch of the Tropical Rainfall Measuring Mission (TRMM) in November 1997 the remote sensing evidence was confined to ocean rainfall because of the very poor sensitivity of the instruments over land. In this paper we summarize the results of a principal component analysis of TRMM's 60-month (January 1998 to December 2002) global land and ocean remote-sensing record of monthly rainfall accumulations. Contrary to the first principal component of the rainfall itself, the first three indices of the anomaly are most sensitive to precipitation over the ocean rather than over the land. With the help of archived surface station data the first TRMM rain anomaly index is extended back several decades. Comparison of the extended index with the Southern Oscillation Index confirms that the first principal component of the rainfall anomaly is strongly correlated with the ENSO indices.

Description: This paper details a novel three-dimensional (3D) hydro-focusing micro cell sorter for micro flow cytometry applications. The unit was microfabricated by means of SU-8 3D lithography. The 3D microstructure for coaxial sheathing was designed, microfabricated, and tested. Three-dimensional hydro-focusing capability was demonstrated with an experiment to sort labeled tanned sheep erythrocytes (red blood cells). This polymer hydro-focusing microstructure is easily microfabricated and integrated with other polymer microfluidic structures.

Description: We apply the GSFC trajectory model with a series of ozonesondes to derive ozone loss rates in the lower stratosphere for the AASE-2/EASOE mission (January - March 1992) and for the SOLVE/THESEO 2000 mission (January - March 2000) in an approach similar to Match. Ozone loss rates are computed by comparing the ozone concentrations provided by ozonesondes launched at the beginning and end of the trajectories connecting the launches. We investigate the sensitivity of the Match results on the various parameters used to reject potential matches in the original Match technique and conclude that only a filter based on potential vorticity changes along the calculated back trajectory seems necessary. Our study also demonstrates that calculated ozone loss rates can vary by up to a factor of two depending upon the precise trajectory paths calculated for each trajectory. As a result an additional systematic error might need to be added to the statistical uncertainties published with previous Match results. The sensitivity to the trajectory path is particularly pronounced in the month of January, the month during which the largest ozone loss rate discrepancies between photochemical models and Match are found. For most of the two study periods, our ozone loss rates agree with those previously published. Notable exceptions are found for January 1992 at 475 K and late February/early March 2000 at 450 K, both periods during which we find less loss than the previous studies. Integrated ozone loss rates in both years compare well with those found in numerous other studies and in a potential vorticity/potential temperature approach shown previously and in this paper. Finally, we suggest an alternate approach to Match using trajectory mapping that appears to more accurately reflect the true uncertainties associated with Match and reduces the dependence upon filters that may bias the results of Match through the rejection of greater than or equal to 80% of the matched sonde pairs and 〉99% of matched observations.

Description: Since the ICRF was generated in 1995, VLBI modeling and estimation, data quality: source position stability analysis, and supporting observational programs have improved markedly. There are developing and potential applications in the areas of space navigation Earth orientation monitoring and optical astrometry from space that would benefit from a refined ICRF with enhanced accuracy, stability and spatial distribution. The convergence of analysis, focused observations, and astrometric needs should drive the production of a new realization in the next few years.

Description: The Observing Strategies Sub-group of IVS's Working Group 3 has been tasked with producing a vision for the following aspects of geodetic VLBI: antenna-network structure and observing strategies; source strength/structure/distribution; frequency bands, RFI; and field system and scheduling. These are high level considerations that have far reaching impact since they significantly influence performance potential and also constrain requirements for a number of other \VG3 sub-groups. The paper will present the status of the sub-group's work on these topics.

Description: Climate models often rely on standard atmospheres to represent various regions; these broadly capture the important physical and radiative characteristics of regional atmospheres, and become benchmarks for simulations by researchers. The high Antarctic plateau is a significant region of the earth for which such standard atmospheres are as yet unavailable. Moreover, representative profiles from atmospheres over other regions of the planet, including &om the northern high latitudes, are not comparable to the atmosphere over the Antarctic plateau, and are therefore only of limited value as substitutes in climate models. Using data from radiosondes, ozonesondes and satellites along with other observations from South Pole station, typical seasonal atmospheric profiles for the high plateau are compiled. Proper representations of rapidly changing ozone concentrations (during the ozone hole) and the effect of surface elevation on tropospheric temperatures are discussed. The differences between standard profiles developed here and the most similar standard atmosphere that already exists - namely, the Arctic Winter profile - suggest that these new profiles will be extremely useful to make accurate representations of the atmosphere over the high plateau.

Description: The Sensors and Electronics Technology Branch of the NASA Glenn Research Center is developing thin-film-based sensors for surface measurement in propulsion system research. Thin-film sensors do not require special machining of the components on which they are mounted, and they are considerably thinner than wire- or foil-based sensors. One type of sensor being advanced is the thin-film thermocouple, specifically for applications in high-temperature combustion environments. Ceramics are being demonstrated as having the potential to meet the demands of thin-film thermocouples in advanced aerospace environments. The maximum-use temperature of noble metal thin-film thermocouples, 1500 C (2700 F), may not be adequate for components used in the increasingly harsh conditions of advanced aircraft and next-generation launch vehicles. Ceramic-based thermocouples are known for their high stability and robustness at temperatures exceeding 1500 C, but are typically in the form of bulky rods or probes. As part of ASTP, Glenn's Sensors and Electronics Technology Branch is leading an in-house effort to apply ceramics as thin-film thermocouples for extremely high-temperature applications as part of ASTP. Since the purity of the ceramics is crucial for the stability of the thermocouples, Glenn's Ceramics Branch and Case Western Reserve University are developing high-purity ceramic sputtering targets for fabricating high-temperature sensors. Glenn's Microsystems Fabrication Laboratory, supported by the Akima Corporation, is using these targets to fabricate thermocouple samples for testing. The first of the materials used were chromium silicide (CrSi) and tantalum carbide (TaC). These refractory materials are expected to survive temperatures in excess of 1500 C. Preliminary results indicate that the thermoelectric voltage output of a thin-film CrSi versus TaC thermocouple is 15 times that of the standard type R (platinum-rhodium versus platinum) thermocouple, producing 20 mV with a 200 C temperature gradient. The photograph on the left shows the CrSi-TaC thermocouple in a test fixture at Glenn, and the resulting output signal is shown on the right. The temperature differential across the sample, from the center of the sample inside the oven to the sample mount outside the oven, is measured using a type R thermocouple on the sample.

Description: The unsteady, incompressible Navier-Stokes equations are used for the direct numerical simulation (DNS) of spatially evolving disturbances in a three-dimensional (3-D) attachment-line boundary layer. Two-dimensional (2-D) disturbances are introduced either by forcing at the in ow or by harmonic-source generators at the wall; 3-D disturbances are introduced by harmonic-source generators at the wall. The DNS results are in good agreement with both 2-D non-parallel theory (for small-amplitude disturbances) and weakly nonlinear theory (for finite-amplitude disturbances), which validates the two theories. The 2-D DNS results indicate that nonlinear disturbance growth occurs near branch II of the neutral stability curve; however, steady suction can be used to stabilize this disturbance growth. For 3-D instabilities that are generated o the attachment line, spreading both toward and away from the attachment line causes energy transfer to the attachment-line and downstream instabilities; suction stabilizes these instabilities. Furthermore, 3-D instabilities are more stable than 2-D or quasi-2-D instabilities.

Description: A fast multigrid solver for the steady incompressible Euler equations is presented. Unlike time-marching schemes, this approach uses relaxation of the steady equations. Application of this method results in a discretization that correctly distinguishes between the advection and elliptic parts of the operator, allowing efficient smoothers to be constructed. Solvers for both unstructured triangular grids and structured quadrilateral grids have been written. Computations for channel flow and flow over a nonlifting airfoil have computed. Using Gauss-Seidel relaxation ordered in the flow direction, textbook multigrid convergence rates of nearly one order-of-magnitude residual reduction per multigrid cycle are achieved, independent of the grid spacing. This approach also may be applied to the compressible Euler equations and the incompressible Navier-Stokes equations.

Description: The structural analysis results for a graphite/epoxy quasi-isotropic circular plate subjected to a forced rotation at the boundary and pressure is presented. The analysis is to support a specialized material characterization test for composite cryogenic tanks. Finite element models were used to ensure panel integrity and determine the pressure necessary to achieve a predetermined equal biaxial strain value. The displacement results due to the forced rotation at the boundary led to a detailed study of the bending stiffness matrix [D]. The variation of the bending stiffness terms as a function of angular position is presented graphically, as well as, an illustrative technique of considering the laminate as an I-beam.

Description: Various artificial compressibility methods for calculating the three-dimensional incompressible Navier-Stokes equations are compared. Each method is described and numerical solutions to test problems are conducted. A comparison based on convergence behavior, accuracy, and robustness is given.

Description: This study provides systematic method for reducing power consumption in reduced gravity systems by adopting minimum velocity required to provide adequate CHF and preclude detrimental effects of reduced gravity . This study proves it is possible to use existing 1 ge flow boiling and CHF correlations and models to design reduced gravity systems provided minimum velocity criteria are met

Description: Recent workshops to define strategic research on critical issues in microgravity fluids and transport phenomena in support of mission orientated needs of NASA and many technical conferences over the years in support of fundamental research targeting NASA's long range missions goal have identified several phase change processes needed to design advanced space and planetary based systems for long duration operations Recommendation noted that phase change processes are profoundly affected by gravitational environment.

Description: We present experimental data on flow pattern transitions, pressure drop and flow characteristics for cocurrent gas-liquid flow through packed columns in microgravity. The flow pattern transition data indicates that the pulse flow regime exists over a wider range of gas and liquid flow rates under microgravity conditions compared to 1-g and the widely used Talmor map in 1-g is not applicable for predicting the transition boundaries. A new transition criterion between bubble and pulse flow in microgravity is proposed and tested using the data. Since there is no static head in microgravity, the pressure drop measured is the true frictional pressure drop. The pressure drop data, which has much smaller scatter than most reported 1-g data clearly shows that capillary effects can enhance the pressure drop (especially in the bubble flow regime) as much as 200% compared to that predicted by the single phase Ergun equation. The pressure drop data are correlated in terms of a two-phase friction factor and its dependence on the gas and liquid Reynolds numbers and the Suratman number. The influence of gravity on the pulse amplitude and frequency is also discussed and compared to that under normal gravity conditions. Experimental work is planned to determine the gas-liquid and liquid-solid mass transfer coefficients. Because of enhanced interfacial effects, we expect the gas-liquid transfer coefficients kLa and kGa (where a is the gas-liquid interfacial area) to be higher in microgravity than in normal gravity at the same flow conditions. This will be verified by gas absorption experiments, with and without reaction in the liquid phase, using oxygen, carbon dioxide, water and dilute aqueous amine solutions. The liquid-solid mass transfer coefficient will also be determined in the bubble as well as the pulse flow regimes using solid benzoic acid particles in the packing and measuring their rate of dissolution. The mass transfer coefficients in microgravity will be compared to those in normal gravity cocurrent flow to determine the mass transfer enhancement and propose new mass transfer correlations for two-phase gas-liquid flows through packed beds in microgravity.

Description: One particular characteristic observed in unsteady shear layers is the phase shift relative to the main flow. In attached boundary layers this will have an effect both on the instantaneous skin friction and heat transfer. In separation bubbles the contribution to the drag is dominated by the pressure distribution. However, the most significant effect appears to be the phase shift on the transition process. Unsteady transition behaviour may determine the bursting of the bubble resulting in an un-recoverable full separation. An early analysis of the phase shift was performed by Stokes for the incompressible boundary layer of an oscillating wall and an oscillating main flow. An amplitude overshoot within the shear layer as well as a phase shift were observed that can be attributed to the relatively slow diffusion of viscous stresses compared to the fast change of pressure. Experiments in a low speed facility with the boundary layer of a flat plate were evaluated in respect to phase shift. A pressure distribution similar to that on the suction surface of a turbomachinery aerofoil was superimposed generating a typical transitional separation bubble. A periodically unsteady main flow in the suction type wind tunnel was introduced via a rotating flap downstream of the test section. The experiments covered a range of the three similarity parameters of momentum-loss-thickness Reynolds-number of 92 to 226 and Strouhal-number (reduced frequency) of 0.0001 to 0.0004 at the separation point, and an amplitude range up to 19 %. The free stream turbulence level was less than 1% .Upstream of the separation point the phase shift in the laminar boundary layer does not appear to be affected significantly bay either of the three parameters. The trend perpendicular to the wall is similar to the Stokes analysis. The problem scales well with the wave velocity introduced by Stokes, however, the lag of the main flow near the wall is less than indicated analytically. The separation point comes closest to the Stokes analysis but the phase is still 20 degrees lower at the wall.

Description: Among the numerous causes for unsteadiness in turbo machinery flows are turbulence and flow environment, wakes from stationary and rotating vanes, boundary layer separation, boundary layer/shear layer instabilities, presence of shock waves and deliberate unsteadiness for flow control purposes. These unsteady phenomena may lead to flow-structure interactions such as flutter and forced vibration as well as system instabilities such as stall and surge. A major issue of unsteadiness relates to the fact that a fundamental understanding of unsteady flow physics is lacking and requires continued attention. Accurate simulations and sufficient high fidelity experimental data are not available. The Glenn Research Center plan for Engine Component Flow Physics Modeling is part of the NASA 21st Century Aircraft Program. The main components of the plan include Low Pressure Turbine National Combustor Code. The goals, technical output and benefits/impacts of each element are described in the presentation. The specific areas selected for discussion in this presentation are blade wake interactions, flow control, and combustor exit turbulence and modeling.

Description: The Effects of elevated free-stream turbulence (FST) on the natural and periodically excited separation bubbles were studied experimentally, due to the relevance of this flow to low-pressure turbine blades at low Reynolds numbers. A bubble was formed at the leading edge of a flat plate and the FST level was altered by placing a grid across the flow at different locations upstream of the plate. The mixing across the separated shear-layer, forming the free boundary of the bubble, increased due to the elevated FST and due to nominally two-dimensional periodic excitation, both flattening and shortening the bubble. Periodic excitation at frequencies that were at least an order of magnitude lower than those associated with the initial shear-layer instability, were very effective at low FST, because the amplitudes of the excitation frequency and its harmonic were amplified over the bubble. High frequency excitation (F+ 3, based on the length of the baseline low FST bubble) had a major effect close to the separation location, while farther downstream the excited fluctuations rapidly decayed in the reattachment region. Low frequency excitation, that generated waves comparable to the length of the unperturbed bubble (F+ 1) were less effective and their magnitude decayed at a slower rate downstream of reattachment. An increase in the level of the FST reduced the net effect of the periodic excitation on the mixing enhancement and subsequent reattachment process, probably due to a destructive interference between the nominally 2D excitation and the random (in space and time) FST, reducing the spanwise coherence and therefore the effectiveness of the current control strategy. However, even at the reduced effectiveness of 2D periodic excitation at elevated FST, it accelerated the reattachment process and the recovery rate of the reattached boundary layer, enhancing the boundary layer resistance to repeat separation and reducing its momentum loss further downstream.

Description: In this brief, we investigate the turbulent condensation of a population of droplets by means of a direct numerical simulation. To that end, a coupled Navier-Stokes/Lagrangian solver is used where each particle is tracked and its growth by water vapor condensation is monitored exactly. The main goals of the study are to find out whether turbulence broadens the droplet size distribution, as observed in in situ measurements. The second issue is to understand if and for how long a correlation between the droplet radius and the local supersaturation exists for the purpose of modeling sub-grid scale microphysics in cloud-resolving codes. This brief is organized as follows. In Section 2 the governing equations are presented, including the droplet condensation model. The implementation of the forcing procedure is described in Section 3. The simulation results are presented in Section 4 together with a sketch of a simple stochastic model for turbulent condensation. Conclusions and the main outcomes of the study are given in Section 5.

Description: Three-dimensional computational techniques, in particular the uncoupled CFD-DSMC of the present study, are available to be applied to problems such as jet interactions with variable density regions ranging from a continuum jet to a rarefied free stream. When the value of the jet to free stream momentum flux ratio approximately greater than 2000 for a sharp leading edge flat plate forward separation vortices induced by the jet interaction are present near the surface. Also as the free stream number density n (infinity) decreases, the extent and magnitude of normalized pressure increases and moves upstream of the nozzle exit. Thus for the flat plate model the effect of decreasing n (infinity) is to change the sign of the moment caused by the jet interaction on the flat plate surface.

Description: InSPACE is a microgravity fluid physics experiment that was operated on the International Space Station (ISS) in the Microgravity Science Glovebox from late March 2003 through early July 2003. (InSPACE is an acronym for Investigating the Structure of Paramagnetic Aggregates From Colloidal Emulsions.) The purpose of the experiment is to obtain fundamental data of the complex properties of an exciting class of smart materials termed magnetorheological (MR) fluids. MR fluids are suspensions, or colloids, comprised of small (micrometer-sized) superparamagnetic particles in a nonmagnetic medium. Colloids are suspensions of very small particles suspended in a liquid. (Examples of other colloids are blood, milk, and paint.) These controllable fluids can quickly transition into a nearly solid state when exposed to a magnetic field and return to their original liquid state when the magnetic field is removed. Controlling the strength of the magnetic field can control the relative stiffness of these fluids. MR fluids can be used to improve or develop new seat suspensions, robotics, clutches, airplane landing gear, and vibration damping systems. The principal investigator for InSPACE is Professor Alice P. Gast of the Massachusetts Institute of Technology (MIT). The InSPACE hardware was developed at the NASA Glenn Research Center. The InSPACE samples were delivered to the ISS in November 2002, on the Space Shuttle Endeavour, on Space Station Utilization Flight UF-2/STS113. Operations began on March 31, 2003, with the processing of three different particle size samples at multiple test parameters. This investigation focused on determining the structural organization of MR colloidal aggregates when exposed to a pulsing magnetic field. On Earth, the aggregates take the shape of footballs with spiky tips. This characteristic shape may be influenced by the pull of gravity, which causes most particles initially suspended in the fluid to sediment, (i.e., settle and collect at the bottom of the cell). In the absence of sedimentation effects on the ISS, the behavior and shape of these MR aggregate structures are dominated exclusively by magnetic and surface tension forces. The microscopic detail of these structures was imaged under two orthogonal camera views. The video was downlinked to the InSPACE team at Glenn's Telescience Support Center and to MIT and also recorded onboard the ISS on videotapes that will be brought back to the ground by the space shuttles. The study examined the effect on the structure formation by varying the magnetic field strength and pulse frequency, and particle size. Fundamental data that characterized the structure formation were obtained. InSPACE completed its last planned test run on July 2, 2003. Operations occurred on 21 days over approximately a 3-month period. Forty-one test points were completed during 26 test runs. During the initial testing, the procedures followed by the crew were modified to maximize the observation of some unexpected and interesting aggregate behavior. As a result Dr. Gast has reported on the formation of aggregate shapes that are more extended and diverse than those observed on the ground. Sheets of magnetic material folded over in a labyrinth pattern and large columnar aggregates with complex interfaces with the surrounding fluid are examples of the interesting structures that have been observed on the ISS. In light of these early findings, the understanding of the fundamental properties of MR fluids on the basis of ground-based observations may need to be reconsidered.The experiments on the ISS have provided a vast amount of video data for analysis. While this analysis is ongoing, plans are being made for additional experimental runs. For this purpose, additional hardware and cells containing samples of different magnetic particles and sizes are being fabricated for a future launch to the ISS. The InSPACE hardware will remain on orbit until this testing is completed.

Description: Leakage and wear are two fundamental problems in all traditional turbine seals that contribute to an engine's inefficiency. The solutions to seal leakage and wear conflict in the conventional design space. Reducing the clearance between the seal and rotating shaft reduces leakage but increases wear because of increased contact incidents. Increasing the clearance to reduce the contact between parts reduces wear but increases parasitic leakage. The goal of this effort is to develop a seal that restricts leakage flow using acoustic pressure while operating in a noncontacting manner, thereby increasing life. In 1996, Dr. Timothy Lucas announced his discovery of a method to produce shock-free high-amplitude pressure waves. For the first time, the formation of large acoustic pressures was possible using dissonant resonators. A pre-prototype acoustic seal developed at the NASA Glenn Research Center exploits this fundamental acoustic discovery: a specially shaped cavity oscillated at the contained fluid's resonant frequency produces high-amplitude acoustic pressure waves of a magnitude approaching those required of today's seals. While the original researchers are continuing their development of acoustic pumps, refrigeration compressors, and electronic thermal management systems using this technology, the goal of researchers at Glenn is to apply these acoustic principles to a revolutionary sealing device. When the acoustic resonator shape is optimized for the sealing device, the flow from a high-pressure cavity to a low-pressure cavity will be restricted by a series of high-amplitude standing pressure waves of higher pressure than the pressure to be sealed. Since the sealing resonator cavity will not touch the adjacent sealing structures, seal wear will be eliminated, improving system life. Under a cooperative agreement between Glenn and the Ohio Aerospace Institute (OAI), an acoustic-based pre-prototype seal was demonstrated for the first time. A pressurized cavity was attached to one end of the resonator while the other end remained open to ambient pressure. Measurements were taken at several values of applied pressure with the assembly stationary, oscillated at an off-resonance frequency, and then oscillated on-resonance. The three cases show that the flow through the conical resonator can be reduced by oscillating the resonator at the resonance frequency of the air contained within the cavity. The results are currently being compared with results obtained from a commercial computational fluid dynamics code. The objective is to improve the design through numerical simulation before fabricating a next-generation prototype sealing device. Future work is aimed at implementing acoustic seal design improvements to further reduce the leakage flow rate through the device and at reducing the device's overall size.

Description: The VLBI antenna (GILCREEK) at Fairbanks, Alaska observes in networks routinely twice a week with operational networks and on additional days with other networks on a more uneven basis. The Fairbanks antenna position is about 150 km north of the Denali fault and from the earthquake epicenter. We examine the transient behavior of the estimated VLBI position during the year following the earthquake to determine how the rate of change of postseismic deformation has changed. This is compared with what is seen in the GPS site position series.

Description: Whereas the surface units of the northern plain of Mars generally exhibit ages ranging from late Hesperian to Amazonian, interpretation of precise topographic measurements indicate that the age of the underlying "basement" is early Noachian, or almost as old as the southern highlands. This suggests that widespread but relatively superficial resurfacing has occurred throughout the northern plains since the end of early heavy bombardment. In this abstract I examine some of the possible implications of the subsurface structure inferred for the Utopia basin from gravity data on the nature of this resurfacing. The large, shallow, circular depression in Utopia Planitia has been identified as a huge impact basin, based on both geological evidence and detailed analysis of MOLA topography. Its diameter (approx. 3000 km) is equivalent to that of the Hellas basin, as is its inferred age (early Noachian). However, whereas Hellas is extremely deep with rough terrain and large slopes, the Utopia basin is a smooth, shallow, almost imperceptible bowl. Conversely, Utopia displays one of the largest (non-Tharsis-related) positive geoid anomalies on Mars, in contrast to a much more subdued negative anomaly over Hellas.

Description: For the years 1999-2003, we estimate the time-varying perennial ice zone (PIZ) coverage and construct the annual cycles of multiyear (MY, including second year) ice coverage of the Arctic Ocean using QuikSCAT backscatter, MY fractions from RADARSAT, and the record of ice export from satellite passive microwave observations. An area balance approach extends the winter MY coverage from QuikSCAT to the remainder of the year. From these estimates, the coverage of MY ice at the beginning of each year is 3774 x 10(exp 3) sq km (2000), 3896 x 10(exp 3) sq km (2001), 4475 x 10(exp 3) sq km (2002), and 4122 x 10(exp 3) sq km (2003). Uncertainties in coverage are approx.150 x 10(exp 3) sq km. In the mean, on 1 January, MY ice covers approx.60% of the Arctic Ocean. Ice export reduces this coverage to approx.55% by 1 May. From the multiple annual cycles, the area of first-year (FY) ice that survives the intervening summers are 1192 x 10(exp 3) sq km (2000), 1509 x 10(exp 3) sq km (2001), and 582 x 10(exp 3) sq km (2002). In order for the MY coverage to remain constant from year to year, these replenishment areas must balance the overall area export and melt during the summer. The effect of the record minimum in Arctic sea ice area during the summer of 2002 is seen in the lowest area of surviving FY ice of the three summers. In addition to the spatial coverage, the location of the PIZ is important. One consequence of the unusual location of the PIZ at the end of the summer of 2002 is the preconditioning for enhanced export of MY ice into the Barents and Kara seas. Differences between the minimums in summer sea ice coverage from our estimates and passive microwave observations are discussed.

Description: Bi-directional reflectances of marine liquid water clouds, as measured by the Multiangle Imaging SpectroRadiometer (MISR), are compared with plane-parallel radiative transfer model calculations.

Description: The Crustal Dynamics Data Information System (CDDIS) has supported the International GPS Service (IGS) as a global data center since 1992. The CDDIS activities within the IGS during 2001 are summarized below; this report also includes any changes or enhancements made to the CDDIS during the past year. General CDDIS background and system information can be found in the CDDIS data center summary included in the IGS 1994 Annual Report (Noll, 1995) as well as the subsequent updates (Noll, 1996, Noll, 1997, Noll, 1998, Noll, 1999, and Noll, 2001).

Description: Temperature-Salinity (T-S) relationship variability in the pycnocline of the eastern equatorial Pacific Ocean (NINO3 region, 5 degrees S ??degrees N, 150 degrees W ?? degrees W) over the last two decades is investigated using observational data and model simulation.

Description: We show that the reconstructed sensitivity of the sea level temperature to long term solar forcing in the Northern Hemisphere is in very good agreement with the empirical temperature pattern corresponding to changes of the North Annular Mode (NAM).

Description: A fluid, mobile atmosphere and oceans surrounds the solid Earth and upon its land surface lays a continually changing distribution of ice, snow, and ground water.

Description: A new technology for reducing turbulent skin friction, called the Microblowing Technique (MBT), is presented. Results from proof-of-concept experiments show that this technology could potentially reduce turbulent skin friction by more than 50% of the skin friction of a solid flat plate for subsonic and supersonic flow conditions. The primary purpose of this review paper is to provide readers with information on the turbulent skin friction reduction obtained from many experiments using the MBT. Although the MBT has a penalty for obtaining the microblowing air associated with it, some combinations of the MBT with suction boundary layer control methods are an attractive alternative for a real application. Several computational simulations to understand the flow physics of the MBT are also included. More experiments and computational fluid dynamics (CFD) computations are needed for the understanding of the unsteady flow nature of the MBT and the optimization of this new technology.

Description: Solve the time-dependent inviscid flow equations for this geometry subject to the specified inflow/outflow mean conditions and the fluctuating inflow velocity distortion. (1) Compute the unsteady solution until periodicity in pressure is achieved by showing that at least two successive periods are identical. Periodicity must be achieved on both the airfoil surface and the inflow/outflow boundaries. (2) Once periodicity is achieved, compute the pressure frequency spectra on the reference airfoil on both the upper and lower surfaces at x=(-0.25c,0.00, +0.25c), on the inflow boundary at (x,y)={1.5c,-0.3c), (-1.5c,0.0),(-1.5c,0.3c)} and on the outflow boundary at (x,y)= {(1.5c,-0.3c),(1.5c,0.0), (1.5c,0.3c)}. Express the spectral results in dB using the standard definition 20 log(P(sub(r.m.s)/P(sub ref), where p(sub ref) == 20 microPa. (3) Extract the harmonic pressure distributions on the inflow and outflow boundaries (i.e., on x= -/+ 1.5c lines) at the fundamental frequency omega and apply a Fourier transform in y direction to identify the spatial (i.e., mode order) structure of the pressure perturbations. Express the result in dB for each mode order. Repeat the process for the frequencies 2omega and 3omega.

Description: Three benchmark problems from the current and previous CAA workshops involving tone noise generated in viscous flows are investigated using the CE/SE finite volume method. The CE/SE method is first briefly reviewed. Then, the benchmark problems, namely, flow past a single cylinder (CAA Workshop II problem), flow past twin cylinders (from the current CAA Workshop IV, Category 5, Problem 1) and flow past a deep cavity with overhang (CAA Workshop III problem) are investigated. Generally good results are obtained in comparison with the experimental data.

Description: The effects of gravity on the bubble formation in an annular jet were studied. The experiments were conducted in the 2.2-second drop tower at the NASA Glenn Research Center. Terrestrial gravity experiments were conducted at the Fluid Dynamics Research Laboratory at the University of Oklahoma. Stainless steel tubing with inner diameters of 1/8" (gas inner annulus) and 5/16" (liquid outer annulus) served as the injector. A rectangular test section, 6" x 6" x 14" tall, made out of half-inch thick Lexan was used. Images of the annular jet were acquired using a high-speed camera. The effects of gravity and varying liquid and gas flow rates on bubble size, wavelength, and breakup length were documented. In general, the bubble diameter was found to be larger in terrestrial gravity than in microgravity for varying Weber numbers (0.05 - 0.16 and 5 - 11) and liquid flow rates (1.5 ft/s - 3.0 ft/s). The wavelength was found to be larger in terrestrial gravity than in microgravity, but remained constant for varying Weber numbers. For low Weber numbers (0.05 - 0.16), the breakup length in microgravity was significantly higher than in terrestrial gravity. Comparison with linear stability analysis showed estimated bubble sizes within 9% of experimental bubble sizes. Bubble size compared to other terrestrial gravity experiments with same flow conditions showed distinct differences in bubble size, which displayed the importance of injector geometry on bubble formation.

Description: The objective of the present study is to develop a two-fluid model formulation with interfacial area transport equation applicable for microgravity conditions. The new model is expected to make a leapfrog improvement by furnishing the constitutive relations for the interfacial interaction terms with the interfacial area transport equation, which can dynamically model the changes of the interfacial structures. In the first year of this three-year project supported by the U.S. NASA, Office of Biological and Physics Research, the primary focus is to design and construct a ground-based, microgravity two-phase flow simulation facility, in which two immiscible fluids with close density will be used. In predicting the two-phase flow behaviors in any two-phase flow system, the interfacial transfer terms are among the most essential factors in the modeling. These interfacial transfer terms in a two-fluid model specify the rate of phase change, momentum exchange, and energy transfer at the interface between the two phases. For the two-phase flow under the microgravity condition, the stability of the fluid particle interface and the interfacial structures are quite different from those under normal gravity condition. The flow structure may not reach an equilibrium condition and the two fluids may be loosely coupled such that the inertia terms of each fluid should be considered separately by use of the two-fluid model. Previous studies indicated that, unless phase-interaction terms are accurately modeled in the two-fluid model, the complex modeling does not necessarily warrant an accurate solution.

Description: Flow regime and pressure drop data was obtained and analyzed. Pulse flow exists at lower liquid flow rates in 0-g compared to 1-g. 1-g flow regime maps do not apply in microgravity. Pressure drop is higher in microgravity (enhanced interfacial effects).

Description: In fiscal year 2003, test cell 23 of the Research Combustion Laboratory (RCL 23) at the NASA Glenn Research Center was upgraded with the addition of gaseous hydrogen as a working propellant and the addition of a 450-psig air-supply system. Test flexibility was further enhanced by upgrades to the facility control systems. RCL 23 can now test with gaseous hydrogen flow rates up to 0.05 lbm/sec and jet fuel flow rates up to 0.62 lbm/sec. Research airflow rates up to 3 lbm/sec are possible with the 450-psig supply system over a range of inlet temperatures. Nonvitiated, heated air is supplied from a shell and tube heat exchanger. The maximum nonvitiated facility air temperature is 1100 F at 1.5 lbm/sec. Research-section exhaust temperatures are limited to 3200 F because of material and cooling capacity limits. A variety of support systems are available depending on the research hardware configuration. Test section ignition can be provided via either a hydrogen air torch system or an electronic spark system. Emissions measurements are obtained with either pneumatically or electromechanically actuated gas sample probes, and the electromechanical system allows for radial measurements at a user-specified axial location for measurement of emissions profiles. Gas analysis data can be obtained for a variety of species, including carbon monoxide (CO), carbon dioxide (CO2), nitrogen oxides (NO and NOx), oxygen (O2), unburnt hydrocarbons, and unburnt hydrogen. Facility control is accomplished with a programmable logic control system. Facility operations have been upgraded to a system based on graphical user interface control screens. A data system is available for real-time acquisition and monitoring of both measurements in engineering units and performance calculations. The upgrades have made RCL 23 a highly flexible facility for research into low emissions gas turbine combustor concepts, and the flame tube configuration inherently allows for a variety of fuel nozzle configurations to be tested in a cost-effective manner. RCL 23 is poised to be a leading facility for developing modern low-emission fuel nozzles for use with jet fuel and alternative fuels.

Description: The operation of a packed bed reactor (PBR) involves gas and liquid flowing simultaneously through a fixed-bed of solid particles. Depending on the application, the particles can be various shapes and sizes but are generally designed to force the two fluid phases through a tortuous route of narrow channels connecting the interstitial space. The PBR is the most common type of reactor in industry because it provides for intimate contact and high rates of transport between the phases needed to sustain chemical or biological reactions. The packing may also serve as either a catalyst or as a support for growing biological material. Furthermore, this type of reactor is relatively compact and requires minimal power to operate. This makes it an excellent candidate for unit operations in support of long-duration human space activities.

Description: Recent research efforts within the Microgravity Fluid Physics Branch of the NASA Glenn Research Center have necessitated the development of a microscope capable of high-resolution, three-dimensional imaging of intracellular structure and tissue morphology. Standard optical microscopy works well for thin samples, but it does not allow the imaging of thick samples because of severe degradation caused by out-of-focus object structure. Confocal microscopy, which is a laser-based scanning microscopy, provides improved three-dimensional imaging and true optical sectioning by excluding the out-of-focus light. However, in confocal microscopy, out-of-focus object structure is still illuminated by the incoming beam, which can lead to substantial photo-bleaching. In addition, confocal microscopy is plagued by limited penetration depth, signal loss due to the presence of a confocal pinhole, and the possibility of live-cell damage. Two-photon microscopy is a novel form of laser-based scanning microscopy that allows three-dimensional imaging without many of the problems inherent in confocal microscopy. Unlike one-photon microscopy, it utilizes the nonlinear absorption of two near-infrared photons. However, the efficiency of two-photon absorption is much lower than that of one-photon absorption because of the nonlinear (i.e., quadratic) electric field dependence, so an ultrafast pulsed laser source must typically be employed. On the other hand, this stringent energy density requirement effectively localizes fluorophore excitation to the focal volume. Consequently, two-photon microscopy provides optical sectioning and confocal performance without the need for a signal-limiting pinhole. In addition, there is a reduction in photo-damage because of the longer excitation wavelength, a reduction in background fluorescence, and a 4 increase in penetration depth over confocal methods because of the reduction in Rayleigh scattering.

Description: Many valuable advances in combustion science have come from observations of microgravity flames. This research is contributing to the improved efficiency and reduced emissions of practical combustors and is benefiting terrestrial and spacecraft fire safety. Unfortunately, difficulties associated with microgravity have prevented many types of measurements in microgravity flames. In particular, temperature measurements in flames are extremely important but have been limited in microgravity. A novel method of measuring temperatures in microgravity flames is being developed in-house at the National Center for Microgravity Research and the NASA Glenn Research Center and is described here. Called thin-filament pyrometry, it involves using a camera to determine the local gas temperature from the intensity of inserted fibers glowing in a flame. It is demonstrated here to provide accurate measurements of gas temperatures in a flame simultaneously at many locations. The experiment is shown. The flame is a laminar gas jet diffusion flame fueled by methane (CH4) flowing from a 14-mm round burner at a pressure of 1 atm. A coflowing stream of air is used to prevent flame flicker. Nine glowing fibers are visible. These fibers are made of silicon carbide (SiC) and have a diameter of 15 m (for comparison, the average human hair is 75 m in diameter). Because the fibers are so thin, they do little to disturb the flame and their temperature remains close to that of the local gas. The flame and glowing filaments were imaged with a digital black-and-white video camera. This camera has an imaging area of 1000 by 1000 pixels and a wide dynamic range of 12 bits. The resolution of the camera and optics was 0.1 mm. Optical filters were placed in front of the camera to limit incoming light to 750, 850, 950, and 1050 nm. Temperatures were measured in the same flame in the absence of fibers using 50-m Btype thermocouples. These thermocouples provide very accurate temperatures, but they generally are not useful in microgravity tests because they measure temperature at only one location at a time. Thermocouple measurements at a height of 11 mm above the burner were used to calibrate the thin-filament pyrometry system at all four wavelengths. This calibration was used to perform thin-filament pyrometry at other heights above the burner. One such profile is shown in this graph; this is for a height of 21 mm. The agreement between the pyrometry measurements and thermocouple results at this height is excellent in the range of 1000 to 2000 K, with an estimated uncertainty of 50 K and an estimated upper limit of 2500 K. Neither the thermocouple nor the thin-filament pyrometry temperatures have been corrected for radiation, but the correction is expected to be nearly the same for both methods. We anticipate that thin-filament pyrometry similar to that developed here will become an important diagnostic for studies of microgravity flames owing to its accuracy and its ability to simultaneously measure finely spaced temperatures.

Description: When discharged into an International Space Station (ISS) payload rack, a carbon dioxide (CO2) portable fire extinguisher (PFE) must extinguish a fire by decreasing the oxygen in the rack by 50 percent within 60 sec. The length of time needed for this oxygen reduction throughout the rack and the length of time that the CO2 concentration remains high enough to prevent the fire from reigniting is important when determining the effectiveness of the response and postfire procedures. Furthermore, in the absence of gravity, the local flow velocity can make the difference between a fire that spreads rapidly and one that self-extinguishes after ignition. A numerical simulation of the discharge of CO2 from PFE into the Combustion Integrated Rack (CIR) in microgravity was performed to obtain the local velocity and CO2 concentration. The complicated flow field around the PFE nozzle exits was modeled by sources of equivalent mass and momentum flux at a location downstream of the nozzle. The time for the concentration of CO2 to reach a level that would extinguish a fire anywhere in the rack was determined using the Fire Dynamics Simulator (FDS), a computational fluid dynamics code developed by the National Institute of Standards and Technology specifically to evaluate the development of a fire and smoke transport. The simulation shows that CO2, as well as any smoke and combustion gases produced by a fire, would be discharged into the ISS cabin through the resource utility panel at the bottom of the rack. These simulations will be validated by comparing the results with velocity and CO2 concentration measurements obtained during the fire suppression system verification tests conducted on the CIR in March 2003. Once these numerical simulations are validated, portions of the ISS labs and living areas will be modeled to determine the local flow conditions before, during, and after a fire event. These simulations can yield specific information about how long it takes for smoke and combustion gases produced by a fire to reach a detector location, how large the fire would be when the detector alarms, and the behavior of the fire until it has been extinguished. This new capability could then be used to optimize the location of fire detectors and fire-suppression ports as well as to evaluate the effectiveness of fire suppressants and response strategies. Numerical data collected from these simulations could also be used to develop a virtual reality fire event for crew training and fire safety awareness. This work is funded by NASA's Bioastronautics Initiative, which has the objective of ensuring and enhancing the health, safety, and performance of humans in space. As part of this initiative, the Microgravity Combustion Science Branch at the NASA Glenn Research Center is conducting spacecraft fire safety research to significantly improve fire safety on inhabited spacecraft.

Description: A fast-running unsteady aerodynamics code, LINFLUX, was previously developed for predicting turbomachinery flutter. This linearized code, based on a frequency domain method, models the effects of steady blade loading through a nonlinear steady flow field. The LINFLUX code, which is 6 to 7 times faster than the corresponding nonlinear time domain code, is suitable for use in the initial design phase. Earlier, this code was verified through application to a research fan, and it was shown that the predictions of work per cycle and flutter compared well with those from a nonlinear time-marching aeroelastic code, TURBO-AE. Now, the LINFLUX code has been applied to real configurations: fans developed under the Energy Efficient Engine (E-cubed) Program and the Quiet Aircraft Technology (QAT) project. The LINFLUX code starts with a steady nonlinear aerodynamic flow field and solves the unsteady linearized Euler equations to calculate the unsteady aerodynamic forces on the turbomachinery blades. First, a steady aerodynamic solution is computed for given operating conditions using the nonlinear unsteady aerodynamic code TURBO-AE. A blade vibration analysis is done to determine the frequencies and mode shapes of the vibrating blades, and an interface code is used to convert the steady aerodynamic solution to a form required by LINFLUX. A preprocessor is used to interpolate the mode shapes from the structural dynamics mesh onto the computational fluid dynamics mesh. Then, LINFLUX is used to calculate the unsteady aerodynamic pressure distribution for a given vibration mode, frequency, and interblade phase angle. Finally, a post-processor uses the unsteady pressures to calculate the generalized aerodynamic forces, eigenvalues, an esponse amplitudes. The eigenvalues determine the flutter frequency and damping. Results of flutter calculations from the LINFLUX code are presented for (1) the E-cubed fan developed under the E-cubed program and (2) the Quiet High Speed Fan (QHSF) developed under the Quiet Aircraft Technology project. The results are compared with those obtained from the TURBO-AE code. A graph of the work done per vibration cycle for the first vibration mode of the E-cubed fan is shown. It can be seen that the LINFLUX results show a very good comparison with TURBO-AE results over the entire range of interblade phase angle. The work done per vibration cycle for the first vibration mode of the QHSF fan is shown. Once again, the LINFLUX results compare very well with the results from the TURBOAE code.

Description: A new in-house test capability has been developed at the NASA Glenn Research Center, where a critical component of the Stirling Radioisotope Generator (SRG) is undergoing extensive testing to aid the development of analytical life prediction methodology and to experimentally aid in verification of the flight-design component's life. The new facility includes two test rigs that are performing creep testing of the SRG heater head pressure vessel test articles at design temperature and with wall stresses ranging from operating level to seven times that (see the following photograph).

Description: The TopMaker technique was developed in an effort to reduce the time required for grid generation in complex numerical studies. Topology generation accounts for much of the man-hours required for structured multiblock grids. With regard to structured multiblock grids, topology refers to how the blocks are arranged and connected. A two-dimensional multiblock topology generation technique has been developed at the NASA Glenn Research Center. Very general configurations can be addressed by the technique. A configuration is defined by a collection of non-intersecting closed curves, which will be referred to as loops. More than a single loop implies that holes exist in the domain, which poses no problem. This technique requires only the medial vertices and the touch points that define each vertex. From the information about the medial vertices, the connectivity between medial vertices is generated. The physical shape of the medial edge is not required. By applying a few simple rules to each medial edge, a multiblock topology can be generated without user intervention. The resulting topologies contain only the level of complexity dictated by the configurations. Grid lines remain attached to the boundary except at sharp concave turns, where a change in index family is introduced as would be desired. Keeping grid lines attached to the boundary is especially important in computational fluid dynamics, where highly clustered grids are used near no-slip boundaries. This technique is simple and robust and can easily be incorporated into the overall grid-generation process.

Description: Over the past thirty years, numerical methods and simulation tools for incompressible flows have been advanced as a subset of the computational fluid dynamics (CFD) discipline. Although incompressible flows are encountered in many areas of engineering, simulation of compressible flow has been the major driver for developing computational algorithms and tools. This is probably due to the rather stringent requirements for predicting aerodynamic performance characteristics of flight vehicles, while flow devices involving low-speed or incompressible flow could be reasonably well designed without resorting to accurate numerical simulations. As flow devices are required to be more sophisticated and highly efficient CFD took become increasingly important in fluid engineering for incompressible and low-speed flow. This paper reviews some of the successes made possible by advances in computational technologies during the same period, and discusses some of the current challenges faced in computing incompressible flows.

Description: Transport of solid-liquid slurries in pipeline transport over short and medium distances is very important in many industries, including mining related processes. The particle image velocimetry technique was successfully utilized to investigate the velocities and kinetic energy fluctuations of slurry particles at the tongue region of an optically-clear centrifugal pump. The experiments were conducted using 500 micron glass beads at volumetric Concentrations of 2.5% and 5% and at pump speeds of 725 rpm and 1000 rpm. The fluctuation kinetic energy increased approximately 200% to 500% as the pump speed was increased from 725 rpm to IO00 rpm. The directional impingement mechanism is more significant at the pressure side of the blade, tongue and the casing. This mechanism becomes more important as the speed increases. This suggests that the impeller; tongue and the casing of the slurry pump can wear out quickly, especially with an increase in speed. In this paper the emphasis is on the tongue region. The random impingement mechanism caused by the fluctuation kinetic energy of the solids can play an important role on the erosion of the tongue area.

Description: The current cache of S/X-band geodetic/astrometric VLBI data accumulated since 1979 is approx. 5 million observations and is increasing by approx. 300,000 observations per year. The long time interval and access to all such VLBI data for re-analysis have contributed to their usefulness for the terrestrial and celestial reference frames, Earth orientation parameters, tidal and non-tidal loading, and troposphere. While data access and integrity have been maintained through the Mark III data base system as storage devices and media have evolved, past transitions have been major projects. A new format and retention concept to ensure eternal archiving and access should make use of self-documentation, generalized media, network connectivity and multiple redundancy. Similarly permanent organizations or sequences of organizations are also necessary.

Description: Solar proton events (SPEs) are known to have caused changes in constituents in the Earth's polar neutral middle atmosphere. The past four years, 2000-2003, have been replete with SPEs and huge fluxes of high energy protons occurred in July and November 2000, September and November 2001, and October 2003. The highly energetic protons produce ionizations, excitations, dissociations, and dissociative ionizations of the background constituents, which lead to the production of HOx (H, OH, HO2) and NOy (N, NO, NO2, NO3, N2O5, HNO3, HO2NO2, ClONO2, BrONO2). The HOx increases lead to short-lived ozone decreases in the polar mesosphere and upper stratosphere due to the short lifetimes of the HOx constituents. Large mesospheric ozone depletions (〉70%) due to the HOx enhancements were observed and modeled as a result of the very large July 2000 SPE. The NOy increases lead to long-lived stratospheric ozone changes because of the long lifetime of the NOy family in this region. Polar total ozone depletions 〉1% were simulated in both hemispheres for extended periods of time (several months) as a result of the NOy enhancements due to the very large SPEs.

Description: I have done detailed petrologic study of Ibitira, nominally classified as a basaltic eucrite. The Fe/Mn ratio of Ibitira pyroxenes with 〈10 mole % wollastonite component is 36.4 0.4, and is well-resolved from those of five basaltic eucrites studied for comparison; 31.2-32.2. Data for the latter completely overlap. Ibitira pyroxenes have lower Fe/Mg than the basaltic eucrite pyroxenes. Thus, the higher Fe/Mn ratio does not reflect a simple difference in oxidation state. Ibitira also has an oxygen isotopic composition, alkali element contents and a Ti/Hf ratio that distinguish it from basaltic eucrites. These differences support derivation from a distinct parent asteroid. Ibitira is the first recognized representative of the fifth known asteroidal basaltic crust.

Description: This Symposium is intended to bring together the often distinct cultures of the Stability and Control (S&C) community and the Computational Fluid Dynamics (CFD) community. The COMSAC program is itself a new effort by NASA Langley to accelerate the application of high end CFD methodologies to the demanding job of predicting stability and control characteristics of aircraft. This talk is intended to set the stage for needing a program like COMSAC. It is not intended to give details of the program itself. The topics include: 1) S&C Challenges; 2) Aero prediction methodology; 3) CFD applications; 4) NASA COMSAC planning; 5) Objectives of symposium; and 6) Closing remarks.

Description: The Shock Compression Laboratory in the Department of Earth and Planetary Sciences at Harvard is a new facility for the study of impact and collisional phenomena. The following describes the experimental capabilities of the laboratory.

Description: The motivation behind the inclusion of unsteady aerodynamics and aeroelastic effects in the computation of stability and control (S&C) derivatives will be discussed as they pertain to aeroelastic and aeroservoelastic analysis. This topic will be addressed in the context of two applications, the first being the estimation of S&C derivatives for a cable-mounted aeroservoelastic wind tunnel model tested in the NASA Langley Research Center (LaRC) Transonic Dynamics Tunnel (TDT). The second application will be the prediction of the nonlinear aeroservoelastic phenomenon known as Residual Pitch Oscillation (RPO) on the B-2 Bomber. Techniques and strategies used in these applications to compute S&C derivatives and perform flight simulations will be reviewed, and computational results will be presented.

Description: COMSAC goals include increasing the acceptance of CFD as a viable tool for S&C predictions, as well as to focus CFD development and improvement towards the needs of the S&C community. We view this as a symbiotic relationship, with increasing improvement of CFD promoting increasing acceptance by the S&C community, and increasing acceptance spurring further improvements. In this presentation we want to provide an overview for the non CFD expert of current CFD strengths and weaknesses, as well as to highlight a few emerging capabilities that we feel will lead toward increased usefulness in S&C applications.

Description: Various killing mechanisms have been suggested to contribute to the mass extinctions at the KT boundary, including severe, global deterioration of the atmosphere and hydrosphere due to SO(x) released from heavily shocked, sulfate-bearing target rocks. The devolatilization of anhydrite is predominantly inferred from thermodynamic considerations and lacks experimental confirmation. To date, the experimentally determined shock behavior of anhydrite is limited to solid-state effects employing X-ray diffraction methods. The present report employs additional methods to characterize experimentally shocked anhydrite.

Description: Osmium isotope data had shown that Ivory Coast tektites contain an extraterrestrial component, but do not allow distinction between chondritic and iron meteorite contamination. PGE abundances of Ivory Coast tektites and impactites and target rocks from the Bosumtwi crater, the source crater of the Ivory Coast tektites, were all relatively high and did not allow to resolve the presence, or identify the nature, of the meteoritic component. However, Cr isotope analyses of an Ivory Coast tektite yielded a distinct 53Cr excess of 0.30+/-0.06, which indicates that the Bosumtwi impactor was an ordinary chondrite.

Description: The Miniature Thermal Emission Spectrometer (Mini-TES) has provided remote measurements of the mineralogy and thermophysical properties of the scene surrounding the Mars Exploration Rovers. The specific scientific objectives of this investigation are to: (1) determine the mineralogy of rocks and soils; (2) determine the thermophysical properties of surface materials; and (3) determine the temperature profile, dust and water-ice opacity, and water vapor abundance in the lower atmospheric boundary layer.

Description: The second MER rover (Opportunity) landed on Meridiani Planum on January 24, 2004 inside a shallow crater. The science rational for the selection of the landing site centered on detection of the mineral hematite from martian orbit by the Mars Global Surveyor Thermal Emission Spectrometer (MGS-TES) [1,2]. Other smaller occurrences of hematite are in Aram Chaos and several isolated spots in Valles Marineris. Proposed formation pathways for martian hematite include both aqueous (e.g., low temperature precipitation of Fe oxides/oxyhydroxides in a lacustrine environment, laterite-style weathering, and precipitation from fluids having a hydrothermal origin) and dry (e.g., oxidation of magnetite rich ash) processes [e.g., 1,2,3]. The crystallographic c-face of martian hematite must be exaggerated to account for the thermal emissions spectra and it must be gray in color so as to account for the absence of the characteristic spectral signature of red hematite at visible wavelengths

Description: Though plausible explanations for the high latitude subsurface hydrogen features on Mars have been put forth, there still lacks a consensus on the nature of the low-latitude hydrogen features found in Arabia Terra and Daedalia Planum. While equivalent water mass fractions in these regions are low enough to potentially be explained by the presence of hydrated minerals, it still remains possible that such features are the remnants of ice deposits left from a previous period of high obliquity and which is now thermally unstable and subliming. In order to explore the thermal stability of putative ice deposits at low latitudes, we use the Geophysical Fluid Dynamics Laboratory (GFDL) Mars GCM with a newly integrated subsurface scheme to trace the deposition and sublimation rates of water ice, adsorbate and vapor across the planet at varying obliquities. In addition, these results help resolve the question of the dominant means of ice emplacement in the near surface, whether such ice is the result of buried surface deposits, or in situ emplaced ice due to vapor diffusion.

Description: The Flynn Creek impact structure is located in Tennessee, USA (36 deg.17 min.N, 85 deg.40 min.W). The structure was first mapped as a crypto-volcanic by Wilson and Born in 1936 [1]. Although they did not properly identify the stratigraphy within the crater or the causal mechanism, they did correctly define the horizontal extent of the crater. More detailed surface and subsurface research by Roddy (1979) accurately described the crater as being an impact structure with a diameter of 3.8 km. It formed around 360 Ma, which corresponds to the interval between the deposition of the Nashville Group and the Chattanooga Shale. Although there is limited rock outcrop in the area, there are exposed surface faults, folds, and large outcrops of impact breccia within the crater.

Description: Angrites constitute a small, but important group of basaltic achondrites showing unusual mineralogy and old crystallization ages. The currently known angrites are divided into two subgroups. Angra dos Reis (ADOR) and LEW86010 show slow cooling histories ("slowly-cooled" angrites) and differ from the later found angrites (LEW87051, Asuka 881371, Sahara 99555, D Orbigny, NWA1670, NWA1298). This second group has textures that suggest rapid cooling histories ("quenched" angrites). The petrogenesis of angrites has been controversial, partly due to the small number of available samples. In this abstract, we suggest a possible parent melt composition for the quenched angrites and its relationship to the partial melts of carbonaceous chondrites.

Description: Ureilites are ultramafic achondrites composed primarily of olivine and pyroxene with intergranular fine-grained metal, sulfides, and silicates. Ureilites contain significant amounts of carbon (up to about 6.5 wt%) as graphite, lonsdaleite, and/or diamond. It has been shown that carbon-silicate redox (i.e. "smelting") reactions are responsible for the negative FeO-MnO (or positive Fe/Mn-Fe/Mg with constant Mn/Mg) trend seen in the mineral and bulk compositions of ureilites and for the positive correlation between modal percent pigeonite and mg#. Carbon redox reactions are strongly exothermic and pressure dependent; so ureilites with the largest mg# are the most reduced, experienced the highest temperatures, and formed at the lowest pressures, i.e. near the surface of the ureilite parent body. Ureilites with the largest mg# have the smallest the delta(sup 18)O and the largest Delta(sup 17)O. To explain this, Singletary and Grove proposed that heterogeneous accretion took place on the ureilite parent body, which lead to a radial distribution of the oxygen isotopes. To further investigate possible relationships, we performed carbon isotope and electron probe measurements on a suite of 27 ureilites in order to see the type of correlation that exists between mg#, oxygen isotopes, and carbon.

Description: When a meteorite enters into the earth's atmosphere, it is immediately subjected to various chemical reactions between its minerals and the elements in our thick atmosphere. These reactions alter the mineralogy of the meteorites. Some examples of the altered minerals in some meteorites are summarized. If meteorites can be linked to their corresponding parent body asteroid in space, we can improve our understanding of the solar system. Spectroscopy has been used in determining the links between asteroids and meteorites. However, meteorite weathering in the Earth s atmosphere also has an effect on meteorite spectra, and altered spectra can easily confuse asteroid connections. The goal of this project is to study the mineral and spectral changes due to Earth's environment and to design a computer simulation to correct meteorite spectra for the effects of terrestrial weathering.

Description: Despite evidence from ground-based data that flow over mountains is a dominant source of gravity waves (GWs) for the Northern Hemisphere winter middle atmosphere, GW-related signals in global limb radiances from the Microwave Limb Sounder (MLS) on the Upper Atmosphere Research Satellite (UARS) have shown little direct evidence of mountain waves.

Description: Absorbing the electromagnetic radiation in several regions of the solar spectrum, CO2 plays an important role in the Earth radiation budget since it produces the greenhouse effect. Many natural processes in the Earth s system add and remove carbon dioxide. Overall, measurements of atmospheric carbon dioxide at different sites around the world show an increased carbon dioxide concentration in the atmosphere. At Mauna Loa Observatory (Hawaii) the measured carbon dioxide increased from 315 to 365 ppm, in the period 1958 2000 [Keeling et al., 2001]. While at the large scale, the relationship between CO2 increase and global warming is established [IPCC, 1996], at the local scale, many studies are still needed to understand regional and local sources of carbon dioxide, such as volcanoes. The volcanic areas are particularly rich in carbon dioxide; this is due to magma degassing in the summit craters region of active volcanoes, and to the presence of fractures and active faults [Giammanco et al., 1998]. Several studies estimate a global flux of volcanic CO2 (34+/-24)10(exp 6) tons/day from effusive volcanic emissions, such as the tropospheric volcanic plume (Table 1) [McClelland et al., 1989]. Plumes are a turbulent mixture of gases, solid particles and liquid droplets, emitted continuously at high temperature from summit craters, fumarolic fields or during eruptive episodes. Inside the plume, water vapour represents 70 90% of the volcanic gases. The main gaseous components are CO2, SO2, HCl, H2, H2S, HF, CO, N2 and CH4. Other plume components are volcanic ash, aqueous and acid droplets and solid sulphur-derived particles [Sparks et al., 1997]. Volcanic gases and aerosols are evidences of volcanic activity [Spinetti et al., 2003] and they have important climatic and environmental effects [Fiocco et al., 1994]. For example, Etna volcano is one of the world s major volcanic gas sources [Allard et al., 1991]. New studies on volcanic gaseous emissions have pointed out that a variation of the gas ratio CO2/SO2 is related to eruptive episodes [Caltabiano et al., 1994]. However, measurements and monitoring of volcanic carbon dioxide are difficult and often hazardous, due to the high background presence of atmospheric CO2 and the inaccessibility of volcanic sites. Hyperspectral remote sensing is a suitable technique to overcome the difficulties of ground measurement. It permits a rapid, comprehensive view of volcanic plumes and their evolution over time, detection of all gases with absorption molecular lines within the sensor s multispectral range and, in general, measurement of all the volatile components evolving from craters. The molecular and particle plume components scatter and absorb incident solar radiation. The integral of the radiation difference composes the signal measured by the remote spectrometer. The inversion technique consists of retrieving the plume component concentrations, hence decomposing the signal into the different contributions. The accuracy of remote sensing techniques depends primarily on the sensor capability and sensitivity.

Description: The high resolution UV capabilities (lamda/delta lambda = 10(sup 5)) of the Hubble Space Telscope (HST) equipped with the Space Telescope Imaging Spectrograph (STIS) reflects a need for high resolution laboratory UV spectral data base for comparison with observation.

Description: The new MLS data are consistent with convective input of H(sub 2)O into the bottom of the TTL followed by slow ascent with a maximum relative amplitude in the seasonal cycle occurring near the tropopause nearly in phase with the tropopause temperature seasonal cycle.

Description: The Hydrosphere State (HYDROS) Mission has been selected for the National Aeronautics and Space Administration (NASA) Earth System Science Pathfinder (ESSP) program. The objectives of HYDROS are to provide frequent, global measurements of surface soil moisture and surface freeze/thaw state. In order to adequately measure these geophysical parameters, a system capable of simultaneously measuring L-Band radiometer brightness temperatures at 40 km resolution and L-Band radar backscatter at 3 km resolution over a very wide swath is required. In addition, these science requirements must be satisfied under the stringent cost-cap imposed on all ESSP missions. As a solution to this challenging set of requirements, a relatively large, six meter, conically-scanning reflector antenna architecture was selected for the mission design. The HYDROS instrument will fly on a General Dynamics SA-200HP spacecraft bus. Although large deployable mesh antennas have been used in communication applications, this will mark the first time such technology is applied in a rotating configuration for high-resolution remote sensing.

Description: Experiments have shown that moderate turbulence levels can nearly double turbine blade stagnation region heat transfer. Data have also shown that heat transfer is strongly affected by the scale of turbulence as well as its level. In addition to the stagnation region, turbulence is often seen to increase pressure surface heat transfer. This is especially evident at low to moderate Reynolds numbers. Vane and rotor stagnation region, and vane pressure surface heat transfer augmentation is often seen in a pre-transition environment. Accurate predictions of transition and relaminarization are critical to accurately predicting blade surface heat transfer. An approach is described which incorporates the effects of both turbulence level and scale into a CFD analysis. The model is derived from experimental data for cylindrical and elliptical leadng edges. Results using this model are compared to experimental data for both vane and rotor geometries. The comparisons are made to illustrate that using a model which includes the effects of turbulence length scale improves agreement with data, and to illustrate where improvements in the modeling are needed.

Description: Studies of several samples of the large Caddo County IAB iron meteorite reveal andesitic material, enriched in Si, Na, Al and Ca, which is essentially unique among meteorites. This material is believed to have formed from a chondritic source by partial melting and to have further segregated by grain coarsening. Such an origin implies extended metamorphism of the IAB parent body. New Ar-39-Ar-40 ages for silicate from three different Caddo samples are consistent with a common age of 4.50-4.51 Gyr ago. Less well defined Ar-Ar degassing ages for inclusions from two other IABs, EET8333 and Udei Station, are approx.4.32 Gyr, whereas the age for Campo del Cielo varies considerably over approx.3.23-4.56 Gyr. New I-129-Xe-129 ages for Caddo County and EET8333 are 4561.9+/-0.1 Myr and 4560- 4563 Myr, respectively, relative to an age of 4566 Myr for Shallowater. Considering all reported Ar-Ar ages for IABs and related winonaites, the range is approx.4.32-4.53 Gyr, but several IABs give similar Ar ages of 4.50-4.52 Gyr. We interpret these older ages to represent cooling after the time of last significant metamorphism on the parent body, and the younger ages to represent later Ar-40 diffusion loss. These older Ar-Ar ages are similar to Sm-Nd and Rb-Sr isochron ages reported in the literature for Caddo County. Considering the possibility that IAB parent body formation was followed by impact disruption, reassembly, and metamorphism (e.g., Benedix et al. 2000), the time of the post-assembly metamorphism may have been as late as approx.4.53 Gyr ago. However, precise I-Xe ages reported for some IABs define a range of ages of approx.4560 to approx.4576 Myr. The older I-Xe ages exceed the oldest precise radiometric ages of meteorites, appear unrealistic, and s,uggest a bias in the calibration of all I-Xe ages. But even with such a bias, the I-Xe ages of IABs cannot easily be reconciled with the much younger Ar-Ar and Sm-Nd ages and with cooling rates deduced from Ni concentration profiles in IAB metal (Herpfer et al., 1994). An explanation for the difference in radiometric ages of IABs may reside in combinations of the following: a) I-Xe ages have very high closure temperatures and were not reset during metamorphism; b) a bias exists in the K-40 decay constants; c) the reported Sm-Nd and Rb-Sr ages for Caddo are in error by amounts equal to or exceeding their reported 2-sigma uncertainties; and 4) the IAB parent body may have experienced a mild metamorphism approx.30 Myr after the initial heating that produced differentiation of Caddo silicate and mixing of silicate and metal.

Description: Computational Fluid Dynamics (CFD) is becoming an important component of injector design in the rocket industry. Injector designers who use CFD in the design process need to understand the accuracy level of the particular code being used for certain aspects of the design. This paper represents a recent effort to demonstrate the accuracy of two CFD codes to calculate the wall heat flux for a single element injector. The FDNS and Loci- CHEM codes currently in use at NASA Marshall Space Flight Center were evaluated. The test case was a single shear coaxial element with gaseous oxygen and hydrogen in a chamber instrumented with coaxial heat flux gauges on the chamber wall down the axis. The data were taken at Penn State University. Measured wall temperatures were used as boundary conditions for the CFD calculations. Calculations were compared to experimental data at chamber pressures of 300,450 600, and 750 psia. The accuracy level of both codes was assessed.

Description: Ice I exhibits a complex rheology at temperature and pressure conditions appropriate for the interiors of the ice I shells of Europa, Ganymede, and Callisto. We use numerical methods and existing parameterizations of the critical Rayleigh number to determine the conditions required to trigger convection in an ice I shell with the stress-, temperature- and grain size- dependent rheology measured in laboratory experiments by Goldsby and Kohlstedt [2001]. The critical Rayleigh number depends on the ice grain size and the amplitude and wavelength of temperature perturbation issued to an initially conductive ice I shell. If the shells have an assumed uniform grain size less than 0.4 mm, deformation during initial plume growth is accommodated by Newtonian volume diffusion. If the ice grain size is between 0.4 mm and 3 cm, deformation during plume growth is accommodated by weakly non-Newtonian grain boundary sliding, where the critical ice shell thickness for convection depends on the amplitude of temperature perturbation to the _0.5 power. If the ice grain size exceeds 2 cm, convection can not occur in the ice I shells of the Galilean satellites regardless of the amplitude or wavelength of temperature perturbation. If the grain size in a convecting ice I shell evolves to effective values greater than 2 cm, convection will cease. If the ice shell has a grain size large enough to permit flow by dislocation creep, the ice is too stiff to permit convection, even in the thickest possible ice I shell. Consideration of the composite rheology implies that estimates of the grain size in the satellites and knowledge of their initial thermal states are required when judging the convective instability of their ice I shells.

Description: Solid-oxide fuel cells (SOFCs) show great potential as a power source for future space exploration missions. Because SOFCs operate at temperatures significantly higher than other types of fuel cells, they can reach overall efficiencies of up to 60% and are able to utilize fossil fuels. The SOFC team at GRC is leading NASA's effort to develop a solid oxide fuel cell with a power density high enough to be used for aeronautics and space applications, which is approximately ten times higher than ground transport targets. layers must be able to operate as a single unit at temperatures upwards of 900'C for at least 40,000 hours with less than ten percent degradation. One key challenge to meeting this goal arises from the thermal expansion mismatch between different layers. The amount a material expands upon heating is expressed by its coefficient of thermal expansion (CTE). If the CTEs of adjacent layers are substantially different, thermal stresses will arise during the cell's fabrication and operation. These stresses, accompanied by thermal cycling, can fracture and destroy the cell. While this is not an issue at the electrolyte-cathode interface, it is a major concern at the electrolyte-anode interface, especially in high power anode-supported systems. electrolyte are nearly identical. Conventionally, this has been accomplished by varying the composition of the anode to match the CTE of the yittria-stabilized zirconia (YSZ) electrolyte (approx.10.8x10(exp -6/degC). A Ni/YSZ composite is typically used as a base material for the anode due to its excellent electrochemical properties, but its CTE is about 13.4x10(exp -6/degC). One potential way to lower the CTE of this anode is to add a small percentage of polycrystalline Al2TiO5, with a CTE of 0.68x10(exp -6/degC, to the Ni/YSZ base. However, Al2TiO5 is thermally unstable and loses its effectiveness as it decomposes to Al2O3 and TiO2 between 750 C and 1280 C. be used as additives to increase the thermal stability of Al2TiO5 in SOFC operating conditions without adversely affecting the electrochemical properties of the SOFC anode. Three candidate materials were chosen through an extensive literature review: MgO, Fe2O3, and ZrTiO4. Although all three have been shown to prevent Al2TiO5 decomposition under various conditions, their effectiveness in the temperature range and atmosphere of the SOFC has not yet been evaluated. Several batches of Al2TiO5 with varying amounts of additives were prepared, exposed to reducing and oxidizing atmospheres at elevated temperatures, and the resulting decomposition of Al2TiO5 was measured. The most promising additives were further evaluated with the goal of ultimately preparing low CTE anodes that are chemically compatible to current systems. Adding minor constituents to stabilize Al2TiO5 could ultimately preserve its low CTE for the life of the fuel cell and improve the cell's long-term performance without a drop in anode conductivity. Further, these low CTE filler additions could allow the use of new sulfur tolerant anode materials, improving the viability of SOFCs for future aeronautics and space applications. Every SOFC consists of a cathode and an anode separated by an electrolyte, These three One way to avoid this problem is to design the cell such that the CTEs of the anode and The objective of this summer research project was to evaluate several materials that could

Description: The Gravity Recovery and Climate Experiment (GRACE), launched on March 17, 2002, represents the state-of-the-art in geodetic observations of the static and time varying components of the Earth's geopotential field. The fundamental measurement used to observe gravity is the inter-satellite range and range rate between two coplanar, low altitude satellites obtained from a K-band ranging (KBR) system. In addition to the K-band ranging system, each satellite possess a Super-STAR Accelerometer, a GPS receiver/antenna package, Star Cameras and a Laser Retro Reflector (LRR) to complete the compliment of science instruments. The GRACE project has now released two years of Level 1B data derived from the science instruments and sensors. An integral component of our time variable gravity research is the reduction, calibration and analyses of these Level 1B data. In particular we have analyzed several months of K-band ranging (KBR1B), accelerometry (ACC1B) and GPS navigation (GNAV1B) data. Accelerometer calibration and KBR data reduction methodology and results will be presented. We discuss the impact of these analyses on the recovery of time variable gravity.

Description: The salar de Uyuni is a massive dry salt lake that lies at the lowest point of an internal/drainage basin in the Bolivian Altiplano. Its topography is remarkable for its extraordinary flatness over almost a full degree of latitude and longitude. We surveyed a 54 x 45 km region of the salar with kinematic GPS in September, 2002 and found a topographic range of only 80 cm over the entire surveyed area. Furthermore, the survey revealed distinct surface features with several dominant wavelengths and orientations. Some of these appear to be aligned with orographic features that intersect the salar, leading us to conjecture that they are the surface expression of high-density mountains that have been buried by low-density basin sediments. Over the oceans, a similar correspondence between basin bathymetry and surface topography is exploited to map the seafloor using sea-surface satellite altimetry measurements, with the sea surface following geoid undulations due to the underwater mass distribution. On the salar, annual flooding creates a shallow lake whose surface also lies on a equipotential surface shaped by the distribution of underlying mass. The link to the actual salar surface is via the dissolution and redeposition of salt by the lake waters, which appears to push the system to an equilibrium of constant water depth and the coincidence of the shapes of the lake surface and bottom. To test our hypothesis about the origin of the surface features on the salar, we compare our GPS survey elevations with the equipotential surface generated from local gravity measurements in conjunction with gravity and potential values from the EGM96 global geopotential model. 50% of the variance of the GPS elevations can be explained by equipotential surface undulations from the EGM96 model alone, and an additional 40% is explained by the shorter-wavelength equipotential surface derived from local gravity. We examine the unexplained 10% of elevation variance from the standpoint of errors in the equipotential surface calculation and possible unmodelled surface processes.

Description: NASA's Geoscience Laser Altimeter System (GLAS) on the Ice, Cloud and land Elevation Satellite (ICESat) is Earth's first polar-orbiting satellite laser altimeter. Its primary purpose is to measure ice sheet elevation change, with a scientific requirement of detecting changes as low as 1.5 centimeters per year. This goal requires precise calibration and validation of the instrument. One approach for validating the GLAS range measurement involves comparison with a land reference target. Dry salt lakes are ideal for this purpose since they are large, stable, easily surveyed using kinematic GPS, and have albedos similar to that of ice. We selected the largest dry salt lake in the world, the salar de Uyuni in Bolivia, as a reference target for GLAS. In September 2002 we surveyed a 54 x 45 kilometer area of the salar with car-mounted kinematic GPS, using the data to construct a DEM of the surface. To date there have been six ICESat passes across our survey area, along two intersecting tracks. The first two passes occurred during "Laser 2A" operation (September to November 2003), the second two during "Laser 2B" operation (February to March 2004), and the final two during "Laser 2C" (June to July 2004). We compare GLAS elevations from all six passes with the salar de Uyuni DEM, showing how differing conditions between passes affected the performance of the GLAS instrument and providing an estimate of the relative and absolute accuracy of the range measurement. To check our results, we also compare the GLAS elevations along the same tracks with the mean sea-surface derived from TOPEX over the ocean on either side of South America.

Description: Salt flats are aptly named: they are composed largely of salt, and are maintained as nearly equipotential surfaces via frequent flooding. The salar de Uyuni, on the Altiplano in southwestern Bolivia, is the largest salt flat on Earth, with an area of 9,800 sq km. Except for a few bedrock islands, it has less than 40 cm of relief. The upper-most salt unit averages 5 m thick and contains 50 cu km of nearly pure halite. It includes most of the salt that was in solution in paleolake Minchin, which attained a maximum area of 60,000 sq km and a maximum depth of 150 m, roughly 15 kyr ago. Despite approx. 10 m of differential isostatic rebound since deposition, the salar surface has been actively maintained as an extraordinarily flat and smooth surface by annual flooding during the rainy season. We have used the strong optical absorption properties of water in the visible band to map spatial variations in water depth during a time when the salar was flooded. As water depth increases, the initially pure white surface appears both darker and bluer. We utilized MISR images taken during the interval from April to November 2001. The red and infra-red bands (672 and 867 nm wavelength) were most useful since the water depth is small and the absorption at those wavelengths is quite strong. Nadir pointed MISR images have 275 m spatial resolution. To aid in our evaluation of water depth variations over the saiar surface, we utilized two sources of direct topographic measurements: several ICESAT altimetry tracks cross the area, and a 40x50 km GPS grid was surveyed to calibrate ICESAT. A difficulty in using these data types is that both give salt surface elevations relative to the ellipsoid, whereas the water surface will, in the absence of wind or tidal disturbances, follow an equipotential surface. Geoid height is not known to the required accuracy of a few cm in the central Andes. As a result, before comparing optical absorption from MISR to salt surface topography from GPS or ICESAT, we removed the longest wavelengths from both.

Description: We have analyzed GRACE Level 1-B data in 2003 and assessed a new approach for extracting time variable gravity that isolates the gravity signal in both time and space. The Level-1B satellite-to-satellite range rate (KBRR) data and accelerometry are processed in daily arcs using the precise orbit products produced by the GRACE team from GPS to calibrate both the accelerometer and KBRR data. We then adjusted select components of the intersatellite baseline vector for each data segment isolated to the region of interest. Herein, we solved for mass anomalies in 45 deg x 45 deg blocks over the Amazon and the nearby Atlantic Ocean and estimate mass flux in units of cm of water over each block. We show with this approach that we can recover mass anomalies on a submonthly basis with 10 to 15 day temporal resolution. We discuss the important issues related to this solution, including the size of the mascon blocks, the weight given to the temporal and spatial constraint used to stabalize the solutions, as well as the optimal correlation in time and distance. We compare the the mascon results with solutions obtained from the more standard approach using spherical harmonics and with independent hydrology models and lake data. This technique demonstrates that sub-monthly medium wavelength mass flux phenomena are well sensed by the hyper-precise line of sight velocity data produced from GRACE.

Description: The purpose of this project was to provide the community access to magnetosheath data near Earth. We provided 27 years of IMP 8 magnetosheath proton velocities, densities, and temperatures with our best (usually 1-min.) time resolution. IMP 8 crosses the magnetosheath twice each 125 day orbit, and we provided magnetosheath data for the roughly 27 years of data for which magnetometer data are also available (which are needed to reliably pick boundaries). We provided this 27 years of IMP 8 magnetosheath data to the NSSDC; this data is now integrated with the IMP 8 solar wind data with flags indicating whether each data point is in the solar wind, magnetosheath, or at the boundary between the two regions. The plasma speed, density, and temperature are provided for each magnetosheath point. These data are also available on the MIT web site ftp://space .mit.edu/pub/plasma/imp/www/imp.html. We provide ASCII time-ordered rows of data giving the observation time, the spacecraft position in GSE, the velocity is GSE, the density and temperature for protons. We also have analyzed and archived on our web site the Wind magnetosheath plasma parameters. These consist of ascii files of the proton and alpha densities, speeds, and thermal speeds. These data are available at ftp://space.mit.edu/pub/plasma/wind/sheath These are the two products promised in the work statement and they have been completed in full.

Description: We present latest result in understanding the potential relationship between tectonic stress, electro-chemical and thermodynamic processes in the Earths crust and atmosphere with an increase in IR flux as a potential signature of electromagnetic (EM) phenomena that are related to earthquake activity, either pre-, co- or post seismic. Thermal infra-red (TIR) surveys performed by the polar orbiting (NOAA/AVHRR MODIS) and geosynchronous weather satellites (GOES, METEOSAT) gave an indication of the appearance (from days to weeks before the event) of "anomalous" space-time TIR transients that are associated with the location (epicenter and local tectonic structures) and time of a number of major earthquakes with M〉5 and focal depths less than 50km. We analyzed broad category of associated pre-earthquake events, which provided evidence for changes in surface temperature, surface latent heat flux, chlorophyll concentrations, soil moisture, brightness temperature, emissivity of surface, water vapour in the atmosphere prior to the earthquakes occurred in Algeria, India, Iran, Italy, Mexico and Japan. The cause of such anomalies has been mainly related to the change of near-surface thermal properties due to complex lithosphere-hydrosphere-atmospheric interactions. As final results we present examples from the most recent (2000-2004) worldwide strong earthquakes and the techniques used to capture the tracks of EM emission mid-IR anomalies and a methodology for practical future use of such phenomena in the early warning systems.

Description: We investigate the time-variable gravity changes in Europe retrieved from the initial GRACE monthly solutions spanning a 18 month duration from April 2002 to October 2003. Gravity anomaly maps are retrieved in Central Europe from the monthly satellite solutions we compare the fields according to various truncation levels (typically between degree 10 and 20) of the initial fields (expressed in spherical harmonics to degree 120). For these different degrees, an empirical orthogonal function (EOF) decomposition of the time-variable gravity field leads us to its main spatial and temporal characteristics. We show that the dominant signal is found to be annual with an amplitude and a phase both in agreement with predictions in Europe modeled using snow and soil-moisture variations from recent hydrology models. We compare these GRACE gravity field changes to surface gravity observations from 6 superconducting gravimeters of the GGP (Global Geodynamics Project) European sub-network, with a special attention to loading corrections. Initial results suggest that all 3 data sets (GRACE, hydrology and GGP) are responding to annual changes in near-surface water in Europe of a few microGal (at length scales of approx.1000 km) that show a high value in winter and a summer minimum. We also point out that the GRACE gravity field evolution seems to indicate that there is a trend in gravity between summer 2002 and summer 2003 which can be related to the 2003 heatwave in Europe and its hydrological consequences (drought). Despite the limited time span of our analysis and the uncertainties in retrieving a regional solution from the network of gravimeters, the calibration and validation aspects of the GRACE data processing based on the annual hydrology cycle in Europe are in progress.