ALBERT

Description: The low-g fluids management group with the Center for Space Construction is engaged in active research on the following topics: gauging; venting; controlling contamination; sloshing; transfer; acquisition; and two-phase flow. Our basic understanding of each of these topics at present is inadequate to design space structures optimally. A brief report is presented on each topic showing the present status, recent accomplishings by our group and our plans for future research. Reports are presented in graphic and outline form.

Description: Experimental and theoretical studies have been conducted to determine critical parameters at the onset of nonlinear counterflow in He II below the lambda point of He-4. Critical temperature differences have been measured in porous media for zero net mass flow and for Darcy permeabilities in the order of magnitude range from 10 to the -10th to 10 to the -8th sq cm. The normalized critical temperature gradients, which covered the liquid temperature range of 1.5 K to the lambda temperature, are found to vary with T proportional to the ratio of the superfluid density to the normal fluid density. This liquid temperature dependence appears to be consistent with duct data which are limited at low temperature by a Reynolds number criterion.

Description: The Multispectral Atmospheric Mapping Sensor (MAMS) is an airborne scanner which collects imagery in eight visible and three IR bands, with ground resolution of 100 m at a nominal aircraft altitude of 20 km. The visible channels provide continuous spectral coverage from 0.42 to 1.05 microns, and the IR channels measure upwelling and reflected radiation at combinations of 3.7, 6.5, 11.1, and 12.5 microns. These include channels nearly identical to the Landsat TM channels 1-4, and similar to those on the VISSR Atmospheric Sounder (VAS) and AVHRR instruments. The high spatial resolution permits investigation of surface thermal and atmospheric water vapor structure at scales not available from current stabilities. Significant modifications have been made since engineering flights in 1985 to improve the radiometric performance. The increased data quality and recent multidisciplinary applications of this data are presented.

Description: In summer 1988 an intensive field experiment was conducted in the vicinity of the Pisgah lava flow in the Mojave Desert. Physical properties such as microtopography, composition, soil moisture, and dielectric constant were measured at five sites representing surfaces with rms heights varying from less than one centimeter to tens of centimeters. In addition, polarimetric radar images at P-band, L-band, and C-band were acquired at three different incidence angles with the NASA/JPL airborne imaging radar polarimeter. Using trihedral corner reflectors deployed in the area prior to imaging, the L- and C-band images were calibrated to provide sigma0 values for each resolution element in the scene. The results of inferring geophysical parameters such as rms surface height and correlation length of the measured surfaces by fitting the observed signatures with those predicted by the small-perturbation model are presented: for smoother surfaces, the rms height values inferred are in good agreement with in situ measurements.

Description: Radar modeling of mangal forest stands in the Sundarbans area of southern Bangladesh is discussed. The modeling uses radar system parameters with forest data on tree height, spacing, biomass, species combinations, and water content (including slightly conductive water) in leaves and trunks of the mangal. For Sundri and Gewa tropical mangal forests, six model components are proposed; the models are required to explain the contributions of various combinations of forest species in the attenuation and scattering of mangal-vegetated nonflooded or flooded surfaces. Statistical data of simulated images have been compared with those of SIR-B images both to refine the modeling procedures and to characterize the model output appropriately. The possibility of delineation of flooded or nonflooded boundaries is discussed.

Description: Surface-air temperature differences are commonly used in a bulk resistance equation for estimating sensible heat flux (H), which is inserted in the one-dimensional energy balance equation to solve for the latent heat flux (LE) as a residual. Serious discrepancies between estimated and measured LE have been observed for partial-canopy-cover conditions, which are mainly attributed to inappropriate estimates of H. To improve the estimates of H over sparse canopies, one- and two-layer resistance models that account for some of the factors causing poor agreement are developed. The utility of the two models is tested with remotely sensed and micrometeorological data for a furrowed cotton field with 20 percent cover and a dry soil surface. It is found that the one-layer model performs better than the two-layer model when a theoretical bluff-body correction for heat transfer is used instead of an empirical adjustment; otherwise, the two-layer model is better.

Description: Experimental treatments in a Douglas-fir forest in NE New Mexico were carried out to determine whether differences in forest canopy chemistry could be detected using data from the Airborne Imaging Spectrometer (AIS-2). Experimental treatments consisted of nitrogen fertilizer additions, sawdust additions, and control plots. After AIS-2 data were collected, the digital number of a given pixel was extracted from each channel, yielding 128 values that were used to form a spectrum. Four spectra were extracted from each treatment plot. Multiple stepwise linear regressions between first and second difference transformations of AIS-2 spectra and the canopy characteristics of biomass, nitrogen concentration, and nitrogen content were performed. The results showed a coefficient of multiple determination of 0.71 between first-difference AIS-2 spectra and measured nitrogen concentration in foliage, indicating that it may be possible to predict nitrogen concentration in Douglas fir using AIS-2 spectra.

Description: The information content in polarimetric SAR images is examined, and the polarimetric image variables containing the information that is important to the classification of terrain features in the images are determined. It is concluded that accurate classification can be done when just over half of the image variables are retained. A reduction in image data dimensionality gives storage savings, and can lead to the improvement of classifier performance. In addition, it is shown that a simplified radar system with only phase-calibrated CO-POL or SINGLE TX channels can give classification performance which approaches that of a fully polarimetric radar.

Description: The spatial and phenological distribution of vegetation was examined in a remote sensing geobotanical study of the East Africa Rift region. Six normalized difference vegetation index (NDVI) AVHRR (advanced very high resolution radiometer) scenes (February, June, and September of both 1984 and 1987) were used as a measure of vegetation presence. NOAA ETOPO-5 elevation and geographic coordinate data were coregistered with the multitemporal NDVI images. Univariate and multivariate statistics indicate that the NDVI values are significantly associated with elevation, latitude, and longitude. This supports the concept of quantifiable, regional gradients that affect large-scale geobotanical studies. A quadratic regression line was fitted to the NDVI, elevation, latitude, and longitudinal data. In this way a regional trend, as affected by elevation, was determined. The deviation of the actual data from this regional trend was displayed as an image, and shows the local climatic and geological influences.

Description: The statistical characteristics of image speckle are reviewed. Existing segmentation techniques that have been used for speckle filtering, edge detection, and texture extraction are sumamrized. The relative effectiveness of each technique is briefly discussed.

Description: Theory and experiments have shown that passive microwave radiometers can be used to measure soil moisture. However, the presence of a vegetative cover alters the measurement that might be obtained under bare conditions. Deterministically accounting for the effect of vegetation and developing algorithms for extracting soil moisture from observations of a vegetable-soil complex present significant obstacles to the practical use of this approach. The presence of a vegetation canopy reduces the sensitivity of passive microwave instruments to soil moisture variations. The reduction in sensitivity, as compared to a bare-soil relationship, increases as microwave frequency increases, implying that the longest wavelength sensors should provide the most information. Sensitivity also decreases as the amount of vegetative wet biomass increases for a given type of vegetation.

Description: Spectral bidirectional reflectances were measured over three natural soil sites using a specially designed radiometer called the Parabola. Two of the sites were bare soils, and the third had a sparse cover of desert scrub. The reflectances were strongly non-Lambertian for all three surfaces, but with markedly different patterns. The measured data were fitted with a quasi-physical reflectance model in which the surface backscattering and forwardscattering are separately formulated. A soil reflectance characterization was obtained by assessing the contributions of the forward, backward, and Lambertian components. This three-parameter characterization produced a satisfactory fit to the measured reflectances and appears promising as a basis for soils categorization.

Description: A radiative-transfer model for the bidirectional reflectance-distribution function (BRDF) shows that snow is moderately anisotropic in the NIR wavelengths. Although the directional-hemispherical albedo of snow decreases as the grains become larger, the forward scattering also increases, with the result that the illumination and viewing geometry must be considered when interpreting physical properties of the surface layer of the snow pack from remote sensing data. Measurements of the BRDF and the transmittance for a variety of snow conditions were made throughout the winter and spring seasons with an SE-590 spectroradiometer, for wavelengths from 0.38 to 1.11 microns. Coincident with these, the surface grain properties were analyzed by stereological methods. A sphere with the same surface-to-volume ratio as the ice grains is used as the equivalent sphere.

Description: Remotely sensed data are being used to estimate foliar chemical content. This paper reviews how stepwise multiple regression and deconvolution have been used to extract chemical information from foliar spectra, and concludes that both methods are useful, but neither is ideal. It is recommended that the focus of research be modeling in the long term and experimentation in the short term. Long-term research should increase our understanding of the interaction between radiation and foliar chemistry so that the focus of research can move from leaf model to canopy model to field experiment. Short-term research should aim to design experiments in which remotely sensed data are used to generate unambiguous and accurate estimates of foliar chemical content.

Description: Surface energy balance components were evaluated by combining satellite-based spectral data with on-site measurements of solar irradiance, air temperature, wind speed, and vapor pressure. Maps of latent heat flux density and net radiant flux density were produced using Landsat TM data for three dates. The TM-based estimates differed from Bowen-ratio and aircraft-based estimates by less than 12 percent over mature fields of cotton, wheat, and alfalfa.

Description: An analysis is made based upon the concept that the velocity fluctuations, and therefore, the Reynolds stresses, driven by the instability of the original flow grow until a new stable state is approached. The Reynolds stresses incorporated into the Orr-Sommerfeld equation are coupled with the main flow such that all the imaginary parts of the complex eigenvalues vanish, i.e., the original instability is eliminated. Using this stabilization principle, it is possible to find the Reynolds stresses as well as the mean velocity for plane Poiseuille flow with the Reynolds number slightly higher than the critical.

Description: Multipolarization P-, L-, and C-band airborne SAR data sets were acquired over a coastal zone and a forested wetland of southern Louisiana. The data sets were used with field-collected surface-parameter data in order to determine the value of SAR systems in assessing and mapping coastal-zone surface features. The coastal-zone surface features in this study are sediments, sediment distribution, and the formation of new isles and banks. Results of the data analysis indicate that the P-band radar with 68-cm wavelength is capable of detecting the submerged sediment if the area is very shallow (i.e., a water depth of less than one meter). The penetration capability of P-band radar is also demonstrated in the forested wetland area. The composition and condition of the ground surface can be detected, as well as the standing water beneath dense tree leaves.

Description: The behavior of the reverse flow ceiling jet against the ventilation flow from 0.58 to 0.87 m/s was investigated in a 1/3 scale model of a wide body aircraft interior. For all tests, strong reverse-flow ceiling jets of hot gases were detected well upstream of the fire. Both thicknesses of the reverse-flow ceiling jet and the smoke layer increased with the fire-crossflow parameter. The thickness of the smoke layer where the smoke flows along the main flow below the reverse-flow ceiling jet was almost twice that of the reverse-flow ceiling jet. Detailed spatial and time-varying temperatures of the gas in the test section were measured, and velocity profiles were also measured using a temperature compensated hot film.

Description: Certain theoretical studies of boundary-layer transition are described, based on high Reynolds numbers and with attention drawn to the various nonlinear interactions and scales present. The article concentrates in particular on theories for which the mean-flow profile is completely altered from its original state. Two- and three-dimensional flow theory and conjectures on turbulent-boundary-layer structures are included. Specific recent findings noted, and in qualitative agreement with experiments, are: nonlinear finite-time break-ups in unsteady interactive boundary layers; strong vortex/wave interactions; and prediction of turbulent boundary-layer displacement- and stress sublayer-thicknesses.

Description: Multiband radiometric data from an airborne imaging thermal scanner are being studied for use in finding buried gravel deposits. The techniques are based on measuring relative differences in the thermal properties between gravel-laden targets and the surrounding gravelless background. These properties are determined from modeling the spectral radiant emittance recorded over both types of surfaces in conjunction with ground measurements of the most significant heat flows above and below the surface. Thermodynamic properties of sampled materials from control sites are determined, and diurnal and annual subsurface heat waves are recorded. Thermal models that account for heat exchange at the surface, as well as varying levels of soil moisture, humidity, and vegetation, are needed for adaptation and modification to simulate the physical and radiative environments of this region.

Description: The advanced solid-state array spectroradiometer (ASAS) is an airborne, off-nadir pointing imaging spectroradiometer used to acquire bidirectional radiance data for terrestrial targets. As its platform aircraft flies over a target the sensor can image the target through a sequence of at least seven fore-to-aft view directions ranging up to 45 deg on either side of nadir. ASAS acquires data for 29 spectral bands in the visible and near-infrared portions of the spectrum with a resolution of 15 nm. ASAS data were recently acquired for a prairie ecosystem and a northern forest ecosystem. The data demonstrate the combined effects of reflectance anisotropy and increased atmospheric path length on off-nadir observations. One result of these effects is a variation in vegetation indices as a function of view direction. Normalized-difference-vegetation-indices for prairie grass, coniferous, and deciduous canopies varied up to 14 percent, 23 percent, and 6 percent, respectively, relative to nadir as a function of view zenith angle along the solar principal plane.

Description: During the First International Land Surface Climatology Project (ISLSCP) Field Experiment (FIFE), 80 gigabytes of image data were generated from a variety of satellite and airborne sensors in a multidisciplinary attempt to study energy and mass exchange between the land surface and the atmosphere. To make these data readily available to researchers with a range of image data handling experience and capabilities, unique image-processing software was designed to perform a variety of nonstandard image-processing manipulations and to derive a set of standard-format image products. The nonconventional features of the software include: (1) adding new layers of geographic coordinates, and solar and viewing conditions to existing data; (2) providing image polygon extraction and calibration of data to at-sensor radiances; and, (3) generating standard-format derived image products that can be easily incorporated into radiometric or climatology models. The derived image products consist of easily handled ASCII descriptor files, byte image data files, and additional per-pixel integer data files (e.g., geographic coordinates, and sun and viewing conditions). Details of the solutions to the image-processing problems, the conventions adopted for handling a variety of satellite and aircraft image data, and the applicability of the output products to quantitative modeling are presented. They should be of general interest to future experiment and data-handling design considerations.

Description: In order to evaluate the importance of surface thermal flux and emissivity variations on surface and boundary layer processes, a technique that uses thermal data from an airborne multispectral scanner to determine the surface skin temperature and thermal emissivity over a regional area has been developed. These values are used to estimate the total flux density emanating from the surface and at the top of the atmosphere. Data from the multispectral atmospheric mapping sensor (MAMS) collected during the First ISLSCP Field Experiment (FIFE) are used to develop the technique, and to show the time and space variability of the flux values. The ground truth data available during FIFE provide a unique resource to evaluate this technique.

Description: A common problem in acquiring ground-truth data for use in microwave interaction modeling is the capture of surface roughness data that are both sampled at distances comparable to a fraction of the wavelength and extensive enough to represent the surface statistics in at least one resolution cell of the microwave remote sensor used. A technique has been developed for acquiring the necessary photogrammetric data using twin 70-mm film cameras mounted on a helicopter boom. The apparatus is described, and the accuracy with which ground surface roughness can be characterized using this device is estimated. Both standard and cross-correlation methods were used for data reduction. Stereogrammetry is compared with a completely automated image-matching technique. Dense disparity images were generated from the helicopter stereo pairs. Using interior orientation parameters supplied by the camera manufacturers, and assuming that exterior orientation parameters remain constant between control target and test field photography, an extremely dense DEM (digital elevation model) for a test field has been derived. Results are compared, and accuracy estimates are presented.

Description: Algorithms to correct for atmospheric-scattering effects in high-spatial resolution land-surface images require the ability to perform rapid and accurate computations of the top-of-atmosphere diffuse radiance field for arbitrarily general surface reflectance distributions (which may be both heterogeneous and non-Lambertian) and atmospheric models. Using three-dimensional radiative transfer (3DRT) theory algorithms are being developed. The methodology used to perform the 3DRT calculations is described. It is shown how these calculations are used to perform atmospheric corrections, and the sensitivity of the retrieved surface reflectances to atmospheric structural parameters is illustrated.

Description: A model for the brightness temperature of a canopy which is based on the physical characteristics of the vegetation is being developed. This model uses the Peake model to relate the scattering coefficient of the vegetation to its radiometric brightness temperature and a discrete scattering model to computer the scattering coefficient. A distorted Born approximation is used to provide a tractable solution to the scattering problem. To verify the theory a comparison has been made with experiments. Radiometric measurements of brightness temperature were made at L- and C-bands over a crop of mature soybeans, along with ground truth measurements. All model input parameters were measured; there were no free parameters to adjust to obtain agreement with the measurements. Results showed that agreement between theory and measurements was reasonable for horizontal polarization at both L- and C-bands with dry underlying soil.

Description: During the spring of 1988, the NASA/JPL multifrequency, multipolarization SAR (synthetic aperture radar) flew in series of calibration experiments over the Goldstone dry lake bed. An array of calibration devices was deployed, including dihedral and trihedral corner reflectors, polarimetric active radar calibrators, passive receivers, and CW (continuous wave) tone generators. The aim of the campaign was to calibrate both amplitude and phase of the resulting SAR images, over long and short time scales. The results of the analysis of the calibration data collected in the spring of 1988 are presented. Trihedral corner reflector signatures and certain image background measures are used to externally calibrate relative amplitude and phase between polarizations at a given frequency, and to calibrate across frequencies. Assessments are made of the calibration accuracy as a function of image frame position for each frequency, and the stability of the radar calibration over long and short time-scales.

Description: An unsupervised clustering algorithm is used to segment multifrequency polarimetric radar data from the NASA/JPL airborne SAR (synthetic aperture radar). Twenty-two parameters are evaluated for their discriminatory capability for each pixel of an image. A clustering analysis is then performed using different subsets of these parameters. This analysis relies on data taken as part of an intensive field experiment during the summer of 1988 in the vicinity of the Pisgah lava flow in the Mojave Desert in southern California. As part of the experiment, extensive ground truth was acquired, including dielectric constant and topography measurements. Segmentation results show good agreement with these measurements.

Description: This paper presents the application of a class of multi-grid methods to the solution of the Navier-Stokes equations for two-dimensional laminar flow problems. The methods consists of combining the full approximation scheme-full multi-grid technique (FAS-FMG) with point-, line- or plane-relaxation routines for solving the Navier-Stokes equations in primitive variables. The performance of the multi-grid methods is compared to those of several single-grid methods. The results show that much faster convergence can be procured through the use of the multi-grid approach than through the various suggestions for improving single-grid methods. The importance of the choice of relaxation scheme for the multi-grid method is illustrated.

Description: Rapid distortion theory is applied to study distortion of homogeneous turbulence subject to two different axisymmetric strain modes: the axisymmetric contraction (AC, nozzle-type flow), and the axisymmetric expansion (AE, diffuser-type flow). The paper explores the differences in effects of the two axisymmetric strain modes on the anisotropy of correlations and structures of turbulence; examines the effect of dilatation on the distortion of turbulence; and provides a theoretical background for turbulence model development. It is found that velocity and vorticity fluctuations are enhanced more efficiently by contraction than by expansion; contraction produces much higher anisotropy in velocity and vorticity than expansion; root-mean-square pressure is slightly reduced during contraction, whereas it increases rapidly during expansion; and vortical structures of rodlike shape develop in a contraction flow, while disklike structures develop in an expansion flow. A simple model that reflects the dependence of turbulence evolution on structural parameters such as the Reynolds-stress anisotropy and total strain is proposed, and is shown to outperform all other models for all cases examined, regardless of the mean strain rate.

Description: An experiment to study the hypothesis that winter cold fronts are more important that violent summer storms in moving sediment in midlatitude coastal river deltas. The experiment is based on Calibrated Airborne Multispectral Scanner (CAMS) data. A new approach for the georeferencing of aircraft imagery has been developed for use in this experiment. The software for the approach uses the aircraft inertial navigation data and is able to automatically map a three-channel 1600 line file through a 90 deg rotation in less than 14 minutes. Results from tests of the technique show rms errors of less than 2/3 of a pixel.

Description: Space photography has been successfully used to extend the space remote sensing data base for environmental monitoring by a decade. In this study of Lake Chad, space photographs were digitized and registered to a topographic base map before water classifications were performed. From 1966 to 1985 over a 21,000-square-kilometer decrease in lake surface area was observed.

Description: Space Shuttle astronauts have used handheld cameras to take about 30,000 photographs of the earth as seen from orbit. These pictures provide valuable, true-color depictions of many geologically significant areas. While the photographs have areal coverages and resolutions similar to the more familiar Landsat MSS and TM images, they differ from the latter in having a wide variety of solar illumination angles and look angles. Astronaut photographs can be used as very small scale aerial photographs for geologic mapping and planning logistical support for field work. Astronaut photography offers unique opportunities, because of the intelligence and training of the on-orbit observer, for documenting dynamic geologic activity such as volcanic eruptions, dust storms, etc. Astronauts have photographed more than 3 dozen volcanic eruption plumes, some of which were not reported otherwise. The stereographic capability of astronaut photography also permits three-dimensional interpretation of geologic landforms which is commonly useful in analysis of structural geology. Astronauts have also photographed about 20 known impact craters as part of project to discover presently unknown examples in Africa, South America, and Australia.

Description: Numerical studies of turbulent flow in an axisymmetric 45-deg-expansion combustor and bifurcated diffuser are presented. The Navier-Stokes equations incorporating a k-epsilon model were solved in a nonorthogonal curvilinear coordinate system. A zonal-grid method, where the flow field was divided into several subsections, was developed. This approach permitted different computational schemes to be used in the various zones. In addition, grid generation was made a more simple task. Boundary overlap and interpolating techniques were used, and an adjustment of the flow variables was required to assure conservation of mass flux. Three finite-differencing methods (hybrid, quadratic upwind, and skew upwind) were used to represent the convection terms. Results were compared with existing experimental data. In general, good agreement between predicted and measured values was obtained.

Description: The environmental, geologic, meteorologic, and oceanographic phenomena documented by earth photography during the Space Shuttle STS-29 mission are reviewed. A map of the nadir point positions of earth-viewing photographs from the mission is given and color photographs of various regions are presented. The mission photographs include atmospheric dust and smoke over parts of Africa and Asia, Sahelian water sites, center pivot irrigation fields in the Middle East, urban smog over Mexico City, isolated burning in the Bolivian Amazon, and various ocean features and cloud formations.

Description: A numerical procedure in which the Navier-Stokes equations are discretized using tightly coupled discretizations of pressure derivatives and continuity equations is used here to extend the range of known terminal velocities of gaseous bubbles in liquids well beyond that in previous investigations. Computations performed for Reynolds numbers up to 2000 and Marangoni numbers up to 1000 show only a modest variation of the scaled bubble velocity between 0.16 and 0.5. The bubble velocity is influenced more by the Marangoni number than by the Reynolds number.

Description: A numerical analysis is performed on thermocapillary buoyancy convection induced by phase change in a liquid droplet. A finite-difference code is developed using an alternating-direction implicit (ADI) scheme. The intercoupling relation between thermocapillary force, buoyancy force, fluid property, heat transfer, and phase change, along with their effects on the induced flow patterns, are disclosed. The flow is classified into three types: thermocapillary, buoyancy, and combined convection. Among the three mechanisms, the combined convection simulates the experimental observations quite well, and the basic mechanism of the observed convection inside evaporating sessile drops is thus identified. It is disclosed that evaporation initiates unstable convection, while condensation always brings about a stable density distribution which eventually damps out all fluid disturbances. Another numerical model is presented to study the effect of boundary recession due to evaporation, and the 'peeling-off' effect (the removal of the surface layer of fluid by evaporation) is shown to be relevant.

Description: Two-dimensional solidification influenced by anisotropic heat conduction has been considered. The interfacial energy balance was derived to account for the heat transfer in one direction (x or y) depending on the temperature gradient in both the x and y directions. A parametric study was made to determine the effects of the Stefan number, aspect ratio, initial superheat, and thermal conductivity ratios on the solidification rate. Because of the imposed boundary conditions, the interface became skewed and sometimes was not a straight line between the interface position at the upper and lower adiabatic walls (spatially nonlinear along the height). This skewness depends on the thermal conductivity ratio k(yy)/k(yx). The nonlinearity of the interface is influenced by the solidification rate, aspect ratio, and k(yy/k(yx).

Description: Remote sensing instruments for monitoring global changes are examined. The use of the earth observing system, a set of instrument platforms in polar, sun-synchoronous orbit that provide coverage of the entire globe, is discussed. The radar and imaging spectrometers utilized to obtain surface measurements are described. Atmospheric data is collected by the atmospheric IR sounder, the tropospheric emission spectrometer, and the stratospheric wind IR limb sounder. Consideration is given to the operation of the microwave limb sounder, the active cavity radiometer, and the TDRSS.

Description: A multiple-time-scale turbulence model of a single point closure and a simplified split-spectrum method is presented. In the model, the effect of the ratio of the production rate to the dissipation rate on eddy viscosity is modeled by use of the multiple-time-scales and a variable partitioning of the turbulent kinetic energy spectrum. The concept of a variable partitioning of the turbulent kinetic energy spectrum and the rest of the model details are based on the previously reported algebraic stress turbulence model. Example problems considered include: a fully developed channel flow, a plane jet exhausting into a moving stream, a wall jet flow, and a weakly coupled wake-boundary layer interaction flow. The computational results compared favorably with those obtained by using the algebraic stress turbulence model as well as experimental data. The present turbulence model, as well as the algebraic stress turbulence model, yielded significantly improved computational results for the complex turbulent boundary layer flows, such as the wall jet flow and the wake boundary layer interaction flow, compared with available computational results obtained by using the standard kappa-epsilon turbulence model.

Description: Multispectral imaging systems provide much of the basic data used by the land and ocean civilian remote-sensing community. There are numerous multispectral imaging systems which have been and are being developed. A common way to compare the radiometric performance of these systems is to examine their noise-equivalent change in reflectance, NE Delta-rho. The NE Delta-rho of a system is the reflectance difference that is equal to the noise in the recorded signal. A comparison is made of the noise equivalent change in reflectance of seven different multispectral imaging systems (AVHRR, AVIRIS, ETM, HIRIS, MODIS-N, SPOT-1, HRV, and TM) for a set of three atmospheric conditions (continental aerosol with 23-km visibility, continental aerosol with 5-km visibility, and a Rayleigh atmosphere), five values of ground reflectance (0.01, 0.10, 0.25, 0.50, and 1.00), a nadir viewing angle, and a solar zenith angle of 45 deg.

Description: Analysis of SAR data of the Greenland ice sheet in summer and winter suggest the use of SAR to monitor the temporal hydrology of ice sheets. Comparisons of each SAR data set with summer Landsat TM imagery show an areal-positive correlation with summer SAR data and a negative correlation with winter SAR data. It is proposed that the summer SAR data are most sensitive to the variable concentrations of free water in the surface snow and that the winter SAR data indicate variations in snow grain size.

Description: The NASA Ames Ecosystem Science and Technology Branch and the U.S. Army Medical Research Institute of Infectious Diseases are conducting research to detect Rift Valley fever (RVF) vector habitats in eastern Africa using active and passive remote-sensing. The normalized difference vegetation index (NDVI) calculated from Landsat TM and SPOT data is used to characterize the vegetation common to the Aedes mosquito. Relationships have been found between the highest NDVI and the 'dambo' habitat areas near Riuru, Kenya on both wet and dry data. High NDVI values, when combined with the vegetation classifications, are clearly related to the areas of vector habitats. SAR data have been proposed for use during the rainy season when optical systems are of minimal use and the short frequency and duration of the optimum RVF mosquito habitat conditions necessitate rapid evaluation of the vegetation/moisture conditions; only then can disease potential be stemmed and eradication efforts initiated.

Description: The possibility of using detailed spatial soil moisture maps as input to a runoff model was investigated. The water balance of a small drainage basin was simulated using a simple storage model. Aircraft microwave measurements of soil moisture were used to construct two-dimensional maps of the spatial distribution of the soil moisture. Data from overflights on different dates provided the temporal changes resulting from soil drainage and evapotranspiration. The study site and data collection are described, and the soil measurement data are given. The model selection is discussed, and the simulation results are summarized. It is concluded that a time series of soil moisture is a valuable new type of data for verifying model performance and for updating and correcting simulated streamflow.

Description: Multiwavelength, multipolarization radar image data can be processed to allow the clear depiction of most information in color composite displays. The first principal component, depicted as variations in image intensity, captures the dominant variation in the data and reduces noise without degradation of spatial detail. Other principal components, depicted primarily as variations in image chromaticity, are derived from spatially filtered data so that image colors can be attributed to significant wavelength and polarization contrasts, and not to speckle. Processing paths can be designed to emphasize wavelength contrasts, polarization contrasts, or overall discrimination of surface features.

Description: An analysis is made of time-dependent covariance of the greenness vegetation index with mapped terrain variables over the Konza Prarie (Kansas) during the 1987 growing season. The analysis was part of an ongoing project to establish appopriate ground-sampling and data-integration strategies for satellite-based monitoring of land surface climate conditions. Greenness images for six dates between May and October were derived from atmospherically corrected thematic mapper (TM) data and coregistered with maps of woody vegetation, fire, and soils. Local variance in greenness peaked in mid-June, falling rapidly until mid-August, and declining gradually thereafter. Greenness images exhibited positive autocorrelation up to distances of 180-210 m, but the dominant scale of pattern occurred at a block size of 60 m by 60 m throughout the growing season. 40-44 percent of total scene variance in July and August was accounted for by the effects of woody vegetation (8.9 percent of the area), prairie burning, and soil type. The effect of these terrain variables was fairly consistent between June and late August and was manifested as additional high-frequency spatial variation in imagery from that period.

Description: The Alaska SAR (synthetic-aperture radar) Facility (ASF) will be capable of receiving, processing, archiving, and producing a variety of SAR image products from three satellite-borne SARs: E-ERS-1 (ESA), J-ERS-1 (NASDA) and Radarsat (Canada). Crucial to the success of the ASF is the Alaska SAR processor (ASP), which will be capable of processing over 200 100-km x 100-km (Seasat-like) frames per day from the raw SAR data, at a ground resolution of about 30 m x 30 m. The processed imagery is of high geometric and radiometric accuracy, and is geolocated to within 500 m. Special-purpose hardware has been designed to execute a SAR processing algorithm to achieve this performance. This hardware is currently undergoing acceptance testing for delivery to the University of Alaska. Particular attention has been devoted to making the operations semi-automated and to providing a friendly operator interface via a computer workstation. The operations and control of the Alaska SAR processor are described.

Description: The problem of estimating areal averages from point measurement has been extensively studied by mining engineers and hydrologists. Its application to satellite measurements has recently been introduced. The semivariaogram has been used in many geostatistical applications to estimate spatial structures of observed properties, such as mineral distributions. An examination is made of snow variations in Colorado from daily snow data collected in 11 SNOTEL stations. The associated semivariogram is estimated. The objective is to estimate the spatial structure of the snow field so that the point data can be used for comparison with, and validation for, satellite measurements.

Description: Spaceborne observations have become increasingly accepted for studies of the hydrological cycle and as a source of information that can be used for management and assessment of water resources over large areas. As data records are extended over long periods of time, using sensors that improve their reliability and quality, the utility of spaceborne observations will become increasingly evident. The Earth Observing System (EOS) concept being studied by NASA holds potential for providing improved means and procedures for hydrological studies, particularly on a global basis, using spaceborne sensors and a comprehensive information system. Continuing emphasis needs to be placed on research exploring the interpretation of remotely sensed observations through the assimilation and use of these data in models of hydrological processes.

Description: A study undertaken to determine SIR-C (Shuttle Imaging Radar-C) science requirements for polarimetric calibration is described. A model describing the effect of miscalibration is presented, followed by a description of specific experiments and the resulting requirements levied on the system. The effects of miscalibration on some commonly used measures which can be extracted from SAR (synthetic aperture radar) data are also discussed. It is concluded that polarimetric calibration requirements are strongly application-dependent, and that a compromise must be made to maximize the domain of application of the data without prohibitive cost. Data summarizing SIR-C polarimetric calibration goals are also presented.

Description: Regression and ratio estimators are used to integrate AVHRR-Global Area Coverage (GAC) and Landsat MSS digital data to estimate forest area in the continental United States. Forestlands are enumerated for the 48 contiguous states using five different AVHRR-GAC data sets. Results indicated that the GAC and MSS forest estimates were not highly correlated. Although the ratio of means and linear regression corrections were, on the average, closer to national U.S. Forest Service forest area estimates, these correction procedures did not consistently improve GAC estimates of forest area. GAC forest area estimates tended to be high in densely forested regions such as the northeast and low in sparsely forested areas.

Description: Consideration is given to the possibility of using Landsat TM data to detect differences and changes in the leaf area index (LAI) of closed canopy pine plantations in central Massachusetts. Field measurements of the LAI were obtained using locally-derived allometric relationships between leaf area and trunk diameter. The relationship between the Landsat TM Band 3/Band 4 ratio and regional differences in LAI is examined. A thinning treatment was conducted to demonstrate the ability of TM to detect a known reduction in LAI. It is concluded that, for the detection of moderate changes and differences in the LAI of closed canopy pine plantations at local scales, the TM sensor is better than field measurements involving allometric equations.

Description: Results are presented from field experiments relating spectral reflectance to intercepted photosynthetically active radiation (PAR) and net CO2 exchange in a natural canopy composed of the marsh cordgrass (Spartina alterniflora). Reflectance measurements made by a hand-held radiometer with Landsat TM spectral wavebands are used to compute remote sensing indices such as the normalized difference vegetation index. Consideration is given to the impact of standing dead canopy material on the relationship between intercepted PAR and spectral vegetation indices and the impact of changes in photosynthetic efficiency on the relationship between vegetation indices and CO2 exchange rates. The results suggest that quantitative remote assessment of photosynthesis and net gas exchange in natural vegetation is feasible, especially if the analysis incorporates information on biological responses to environmental variables.

Description: A method to detect plant water stress by remote sensing is proposed using indices of near-IR and mid-IR wavelengths. The ability of the Leaf Water Content Index (LWCI) to determine leaf relative water content (RWC) is tested on species with different leaf morphologies. The way in which the Misture Stress Index (MSI) varies with RWC is studied. On test with several species, it is found that LWCI is equal to RWC, although the reflectances at 1.6 microns for two different RWC must be known to accurately predict unknown RWC. A linear correlation is found between MSI and RWC with each species having a different regression equation. Also, MSI is correlated with log sub 10 Equivalent Water Thickness (EWT) with data for all species falling on the same regression line. It is found that the minimum significant change of RWC that could be detected by appying the linear regression equation of MSI to EWT is 52 percent. Because the natural RWC variation from water stress is about 20 percent for most species, it is concluded that the near-IR and mid-IR reflectances cannot be used to remotely sense water stress.

Description: The emission from a gray radiating medium is analyzed for transient cooling in surroundings at a low temperature. The medium is rectangular with no variations in the direction normal to the cross section. The integral equation for the transient temperature distribution is solved numerically using a two-dimensional Gaussian integration subroutine. The emissive ability for a rectangle at uniform temperature is compared with that for transient cooling where the temperature distribution of the region has reached a fully developed shape, as shown by a separation of variables solution. The two solutions provide the upper and lower bounds for the emittance of a rectangle during transient cooling. The emittances for various aspect ratios are presented as a function of the mean length of the rectangle and are compared with results for a plane layer and a cylinder.

Description: Calculations for n-decane drops evaporating in a spherical cluster surrounded by unvitiated ambient air at atmospheric pressure were performed using two previously proposed cluster models. Both cluster models predict that turbulent transport effects are more important in the case of small clusters. This is due to the smaller volume to surface ratio and thus to the greater transport of hot unvitiated gas to the drops in order to promote evaporation. The results obtained are compared with those of two turbulent models for each one of the 'trapping factors' and similarity models.

Description: The effects of critical layer nonlinearity are considered on spatially growing oblique instability waves on nominally two-dimensional shear layers between parallel streams. The analysis shows that three-dimensional effects cause nonlinearity to occur at much smaller amplitudes than it does in two-dimensional flows. The nonlinear instability wave amplitude is determined by an integro-differential equation with cubic type nonlinearity. The numerical solutions to this equation are worked out and discussed in some detail. The numerical solutions always end in a singularity at a finite downstream distance.

Description: The effects of critical-layer nonlinearity on spatially growing oblique instability waves on compressible shear layers between two parallel streams are considered. The analysis shows that mean temperature nonuniformities cause nonlinearity to occur at much smaller amplitudes than it does when the flow is isothermal. The nonlinear instability wave growth rate effects are described by an integrodifferential equation which bears some resemblance to the Landau equation, in that it involves a cubic-type nonlinearity. The numerical solutions to this equation are worked out and discussed in some detail. Inviscid solutions always end in a singularity at a finite downstream distance, but viscosity can eliminate this singularity for certain parameter ranges.

Description: A simulation is performed of a passive scalar field convected by a rapidly fluctuating velocity field whose correlation time approaches zero. By using a code proposed in a previous study (Chasnov et al., 1988), the turbulence spectrum of the passive temperature field in the conductive subrange is determined. A theoretical model is proposed which explains the result obtained by representing the transfer of scalar variance by an eddy conductivity, whose correlation time is limited by the correlation time of the velocity field.

Description: Long's self-similar vortex is known to have two solutions for each supercritical value of the flow force. Each of those solutions is shown to have a double structure if the flow force is large. The inertial instabilities of one of those large-flow-force limit solutions are investigated, showing that they are related to the instabilities of the Bickley jet in one regime. However, the swirl in the vortex becomes important for long waves, very strongly modifying the sinuous and varicose, Bickley modes. The asymptotic results obtained agree well with the numerical solutions for the sinuous mode, but not for the varicose mode, the difficulty in the latter case being apparently due to mode jumping.

Description: Numerical techniques are developed to solve the Navier-Stokes equations for unsteady incompressible flow. The extension of the finite-difference Galerkin (FDG) method of Stephens et al. (1984) to the continuous-time case in two or three space dimensions is explained, and the numerical implementation of the method is discussed with particular attention to the staggered-MAC-grid primitive-variable discretization, the application of discrete mass balance to avoid problems inherent in FDG schemes, the direct interpretation of the FDG expansion variables as a discrete streamfunction, and a mass-balance approach to two-dimensional problems with throughflow or obstacles. Numerical results are presented graphically for the evolution of asymptotic steady flow in a driven cavity at Reynolds number 400, 1000, or 3200; good agreement with published experimental data is demonstrated, with accurate predictions of secondary-vortex formation from wall bubble recirculations at Reynolds number 1000.

Description: Flow-field measurements of unsteady turbulent flow downstream of a rotating spoked-wheel wake generator were performed in a short-duration light-piston tunnel, and the instantaneous-velocity data were phase averaged based on a signal synchronized with the bar-passing frequency. Mean axial velocities were found to agree well with those obtained from measurements behind a stationary cylinder and to be independent of both Reynolds and bar-passing Strouhal numbers. Reynolds stresses were found to be consistent with related cylinder-wake measurements, but were significantly higher than corresponding measurements obtained in large-scale research turbomachines. Phase-averaged triple velocity correlations were calculated from the digital velocity records, revealing the sign and the magnitude of skewness in the velocity probability density distributions for the two components.

Description: The bifurcation diagram corresponding to the Eckhaus stability curve has been constructed for the one-dimensional Swift-Hohenberg equation in a finite domain. Finite-amplitude solutions with particular spatial wavelength recover linear stability, as predicted by the Eckhaus curve, after a sequence of secondary bifurcations from the branch of solutions with this wavelength. No connectivity between the primary-solution branches is admissible if the stability predicted by this bifurcation diagram is to correspond to the prediction of the Eckhaus analysis. The Eckhaus curve does not exist if nonlinear couplings destroy this pattern. This is demonstrated by analysis of a coupled pair of Swift-Hohenberg equations.

Description: The Landsat Thematic Mapper (TM) imagery of the Basin and Range province of southern Nevada was analyzed to identify and map volcanic rock assemblages at three Tertiary calderas. It was found that the longer-wavelength visible and the NIR TM Bands 3, 5, and 7 provide more effective lithologic discrimination than the shorter-wavelength bands, due partly to deeper penetration of the longer-wavelength bands, resulting in more lithologically driven radiances. Shorter-wavelength TM Bands 1 and 2 are affected more by surficial weathering products including desert varnish which may or may not provide an indirect link to lithologic identity. Guidelines for lithologic analysis of volcanic terrains using Landsat TM imagery are outlined.

Description: Forest stand structure and biomass data were collected using conventional forest inventory techniques in tropical, subtropical, and warm temperate forest biomes. The feasibility of detecting tropical forest successional age class and total biomass differences using Landsat-Thematic mapper (TM) data, was evaluated. The Normalized Difference Vegetation Index (NDVI) calculated from Landsat-TM data were not significantly correlated with forest regeneration age classes in the mountain terrain of the Luquillo Experimental Forest, Puerto Rico. The low sun angle and shadows cast on steep north and west facing slopes reduced spectral reflectance values recorded by TM orbital altitude. The NDVI, calculated from low altitude aircraft scanner data, was significatly correlated with forest age classes. However, analysis of variance suggested that NDVI differences were not detectable for successional forests older than approximately 15-20 years. Also, biomass differences in young successional tropical forest were not detectable using the NDVI. The vegetation index does not appear to be a good predictor of stand structure variables (e.g., height, diameter of main stem) or total biomass in uneven age, mixed broadleaf forest. Good correlation between the vegetation index and low biomass in even age pine plantations were achieved for a warm temperate study site. The implications of the study for the use of NDVI for forest structure and biomass estimation are discussed.

Description: Thematic Mapper data were analyzed with respect to its capability for mapping the complex structure and dynamics of suspended matter distribution in the coastal area of the German Bight (North Sea). Three independent pieces of information were found by factor analysis of all seven TM channels: suspended matter concentration, atmospheric scattering, and sea surface temperature. For the required atmospheric correction, the signal-to-noise ratios of Channels 5 and 7 have to be improved by averaging over 25 x 25 pixels, which also makes it possible to monitor the aerosol optical depth and aerosol type over cloud-free water surfaces. Near-surface suspended matter concentrations may be detected with an accuracy of factor less than 2 by using an algorithm derived from radiative transfer model calculation. The patchiness of suspended matter and its relation to underwater topography was analyzed with autocorrelation and cross-correlation.

Description: Reflectances calculated from TM data and corrected for atmospheric effects correspond with in situ measured reflectances in the nadir-viewing mode, and are shown to be related to a glacier's mass balance if measured over a period of years. A reflectance of 0.895 for a test site in the Wrangell Mountains, Alaska, was calculated from TM Band 4 (0.76 - 0.90 micron) data and corrected for atmospheric effects. This value was comparable to the in situ reflectance of 0.90 measured in the same 0.76 - 0.90 micron wavelength region. For the same site, a reflectance value of 0.79 derived from integrating over most (0.40 - 3.0 micron) of the reflective portion of the electromagnetic spectrum was quite different from the integrated reflectance of 0.95 calculated for the spectral range 0.40 - 1.0 micron. This demonstrates the importance of using the full reflective energy spectrum for calculating the albedo of snow, and for obtaining a meaningful computation of a glacier's energy and mass balance change.

Description: This paper's studies were initiated under the NASA program for the purpose of conducting the earth sciences research using the Landsat Thematic Mapper. The goals of the program include studies of the factors influencing the growth, health, condition, and distribution of vegetation on the earth; the processes controlling the evolution of the earth's crust; the earth's water budget and the hydrologic processes that operate at local, regional, and global scales; the physical and chemical interaction between different types of surficial materials; and the interaction between the earth's surface and its atmosphere. Twenty-seven domestic and five foreign investigations were initiated in 1985, with the results from most of them already published (one study was terminated due to the delay in the TDRSS). Twelve of the studies addressed hydrology, snow and ice, coastal processes, and near-shore oceanographic phenomena; seven addressed vegetation, soils, or animal habitat; and twelve addressed geologic subjects.

Description: The effects of asymmetry in furnace temperature profile and pulling velocity on the crystal interface shape are demonstrated while neglecting the latent heat of solidification. It is seen that the furnace temperature profile may be varied in order to influence the shape of the melt-crystal interface. An exact thermal analysis is then performed on the Bridgman technique by including the latent heat of solidification as a source term. The exact temperature field required for yielding a flat melt-crystal interface is obtained. The earlier observation regarding the influence of furnace temperature profile on the interface shape is confirmed and a criterion for achieving a flat interface is obtained. Various furnace temperature profiles are selected and their corresponding melt-crystal interface results are presented.

Description: Transient cooling by radiation is analyzed for a cylindrical region filled with axially flowing streams of drops that are becoming solidified. This is of interest for the dissipation of waste heat from orbiting power system in space. The drops absorb, emit, and scatter radiation, and the surroundings are at a lower uniform temperature. The radiative properties are assumed gray, and the scattering is isotropic. The radiating region is a two-phase mixture that remains at the melting temperature of the drops. Its temperature uniformity maintains a high emissive power as energy is lost. This is an advantage over a sensible heat radiator in which the temperature decreases, thereby reducing the emissive power. The results provide the axial length that remains two-phase and the fraction of energy dissipated within this length in which the emissive power has not decreased because of sensible cooling. It is also shown how the radial distribution of the axial velocity of the drops can be modified to increase this energy fraction.

Description: A number of successful applications of a spectral collocation method extended by a multi-domain patching technique are shown. The multi-domain technique can be used to improve resolution for problems with widely disparate scales, and to reduce the ill-conditioning of the spectral operators for problems in which a large number of points are required for distributed resolution. A new nonreflecting outflow boundary treatment for unsteady transition-to-turbulence simulations is also presented, which relies on the multi-domain technique. The role of multi-domain in improving the efficiency of such calculations is discussed.

Description: Spectral element methods are high-order weighted residual techniques based on spectral expansions of variables and geometry for the Navier-Stokes (NS) and transport equations. Here, practical aspects of these methods and their efficient implementation are examined, and several examples of flows in truly complex geometries are presented. The spectral element discretization for NS equations is introduced, and the convergence of the method is addressed. An efficient data management scheme is discussed in the context of parallel processing computations. The method is validated by comparing the spectral element solutions with the exact eigensolutions for the Orr-Sommerfeld equations in two and three dimensions. Computer-aided flow visualizations are presented for an impulsive flow past a sharp edge wedge. Three-dimensional states of channel flow disrupted by an array of cylindrical eddy promoters are studied, and the results of a direct simulation of the turbulent flow in a plane channel are presented.

Description: Comparison between the microwave polarized difference temperature (MPDT) derived from 37 GHz band data and the normalized difference vegetation index (NDVI) derived from near-infrared and red bands, from several empirical investigations are summarized. These indicate the complementary character of the two measures in environmental monitoring. Overall the NDVI is more sensitive to green leaf activity, whereas the MPDT appears also to be related to other elements of the above-ground biomass. Monitoring of hydrological phenomena is carried out much more effectively by the MPDT. Further work is needed to explain spectral and temporal variation in MPDT both through modelling and field experiments.

Description: Coincident Scanning Microwave Multi-channel Radiometer 37 GHz and Advanced Very High Resolution Radiometer normalized difference vegetation index satellite data have been compared from drought-affected regions of sub-Saharan Africa and northeastern Brazil for the time period of 1980-1985. Although the two satellite data types can be highly correlated, differences between them were found for the Sahel zone in 1985 and for northeastern Brazil from 1984-1985. These findings suggest that scattering or surface roughness contributions may be greater than previously assumed for the 37 GHz microwave data.

Description: The microwave polarization difference temperature (MPDT) prepared from 37 GHz Scanning Multifrequnecy Microwave Radiometer (SMMR) data from the Nimbus-7 satellite allows a unique vision of hydrological features, especially of humid areas. From the series of monthly images from 1979 to 1985 the major hydrological features of four major river basins of South America are briefly examined. These include rivers of the Amazon, La Plata, Orinoco and Sao Francisco, wetlands in the La Plata and Amazon basins and floods in the La Plata basin. Many hydrological features appear much as they do on conventional maps. An interesting prespective can be derived from the ways in which they differ from conventional maps.

Description: Monthly 37 GHz microwave polarization difference temperatures (MPDT) derived from the Nimbus-7 scanning multichannel microwave radiometer (SMMR) for southern Africa from 1979 to 1985 are compared with rainfall and Advanced Very High Resolution Radiometer (AVHRR) normalized difference vegetation index (NDVI) data. MPDT rose sharply during a drought episode which occurred within the period included in the data. The rise was seen not only in the growing season, but also in the dry season MPDT when no actively photosynthetic, water-containing leaves are present. The results suggest that scattering of the emitted microwave radiation by dead and living vegetation is a more important factor than has previously been recognized.

Description: Simplified radiative transfer models for 37 GHz brightness temperatures are developed and applied to observations taken with the Nimbus-7 scanning multichannel microwave radiometer. It is found that, as the vegetation density increases, polarization difference decreases. The polarization difference also decreases with increasing surface roughness. Empirical relations are established between the satellite observations and annual rainfall over Africa and Australia, primary productivity, and actual evaporation. Color-coded maps of primary productivity and evaporation are presented.

Description: A description of the Axial Flow Turbine Research Facility (AFTRF) being built at the Turbomachinery Laboratory of the Pennsylvania State University is presented. The purpose of the research to be performed in this facility is to obtain a better understanding of the rotor/stator interaction, three dimensional viscous flow field in nozzle and rotor blade passages, spanwise mixing and losses in these blade rows, transport of wake through rotor passage, and unsteady aerodynamics and heat transfer of rotor blade row. The experimental results will directly feed and support the analytical and the computational tool development. This large scale low speed facility is heavily instrumented with pressure and temperature probes and has provision for flow visualization and laser Doppler anemometer measurement. The facility design permits extensive use of the high frequency response instrumentation on the stationary vanes and more importantly on the rotating blades. Furthermore it facilitates detailed nozzle wake, rotor wake, and boundary layer surveys. The large size of the rig also has the advantage of operating at Reynolds numbers representative of the engine environment.

Description: A quasi-three-dimensional analysis has been developed for unsteady rotor-stator interaction in turbomachinery. The analysis solves the unsteady Euler or thin-layer Navier-Stokes equations in a body-fitted coordinate system. It accounts for the effects of rotation, radius change, and stress-surface thickness. The Baldwin-Lomax eddy-viscosity model is used for turbulent flows. The equations are integrated in time using an explicit four-stage Runge-Kutta scheme with a constant time step. Implicit residual smoothing is used to increase the stability limit of the time-accurate computations. The scheme is described, and stability and accuracy analyses are given.

Description: A 3-D model was developed for simulating multistage turbomachinery flows using supercomputers. This average passage flow model described the time averaged flow field within a typical passage of a bladed wheel within a multistage configuration. To date, a number of inviscid simulations were executed to assess the resolution capabilities of the model. Recently, the viscous terms associated with the average passage model were incorporated into the inviscid computer code along with an algebraic turbulence model. A simulation of a stage-and-one-half, low speed turbine was executed. The results of this simulation, including a comparison with experimental data, is discussed.

Description: A multiphase turbulence closure model is presented which employs one transport equation, namely the turbulence kinetic energy equation. The proposed form of this equation is different from the earlier formulations in some aspects. The power spectrum of the carrier fluid is divided into two regions, which interact in different ways and at different rates with the suspended particles as a function of the particle-eddy size ratio and density ratio. The length scale is described algebraically. A mass/time averaging procedure for the momentum and kinetic energy equations is adopted. The resulting turbulence correlations are modeled under less retrictive assumptions comparative to previous work. The closures for the momentum and kinetic energy equations are given. Comparisons of the predictions with experimental results on liquid-solid jet and gas-solid pipe flow show satisfactory agreement.

Description: Generally, two types of theory are used to describe the field equations for suspensions. The so-called postulated equations are based on the kinetic theory of mixtures, which logically should give reasonable equations for solutions. The basis for the use of such theory for suspensions is tenuous, though it at least gives a logical path for mathematical arguments. It has the disadvantage that it leads to a system of equations which is underdetermined, in a sense that can be made precise. On the other hand, the so-called averaging theory starts with a determined system, but the very process of averaging renders the resulting system underdetermined. A third type of theory is proposed in which the kinetic theory of gases is used to motivate continuum equations for the suspended particles. This entails an interpretation of the stress in the particles that is different from the usual one. Classical theory is used to describe the motion of the suspending medium. The result is a determined system for a dilute suspension. Extension of the theory to more concentrated systems is discussed.

Description: A coherent HF radar system developed by Barrick has successfully measured ocean surface currents near shore. This innovative system, called CODAR, can map the current vector for coastal areas as large as 10,000 sq km. CODAR's range is limited owing to the strong attenuation suffered by HF ground waves. An alternate technique was proposed by Schuler, in which the cross-product power spectrum of two (different frequency) microwave signals is processed. The frequency of the resonant peak corresponds close by to the Doppler shift of an ocean gravity wave traveling toward the radar at the phase velocity, v(sub p). The slight difference between the frequency of the measured resonant delta K peak and the Doppler frequency shift caused by the motion of the gravity wave is attributed to be the current velocity in the pointing direction of the radar. The Microwave Remote Sensing Laboratory (MIRSL) has considered the feasibility of using this technique to measure ocean surface currents from geostationary satellite platforms. Problems are discussed that must be overcome if a satellite current measurement system is to be realized. MIRSL research activities that address some of these problem areas are discussed. Current measurements are presented that were made using a specially-designed C-Band, step-frequency delta K radar. These measurements suggest that progress is being achieved in detecting ocean surface current motion for a wide variety of ocean surface conditions.

Description: The vaporization of a droplet, interacting with its neighbors in a non-dilute spray environment is examined as well as a vaporization scaling law established on the basis of a recently developed theory of renormalized droplet. The interacting droplet consists of a centrally located droplet and its vapor bubble which is surrounded by a cloud of droplets. The distribution of the droplets and the size of the cloud are characterized by a pair-distribution function. The vaporization of a droplet is retarded by the collective thermal quenching, the vapor concentration accumulated in the outer sphere, and by the limited percolative passages for mass, momentum and energy fluxes. The retardation is scaled by the local collective interaction parameters (group combustion number of renormalized droplet, droplet spacing, renormalization number and local ambient conditions). The numerical results of a selected case study reveal that the vaporization correction factor falls from unity monotonically as the group combustion number increases, and saturation is likely to occur when the group combustion number reaches 35 to 40 with interdroplet spacing of 7.5 diameters and an environment temperature of 500 K. The scaling law suggests that dense sprays can be classified into: (1) a diffusively dense cloud characterized by uniform thermal quenching in the cloud; (2) a stratified dense cloud characterized by a radial stratification in temperature by the differential thermal quenching of the cloud; or (3) a sharply dense cloud marked by fine structure in the quasi-droplet cloud and the corresponding variation in the correction factor due to the variation in the topological structure of the cloud characterized by a pair-distribution function of quasi-droplets.

Description: The capability of accurate nonlinear flow analysis of resonance systems is essential in many problems, including combustion instability. Classical numerical schemes are either too diffusive or too dispersive especially for transient problems. In the last few years, significant progress has been made in the numerical methods for flows with shocks. The objective was to assess advanced shock capturing schemes on transient flows. Several numerical schemes were tested including TVD, MUSCL, ENO, FCT, and Riemann Solver Godunov type schemes. A systematic assessment was performed on scalar transport, Burgers' and gas dynamic problems. Several shock capturing schemes are compared on fast transient resonant pipe flow problems. A system of 1-D nonlinear hyperbolic gas dynamics equations is solved to predict propagation of finite amplitude waves, the wave steepening, formation, propagation, and reflection of shocks for several hundred wave cycles. It is shown that high accuracy schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time.

Description: An explicit multistage Runge-Kutta type of time-stepping scheme is used for solving transonic flow past a transport type wing/fuselage configuration. Solutions for both Euler and Navier-Stokes equations are obtained for quantitative assessment of boundary layer interaction effects. The viscous solutions are obtained on both a medium resolution grid of approximately 270,000 points and a find grid of 460,000 points to assess the effects of grid density on the solution. Computed pressure distributions are compared with the experimental data.

Description: Some considerations toward developing numerical procedures for simulating viscous compressible flows are discussed. Both Navier-Stokes and boundary layer field methods are considered. Because efficient viscous-inviscid interaction methods have been difficult to extend to complex 3-D flow simulations, Navier-Stokes procedures are more frequently being utilized even though they require considerably more work per grid point. It would seem a mistake, however, not to make use of the more efficient approximate methods in those regions in which they are clearly valid. Ideally, a general purpose compressible flow solver that can optionally take advantage of approximate solution methods would suffice, both to improve accuracy and efficiency. Some potentially useful steps toward this goal are described: a generalized 3-D boundary layer formulation and the fortified Navier-Stokes procedure.