ALBERT

Description: The knowledge of vegetation dielectric behavior is important in studying the scattering properties of the vegetation canopy and radar backscatter modelling. Until now, a limited number of studies have been published on the dielectric properties in the boreal forest context. This paper presents the results of the dielectric constant as a function of depth in the trunks of two common boreal forest species: black spruce and trembling aspen, obtained from field measurements. The microwave penetration depth for the two species is estimated at C, L, and P bands and used to derive the equivalent dielectric constant for the trunk as a whole. The backscatter modelling is carried out in the case of black spruce and the results are compared with the JPL AIRSAR data. The sensitivity of the backscatter coefficient to the dielectric constant is also examined.

Description: In the present work, the high Reynolds number flow past an inclined plate with a splitter plate placed in its wake is considered numerically. A numerical conformal mapping technique is employed to transform the two-plate system into the same number of cylinders: the flow field is assumed to be two-dimensional. The vortex shedding from the inclined plate is modelled using the discrete vortex method. It is shown that the splitter plate has a profound effect on the development of the flow over a range of values of a suitably defined offset parameter and for a range of positions of the leading edge of the splitter plate. The acoustic field is also calculated and the spectrum reflects the flow results.

Description: Because changes in the Earth's environment have become major global issues, continuous, longterm scientific information is required to assess global problems such as deforestation, desertification, greenhouse effects and climate variations. Global change studies require understanding of interactions of complex processes regulating the Earth system. Space-based Earth observation is an essential element in global change research for documenting changes in Earth environment. It provides synoptic data for conceptual predictive modeling of future environmental change. This paper provides a brief overview of remote sensing technology from the perspective of global change research.

Description: This paper presents the results of the correlation analysis of the Skylab S-193 13.9 GHz Radiometer/Scatterometer data. Computer analysis of the S-193 data shows more than 50 percent of the radiometer and scatterometer data are uncorrelated. The correlation coefficients computed for the data gathered over various ground scenes indicates the desirability of using both active and passive sensors for the determination of various Earth phenomena.

Description: Four radiometric correction methods for the reduction of slope-aspect effects in a Landsat TM data set are tested in a mountainous test site with regard to their physical soundness and their influence on forest classification, as well as on the visual appearance of the scene. Excellent ground reference information and a fine-resolution DEM allowed precise assessment of the applicability of the methods under investigation. The results of the study presented here demonstrate the weakness of the classical cosine correction method for radiometric correction in rugged terrain. The statistical, Minnaert and C-correction approaches, however, yielded an improvement of the forest classification and an impressive reduction of the visual topography effect.

Description: Energy of a high harmonic standing wave tends to be distributed equally over the whole wave even in a stratified medium where the wave's peak amplitude can be much larger near the upper boundary than the lower one. This fact is generalized to the many diverse physical problems which solve second-order differential equations of Sturm-Liouville type. For any such solution y(z) whose sign fluctuates along the z-axis, quantities are found which have the same value between any two neighboring zeros of y. One of the equidistributed quantities for an oscillating fluid sphere is similar to kinetic energy but is identical only in limiting cases. The acoustic midpoint of a cavity can be a unique place where some nonlinear perturbations have extra strength. This may apply to the puzzling solar phenomenon called supergranulation.

Description: Development of a first-order radiative transfer model for predicting backscatter from tree canopies has been underway at the University of Michigan Radiation Laboratory for some time. This model is known as the Michigan Microwave Canopy Scattering (MIMICS) model. This article presents the second-generation MIMICS model (MIMICS II) which accounts for canopies with discontinuous (open) crown layer geometries. MIMICS II models open crown layers by treating the location, size, and shape of the individual tree crowns as random variables. The backscattering coefficients for the canopy are then determined by introducing statistics derived from these parameters into the radiative transfer solution. Application of the radiative transfer equations to the discontinuous canopy geometry is presented. The resulting model is a robust fully polarimetric solution that is applicable over a wide variety of canopy architectures. Model simulations are compared to results generated with the continuous canopy model. The effect of the open crown geometry is found to be most significant at shallow incidence angles and at high frequencies for trees with well-developed crowns. Under these conditions, the gaps in the crown layer give rise to a notable increase in crown layer transmissivity which allows the radar to see through to the lower layers of the canopy more easily, thereby directly affecting the backscatter contribution of the trunks and ground.

Description: Experimental studies are performed on some coniferous trees (Austrian pine, Nordmann spruce, and Norway spruce) to investigate the relation between the tree architecture and radar signal at X-band. For a single tree, the RCS is measured as a function of the scatterer location at 90 deg incidence. It is found that the main scatterers are the leafy branches and the difference between sigma(vv) and sigma(hh) is significant at the upper portion of the tree. At the lower portion of the tree, sigma(vv) and sigma(hh) have almost the same level. For a group of trees the angular trends of sigma(vv) and sigma(hh) are measured. It is found that the levels of sigma(vv) and sigma(hh) are of the same order, but their angular trends vary from one tree species to the other depending on the tree species structure. The interpretation of these experimental results is carried out with the help of a theoretical model which accounts for the structure of the tree. According to this theoretical study, the major scattering trend is due to the leaves, while the perturbation to the angular trend and the level difference between sigma(vv) and sigma(hh) are due to the branch orientation distributions (i.e., the tree architecture).

Description: A computational code for the vorticity-potential method is developed for a three-dimensional bounded vorticity field. The evaluation of the boundary data for the vector potential in the code is improved so that the numerical solution simulates that in an unbounded domain to a high order. The time evolution of two vortex rings and that of an elliptic ring are investigated with this code. The cut-and-connect phenomena of vortex rings are successfully captured. The results are compared with those of asymptotic theory and the experiment. They also highlight the need for additional theoretical and numerical investigations.

Description: The feasibility of using imaging spectrometry in studies of playa evaporites is demonstrated by mapping efflorescent salt crusts in Death Valley (California), using Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) data and a recently developed least-squares spectral band-fitting algorithm. It is shown that it was possible to remotely identify eight different saline minerals, including three borates that have not been previously reported for the Death Valley efflorescent crusts: hydroboracite, pinnoite, and rivadavite. The three borates are locally important phases in the crusts; at least one of them, rivadavite, appears to be forming directly from brine.

Description: Data on chlorophyll content and bathymetry of Lake Tahoe obtained on August 9, 1990 by the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) are compared to concurrent in situ surface and in-water measurements. Measured parameters included profiles of percent transmission of monochromatic light, stimulated chlorophyll fluorescence, photosynthetically available radiation, spectral upwelling and downwelling irradiance, and upwelling radiance. Several analyses were performed illustrating the utility of the AVIRIS over a dark water scene. Image-derived chlorophyll concentration compared extremely well with that measured with bottle samples. A bathymetry map of the shallow parts of the lake was constructed which compares favorably with published lake soundings.

Description: Remotely sensed optical and microwave data can be synergistically used to infer land surface properties. Optical data can be used to estimate surface albedo, radiation absorption by vegetation canopies and their photosynthetic efficiencies. Vegetation canopy reflectance at red and near-infrared wavelengths can be used to correct for vegetation effect on microwave emissivities at low frequencies for estimating soil moisture. Optical data can also provide information about surface and air temperatures, precipitable water vapor, cloud top temperature and its water content. This information can be utilized to correct microwave data for atmospheric effects. These points are illustrated with theoretical analyses and by application to satellite data. The basic physical mechanisms operative at the various wavelengths are also discussed.

Description: Quantitative use of remote multispectral measurements to study and map land surface evapotranspiration has been a challenging issue for the past 20 years. Past work is reviewed against process physics. A simple two-layer combination-type model is used which is applicable to both vegetation and bare soil. The theoretic analysis is done to show which land surface properties are implicitly defined by such evaporation models and to assess whether they are measurable as a matter of principle. Conceptual implications of the spatial correlation of land surface properties, as observed by means of remote multispectral measurements, are illustrated with results of work done in arid zones. A normalization of spatial variability of land surface evaporation is proposed by defining a location-dependent potential evaporation and surface temperature range. Examples of the application of remote based estimates of evaporation to hydrological modeling studies in Egypt and Argentina are presented.

Description: Examples are presented of applications of a fast Fourier transform algorithm to analyze time series of images of Normalized Difference Vegetation Index values. The results obtained for a case study on Zambia indicated that differences in vegetation development among map units of an existing agroclimatic map were not significant, while reliable differences were observed among the map units obtained using the Fourier analysis.

Description: This paper discusses a multisensor satellite approach for the study of hydrological applications. Spectral as well as spatial and temporal characteristics of specific operational and planned instruments applicable to hydrology are presented. A hydrology specific series of sensors are proposed to fill the gaps not covered by the current and planned systems. We have called this hypothetical platform HYDROSAT. In addition, the trade-offs between a geostationary satellite and a polar orbiter are explored.

Description: In the present grid-independent approximate Riemann solver for 2D and 3D flows that are governed by the Euler or Navier-Stokes equations, fluxes on grid faces are obtained by wave decomposition; the assumption of information-propagation in the velocity-difference directions leads to a more accurate resolution of shear and shock waves, when these are are oblique to the grid. The model, which yields significantly greater accuracy in both supersonic and subsonic first-order spatially accurate computations, describes the difference in states at each grid interface by the action of five waves.

Description: The action of weak, streamwise vortices on a plane, incompressible, steady mixing layer is examined in the large Reynolds number limit. The outer, inviscid region is bounded by a vortex sheet to which the viscous region is confined. It is shown that the local linear analysis becomes invalid at streamwise distances O(epsilon sup -1), where (epsilon much less than 1) is the crossflow amplitude, and a new nonlinear analysis is constructed for this region. Numerical solutions of the nonlinear problem show that the vortex sheet undergoes an O(1) change in position and that the solution is ultimately terminated by a breakdown in the numerical procedure. The corresponding viscous layer shows downstream thickening, but appears to remain well behaved up to the terminal location.

Description: This paper investigates the linear stability of confined mixing layers with special emphasis on the effects of heat release and compressibility. The results show that reflection of supersonic disturbances by the walls makes the confined supersonic mixing layer more unstable than the unconfined free shear layer. Decreasing the distance between the walls makes the flow more unstable. However, subsonic disturbances are relatively unaffected by the walls. Heat release and Mach number hardly change the growth rates of supersonic disturbances. The most unstable supersonic disturbances are two-dimensional in rectangular channel flows, but three-dimensional in partially confined flows. Finally, the reactants are not strongly mixed by supersonic instabilities, which mainly disturb one side of the layer.

Description: Spectral local isotropy tests are applied to direct numerical simulation data, mainly at the centerline of a fully developed turbulent channel flow. Despite the small Reynolds number of the simulation, the high-wavenumber behavior of velocity and vorticity spectra is consistent with local isotropy. This consistency is verified by the relationship between streamwise wavenumber spectra and spanwise wavenumber spectra. The high-wavenumber behavior of the pressure spectrum is also consistent with local isotropy and compares favorably with the calculation of Batchelor (1951), which assumes isotropy and joint normality of the velocity field at two points in space. The latter assumption is validated by the shape but not the magnitude of the quadruple correlation of the streamwise velocity fluctuation at small separations. There is only partial support for local spectral isotropy away from the centerline as the magnitude of the mean strain rate increases.

Description: A previously derived analytical model for the small-scale structure of turbulence is reformulated in such a way that the energy spectrum may be computed. The model is an ensemble of two-dimensional (2D) vortices with internal spiral structure, each stretched by an axially symmetric strain flow. Stretching and differential rotation produce an energy cascade to smaller scales in which the stretching represents the effect of instabilities and the spiral structure is the source of dissipation at the end of the cascade. The energy spectrum of the resulting flow may be expressed as a time integration involving only the enstrophy spectrum of the time evolving 2D cross section flow, which may be obtained numerically. Examples are given in which a k exp -5/3 spectrum is obtained by this method. The k exp -5/3 inertial range spectrum is shown to be related to the existence of a self-similar enstrophy preserving range in the 2D enstrophy spectrum. The results are found to be insensitive to time dependence of the strain rate, including even intermittent on-or-off strains.

Description: A k-epsilon model is proposed for wall bonded turbulent flows. In this model, the eddy viscosity is characterized by a turbulent velocity scale and a turbulent time scale. The time scale is bounded from below by the Kolmogorov time scale. The dissipation equation is reformulated using this time scale and no singularity exists at the wall. The damping function used in the eddy viscosity is chosen to be a function of R(sub y) = (k(sup 1/2)y)/v instead of y(+). Hence, the model could be used for flows with separation. The model constants used are the same as in the high Reynolds number standard k-epsilon model. Thus, the proposed model will be also suitable for flows far from the wall. Turbulent channel flows at different Reynolds numbers and turbulent boundary layer flows with and without pressure gradient are calculated. Results show that the model predictions are in good agreement with direct numerical simulation and experimental data.

Description: This paper shows how the radar scattering from vegetated areas is affected by the topography of the surface underneath the vegetation. It is shown, using a discrete scatterer model, that the dominant scattering mechanism may change drastically when the ground surface is tilted relative to the horizontal. In the case of a horizontal ground surface, total scattering may be dominated by scattering off the tree trunks, followed by a reflection off the ground surface. For a relatively small tilt in the ground surface (about 2 deg from horizontal), the ground-trunk interaction term may be replaced by scattering from the branches alone as the dominant scattering mechanism. We also show that the effect of the topography is more pronounced for scattering by longer wavelengths, and discuss the implications on algorithms designed to infer forest woody biomass and soil and vegetation moisture using polarimetric SAR data. The effect of the topography on the scattering behavior from forested areas is illustrated with images acquired by the NASA/JPL three-frequency polarimetric SAR over the Black Forest in Germany.

Description: Results are presented of an experimental and numerical study of a turbulent boundary layer with pressure gradients conducted using the recent 'fringe method' with its numerical advantages and good inflow quality. After an inflow transient good agreement is observed; the differences, of up to 13 percent, are discussed. Moderate deviations from the law of the wall are found in the velocity profiles of the simulation. They are fully correlated with the pressure gradient, are in fair quantitative agreement with the experimental results of Nagano et al. (1992), and are roughly the opposite of uncorrected mixing-length-model predictions. Large deviations from the wall scaling are observed for other quantities, notably for the turbulence dissipation rate. The a(1) structure parameter drops mildly in the upper layer with adverse pressure gradient.

Description: The paper formulates a tensorially consistent near-wall second-order closure model. Redistributive terms in the Reynolds stress equations are modeled by an elliptic relaxation equation in order to represent strongly nonhomogeneous effects produced by the presence of walls; this replaces the quasi-homogeneous algebraic models that are usually employed, and avoids the need for ad hoc damping functions. The model is solved for channel flow and boundary layers with zero and adverse pressure gradients. Good predictions of Reynolds stress components, mean flow, skin friction, and displacement thickness are obtained in various comparisons to experimental and direct numerical simulation data. The model is also applied to a boundary layer flowing along a wall with a 90-deg, constant-radius, convex bend.

Description: A computer code for solving the Reynolds-averaged full Navier-Stokes equations has been developed and applied using H- and C-type grids. The Baldwin-Lomax eddy-viscosity model is used for turbulence closure. The integration in time is based on an explicit four-stage Runge-Kutta scheme. Local time stepping, variable coefficient implicit residual smoothing, and a full multigrid method have been implemented to accelerate steady-state calculations. A grid independence analysis is presented for a transonic rotor blade. Comparisons with experimental data show that the code is an accurate viscous solver and can give very good blade-to-blade predictions for engineering applications.

Description: Land surface hydrologic-atmospheric interactions in humid and semi-arid watersheds were investigated. Active and passive microwave sensors were used to estimate the spatial and temporal distribution of soil moisture at the catchment scale in four areas. Results are presented and discussed. The eventual use of this information in the analysis and prediction of associated hydrologic processes is examined.

Description: Regional ecological studies are considered in the context of the global change problem. The Kursk-91 international experiment is used to illustrate applications of remote sensing data and data bases of field experiments for assessment of parameters of the state of the soil and vegetative cover and subsequent study of biospheric stability on the basis of regular satellite observations.

Description: A numerical study was conducted to investigate how bleed through a two-dimensional slot affects shock-wave induced, boundary-layer separation on a flat plate. This study is based on the ensemble-averaged, compressible, Navier-Stokes equations closed by the Baldwin-Lomax, algebraic turbulence model. The algorithm used to obtain solutions was the implicit, partially split, two-factored scheme of Steger. This study examined the effects of the following parameters in controlling shock-wave induced flow separation: location of slot in relation to where the incident shock wave impinged on the boundary layer, size of slot in relation to the boundary-layer thickness, number of slots, spacings between slots, and strength of the incident shock wave. This study also showed the nature of the very complex flowfield about the slot or slots and how the plenum affects the bleed process. The results of this study are relevant to problems where bleed is used to control shock-wave induced, boundary-layer separation (e.g., inside jet engine inlets and wind tunnels).

Description: Turbulent flow in a 2D channel with repeated rectangular rib roughness was numerically simulated using a low Reynolds number form of the k-epsilon turbulence model. Friction factors and average Stanton numbers were calculated for various pitch to rib height ratios and bulk Reynolds numbers. Comparisons with experiment were generally adequate, with the predictions of friction superior to those for heat transfer. The effect of variable properties for channel flow was investigated, and the results showed a greater effect for friction than for heat transfer. Comparison with experiment yielded no clear conclusions. The turbulence model was also validated for a related problem, that of flow downstream of an abrupt pipe expansion.

Description: The purpose of this note is to construct a local solution that eliminates a residual velocity discontinuity in the inviscid portion of a solution obtained in a recent paper by Goldstein, Leib and Cowley (1992). This result is of importance because it shows that the solution obtained in that paper is entirely non-singular outside the viscous wall boundary layer and that any singularity in the problem will have to arise in the usual way through a breakdown in the viscous boundary layer.

Description: We use a set of equations, sometimes referred to as the 'delta-equations', to approximate the two-dimensional inviscid motion of an initially circular vortex sheet released from rest in a cross-flow. We present numerical solutions of these equations for the case with delta-square = 0 (for which the equations are exact) and for delta-square greater than 0. For small values of the smoothing parameter delta, a spectral filter must be used to eliminate spurious instabilities due to round-off error. Two singularities appear simultaneously in the vortex sheet when delta-square = 0 at a critical time t(c). After t(c), the solutions do not converge as the computational mesh is refined. With delta-square greater than 0, converged solutions were found for all values of delta-square when t is less than t(c), and for all but the two smallest values of delta-square used when t is greater than t(c). Our results show that, when delta-square is greater than 0, the vortex sheet deforms into two doubly branched spirals some time after t(c). The limiting solution as delta approaching 0 clearly exists and equals the delta = 0 solution when t is less than t(c).

Description: The evolution of three-dimensional temporally evolving plane mixing layers through as many as three pairings has been simulated numerically. All simulations were begun from a few low-wavenumber disturbances, usually derived from linear stability theory, in addition to the mean velocity. Three-dimensional perturbations were used with amplitudes ranging from infinitesimal to large enough to trigger a rapid transition to turbulence. Pairing is found to inhibit the growth of infinitesimal three-dimensional disturbances, and to trigger the transition to turbulence in highly three-dimensional flows. The mechanisms responsible for the growth of three-dimensionality and onset of transition to turbulence are described. The transition to turbulence is accompanied by the formation of thin sheets of spanwise vorticity, which undergo secondary rollups. The post-transitional simulated flow fields exhibit many properties characteristic of turbulent flows.

Description: The visible bands of the Landsat Thematic Mapper (TM) sensor were used in an empirical assessment of seagrass biomass on shallow banks near Lee Stocking Island in the Bahamas. The TM bands were transformed to minimize the depth-dependent variance in the bottom reflectance signal. Regression analyses were performed between the transformed bands and field measurements of seagrass standing crop (above-ground biomass). Regression equations using spectral data accounted for up to 80 per cent of the variability in seagrass biomass. The unexplained variance was ascribed to variations in bottom sediment color.

Description: An analytical solution was obtained by Siegel (1991, 1992) for local boundary heat fluxes by a radiating medium at uniform temperature in a 2D rectangular region. It is shown here that, after local fluxes from the medium to the walls have been evaluated, it is very easy to compute local fluxes arriving from the adjacent and opposite walls. This extends the previous analysis and provides convenient relations to include radiation from a black boundary, each side of the rectangle being at a different uniform temperature. The final expressions are helpful in performing spectral calculations that must be made for many spectral bands.

Description: Effects of helicity in three-dimensional incompressible inhomogeneous turbulence are examined with the aid of a two-scale direct-interaction approximation (DIA). The turbulent helicity gives a measure of the reflectional asymmetry in a turbulent flow and its inhomogeneity contributes to the sustainment of large-scale vorticity field in a three-dimensional mean flow. The importance of helicity effects is discussed in the context of flows in a rotating system and swirling flows in a pipe. A three-equation model with the turbulent helicity incorporated is proposed using the theoretical results. The validity of the model is confirmed quantitatively through the application to a decaying swirling flow in a pipe.

Description: The single scattering albedo and optical depth of typical savanna vegetation in Botswana (Africa) have been determined by inverse modelling using satellite observed microwave signatures and surface soil moisture. Soil emissivity was modelled using a multi-layer radiative transfer model. The study is based on large scale surface moisture data and Nimbus/SMMR 6-6 GHz and 37 GHz dual polarized brightness temperatures over a 3-year period. As compared to the optical depths, the derived single scattering albedos displayed only minor seasonal variations, whereas the values fit well within the range reported in the literature from laboratory and field experiments. Both 6-6 and 37GHz optical depths were found to be significantly related to NDVI-values derived from NOAA/AVHRR.

Description: The hypothesis tested was that some part of the ecosystem-dependent variability of vegetation indices was attributable to the effects of light specularly reflected by leaves. 'Minus specular' indices were defined excluding effects of specular light which contains no cellular pigment information. Results, both empirical and theoretical, show that the 'minus specular' indices, when compared to the traditional vegetation indices, potentially provide better estimates of the photosynthetic activity within a canopy - and therefore canopy primary production - specifically as a function of sun and view angles.

Description: A linear mixing model was applied to coarse spatial resolution data from the NOAA Advanced Very High Resolution Radiometer. The reflective component of the 3.55-3.95 micron channel was used with the two reflective channels 0.58-0.68 micron and 0.725-1.1 micron to run a constrained least squares model to generate fraction images for an area in the west central region of Brazil. The fraction images were compared with an unsupervised classification derived from Landsat TM data acquired on the same day. The relationship between the fraction images and normalized difference vegetation index images show the potential of the unmixing techniques when using coarse spatial resolution data for global studies.

Description: Recently, a number of studies have investigated the use of the 37 GHz channels of the Nimbus-7 Scanning Multichannel Microwave Radiometer (SMMR) for vegetation monitoring and for studying synergisms between the SMMR and the NOAA Advanced Very High Resolution Radiometer (AVHRR). The approaches are promising but raise a number of issues concerning interpretation of the results, specifically on the relative effects of vegetation and other surface and atmospheric characteristics on the observed signal. This article analyzes the 37 GHz Microwave Polarization Difference Temperature (MPDT) in terms of its sensitivity to surface and atmospheric parameters. For this, a radiative transfer model is used which indicates some limitations of the MPDT index and suggests the importance of accounting for atmospheric effects in the data analysis. An alternative approach to the MPDT, including lower SMMR frequencies than 37 GHz, is discussed.

Description: A Runge-Kutta formula in time is presently used to advance schemes in which central differences are used to solve the time-dependent Euler equations; a second difference is added near shocks as an artificial viscosity to reduce the given scheme to a first-order upwind one at shocks. A matrix-valued dissipation is introduced and compared with the scalar viscosity; a connection is shown between this artificial viscosity and flux limiters. The use of various flux limiters for this central difference scheme is compared.

Description: Spaceborne laser altimeter systems intended to operate at lunar and Martian orbits are reviewed. Laser altimeter systems capable of long lifetimes with centimeter precision ranging electronics are considered to be essential components of NASA's EOS.

Description: The Visiting Investigator Program (VIP) developed at NASA-Stennis' Science and Technology Laboratory (STL) allows U.S. industry to use the specialized resources of STL in the fields of remote sensing and GIS, with a view to the development of new commercial processes and improved services. Attention is given to the novel agreement mechanisms developed by NASA to implement VIP. These agreements encompass a memorandum of understanding, a technical exchange agreement, a sponsored-transfer agreement, a proprietary work agreement, and a joint endeavor agreement.

Description: An efficient method for calculating unsteady flows is presented, with emphasis on a modified version of the thin-layer Navier-Stokes equations. Fourier stability analysis is used to illustrate the effect of treating the source term implicitly instead of explicity, as well as to illustrate other algorithmic choices. A 2D circular cylinder (with a Reynolds number of 1200 and a Mach number of 0.3) is calculated. The present scheme requires only about 10 percent of the computer time required by global minimum time stepping.

Description: Flow solutions were analyzed using three visualization tools, FAST, UFAT, and Visual3. The simulation models axisymmetric unsteady flow inside a closed circular cylinder with a rotating lid. The capabilities and limitations of these visualization packages are presented. The versatility of these tools enhances scientific study and presentation of numerical results.

Description: The present time-accurate coupled-solution procedure addresses the chemical nonequilibrium Navier-Stokes equations over a wide Mach-number range uses, in conjunction with the strong conservation form of the governing equations, five unknown primitive variables. The numerical tests undertaken address steady convergent-divergent nozzle flows with air dissociation/recombination, dump combustor flows with n-pentane/air chemistry, and unsteady nonreacting cavity flows.

Description: A generalized anisotropic porous medium approach is developed for modelling the flow, heat and mass transport processes during binary mixture solidification. Transient predictions are obtained using FEM, coupled with an implicit time-marching scheme, for solidification inside a two-dimensional rectangular enclosure. A parametric study focusing attention on the effects of solutal buoyancy and thermal buoyancy is presented. It is observed that three parameters, namely the thermal Rayleigh number, the solutal Rayleigh number, and the relative density change parameter, significantly alter the flow fields in the liquid and the mushy regions. Depending upon the nature of these flow fields, the solute enrichment caused by macrosegregation may occur in the top or the bottom region of the enclosure.

Description: The constant C sub mu and the near-wall damping function f sub mu in the eddy-viscosity relation of the k-epsilon model are evaluated from direct numerical simulation (DNS) data for developed channel and boundary layer flow at two Reynolds numbers each. Various existing f sub mu model functions are compared with the DNS data, and a new function is fitted to the high-Reynolds-number channel flow data. The epsilon-budget is computed for the fully developed channel flow. The relative magnitude of the terms in the epsilon-equation is analyzed with the aid of scaling arguments, and the parameter governing this magnitude is established. Models for the sum of all source and sink terms in the epsilon-equation are tested against the DNS data, and an improved model is proposed.

Description: Knowledge of the amount of heat that is conducted, advected and radiated between an ampoule and the furnace is important for understanding vertical Bridgman crystal growth. This heat transfer depends on the temperature, emissivities and geometries of both the furnace and ampoule, as well as the choice of ambient gas inside the furnace. This paper presents a method which directly measures this heat transfer without the need to know any physical properties of the furnace, the ampoule, or the gaseous environment. Data are given for one specific furnace in which this method was used.

Description: Topographic profiles at 25- and 5-cm horizontal resolution for three sites along a lava flow on Kilauea Volcano are presented, and these data are used to illustrate techniques for surface roughness analysis. Height and slope distributions and the height autocorrelation function are evaluated as a function of varying lowpass filter wavelength for the 25-cm data. Rms slopes are found to increase rapidly with decreasing topographic scale and are typically much higher than those found by modeling of Magellan altimeter data for Venus. A more robust description of the surface roughness appears to be the ratio of rms height to surface height correlation length. For all three sites this parameter falls within the range of values typically found from model fits to Magellan altimeter waveforms. The 5-cm profile data are used to estimate the effect of small-scale roughness on quasi-specular scattering.

Description: A comparative description is presented for the least-squares FEM (LSFEM) for 2D steady-state pure convection problems. In addition to exhibiting better control of the streamline derivative than the streamline upwinding Petrov-Galerkin method, numerical convergence rates are obtained which show the LSFEM to be virtually optimal. The LSFEM is used as a framework for an iteratively reweighted LSFEM yielding nonoscillatory and nondiffusive solutions for problems with contact discontinuities; this method is shown to convect contact discontinuities without error when using triangular and bilinear elements.

Description: The directions of the triadic energy transfers assumed in a previous paper agree with the spectral closures, in a similarity range, if the exponent of the power-law energy spectrum is less than unity. The helical interactions showing a strong local energy transfer when the triad is nonlocal, sum up to a reverse cascade unless the spectrum falls off faster than a -7/3 power of the wave number. The energy cascades from each type of helical interaction are calculated for a -5/3 inertial range using the eddy damped quasinormal Markovian model. One type of interaction is responsible for 86 percent of the cascade. The contributions of the two classes of helical interactions to the subgrid-scale eddy viscosity are presented, together with the contributions from the forward and reverse cascading interactions. An application of the assumption on the triadic energy transfers to turbulence under rapid rotation gives a simple argument for the tendency toward nonlinear two-dimensionalization of the flow.

Description: A database obtained from direct numerical simulation of a turbulent channel flow is analyzed to extract the streamwise component of the propagation velocity V of velocity, vorticity, and pressure fluctuations from their space-time correlations. A surprising result is that V is approximately the same as the local mean velocity for most of the channel, except for the near-wall region. For y(+) less than 15, V is virtually constant, implying that perturbations of all flow variables propagate like waves near the wall. In this region V is 55 percent of the centerline velocity Uc for velocity and vorticity perturbations and 75 percent of U sub c for pressure perturbations. This is equal to U at y(+) = 15 for velocity and vorticity perturbations, and equal to U at y(+) = 20 for pressure perturbations, indicating that the dynamics of the nearwall turbulence is controlled by turbulence structures present near y(+) about 15-20. Scale dependence of V is also examined by analyzing the bandpass-filtered flow fields. This paper contains comprehensive documentation on the propagation velocities, which should prove useful in the evaluation of Taylor's hypothesis.

Description: Direct numerical simulation data for the lateral velocity derivative delta(u)/delta(y) at the centerline of a fully developed turbulent channel flow provide reasonable support for Wyngaard's analysis of the error involved in measuring this quantity using parallel hot wires. Numerical data in the wall region of the channel flow also provide a useful indication of how to select the separation between the wires. Justification for this choice is obtained by comparing several measured statistics of delta(u)/delta(y) with the corresponding numerical data.

Description: The nonlinear instability of the boundary layer on a heated flat plate placed in an oncoming flow is investigated. Such flows are unstable to stationary vortex instabilities and inviscid travelling wave disturbances governed by the Taylor-Goldstein equation. For small temperature differences the Taylor-Goldstein equation reduces to Rayleigh's equation. When the temperature difference between the wall and free stream is small the preferred mode of instability is a streamwise vortex. It is shown in this case that the vortex, assumed to be of small wavelength, restructures the underlying mean flow to produce a profile which can be massively unstable to inviscid travelling waves. The mean state is shown to be destabilized if the Prandtl number is less than unity.

Description: A new era is dawning in the ability to predict convection heat transfer in the turbine gas path. We feel that the technical community now has the capability to mount a major assault on this problem, which has eluded significant progress for a long time. In this paper we hope to make a case for this bold statement by reviewing the state of the art in three major and related areas, which we believe are indispensable to the understanding and accurate prediction of turbine gas path heat transfer: configuration-specific experiments, fundamental physics and model development, and code development. We begin our review with the configuration-specific experiments, whose data have provided the big picture and guided both the fundamental modeling research and the code development. Following that, we examine key modeling efforts and comment on what will be needed to incorporate them into the codes. In this region we concentrate on bypass transition, 3D endwalls, and film cooling. We then review progress and directions in the development of computer codes to predict turbine gas path heat transfer. Finally, we cite examples and make observations on the more recent efforts to do all this work in a simultaneous, interactive, and more synergistic manner. We conclude with an assessment of progress, suggestions for how to use the current state of the art, and recommendations for the future.

Description: An improved method of estimating fractal surface dimensions has been developed. The accuracy of this method is illustrated using artificially generated fractal surfaces. A slightly different from usual concept of linear dimension is developed, allowing a direct link between that and the corresponding surface dimension estimate. These methods are applied to a series of images of lava flows, representing a variety of physical and chemical conditions. These include lavas from California, Idaho, and Hawaii, as well as some extraterrestrial flows. The fractal surface dimension estimations are presented, as well as the fractal line dimensions where appropriate.

Description: Three AVHRR-LAC data sets acquired in September 1990 and January 1991 were used to map the forest resources of Madagascar. The island was partitioned into four strata to include: (1) the western hardwoods, (2) the central grasslands, (3) the eastern rainforest, and (4) spiny forest. Each stratum was classified separately using AVHRR-LAC data in conjunction with 1984-1988 Landsat-MSS photoproducts. The results of AVHRR classification indicate that approximately 11 percent of the island is covered by forest. Estimates of forest area, by stratum, are as follows: western hardwoods, 6697 sq km; central grasslands, 2830 sq km; eastern rainforest 34,167 sq km; and spiny forest, 17,224 sq km. The total forest area on the 587,041 sq km island is estimated to be 60,918 sq km. The AVHRR forest map was compared to a mid-1970s land cover map which was developed using Landsat-MSS photoproducts. The average class agreement between the mid 1970s ground reference map and the 1990 AVHRR-LAC map was 78.2 percent, the overall accuracy was 81.1 percent. Much of the per-pixel disagreement between the ground reference and AVHRR maps involved areas identified as forest in the 1970s and as nonforest in 1990.

Description: We describe an application of a scale-space clustering algorithm to the classification of a multispectral and polarimetric SAR image of an agricultural site. After the initial polarimetric and radiometric calibration and noise cancellation, we extracted a 12-dimensional feature vector for each pixel from the scattering matrix. The clustering algorithm was able to partition a set of unlabeled feature vectors from 13 selected sites, each site corresponding to a distinct crop, into 13 clusters without any supervision. The cluster parameters were then used to classify the whole image. The classification map is much less noisy and more accurate than those obtained by hierarchical rules. Starting with every point as a cluster, the algorithm works by melting the system to produce a tree of clusters in the scale space. It can cluster data in any multidimensional space and is insensitive to variability in cluster densities, sizes and ellipsoidal shapes. This algorithm, more powerful than existing ones, may be useful for remote sensing for land use.

Description: Results are presented of an investigation to determine the degree to which digitally processed Landsat TM imagery can be used to discriminate among vegetated lava flows of different ages in the Menengai Caldera, Kenya. A selective series of five images, consisting of a color-coded Landsat 5 classification and four color composites, are compared with geologic maps. The most recent of more than 70 postcaldera flows within the caldera are trachytes, which are variably covered by shrubs and subsidiary grasses. Soil development evolves as a function of time, and as such supports a changing plant community. Progressively older flows exhibit the increasing dominance of grasses over bushes. The Landsat images correlated well with geologic maps, but the two mapped age classes could be further subdivided on the basis of different vegetation communities. It is concluded that field maps can be modified, and in some cases corrected by use of such imagery, and that digitally enhanced Landsat imagery can be a useful aid to field mapping in similar terrains.

Description: The length scales appearing in the relations for the eddy viscosity and dissipation rate in one-equation models were evaluated from direct numerical (DNS) simulation data for developed channel and boundary-layer flow at two Reynolds numbers each. To prepare the ground for the evaluation, the distribution of the most relevant mean-flow and turbulence quantities is presented and discussed, also with respect to Reynolds-number influence and to differences between channel and boundary-layer flow. An alternative model is tested as near wall component of a two-layer model by application to developed-channel, boundary-layer and backward-facing-step flows.

Description: Measurements of the C-band (wavelength = 5 cm) radar cross section of an area in the Brooks Range foothills on the North Slope of Alaska using images from the ERS-1 satellite show significant temporal changes. These changes are strongly correlated with elevation and hillslope orientation and are greatest on some of the elevated areas and weaker in river drainages. By constructing 'difference images' using various image pairs, and by analyzing climatological and hydrological data from the site, we conclude that the radar backscatter changes are largely due to changes in soil and vegetation liquid water content induced by freeze/thaw events. The correlation with topography in the difference images arises from the dependence of vegetation, organic layer thickness, and volumetric water content on hillslope position and orientation. These results demonstrate the viability of radar backscatter intensity comparisons using repeat-pass images as a means of change detection.

Description: Measured raw transfer interactions from which local energy transfer is argued to result are summed in a way that directly indicates the scale disparity (s) of contributions to the net energy flux across the spectrum. It is found that the dependence upon s closely follows the s exp -4/3 form predicted by classical arguments. As a result, it is concluded that direct numerical simulation measurements lend support to the classical Kolmogorov phenomenology of local interactions and local transfer in an inertial range.

Description: The present numerical study of unsteady, low Reynolds number flow past a 2D airfoil attempts to ascertain the bifurcation sequence leading from simple periodic to complex aperiodic flow with rising Reynolds number, as well as to characterize the degree of chaos present in the aperiodic flow and assess the role of numerics in the modification and control of the observed bifurcation scenario. The ARC2D Navier-Stokes code is used in an unsteady time-accurate mode for most of these computations. The system undergoes a period-doubling bifurcation to chaos as the Reynolds number is increased from 800 to 1600; its chaotic attractors are characterized by estimates of the fractal dimension and partial Liapunov exponent spectra.

Description: Direct numerical simulations of temporally evolving plane mixing layers undergoing as many as three pairings have been examined for evidence of spanwise scale change. All simulations were begun from a few low-wavenumber disturbances, usually derived from linear stability theory, in addition to the mean velocity. The amplitude of the initial three-dimensional disturbances varied from infinitesimal to large enough to trigger a rapid transition to turbulence. The time required for a change of characteristic spanwise scale with infinitesimal three-dimensional disturbances was found to be very long, requiring three or more pairings to complete a doubling of the spanwise scale. Stronger three-dimensionality can produce more rapid scale changes, but it is also likely to trigger transition to turbulence. No evidence was found for a change from an organized array of rib vortices at one spanwise scale to a similar array at a larger spanwise scale.

Description: A method is presented which allows diffusion terms to be treated implicitly in the lower-upper (LU) algorithm (which is a commonly used method for solving 'compressible' Euler and Navier-Stokes equations) so that the algorithm's good stability properties will not be impaired. The new method generalizes the concept of LU factorization from that associated with the sign of eigenvalues to that associated with backward- and forward-difference operators without regard to eigenvalues. The method is verified in a turbulent boundary layer study.

Description: We consider the effects of a small-amplitude, steady, streamwise vorticity field on the flow over an infinitely thin flat plate in an otherwise uniform stream. We show how the initially linear perturbation, ultimately leads to a small-amplitude but nonlinear cross flow far downstream from the leading edge. This motion is imposed on the boundary-layer flow and eventually causes the boundary layer to separate. The streamwise velocity profiles within the boundary layer become inflexional in localized spanwise regions just upstream of the separation point. The flow in these regions is therefore susceptible to rapidly growing inviscid instabilities.

Description: This research attempts to map small-scale vegetation changes in Mexico. Forty-eight weeks of coarse resolution Advanced Very High Resolution Radiometer Normalized Difference Vegetation Index (NDVI), a digitized climax vegetation map, land cover samples from space shuttle photographs and actual vegetation samples were used as inputs. Principal components analyses and a clustering algorithm were applied to the NDVI data to generate a single layer that was stratified by the climax vegetation zones map. The purpose is to create a new layer that differentiates climax vegetation (hypothesized potential vegetation) from non-climax vegetation land covers. One of the keys to developing a present-day vegetation map was differentiating intrazone land covers based on the stratification; as great as 75% of the sampled land cover types differed from the climax vegetation. The present-day vegetation map achieved 80% classification accuracy when calculated from available ground reference data. About 55% of the temperate zones and 37% of the tropical zones were found to contain original climax vegetation. Most changes coincide with areas of major agricultural activity.

Description: The photographic and videographic documentation of the Earth during STS-46 mission has enhanced the Space Shuttle Earth Observations Project (SSEOP) database. Increasing numbers of scientists are using this database; many are downloading the imagery from our electronic database for specific scientific analyses. We believe the scientific returns of the Earth Observations photography from this mission will add to the global change databases and will contribute to the better understanding of our home planet. The use of manned space flights in understanding the global process first hand is a vital component in NASA's Mission to Planet Earth. The following are discussed along with photographs from the mission: landforms and geologic observation; environmental observations; meteorological/atmospheric observation; and oceanographic observations.

Description: The advantages of the astronaut photography during Space Shuttle missions are briefly examined, and the scope and applications of the Space Shuttle earth observations photography database are discussed. The global and multidisciplinary nature of the data base is illustrated by several examples of geologic applications. These include the eruption of Mount Pinatubo (Philippine Islands), heat flow and ice cover on Lake Baikal in Siberia (Russia), and windblown dust in South America. It is noted that hand-held photography from the U.S. Space Shuttle provides unique remotely-sensed data for geologic applications because of the combination of varying perspectives, look angles, and illumination, and changing resolution resulting from different lenses and altitudes.

Description: A review of the imagery acquired during the STS 50 mission of the Space Shuttle is presented. The earth viewing photography from this flight includes photos of dust plumes over several portions of the Red Sea, Arabian Sea, Persian Gulf, the Mediterranean Sea, and the Atlantic Ocean. Over land, prominent dust plumes were seen over Iraq, North Africa, Sudan, and West Africa. The color infrared photography includes images of the tropical rain forests of South America and South and Southeast Asia. Other examples include photographs of floods in Argentina, photos of Lake Chad in Africa, Coastal Madagascar, the Aswan dam and the Nile, geologic features of North Africa, the center pivot irrigation land areas of Saudi Arabia, flooding in Asian rivers, and sediment plumes of South American and South and Southeast Asian coasts.

Description: The paper presents a critical comment on the work of Henkes and Hoogendoorn. In particular, it is argued that the model of To and Humphrey (TH) (1986) has been implemented incorrectly, rendering invalid the conclusion that the TH model considerably deviates from the other models evaluated. In their response to the criticism, the authors present arguments in support of the correctness of their implementation of the TH model and of the validity of the conclusions reached in their work.

Description: The present extention of a Lagrangian approach of the Riemann solution procedure, which was originally proposed for perfect gases, to real gases, is nontrivial and requires the development of an exact real-gas Riemann solver for the Lagrangian form of the conservation laws. Calculations including complex wave interactions of various types were conducted to test the accuracy and robustness of the approach. Attention is given to the case of 2D oblique waves' capture, where a slip line is clearly in evidence; the real gas effect is demonstrated in the case of a generic engine nozzle.

Description: Implicit difference operator spectra are presently computed by applying eigensystem analysis techniques to finite-difference formulations of 2D Euler and Navier-Stokes equations, and attention is given to these iterative methods' convergence and stability characteristics by taking into account the effects of grid geometry, time-step, numerical viscosity, and boundary conditions. On the basis of the eigenvalue distributions for various flow configurations, the feasibility of applying such convergence-acceleration techniques as eigenvalue annihilation and relaxation is discussed. Spectrum-shifting is applied to NASA-Ames' ARC2D flow code, achieving a 20-33 percent efficiency.

Description: An Automated Bidirectional Reflection Acquisition Measurement System (ABRAMS) has been constructed to facilitate measurement of bidirectional reflectance from soil and vegetative samples in the laboratory. The system illuminates a sample with linearly-polarized laser light, lambda = 632.8 nm, and measures the like- and cross-polarized scattered intensities over half a hemisphere. System design and polarized bidirectional reflectance measurements from a soil sample and SiO2 spherical particles are discussed in this work. It is shown that polarization information in the plane of incidence is useful for identifying certain scattering mechanisms associated with soil reflectance. This is because the like-polarized intensity, Ivv, is influenced by single-scattered light and the cross-polarized intensity, IHv, is strongly influenced by multiple-scattered light. For example, comparable levels of Ivv and IHv indicated that the reflectance of soils, is dominated by significant multiple scattering because single scattering causes minimal depolarization in the plane of incidence.

Description: Study of turbulent flows in rotating reference frames has long been an area of considerable scientific and engineering interest. Because of its importance, the subject of turbulence in rotating reference frames has motivated over the years a large number of theoretical, experimental, and computational studies. The bulk of these previous works has served to demonstrate that the effect of system rotation on turbulence is subtle and remains exceedingly difficult to predict. A rotating flow of particular interest in many studies, including the present work, is examination of the effect of solid-body rotation on an initially isotropic turbulent flow. One of the principal reasons for the interest in this flow is that it represents the most basic turbulent flow whose structure is altered by system rotation but without the complicating effects introduced by mean strains or flow inhomogeneities. The assumption of statistical homogeneity considerably simplifies analysis and computation. The principal objective of the present study has been to examine the asymptotic state of solid-body rotation applied to an initially isotropic, high Reynolds number turbulent flow. Of particular interest has been to determine the degree of two-dimensionalization and the existence of asymptotic self-similar states in homogeneous rotating turbulence.

Description: The objective here is to test the Dynamic Localization (DL) model in a wall-bounded channel flow for numerical stability and accuracy of results. Algebraic stress models suggest that the model for the residual subgrid-scale (SGS) Reynolds stress and scalar flux should generally have terms comprising most of the unique products of the resolved strain (S) and rotation (R) tensors with S and the resolved scalar gradient. The standard dynamic SGS model uses a simple (Smagorinsky) base model for the residual Reynolds stress, which is made proportional to S, and down-gradient base models for residual scalar fluxes; these correspond to the lowest, 'first-order' terms in algebraic stress models. Temporal scaling terms in these base models are formed from the magnitude of the resolved strain rate. While this is appropriate for simple shear flows, it may not be appropriate for more complicated flows (relevant to geophysical and astrophysical problems) that include any combination of shear, rotation, buoyancy, etc. On the other hand, the coefficient in the dynamic SGS model readily adjusts itself to different flow conditions and may adequately take account of these effects without the need for more complicated base models. Cabot (1993) has begun to test the dynamic SGS model in buoyant flows (Rayleigh-Benard and internally heated convection) with and without buoyancy terms explicitly included in the scaling terms of the base model; no great differences were found in large eddy simulation (LES) results for the different base model scalings. The second objective in this work is to test base models with additional, 'second-order' terms (e.g., S(sup 2) and RS for the residual Reynolds stress). These terms have been found to improve large-scale flow predictions by kappa-epsilon models in the presence of rotation and shear. Second-order base models will be tested here in the LES of channel flow with and without solid-body rotation and compared with results from the standard first-order base models to determine if there are significant differences or improvements in results that would warrant the added complexity of the second-order base models.

Description: The present study is a continuation of the examination by Hamilton et al. of the regeneration mechanisms of near-wall turbulence and an attempt to investigate the conjecture of Waleffe et al. The basis of this study is an extension of the 'minimal channel' approach of Jimenez and Moin that emphasizes the near-wall region and reduces the complexity of the turbulent flow by considering a plane Couette flow of near minimum Reynolds number and stream-wise and span-wise extent. Reduction of the flow Reynolds number to the minimum value which will allow turbulence to be sustained has the effect of reducing the ratio of the largest scales to the smallest scales or, equivalently, of causing the near-wall region to fill more of the area between the channel walls. A plane Couette flow was chosen for study since this type of flow has a mean shear of a single sign, and at low Reynolds numbers, the two wall regions are found to share a single set of structures.

Description: Free shear flows, like those of mixing layers, are encountered in aerodynamics, in the atmosphere, and in the ocean as well as in many industrial applications such as flow reactors or combustion chambers. It is, therefore, crucial to understand the mechanisms governing the process of transition to turbulence in order to predict and control the evolution of the flow. Delaying transition to turbulence as far downstream as possible allows a gain in energy expenditure while accelerating the transition can be of interest in processes where high mixing is desired. Various methods, including the use of polymer additives, can be effective in controlling fluid flows. The drag reduction obtained by the addition of small amounts of high polymers has been an active area of research for the last three decades. It is now widely believed that polymer additives can affect the stability of a large variety of flows and that dilute solutions of these polymers have been shown to produce drag reductions of over 80 percent in internal flows and over 60 percent in external flows under a wide range of conditions. The major thrust of this work is to study the effects of polymer additives on the stability of the incompressible mixing layer through large scale numerical simulations. In particular, we focus on the two dimensional flow and examine how the presence of viscoelasticity may affect the typical structures of the flow, namely roll-up and pairing of vortices.

Description: Although the different regimes of premixed combustion are not well defined, most of the recent developments in turbulent combustion modeling are led in the so-called flamelet regime. The goal of these models is to give a realistic expression to the mean reaction rate (w). Several methods can be used to estimate (w). Bray and coworkers (Libby & Bray 1980, Bray 1985, Bray & Libby 1986) express the instantaneous reaction rate by means of a flamelet library and a frequency which describes the local interaction between the laminar flamelets and the turbulent flowfield. In another way, the mean reaction rate can be directly connected to the flame surface density (Sigma). This quantity can be given by the transport equation of the coherent flame model initially proposed by Marble & Broadwell 1977 and developed elsewhere. The mean reaction rate, (w), can also be estimated thanks to the evolution of an arbitrary scalar field G(x, t) = G(sub O) which represents the flame sheet. G(x, t) is obtained from the G-equation proposed by Williams 1985, Kerstein et al. 1988 and Peters 1993. Another possibility proposed in a recent study by Mantel & Borghi 1991, where a transport equation for the mean dissipation rate (epsilon(sub c)) of the progress variable c is used to determine (w). In their model, Mantel & Borghi 1991 considered a medium with constant density and constant diffusivity in the determination of the transport equation for (epsilon(sub c)). A comparison of different flamelet models made by Duclos et al. 1993 shows the realistic behavior of this model even in the case of constant density. Our objective in this present report is to present preliminary results on the study of this equation in the case of variable density and variable diffusivity. Assumptions of constant pressure and a Lewis number equal to unity allow us to significantly simplify the equation. A systematic order of magnitude analysis based on adequate scale relations is performed on each term of the equation. As in the case of constant density and constant diffusivity, the effects of stretching of the scalar field by the turbulent strain field, of local curvature, and of chemical reactions are predominant. In this preliminary work, we suggest closure models for certain terms, which will be validated after comparisons with DNS data.

Description: The recent development of the dynamic subgrid-scale (SGS) model has provided a consistent method for generating localized turbulent mixing models and has opened up great possibilities for applying the large eddy simulation (LES) technique to real world problems. Given the fact that the direct numerical simulation (DNS) can not solve for engineering flow problems in the foreseeable future (Reynolds 1989), the LES is certainly an attractive alternative. It seems only natural to bring this new development in SGS modeling to bear on the reacting flows. The major stumbling block for introducing LES to reacting flow problems has been the proper modeling of the reaction source terms. Various models have been proposed, but none of them has a wide range of applicability. For example, some of the models in combustion have been based on the flamelet assumption which is only valid for relatively fast reactions. Some other models have neglected the effects of chemical reactions on the turbulent mixing time scale, which is certainly not valid for fast and non-isothermal reactions. The probability density function (PDF) method can be usefully employed to deal with the modeling of the reaction source terms. In order to fit into the framework of LES, a new PDF, the large eddy PDF (LEPDF), is introduced. This PDF provides an accurate representation for the filtered chemical source terms and can be readily calculated in the simulations. The details of this scheme are described.

Description: Application of the method of Large Eddy Simulation (LES) to a turbulent flow consists of three separate steps. First, a filtering operation is performed on the Navier-Stokes equations to remove the small spatial scales. The resulting equations that describe the space time evolution of the 'large eddies' contain the subgrid-scale (sgs) stress tensor that describes the effect of the unresolved small scales on the resolved scales. The second step is the replacement of the sgs stress tensor by some expression involving the large scales - this is the problem of 'subgrid-scale modeling'. The final step is the numerical simulation of the resulting 'closed' equations for the large scale fields on a grid small enough to resolve the smallest of the large eddies, but still much larger than the fine scale structures at the Kolmogorov length. In dividing a turbulent flow field into 'large' and 'small' eddies, one presumes that a cut-off length delta can be sensibly chosen such that all fluctuations on a scale larger than delta are 'large eddies' and the remainder constitute the 'small scale' fluctuations. Typically, delta would be a length scale characterizing the smallest structures of interest in the flow. In an inhomogeneous flow, the 'sensible choice' for delta may vary significantly over the flow domain. For example, in a wall bounded turbulent flow, most statistical averages of interest vary much more rapidly with position near the wall than far away from it. Further, there are dynamically important organized structures near the wall on a scale much smaller than the boundary layer thickness. Therefore, the minimum size of eddies that need to be resolved is smaller near the wall. In general, for the LES of inhomogeneous flows, the width of the filtering kernel delta must be considered to be a function of position. If a filtering operation with a nonuniform filter width is performed on the Navier-Stokes equations, one does not in general get the standard large eddy equations. The complication is caused by the fact that a filtering operation with a nonuniform filter width in general does not commute with the operation of differentiation. This is one of the issues that we have looked at in detail as it is basic to any attempt at applying LES to complex geometry flows. Our principal findings are summarized.

Description: An unified approach for solving both compressible and incompressible flows was investigated in this study. The difference in CFD code development between incompressible and compressible flows is due to the mathematical characteristics. However, if one can modify the continuity equation for incompressible flows by introducing pseudocompressibility, the governing equations for incompressible flows would have the same mathematical characters as compressible flows. The application of a compressible flow code to solve incompressible flows becomes feasible. Among numerical algorithms developed for compressible flows, the Centered Total Variation Diminishing (CTVD) schemes possess better mathematical properties to damp out the spurious oscillations while providing high-order accuracy for high speed flows. It leads us to believe that CTVD schemes can equally well solve incompressible flows. In this study, the governing equations for incompressible flows include the continuity equation and momentum equations. The continuity equation is modified by adding a time-derivative of the pressure term containing the artificial compressibility. The modified continuity equation together with the unsteady momentum equations forms a hyperbolic-parabolic type of time-dependent system of equations. The continuity equation is modified by adding a time-derivative of the pressure term containing the artificial compressibility. The modified continuity equation together with the unsteady momentum equations forms a hyperbolic-parabolic type of time-dependent system of equations. Thus, the CTVD schemes can be implemented. In addition, the boundary conditions including physical and numerical boundary conditions must be properly specified to obtain accurate solution. The CFD code for this research is currently in progress. Flow past a circular cylinder will be used for numerical experiments to determine the accuracy and efficiency of the code before applying this code to more specific applications.

Description: We describe a multigrid multiblock method for compressible turbulent flow simulations and present results obtained from calculations on a two-element airfoil. A vertex-based spatial discretization method and explicit multistage Runge-Kutta time-stepping are used. The slow convergence of a single grid method makes the multigrid method, which yields a speed factor of about 20, indispensable. The numerical predictions are in good agreement with experimental results. It is shown that the convergence of the multigrid process depends considerably on the ordering of the various loops. If the block loop is put inside the stage loop, the process converges more rapidly than if the block loop is situated outside the stage loop when a three-stage Runge-Kutta method is used. If a five-stage scheme is used, the process does not converge in the latter block ordering. Finally, the process based on the five-stage method is about 60 percent more efficient than with the three-stage method, if the block loop is inside the stage loop.

Description: Grid staggering is a well known remedy for the problem of velocity/pressure coupling in incompressible flow calculations. Numerous inconveniences occur, however, when staggered grids are implemented, particularly when a general-purpose code, capable of handling irregular three-dimensional domains, is sought. In several non-staggered grid numerical procedures proposed in the literature, the velocity/pressure coupling is achieved by either pressure or velocity (momentum) averaging. This approach is not convenient for simultaneous (block) solvers that are preferred when using multigrid methods. A new method is introduced in this paper that is based upon non-staggered grid formulation with a set of virtual cell face velocities used for pressure/velocity coupling. Instead of pressure or velocity averaging, a momentum balance at the cell face is used as a link between the momentum and mass balance constraints. The numerical stencil is limited to 9 nodes (in 2D) or 27 nodes (in 3D), both during the smoothing and inter-grid transfer, which is a convenient feature when a block point solver is applied. The results for a lid-driven cavity and a cube in a lid-driven cavity are presented and compared to staggered grid calculations using the same multigrid algorithm. The method is shown to be stable and produce a smooth (wiggle-free) pressure field.

Description: The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) is a facility instrument selected for launch in 1998 on the first in a series of spacecraft for NASA's Earth Observing System (EOS). The ASTER instrument is being sponsored and built in Japan. It is a three telescope, high spatial resolution imaging instrument with 15 spectral bands covering the visible through to the thermal infrared. It will play a significant role within EOS providing geological, biological, land hydrological information necessary for intense study of the Earth. The operational capabilities for ASTER, including the necessary interfaces and operational collaborations between the US and Japanese participants, are under development. EOS operations are the responsibility of the EOS Project at NASA's Goddard Space Flight Center (GSFC). Although the primary EOS control center is at GSFC, the ASTER control facility will be in Japan. Other aspects of ASTER are discussed.

Description: The study of a thin, incompressible Newtonian fluid layer trapped between two almost parallel, sliding surfaces has been actively pursued in the last decades. This subject includes lubrication applications such as slider bearings or the sealing of non-pressurized fluids with rubber rotary shaft seals. In the present work we analyze numerically the flow of lubricant fluid through a micro-gap of sealing devices. The first stage of this study is carried out assuming that a 'small-gap' parameter delta attains an extreme value in the Navier-Stokes equations. The precise meaning of small-gap is achieved by the particular limit delta = 0 which, within the bounds of the hypotheses, predicts transport of lubricant through the sealed area by centrifugal instabilities. Numerical results obtained with the penalty function approximation in the finite element method are presented. In particular, the influence of inflow and outflow boundary conditions, and their impact in the simulated flow, are discussed.

Description: The development of improved performance models for the Space Shuttle Main Engine (SSME) is an important, ongoing program at NASA MSFC. These models allow prediction of overall system performance, as well as analysis of run-time anomalies which might adversely affect engine performance or safety. Due to the complexity of the flow fields associated with the SSME, NASA has increasingly turned to Computational Fluid Dynamics (CFD) techniques as modeling tools. An important component of the SSME system is the fuel preburner, which consists of a cylindrical chamber with a plate containing 264 coaxial injector elements at one end. A fuel rich mixture of gaseous hydrogen and liquid oxygen is injected and combusted in the chamber. This process preheats the hydrogen fuel before it enters the main combustion chamber, powers the hydrogen turbo-pump, and provides a heat dump for nozzle cooling. Issues of interest include the temperature and pressure fields at the turbine inlet and the thermal compatibility between the preburner chamber and injector plate. Performance anomalies can occur due to incomplete combustion, blocked injector ports, etc. The performance model should include the capability to simulate the effects of these anomalies. The current approach to the numerical simulation of the SSME fuel preburner flow field is to use a global model based on the MSFC sponsored FNDS code. This code does not have the capabilities of modeling several aspects of the problem such as detailed modeling of the coaxial injectors. Therefore, an effort has been initiated to develop a detailed simulation of the preburner coaxial injectors and provide gas phase boundary conditions just downstream of the injector face as input to the FDNS code. This simulation should include three-dimensional geometric effects such as proximity of injectors to baffles and chamber walls and interaction between injectors. This report describes an investigation into the numerical simulation of GH2/LOX coaxial injectors. The following sections will discuss the physical aspects of injectors, the CFD code employed, and preliminary results of a simulation of a single coaxial injector for which experimental data is available. It is hoped that this work will lay the foundation for the development of a unique and useful tool to support the SSME program.

Description: As the workstation and personal computer become more popular than a centralized mainframe to perform thermal analysis, the methods for space vehicle thermal analysis will change. Already, many thermal analysis codes are now available for workstations, which were not in existence just five years ago. As these changes occur, some organizations will adopt the new codes and analysis techniques, while others will not. This might lead to misunderstandings between thermal shops in different organizations. If thermal analysts make an effort to understand the major differences between the new and old methods, a smoother transition to a more efficient and more versatile thermal analysis environment will be realized.

Description: Thermal loads on current and future aircraft are increasing and as a result are stressing the energy collection, control, and dissipation capabilities of current thermal management systems and technology. The thermal loads for hypersonic vehicles will be no exception. In fact, with their projected high heat loads and fluxes, hypersonic vehicles are a prime example of systems that will require thermal management systems (TMS) that have been optimized and integrated with the entire vehicle to the maximum extent possible during the initial design stages. This will not only be to meet operational requirements, but also to fulfill weight and performance constraints in order for the vehicle to takeoff and complete its mission successfully. To meet this challenge, the TMS can no longer be two or more entirely independent systems, nor can thermal management be an after thought in the design process, the typical pervasive approach in the past. Instead, a TMS that was integrated throughout the entire vehicle and subsequently optimized will be required. To accomplish this, a method that iteratively optimizes the TMS throughout the vehicle will not only be highly desirable, but advantageous in order to reduce the manhours normally required to conduct the necessary tradeoff studies and comparisons. A thermal management engineering computer code that is under development and being managed at Wright Laboratory, Wright-Patterson AFB, is discussed. The primary goal of the code is to aid in the development of a hypersonic vehicle TMS that has been optimized and integrated on a total vehicle basis.