ALBERT

Description: This volume and its accompanying CD-ROM contain materials presented at the Minnowbrook III-2000 Workshop on Boundary Layer Transition and Unsteady Aspects of Turbomachinery Flows held at the Syracuse University Minnowbrook Conference Center, Blue Mountain Lake, New York, August 20-23, 2000. Workshop organizers were John E. LaGraff (Syracuse University), Terry V Jones (Oxford University), and J. Paul Gostelow (University of Leicester). The workshop followed the theme, venue, and informal format of two earlier workshops: Minnowbrook I (1993) and Minnowbrook II (1997). The workshop was focused on physical understanding the late stage (final breakdown) boundary layer transition, separation, and effects of unsteady wakes with the specific goal of contributing to engineering application of improving design codes for turbomachinery. The workshop participants included academic researchers from the USA and abroad, and representatives from the gas-turbine industry and government laboratories. The physical mechanisms discussed included turbulence disturbance environment in turbomachinery, flow instabilities, bypass and natural transition, turbulent spots and calmed regions, wake interactions with attached and separated boundary layers, turbulence and transition modeling and CFD, and DNS. This volume contains abstracts and copies of the viewgraphs presented, organized according to the workshop sessions. The viewgraphs are included on the CD-ROM only. The workshop summary and the plenary-discussion transcripts clearly highlight the need for continued vigorous research in the technologically important area of transition, separated and unsteady flows in turbomachines.

Description: Microgravity research at NASA has been an undertaking that has included both science and commercial approaches since the late 80s and early 90s. The Fluid Physics and Transport Phenomena community has been developed, through NASA's science grants, into a valuable base of expertise in microgravity science. This was achieved through both ground and flight scientific research. Commercial microgravity research has been primarily promoted thorough NASA sponsored Centers for Space Commercialization which develop cost sharing partnerships with industry. As an example, the Center for Advanced Microgravity Materials Processing (CAMMP)at Northeastern University has been working with cost sharing industry partners in developing Zeolites and zeo-type materials as an efficient storage medium for hydrogen fuel. Greater commercial interest is emerging. The U.S. Congress has passed the Commercial Space Act of 1998 to encourage the development of a commercial space industry in the United States. The Act has provisions for the commercialization of the International Space Station (ISS). Increased efforts have been made by NASA to enable industrial ventures on-board the ISS. A Web site has been established at http://commercial/nasa/gov which includes two important special announcements. One is an open request for entrepreneurial offers related to the commercial development and use of the ISS. The second is a price structure and schedule for U.S. resources and accommodations. The purpose of the presentation is to make the Fluid Physics and Transport Phenomena community, which understands the importance of microgravity experimentation, aware of important aspects of ISS commercial development. It is a desire that this awareness will be translated into a recognition of Fluid Physics and Transport Phenomena application opportunities coordinated through the broad contacts of this community with industry.

Description: The performance of a heat pipe system is greatly improved by the use of a dilute aqueous solution of about 0.0005 and about 0.005 moles per liter of a long chain alcohol as the working fluid. The surface tension-temperature gradient of the long-chain alcohol solutions turns positive as the temperature exceeds a certain value, for example about 40.degree. C. for n-heptanol solutions. Consequently, the Marangoni effect does not impede, but rather aids in bubble departure from the heating surface. Thus, the bubble size at departure is substantially reduced at higher frequencies and, therefore, increases the boiling limit of heat pipes. This feature is useful in microgravity conditions. In addition to microgravity applications, the heat pipe system may be used for commercial, residential and vehicular air conditioning systems, micro heat pipes for electronic devices, refrigeration and heat exchangers, and chemistry and cryogenics.

Description: The vision document provides an overview of the Climate Change Science Program (CCSP) long-term strategic plan to enhance scientific understanding of global climate change.This document is a companion to the comprehensive Strategic Plan for the Climate Change Science Program. The report responds to the Presidents direction that climate change research activities be accelerated to provide the best possible scientific information to support public discussion and decisionmaking on climate-related issues.The plan also responds to Section 104 of the Global Change Research Act of 1990, which mandates the development and periodic updating of a long-term national global change research plan coordinated through the National Science and Technology Council.This is the first comprehensive update of a strategic plan for U.S. global change and climate change research since the origal plan for the U.S. Global Change Research Program was adopted at the inception of the program in 1989.

Description: The purpose of this Science and Implementation Plan is to describe GAPP science objectives and the activities required to meet these objectives, both specifically for the near-term and more generally for the longer-term. The GEWEX Americas Prediction Project (GAPP) is part of the Global Energy and Water Cycle Experiment (GEWEX) initiative that is aimed at observing, understanding and modeling the hydrological cycle and energy fluxes at various time and spatial scales. The mission of GAPP is to demonstrate skill in predicting changes in water resources over intraseasonal-to-interannual time scales, as an integral part of the climate system.

Description: During its maiden voyage in May 1962, a Centaur upper stage rocket, mated to an Atlas booster, exploded 54 seconds after launch, engulfing the rocket in a huge fireball. Investigation revealed that Centaur's light, stainless-steel tank had split open, spilling its liquid-hydrogen fuel down its sides, where the flame of the rocket exhaust immediately ignited it. Coming less than a year after President Kennedy had made landing human beings on the Moon a national priority, the loss of Centaur was regarded as a serious setback for the National Aeronautics and Space Administration (NASA). During the failure investigation, Homer Newell, Director of Space Sciences, ruefully declared: "Taming liquid hydrogen to the point where expensive operational space missions can be committed to it has turned out to be more difficult than anyone supposed at the outset." After this failure, Centaur critics, led by Wernher von Braun, mounted a campaign to cancel the program. In addition to the unknowns associated with liquid hydrogen, he objected to the unusual design of Centaur. Like the Atlas rocket, Centaur depended on pressure to keep its paper-thin, stainless-steel shell from collapsing. It was literally inflated with its propellants like a football or balloon and needed no internal structure to give it added strength and stability. The so-called "pressure-stabilized structure" of Centaur, coupled with the light weight of its high- energy cryogenic propellants, made Centaur lighter and more powerful than upper stages that used conventional fuel. But, the critics argued, it would never become the reliable rocket that the United States needed.

Description: An integral, lightweight combustion chamber/nozzle assembly for a rocket engine has a refractory metal shell defining a chamber of generally frusto-conical contour. The shell communicates at its smaller end with a rocket body, and terminates at its larger end in a generally contact contour, which is open at its terminus and which serves as a nozzle for the rocket engine. The entire inner surface of the refractory metal shell has a thermal and oxidation barrier layer applied thereto. An ablative silica phenolic insert is bonded to the exposed surface of the thermal and oxidation barrier layer. The ablative phenolic insert provides a chosen inner contour for the combustion chamber and has a taper toward the open terminus of the nozzle. A process for fabricating the integral, lightweight combustion chamber/nozzle assembly is simple and efficient, and results in economy in respect of both resources and time.

Description: An injector for liquid fueled rocket engines wherein a generally flat core having a frustoconical dome attached to one side of the core to serve as a manifold for a first liquid, with the core having a generally circular configuration having an axis. The other side of the core has a plurality of concentric annular first slots and a plurality of annular concentric second slots alternating with the first slots, the second slots having a greater depth than said first slots. A bore extends through the core for inletting a second liquid into said core, the bore intersecting the second slots to feed the second liquid into the second slots. The core also has a plurality of first passageways leading from the manifold to the first annular slots for feeding the first liquid into said first slots. A faceplate brazed to said other side of the core is provided with apertures extending from the first and second slots through said face plate, these apertures being positioned to direct fuel and liquid oxygen into contact with each other in the combustion chamber. The first liquid may be liquid oxygen and the second liquid may be kerosene or liquid hydrogen.

Description: Much of NASA's Arctic Research is run through its Cryospheric Sciences Program. Arctic research efforts to date have focused primarily on investigations of the mass balance of the largest Arctic land-ice masses and the mechanisms that control it, interactions among sea ice, polar oceans, and the polar atmosphere, atmospheric processes in the polar regions, energy exchanges in the Arctic. All of these efforts have been focused on characterizing, understanding, and predicting, changes in the Arctic. NASA's unique vantage from space provides an important perspective for the study of these large scale processes, while detailed process information is obtained through targeted in situ field and airborne campaigns and models. An overview of NASA investigations in the Arctic will be presented demonstrating how the synthesis of space-based technology, and these complementary components have advanced our understanding of physical processes in the Arctic.

Description: This paper presents viewgraphs of physical sciences research priorities and plans at the Office of Biological and Physical Sciences Research (OBPR). The topics include: 1) Sixth Microgravity Fluid Physics and Transport Phenomena Conference; 2) Beneficial Characteristics of the Space Environment; 3) Windows of Opportunity for Research Derived from Microgravity; 4) Physical Sciences Research Program; 5) Fundamental Research: Space-based Results and Ground-based Applications; 6) Nonlinear Oscillations; and 7) Fundamental Research: Applications to Mission-Oriented Research.

Description: NASA's In-Space Propulsion (ISP) Program is designed to develop advanced propulsion technologies that can enable or greatly enhance near and mid-term NASA science missions by significantly reducing cost, mass, and/or travel times. These technologies include: Solar Electric Propulsion, Aerocapture, Solar Sails, Momentum Exchange Tethers, Plasma Sails and other technologies such as Advanced Chemical Propulsion. The ISP Program intends to develop cost-effective propulsion technologies that will provide a broad spectrum of mission possibilities, enabling NASA to send vehicles on longer, more useful voyages and in many cases to destinations that were previously unreachable using conventional means. The ISP approach to identifying and prioritizing these most promising technologies is to use mission and system analysis and subsequent peer review. The ISP program seeks to develop technologies under consideration to Technology Readiness Level (TRL) -6 for incorporation into mission planning within 3-5 years of initiation. The NASA TRL 6 represents a level where a technology is ready for system level demonstration in a relevant environment, usually a space environment. In addition, maximum use of open competition is encouraged to seek optimum solutions under ISP. Several NASA Research Announcements (NRA's) have been released asking industry, academia and other organizations to propose propulsion technologies designed to improve our ability to conduct scientific study of the outer planets and beyond. The ISP Program is managed by NASA Headquarters Office of Space Science and implemented by the Marshall Space Flight Center in Huntsville, Alabama.

Description: This paper presents viewgraphs of NASA's strategic and fundamental research program at the Office of Biological and Physical Research (OBPR). The topics include: 1) Colloid-Polymer Samples; 2) Pool Boiling Experiment; 3) The Dynamics of Miscible Interfaces: A Space Flight Experiment (MIDAS); and 4) ISS and Ground-based Facilities.

Description: This overview presents in viewgraph form, the NASA Program organization regarding fluid physics, physical sciences research in space and the connection to biology, the dual thrust of the fluid physics program, and the immediate and future plans of the physical science research division.

Description: Aeroassist technology development is a vital part of the NASA In-Space Propulsion Program (ISP), which is managed by the NASA Headquarters Office of Space Science, and implemented by the Marshall Space Flight Center in Huntsville, Alabama. Aeroassist is the general term given to various techniques to maneuver a space vehicle within an atmosphere, using aerodynamic forces in lieu of propulsive fuel. Within the ISP, the current aeroassist technology development focus is aerocapture. The objective of the ISP Aerocapture Technology Project (ATP) is to develop technologies that can enable and/or benefit NASA science missions by significantly reducing cost, mass, and/or travel times. To accomplish this objective, the ATP identifies and prioritizes the most promising technologies using systems analysis, technology advancement and peer review, coupled with NASA Headquarters Office of Space Science target requirements. Plans are focused on developing mid-Technology Readiness Level (TRL) technologies to TRL 6 (ready for technology demonstration in space).

Description: NASA conducted the Cirrus Regional Study of Tropical Anvils and Cirrus Layers (CRYSTAL) Florida Area Cirrus Experiment (FACE) during July 2002 for improved understanding of tropical cirrus. One of the goals was to improve the understanding of cirrus generation by convective updrafts. The reasons why some convective storms produce extensive cirrus anvils is only partially related to convective instability and the vertical transport ice mass by updrafts. Convective microphysics must also have an important role on cirrus generation, for example, there are hypotheses that homogeneous nucleation in convective updrafts is a major source of anvil ice particles. In this paper, we report on one intense CRYSTAL- FACE convective case on 16 July 2002 that produced extensive anvil.

Description: The capability for massively parallel programming (MPP) using a message passing interface (MPI) has been implemented into a three-dimensional version of the Goddard Cumulus Ensemble (GCE) model. The design for the MPP with MPI uses the concept of maintaining similar code structure between the whole domain as well as the portions after decomposition. Hence the model follows the same integration for single and multiple tasks (CPUs). Also, it provides for minimal changes to the original code, so it is easily modified and/or managed by the model developers and users who have little knowledge of MPP. The entire model domain could be sliced into one- or two-dimensional decomposition with a halo regime, which is overlaid on partial domains. The halo regime requires that no data be fetched across tasks during the computational stage, but it must be updated before the next computational stage through data exchange via MPI. For reproducible purposes, transposing data among tasks is required for spectral transform (Fast Fourier Transform, FFT), which is used in the anelastic version of the model for solving the pressure equation. The performance of the MPI-implemented codes (i.e., the compressible and anelastic versions) was tested on three different computing platforms. The major results are: 1) both versions have speedups of about 99% up to 256 tasks but not for 512 tasks; 2) the anelastic version has better speedup and efficiency because it requires more computations than that of the compressible version; 3) equal or approximately-equal numbers of slices between the x- and y- directions provide the fastest integration due to fewer data exchanges; and 4) one-dimensional slices in the x-direction result in the slowest integration due to the need for more memory relocation for computation.

Description: A detailed analysis of available in situ and remotely sensed N2O and CH4 data measured in the 1999/2000 winter Arctic vortex has been performed in order to quantify the temporal evolution of vortex descent. Differences in potential temperature (theta) among balloon and aircraft vertical profiles (an average of 19-23 K on a given N2O or CH4 isopleth) indicated significant vortex inhomogeneity in late fall as compared with late winter profiles. A composite fall vortex profile was constructed for 26 November 1999, whose error bars encompassed the observed variability. High-latitude extravortex profiles measured in different years and seasons revealed substantial variability in N2O and CH4 on theta surfaces, but all were clearly distinguishable from the first vortex profiles measured in late fall 1999. From these extravortex-vortex differences we inferred descent prior to 26 November: as much as 397 plus or minus 15 K (lsigma) at 30 ppbv N2O and 640 ppbv CH4, and falling to 28 plus or minus 13 K above 200 ppbv N2O and 1280 ppbv CH4. Changes in theta were determined on five N2O and CH4 isopleths from 26 November through 12 March, and descent rates were calculated on each N2O isopleth for several time intervals. The maximum descent rates were seen between 26 November and 27 January: 0.82 plus or minus 0.20 K/day averaged over 50- 250 ppbv N2O. By late winter (26 February to 12 March), the average rate had decreased to 0.10 plus or minus 0.25 K/day. Descent rates also decreased with increasing N2O; the winter average (26 November to 5 March) descent rate varied from 0.75 plus or minus 0.10 K/day at 50 ppbv to 0.40 plus or minus 0.11 K/day at 250 ppbv. Comparison of these results with observations and models of descent in prior years showed very good overall agreement. Two models of the 1999/2000 vortex descent, SLIMCAT and REPROBUS, despite theta offsets with respect to observed profiles of up to 20 K on most tracer isopleths, produced descent rates that agreed very favorably with the inferred rates from observation.

Description: The stability of cylindrical liquid bridges in reduced gravity is affected by ambient vibrations of the spacecraft. Such vibrations are expected to excite capillary modes of the bridge. The lowest-order unstable mode is particularly susceptible to vibration as the length of the bridge approaches the stability limit. This low-order mode is known as the (2,0) mode and is an axisymmetric varicose mode of one wavelength in the axial direction. In this work, an optical system is used to detect the (2,0)-mode amplitude. The derivative of the error signal produced by this detector is used to produce the appropriate voltages on a pair of ring electrodes which are concentric with the bridge. A mode-coupled Maxwell stress profile is thus generated in proportional to the modal velocity. Depending on the sign of the gain, the damping of the capillary oscillation can be either increased or decreased. This effect has been demonstrated in Plateau-tank experiments. Increasing the damping of the capillary modes on free liquid surfaces in space could be beneficial for containerless processing and other novel technologies.

Description: The Goddard Cumulus Ensemble (GCE) model is used to examine the sensitivities of surface fluxes, explicit radiation, and ice microphysical processes on multi-day simulations of deep tropical convection over the Tropical Ocean Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA COARE). The simulations incorporate large-scale advective temperature and moisture forcing, as well as large-scale momentum, that are updated every time step on a periodic lateral boundary grid. This study shows that when surface fluxes are eliminated, the mean atmosphere is much cooler and drier, convection and CAPE are much weaker, precipitation is less, and cloud coverage in stratiform regions much greater. Surface fluxes using the TOGA COARE flux algorithm are weaker than with the aerodynamic formulation, but closer to the observed fluxes. In addition, similar trends noted above for the case without surface fluxes are produced for the TOGA flux case, albeit to a much lesser extent. The elimination of explicit shortwave and longwave radiation is found to have only minimal effects on the mean thermodynamics, convection, and precipitation. However explicit radiation does have a significant impact on cloud temperatures and structure above 200 mb and on the overall mean vertical circulation. The removal of ice processes produces major changes in the structure of the cloud. Much of the liquid water is transported aloft and into anvils above the melting layer (600 mb), leaving narrow, but intense bands of rainfall in convective regions. The elimination of melting processes leads to greater hydrometeor mass below the melting layer, and produces a much warmer and moister boundary layer, leading to a greater mean CAPE. Finally, the elimination of the graupel species has only a small impact on mean total precipitation, thermodynamics, and dynamics of the simulation, but does produce much greater snow mass just above the melting layer. Some of these results differ from previous CRM studies of tropical systems, which is likely due to the type of simulated system, total time integration, and model setup.

Description: Observing system simulation experiments (OSSE's) provide an effective means to evaluate the potential impact of a proposed observing system, as well as to determine tradeoffs in their design, and to evaluate data assimilation methodology. Great care must be taken to ensure realism of the OSSE's, and in the interpretation of OSSE results. All of the OSSE's that have been conducted to date have demonstrated tremendous potential for space-based wind profile data to improve atmospheric analyses, forecasts, and research. This has been true for differing data assimilation systems, analysis methodology, and model resolutions. OSSE's clearly show much greater potential for observations of the complete wind profile than for single-level wind data or observations of the boundary layer alone.

Description: A mixing chamber used in rocket engine testing at the NASA Stennis Space Center is modelled by a system of two nonlinear ordinary differential equations. The mixer is used to condition the thermodynamic properties of cryogenic liquid propellant by controlled injection of the same substance in the gaseous phase. The three inputs of the mixer are the positions of the valves regulating the liquid and gas flows at the inlets, and the position of the exit valve regulating the flow of conditioned propellant. Mixer operation during a test requires the regulation of its internal pressure, exit mass flow, and exit temperature. A mathematical model is developed to facilitate subsequent controller designs. The model must be simple enough to lend itself to subsequent feedback controller design, yet its accuracy must be tested against real data. For this reason, the model includes function calls to thermodynamic property data. Some structural properties of the resulting model that pertain to controller design, such as uniqueness of the equilibrium point, feedback linearizability and local stability are shown to hold under conditions having direct physical interpretation. The existence of fixed valve positions that attain a desired operating condition is also shown. Validation of the model against real data is likewise provided.

Description: I will discuss recent experiments from my lab, which use surface templates to induce ordered colloidal structures. Particle assembly driven by entropic depletion, fluid convection, and sedimentation will be described. Confocal microscopy was used to visualize most of these samples.

Description: The widespread use of electro-hydrodynamic devices and processes emphasizes a critical need for developing a comprehensive predictive theory capable of improving our fundamental understanding of the behavior of a suspension subject to an AC electric field and shear, and of facilitating the design and optimization of such devices. The currently favored approach to the qualitative interpretation of the AC field driven manipulation of suspensions is based on a model which considers only the force exerted on a single particle by an external field and neglects the field-induced and hydrodynamic interparticle interactions both being inversely proportional to the interparticle distance raised to the power three. On the other hand, the purpose of the field-induced separation is to concentrate particles in certain regions of a device. This clearly raises the fundamental question regarding the extent to which we can neglect these slow decaying electrical and hydrodynamic collective interactions and rely on the predictions of a single-particle model. Another important issue that still remains open is how to characterize the polarization of a particle exposed to a strong electric field. The presentation will address both these questions. Experiments were conducted in a parallel-plate channel in which a 10(exp -3) (v/v) suspension of heavy, positively polarized Al2O3 spheres was exposed to an AC field under conditions such that the field lines were arranged in the channel cross-section perpendicular to the streamlines of the main flow. To reduce the effects of the gravitational settling of the particles, the channel was slowly rotated (4 rpm) around a horizontal axis. Following the application of a high-gradient strong AC field (approx. kV/mm), the particles were found to move towards both the high-voltage (HV) and grounded (GR) electrodes and to form 'bristles' along their edges.

Description: Spaceborne scatterometers are active microwave radar instruments designed to acquire near-simultaneous, spatially collocated measurements of the normalized radar backscattering cross section (sigma0) of the global surface from several azimuth and/or incidence angles. The primary objective of the scatterometer mission is to measure the near-surface wind speed and direction over the global ocean using sigma0 measurements together with a wind geophysical model function. However, since sigma0 measurements are collected globally all the time, sigma0 data can also be used for global land and ice applications. In this paper, we will first present the objectives of the QSCAT mission, the instrument design, and the unique features of the Ku-band scatterometer currently in operation, called SeaWinds on QuikSCAT (QSCAT). We will then present some emerging land and ocean applications of the QSCAT data, which include (1) global snow detection and monitoring, (2) melt region mapping on the Greenland ice sheet, (3) Monsoon flood detection and monitoring, (4) soil wetness application at large scale, and (5) hurricane monitoring and tracking.

Description: This paper studies the performance of a spaceborne precipitation radar in measuring vertical Doppler velocity of rainfall. As far as a downward pointing precipitation radar is concerned, one of the major problems affecting Doppler measurement at the nadir direction arises from the Non-Uniform Beam-Filling effect (NUBF). That is, when significant variation in rain rate is present within the radar IFOV (Instrument Field of View) in the along track direction. the Doppler shift caused by the radial component of the horizontal speed of the satellite is weighted differently among the portions of IFOV. The effects of this non-uniform weighting may dominate any other contribution. Under this condition, shape, average value and width of the Doppler spectrum may not be directly correlated with the vertical velocity of the precipitating particles. However, by using an inversion technique which over-samples the radar measurements in the along track direction, we show that the shift due to NUBF can be evaluated, and that the NUBF induced errors on average fall speed can be reduced.

Description: The conversion of cloud ice to snow by depositional growth, designated P(sub SFI), in the Goddard Cumulus Ensemble Model cloud physics parameterization is examined. The original formulation of P(sub SFI) is shown to produce excessive conversion of cloud ice to snow because of an implicit assumption that the relative humidity is 100% with respect to water even though the air may actually be quite less humid. Two possible corrections to this problem are proposed, the first involving application of a relative humidity dependent correction factor to the original formulation of P(sub SFI), and the second involving a new formulation of P(sub SFI) based on the equation for depositional growth of cloud ice.

Description: This project, supported by the NASA New Investigator Program, has primarily been funding a graduate student, Darren McKague. Since August 1999 Darren has been working part time at Raytheon, while continuing his PhD research. Darren is planning to finish his thesis work in May 2001, thus some of the work described here is ongoing. The proposed research was to use GOES visible and infrared imager data and SSM/I microwave data to obtain joint distributions of cirrus cloud ice mass and precipitation for a study region in the Eastern Tropical Pacific. These joint distributions of cirrus cloud and rainfall were to be compared to those from the CSU general circulation model to evaluate the cloud microphysical amd cumulus parameterizations in the GCM. Existing algorithms were to be used for the retrieval of cloud ice water path from GOES (Minnis) and rainfall from SSM/I (Wilheit). A theoretical study using radiative transfer models and realistic variations in cloud and precipitation profiles was to be used to estimate the retrieval errors. Due to the unavailability of the GOES satellite cloud retrieval algorithm from Dr. Minnis (a co-PI), there was a change in the approach and emphasis of the project. The new approach was to develop a completely new type of remote sensing algorithm - one to directly retrieve joint probability density functions (pdf's) of cloud properties from multi-dimensional histograms of satellite radiances. The usual approach is to retrieve individual pixels of variables (i.e. cloud optical depth), and then aggregate the information. Only statistical information is actually needed, however, and so a more direct method is desirable. We developed forward radiative transfer models for the SSM/I and GOES channels, originally for testing the retrieval algorithms. The visible and near infrared ice scattering information is obtained from geometric ray tracing of fractal ice crystals (Andreas Macke), while the mid-infrared and microwave scattering is computed with Mie scattering. The radiative transfer is performed with the Spherical Harmonic Discrete Ordinate Method (developed by the PI), and infrared molecular absorption is included with the correlated k-distribution method. The SHDOM radiances have been validated by comparison to version 2 of DISORT (the community "standard" discrete-ordinates radiative transfer model), however we use SHDOM since it is computationally more efficient.

Description: The microphysical parameterization of clouds and rain-cells plays a central role in atmospheric forward radiative transfer models used in calculating passive microwave brightness temperatures. The absorption and scattering properties of a hydrometeor-laden atmosphere are governed by particle phase, size distribution, aggregate density., shape, and dielectric constant. This study identifies the sensitivity of brightness temperatures with respect to the microphysical cloud parameterization. Cloud parameterizations for wideband (6-410 GHz observations of baseline brightness temperatures were studied for four evolutionary stages of an oceanic convective storm using a five-phase hydrometeor model in a planar-stratified scattering-based radiative transfer model. Five other microphysical cloud parameterizations were compared to the baseline calculations to evaluate brightness temperature sensitivity to gross changes in the hydrometeor size distributions and the ice-air-water ratios in the frozen or partly frozen phase. The comparison shows that, enlarging the rain drop size or adding water to the partly Frozen hydrometeor mix warms brightness temperatures by up to .55 K at 6 GHz. The cooling signature caused by ice scattering intensifies with increasing ice concentrations and at higher frequencies. An additional comparison to measured Convection and Moisture LA Experiment (CAMEX 3) brightness temperatures shows that in general all but, two parameterizations produce calculated T(sub B)'s that fall within the observed clear-air minima and maxima. The exceptions are for parameterizations that, enhance the scattering characteristics of frozen hydrometeors.

Description: A persistent, mesoscale region of intense eyewall convection contained within Hurricane Bonnie on 23 August 1998 is examined from multiple observations synthesized from the NASA ER-2 and DC-8 aircraft. The intense convection occurred late in the day as Bonnie was attaining its minimum central pressure and during a stage when the inner core featured a markedly asymmetric structure. The main purpose of this paper is to describe the internal structure of a convective burst and its relationship to the warm core using a synthesis of high-resolution satellite, aircraft radar, and in situ data. An exceptionally vigorous eyewall tower penetrating to nearly 18 km is described. A second intense eyewall tower, adjacent to the eye, is shown to be associated with a mesoscale subsiding current of air that extends horizontally nearly 20 km into the eye interior. The subsidence occupies a substantial depth within the eye and appears to be a much larger scale feature than the convectively-induced, symmetric overturning which commonly occurs on the upper-level flanks of convective towers in other tropical environments.

Description: Ensemble sets of simulation experiments were conducted with a single column model (SCM) using the Goddard GEOS II GCM physics containing a recent version of the Cumulus Scheme (McRAS) and a biosphere based land-fluxes scheme (SSiB). The study used the 18 July to 5 August 1995 ARM-CART (Atmospheric Radiation Measurement-Cloud Atmospheric Radiation Test-bed) data, which was collected at the ARM-CART site in the mid-western United States and analyzed for single column modeling (SCM) studies. The new findings affirm the earlier findings that the vegetation, which increases the solar energy absorption at the surface together with soil and soil-moisture dependent processes, which modulate the surface, fluxes (particularly evapotranspiration) together help to increase the local rainfall. In addition, the results also show that for the particular study period roughly 50% of the increased evaporation over the ARM-CART site would be converted into rainfall with the Column, while the remainder would be advected out to the large-scale. Notwithstanding the limitations of only one-way interaction (i.e., the large-scale influencing the regional physics and not vice versa), the current SCM simulations show a very robust relationship. The evaporation-precipitation relationship turns out to be independent of the soil types, and soil moisture; however, it is weakly dependent on the vegetation cover because of its surface-albedo effect. Clearly, these inferences are prone to weaknesses of the SCM physics, the assumptions of the large-scale being unaffected by gridscale (SCM-scale) changes in moist processes, and other limitations of the evaluation procedures.

Description: Multiple, new technologies for chemical systems are becoming available and include high temperature rockets, very light propellant tanks and structures, new bipropellant and monopropellant options, lower mass propellant control components, and zero boil off subsystems. Such technologies offer promise of increasing the performance of in-space chemical propulsion for energetic space missions. A mass model for pressure-fed, Earth and space-storable, advanced chemical propulsion systems (ACPS) was developed in support of the NASA MSFC In-Space Propulsion Program. Data from flight systems and studies defined baseline system architectures and subsystems and analyses were formulated for parametric scaling relationships for all ACPS subsystem. The paper will first provide summary descriptions of the approaches used for the systems and the subsystems and then present selected analyses to illustrate use of the model for missions with characteristics of current interest.

Description: A response is provided to Diem and colleagues discussion on the whether urbanization-enhanced precipitation should be maximized in the South-Southeast of Atlanta, Georgia as documented in Shepherd et al. (2003). Diem and colleagues have offered a critical response to Shepherd et al. (2002). The response herein offers both general and specific responses to issues raised by Diem and colleagues.

Description: Performance expectations of closed-Brayton-cycle heat exchangers to be used in 100-kWe nuclear space power systems were forecast. Proposed cycle state points for a system supporting a mission to three of Jupiter s moons required effectiveness values for the heat-source exchanger, recuperator and rejection exchanger (gas cooler) of 0.98,0.95 and 0.97, respectively. Performance parameters such as number of thermal units (Nm), equivalent thermal conductance (UA), and entropy generation numbers (Ns) varied from 11 to 19,23 to 39 kWK, and 0.019 to 0.023 for some standard heat exchanger configurations. Pressure-loss contributions to entropy generation were significant; the largest frictional contribution was 114% of the heat-transfer irreversibility. Using conventional recuperator designs, the 0.95 effectiveness proved difficult to achieve without exceeding other performance targets; a metallic, plate-fin counterflow solution called for 15% more mass and 33% higher pressure-loss than the target values. Two types of gas-coolers showed promise. Single-pass counterflow and multipass cross-counterflow arrangements both met the 0.97 effectiveness requirement. Potential reliability-related advantages of the cross-countefflow design were noted. Cycle modifications, enhanced heat transfer techniques and incorporation of advanced materials were suggested options to reduce system development risk. Carbon-carbon sheeting or foam proved an attractive option to improve overall performance.

Description: The paper presents an order-of-magnitude analysis of the physical processes occurring during the pulsations of a vapor bubble subject to a sound field and shows several numerical examples relating to vapor bubbles in water with and without a translational velocity relative to the liquid. Finally, the growth and collapse of a bubble in a small tube under the action of a heat pulse is considered and it is pointed out that, in suitable conditions, a potentially useful pumping effect without mechanical moving parts can be achieved.

Description: Measurements from the Halogen Occultation Experiment (HALOE) revealed the infrared signature of polar mesospheric clouds (PMCs), for the first time, HALOE PMC observations at eight wavelengths (2.45 - 10 microns) show remarkable agreement with model PMC spectra based on ice particle extinction, and thus provide the first confirmation that water ice is the primary component of PMCs. Because PMCs respond to changes in temperature and water vapor, they are considered an indicator of global climate change. We propose to further the understanding of PMCs using a decade of infrared measurements form HALOE. This effort will characterize PMC spectral properties, extinction profiles, and size distributions. Using this information, HALOE measurements will be used to make simultaneous retrievals of H2O3, and temperature, in the presence of PMCs. The simultaneous retrievals of particle properties, H2O3, and temperature will be used with HALOE NO data to provide a significant step forward in the knowledge of PMC characteristics and formation conditions. We will challenge fundamental theories of PMC formation, and investigate changes in PMC properties and related conditions over the length of the HALOE measurement record. HALOE has been operating without flaw since it was launched on October 11, 1991. Consequently, ten southern and ten northern PMC seasons have been observed thus far, providing a wealth of data for the study of PMC and related parameters.

Description: Heavy rainfall occurred over the western side of Taiwan's complex terrain from August 10 to 13, 1994 after Typhoon Doug moved northward from the East China Sea into Taiwan and on towards the Yellow Sea. On August 10, most of the rainfall fell over sloped areas. The heaviest daily rainfall totals were in excess of 200 mm over southwestern as well as central Taiwan. However, not much rainfall occurred over northern Taiwan. The lack of rainfall over northern Taiwan also occurred on August 11, 12 and 13. The larger rainfall amounts shifted westward from the sloped areas on August 10 toward lower terrain on August 11. On August 12 and 13, most of the higher rainfall amounts were found over the coastal area in southwestern Taiwan. Notably, about 300 to 400 mm per day fell over the coastal area in southwest Taiwan on August 12 and 13. The distribution of rainfall amount was different on August 10 and 11 (termed as Case 1) compared to August 12 and 13 (termed as Case 2). The environmental situation and precipitation characteristics are analyzed using EC/TOGA data, ground-based radar data, surface rainfall patterns, surface wind data, and upper air soundings. Chen at al. (2001) also categorized the precipitation pattern into two types, propagating and quasi-stationary. For the propagating type of precipitation, rainrates increased or remained the same as systems went from the plains to mountainous regions. With the quasi-stationary type of precipitation, however, rainrates decreased as precipitation propagated across the plains and into the mountains. The focus of this study is to understand what causes the h1aher amounts of rainfall over Taiwan, and what factors influence where the higher amounts of rainfall will occur, over sloped areas or over coastal areas.

Description: An observing system comprised of two lidars in geosynchronous orbit would enable the synoptic and meso-scale measurement of atmospheric winds and moisture, both of which are key first-order variables of the Earth's weather equation. Simultaneous measurement of these parameters at fast revisit rates promises large advancements in our weather prediction skills. Such capabilities would be unprecedented and yield greatly improved and finer resolution initial conditions for models, make existing costly and cumbersome measurement approaches obsolete, and obviate the use of numerical techniques needed to correct data obtained using present observing systems. Additionally, simultaneous synoptic wind and moisture observations would lead to improvements in model parameterizations, and in our knowledge of small-scale weather processes. Technology and science data product assessments are ongoing. Results will be presented during the conference.

Description: On a twelve-month voyage to Mars, one astronaut will require at least two tons of potable water and two tons of pure oxygen. Efficient, reliable fluid reclamation is therefore necessary for manned space exploration. Space habitats require a compact, flexible, and robust apparatus capable of solid-fluid mechanical separation over a wide range of fluid and particle densities and particle sizes. In space, centrifugal filtration, where particles suspended in fluid are captured by rotating fixed-fiber mat filters, is a logical candidate for mechanical separation. Non-colloidal particles are deposited on the fibers due to inertial impaction or direct interception. Since rotation rates are easily adjustable, inertial effects are the most practical way to control separation rates for a wide variety of multiphase mixtures in variable gravity environments. Understanding how fluid inertia and differential fluid-particle inertia, characterized by the Reynolds and Stokes numbers, respectively, affect deposition is critical in optimizing filtration in a microgravity environment. This work will develop non-intrusive optical diagnostic techniques for directly visualizing where and when non-colloidal particles deposit upon, or contact, solid surfaces: 'particle proximity sensors'. To model particle deposition upon a single filter fiber, these sensors will be used in ground-based experiments to study particle dynamics as in the vicinity of a large (compared with the particles) cylinder in a simply sheared (i.e., linearly-varying, zero-mean velocity profile) neutrally-buoyant, refractive-index matched solid-liquid suspension.

Description: The development of a satellite infrared technique for estimating convective and stratiform rainfall and its application in studying the diurnal variability of rainfall in Amazonia are presented. The Convective-Stratiform. Technique, calibrated by coincident, physically retrieved rain rates from the Tropical Rain Measuring Mission (TRMM) Microwave Imager (TMI), is applied during January to April 1999 over northern South America. The diurnal cycle of rainfall, as well as the division between convective and stratiform rainfall is presented. Results compare well (a one-hour lag) with the diurnal cycle derived from Tropical Ocean-Global Atmosphere (TOGA) radar-estimated rainfall in Rondonia. The satellite estimates reveal that the convective rain constitutes, in the mean, 24% of the rain area while accounting for 67% of the rain volume. The effects of geography (rivers, lakes, coasts) and topography on the diurnal cycle of convection are examined. In particular, the Amazon River, downstream of Manaus, is shown to both enhance early morning rainfall and inhibit afternoon convection. Monthly estimates from this technique, dubbed CST/TMI, are verified over a dense rain gage network in the state of Ceara, in northeast Brazil. The CST/TMI showed a high bias equal to +33% of the gage mean, indicating that possibly the TMI estimates alone are also high. The root mean square difference (after removal of the bias) equaled 36.6% of the gage mean. The correlation coefficient was 0.77 based on 72 station-months.

Description: Global rainfall is the primary redistributor of earth's energy by the process of latent heat release. This forms the main driving force for the tropical circulation, which in turn impacts the global circulation .through transient events-such as El Nino. Hence, more precise and long-term time series of the rainfall and its variability is crucial to the understanding and prediction of the global climate and climate change. The Precipitation Radar (PR) abroad the US/Japan Tropical Rainfall Measuring Mission (TRMM) is the first radar ever launched into space that measures detailed vertical profiles of rain intensity over the tropics. One of the challenges in estimating rainfall from spaceborne radars is the presence of attenuation at frequencies, such as 14 GHz of the TRMM PR and future planned systems at this and higher frequencies. A common approach in current rainfall retrieval algorithms is to employed the path integrated attenuation (PIA) as a constraint to the retrieval, and hence overcome errors in the radar calibration or in the assumed rainfall parameters. PIA can either be derived from a radiometer or from the surface reference technique, in which a clear air radar measurement is compared with the measurement in the raining area. The current TRMM 2A21 PIA data product makes use of both a temporal and spatial clear air database for comparison to rainy measurements. In this paper we present results from analysis of TRMM surface backscatter cross-section (sigmaO) measurements from Nov 97-Feb99, and a comparison with sigmaO measurements obtained by the NASA Scatterometer (NSCAT) between Sept96-June97. Measurements for a given month from both instruments are compiled on a 1 deg. (lat.) x 1 deg. (lon.) x 1 hr. grid. This enables TRMM--NSCAT comparison and the investigation of seasonal and diurnal trends in both data sets. From preliminary analysis of TRMM sigmaO's we have decided not to treat the ocean as a single homogeneous region but to select a number of ocean sub-regions and individually analyze their trends. Likewise, and in a similar approach to previous studies of Seasat over-land data, we have selected a number of over-land regions for study. From said sigmaO maps and regional trend analysis we investigate possible sources of trends and variability. In addition, we study the effects of TRMM PR sensitivity through the PR "possible rain" class. Given NSCAT's inability to flag rain contaminated measurements we are able to gauge the impact of rain contamination on NSCAT monthly sigmaO maps, using TRMM measurements. The research described in this paper was carried out by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration (NASA), U.S.A.

Description: The NASA ER-2 and DC-8 aircraft collected remote sensing and in situ data sets from Hurricane Bonnie (23, 24, and 26 August 1998) during the Convection And Moisture Experimental-3 (CAMEX-3). Bonnie was an exceptional case where NASA and NOAA had five aircraft sampling both upper levels and lower altitudes. The ER-2 was instrumented with the ER-2 Doppler XBand radar (EDOP) and several radiometers ranging from visible to lower frequency microwaves. EDOP is a fixed dual-beam radar (nadir and forward-looking beams) which allows computation of both vertical and alongtrack horizontal winds. The hurricane secondary circulation is typically difficult to measure at upper levels due to aircraft altitude limitations and sensitivity of the lower altitude airborne radars. EDOP is in principle, well suited to measure these components of the wind. When ER-2 flies across the approximate center of the hurricane circulation, the along-track winds derived from EDOP, are approximately equal to the hurricane radial flow comprising the secondary circulation. Assuming that the hydrometeor fallspeeds can be approximated, the radial and vertical wind components of the secondary circulation can be measured. Since the hydrometeor motions can be estimated with more confidence in the higher altitude ice regions (i.e., graupel and mixed phase are complicated at lower altitudes), the derived radial and vertical winds have higher accuracy at upper levels. On the other hand, the reflectivities are extremely low at higher altitudes, resulting in fewer Doppler velocity estimates.

Description: NASA's In-Space Propulsion Technology Program is investing in technologies that have the potential to revolutionize the robotic exploration of deep space. For robotic exploration and science missions, increased efficiencies of future propulsion systems are critical to reduce overall life-cycle costs and, in some cases, enable missions previously considered impossible. Continued reliance on conventional chemical propulsion alone will not enable the robust exploration of deep space - the maximum theoretical efficiencies have almost been reached and they are insufficient to meet needs for many ambitious science missions currently being considered. The In-Space Propulsion Technology Program's technology portfolio includes many advanced propulsion systems. From the next-generation ion propulsion system operating in the 5- to 10-kW range to aerocapture and solar sails, substantial advances in - spacecraft propulsion performance are anticipated. Some of the most promising technologies for achieving these goals use the environment of space itself for energy and propulsion and are generically called 'propellantless' because they do not require onboard fuel to achieve thrust. Propellantless propulsion technologies include scientific innovations such as solar sails, electrodynamic and momentum transfer.tethers, aeroassist and aerocapture. This paper will provide an overview of both propellantless and propellant-based advanced propulsion technologies, as well as NASA's plans for advancing them as part of the In-Space Propulsion Technology Program.

Description: The onset of the southeast Asian monsoon during 1997 and 1998 was simulated by coupling a mesoscale atmospheric model (MM5) and a detailed, land surface model, PLACE (the Parameterization for Land-Atmosphere-Cloud Exchange). The rainfall results from the simulations were compared with observed satellite data from the TRMM (Tropical Rainfall Measuring Mission) TMI (TRMM Microwave Imager) and GPCP (Global Precipitation Climatology Project). The control simulation with the PLACE land surface model and variable sea surface temperature captured the basic signatures of the monsoon onset processes and associated rainfall statistics. Sensitivity tests indicated that simulations were significantly improved by including the PLACE land surface model. The mechanism by which the land surface processes affect the moisture transport and the convection during the onset of the southeast Asian monsoon were analyzed. The results indicated that land surface processes played an important role in modifying the low-level wind field over two major branches of the circulation: the southwest low-level flow over the Indo-China peninsula and the northern, cold frontal intrusion from southern China. The surface sensible and latent heat fluxes modified the low-level temperature distribution and merit, and therefore the low-level wind due to the thermal wind effect. The more realistic forcing of the sensible and latent heat fluxes from the detailed, land surface model improved the low-level wind simulation and associated moisture transport and convection.

Description: A cloud-resolving model is used to study sensitivities of two different microphysical schemes, one is the bulk type, and the other is an explicit bin scheme, in simulating a mid-latitude squall line case (PRE-STORM, June 10-11, 1985). Simulations using different microphysical schemes are compared with each other and also with the observations. Both the bulk and bin models reproduce the general features during the developing and mature stage of the system. The leading convective zone, the trailing stratiform region, the horizontal wind flow patterns, pressure perturbation associated with the storm dynamics, and the cool pool in front of the system all agree well with the observations. Both the observations and the bulk scheme simulation serve as validations for the newly incorporated bin scheme. However, it is also shown that, the bulk and bin simulations have distinct differences, most notably in the stratiform region. Weak convective cells exist in the stratiform region in the bulk simulation, but not in the bin simulation. These weak convective cells in the stratiform region are remnants of the previous stronger convections at the leading edge of the system. The bin simulation, on the other hand, has a horizontally homogeneous stratiform cloud structure, which agrees better with the observations. Preliminary examinations of the downdraft core strength, the potential temperature perturbation, and the evaporative cooling rate show that the differences between the bulk and bin models are due mainly to the stronger low-level evaporative cooling in convective zone simulated in the bulk model. Further quantitative analysis and sensitivity tests for this case using both the bulk and bin models will be presented in a companion paper.

Description: Despite the obvious notion that the presence of hail or graupel is a good indication of convection, the model results show this does not provide an objective benchmark partly due to the unrealistic presence of small amounts of hail or graupel throughout the anvil in the model but mainly because of the significant amounts of hail or graupel, especially in the tropical TOGA COARE simulation, in the transition zone. Without use of a "transition" category, it is open to debate as how this region should best be defined, as stratiform or as convective. So, the presence of significant hail or graupel contents in this zone significantly degrades its use an objective benchmark for convection. The separation algorithm comparison was done in the context of a cloud-resolving model. These models are widely used and serve a variety of purposes especially with regard to retrieving information that cannot be directly measured by providing synthetic data sets that are consistent and complete. Separation algorithms are regularly applied in these models. However, as with any modeling system, these types 'of models are constantly being improved to overcome any known deficiencies and make them more accurate representations of observed systems. The presence of hail and graupel in the anvil and the bias towards heavy rainfall rates are two such examples of areas that need improvement. Since, both of these can effect the perceived performance of the separation algorithms, the Lang et al. (2003) study did not want to overstate the relative performance of any specific algorithms.

Description: In this study, we investigate the seasonal variations in surface rainfall and associated large-scale processes in the tropical eastern Atlantic and West African region. The 5-yr (1998-2002) high-quality TRMM rainfall, sea surface temperature (SST), water vapor and cloud liquid water observations are applied along with the NCEP/NCAR reanalysis wind components and a 3-yr (2000-2002) Quickscat satellite-observed surface wind product. Major mean rainfall over West Africa tends to be concentrated in two regions and is observed in two different seasons, manifesting an abrupt shift of the mean rainfall zone during June-July. (i) Near the Gulf of Guinea (about 5 degN), intense convection and rainfall are seen during April-June and roughly follow the seasonality of SST in the tropical eastern Atlantic. (ii) Along the latitudes of about 10 deg. N over the interior West African continent, a second intense rain belt begins to develop from July and remains there during the later summer season. This belt co-exists with a northwardmoved African Easterly Jet (AEJ) and its accompanying horizonal and vertical shear zones, the appearance and intensification of an upper tropospheric Tropical Easterly Jet (TEJ), and a strong low-level westerly flow. Westward-propagating wave signals [ i e . , African easterly waves (AEWs)] dominate the synoptic-scale variability during July-September, in contrast to the evident eastward-propagating wave signals during May- June. The abrupt shift of mean rainfall zone thus turns out to be a combination of two different physical processes: (i) Evident seasonal cycles in the tropical eastern Atlantic ocean which modulate convection and rainfall in the Gulf of Guinea by means of SST thermal forcing and SST-related meridional gradient; (ii) The interaction among the AEJ, TEJ, low-level westerly flow, moist convection and AEWs during July-September which modulates rainfall variability in the interior West Africa, primarily within the ITCZ rain band. Evident seasonality in synoptic-scale wave signals is shown to be a good evidence for this seasonal evolution.

Description: Onboard radioisotope power systems being developed to support future NASA exploration missions require reliable design lifetimes of up to 14 yr and beyond. The structurally critical heater head of the high-efficiency developmental Stirling power converter has undergone extensive computational analysis of operating temperatures (up to 650 C), stresses, and creep resistance of the thin-walled Inconel 718 bill of material. Additionally assessment of the effect of uncertainties in the creep behavior of the thin-walled heater head, the variation in the manufactured thickness, variation in control temperature, and variation in pressure on the durability and reliability were performed. However, it is possible for the heater head to experience rare incidences of random temperature spikes (excursions) of short duration. These incidences could occur randomly with random magnitude and duration during the desired mission life. These rare incidences could affect the creep strain rate and therefore the life. The paper accounts for these uncertainties and includes the effect of such rare incidences, random in nature, on the reliability. The sensitivities of variables affecting the reliability are quantified and guidelines developed to improve the reliability are outlined. Furthermore, the quantified reliability is being verified with test data from the accelerated benchmark tests being conducted at the NASA Glenn Research Center.

Description: We seek funding from NASA for the third year (2005) of the four-year period January 1, 2003 - December 31, 2006 for continued support of the MIT contributions to the multi-national global atmospheric trace species measurement program entitled Advanced Global Atmospheric Gases Experiment (AGAGE). The case for real-time high-frequency measurement networks like AGAGE is very strong and the observations and their interpretation are widely recognized for their importance to ozone depletion and climate change studies and to verification issues arising from the Montreal Protocol (ozone) and Kyoto Protocol (climate). The proposed AGAGE program is distinguished by its capability to measure over the globe at high frequency almost all of the important species in the Montreal Protocol and almost all of the significant non-CO2 gases in the Kyoto Protocol.

Description: Stirling power conversion is being considered for use in a Radioisotope Power System for deep-space science missions because it offers a multifold increase in the conversion efficiency of heat to electric power. Quantifying the reliability of a Radioisotope Power System that utilizes Stirling power conversion technology is important in developing and demonstrating the capability for long-term success. A description of the Stirling power convertor is provided, along with a discussion about some of the key components. Ongoing efforts to understand component life, design variables at the component and system levels, related sources, and the nature of uncertainties is discussed. The requirement for reliability also is discussed, and some of the critical areas of concern are identified. A section on the objectives of the performance model development and a computation of reliability is included to highlight the goals of this effort. Also, a viable physics-based reliability plan to model the design-level variable uncertainties at the component and system levels is outlined, and potential benefits are elucidated. The plan involves the interaction of different disciplines, maintaining the physical and probabilistic correlations at all the levels, and a verification process based on rational short-term tests. In addition, both top-down and bottom-up coherency were maintained to follow the physics-based design process and mission requirements. The outlined reliability assessment approach provides guidelines to improve the design and identifies governing variables to achieve high reliability in the Stirling Radioisotope Generator design.

Description: Computational phenomena (i.e., spurious supersaturation and negative mixing ratio of cloud water) usually exist in cloud-resolving models when the time step for explicit integration is larger than a microphysical timescale in clouds. In this paper, the microphysical timescales in clouds are studied, showing that the timescale of water vapor condensation (or cloud water evaporation) is smaller than 10 s - the order of a typical time step for cloud-resolving models. To avoid spurious computational phenomena in cloud-resolving modeling, it is suggested that moist entropy be used as a prognostic thermodynamic variable, and temperature be diagnosed from that and other prognostic variables. A simple numerical model with moist entropy as a prognostic variable, for example, is presented to show that spurious computational phenomena are removed when moist entropy is used as a prognostic variable.

Description: Recently Grabowski and Khairoutdinov and Randall have proposed and demonstrated, the feasibility of using 2D CRMs (these did not include detailed land processes and used periodic lateral boundary conditions) as 'super-parameterizations or a multi-scale modeling framework MMF' for cloud processes within atmospheric generation circulation models (GCMs). In the MMF, a fine-resolution 2D CRM takes the place of the single-column parameterization that is used in conventional GCMs. An overview of this approach is given by Randall et al.. Since a CRM can explicitly simulate cloud processes at the natural space and time scales of clouddynamical processes, it's possible to compute statistics of cloud fields, including radiances and radar backscatters, directly simulated/produced by the CRM. A seed fund is available at NASA Goddard to build a MMF based on a 2D GCE model and the Goddard fvGCM. A prototype MMF will be developed at end of 2004. The major objectives are to test and evaluate its performance using the GCE model as the super parameterization and compare its results with other super parameterizations, conventional GCMs and observations. Recently, a detailed spectral-bin microphysical scheme was implemented into the Goddard Cumulus Ensemble (WE) model. The GCE model is a CRM and it has been improved and applied to study the precipitation processes in the past two decades. Atmospheric aerosols are also described using number density size-distribution functions. The model is tested by studying the evolution of deep tropical clouds in the west Pacific warm pool region and summertime convection over a mid-latitude continent with different concentrations of CCN: a low 'clean' concentration and a high 'dirty' concentration. The impact of atmospheric aerosol concentration on cloud and precipitation will be investigated. In addition, the similarities and differences between bulk- and spectral-bin microphysics will be presented,

Description: The Electric Particulate Suspension is a fire safety ignition test system being developed at Iowa State University with NASA support for evaluating combustion properties of powders, powder-gas mixtures, and pure gases in microgravity and gravitational atmospheres (quenching distance, ignition energy, flammability limits). A separate application is the use of EPS technology to control heat transfer in vacuum and space environment enclosures. In combustion testing, ignitable powders (aluminum, magnesium) are introduced in the EPS test cell and ignited by spark, while the addition of inert particles act as quenching media. As a combustion research tool, the EPS method has potential as a benchmark design for quenching powder flames that would provide NASA with a new fire safety standard for powder ignition testing. The EPS method also supports combustion modeling by providing accurate measurement of flame-quenching distance as an important parameter in laminar flame theory since it is closely related to characteristic flame thickness and flame structure. In heat transfer applications, inert powder suspensions (copper, steel) driven by electric fields regulate heat flow between adjacent surfaces enclosures both in vacuum (or gas) and microgravity. This simple E-field control can be particularly useful in space environments where physical separation is a requirement between heat exchange surfaces.

Description: Fast amplitude variations due to atmospheric scintillation are the main concerns for the Deep Space Network (DSN) Ka-band downlink under clear weather conditions. A theoretical study of the amplitude scintillation variances for a finite aperture antenna is presented. Amplitude variances for weak scattering scenarios are examined using turbulence theory to describe atmospheric irregularities. We first apply the Kolmogorov turbulent spectrum to a point receiver for three different turbulent profile models, especially for an exponential model varying with altitude. These analytic solutions then are extended to a receiver with a finite aperture antenna for the three profile models. Smoothing effects of antenna aperture are expressed by gain factors. A group of scaling factor relations is derived to show the dependences of amplitude variances on signal wavelength, antenna size, and elevation angle. Finally, we use these analytic solutions to estimate the scintillation intensity for a DSN Goldstone 34-m receiving station. We find that the (rms) amplitude fluctuation is 0.13 dB at 20-deg elevation angle for an exponential model, while the fluctuation is 0.05 dB at 90 deg. These results will aid us in telecommunication system design and signal-fading prediction. They also provide a theoretical basis for further comparison with other measurements at Ka-band.

Description: A key aspect of this project is the establishment of a causal link from circulation anomalies in the lower mesosphere and stratopause region downward through the stratosphere to the troposphere. The observational link for stratospheric sudden warmings and surface climate is fairly clear. However, our understanding of the dynamics is incomplete. We have been making significant progress in the area of dynamical mechanisms by which circulation anomalies in the stratosphere affect the troposphere. We are trying to understand the details and sequence of events that occur when a middle atmosphere (wind) anomaly propagates downward to near the tropopause. The wind anomaly could be caused by a warming or solar variations in the low-latitude stratopause region, or could have other causes. The observations show a picture that is consistent with a circulation anomaly that descends to the tropopause region, and can be detected as low as the mid-troposphere. Processes near the stratopause in the tropics appear to be important precursors to the wintertime development of the northern polar vortex. This may affect significantly our understanding of the process by which low-latitude wind anomalies in the low mesosphere and upper stratosphere evolve through the winter and affect the polar vortex.

Description: The stability of cylindrical liquid bridges in reduced gravity is affected by ambient vibrations of the spacecraft. Such vibrations are expected to excite capillary modes of the bridge. The lowest-order unstable mode is particularly susceptible to vibration as the length of the bridge approaches the stability limit. This low-order mode is known as the (2,0) mode and is an axisymmetric varicose mode of one wavelength in the axial direction. In this work, an optical system is used to detect the (2,0)-mode amplitude. The derivative of the error signal produced by this detector is used to produce the appropriate voltages on a pair of ring electrodes which are concentric with the bridge. A mode-coupled Maxwell stress profile is thus generated in proportional to the modal velocity. Depending on the sign of the gain, the damping of the capillary oscillation can be either increased or decreased. This effect has been demonstrated in Plateau-tank experiments. Increasing the damping of the capillary modes on free liquid surfaces in space could be beneficial for containerless processing and other novel technologies. [work supported by NASA]

Description: Spray cooling has high potential in thermal management and life support systems by overcoming the deleterious effect of microgravity upon two-phase heat transfer. In particular spray cooling offers several advantages in heat flux removal that include the following: 1. By maintaining a wetted surface, spray droplets impinge upon a thin fluid film rather than a dry solid surface 2. Most heat transfer surfaces will not be smooth but rough. Roughness can enhance conductive cooling, aid liquid removal by flow channeling. 3. Spray momentum can be used to a) substitute for gravity delivering fluid to the surface, b) prevent local dryout and potential thermal runaway and c) facilitate liquid and vapor removal. Yet high momentum results in high We and Re numbers characterizing the individual spray droplets. Beyond an impingement threshold, droplets splash rather than spread. Heat flux declines and spray cooling efficiency can markedly decrease. Accordingly we are investigating droplet impingement upon a) dry solid surfaces, b) fluid films, c) rough surfaces and determining splashing thresholds and relationships for both dry surfaces and those covered by fluid films. We are presently developing engineering correlations delineating the boundary between splashing and non-splashing regions.

Description: The Pennsylvania State University-National Center for Atmospheric Research mesoscale model MM5 is used to simulate Hurricane Bonnie at high resolution (2-km spacing) in order to examine budgets of water vapor, cloud condensate, and precipitation. Virtually all budget terms are derived directly from the model (except for the effects of storm motion). The water vapor budget reveals that a majority of the condensation in the eyewall occurs in convective hot towers, while outside of the eyewall most of the condensation occurs in weak updrafts, indicative of a dominance of stratiform precipitation processes. The ocean source of water vapor in the eyewall region is only a very small fraction of that transported inward in the boundary layer inflow or that condensed in the updrafts. In contrast, in the outer regions, the ocean vapor source is larger owing to the larger area, counters the drying effect of low-level subsidence, and enhances the moisture transported in toward the eyewall. In this mature storm, cloud condensate is consumed as rapidly as it is produced. Cloud water peaks at the top of the boundary layer and within the melting layer, where. cooling from melting enhances condensation. Unlike squall lines, in the hurricane, very little condensate produced in the eyewall convection is transported outward into the surrounding precipitation area. Most of the mass ejected outward is likely in the form of small snow particles that seed the outer regions and enhance stratiform precipitation development through additional growth by vapor deposition and aggregation. In addition, artificial source terms for cloud and precipitation mass associated with setting negative mixing ratios, produced as a result of numerical advection errors, to zero are described. Although small at any given point and time, the cumulative effect of these terms contributes an amount of mass equivalent to 13% of the total condensation. Thus, these terms must be accounted for to balance the model budgets and the results suggest the need for improved model numerics.

Description: The present work analyses the dynamics of a suspension of heavy particles in shear flow. The magnitude of the particle inertia is given by the Stokes number St = m(gamma/6(pi)a, which is the ratio of the viscous relaxation time of a particle tau(sub p) = m=6pi(eta)a to the flow time gamma(sup -1). Here, m is the mass of the particle, a is its size, eta is the viscosity of the suspending fluid and gamma is the shear rate. The ratio of the Stokes number to the Reynolds number, Re = (rho)f(gamma)a(exp 2)/eta, is the density ratio rho(sub p)/rho(sub f). Of interest is to understand the separate roles of particle (St) and fluid (Re) inertia in the dynamics of suspensions. In this study we focus on heavy particles, rho(sub p)/rho(sub f) much greater than 1, for which the Stokes number is finite, but the Reynolds number is sufficiently small for inertial forces in the fluid to be neglected; thus, the fluid motion is governed by the Stokes equations. On the other hand, the probability density governing the statistics of the suspended particles satisfies a Fokker-Planck equation that accounts for both configuration and momentum coordinates, the latter being essential for finite St. The solution of the Fokker-Planck equation is obtained to O(St) via a Chapman-Enskog type-procedure, and the conditional velocity distribution so obtained is used to derive a configuration-space Smoluchowski equation with inertial corrections. The inertial effects are responsible for asymmetry in the relative trajectories of two spheres in shear flow, in contrast to the well known symmetric structure in the absence of inertia. Finite St open trajectories in the plane of shear suffer a downward lateral displacement resulting from the inability of a particle of finite mass to follow the curvature of the zero-Stokes-number pathlines. In addition to the induced asymmetry, the O(St) inertial perturbation dramatically alters the nature of the near-field trajectories. The stable closed orbits (for St = 0) in the plane of shear now spiral in, approaching particle-particle contact in the limit. All trajectories starting from an initial offset of O(St(sup 1/2) or less (which remain open for St = 0) also spiral in. The asymmetry of the trajectories leads to a non-Newtonian rheology and diffusive behavior. The latter because a given particle (moving along a finite St open trajectory) suffers a net displacement in the transverse direction after a single interaction. A sequence of such uncorrelated displacements leads to the particle executing a random walk. The inertial diffusivity tensor is anisotropic on account of differing strengths of interaction in the gradient and vorticity directions. Since the entire region (constituting an in finite area) of closed orbits in the plane of shear spirals onto contact for #finite St, the latter represents a singular surface for the pair-distribution function. The exact form of the pair-distribution function at contact is still, however, indeterminate in the absence of non-hydrodynamic effects. It should also be noted that finite St non-rectilinear flows do not support a spatially uniform number density owing to the cross-streamline inertial migration of particles.

Description: The notion that the continental-scale land-sea contrast is the main reason that monsoon circulation exists has been a long-held belief. The purpose of this paper is to point out that this notion should be substantially modified. The central idea of this notion states that in summer, radiative heating of the continent, say Asia, gives rise to a continental-scale thermal low and surrounding the thermal low in its southeast direction the low level wind flows in from south-west. This low-level inflow creates a convergence of moisture, which maintains the cumulus convection. And in winter, radiative cooling of continent gives rise to a thermal high and to its southeast the low-level wind is from northeast. The mechanism in this interpretation does undoubtedly exist. However, this mechanism, though believed to be the main driving force of monsoon, has not been tested in numerical experiments. There has been an increasing recognition in the recent years that monsoon is inextricably tied to the heating in the intertropical convergence zone (ITCZ). We propose that the main cause of monsoon is ITCZ's being substantially away from the equator. A brief qualitative explanation of why the ITCZ can be a source of monsoon circulation can be offered based on the circulation field forced by the ITCZ heating. The existence of the ITCZ's does not always have to rely on land-sea contrast on the continental scale. This is hinted in the fact that in February the ITCZ close to Australia (and its associated monsoon circulation) covers a longitudinal range several times as long as that of Australia and thus cannot possibly be caused mainly by the land-sea contrast associated with Australia. Yet, this cannot be used as a proof that the ITCZ in the Asian summer monsoon is not mainly due to land-sea contrast. One of the purposes of this work is to provide a convincing proof. In this work the role of land-sea contrast in the origin of monsoon is examined through numerical simulation with the Goddard general circulation model. The Asian and Australian monsoon circulations are obtained in a four-year integration and then the integration is repeated with Asia, the maritime continent, and Australia replaced by ocean. The sea surface temperature (SST) at each affected grid is specified as the SST at the first grid to the east that is an ocean grid in the first experiment. The latter integration shows that the monsoon circulation pattern over where south Asia and Australia were and the surrounding region has largely remained. The results discount land-sea contrast as the main cause of Asian monsoon. A third experiment is the same as the first except that the topography of Asia, the maritime continent, and Australia is reduced to zero. This experiment reveals that the difference between the first two experiments is due more to the removal of topography than to the removal of land-sea contrast. August precipitation is shown averaged over the last three years of each of the three experiments. They show that the Asian monsoon rainy region is largely intact in the second experiment and the difference between the second and the third experiment is mainly in the longitudinal location of the maximum precipitation. Additionally, in Asian and Australian winter monsoons land-sea contrast also plays only a modifying role. Although land-sea contrast plays only a modifying role in Asian and Australian (and Central American including Mexican) monsoons, it is the main reason that ITCZ (and thus monsoon) exists in Africa and South America. Thus, monsoons can be classified into two groups depending on whether land-sea contrast plays a major role.

Description: In dynamic meteorology, singular vectors (SVs) are the structures that maximize a given norm of a forecast perturbation given a tangent linear model and a quadratic constraint on the initial perturbation. That constraint is a prescription of the value of either the same or a different norm applied to the initial perturbations. In the sense that SVs maximize the forecast perturbations according to a specified measure, they may be considered as optimal perturbations. SVs are used to characterize predictability, to identify and correct possible initial errors given forecast errors, to create a set of significant perturbations for ensemble forecasting, or to determine locations for observation targeting. They are a specific application of generally defined singular vectors in mathematics.

Description: Land-atmosphere feedback, by which precipitation-induced soil moisture anomalies affect subsequent precipitation, may be an important element of Earth's climate system, but its very existence has never been demonstrated conclusively at regional to continental scales. Evidence for the feedback is sought in a 50-year observational precipitation dataset covering the United States. The precipitation variance and autocorrelation fields are characterized by features that agree (in structure, though not in magnitude) with those produced by an atmospheric general circulation model (AGCM). Because the model-generated features are known to result from land-atmosphere feedback alone, the observed features are highly suggestive of the existence of feedback in nature.

Description: The focus of this project was to study the physical processes that govern tachocline dynamics and structure. Specific features explored included stratification, shear, waves, and toroidal and poloidal background fields. In order to address recent theoretical work on anisotropic mixing and dynamics in the tachocline, we were particularly interested in such anisotropic mixing for the specific tachocline processes studied. Transition to turbulence often shapes the largest-scale features that appear spontaneously in a flow during the development of turbulence. The resulting large-scale straining field can control the subsequent dynamics; therefore, anticipation of the large-scale straining field that results for individual realizations of the transition to turbulence can be important for subsequent dynamics, flow morphology, and transport characteristics. As a result, we paid particular attention to the development of turbulence in the stratified and sheared environment of the tachocline. This is complicated by the fact that the linearly stability of sheared MHD flows is non-self-adjoint, implying that normal asymptotic linear stability theory may not be relevant.

Description: In this paper we consider the robust control of a thermal mixer using multivariable Sliding Mode Control (SMC). The mixer consists of a mixing chamber, hot and cold fluid valves, and an exit valve. The commanded positions of the three valves are the available control inputs, while the controlled variables are total mass flow rate, chamber pressure and the density of the mixture inside the chamber. Unsteady thermodynamics and linear valve models are used in deriving a 5th order nonlinear system with three inputs and three outputs, An SMC controller is designed to achieve robust output tracking in the presence of unknown energy losses between the chamber and the environment. The usefulness of the technique is illustrated with a simulation.

Description: The South China Sea Monsoon Experiment (SCSMEX) was conducted in May-June 1998. One of its major objectives is to better understand the key physical processes for the onset and evolution of the summer monsoon over Southeast Asia and southern China. Multiple observation platforms (e.g., upper-air soundings, Doppler radar, ships, wind profilers, radiometers, etc.) during SCSMEX provided a first attempt at investigating the detailed characteristics of convective storms and air pattern changes associated with monsoons over the South China Sea region. SCSMEX also provided rainfall estimates which allows for comparisons with those obtained from the Tropical Rainfall Measuring Mission (TRMM), a low earth orbit satellite designed to measure rainfall from space. The Goddard Cumulus Ensemble (GCE) model (with 1-km grid size) is used to understand and quantify the precipitation processes associated with the summer monsoon over the South China Sea. This is the first (loud-resolving model used to simulate precipitation processes in this particular region. The GCE-model results captured many of the observed precipitation characteristics because it used a fine grid size. For example, the temporal variation of the simulated rainfall compares quite well to the sounding-estimated rainfall variation. The time and domain-averaged temperature (heating/cooling) and water vapor (drying/ moistening) budgets are in good agreement with observations. The GCE-model-simulated rainfall amount also agrees well with TRMM rainfall data. The results show there is more evaporation from the ocean surface prior to the onset of the monsoon than after the on-~et of monsoon when rainfall increases. Forcing due to net radiation (solar heating minus longwave cooling) is responsible for about 25% of the precipitation in SCSMEX The transfer of heat from the ocean into the atmosphere does not contribute significantly to the rainfall in SCSMEX. Model sensitivity tests indicated that total rain production is reduced 17-18% in runs neglecting the ice phase. The SCSMEX results are compared to other GCE-model-simulated weather systems that developed during other field campaigns (i.e., west Pacific warm pool region, eastern Atlantic region and central USA). Large-scale forcing vie temperature and water vapor tendency, is the major energy source for net condensation in the tropical cases. The effects of large-scale cooling exceed that of large-scale moistening in the west pacific warm pool region and eastern Atlantic region. For SCSMEX, however, the effects of large-scale moistening predominate. Net radiation and sensible and latent hc,it fluxes play a much more important role in the central USA.

Description: CloudSat is a satellite experiment designed to measure the vertical structure of clouds from space. The expected launch of CloudSat is planned for 2004 and, once launched, CloudSat will orbit in formation as part of a constellation of satellites including NASA's Aqua and Aura satellites, a NASA-CNES lidar satellite (P-C) and a CNES satellite carrying a polarimeter (PARASOL). A unique feature that CloudSat brings to this constellation is the ability to fly a precise orbit enabling the fields of view of the CloudSat radar to be overlapped with the P-C lidar footprint and the other measurements of the EOS constellation. The precision of this overlap creates a unique multi-satellite observing system for studying the atmospheric processes essential to the hydrological cycle. The vertical profile of cloud properties provided by CloudSat fills a critical gap in the investigation of feedback mechanisms linking clouds to climate. Measuring the vertical profile of cloud properties requires a combination of active and passive instruments, and this will be achieved by combining the radar data of CloudSat with active and passive data from other sensors of the constellation. This paper describes the underpinning science, and gives an overview of the mission, and provides some idea of the expected products and anticipated application of these products. Notably, the CloudSat mission is expected to provide new knowledge about global cloudiness, stimulating new areas of research on clouds including data assimilation and cloud parameterization. The mission also provides an important opportunity to demonstrate active sensor technology for future scientific and tactical applications. The CloudSat mission is a partnership between NASA/JPL, the Canadian Space Agency, Colorado State University, the US Air Force, and the US Department of Energy.

Description: Several heavy precipitation episodes occurred over Taiwan from August 10 to 13, 1994. Precipitation patterns and characteristics are quite different between the precipitation events that occurred from August 10 and I I and from August 12 and 13. In Part I (Chen et al. 2001), the environmental situation and precipitation characteristics are analyzed using the EC/TOGA data, ground-based radar data, surface rainfall patterns, surface wind data, and upper air soundings. In this study (Part II), the Penn State/NCAR Mesoscale Model (MM5) is used to study the precipitation characteristics of these heavy precipitation events. Various physical processes (schemes) developed at NASA Goddard Space Flight Center (i.e., cloud microphysics scheme, radiative transfer model, and land-soil-vegetation surface model) have recently implemented into the MM5. These physical packages are described in the paper, Two way interactive nested grids are used with horizontal resolutions of 45, 15 and 5 km. The model results indicated that Cloud physics, land surface and radiation processes generally do not change the location (horizontal distribution) of heavy precipitation. The Goddard 3-class ice scheme produced more rainfall than the 2-class scheme. The Goddard multi-broad-band radiative transfer model reduced precipitation compared to a one-broad band (emissivity) radiation model. The Goddard land-soil-vegetation surface model also reduce the rainfall compared to a simple surface model in which the surface temperature is computed from a Surface energy budget following the "force-re store" method. However, model runs including all Goddard physical processes enhanced precipitation significantly for both cases. The results from these runs are in better agreement with observations. Despite improved simulations using different physical schemes, there are still some deficiencies in the model simulations. Some potential problems are discussed. Sensitivity tests (removing either terrain or radiative processes) are performed to identify the physical processes that determine the precipitation patterns and characteristics for heavy rainfall events. These sensitivity tests indicated that terrain can play a major role in determining the exact location for both precipitation events. The terrain can also play a major role in determining the intensity of precipitation for both events. However, it has a large impact on one event but a smaller one on the other. The radiative processes are also important for determining, the precipitation patterns for one case but. not the other. The radiative processes can also effect the total rainfall for both cases to different extents.

Description: Experiments as well as accompanying simulations are described that serve in preparation of a space flight experiment to study the dynamics of miscible interfaces. The investigation specifically addresses the importance of both nonsolenoidal effects as well as nonconventional Korteweg stresses in flows that give rise to steep but finite concentration gradients. The investigation focuses on the flow in which a less viscous fluid displaces one of higher viscosity and different density within a narrow capillary tube. The fluids are miscible in all proportions. An intruding finger forms that occupies a fraction of the total tube diameter. Depending on the flow conditions, as expressed by the Peclet number, a dimensionless viscosity ratio, and a gravity parameter, this fraction can vary between approximately 0.9 and 0.2. For large Pe values, a quasi-steady finger forms, which persists for a time of O(Pe) before it starts to decay, and Poiseuille flow and Taylor dispersion are approached asymptotically. Depending on the specific flow conditions, we observe a variety of topologically different streamline patterns, among them some that leak fluid from the finger tip. For small Pe values, the flow decays from the start and asymptotically reaches Taylor dispersion after a time of O(Pe). Comparisons between experiments and numerical simulations based on the 'conventional' assumption of solenoidal velocity fields and without Korteweg stresses yield poor agreement as far as the Pe value is concerned that distinguishes these two regimes. As one possibility, we attribute this lack of agreement to the disregard of these terms. An attempt is made to use scaling arguments in order to evaluate the importance of the Korteweg stresses and of the assumption of solenoidality. While these effects should be strongest in absolute terms when steep concentration fronts exist, i.e., at large Pe, they may be relatively most important at lower values of Pe. We subsequently compare these conventional simulations to more complete simulations that account for nonvanishing divergence as well as Korteweg stresses. While the exact value of the relevant stress coefficients are not known, ballpark numbers do exist, and their use in the simulations indicates that these stresses may indeed be important. We plan to evaluate these issues in detail by means of comparing a space experiment with corresponding simulations, in order to extract more accurate Korteweg stress coefficients, and to confirm or deny the importance of such stresses.

Description: Six different convective-stratiform separation techniques, including a new technique that utilizes the ratio of vertical and terminal velocities, are compared and evaluated using two-dimensional numerical simulations of a tropical [Tropical Ocean Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA COARE)] and midlatitude continental [Preliminary Regional Experiment for STORM-Central (PRESTORM)] squall line. The simulations are made using two different numerical advection schemes: 4th order and positive definite advection. Comparisons are made in terms of rainfall, cloud coverage, mass fluxes, apparent heating and moistening, mean hydrometeor profiles, CFADs (Contoured Frequency with Altitude Diagrams), microphysics, and latent heating retrieval. Overall, it was found that the different separation techniques produced results that qualitatively agreed. However, the quantitative differences were significant. Observational comparisons were unable to conclusively evaluate the performance of the techniques. Latent heating retrieval was shown to be sensitive to the use of separation technique mainly due to the stratiform region for methods that found very little stratiform rain. The midlatitude PRESTORM simulation was found to be nearly invariant with respect to advection type for most quantities while for TOGA COARE fourth order advection produced numerous shallow convective cores and positive definite advection fewer cells that were both broader and deeper penetrating above the freezing level.

Description: The vibratory conveyor, routinely employed for normal-gravity transport of granular materials, usually consists of a continuous open trough vibrated sinusoidally to induce axial movement of a granular material. Motivated in part by a hypothetical application in zero gravity, we propose a novel modification of the vibratory conveyor based on a closed 2d trough operating in a "slide-conveying" mode, with the granular mass remaining permanently in contact with the trough walls. We present a detailed analysis of the mechanics of transport, based on a rigid-slab model for the granular mass with frictional (Coulomb) slip at the upper and lower walls. The form of the vibration cycle plays a crucial role, and the optimal conveying cycle is not the commonly assumed rectilinear sinusoidal motion. The conveying efficiency for the novel slide conveyor will be presented for several simple vibration cycles, including one believed to represent the theoretical optimum.

Description: The gravity-driven flow of non-neutrally buoyant suspensions is shown to be unstable to spanwise perturbations when the shearing motion generates a density profile that increases with height. The instability is simply due to having heavier material over light. The wavelength of the perturbation is found to be on the order of the thickness of the suspension layer. The parameters important to the problem are the angle of inclination of the layer relative to gravity, the relative density difference between the particles and fluid, the ratio of the particle size to the suspension layer, and the bulk volume fraction of particles. An example showing the growth rate as a function of wave number is shown.

Description: Results of an ongoing effort to quantify the role turbulence in scattering sound in jets are reported. Using a direct numerical simulation database to provide the flow data, ray paths traced through the mean flow are compared with those traced through the actual time evolving turbulent flow. Significant scattering by the turbulence is observed. The most notable effect is that upstream traveling waves that are trapped in the potential core by the mean flow, which acts as a wave guide, easily escape in the turbulent flow. A crude statistical estimate based on ray number density suggests that directivity is modified by the turbulence, but no rigorous treatment of non-uniformities in the high-frequency approximation is attempted.

Description: This study used a two-dimensional coupled land/atmosphere (cloud-resolving) model to investigate the influence of land cover on the water budgets of squall lines in the Sahel. Study simulations used the same initial sounding and one of three different land covers, a sparsely vegetated semi-desert, a grassy savanna, and a dense evergreen broadleaf forest. All simulations began at midnight and ran for 24 hours to capture a full diurnal cycle. In the morning, the latent heat flux, boundary layer mixing ratio, and moist static energy in the boundary layer exhibited notable variations among the three land covers. The broadleaf forest had the highest latent heat flux, the shallowest, moistest, slowest growing boundary layer, and significantly more moist static energy per unit area than the savanna and semi-desert. Although all simulations produced squall lines by early afternoon, the broadleaf forest had the most intense, longest-lived squall lines with 29% more rainfall than the savanna and 37% more than the semi-desert. The sensitivity of the results to vegetation density, initial sounding humidity, and grid resolution was also assessed. There were greater differences in rainfall among land cover types than among simulations of the same land cover with varying amounts of vegetation. Small changes in humidity were equivalent in effect to large changes in land cover, producing large changes in the condensate and rainfall. Decreasing the humidity had a greater effect on rainfall volume than increasing the humidity. Reducing the grid resolution from 1.5 km to 0.5 km decreased the temperature and humidity of the cold pools and increased the rain volume.

Description: A numerical model of the tropical Atlantic ocean is used to investigate the upper layer pathways of the Meridional Overturning Circulation (MOC) in the tropical Atlantic. The main focus of this thesis is on those parts of the tropical circulation that are thought to be important for the MOC return flow, but whose dynamics have not been understood yet. It is shown how the particular structure of the tropical gyre and the MOO act to inhibit the flow of North Atlantic water into the equatorial thermocline. As a result, the upper layers of the tropical Atlantic are mainly fed by water from the South Atlantic. The processes that carry the South Atlantic water across the tropical Atlantic into the North Atlantic as part of the MOO are described here, and three processes that were hitherto not understood are explained as follows: The North Brazil Current rings are created as the result of the reflection of Rossby waves at the South American coast. These Rossby waves are generated by the barotropically unstable North Equatorial Countercurrent. The deep structure of the rings can be explained by merger of the wave's anticyclones with the deeper intermediate eddies that are generated as the intermediate western boundary current crosses the equator. The bands of strong zonal velocity in intermediate depths along the equator have hitherto been explained as intermediate currents. Here, an alternative interpretation of the observations is offered: The Eulerian mean flow along the equator is negligible and the observations are the signature of strong seasonal Rossby waves. The previous interpretation of the observations can then be explained as aliasing of the tropical wave field. The Tsuchyia Jets are driven by the Eliassen-Palm flux of the tropical instability waves. The equatorial current system with its strong shears is unstable and generates tropical instability waves.

Description: The magnetic Kelvin force has been proposed as an artificial gravity to control the orientation of paramagnetic liquid propellants such as liquid oxygen in a microgravity environment. This paper reports experiments performed in the NASA "Weightless Wonder" KC-135 aircraft, through the Reduced Gravity Student Flight Opportunities Program. The aircraft flies through a series of parabolic arcs providing about 25 s of microgravity in each arc. The experiment was conceived, designed, constructed, and performed by the undergraduate student team and their two faculty advisors. Two types of tanks were tested: square-base prismatic tanks 5 cm x 5 cm x 8.6 cm and circular cylinders 5 cm in diameter and 8.6 cm tall. The paramagnetic liquid was a 3.3 molar solution of MnCl2 in water. Tests were performed with each type of tank filled to depths of 1 cm and 4 cm. Each test compared a pair of tanks that were identical except that the base of one was a pole face of a 0.6 Tesla permanent magnet. The Kelvin force attracts paramagnetic materials toward regions of higher magnetic field. It was hypothesized that the Kelvin force would hold the liquid in the bottom of the tanks during the periods of microgravity. The tanks were installed in a housing that could slide on rails transverse to the flight direction. By manually shoving the housing, an identical impulse could be provided to each tank at the beginning of each period of microgravity. The resulting fluid motions were videotaped for later analysis.

Description: An extensive body of research this quarter is documented. Further methodical analysis of temperature residuals in Cryogenic Limb Array Etalon Spectrometer (CLAES) Version 8 level 3AT data show signatures during December 1992 at middle and high northern latitudes that, when compared to Naval Research Laboratory/Mountain Wave Forecast Model (NRL)/(MWFM) mountain wave hindcasts, reveal evidence of long mountain waves in these data over Eurasia, Greenland, Scandinavia and North America. The explicit detection of gravity waves in limb-scanned Cryogenic Infrared Spectrometers and Telescopes for the Atmosphere (CRISTA) temperatures is modeled at length, to derive visibility functions. These insights are used to convert CRISTA gravity wave temperature residuals into data that more closely resemble gravity wave fluctuations detected in data from other satellite instruments, such as Microwave Limb Sounder (MLS), Limb Infrared Monitor of the Stratosphere (LIMS) and Global Positioning System/Meteorology (GPS)/(MET). Finally, newly issued mesospheric temperatures from inversion of CRISTA 15gin emissions are analyzed using a new method that uses separate Kalman fits to the ascending and descending node data. This allows us to study global gravity wave amplitudes at two local times, 12 hours apart. In the equatorial mesosphere, where a large diurnal tidal temperature signal exists, we see modulations of gravity wave activity that are consistent with gravity wave-tidal interactions produced by tidal temperature variability.

Description: Since the advent of meteorological satellites in the 1960's, numerous experiments have been conducted in order to evaluate the impact of these and other data on atmospheric analysis and prediction. Such studies have included both OSE'S and OSSE's. The OSE's were conducted to evaluate the impact of specific observations or classes of observations on analyses and forecasts. Such experiments have been performed for selected types of conventional data and for various satellite data sets as they became available. (See for example the 1989 ECMWF/EUMETSAT workshop proceedings on "The use of satellite data in operational numerical weather prediction" and the references contained therein.) The ODYSSEY were conducted to evaluate the potential for future observing systems to improve Numerical Weather Prediction NWP and to plan for the Global Weather Experiment and more recently for EVANS (Atlas et al., 1985a; Arnold and Day, 1986; Hoffman et al., 1990). In addition, OSSE's have been run to evaluate trade-offs in the design of observing systems and observing networks (Atlas and Emmitt, 1991; Rohaly and Krishnamurti, 1993), and to test new methodology for data assimilation (Atlas and Bloom, 1989).

Description: A dynamic model for a free-piston Stirling convertor is being developed at the NASA Glenn Research Center. The model is an end-to-end system model that includes the cycle thermodynamics, the dynamics, and electrical aspects of the system. The subsystems of interest are the heat source, the springs, the moving masses, the linear alternator, the controller, and the end-user load. The envisioned use of the model will be in evaluating how changes in a subsystem could affect the operation of the convertor. The model under development will speed the evaluation of improvements to a subsystem and aid in determining areas in which most significant improvements may be found. One of the first uses of the end-toend model will be in the development of controller architectures. Another related area is in evaluating changes to details in the linear alternator.

Description: The Department of Energy, Stirling Technology Company (STC), and NASA Glenn Research Center (NASA Glenn) are developing a free-piston Stirling convertor for a high-efficiency Stirling Radioisotope Generator (SRG) for NASA Space Science missions. The SRG is being developed for multimission use, including providing electric power for unmanned Mars rovers and deep space missions. NASA Glenn is conducting an in-house technology project to assist in developing the convertor for space qualification and mission implementation. Recent testing, of 55-We Technology Demonstration Convertors (TDC's) built by STC includes mapping, of a second pair of TDC's, single TDC testing, and TDC electromagnetic interference and electromagnetic compatibility characterization on a nonmagnetic test stand. Launch environment tests of a single TDC without its pressure vessel to better understand the convertor internal structural dynamics and of dual-opposed TDC's with several engineering mounting structures with different natural frequencies have recently been completed. A preliminary life assessment has been completed for the TDC heater head, and creep testing of the IN718 material to be used for the flight convertors is underway. Long-term magnet aging tests are continuing to characterize any potential aging in the strength or demagnetization resistance of the magnets used in the linear alternator (LA). Evaluations are now beginning on key organic materials used in the LA and piston/rod surface coatings. NASA Glenn is also conducting finite element analyses for the LA, in part to look at the demagnetization margin on the permanent magnets. The world's first known integrated test of a dynamic power system with electric propulsion was achieved at NASA Glenn when a Hall-effect thruster was successfully operated with a free-piston Stirling power source. Cleveland State University is developing a multidimensional Stirling computational fluid dynamics code to significantly improve Stirling loss predictions and assist in identifying convertor areas for further improvements. This paper will update the status and results for these efforts.

Description: Six different convective-stratiform separation techniques, including a new technique that utilizes the ratio of vertical and terminal velocities, are compared and evaluated using two-dimensional numerical simulations of a tropical [Tropical Ocean Global Atmosphere Coupled Ocean Atmosphere Response Experiment (TOGA COARE)] and midlatitude continental [Preliminary Regional Experiment for STORM-Central (PRESTORM)] squall line. Comparisons are made in terms of rainfall, cloud coverage, mass fluxes, apparent heating and moistening, mean hydrometeor profiles, CFADs (Contoured Frequency with Altitude Diagrams), microphysics, and latent heating retrieval. Overall, it was found that the different separation techniques produced results that qualitatively agreed. However, the quantitative differences were significant. Observational comparisons were unable to conclusively evaluate the performance of the techniques. Latent heating retrieval was shown to be sensitive to the use of separation technique mainly due to the stratiform region for methods that found very little stratiform rain.

Description: Hurricane Earl developed from a tropical wave that moved into the Gulf of Mexico, which triggered abundant convection. On 1 Sept. 1998, the wave was upgraded directly to a tropical storm. Earl reached hurricane status the next morning. The system moved erratically as it interacted with an upper level short wave trough rotating around a long wave trough to the northeast. The storm made landfall near 0600 UTC on 3 September near Panama City, FL. During August and September 1998, NASA conducted the Third Convection and Moisture Experiment (CAMEX-3). It focused on studying the intensity, track, and impacts at landfall of hurricanes. On the afternoon of 2 September 1998, the NASA ER2 high-altitude aircraft flying at 65,000 feet in tandem with the NASA DC-8 flying at 35,000 feet flew over and through, respectively, the eastern rainbands of Earl near the Florida Panhandle as the storm neared landfall in the region. Two approaches to studying Earl are undertaken here: first, an examination of the source and height of the dry air region using GOES-8 water vapor data and, second, a look into the impact of the dry air entrainment on the system using aircraft remote sensing data.

Description: A pulse thermal loop heat transfer system includes a means to use pressure rises in a pair of evaporators to circulate a heat transfer fluid. The system includes one or more valves that iteratively, alternately couple the outlets the evaporators to the condenser. While flow proceeds from one of the evaporators to the condenser, heating creates a pressure rise in the other evaporator, which has its outlet blocked to prevent fluid from exiting the other evaporator. When the flow path is reconfigured to allow flow from the other evaporator to the condenser, the pressure in the other evaporator is used to circulate a pulse of fluid through the system. The reconfiguring of the flow path, by actuating or otherwise changing the configuration of the one or more valves, may be triggered when a predetermined pressure difference between the evaporators is reached.

Description: The penetration of the MARSIS radar signal into the polar ice mass is modeled to determine the capability of the instrument to locate sub-glacial aquifers. As a ground penetrating radar, the orbiting MARSIS transmits a signal greater than 1 W between 1-5 MHz. In this work we will investigate the effect of ice conductive losses on the radar-detection of subsurface aquifers. Based on wave propagation analysis, it is found that for a bulk ice conductivity below 10-5 S/m, conductive losses in the medium are not significant. However, if the bulk ice conductivity is relatively large (greater than 10-5 S/m), the reflected signal from any deep aquifer will be absorbed as it propagates in the lossy ice medium limiting the probing depth.

Description: Summer tropical synoptic-scale waves over West Africa are quantified by the 850 mb meridional wind component from the NCEP/NCAR reanalysis project. Their relationships with surface precipitation patterns are further explored by applying the data from the Tropical Rainfall Measuring Mission (TRMM) satellite in combination with other satellite observations during 1998-2002. Evident wavelet spectral power peaks are seen within a period of 2.5 - 6 days in both meridional wind and precipitation. The most intense wave signals in meridional wind are concentrated along 15 deg N- 25 deg N. Wave signals in precipitation and corresponding wavelet cross-spectral signals between these two variables, however, are primarily located at 5 deg N- 15 deg N, the latitudes of major summer rain events. There is a tendency for the perturbations in meridional wind component to lag (lead) precipitation signals south (north) of 15 deg N. In some cases, either an in-phase or out-of-phase relationship can even be found between these two variables, suggesting a latitude-dependent horizontal structure for these waves and probably implying two distinct wave-convective coupling mechanisms. Moreover, the lagging relationship (and/or the out-of-phase tendency) is only observed south of 15 deg N during July-September, indicating a strong seasonal preference. This phase relationship is generally consistent with the horizontal wave structures from a composite analysis.

Description: The development of a ground based direct detection Doppler lidar based on the recently described aerosol double edge technique is reported. A pulsed, injection seeded Nd:YAG laser operating at 1064 nm is used to make range resolved measurements of atmospheric winds in the free troposphere. The wind measurements are determined by measuring the Doppler shift of the laser signal backscattered from atmospheric aerosols. The lidar instrument and double edge method are described and initial tropospheric wind profile measurements are presented. Wind profiles are reported for both day and night operation. The measurements extend to altitudes as high as 14 km and are compared to rawinsonde wind profile data from Dulles airport in Virginia. Vertical resolution of the lidar measurements is 330 m and the rms precision of the measurements is a low as 0.6 m/s.

Description: The objectives of this program are to 1) Assess viability of using lithium-ion technology for future NASA applications, with emphasis upon Mars landers and rovers which will operate on the planetary surface; 2) Support the JPL 2003 Mars Exploration Rover program to assist in the delivery and testing of a 8 AHr Lithium-Ion battery (Lithion/Yardney) which will power the rover; 3) Demonstrate applicability of using lithium-ion technologyfor future Mars applications: Mars 09 Science Laboratory (Smart Lander) and Future Mars Surface Operations (General). Mission simulation testing was carried out for cells and batteries on the Mars Surveyor 2001 Lander and the 2003 Mars Exploration Rover.

Description: The SeaWinds on QuikSCAT scatterometer was developed by NASA JPL to measure the speed and direction of ocean surface winds. Simulations performed to estimate the performance of the instrument prior to its launch have indicated that the mid-swath accuracy is worse than that of the rest of the swath. This behavior is a general characteristic of scanning pencil beam scatterometers. For SeaWinds, the accuracy of the rest of the swath, and the size of the swath are such that the instrument meets its science requirements despite mid-swath shortcomings. However, by understanding the problem at mid-swath, we can improve the performance there as well. We discuss the underlying causes of the problem in detail and propose a new wind retrieval algorithm which improves mid-swath performance. The directional discrimination ability of the instrument varies with cross track distance wind speed, and direction. By estimating the range of likely wind directions for each measurement cell, one can optimally apply information from neighboring cells where necessary in order to reduce random wind direction errors without significantly degrading the resolution of the resultant wind field. In this manner we are able to achieve mid-swath RMS wind direction errors as low as 15 degrees for low winds and 10 degrees for moderate to high winds, while at the same time preserving high resolution structures such as cyclones and fronts.

Description: A preliminary propulsion system design for a 10-100 kg. micro/nanospacecraft for generic scientific missions is described. The design drivers resulting in this proposed system are discussed as well as current technology developments.

Description: Ram Burn Observations (RAMBO) is a Department of Defense experiment that observes shuttle Orbital Maneuvering System engine burns for the purpose of improving plume models. On STS-107 the appropriate sensors will observe selected rendezvous and orbit adjust burns.

Description: The onset of the southeast Asian monsoon during 1997 and 1998 was simulated with a coupled mesoscale atmospheric model (MM5) and a detailed land surface model. The rainfall results from the simulations were compared with observed satellite data fiom the TRMM (Tropical Rainfall Measuring Mission) TMI (TRMM Microwave Imager) and GPCP (Global Precipitation Climatology Project). The simulation with the land surface model captured basic signatures of the monsoon onset processes and associated rainfall statistics. The sensitivity tests indicated that land surface processes had a greater impact on the simulated rainfall results than that of a small sea surface temperature change during the onset period. In both the 1997 and 1998 cases, the simulations were significantly improved by including the land surface processes. The results indicated that land surface processes played an important role in modifying the low-level wind field over two major branches of the circulation; the southwest low-level flow over the Indo- China peninsula and the northern cold front intrusion from southern China. The surface sensible and latent heat exchange between the land and atmosphere modified the lowlevel temperature distribution and gradient, and therefore the low-level. The more realistic forcing of the sensible and latent heat from the detailed land surface model improved the monsoon rainfall and associated wind simulation.

Description: The Department of Energy, Lockheed Martin (LM), Stirling Technology Company, and NASA Glenn Research Center (GRC) are developing a high-efficiency Stirling Radioisotope Generator (SRG) for potential NASA Space Science missions. The SRG is being developed for multimission use, including providing spacecraft onboard electric power for NASA deep space missions and power for unmanned Mars rovers. NASA GRC is conducting an in- house supporting technology project to assist in developing the Stirling convertor for space qualification and mission implementation. Preparations are underway for a thermalhacuum system demonstration and unattended operation during endurance testing of the 55-We Technology Demonstration Convertors. Heater head life assessment efforts continue, including verification of the heater head brazing and heat treatment schedules and evaluation of any potential regenerator oxidation. Long-term magnet aging tests are continuing to characterize any possible aging in the strength or demagnetization resistance of the permanent magnets used in the linear alternator. Testing of the magnet/lamination epoxy bond for performance and lifetime characteristics is now underway. These efforts are expected to provide key inputs as the system integrator, LM, begins system development of the SRG. GRC is also developing advanced technology for Stirling convertors. Cleveland State University (CSU) is progressing toward a multi-dimensional Stirling computational fluid dynamics code, capable of modeling complete convertors. Validation efforts at both CSU and the University of Minnesota are complementing the code development. New efforts have been started this year on a lightweight convertor, advanced controllers, high-temperature materials, and an end-to-end system dynamics model. Performance and mass improvement goals have been established for second- and third-generation Stirling radioisotope power systems.

Description: A NASA grant has been awarded to Cleveland State University (CSU) to develop a multi-dimensional (multi-D) Stirling computer code with the goals of improving loss predictions and identifying component areas for improvements. The University of Minnesota (UMN) and Gedeon Associates are teamed with CSU. Development of test rigs at UMN and CSU and validation of the code against test data are part of the effort. The one-dimensional (1-D) Stirling codes used for design and performance prediction do not rigorously model regions of the working space where abrupt changes in flow area occur (such as manifolds and other transitions between components). Certain hardware experiences have demonstrated large performance gains by varying manifolds and heat exchanger designs to improve flow distributions in the heat exchangers. 1-D codes were not able to predict these performance gains. An accurate multi-D code should improve understanding of the effects of area changes along the main flow axis, sensitivity of performance to slight changes in internal geometry, and, in general, the understanding of various internal thermodynamic losses. The commercial CFD-ACE code has been chosen for development of the multi-D code. This 2-D/3-D code has highly developed pre- and post-processors, and moving boundary capability. Preliminary attempts at validation of CFD-ACE models of MIT gas spring and "two space" test rigs were encouraging. Also, CSU's simulations of the UMN oscillating-flow fig compare well with flow visualization results from UMN. A complementary Department of Energy (DOE) Regenerator Research effort is aiding in development of regenerator matrix models that will be used in the multi-D Stirling code. This paper reports on the progress and challenges of this

Description: A liquid drop present on a solid surface can move because of a gradient in wettability along the surface, as manifested by a gradient in the contact angle. The contact angle at a given point on the contact line between a solid and a liquid in a gaseous medium is the angle between the tangent planes to the liquid and the solid surfaces at that point and is measured within the liquid side, by convention. The motion of the drop occurs in the direction of increasing wettability. The cause of the motion is the net force exerted on the drop by the solid surface because of the variation of the contact angle around the periphery. This force causes acceleration of an initially stationary drop, and leads to its motion in the direction of decreasing contact angle. The nature of the motion is determined by the balance between the motivating force and the resisting hydrodynamic force from the solid surface and the surrounding gaseous medium. A wettability gradient can be chemically induced as shown by Chaudhury and Whitesides who provided unambiguous experimental evidence that drops can move in such gradients. The phenomenon can be important in heat transfer applications in low gravity, such as when condensation occurs on a surface. Daniel et al have demonstrated that the velocity of a drop on a surface due to a wettability gradient in the presence of condensation can be more than two orders of magnitude larger than that observed in the absence of condensation. In the present research program, we have begun to study the motion of a drop in a wettability gradient systematically using a model system. Our initial efforts will be restricted to a system in which no condensation occurs. The experiments are performed as follows. First, a rectangular strip of approximate dimensions 10 x 20 mm is cut out of a silicon wafer. The strip is cleaned thoroughly and its surface is exposed to the vapor from an alkylchlorosilane for a period lasting between one and two minutes inside a desiccator. This is done using an approximate line source of the vapor in the form of a string soaked in the alkylchlorosilane. Ordinarily, many fluids, including water, wet the surface of silicon quite well. This means that the contact angle is small. But the silanized surface resists wetting, with contact angles that are as large as 100 degs. Therefore, a gradient of wettability is formed on the silicon surface. The region near the string is highly hydrophobic, and the contact angle decreases gradually toward a small value at the hydrophilic end away from this region. The change in wettability occurs over a distance of several mm. The strip is placed on a platform within a Plexiglas cell. Drops of a suitable liquid are introduced on top of the strip near the hydrophobic end. An optical system attached to a video camera is trained on the drop so that images of the moving drop can be captured on videotape for subsequent analysis. We have performed preliminary experiments with water as well as ethylene glycol drops. Results from these experiments will be presented in the poster. Future plans include the refinement of the experimental system so as to permit images to be recorded from the side as well as the top, and the conduct of a systematic study in which the drop size is varied over a good range. Experiments will be conducted with different fluids so as to obtain the largest possible range of suitably defined Reynolds and Capillary numbers. Also, an effort will be initiated on theoretical modeling of this motion. The challenges in the development of the theoretical description lie in the proper analysis of the region in the vicinity of the contact line, as well as in the free boundary nature of the problem. It is known that continuum models assuming the no slip condition all the way to the contact line fail by predicting that the stress on the solid surface becomes singular as the contact line is approached. One approach for dealing with this issue has been to relax the no-slip boundary condition using the Navier model. Molecular dynamics simulations of the contact line region show that for a non-polar liquid on a solid surface, the no-slip boundary condition is in fact incorrect near the contact line. Furthermore, the same simulations also show that the usual relationship between stress and the rate of deformation breaks down in the vicinity of the contact line. In developing continuum theoretical models of the system, we shall accommodate this knowledge to the extent possible.

Description: This research addresses turbulent gas flows laden with fine solid particles at sufficiently large mass loading that strong two-way coupling occurs. By two-way coupling we mean that the particle motion is governed largely by the flow, while the particles affect the gas-phase mean flow and the turbulence properties. Our main interest is in understanding how the particles affect the turbulence. Computational techniques have been developed which can accurately predict flows carrying particles that are much smaller than the smallest scales of turbulence. Also, advanced computational techniques and burgeoning computer resources make it feasible to fully resolve very large particles moving through turbulent flows. However, flows with particle diameters of the same order as the Kolmogorov scale of the turbulence are notoriously difficult to predict. Some simple flows show strong turbulence attenuation with reductions in the turbulent kinetic energy by up to a factor of five. On the other hand, some seemingly similar flows show almost no modification. No model has been proposed that allows prediction of when the strong attenuation will occur. Unfortunately, many technological and natural two-phase flows fall into this regime, so there is a strong need for new physical understanding and modeling capability. Our objective is to study the simplest possible turbulent particle-laden flow, namely homogeneous, isotropic turbulence with a uniform dispersion of monodisperse particles. We chose such a simple flow for two reasons. First, the simplicity allows us to probe the interaction in more detail and offers analytical simplicity in interpreting the results. Secondly, this flow can be addressed by numerical simulation, and many research groups are already working on calculating the flow. Our detailed data can help guide some of these efforts. By using microgravity, we can further simplify the flow to the case of no mean velocity for either the turbulence or the particles. In fact the addition of gravity as a variable parameter may help us to better understand the physics of turbulence attenuation. The experiments are conducted in a turbulence chamber capable of producing stationary or decaying isotropic turbulence with nearly zero mean flow and Taylor microscale Reynolds numbers up to nearly 500. The chamber is a 410 mm cubic box with the corners cut off to make it approximately spherical. Synthetic jet turbulence generators are mounted in each of the eight corners of the box. Each generator consists of a loudspeaker forcing a plenum and producing a pulsed jet through a 20 mm diameter orifice. These synthetic jets are directed into ejector tubes pointing towards the chamber center. The ejector tubes increase the jet mass flow and decrease the velocity. The jets then pass through a turbulence grid. Each of the eight loudspeakers is forced with a random phase and frequency. The resulting turbulence is highly Isotropic and matches typical behavior of grid turbulence. Measurements of both phases are acquired using particle image velocimetry (PIV). The gas is seeded with approximately 1 micron diameter seeding particles while the solid phase is typically 150 micron diameter spherical glass particles. A double-pulsed YAG laser and a Kodak ES-1.0 10-bit PIV camera provide the PIV images. Custom software is used to separate the images into individual images containing either gas-phase tracers or large particles. Modern high-resolution PIV algorithms are then used to calculate the velocity field. A large set of image pairs are acquired for each case, then the results are averaged both spatially and over the ensemble of acquired images. The entire apparatus is mounted in two racks which are carried aboard NASA's KC-135 Flying Microgravity Laboratory. The rack containing the turbulence chamber, the laser head, and the camera floats freely in the airplane cabin (constrained by competent NASA personnel) to minimize g-jitter.

Description: This quarter was largely devoted to a detailed study of temperature data acquired by the Cryogenic Limb Array Etalon Spectrometer (CLAES) on UARS. Our analysis used the same sequence of methods that have been developed, tested and refined on a more limited subset of temperature data acquired by the CRISTA instrument. We focused on a limited subset of our reasoning that geographical and vertical trends in the small-scale temperature variability could be compared with similar trends observed in November 1994 by the CRISTA-SPAS satellite. Results, backed up with hindcasts from the Mountain Wave Forecast Model (MWFM), reveal strong evidence of mountain waves, most persuasively in the Himalayas on 16-17 November, 1992. These CLAES results are coherent over the 30-50 km range and compare well with MWFM hindcasts for the same period. This constitutes, we believe, the first clear evidence that CLAES explicitly resolved long wavelength gravity waves in its CO2 temperature channel. A series of other tasks, related to mesoscale modeling of mountain waves in CRISTA data and fitting of ground-based and HRDI data on global scales, were seen through to publication stage in peer-reviewed journals.

Description: Ideal cloud-resolving models contain little-accumulative errors. When their domain is so large that synoptic large-scale circulations are accommodated, they can be used for the simulation of the interaction between convective clouds and the large-scale circulations. This paper sets up a framework for the models, using moist entropy as a prognostic variable and employing conservative numerical schemes. The models possess no accumulative errors of thermodynamic variables when they comply with a discrete constraint on entropy conservation and sound waves. Alternatively speaking, the discrete constraint is related to the correct representation of the large-scale convergence and advection of moist entropy. Since air density is involved in entropy conservation and sound waves, the challenge is how to compute sound waves efficiently under the constraint. To address the challenge, a compensation method is introduced on the basis of a reference isothermal atmosphere whose governing equations are solved analytically. Stability analysis and numerical experiments show that the method allows the models to integrate efficiently with a large time step.

Description: Reasor et al. (2004) examined the tilting of a vortex by shear and found that a damping mechanism intrinsic to the dry adiabatic dynamics suppresses departures from an upright state. This realignment occurs through projection of the tilt asymmetry onto two types of vortex Rossby waves: sheared vortex Rossby waves, in which the radial shear of the swirling flow axisymmetrizes tilt asymmetries, and a quasi-mode, or discrete, vortex Rossby wave that, in the absence of damping, causes precession of the upper vortex until it realigns with the lower one (thereafter undergoing repeated cycles of tilting and realignment). With damping, the vortex achieves a downshear-left equilibrium tilt. This paper examines the role of a damped quasi-mode in producing asymmetric vertical motions in a high-resolution simulation of Hurricane Bonnie (1998).

Description: Land-atmosphere feedback, by which (for example) precipitation-induced moisture anomalies at the land surface affect the overlying atmosphere and thereby the subsequent generation of precipitation, has been examined and quantified with many atmospheric general circulation models (AGCMs). Generally missing from such studies, however, is an indication of the extent to which the simulated feedback strength is model dependent. Four modeling groups have recently performed a highly controlled numerical experiment that allows an objective inter-model comparison of land-atmosphere feedback strength. The experiment essentially consists of an ensemble of simulations in which each member simulation artificially maintains the same time series of surface prognostic variables. Differences in atmospheric behavior between the ensemble members then indicates the degree to which the state of the land surface controls atmospheric processes in that model. A comparison of the four sets of experimental results shows that feedback strength does indeed vary significantly between the AGCMs.

Description: The present paper studies the numerical simulation of flows with shock/boundary-layer upstream interaction, under conditions of symmetry in geometry, boundary conditions, and grid. For this purpose, a series of two- and three-dimensional numerical test-cases were carried out. The tests showed that standard numerical schemes, which appear to be symmetry preserving under most flow configurations, produce nonsymmetric perturbations when large separated regions are present. These perturbations are amplified when the core flow is under compression. If the flow-blockage due to separation is sufficiently large, the symmetry of the flow may collapse altogether. Experimental evidence of this numerical behavior is also considered.

Description: The One-Degree Daily (1DD) technique is described for producing globally complete daily estimates of precipitation on a 1 deg x 1 deg lat/long grid from currently available observational data. Where possible (40 deg N-40 deg S), the Threshold-Matched Precipitation Index (TMPI) provides precipitation estimates in which the 3-hourly infrared brightness temperatures (IR T(sub b)) are thresholded and all "cold" pixels are given a single precipitation rate. This approach is an adaptation of the Geostationary Operational Environmental Satellite (GOES) Precipitation Index (GPI), but for the TMPI the IR Tb threshold and conditional rain rate are set locally by month from Special Sensor Microwave/Imager (SSM/I)-based precipitation frequency and the Global Precipitation Climatology Project (GPCP) satellite-gauge (SG) combined monthly precipitation estimate, respectively. At higher latitudes the 1DD features a rescaled daily Television Infrared Observation Satellite (TIROS) Operational Vertical Sounder (TOVS) precipitation. The frequency of rain days in the TOVS is scaled down to match that in the TMPI at the data boundaries, and the resulting non-zero TOVS values are scaled locally to sum to the SG (which is a globally complete monthly product). The time series of the daily 1DD global images shows good continuity in time and across the data boundaries. Various examples are shown to illustrate uses. Validation for individual grid -box values shows a very high root-mean-square error but, it improves quickly when users perform time/space averaging according to their own requirements.

Description: A two-dimensional version of the Goddard Cumulus Ensemble (GCE) Model is used to simulate convective systems that developed in various geographic locations. Observed large-scale advective tendencies for potential temperature, water vapor mixing ratio, and horizontal momentum derived from field campaigns are used as the main forcing. By examining the surface energy budgets, the model results show that the two largest terms are net condensation (heating/drying) and imposed large-scale forcing (cooling/moistening) for tropical oceanic cases. These two terms arc opposite in sign, however. The contributions by net radiation and latent heat flux to the net condensation vary in these tropical cases, however. For cloud systems that developed over the South China Sea and eastern Atlantic, net radiation (cooling) accounts for about 20% or more of the net condensation. However, short-wave heating and long-wave cooling are in balance with each other for cloud systems over the West Pacific region such that the net radiation is very small. This is due to the thick anvil clouds simulated in the cloud systems over the Pacific region. Large-scale cooling exceeds large-scale moistening in the Pacific and Atlantic cases. For cloud systems over the South China Sea, however, there is more large-scale moistening than cooling even though the cloud systems developed in a very moist environment. though For three cloud systems that developed over a mid-latitude continent, the net radiation and sensible and latent heat fluxes play a much more important role. This means the accurate measurement of surface fluxes and radiation is crucial for simulating these mid-latitude cases.