ALBERT

Description: An explicit time filter is applied to the Navier-Stokes equation prior to a space filter. The time filter is supposed to be smooth, and an exact expansion depending on the time derivatives of the velocity is derived for the associated stress tensor. On the contrary, the effect of the space filter is treated as usual and an eddy viscosity model is introduced in the LES equation. The total stress is thus represented using a new class of mixed models combining time and space derivatives of the LES field.

Description: As pointed out by Rodi standard integral solutions for jets and plumes developed for discharge into infinite, quiescent ambient are difficult to extend to complex situations, particularly in the presence of boundaries such as the sea floor or ocean surface. In such cases the assumption of similarity breaks down and it is impossible to find a suitable entrainment coefficient. The models are also incapable of describing any but the most slowly varying unsteady motions. There is therefore a need for full time-dependent modeling of the flow field for which there are three main approaches: (1) Reynolds averaged numerical simulation (RANS), (2) large eddy simulation (LES), and (3) direct numerical simulation (DNS). Rodi applied RANS modeling to both jets and plumes with considerable success, the test being a match with experimental data for time-averaged velocity and temperature profiles as well as turbulent kinetic energy and rms axial turbulent velocity fluctuations. This model still relies on empirical constants, some eleven in the case of the buoyant jet, and so would not be applicable to a partly laminar plume, may have limited use in the presence of boundaries, and would also be unsuitable if one is after details of the unsteady component of the flow (the turbulent eddies). At the other end of the scale DNS modeling includes all motions down to the viscous scales. Boersma et al. have built such a model for the non-buoyant case which also compares well with measured data for mean and turbulent velocity components. The model demonstrates its versatility by application to a laminar flow case. As its name implies, DNS directly models the Navier-Stokes equations without recourse to subgrid modeling so for flows with a broad spectrum of motions (high Re) the cost can be prohibitive - the number of required grid points scaling with Re(exp 9/4) and the number of time steps with Re(exp 3/4). The middle road is provided by LES whereby the Navier-Stokes equations are formally filtered with the filter chosen to only exclude the smallest turbulent motions. If successful, LES should provide much of the detail available to DNS but at more bearable cost. Fatica et al. in comparing LES with DNS for a low Reynolds number jet showed that the LES could simulate the temporally evolving behavior including growth of the jet thickness. It is the intention of this report to explore the application of an LES model to jets and plumes. As always, before tackling complex situations, the model must be tested for the simplest of cases and so we address only two, a non-buoyant axisymmetric jet issuing steadily from an orifice into a semi-infinite stationary environment and a buoyant jet in the same environment. The work is a continuation of Basu and Mansour.

Description: This paper compares the filtering used in Coherent Vortex Simulation (CVS) decomposition with an orthogonal wavelet basis, with the Proper Orthogonal Decomposition (POD) or Fourier filtering. Both methods are applied to a field of Direct Numerical Simulation (DNS) data of 3D forced homogeneous isotropic turbulence at microscale Reynolds number R(sub lambda) = 168. We show that, with only 3%N retained modes, CVS filtering separates the coherent vortex tubes from the incoherent background flow. The latter is structureless, has an equipartition energy spectrum, and has a Gaussian velocity probability distribution function (PDF) and an exponential vorticity PDF. On the other hand, the Fourier basis does not extract the coherent vortex tubes cleanly and leaves organized structures in the residual high wavenumber modes whose PDFs are stretched exponentials for both the velocity and the vorticity.

Description: The objective of this study is to investigate the use of the second generation bi-orthogonal wavelet transform for the field decomposition in the Coherent Vortex Simulation of turbulent flows. The performances of the bi-orthogonal second generation wavelet transform and the orthogonal wavelet transform using Daubechies wavelets with the same number of vanishing moments are compared in a priori tests using a spectral direct numerical simulation (DNS) database of isotropic turbulence fields: 256(exp 3) and 512(exp 3) DNS of forced homogeneous turbulence (Re(sub lambda) = 168) and 256(exp 3) and 512(exp 3) DNS of decaying homogeneous turbulence (Re(sub lambda) = 55). It is found that bi-orthogonal second generation wavelets can be used for coherent vortex extraction. The results of a priori tests indicate that second generation wavelets have better compression and the residual field is closer to Gaussian. However, it was found that the use of second generation wavelets results in an integral length scale for the incoherent part that is larger than that derived from orthogonal wavelets. A way of dealing with this difficulty is suggested.

Description: The main goal of this presentation is to give some of the objectives of the testing program. This includes: develop jet noise data base for separate flow nozzles with bypass ratio's 5 to 14; evaluate effect of pylon on noise; develop low performance impact noise suppression concepts; and evaluate potential for active control of jet noise.

Description: Coherent Vortex Simulation (CVS) filtering has been applied to Direct Numerical Simulation (DNS) data of forced and unforced time-developing turbulent mixing layers. CVS filtering splits the turbulent flow into two orthogonal parts, one corresponding to coherent vortices and the other to incoherent background flow. We have shown that the coherent vortices can be represented by few wavelet modes and that these modes are sufficient to reproduce the vorticity probability distribution function (PDF) and the energy spectrum over the entire inertial range. The remaining incoherent background flow is homogeneous, has small amplitude, and is uncorrelated. These results are compared with those obtained for the same compression rate using large eddy simulation (LES) filtering. In contrast to the incoherent background flow of CVS filtering, the LES subgrid scales have a much larger amplitude and are correlated, which makes their statistical modeling more difficult.

Description: This presentation will discuss Ir for online diagnostics with acoustics; Ptot, Ttot, Pstat rake surveys for mean flow measurements; focused Schleiren for density and some turbulence structure; Laser sheet visualization for near-nozzle diagnostics; and Two-point hotwire measurements for turbulence models.

Description: The Lagrangian-averaged Navier-Stokes (LANS) equations are numerically evaluated as a turbulence closure. They are derived from a novel Lagrangian averaging procedure on the space of all volume-preserving maps and can be viewed as a numerical algorithm which removes the energy content from the small scales (smaller than some a priori fixed spatial scale alpha) using a dispersive rather than dissipative mechanism, thus maintaining the crucial features of the large scale flow. We examine the modeling capabilities of the LANS equations for decaying homogeneous turbulence, ascertain their ability to track the energy spectrum of fully resolved direct numerical simulations (DNS), compare the relative energy decay rates, and compare LANS with well-accepted large eddy simulation (LES) models.

Description: We have completed a new generation of water vapor radiometers (WVR), the A- series, in order to support radio science experiments with the Cassini spacecraft. These new instruments sense three frequencies in the vicinity of the 22 GHz emission line of atmospheric water vapor within a 1 degree beamwidth from a clear aperture antenna that is co-pointed with the radio telescope down to 10 degree elevation. The radiometer electronics features almost an order of magnitude improvement in temperature stability compared with earlier WVR designs. For many radio science experiments, the error budget is likely to be dominated by path delay fluctuations due to variable atmospheric water vapor along the line-of-sight to the spacecraft. In order to demonstrate the performance of these new WVRs we are attempting to calibrate the delay fluctuations as seen by a radio interferometer operating over a 21 km baseline with a WVR near each antenna. The characteristics of these new WVRs will be described and the results of our preliminary analysis will be presented indicating an accuracy of 0.2 to 0.5 mm in tracking path delay fluctuations over time scales of 10 to 10,000 seconds.

Description: Science education is taking the teaching of science from a traditional (lecture) approach to a multidimensional sense-making approach which allows teachers to support students by providing exploratory experiences. Using projects is one way of providing students with opportunities to observe and participate in sense-making activity. We created a learning environment that fostered inquiry-based learning. Students were engaged in a variety of Inquiry activities that enabled them to work in cooperative planning teams where respect for each other was encouraged and their ability to grasp, transform and transfer information was enhanced. Summer, 1998: An air pollution workshop was conducted for high school students in the Medgar Evers College/Middle College High School Liberty Partnership Summer Program. Students learned the basics of meteorology: structure and composition of the atmosphere and the processes that cause weather. The highlight of this workshop was the building of hand-held sunphotometers, which measure the intensity of the sunlight striking the Earth. Summer, 1999: high school students conducted a research project which measured the mass and size of ambient particulates and enhanced our ability to observe through land based measurements changes in the optical depth of ambient aerosols over Brooklyn. Students used hand held Sunphotometers to collect data over a two week period and entered it into the NASA GISS database by way of the internet.

Description: Rapid climate change characterizes numerous terrestrial sediment records during and since the last glaciation. Vegetational response is best expressed in terrestrial records near ecotones, where sensitivity to climate change is greatest, and response times are as short as decades.

Description: Geological, geophysical, and geochemical data support a theory that Earth experienced several intervals of intense, global glaciation ("snowball Earth" conditions) during Precambrian time. This snowball model predicts that postglacial, greenhouse-induced warming would lead to the deposition of banded iron formations and cap carbonates. Although global glaciation would have drastically curtailed biological productivity, melting of the oceanic ice would also have induced a cyanobacterial bloom, leading to an oxygen spike in the euphotic zone and to the oxidative precipitation of iron and manganese. A Paleoproterozoic snowball Earth at 2.4 Giga-annum before present (Ga) immediately precedes the Kalahari Manganese Field in southern Africa, suggesting that this rapid and massive change in global climate was responsible for its deposition. As large quantities of O(2) are needed to precipitate this Mn, photosystem II and oxygen radical protection mechanisms must have evolved before 2.4 Ga. This geochemical event may have triggered a compensatory evolutionary branching in the Fe/Mn superoxide dismutase enzyme, providing a Paleoproterozoic calibration point for studies of molecular evolution.

Description: The Texas A&M monthly total oceanic rainfall retrieval algorithm is based on radiative transfer models and can only be modified on a physically sound basis. Within this constraint we have examined some improvements to the algorithm and it appears that it can be made significantly better. In particular, it appears that by proper use of the range of frequencies available on TMI (TRMM Microwave Imager) and AMSR that the need for the log-normal fit can be eliminated.

Description: We evaluated the performance of the Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI) at-launch algorithm for monthly oceanic rain rate using two years (January 1998 - December 1999) of TMI data. The TMI at-launch algorithm is based on Wilheit et al.'s technique for estimating monthly oceanic rainfall that relies on histograms of multichannel microwave measurements. Comparisons with oceanic monthly rain rates derived from the Defense Meteorological Satellite Program (DMSP) F-13 and F-14 Special Sensor Microwave Imager (SSM/I) data show the average rain rates over the TRMM region (between 400S and 40N) are 3.0, 2.85 and 2.89 mm/day, respectively for F-13, F-14 and TMI. Based on the latest version of TB data (version 5), both rainrate and freezing height derived from TMI are similar to those from the F-13 and F-14 SSM/I data. However, regionally the differences are statistically significant at the 95% confidence. Three hourly monthly rainrates are also computed from 3-hourly TB histograms to examine the diurnal cycle of precipitation. Over most of the oceanic TRMM area, a distinct early morning rainfall peak is found. A harmonic analysis shows that the amplitude of the 12h harmonic is significant and comparable to that of the 24h harmonic.

Description: A one-week in situ intercomparison campaign was completed on the Rice University campus for measuring HCHO using three different techniques, including a novel optical sensor based on difference frequency generation (DFG) operating at room temperature. Two chemical derivatization methods, 2,4-dinitrophenylhydrazine (DNPH) and o-(2,3,4,5,6-pentafluorobenzyl) hydroxylamine (PFBHA), were deployed during the daylight hours for three- to four-hour time-integrated samples. A real-time optical sensor based on laser absorption spectroscopy was operated simultaneously, including nighttime hours. This tunable spectroscopic source based on difference frequency mixing of two fiber-amplified diode lasers in periodically poled LiNb03 (PPLN) was operated at 3.5315 micrometers (2831.64 cm 1) to access a strong HCHO ro-vibrational transition free of interferences from other species. The results showed a bias of -1.7 and -1.2 ppbv and a gross error of 2.6 and 1.5 ppbv for DNPH and PFBHA measurements, respectively, compared with DFG measurements. These results validate the DFG sensor for time-resolved measurements of HCHO in urban areas.

Description: Despite the fact that supercritical fluids occur both in nature and in industrial situations, the fundamentals of their behavior is poorly understood because supercritical fluids combine the characteristics of both liquids and gases, and therefore their behavior is not intuitive. There are several specific reasons for the lack of understanding: First, data from (mostly optical) measurements can be very misleading because regions of high density thus observed are frequently identified with liquids. A common misconception is that if in an experiment one can optically identify "drops" and "ligaments", the observed fluid must be in a liquid state. This inference is incorrect because in fact optical measurements detect any large change (i.e. gradients) in density. Thus, the density ratio may be well below Omicron(10(exp 3)) that characterizes its liquid/gas value, but the measurement will still identify a change in the index of refraction providing that the change is sudden (steep gradients). As shown by simulations of supercritical fluids, under certain conditions the density gradients may remain large during the supercritical binary fluids mixing, thus making them optically identifiable. Therefore, there is no inconsistency between the optical observation of high density regions and the fluids being in a supercritical state. A second misconception is that because a fluid has a liquid-like density, it is appropriate to model it as a liquid. However, such fluids may have liquid-like densities while their transport properties differ from those of a liquid. Considering that the critical pressure of most fuel hydrocarbons used in Diesel and gas turbine engines is in the range of 1.5 - 3 MPa, and the fact that the maximum pressure attained in these engines is about 6 Mps, it is clear that the fuel in the combustion chamber will experience both subcritical and supercritical conditions. Studies of drop behavior over a wide range of pressures were performed in the past, however none of these studies identified the crucial differences between the subcritical and supercritical behavior. In fact, in two of these studies, it was found that the subcritical and supercritical behavior is similar as the drop diameter decreased according to the classical d(exp 2)-law over a wide range of pressures and drop diameters. The present study is devoted to the exploration of differences in fluid-behavior characteristics under subcritical and supercritical conditions in the particular case of heptane fluid drops in nitrogen; these substances were selected because of the availability of experimental observations for model validation.

Description: Submillimeter-wave cloud ice radiometry is an innovative technique for determining the amount of ice present in cirrus clouds, measuring median crystal size, and constraining crystal shape. The radiometer described in this poster is being developed to acquire data to validate radiometric retrievals of cloud ice at submillimeter wavelengths. The goal of this effort is to develop a technique to enable spaceborne characterization of cirrus, meeting key climate modeling and NASA measurement needs.

Description: A summary is presented of basic lightning characteristics/criteria applicable to current and future aerospace vehicles. The paper provides estimates on the probability of occurrence of a 200 kA peak lightning return current, should lightning strike an aerospace vehicle in various operational phases, i.e., roll-out, on-pad, launch, reenter/land, and return-to-launch site. A literature search was conducted for previous work concerning occurrence and measurement of peak lighting currents, modeling, and estimating the probabilities of launch vehicles/objects being struck by lightning. This paper presents a summary of these results.

Description: There are three kinds of electrohydrodynamics (EHD) pumping based on Coulomb force: induction pumping, ion-drag pumping, and pure conduction pumping. EHD induction pumping relies on the generation of induced charges. This charge induction in the presence of an electric field takes place due to a non-uniformity in the electrical conductivity of the fluid which can be caused by a non-uniform temperature distribution and/or an inhomogeneity of the fluid (e.g. a two-phase fluid). Therefore, induction pumping cannot be utilized in an isothermal homogeneous liquid. In order to generate Coulomb force, a space charge must be generated. There are two main mechanisms for generating a space charge in an isothermal liquid. The first one is associated with the ion injection at a metal/liquid interface and the related pumping is referred to as ion-drag pumping. Ion-drag pumping is not desirable because it can deteriorate the electrical properties of the working fluid. The second space charge generation mechanism is associated with the heterocharge layers of finite thickness in the vicinity of the electrodes. Heterocharge layers result from dissociation of the neutral electrolytic species and recombination of the generated ions. This type of pumping is referred to as pure conduction pumping. This project investigates the EHD pumping through pure conduction phenomenon. Very limited work has been conducted in this field and the majority of the published papers in this area have mistakenly assumed that the electrostriction force was responsible for the net flow generated in an isothermal liquid. The main motivation behind this study is to investigate an EHD conduction pump for a two-phase loop to be operated in the microgravity environment. The pump is installed in the liquid return passage (isothermal liquid) from the condenser section to the evaporator section. Unique high voltage and ground electrodes have been designed that generate sufficient pressure heads with very low electric power requirements making the EHD conduction pumping attractive to applications such as two-phase systems (e.g. capillary pumped loops and heat pipes). Currently, the EHD conduction pump performance is being tested on a two-phase loop under various operating conditions in the laboratory environment. The simple non-mechanical and lightweight design of the EHD pump combined with the rapid control of performance by varying the applied electric field, low power consumption, and reliability offer significant advantages over other pumping mechanisms; particularly in reduced gravity applications.

Description: In earlier investigations, Marangoni effects were observed to be the dominant mechanism of boiling transport in 2-propanol/water mixtures under reduced gravity conditions. In this investigation we have examined the mechanisms of binary mixture boiling by exploring the transport near a single bubble generated in a binary mixture between a heated surface and cold surface. The temperature field created in the liquid around the bubble produces vaporization over the portion of its interface near the heated surface and condensation over portions of its interface near the cold surface. Experiments were conducted using different mixtures of water and 2-propanol under 1g conditions and under reduced gravity conditions aboard the KC135 aircraft. Since 2-propanol is more volatile than water, there is a lower concentration of 2-propanol near the hot surface and a higher concentration of 2-propanol near the cold plate relative to the bulk quantity. This difference in interface concentration gives rise to strong Marangoni effects that move liquid toward the hot plate in the near bubble region for 2-propanol and water mixtures. In the experiments in this study, the pressure of the test system was maintained at about 5 kPa to achieve the full spectrum of boiling behavior (nucleate boiling, critical heat flux and film boiling) at low temperature and heat flux levels. Heat transfer data and visual documentation of the bubble shape were extracted from the experimental results. In the 1-g experiments at moderate to high heat flux levels, the bubble was observed to grow into a mushroom shape with a larger top portion near the cold plate due to the buoyancy effect. The shape of the bubble was somewhat affected by the cold plate subcooling and the superheat of the heated surface. At low superheat levels for the heated surface, several active nucleation sites were observed, and the vapor stems from them merged to form a larger bubble. The generation rate of vapor is moderate in this regime and the bubble shape is cylindrical in appearance. In some instances, the bubble interface appeared to oscillate. At higher applied heat flux levels, the top of the bubble became larger, apparently to provide more condensing interface area adjacent to the cold plate. Increasing the applied heat flux ultimately led to dry-out of the heated surface, with conditions just prior to dryout corresponding to the maximum heat flux (CHF). A more stable bubble was observed when the system attained the minimum heat flux (for film boiling). In this regime, most of the surface under the bottom of the bubble was dry with nucleate boiling sometimes occuring around the contact perimeter of the bubble at heated surface. Different variations (e.g. gap between two plates, molar concentration of the liquid mixture) of the experiments were examined to determine parametric effects on the boiling process and to determine the best conditions for the KC135 reduced gravity tests. Variation of the gap was found to have a minor impact on the CHF. However, reducing the gap between the hot and cold surface was observed to significantly reduce the minimum heat flux for fixed molar concentration of 2-propanol. In the reduced gravity experiments aboard the KC135 aircraft, the bubble formed in the 6.4 mm gap was generally cylindrical or barrel shaped and it increased its extent laterally as the surface superheat increased. In reduced gravity experiments, dryout of the heated surface under the bubble was observed to occur at a lower superheated temperature than for 1g conditions. Observed features of the boiling process and heat transfer data under reduced gravity will be discussed in detail. The results of the reduced gravity experiments will also be compared to those obtained in comparable 1g experiments. In tandem with the experiments we are also developing a computational model of the transport in the liquid surrounding the bubble during the boiling process. The computational model uses a level set method to model motion of the interface. It will incorporate a macroscale treatment of the transport in the liquid gap between the surfaces and a microscale treatment of transport in the regions between the bubble interface and the solid surfaces. The features of the model will be described in detail. Future research directions suggested by the results to date will also be discussed.

Description: The goals of this study are the evaluation of current fast radiative transfer models (RTMs) and line-by-line (LBL) models. The intercomparison focuses on the modeling of 11 representative sounding channels routinely used at numerical weather prediction centers: seven HIRS (High-resolution Infrared Sounder) and four AMSU (Advanced Microwave Sounding Unit) channels. Interest in this topic was evidenced by the participation of 24 scientists from 16 institutions. An ensemble of 42 diverse atmospheres was used and results compiled for 19 infrared models and 10 microwave models, including several LBL RTMs. For the first time, not only radiances, but also Jacobians (of temperature, water vapor, and ozone) were compared to various LBL models for many channels. In the infrared, LBL models typically agree to within 0.05-0.15 K (standard deviation) in terms of top-of-the-atmosphere brightness temperature (BT). Individual differences up to 0.5 K still exist, systematic in some channels, and linked to the type of atmosphere in others. The best fast models emulate LBL BTs to within 0.25 K, but no model achieves this desirable level of success for all channels. The ozone modeling is particularly challenging. In the microwave, fast models generally do quite well against the LBL model to which they were tuned. However significant differences were noted among LBL models. Extending the intercomparison to the Jacobians proved very useful in detecting subtle and more obvious modeling errors. In addition, total and single gas optical depths were calculated, which provided additional insight on the nature of differences. Recommendations for future intercomparisons are suggested.

Description: Nearly three years of Tropical Rainfall Measuring Mission Satellite (TRMM Satellite) monthly estimates of tropical surface rainfall are analyzed to document and understand the differences among the TRMM-based estimates and how these differences relate to the pre-TRMM estimates and current operational analyses. Variation among the TRMM estimates is shown to be considerably smaller than among a pre-TRMM collection of passive microwave-based products. Use of both passive and active microwave techniques in TRMM should lead to increased confidence in converged estimates. Current TRMM estimates are shown to have a range of about 20% for the tropical ocean as a whole, with variations in heavily raining ocean areas of the ITCZ and SPCZ having differences over 30%. In mid-latitude ocean areas the differences are smaller. Over land there is a distinct difference between the tropics and mid-latitude with a reversal between some of the products as to which tends to be relatively high or low. Comparisons of TRMM estimates with ocean atoll and land gauge information point to products that might have significant regional biases. The radar-based product is significantly low biased compared with atoll raingauge data, while the passive microwave product is significantly high compared to raingauge data in the deep tropics. The evolution of rainfall patterns during the recent change from intense El Nino to a long period of La Nina and then a gradual return to near neutral conditions is described using TRMM. The time history of integrated rainfall over the tropical oceans (and land) during this period differs among the passive and active microwave TRMM estimates.

Description: Observations made by the Precipitation Radar (PR) and the Microwave Imager (TMI) radiometer on board the Tropical Rainfall Measuring Mission (TRMM) satellite help us to show the significance of the 85 GHz polarization difference, PD85, measured by TMI. Rain type, convective or stratiform, deduced from the PR allows us to infer that PD85 is generally positive in stratiform rain clouds, while PD85 can be markedly negative in deep convective rain clouds. Furthermore, PD85 increases in a gross manner as stratiform rain rate increases. On the contrary, in a crude fashion PD85 decreases as convective rain rate increases. From the observations of TMI and PR, we find that PD85 is a weak indicator of rain rate. Utilizing information from existing polarimetric radar studies, we infer that negative values of PD85 are likely associated with vertically-oriented small oblate or wet hail that are found in deep convective updrafts.

Description: Major droughts and floods over the U.S. continent may be related to a far field energy source in the Asian Pacific. This is illustrated by two climate patterns associated with summertime rainfall over the U.S. and large-scale circulation on interannual timescale. The first shows an opposite variation between the drought/flood over the Midwest and that over eastern and southeastern U.S., coupled to a coherent wave pattern spanning the entire East Asia-North Pacific-North America region related to the East Asian jetstream. The second shows a continental-scale drought/flood in the central U.S., coupled to a wavetrain linking Asian/Pacific monsoon region to North America.

Description: A set of global, monthly rainfall products has been intercompared to understand the quality and utility of the estimates. The products include 25 observational (satellite-based), four model and two climatological products. The results of the intercomparison indicate a very large range (factor of two or three) of values when all products are considered. The range of values is reduced considerably when the set of observational products is limited to those considered quasi-standard. The model products do significantly poorer in the tropics, but are competitive with satellite-based fields in mid-latitudes over land. Over ocean, products are compared to frequency of precipitation from ship observations. The evaluation of the observational products point to merged data products (including rain gauge information) as providing the overall best results.

Description: Observational and modeling studies have described the relationships between convective/stratiform rain proportion and the vertical distributions of vertical motion, latent heating, and moistening in mesoscale convective systems. Therefore, remote sensing techniques which can quantify the relative areal proportion of convective and stratiform, rainfall can provide useful information regarding the dynamic and thermodynamic processes in these systems. In the present study, two methods for deducing the convective/stratiform areal extent of precipitation from satellite passive microwave radiometer measurements are combined to yield an improved method. If sufficient microwave scattering by ice-phase precipitating hydrometeors is detected, the method relies mainly on the degree of polarization in oblique-view, 85.5 GHz radiances to estimate the area fraction of convective rain within the radiometer footprint. In situations where ice scattering is minimal, the method draws mostly on texture information in radiometer imagery at lower microwave frequencies to estimate the convective area fraction. Based upon observations of ten convective systems over ocean and nine systems over land, instantaneous 0.5 degree resolution estimates of convective area fraction from the Tropical Rainfall Measuring Mission Microwave Imager (TRMM TMI) are compared to nearly coincident estimates from the TRMM Precipitation Radar (TRMM PR). The TMI convective area fraction estimates are slightly low-biased with respect to the PR, with TMI-PR correlations of 0.78 and 0.84 over ocean and land backgrounds, respectively. TMI monthly-average convective area percentages in the tropics and subtropics from February 1998 exhibit the greatest values along the ITCZ and in continental regions of the summer (southern) hemisphere. Although convective area percentages. from the TMI are systematically lower than those from the PR, monthly rain patterns derived from the TMI and PR rain algorithms are very similar. TMI rain depths are significantly higher than corresponding rain depths from the PR in the ITCZ, but are similar in magnitude elsewhere.

Description: We built a direct detection Doppler lidar based on the double-edge molecular technique and made the first molecular based wind measurements using the eyesafe 355 nm wavelength. Three etalon bandpasses are obtained with Step etalons on a single pair of etalon plates. Long-term frequency drift of the laser and the capacitively stabilized etalon is removed by locking the etalon to the laser frequency. We use a low angle design to avoid polarization effects. Wind measurements of 1 to 2 m/s accuracy are obtained to 10 km altitude with 5 mJ of laser energy, a 750s integration, and a 25 cm telescope. Good agreement is obtained between the lidar and rawinsonde measurements.

Description: We use clear sky heating rates to show that convective outflow in the tropics decreases rapidly with height between the 350 K and 360 K potential temperature surfaces (or between roughly 13 and 15 km). There is also a rapid fall-off in the pseudoequivalent potential temperature probability distribution of near surface air parcels between 350 K and 360 K. This suggests that the vertical variation of convective outflow in the upper tropical troposphere is to a large degree determined by the distribution of sub cloud layer entropy.

Description: A symposium celebrating the first 50 years of Dr. Joanne Simpson's career took place at the NASA/Goddard Space Flight Center from December 1 - 3, 1999. This symposium consisted of presentations that focused on: historical and personal points of view concerning Dr. Simpson's research career, her interactions with the American Meteorological Society, and her leadership in TRMM; scientific interactions with Dr. Simpson that influenced personal research; research related to observations and modeling of clouds, cloud systems and hurricanes; and research related to the Tropical Rainfall Measuring Mission (TRMM). There were a total of 36 presentations and 103 participants from the US, Japan and Australia. The specific presentations during the symposium are summarized in this paper.

Description: We present results of simulations of the distribution of 1809 keV radiation from the decay of Al-26 in the Galaxy. Recent observations of this emission line using the Gamma Ray Imaging Spectrometer (GRIS) have indicated that the bulk of the AL-26 must have a velocity of approx. 500 km/ s. We have previously shown that a velocity this large could be maintained over the 10(exp 6) year lifetime of the Al-26 if it is trapped in dust grains that are reaccelerated periodically in the ISM. Here we investigate whether a dust grain velocity of approx. 500 km/ s will produce a distribution of 1809 keV emission in latitude that is consistent with the narrow distribution seen by COMPTEL. We find that dust grain velocities in the range 275 - 1000 km/ s are able to reproduce the COMPTEL 1809 keV emission maps reconstructed using the Richardson-Lucy and Maximum Entropy image reconstruction methods while the emission map reconstructed using the Multiresolution Regularized Expectation Maximization algorithm is not well fit by any of our models. The Al-26 production rate that is needed to reproduce the observed 1809 keV intensity yields in a Galactic mass of Al-26 of approx. 1.5 - 2 solar mass which is in good agreement with both other observations and theoretical production rates.

Description: This study examines the uncertainty in forecasts of the January-February-March (JFM) mean extratropical circulation, and how that uncertainty is modulated by the El Nino/Southern Oscillation (ENSO). The analysis is based on ensembles of hindcasts made with an Atmospheric General Circulation Model (AGCM) forced with sea surface temperatures observed during; the 1983 El Nino and 1989 La Nina events. The AGCM produces pronounced interannual differences in the magnitude of the extratropical seasonal mean noise (intra-ensemble variability). The North Pacific, in particular, shows extensive regions where the 1989 seasonal mean noise kinetic energy (SKE), which is dominated by a "PNA-like" spatial structure, is more than twice that of the 1983 forecasts. The larger SKE in 1989 is associated with a larger than normal barotropic conversion of kinetic energy from the mean Pacific jet to the seasonal mean noise. The generation of SKE due to sub-monthly transients also shows substantial interannual differences, though these are much smaller than the differences in the mean flow conversions. An analysis of the Generation of monthly mean noise kinetic energy (NIKE) and its variability suggests that the seasonal mean noise is predominantly a statistical residue of variability resulting from dynamical processes operating on monthly and shorter times scales. A stochastically-forced barotropic model (linearized about the AGCM's 1983 and 1989 base states) is used to further assess the role of the basic state, submonthly transients, and tropical forcing, in modulating the uncertainties in the seasonal AGCM forecasts. When forced globally with spatially-white noise, the linear model generates much larger variance for the 1989 base state, consistent with the AGCM results. The extratropical variability for the 1989 base state is dominanted by a single eigenmode, and is strongly coupled with forcing over tropical western Pacific and the Indian Ocean, again consistent with the AGCM results. Linear calculations that include forcing from the AGCM variance of the tropical forcing and submonthly transients show a small impact on the variability over the Pacific/North American region compared with that of the base state differences.

Description: Using global rainfall and sea surface temperature (SST) data for the past two decades (1979-1998), we have investigated the intrinsic modes of Asian summer monsoon (ASM) and ENSO co-variability. Three recurring ASM rainfall-SST coupled modes were identified. The first is a basin scale mode that features SST and rainfall variability over the entire tropics (including the ASM region), identifiable with those occurring during El Nino or La Nina. This mode is further characterized by a pronounced biennial variation in ASM rainfall and SST associated with fluctuations of the anomalous Walker circulation that occur during El Nino/La Nina transitions. The second mode comprises mixed regional and basin-scale rainfall and SST signals, with pronounced intraseasonal and interannual variabilities. This mode features a SST pattern associated with a developing La Nina, with a pronounced low level anticyclone in the subtropics of the western Pacific off the coast of East Asia. The third mode depicts an east-west rainfall and SST dipole across the southern equatorial Indian Ocean, most likely stemming from coupled ocean-atmosphere processes within the ASM region. This mode also possesses a decadal time scale and a linear trend, which are not associated with El Nino/La Nina variability. Possible causes of year-to-year rainfall variability over the ASM and sub-regions have been evaluated from a reconstruction of the observed rainfall from singular eigenvectors of the coupled modes. It is found that while basin-scale SST can account for portions of ASM rainfall variability during ENSO events (up to 60% in 1998), regional processes can accounts up to 20-25% of the rainfall variability in typical non-ENSO years. Stronger monsoon-ENSO relationship tends to occur in the boreal summer immediately preceding a pronounced La Nina, i.e., 1998, 1988 and 1983. Based on these results, we discuss the possible impacts of the ASM on ENSO variability via the west Pacific anticyclone and articulate a hypothesis that anomalous wind forcings derived from the anticyclone may be instrumental in inducing a strong biennial modulation to natural ENSO cycles.

Description: In general, there are two broad scientific objectives when using cloud resolving models (CRMs or cloud ensemble models-CEMs) to study tropical convection. The first one is to use them as a physics resolving models to understand the dynamic and microphysical processes associated with the tropical water and energy cycles and their role in the climate system. The second approach is to use the CRMs to improve the representation of moist processes and their interaction with radiation in large-scale models. In order to improve the credibility of the CRMs and achieve the above goals, CRMs using identical initial conditions and large-scale influences need to produce very similar results. Two CRMs produced different statistical equilibrium (SE) states even though both used the same initial thermodynamic and wind conditions. Sensitivity tests to identify the major physical processes that determine the SE states for the different CRM simulations were performed. Their results indicated that atmospheric horizontal wind is treated quite differently in these two CRMs. The model that had stronger surface winds and consequently larger latent and sensible heat fluxes from the ocean produced a warmer and more humid modeled thermodynamic SE state. In addition, the domain mean thermodynamic state is more unstable for those experiments that produced a warmer and more humid SE state. Their simulated wet (warm and humid) SE states are thermally more stable in the lower troposphere (from the surface to 4-5 km in altitude). The large-scale horizontal advective effects on temperature and water vapor mixing ratio are needed when using CRMs to perform long-term integrations to study convective feedback under specified large-scale environments. In addition, it is suggested that the dry and cold SE state simulated was caused by enhanced precipitation but not enough surface evaporation. We find some problems with the interpretation of these three phenomena.

Description: The effectiveness of techniques for creating "bogus" vortices in numerical simulations of hurricanes is examined by using the Penn State/NCAR nonhydrostatic mesoscale model (MM5) and its adjoint system. A series of four-dimensional variational data assimilation (4-D VAR) experiments is conducted to generate an initial vortex for Hurricane Georges (1998) in the Atlantic Ocean by assimilating bogus sea-level pressure and surface wind information into the mesoscale numerical model. Several different strategies are tested for improving the vortex representation. The initial vortices produced by the 4-D VAR technique are able to reproduce many of the structural features of mature hurricanes. The vortices also result in significant improvements to the hurricane forecasts in terms of both intensity and track. In particular, with assimilation of only bogus sea-level pressure information, the response in the wind field is contained largely within the divergent component, with strong convergence leading to strong upward motion near the center. Although the intensity of the initial vortex seems to be well represented, a dramatic spin down of the storm occurs within the first 6 h of the forecast. With assimilation of bogus surface wind data only, an expected dominance of the rotational component of the wind field is generated, but the minimum pressure is adjusted inadequately compared to the actual hurricane minimum pressure. Only when both the bogus surface pressure and wind information are assimilated together does the model produce a vortex that represents the actual intensity of the hurricane and results in significant improvements to forecasts of both hurricane intensity and track.

Description: This paper presents viewgraphs on Thermal Vacuum Tests of the GLAS (Geoscience Laser Altimeter System) Propylene Loop Heat Pipe Development Model. The topics include: 1) Flight LHP System (Laser); 2) Test Design and Objectives; 3) DM (Development Model) LHP (Loop Heat Pipe) Test Design; 4) Starter Heater and Coupling Blocks; 5) CC Control Heaters and PRT; 6) Heater Plates (Shown in Reflux Mode); 7) Startup Tests; 8) CC Control Heater Power Tests for CC Temperature Control; and 9) Control Temperature Stability.

Description: This paper represents the first attempt to use TRMM rainfall information to estimate the four dimensional latent heating structure over the global tropics for February 1998. The mean latent heating profiles over six oceanic regions (TOGA COARE IFA, Central Pacific, S. Pacific Convergence Zone, East Pacific, Indian Ocean and Atlantic Ocean) and three continental regions (S. America, Central Africa and Australia) are estimated and studied. The heating profiles obtained from the results of diagnostic budget studies over a broad range of geographic locations are used to provide comparisons and indirect validation for the heating algorithm estimated heating profiles. Three different latent heating algorithms, the Goddard Convective-Stratiform (CSH) heating, the Goddard Profiling (GPROF) heating, and the Hydrometeor heating (HH) are used and their results are intercompared. The horizontal distribution or patterns of latent heat release from the three different heating retrieval methods are quite similar. They all can identify the areas of major convective activity (i.e., a well defined ITCZ in the Pacific, a distinct SPCZ) in the global tropics. The magnitude of their estimated latent heating release is also not in bad agreement with each other and with those determined from diagnostic budget studies. However, the major difference among these three heating retrieval algorithms is the altitude of the maximum heating level. The CSH algorithm estimated heating profiles only show one maximum heating level, and the level varies between convective activity from various geographic locations. These features are in good agreement with diagnostic budget studies. By contrast, two maximum heating levels were found using the GPROF heating and HH algorithms. The latent heating profiles estimated from all three methods can not show cooling between active convective events. We also examined the impact of different TMI (Multi-channel Passive Microwave Sensor) and PR (Precipitation Radar) rainfall information on latent heating structures.

Description: The NASA/GSFC Scanning Raman Lidar (SRL) was stationed on Andros Island in the Bahamas during August - September, 1998 as a part of the third Convection and Moisture Experiment (CAMEX-3) which focussed on hurricane development and tracking. During the period August 21 - 24, hurricane Bonnie passed near Andros Island and influenced the water vapor and cirrus cloud measurements acquired by the SRL. Two drying signatures related to the hurricane were recorded by the SRL (Scanning Raman Lidar) and other sensors. Cirrus cloud optical depths (at 351 nm) were also measured during this period. Optical depth values ranged from approximately 0.01 to 1.4. The influence of multiple scattering on these optical depth measurements was studied with the conclusion that the measured values of optical depth are less than the actual value by up to 20% . The UV/IR cirrus cloud optical depth ratio was estimated based on a comparison of lidar and GOES measurements. Simple radiative transfer model calculations compared with GOES satellite brightness temperatures indicate that satellite radiances are significantly affected by the presence of cirrus clouds if IR optical depths are approximately 0.02 or greater. This has implications for satellite cirrus detection requirements.

Description: Central Florida is the ideal test laboratory for studying convergence zone-induced convection. The region regularly experiences sea breeze fronts and rainfall-induced outflow boundaries. The focus of this study is the common yet poorly-studied convergence zone established by the interaction of the sea breeze front and an outflow boundary. Previous studies have investigated mechanisms primarily affecting storm initiation by such convergence zones. Few have focused on rainfall morphology yet these storms contribute a significant amount precipitation to the annual rainfall budget. Low-level convergence and mid-tropospheric moisture have both been shown to correlate with rainfall amounts in Florida. Using 2D and 3D numerical simulations, the roles of low-level convergence and mid-tropospheric moisture in rainfall evolution are examined. The results indicate that time-averaged, vertical moisture flux (VMF) at the sea breeze front/outflow convergence zone is directly and linearly proportional to initial condensation rates. This proportionality establishes a similar relationship between VMF and initial rainfall. Vertical moisture flux, which encompasses depth and magnitude of convergence, is better correlated to initial rainfall production than surface moisture convergence. This extends early observational studies which linked rainfall in Florida to surface moisture convergence. The amount and distribution of mid-tropospheric moisture determines how rainfall associated with secondary cells develop. Rainfall amount and efficiency varied significantly over an observable range of relative humidities in the 850- 500 mb layer even though rainfall evolution was similar during the initial or "first-cell" period. Rainfall variability was attributed to drier mid-tropospheric environments inhibiting secondary cell development through entrainment effects. Observationally, 850-500 mb moisture structure exhibits wider variability than lower level moisture, which is virtually always present in Florida. A likely consequence of the variability in 850-500 moisture is a stronger statistical correlation to rainfall, which observational studies have noted. The study indicates that vertical moisture flux forcing at convergence zones is critical in determining rainfall in the initial stage of development but plays a decreasing role in rainfall evolution as the system matures. The mid-tropospheric moisture (e.g. environment) plays an increasing role in rainfall evolution as the system matures. This suggests the need to improve measurements of magnitude/depth of convergence and mid-tropospheric moisture distribution. It also highlights the need for better parameterization of entrainment and vertical moisture distribution in larger-scale models.

Description: Chao's numerical and theoretical work on multiple quasi-equilibria of the intertropical convergence zone (ITCZ) and the origin of monsoon onset is extended to solve two additional puzzles. One is the highly nonlinear dependence on latitude of the "force" acting on the ITCZ due to earth's rotation, which makes the multiple quasi-equilibria of the ITCZ and monsoon onset possible. The other is the dramatic difference in such dependence when different cumulus parameterization schemes are used in a model. Such a difference can lead to a switch between a single ITCZ at the equator and a double ITCZ, when a different cumulus parameterization scheme is used. Sometimes one of the double ITCZ can diminish and only the other remain, but still this can mean different latitudinal locations for the single ITCZ. A single idea based on two off-equator attractors for the ITCZ, due to earth's rotation and symmetric with respect to the equator, and the dependence of the strength and size of these attractors on the cumulus parameterization scheme solves both puzzles. The origin of these rotational attractors, explained in Part I, is further discussed. The "force" acting on the ITCZ due to earth's rotation is the sum of the "forces" of the two attractors. Each attractor exerts on the ITCZ a "force" of simple shape in latitude; but the sum gives a shape highly varying in latitude. Also the strength and the domain of influence of each attractor vary, when change is made in the cumulus parameterization. This gives rise to the high sensitivity of the "force" shape to cumulus parameterization. Numerical results, of experiments using Goddard's GEOS general circulation model, supporting this idea are presented. It is also found that the model results are sensitive to changes outside of the cumulus parameterization. The significance of this study to El Nino forecast and to tropical forecast in general is discussed.

Description: The 1997-1999 ENSO period was very powerful, but also well observed. Multiple satellite rainfall estimates combined with gauge observations allow for a quantitative analysis of precipitation anomalies in the tropics and elsewhere accompanying the 1997-99 ENSO cycle. An examination of the evolution of the El Nino and accompanying precipitation anomalies revealed that a dry Maritime Continent preceded the formation of positive SST anomalies in the eastern Pacific Ocean. 30-60 day oscillations in the winter of 1996/97 may have contributed to this lag relationship. Furthermore, westerly wind burst events may have maintained the drought over the Maritime Continent. The warming of the equatorial Pacific was then followed by an increase in convection. A rapid transition from El Nino to La Nina occurred in May 1998, but as early as October-November 1997 precipitation indices captured substantial changes in Pacific rainfall anomalies. The global precipitation patterns for this event were in good agreement with the strong consistent ENSO-related precipitation signals identified in earlier studies. Differences included a shift in precipitation anomalies over Africa during the 1997-98 El Nino and unusually wet conditions over northeast Australia during the later stages of the El Nino. Also, the typically wet region in the north tropical Pacific was mostly dry during the 1998-99 La Nina. Reanalysis precipitation was compared to observations during this time period and substantial differences were noted. In particular, the model had a bias towards positive precipitation anomalies and the magnitudes of the anomalies in the equatorial Pacific were small compared to the observations. Also, the evolution of the precipitation field, including the drying of the Maritime Continent and eastward progression of rainfall in the equatorial Pacific was less pronounced for the model compared to the observations.

Description: The predictability of the 1997 and 1998 south Asian summer monsoon winds is examined from an ensemble of 10 Atmospheric General Circulation Model (AGCM) simulations with prescribed sea surface temperatures (SSTs) and soil moisture, The simulations are started in September 1996 so that they have lost all memory of the atmospheric initial conditions for the periods of interest. The model simulations show that the 1998 monsoon is considerably more predictable than the 1997 monsoon. During May and June of 1998 the predictability of the low-level wind anomalies is largely associated with a local response to anomalously warm Indian Ocean SSTs. Predictability increases late in the season (July and August) as a result of the strengthening of the anomalous Walker circulation and the associated development of easterly low level wind anomalies that extend westward across India and the Arabian Sea. During these months the model is also the most skillful with the observations showing a similar late-season westward extension of the easterly CD wind anomalies. The model shows little predictability or skill in the low level winds over southeast Asia during, 1997. Predictable wind anomalies do occur over the western Indian Ocean and Indonesia, however, over the Indian Ocean they are a response to SST anomalies that were wind driven and they show no skill. The reduced predictability in the low level winds during 1997 appears to be the result of a weaker (compared with 1998) simulated anomalous Walker circulation, while the reduced skill is associated with pronounced intraseasonal activity that is not well captured by the model. Remarkably, the model does produce an ensemble mean Madden-Julian Oscillation (MJO) response that is approximately in phase with (though weaker than) the observed MJ0 anomalies. This is consistent with the idea that SST coupling may play an important role in the MJO.

Description: This article is partly a review and partly a new research paper on monsoon-ENSO relationship. The paper begins with a discussion of the basic relationship between the Indian monsoon and ENSO dating back to the work of Sir Gilbert Walker up to research results in more recent years. Various factors that may affect the monsoon-ENSO, relationship, including regional coupled ocean-atmosphere processes, Eurasian snow cover, land-atmosphere hydrologic feedback, intraseasonal oscillation, biennial variability and inter-decadal variations, are discussed. The extreme complex and highly nonlinear nature of the monsoon-ENSO relationship is stressed. We find that for regional impacts on the monsoon, El Nino and La Nina are far from simply mirror images of each other. These two polarities of ENSO can have strong or no impacts on monsoon anomalies depending on the strength of the intraseasonal oscillations and the phases of the inter-decadal variations. For the Asian-Australian monsoon (AAM) as a whole, the ENSO impact is effected through a east-west shift in the Walker Circulation. For rainfall anomalies over specific monsoon areas, regional processes play important roles in addition to the shift in the Walker Circulation. One of the key regional processes identified for the boreal summer monsoon is the anomalous West Pacific Anticyclone (WPA). This regional feature has similar signatures in interannual and intraseasonal time scales and appears to determine whether the monsoon-ENSO relationship is strong or weak in a given year. Another important regional feature includes a rainfall and SST dipole across the Indian Ocean, which may have strong impact on the austral summer monsoon. Results are shown indicating that monsoon surface wind forcings may induce a strong biennial signal in ENSO and that strong monsoon-ENSO coupling may translate into pronounced biennial variability in ENSO. Finally, a new paradigm is proposed for the study of monsoon variability. This paradigm provides a unified framework in which monsoon predictability, the role of regional vs. basin-scale processes, its relationship with different climate subsystems, and causes of secular changes in monsoon-ENSO relationship can be investigated.

Description: Abstract A 1999 study reports an advancement of spring in Europe by 0.2 days per year in the 30 years since 1960. Our analysis indicates that this trend results directly from a change in the late-winter surface winds over the eastern North Atlantic: the southwesterly direction became more dominant, and the speed of these southwesterlies increased slightly. Splitting the 52-year NCEP reanalysis dataset into the First Half, FH (1948-1973)), and the Second Half, SH (1974-1999), we analyze the wind direction for the February mean at three sites at 45N: site A at 30W, site B at 20W, and site C at 10W. The incidence (number of years) of the southwesterlies in SH Vs. (FH) at these sites respectively increased in SH as follows: 24(18), 19(12), 14(l 1); whereas the incidence of northeasterlies decreased: 0(2), 1(2), and 1(6). When the February mean wind is southwesterly, the monthly mean sensible heat flux from the ocean at these sites takes zero or slightly negative values, that is, the surface air is warmer than the ocean. Analyzing the scenario in the warm late winter 1990, we observe that the sensible heat flux from the ocean surface in February 1990 shows a "tongue" of negative values extending southwest from southern England to 7N. This indicates that the source of the maritime air advected into Europe lies to the south of the "tongue." Streamline analysis suggests that the Southwestern or southcentral North Atlantic is the source. For February 1990, we find strong, ascending motions over Europe at 700 mb, up to -0.4 Pa/s as monthly averages. Associated with the unstable low-levels of the troposphere are positive rain and cloud anomalies. Thus, positive in situ feedback over land in late winter (when shortwave absorption is not significant) apparently further enhances the surface temperature through an increase in the greenhouse effect due to increased water vapor and cloudiness.

Description: It is a long-held fundamental belief that the basic cause of a monsoon is land-sea thermal contrast on the continental scale. Through general circulation model experiments we demonstrate that this belief should be changed. The Asian and Australian summer monsoon circulations are largely intact in an experiment in which Asia, maritime continent, and Australia are replaced by ocean. It is also shown that the change resulting from such replacement is in general due more to the removal of topography than to the removal of land-sea contrast. Therefore, land-sea contrast plays only a minor modifying role in Asian and Australian summer monsoons. This also happens to the Central American summer monsoon. However, the same thing cannot be said of the African and South American summer monsoons. In Asian and Australian winter monsoons land-sea contrast also plays only a minor role. Our interpretation for the origin of monsoon is that the summer monsoon is the result of ITCZ's (intertropical convergence zones) peak being substantially (more than 10 degrees) away from the equator. The origin of the ITCZ has been previously interpreted by Chao. The circulation around thus located ITCZ, previously interpreted by Chao and Chen through the modified Gill solution and briefly described in this paper, explains the monsoon circulation. The longitudinal location of the ITCZs is determined by the distribution of surface conditions. ITCZ's favor locations of higher SST as in western Pacific and Indian Ocean, or tropical landmass, due to land-sea contrast, as in tropical Africa and South America. Thus, the role of landmass in the origin of monsoon can be replaced by ocean of sufficiently high SST. Furthermore, the ITCZ circulation extends into the tropics in the other hemisphere to give rise to the winter monsoon circulation there. Also through the equivalence of land-sea contrast and higher SST, it is argued that the basic monsoon onset mechanism proposed by Chao is valid for all monsoons.

Description: We evaluated the impact of several newly available sources of meteorological data on mesoscale model forecasts of precipitation produced by the extra-tropical cyclone that struck Florida on February 2, 1998. Precipitation distributions of convective rainfall events were derived from Special Sensor Microwave Imager (SSM/I) and Multi-Channel Passive Microwave Sensor (TMI) microwave radiometric data by means of the Goddard PROFiling (GPROF) algorithm. Continuous lightning distributions were obtained from sferics measurements obtained from a network of VLF radio receivers. Histograms of coincident sferics frequency distributions were matched to those of precipitation to derive bogus convective rainfall rates from the continuously available sferics measurements. SSM/I and TMI microwave data were used to derive Integrated Precipitable Water (IPW) distributions. The TMI also provided sea surface temperatures (SSTS) of the Loop Current and Gulf Stream with improved structural detail. A series of experiments assimilated IPW and latent heating from the bogus convective rainfall for six-hours in the MM5 mesoscale forecast model to produce nine-hour forecasts of all rainfall as well as other weather parameters. Although continuously assimilating latent heating only slightly improved the surface pressure distribution forecast, it significantly improved the precipitation forecasts. Correctly locating convective rainfall was found critical for assimilating latent heating in the forecast model, but measurement of the rainfall intensity proved to be less important. The improved SSTs also had a positive impact on rainfall forecasts for this case. Assimilating bogus rainfall in the model produced nine-hour forecasts of radar reflectivity distributions that agreed well with coincident observations from the TRMM spaceborne precipitation radar, ground based radar and spaceborne microwave measurements.

Description: The NASA/GSFC Scanning Raman Lidar (SRL) was stationed on Andros Island in the Bahamas during August - September, 1998 as a part of the third Convection and Moisture Experiment (CAMEX-3) which focussed on hurricane development and tracking. During the period August 21 - 24, hurricane Bonnie passed near Andros Island and influenced the water vapor and cirrus cloud measurements acquired by the SRL. Two drying signatures related to the hurricane were recorded by the SRL and other sensors. Cirrus cloud optical depths (at 351 nm) were also measured during this period. Optical depth values ranged from less than 0.01 to 1.5. The influence of multiple scattering on these optical depth measurements was studied. A correction technique is presented which minimizes the influences of multiple scattering and derives information about cirrus cloud optical and physical properties. The UV/IR cirrus cloud optical depth ratio was estimated based on a comparison of lidar and GOES measurements. Simple radiative transfer model calculations compared with GOES satellite brightness temperatures indicate that satellite radiances are significantly affected by the presence of cirrus clouds if IR optical depths are approximately 0.005 or greater. Using the ISCCP detection threshold for cirrus clouds on the GOES data presented here, a high bias of up to 40% in the GOES precipitable water retrieval was found.

Description: Quantitative use of satellite-derived maps of monthly rainfall requires some measure of the accuracy of the satellite estimates. The rainfall estimate for a given map grid box is subject to both remote-sensing error and, in the case of low-orbiting satellites, sampling error due to the limited number of observations of the grid box provided by the satellite. A simple model of rain behavior predicts that Root-mean-square (RMS) random error in grid-box averages should depend in a simple way on the local average rain rate, and the predicted behavior has been seen in simulations using surface rain-gauge and radar data. This relationship was examined using satellite SSM/I data obtained over the western equatorial Pacific during TOGA COARE. RMS error inferred directly from SSM/I rainfall estimates was found to be larger than predicted from surface data, and to depend less on local rain rate than was predicted. Preliminary examination of TRMM microwave estimates shows better agreement with surface data. A simple method of estimating rms error in satellite rainfall estimates is suggested, based on quantities that can be directly computed from the satellite data.

Description: The tropical cyclone rainfall climatology study that was performed for the North Pacific was extended to the North Atlantic. Similar to the North Pacific tropical cyclone study, mean monthly rainfall within 444 km of the center of the North Atlantic tropical cyclones (i.e., that reached storm stage and greater) was estimated from passive microwave satellite observations during, an eleven year period. These satellite-observed rainfall estimates were used to assess the impact of tropical cyclone rainfall in altering the geographical, seasonal, and inter-annual distribution of the North Atlantic total rainfall during, June-November when tropical cyclones were most abundant. The main results from this study indicate: 1) that tropical cyclones contribute, respectively, 4%, 3%, and 4% to the western, eastern, and entire North Atlantic; 2) similar to that observed in the North Pacific, the maximum in North Atlantic tropical cyclone rainfall is approximately 5 - 10 deg poleward (depending on longitude) of the maximum non-tropical cyclone rainfall; 3) tropical cyclones contribute regionally a maximum of 30% of the total rainfall 'northeast of Puerto Rico, within a region near 15 deg N 55 deg W, and off the west coast of Africa; 4) there is no lag between the months with maximum tropical cyclone rainfall and non-tropical cyclone rainfall in the western North Atlantic, while in the eastern North Atlantic, maximum tropical cyclone rainfall precedes maximum non-tropical cyclone rainfall; 5) like the North Pacific, North Atlantic tropical cyclones Of hurricane intensity generate the greatest amount of rainfall in the higher latitudes; and 6) warm ENSO events inhibit tropical cyclone rainfall.

Description: The mechanism of the quasi-biennial tendency in El Nino Southern Oscillation (ENSO)-monsoon coupled system is investigated using an intermediate coupled model. The monsoon wind forcing is prescribed as a function of Sea Surface Temperature (SST) anomalies based on the relationship between zonal wind anomalies over the western Pacific to sea level change in the equatorial eastern Pacific. The key mechanism of quasi-biennial tendency in El Nino evolution is found to be in the strong coupling of ENSO to monsoon wind forcing over the western Pacific. Strong boreal summer monsoon wind forcing, which lags the maximum SST anomaly in the equatorial eastern Pacific approximately 6 months, tends to generate Kelvin waves of the opposite sign to anomalies in the eastern Pacific and initiates the turnabout in the eastern Pacific. Boreal winter monsoon forcing, which has zero lag with maximum SST in the equatorial eastern Pacific, tends to damp the ENSO oscillations.

Description: Idealized numerical simulations are performed with a coupled atmosphere/land-surface model to identify the roles of initial soil moisture, coastline curvature, and land breeze circulations on sea breeze initiated precipitation. Data collected on 27 July 1991 during the Convection and Precipitation Electrification Experiment (CAPE) in central Florida are used. The 3D Goddard Cumulus Ensemble (GCE) cloud resolving model is coupled with the Goddard Parameterization for Land-Atmosphere-Cloud Exchange (PLACE) land surface model, thus providing a tool to simulate more realistically land-surface/atmosphere interaction and convective initiation. Eight simulations are conducted with either straight or curved coast-lines, initially homogeneous soil moisture or initially variable soil moisture, and initially homogeneous horizontal winds or initially variable horizontal winds (land breezes). All model simulations capture the diurnal evolution and general distribution of sea-breeze initiated precipitation over central Florida. The distribution of initial soil moisture influences the timing, intensity and location of subsequent precipitation. Soil moisture acts as a moisture source for the atmosphere, increases the connectively available potential energy, and thus preferentially focuses heavy precipitation over existing wet soil. Strong soil moisture-induced mesoscale circulations are not evident in these simulations. Coastline curvature has a major impact on the timing and location of precipitation. Earlier low-level convergence occurs inland of convex coastlines, and subsequent precipitation occurs earlier in simulations with curved coastlines. The presence of initial land breezes alone has little impact on subsequent precipitation. however, simulations with both coastline curvature and initial land breezes produce significantly larger peak rain rates due to nonlinear interactions.

Description: The numerical simulation of precipitation helps scientists understand the complex mechanisms that determine how and why rainfall is distributed across the globe. Simulation aids in the development of forecastin,g efforts that inform policies regarding the management of water resources. Precipitation modeling also provides short-term warnings, for emergencies such as flash floods and mudslides. Just as precipitation modeling can warn of an impending abundance of rainfall, it can help anticipate the absence of rainfall in drought. What constitutes a drought? A meteorological drought simply means that an area is getting a significantly lower amount of rain than usual over a prolonged period of time and an agricultural drought is based on the level of soil moisture.

Description: This presentation describes the state of transonic flow simulation using nonlinear potential methods for external aerodynamic applications. The presentation begins with a review of the various potential equation forms (with emphasis on the full potential equation) and includes a discussion of pertinent mathematical characteristics and all derivation assumptions. Impact of the derivation assumptions on simulation accuracy, especially with respect to shock wave capture, is discussed. Key characteristics of all numerical algorithm types used for solving nonlinear potential equations, including steady, unsteady, space marching, and design methods, are described. Both spatial discretization and iteration scheme characteristics are examined. Numerical results for various aerodynamic applications are included throughout the presentation to highlight key discussion points. The presentation ends with concluding remarks and recommendations for future work. Overall. nonlinear potential solvers are efficient, highly developed and routinely used in the aerodynamic design environment for cruise conditions. Published by Elsevier Science Ltd. All rights reserved.

Description: Subgrid analysis of a transitional temporal mixing layer with evaporating droplets has been performed using a direct numerical simulation (DNS) database. The DNS is for a Reynolds number (based on initial vorticity thickness) of 600, with droplet mass loading of 0.2. The gas phase is computed using a Eulerian formulation, with Lagrangian droplet tracking. Since Large Eddy Simulation (LES) of this flow requires the computation of unfiltered gas-phase variables at droplet locations from filtered gas-phase variables at the grid points, it is proposed to model these by assuming the gas-phase variables to be given by the filtered variables plus a correction based on the filtered standard deviation, which can be computed from the sub-grid scale (SGS) standard deviation. This model predicts unfiltered variables at droplet locations better than simply interpolating the filtered variables. Three methods are investigated for modeling the SGS standard deviation: Smagorinsky, gradient and scale-similarity. When properly calibrated, the gradient and scale-similarity methods give results in excellent agreement with the DNS.

Description: Future generations of ultra high bypass-ratio jet engines will require far higher pressure ratios and operating temperatures than those of current engines. For the foreseeable future, engine materials will not be able to withstand the high temperatures without some form of cooling. In particular the turbine blades, which are under high thermal as well as mechanical loads, must be cooled. Cooling of turbine blades is achieved by bleeding air from the compressor stage of the engine through complicated internal passages in the turbine blades (internal cooling, including jet-impingement cooling) and by bleeding small amounts of air into the boundary layer of the external flow through small discrete holes on the surface of the blade (film cooling and transpiration cooling). The cooling must be done using a minimum amount of air or any increases in efficiency gained through higher operating temperature will be lost due to added load on the compressor stage. Turbine cooling schemes have traditionally been based on extensive empirical data bases, quasi-one-dimensional computational fluid dynamics (CFD) analysis, and trial and error. With improved capabilities of CFD, these traditional methods can be augmented by full three-dimensional simulations of the coolant flow to predict in detail the heat transfer and metal temperatures. Several aspects of turbine coolant flows make such application of CFD difficult, thus a highly effective CFD methodology must be used. First, high resolution of the flow field is required to attain the needed accuracy for heat transfer predictions, making highly efficient flow solvers essential for such computations. Second, the geometries of the flow passages are complicated but must be modeled accurately in order to capture all important details of the flow. This makes grid generation and grid quality important issues. Finally, since coolant flows are turbulent and separated the effects of turbulence must be modeled with a low Reynolds number turbulence model to accurately predict details of heat transfer.

Description: Aerosol is any small particle of matter that rests suspended in the atmosphere. Natural sources, such as deserts, create some aerosols; consumption of fossil fuels and industrial activity create other aerosols. All the microscopic aerosol particles add up to a large amount of material floating in the atmosphere. You can see the particles in the haze that floats over polluted cities. Beyond this visible effect, aerosols can actually lower temperatures. They do this by blocking, or scattering, a portion of the sun's energy from reaching the surface. Because of this influence, scientists study the physical properties of atmospheric aerosols. Reliable numerical models for atmospheric aerosols play an important role in research.

Description: During the TEFLUN-B (Texas-Florida under-flights for TRMM) field experiment of August-September, 1998, a number of ER-2 aircraft flights with a host of microwave instruments were conducted over many convective storms, including some hurricanes, in the coastal region of Florida and Texas. These instruments include MIR (Millimeter-wave Imaging Radiometer), AMPR (Advanced Microwave Precipitation Radiometer), and EDOP (ER-2 Doppler Radar). EDOP is operated at the frequency of 9.7 GHz, while the AMPR and the MIR together give eleven channels of radiometric measurements in the frequency range of 10-340 GHz. The concurrent measurements from these instruments provide unique data sets for studying the details of the microphysics of hydrometeors. Preliminary examination of these data sets shows features that are generally well understood; i.e., radiometric measurements at frequencies less than or equal to 37 GHz mainly respond to rain, while those at frequencies greater than or equal to 150 GHz, to ice particles above the freezing level. Model calculations of brightness temperature and radar reflectivity are performed and results compared with these measurements. For simplicity the analysis is limited to the anvil region of the storms where hydrometeors are predominantly frozen. Only one ice particle size distribution is examined in the calculations of brightness temperature and radar reflectivity in this initial study. Estimation of ice water path is made based on the best agreement between the measurements and calculations of brightness temperature and reflectivity. Problems associated with these analyses and measurement accuracy will be discussed.

Description: GLOW (Goddard Lidar Observatory for Winds) is a mobile Doppler lidar system which uses direct detection Doppler lidar techniques to measure wind profiles from the surface into the lower stratosphere. The system is contained in a modified van to allow deployment in field operations. The lidar system uses a Nd:YAG laser transmitter to measure winds using either aerosol backscatter at 1064 nm or molecular backscatter at 355 nm. The receiver telescope is a 45 cm Dall-Kirkham which is fiber coupled to separate Doppler receivers, one optimized for the aerosol backscatter wind measurement and another optimized for the molecular backscatter wind measurement. The receivers are implementations of the 'double edge' technique and use high spectral resolution Fabry-Perot etalons to measure the Doppler shift. A 45 cm aperture azimuth-over-elevation scanner is mounted on the roof of the van to allow full sky access and a variety of scanning options. GLOW is intended to be used as a deployable field system for studying atmospheric dynamics and transport and can also serve as a testbed to evaluate candidate technologies developed for use in future spaceborne systems. In addition, it can be used for calibration/validation activities following launch of spaceborne wind lidar systems. A description of the mobile system is presented along with the examples of lidar wind profiles obtained with the system.

Description: The emphasis of combustion research efforts at NASA Glenn Research Center (GRC) is on collaborating with industry to design and test gas-turbine combustors and subcomponents for both sub- and supersonic applications. These next-generation aircraft combustors are required to meet strict international environmental restrictions limiting emissions. To meet these goals, innovative combustor concepts require operation at temperatures and pressures far exceeding those of current designs. New and innovative diagnostic tools are necessary to characterize these flow streams since existing methods are inadequate. The combustion diagnostics team at GRC has implemented a suite of highly sensitive, nonintrusive optical imaging methods to diagnose the flowfields of these new engine concepts. By using optically accessible combustors and flame-tubes, imaging of fuel and intermediate combustion species via planar laser-induced fluorescence (PLIF) at realistic pressures are now possible. Direct imaging of the fuel injection process through both planar Mie scattering and PLIF methods is also performed. Additionally, a novel combination of planar fuel fluorescence imaging and computational analysis allows a 3-D examination of the flowfield, resulting in spatially and temporally resolved fuel/air volume distribution maps. These maps provide detailed insight into the fuel injection process at actual conditions, thereby greatly enhancing the evaluation of fuel injector performance and other combustion phenomena. Stable species such as CO2, O2, N2, H2O, and hydrocarbons are also investigated by a newly demonstrated 1-D, spontaneous Raman spectroscopic method. This visible wavelength Raman technique allows the acquisition of quantitative, stable species concentration measurements from the flow.

Description: Supernovae play a key role in the dynamics, structure, and chemical evolution of galaxies. The massive stars that end their lives as supernovae live for short enough times that many are still associated with dusty star formation regions when they explode, making them difficult to observe at visible wavelengths. In active star forming regions (galactic nuclei and starburst regions), dust extinction is especially severe. Thus, determining the supernova rate in active star forming regions of galaxies, where the supernova rate can be one or two orders of magnitude higher than the average, has proven to be difficult. From observations of SN1987A, we know that the [NiII] 6.63 micrometer emission line was the strongest line in the infrared spectrum for a period of a year and half after th explosion. Since dust extinction is much less at 6.63 micrometers than at visible wavelengths (A(sub 6.63)/A(sub V) = 0.025), the [NiII] line can be used as a sensitive probe for the detection of recent supernovae. We have observed a sample of starburst galaxies at 6.63 micrometers using ISOCAM to search for the [NiII] emission line characteristic of recent supernovae. We did not detect any [NiII] line emission brighter than a 5-sigma limit of 5 mJy. We can set upper limits to the supernova rate in our sample, scaled ot the rate in M82, of less than 0.3 per year at the 90% confidence level using Bayesian methods. Assuming that a supernova would have a [NiII] line with the same luminosity as observed in SN1987A, we find less than 0.09 and 0.15 per year at the 50% and 67% confidence levels. These rates are somewhat less if a more normal type II supernovae has a [NiII] line luminosity greater than the line in SN1987A.

Description: This document is a presentation in viewgraph format which reviews the laboratory facilities and their construction for the International Space Station(ISS). Graphic displays of the ISS are included, with special interest in the facilities available on the US Destiny module and other modules which will be used in the study of fluid physics on the ISS. There are also pictures and descriptions of various components of the Fluids and Combustion Facility.

Description: Application of gas-liquid two-phase flow systems for space-based thermal management and for the HEDS program demands a precise control of bubble size distribution in liquid. The necessity of bulk liquid motion for controlling bubble size and frequency in the space environment has been suggested by recent studies on pool, forced convection boiling and bubble formation in flowing liquid. The present work, consisting of two parts, explores bubble generation at wall in a cross-flowing liquid, i.e., in a forced convection boiling configuration. A schematic is shown. The first part looks into the bubble formation process under isothermal conditions in a reduced gravity environment, by injecting gas through a hole in the wall of a flowing liquid channel. In the latter part with channel wall heating, flow and temperature fields near a single bubble are studied under normal (1-g) and micro-gravity (mu-g) conditions. The objective of the isothermal experiments is to experimentally investigate the effects of liquid cross-flow velocity, gas flow rate, and orifice diameter on bubble formation. Data were taken mainly under reduced gravity conditions but some data were taken in normal gravity for comparison. The reduced gravity experiment was conducted aboard the NASA DC-9 Reduced Gravity Aircraft. The results show that the process of bubble formation and detachment depends on gravity, the orifice diameter (D(sub N)), the gas flow rate (Q(sub g)), and the liquid cross-flow velocity (U(sub L)). The reduced gravity data are shown. The data are analyzed based on a force balance, and two different detachment mechanisms are identified. When the gas momentum is large, the bubble detaches from the injection orifice as the gas momentum overcomes the attaching effects of liquid drag and inertia. The surface tension force is much reduced because a large part of the bubble pinning edge at the orifice is lost as the bubble axis is tilted by the liquid flow. When the gas momentum is small, the force balance in the liquid flow direction is important, and the bubble detaches when the bubble axis inclination exceeds a certain angle. With wall heating, liquid motion around an air bubble in cross-flowing 2cs silicone oil is experimentally investigated in 1-g. A spectral element based steady 2D numerical model is also developed. The traces of particles from experimental flow visualization and the corresponding computed streamlines are shown. At the upstream side of the bubble facing the cross-flow, thermocapillary and forced convection create liquid motion away from the wall, up along the surface. At the downstream side, a competing interaction between the two creates a recirculation cell, causing the bulk liquid to stagnate on the surface and separate thereafter. The important dimensionless parameters are - Surface tension and local cross-flow Reynolds numbers R(sub sigma)(U(sub ref)a/v) and Re(sub loc)(U(sub L)a/v), respectively based on reference thermocapillary U(sub ref)(sigma(sub T)Delta.T/mu, Delta T=T(sub wall)-T(sub liquid)) and local cross-flow velocity U(sub L), Prandtl number P(sub r) and Grashoff number Gr(rho.g.beta.DeltaTa(sup 3)/mu.v). Variation of the stagnation point with R(sub sigma) and Re(sub loc) is shown. Also shown are good agreement between experimental and numerical results in 1-g. The computational model is extended to mu-g condition to investigate temperature and velocity on the bubble surface, stagnation and reattachment points of the recirculation cell and wall heat transfer. It is observed that wall heating significantly alters the flow field around the bubble and thus the forces acting on the bubble, which control its detachment. Thus a combination of heating and liquid cross-flow can be utilized to precisely control bubble formation in a mu-g environment.

Description: Solid oxide electrolyzers, such as electrolysis cells utilizing yttria-stabilized zirconia, can produce oxygen from Mars atmospheric carbon dioxide and reject carbon monoxide and unreacted carbon dioxide in a separate stream. The oxygen-production process has been shown to be far more efficient if the high-pressure, unreacted carbon dioxide can be separated and recycled back into the feed stream. Additionally, the mass of the adsorption compressor can be reduced. Also, the carbon monoxide by-product is a valuable fuel for space exploration and habitation, with applications from fuel cells to production of hydrocarbons and plastics. In our research, we will design, construct, and test an innovative, robust, low mass, low power separation device that can recover carbon dioxide and carbon monoxide for Mars ISRU. Such fundamental process technology, involving gas-solid phase separation in a reduced gravitational environment, will help to enable Human Exploration and Development of Space. The separation device will be scaled to operate with a CO2 sorption compressor and a zirconia electrolysis device built at the NASA Ames Research Center and the University of Arizona, respectively. In our research, we will design, construct, and test an innovative, robust, low mass, low power separation device that can recover carbon dioxide and carbon monoxide for Mars ISRU, Such fundamental process technology, involving gas-solid phase separation in a reduced gravitational environment, will help to enable Human Exploration and Development of Space. The separation device will be scaled to operate with a CO2 sorption compressor and a zirconia electrolysis device built at the NASA Ames Research Center and the University of Arizona, The separation device will be scaled to operate with a CO2 sorption compressor and a zirconia electrolysis device built at the NASA Ames Research Center and the University of Arizona, Research needs for the design shown are as follows: (1) The best adsorbent for the process must be determined. (2) Adsorption isotherms must be measured, both for pure components and mixtures. (3) Mathematical modeling must be performed to provide a solid framework for design. (4) The separation system must be constructed and tested. (5) System integration must be studied.

Description: Aggregates were observed to form very suddenly in a lab-contained dust cloud, transforming (within seconds) an opaque monodispersed cloud into a clear volume containing rapidly-settling, long hair-like aggregates. The implications of such a "phase change" led to a series of experiments progressing from the lab, to KC-135, followed by micro-g flights on USML-1 and USML-2, and now EGM slated for Space Station. We attribute the sudden "collapse" of a cloud to the effect of dipoles. This has significant ramifications for all types of cloud systems, and additionally implicates dipoles in the processes of cohesion and adhesion of granular matter. Notably, there is the inference that like-charged grains need not necessarily repel if they are close enough together: attraction or repulsion depends on intergranular distance (the dipole being more powerful at short range), and the D/M ratio for each grain, where D is the dipole moment and M is the net charge. We discovered that these ideas about dipoles, the likely pervasiveness of them in granular material, the significance of the D/M ratio, and the idea of mixed charges on individual grains resulting from tribological processes --are not universally recognized in electrostatics, granular material studies, and aerosol science, despite some early seminal work in the literature, and despite commercial applications of dipoles in such modern uses as "Krazy Glue", housecleaning dust cloths, and photocopying. The overarching goal of EGM is to empirically prove that (triboelectrically) charged dielectric grains of material have dipole moments that provide an "always attractive" intergranular force as a result of both positive and negative charges residing on the surfaces of individual grains. Microgravity is required for this experiment because sand grains can be suspended as a cloud for protracted periods, the grains are free to rotate to express their electrostatic character, and Coulombic forces are unmasked. Suspended grains will be "interrogated" by applied electrical fields. In one module, grains will be immersed in an inhomogeneous electric field and allowed to be attracted towards or repelled from the central electrode of the module: part of the grain's speed will be a function of its net charge (monopole), part will be a function of the dipole. Observed grain position vs. time will provide a curve that can be deconvolved into the dipole and monopole forces responsible, since both have distinctive radial dependencies. In a second approach, the inhomogeneous field will be alternated at low frequency (e.g., every 5-10 seconds) so that the grains are alternately attracted and repelled from the center of the field. The resulting "zigzag" grain motion will gradually drift inwards, then suddenly change to a unidirectional inward path when a critical radial distance is encountered (a sort of "Coulombic event horizon") at which the dipole strength supersedes the monopole strength --thus proving the presence of a dipole, while also quantifying the D/M ratio. In a second module, an homogeneous electric field eliminates dipole effects (both Coulombic and induced) to provide calibration of the monopole and to more readily evaluate net charge statistical variance. In both modules, the e-fields will be exponentially step-ramped in voltage during the experiment, so that the field "nominalizes" grain speed while spreading the response time --effectively forcing each grain to "wait its turn" to be measured. In addition to rigorously quantifying M, D, and the D/M ratio for many hundreds of grains, the experiment will also observe gross electrometric and RF discharge phenomena associated with grain activity. The parameter space will encompass grain charging levels (via intentional triboelectrification), grain size, cloud density, and material type. Results will prove or disprove the dipole hypothesis. In either case, light will be shed on the role of electrostatic forces in governing granular systems. Knowledge so gained can be applied to natural clouds such as protostellar and protoplanetary dust and debris systems, planetary rings, planetary dust palls and aerosols created by volcanic, impact, aeolian, firestorm, or nuclear winter processes. The data are also directly applicable to adhesion, cohesion, transport, dispersion, and collection of granular materials in industrial, agricultural, pharmaceutical applications, and in fields as diverse as dust contamination of space suits on Mars and crop spraying on Earth.

Description: This document presents in viewgraph form the opportunities, challenges, core capabilities and technologies involved in human deep space exploration.

Description: The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work have been limited. Here, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. Boiling involves both fluid flow and heat transfer and thus requires the solution of the Navier-Stokes and the energy equations. The numerical method is based on writing one set of governing transport equations which is valid in both the liquid and vapor phases. This local, single-field formulation incorporates the effect of the interface in the governing equations as source terms acting only at the interface. These sources account for surface tension and latent heat in the equations for conservation of momentum and energy as well as mass transfer across the interface due to phase change. The single-field formulation naturally incorporates the correct mass, momentum and energy balances across the interface. Integration of the conservation equations across the interface directly yields the jump conditions derived in the local instant formulation for two-phase systems. In the numerical implementation, the conservation equations for the whole computational domain (both vapor and liquid) are solved using a stationary grid and the phase boundary is followed by a moving unstructured two-dimensional grid. While two-dimensional simulations have been used for preliminary studies and to examine the resolution requirement, the focus is on fully three-dimensional simulations. The numerical methodology, including the parallelization and grid refinement strategy is discussed, and preliminary results shown. For buoyancy driven flow, the heat transfer is in good agreement with experimental correlations. The changes when gravity is turned off and/or fluid shear is added are discussed, as well as the difference between simulations of a layer freely releasing bubbles versus simulations using only one wavelength initial perturbation. Figure 1 shows the early stages of the formation of a three-dimensional bubble from a thin vapor layer. The boundary conditions are periodic in the x and y direction, the bottom is a hot and the top allows a free outflow. The jagged edge of the surface close to the bottom of the computational domain is due to some of the surface elements being on the other side of the domain and some elements not plotted by our plotting routine. In the second figure, we show the temperature distribution through two perpendicular planes.

Description: This talk will present recent results from ground-based research to support the "Physics of Colloids in Space" project which is scheduled to fly in the ISS approximately one year from now. In addition, results supporting future planned flights will be discussed.

Description: Experimental studies of growth and detachment processes of a single bubble and multiple bubbles formed on a heated surface have been conducted in the parabola flights of KC-135 aircraft. Distilled water and PF5060 were used as the test liquids. A micro-fabricated test surface was designed and built. Artificial cavities of diameters 10 microns, 7 microns and 4 microns were made on a thin polished Silicon wafer that was electrically heated by a number of small heating elements on the back side in order to control the surface superheat. Bubble growth period, bubble size and shape from nucleation to departure were measured under subcooled and saturation conditions. Significantly larger bubble departure diameters and bubble growth periods than those at earth normal gravity were observed. Bubble departure diameters as large as 20 mm for water and 6 mm for PF5060 were observed as opposed to about 3 mm for water and less than 1 mm for PF5060 at earth normal gravity respectively. It is found that the bubble departure diameter can be approximately related to the gravity level through the relation D(sub d) proportional 1/g(exp 1/2). For water,the effect of wall superheat and liquid subcooling on bubble departure diameter is found to be small.The growth periods are found to be very sensitive to liquid subcooling at a given wall superheat. However,the preliminary results of single bubble dynamics using PF5060 showed that the departure diameter increases when wall superheat is elevated at the same gravity and subcooling. Growth period of single bubbles in water has been found to vary as t(sub g) proportional g(exp -.93). For water, when the magnitude of horizontal gravitational components was comparable to that of gravity normal to the surface, single bubbles slid along the heater surface and departed with smaller diameter at the same gravity level in the direction normal to the surface. For PF5060, even a very small horizontal gravitational component caused the sliding of bubble along the surface. The numerical simulation has been carried out by solving under the condition of axisymmetry, the mass, momentum, and energy equations for the vapor and the liquid phases. In the model the contribution of micro-layer has been included and instantaneous shape of the evolving vapor-liquid interface is determined from the analysis. Consistent with the experimental results, it is found that effect of reduced gravity is to stretch the growth period and bubble diameter It is found that effect of reduced gravity is to stretch the growth period and bubble diameter at departure. The numerical simulations are in good agreement with the experimental data for both the departure diameters and the growth periods. In the study on dynamics of multiple bubbles, horizontal merger of 2,3 4,and 5 bubbles was observed. It is found that after merger of 2 and 3 bubbles the equivalent diameter of the detached bubble is smaller than that of a single bubble departing at the same gravity level. During and after bubble merger, liquid still fills the space between the vapor stems so as to form mushroom type bubbles. The experimental and numerical studies conducted so far have brought us a step closer to prediction of nucleate boiling heat fluxes under low gravity conditions. Preparations for a space flight are continuing.

Description: In the previous lecture, the issues related to arcjet flow modeling were introduced, and the limitations of conventional instrumentation in addressing these issues were discussed. The general level of understanding of the arcjet flows was seen to preclude the use of arcjets as aerothermodynamic test facilities beyond the current role in aerothermal material testing, despite their long test duration capability. In this section, the focus will be on new developments in spectroscopic instrumentation and techniques that can be brought to bear on the fundamental problem of arcjet stream characterization. Although a wide selection of arcjet facilities were introduced in the previous section, the discussion of nonintrusive diagnostic instrumentation will be restricted to the large-scale, segmented, constricted-arc heater facilities that are most widely used in thermal protection material testing for aerospace applications. After a brief review of the important features of arcjet flows, the topic of nonintrusive, optical diagnostics is introduced with a discussion of some of the basic aspects of radiative transitions. The lecture is then organized into two sections covering emission measurements and laser-induced fluorescence measurements. Emission measurements are presented next for different regions of arcjet flows, while the fluorescence measurements are presented for the free stream region only. Summaries are given for each of the two main sections, and observations on arcjet characterization by optical diagnostics in general are given at the end.

Description: Current efforts to develop new reusable launch vehicles and to pursue low-cost robotic planetary missions have led to a renewed interest in understanding arc-jet flows. Part of this renewed interest is concerned with improving the understanding of arc-jet test results and the potential use of available computational-fluid- dynamic (CFD) codes to aid in this effort. These CFD codes have been extensively developed and tested for application to nonequilibrium, hypersonic flow modeling. It is envisioned, perhaps naively, that the application of these CFD codes to the simulation of arc-jet flows would serve two purposes: first. the codes would help to characterize the nonequilibrium nature of the arc-jet flows; and second. arc-jet experiments could potentially be used to validate the flow models. These two objectives are, to some extent, mutually exclusive. However, the purpose of the present discussion is to address what role CFD codes can play in the current arc-jet flow characterization effort, and whether or not the simulation of arc-jet facility tests can be used to eva1uate some of the modeling that is used to formu1ate these codes. This presentation is organized into several sections. In the introductory section, the development of large-scale, constricted-arc test facilities within NASA is reviewed, and the current state of flow diagnostics using conventional instrumentation is summarized. The motivation for using CFD to simulate arc-jet flows is addressed in the next section, and the basic requirements for CFD models that would be used for these simulations are briefly discussed. This section is followed by a more detailed description of experimental measurements that are needed to initiate credible simulations and to evaluate their fidelity in the different flow regions of an arc-jet facility. Observations from a recent combined computational and experiment.al investigation of shock-layer flows in a large-scale arc-jet facility are then used to illustrate the current state of development of diagnostic instrumentation, CFD simulations, and general knowledge in the field of arc-jet characterization. Finally, the main points are summarized and recommendations for future efforts are given.

Description: We are developing particle methods oriented at improving computational modeling capabilities of multiscale physical phenomena in : (i) high Reynolds number unsteady vortical flows, (ii) particle laden and interfacial flows, (iii)molecular dynamics studies of nanoscale droplets and studies of the structure, functions, and evolution of the earliest living cell. The unifying computational approach involves particle methods implemented in parallel computer architectures. The inherent adaptivity, robustness and efficiency of particle methods makes them a multidisciplinary computational tool capable of bridging the gap of micro-scale and continuum flow simulations. Using efficient tree data structures, multipole expansion algorithms, and improved particle-grid interpolation, particle methods allow for simulations using millions of computational elements, making possible the resolution of a wide range of length and time scales of these important physical phenomena.The current challenges in these simulations are in : [i] the proper formulation of particle methods in the molecular and continuous level for the discretization of the governing equations [ii] the resolution of the wide range of time and length scales governing the phenomena under investigation. [iii] the minimization of numerical artifacts that may interfere with the physics of the systems under consideration. [iv] the parallelization of processes such as tree traversal and grid-particle interpolations We are conducting simulations using vortex methods, molecular dynamics and smooth particle hydrodynamics, exploiting their unifying concepts such as : the solution of the N-body problem in parallel computers, highly accurate particle-particle and grid-particle interpolations, parallel FFT's and the formulation of processes such as diffusion in the context of particle methods. This approach enables us to transcend among seemingly unrelated areas of research.

Description: The functionality of the overset, static aeroelasticity, Navier-Stokes flow solver OVERAERO was increased by adding capability to the flow solver and enhancing code performance. Improvements were made to the fluids/structure interface, an MLP version of the parallel OVERAERO code was developed, and the OVERAERO-MPI code was ported to the Cray T3E. The OVERFLOW-MPI and OVERAERO-MPI codes were tested successfully on the IPG testbed and a means of reducing communication overhead within OVERFLOW-MPI was investigated. To solve an aeroelastic problem computationally, a structures grid surface definition and a fluids grid surface definition are required. Typically, the structures grid surface has a lower fidelity than the fluids grid surface. Thus, the methods developed to transfer data between the two grid systems are vital to the accuracy and efficiency of the aeroelasticity code. The fluids/structures interface developed for the OVERAERO code was improved to more accurately treat fluids surfaces that bridge between two different structural surfaces. For example, the method allowed the forward portion of a flap track fairing to deform with the wing and the aft end of the fairing to deform with the flap. A tightly-coupled version of the code based on OVERFLOW-MLP was developed to improve code performance on the SGI Origin 2000. This required a new parallelization strategy to couple the fluids and structures codes. The OVERAERO-MPI code was ported to the Cray T3E to extend the usability of the code. The port required extensive use of dynamic memory management techniques to fit large problems within the memory limitations of the T3E. The OVERFLOW-MPI and OVERAERO-MPI codes were tested on the IPG testbed being developed within NASA. For small problems with minimal data transfer between grids, there was little to no performance penalty spreading the computation across two machines. For very large problems, methods were developed to minimize intermachine communication via the grid partitioning scheme. By minimizing the intermachine communication requirements of the problem, it may still be beneficial to run a tightly-coupled flow solver across two machines within the IPG.

Description: Based on work begun in the 1970s, Analytical Methods, Inc. (AMI) developed computational fluid dynamics (CFD) technology called VSAERO. VSAERO can image aerodynamic flows and has been used in a number of unique designs such as the U.S. sailboat used in the 1996 America's Cup, solar automobiles, and solar airplaces. VSAERO is applicable to fluid flow problems in aerospace, aeronautical, automotive, and marine engineering.

Description: One of the propulsion system concepts to be considered for the High-Speed Civil Transport (HSCT) is an underwing, dual-propulsion, pod-per-wing installation. Adverse transient phenomena such as engine compressor stall and inlet unstart could severely degrade the performance of one of these propulsion pods. The subsequent loss of thrust and increased drag could cause aircraft stability and control problems that could lead to a catastrophic accident if countermeasures are not in place to anticipate and control these detrimental transient events. Aircraft system engineers must understand what happens during an engine compressor stall and inlet unstart so that they can design effective control systems to avoid and/or alleviate the effects of a propulsion pod engine compressor stall and inlet unstart. The objective of the Inlet Unstart Propulsion Airframe Integration test program was to assess the underwing flow field of a High-Speed Civil Transport propulsion system during an engine compressor stall and subsequent inlet unstart. Experimental research testing was conducted in the 10- by 10-Foot Supersonic Wind Tunnel at the NASA Glenn Research Center at Lewis Field. The representative propulsion pod consisted of a two-dimensional, bifurcated inlet mated to a live turbojet engine. The propulsion pod was mounted below a large flat plate that acted as a wing simulator. Because of the plate s long length (nominally 10-ft wide by 18-ft long), realistic boundary layers could form at the inlet cowl plane. Transient instrumentation was used to document the aerodynamic flow-field conditions during an unstart sequence. Acquiring these data was a significant technical challenge because a typical unstart sequence disrupts the local flow field for about only 50 msec. Flow surface information was acquired via static pressure taps installed in the wing simulator, and intrusive pressure probes were used to acquire flow-field information. These data were extensively analyzed to determine the impact of the unstart transient on the surrounding flow field. This wind tunnel test program was a success, and for the first time, researchers acquired flow-field aerodynamic data during a supersonic propulsion system engine compressor stall and inlet unstart sequence. In addition to obtaining flow-field pressure data, Glenn researchers determined other properties such as the transient flow angle and Mach number. Data are still being reduced, and a comprehensive final report will be released during calendar year 2000.

Description: Under certain conditions, such as in thin liquid films or microgravity, surface tension variations along a free surface can induce convection. Convection onset due to surface tension variation is important to many terrestrial technological processes in addition to microgravity materials processing applications. Examples include coating, drying crystallization, solidification, liquid surface contamination, and containerless processing. In double-diffusive and multicomponent systems, the spatial variations of surface tension are associated with two or more stratifying agencies, respectively. For example, both temperature and species (concentration) gradients are associated with convection in the solidification of binary alloys or salt ponds. The direction of the two (or more) gradients has a profound effect on the nature of the flow at or slightly beyond the onset of convection. Our recent work at the NASA Lewis Research Center focused on characterizing surface-tension-induced onset of convection, often referred to as Marangoni-Benard convection. Exact solutions for the stationary neutral stability of multicomponent fluid layers with interfacial deformation were derived. These solutions also permit the computation of a boundary curve that separates the long and finite wavelength instabilities. Computing points along this boundary using the exact solution (when possible) is more efficient than the typical numerical approaches, such as finite difference or spectral methods. Above the curve, a long wavelength instability was predicted, suggesting that convection would occur principally through one large flow cell in the layer, whereas below the curve, finite wavelength instabilities occur which suggest multiple finite-sized circulation cells. For many common liquids with layer depths greater than 100 mm, finite wave instability is predicted under terrestrial conditions; however, with little exception, long wavelength instability is predicted in microgravity for the identical fluid systems.

Description: SHDOM is a general purpose, publicly available, three-dimensional atmospheric radiative transfer model. SHDOM is an explicit method, which means it solves for the whole radiation field, as distinct from Monte Carlo methods which solve for particular radiative outputs. SHDOM is particularly well suited for remote sensing applications, where it can compute outgoing radiances at many angles from a cloud field at virtually no extra cost. SHDOM is not appropriate for calculating domain average quantities for which Monte Carlo methods excel. The I3RC intercomparison offers an opportunity to explore the pros and cons of SHDOM and Monte Carlo models on some real world inhomogeneous cloud fields. Specifically, we wish to determine the computer resources required to achieve a particular accuracy for a certain number of outputs using SHDOM and Monte Carlo models. This will help guide modelers on the appropriate choice of SHDOM or Monte Carlo for their applications. To emphasize the importance of this accuracy versus CPU time tradeoff, we are submitting two SHDOM entries (low and high resolution) in the I3RC.

Description: The Students' Cloud Observations On-Line (S'COOL) Project involved students in K-16 as ground truth observers for a NASA Earth-Observing satellite instrument. The Clouds and Earth's Radiant Energy System (CERES) instrument allows scientists to study the Earth's energy budget and how clouds affect it. Student reports of cloud conditions help scientists verify their algorithms and allow students to be involved in obtaining and analyzing real scientific data. The presentation contains 23 slides.

Description: Researchers at the NASA Glenn Research Center at Lewis Field successfully tested a variable cowl lip inlet at simulated takeoff conditions in Glenn s 10- by 10-Foot Supersonic Wind Tunnel (10x10 SWT) as part of the High-Speed Research Program. The test was a follow-on to the Two-Dimensional Bifurcated (2DB) Inlet/Engine test. At the takeoff condition for a High-Speed Civil Transport aircraft, the inlet must provide adequate airflow to the engine with an acceptable distortion level and high-pressure recovery. The test was conducted to study the effectiveness of installing two rotating lips on the 2DB Inlet cowls to increase mass flow rate and eliminate or reduce boundary layer flow separation near the lips. Hardware was mounted vertically in the test section so that it extended through the tunnel ceiling and that the 2DB Inlet was exposed to the atmosphere above the test section. The tunnel was configured in the aerodynamic mode, and exhausters were used to pump down the tunnel to vacuum levels and to provide a maximum flow rate of approximately 58 lb/sec. The test determined the (1) maximum flow in the 2DB Inlet for each variable cowl lip, (2) distortion level and pressure recovery for each lip configuration, (3) boundary layer conditions near variable lips inside the 2DB Inlet, (4) effects of a wing structure adjacent to the 2DB Inlet, and (5) effects of different 2DB Inlet exit configurations. It also employed flow visualization to generate enough qualitative data on variable lips to optimize the variable lip concept. This test was a collaborative effort between the Boeing Company and Glenn. Extensive inhouse support at Glenn contributed significantly to the progress and accomplishment of this test.

Description: Two current techniques exist for the measurement of planar, three-component velocity fields. Both techniques require multiple views of the illumination plane in order to extract all three velocity components. Particle image velocimetry (PIV) is a high-resolution, high accuracy, planar velocimetry technique that provides valuable instantaneous velocity information in aeropropulsion test facilities. PIV can provide three-component flow-field measurements using a two-camera, stereo viewing configuration. Doppler global velocimetry (DGV) is another planar velocimetry technique that can provide three component flow-field measurements; however, it requires three detector systems that must be located at oblique angles from the measurement plane. The three-dimensional configurations of either technique require multiple (DGV) or at least large (stereo PIV) optical access ports in the facility in which the measurements are being conducted. Optical access is extremely limited in aeropropulsion test facilities. In many cases, only one optical access port is available. A hybrid measurement technique has been developed at the NASA Glenn Research Center, planar particle image and Doppler velocimetry (PPIDV), which combines elements from both the PIV and DGV techniques into a single detection system that can measure all three components of velocity across a planar region of a flow field through a single optical access port. In the standard PIV technique, a pulsed laser is used to illuminate the flow field at two closely spaced instances in time, which are recorded on a "frame-straddling" camera, yielding a pair of single-exposure image frames. The PIV camera is oriented perpendicular to the light sheet, and the processed PIV data yield the two-component velocity field in the plane of the light sheet. In the standard DGV technique, an injection-seeded Nd:YAG pulsed laser light sheet illuminates the seeded flow field, and three receiver systems are used to measure three components of velocity. The receiver systems are oriented at oblique angles to the light sheet in order to accurately resolve the three-component velocity. Each DGV receiver system contains two cameras, which share a common view of the illuminated flow through a beam-splitting cube. One camera views the illuminated flow directly (reference camera) and the second camera images the illuminated flow through an iodine vapor cell (signal camera). The laser frequency (wavelength) is adjusted so that the Doppler-shifted light from particles in the flow falls on an iodine absorption feature, see the following graph. The iodine vapor cell acts as a frequency-to-velocity filter by modulating the intensity of the transmitted light as a function of the flow velocity (Doppler shift). The ratio of the signal and reference images yields the component of the flow velocity along the bisector of the laser sheet propagation direction and the receiver system observation direction. The hybrid system employs a single-component DGV receiver system configured to simultaneously acquire PIV image data, as shown in the following diagram. The cameras used in the DGV receiver are replaced with PIV frame-straddling cameras, and the receiver system views the illuminated light sheet plane at 90 (as in the standard PIV configuration).

Description: A steady-state laser heat flux technique has been developed at the NASA Glenn Research Center at Lewis Field to obtain critical thermal conductivity data of ceramic thermal barrier coatings under the temperature and thermal gradients that are realistically expected to be encountered in advanced engine systems. In this study, thermal conductivity change kinetics of a plasma-sprayed, 254-mm-thick ZrO2-8 wt % Y2O3 ceramic coating were obtained at high temperatures. During the testing, the temperature gradients across the coating system were carefully measured by the surface and back pyrometers and an embedded miniature thermocouple in the substrate. The actual heat flux passing through the coating system was determined from the metal substrate temperature drop (measured by the embedded miniature thermocouple and the back pyrometer) combined with one-dimensional heat transfer models.

Description: This report contains the 2000 annual progress reports of the postdoctoral Fellows and visiting scholars of the Center for Turbulence Research (CTR). It summarizes the research efforts undertaken under the core CTR program. Last year, CTR sponsored sixteen resident Postdoctoral Fellows, nine Research Associates, and two Senior Research Fellows, hosted seven short term visitors, and supported four doctoral students. The Research Associates are supported by the Departments of Defense and Energy. The reports in this volume are divided into five groups. The first group largely consists of the new areas of interest at CTR. It includes efficient algorithms for molecular dynamics, stability in protoplanetary disks, and experimental and numerical applications of evolutionary optimization algorithms for jet flow control. The next group of reports is in experimental, theoretical, and numerical modeling efforts in turbulent combustion. As more challenging computations are attempted, the need for additional theoretical and experimental studies in combustion has emerged. A pacing item for computation of nonpremixed combustion is the prediction of extinction and re-ignition phenomena, which is currently being addressed at CTR. The third group of reports is in the development of accurate and efficient numerical methods, which has always been an important part of CTR's work. This is the tool development part of the program which supports our high fidelity numerical simulations in such areas as turbulence in complex geometries, hypersonics, and acoustics. The final two groups of reports are concerned with LES and RANS prediction methods. There has been significant progress in wall modeling for LES of high Reynolds number turbulence and in validation of the v(exp 2) - f model for industrial applications.

Description: A thin free shear layer containing an inflection point in the mean velocity profile is inherently unstable. Disturbances in the flow field can excite the unstable behavior of a shear layer, if the appropriate combination of frequencies and shear layer thicknesses exists, causing instability waves to grow. For other combinations of frequencies and thicknesses, these instability waves remain neutral in amplitude or decay in the downstream direction. A growing instability wave radiates noise when its phase velocity becomes supersonic relative to the ambient speed of sound. This occurs primarily when the mean jet flow velocity is supersonic. Thus, the small disturbances in the flow, which themselves may generate noise, have generated an additional noise source. It is the purpose of this problem to test the ability of CAA to compute this additional source of noise. The problem is idealized such that the exciting disturbance is a fixed known acoustic source pulsating at a single frequency. The source is placed inside of a 2D jet with parallel flow; hence, the shear layer thickness is constant. With the source amplitude small enough, the problem is governed by the following set of linear equations given in dimensional form.

Description: Because grid generation can consume 70 percent of the total analysis time for a typical three-dimensional viscous flow simulation for a practical engineering device, payoffs from research and development could reduce costs and increase throughputs considerably. In this study, researchers at the NASA Glenn Research Center at Lewis Field developed a new hybrid grid approach with the advantages of flexibility, high-quality grids suitable for an accurate resolution of viscous regions, and a low memory requirement. These advantages will, in turn, reduce analysis time and increase accuracy. They result from an innovative combination of structured and unstructured grids to represent the geometry and the computation domain. The present approach makes use of the respective strengths of both the structured and unstructured grid methods, while minimizing their weaknesses. First, the Chimera grid generates high-quality, mostly orthogonal meshes around individual components. This process is flexible and can be done easily. Normally, these individual grids are required overlap each other so that the solution on one grid can communicate with another. However, when this communication is carried out via a nonconservative interpolation procedure, a spurious solution can result. Current research is aimed at entirely eliminating this undesired interpolation by directly replacing arbitrary grid overlapping with a nonstructured grid called a DRAGON grid, which uses the same set of conservation laws over the entire region, thus ensuring conservation everywhere. The DRAGON grid is shown for a typical film-cooled turbine vane with 33 holes and 3 plenum compartments. There are structured grids around each geometrical entity and unstructured grids connecting them. In fiscal year 1999, Glenn researchers developed and tested the three-dimensional DRAGON grid-generation tools. A flow solver suitable for the DRAGON grid has been developed, and a series of validation tests are underway.

Description: Responses of the nightside magnetosphere and auroral zone to interplanetary shocks are studies using WIND solar wind data and POLAR UV imaging data.

Description: The Cassini spacecraft flew past the Earth in a trajectory almost along the Sun-Earth line, giving a unique perspective of low frequency waves in geospace.

Description: The turbulent latent and sensible heat fluxes are necessary to study heat budget of the upper ocean or initialize ocean general circulation models. In order to retrieve the latent heat flux from satellite observations authors mostly use a bulk approximation of the flux whose parameters are derived from different instrument. In this paper, an approach based on artificial neural networks is proposed and compared to the bulk method on a global data set and 3 local data sets.

Description: This paper describes the mission and system design for a viable, promising first mission application of solar sail technology.

Description: This paper summarizes the results of testing and gives a detailed report on the spectral, radiometric and spatial performance of the AIRS.

Description: The 1999-2000 Arctic stratospheric vortex was unusually cold, especially in the early winter lower stratosphere, with a larger area near polar stratospheric cloud formation temperatures in Dec and Jan, and much lower temperatures averaged over Nov-Jan, than any previously observed Arctic winter.

Description: The solar fields and particles environment and its interaction with planetary magnetospheres are not only of considerable scientific interest, they can impact human endeavors as well.

Description: This work shows quantitatively the pronounced differences between the density equilibration of very compressible dense fluids in Earth's gravity and those in microgravity. The work was performed onsite at the NASA Glenn Research Center at Lewis Field and is complete. Full details are given in references 1 and 2. Liquid-vapor critical fluids (e.g., water) at their critical temperature and pressure, are very compressible. They collapse under their own weight in Earth's gravity, allowing only a thin meniscus-like layer with the critical pressure to survive. This critical layer, however, greatly slows down the equilibration process of the entire sample. A complicating feature is the buoyancy-driven slow flows of layers of heavier and lighter fluid. This work highlights the incomplete understanding of the hydrodynamics involved in these fluids.