Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distribution grid”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Wall-Modeled Large Eddy Simulation Method for Unstructured-Grid Navier-Stokes Solvers

This paper reports on the implementation and assessment of a Wall-Modeled Large-Eddy Simulation (WMLES) methodology in an unstructured-grid, node-centered flow solver, FUN3D that is developed and supported at the NASA Langley Research Center. Finite-volume (FV) and finite-element (FE) discretization schemes considered in the study provide formal second-order spatial accuracy. Large-Eddy Simulations (LES) resolve large-scale turbulent-flow features and filter out small-scale effects using the Vreman subgrid-scale model. At solid-wall boundaries, a shear-stress model is employed to provide a proper boundary-flux closure. The nonlinear equations are integrated in time using either an optimized backward difference formula or an implicit multistage Runge-Kutta temporal scheme. The implicit equations at each time step are solved by strong nonlinear iteration schemes. WMLES demonstrations are shown for two high-lift configurations, namely, the McDonnell Douglas 30P30N multielement airfoil and a NASA High-Lift Common Research Model. Results show that the WMLES approaches implemented in the FV and FE discretization methods produce consistent solutions and are capable of capturing key aerodynamic characteristics and flow structures for high-lift configurations at a wide range of angles of attack including maximum-lift conditions. In the 30P30N example, correct trends in the variations of integrated aerodynamic forces and moments, surface pressure distributions, and boundary-layer profiles are captured as the Reynolds number is increased.

CFD; turbulence modeling; High-Lift flow simulatio↗
NASA NTRS2008

Application of the CloudSat and NEXRAD Radars Toward Improvements in High Resolution Operational Forecasts

As computational power increases, operational forecast models are performing simulations with higher spatial resolution allowing for the transition from sub-grid scale cloud parameterizations to an explicit forecast of cloud characteristics and precipitation through the use of single- or multi-moment bulk water microphysics schemes. investments in space-borne and terrestrial remote sensing have developed the NASA CloudSat Cloud Profiling Radar and the NOAA National Weather Service NEXRAD system, each providing observations related to the bulk properties of clouds and precipitation through measurements of reflectivity. CloudSat and NEXRAD system radars observed light to moderate snowfall in association with a cold-season, midlatitude cyclone traversing the Central United States in February 2007. These systems are responsible for widespread cloud cover and various types of precipitation, are of economic consequence, and pose a challenge to operational forecasters. This event is simulated with the Weather Research and Forecast (WRF) Model, utilizing the NASA Goddard Cumulus Ensemble microphysics scheme. Comparisons are made between WRF-simulated and observed reflectivity available from the CloudSat and NEXRAD systems. The application of CloudSat reflectivity is made possible through the QuickBeam radiative transfer model, with cautious application applied in light of single scattering characteristics and spherical target assumptions. Significant differences are noted within modeled and observed cloud profiles, based upon simulated reflectivity, and modifications to the single-moment scheme are tested through a supplemental WRF forecast that incorporates a temperature dependent snow crystal size distribution.

Molthan, A. L.↗
NASA NTRSDate not supplied

Are You In or Out: the Influence of Sea Ice Drift on Sea Ice Survivability

The state of Arctic sea ice influences aspects of many global systems including ecosystems, economies, geopolitics, and climate. And, it is in severe decline. Numerous studies have considered the factors that influence Arctic sea ice decline using monthly mean or gridded data. However, these methods struggle to cleanly separate the influences on sea ice melt because the sea ice cover is always moving; thus, different sea ice floes and sea ice regimes constitute the monthly average sea ice state within the Eulerian frame. This study uses a LaGrangian sea ice parcel tracking satellite database to investigate the factors that influence sea ice parcel survivability. The concept of sea ice survivability is defined as the likelihood that a sea ice parcel will last through the summer melt season. Survivability can be stratified in many ways to analyze how it changes as a function of sea ice regime, region, and sea ice characteristics. How is the survivability of first year (FY) and multi-year (MY) ice classes influenced by sea ice motion? This is the central question addressed in this presentation. Specifically, this analysis determines sea ice survivability as a function of distance traveled by sea ice parcels and compares the survivability statistics for sea ice parcels that stay in their starting region versus those that leave their starting region. The preliminary results indicate that the influence of sea ice movement on survivability depends strongly on which region the sea ice resides in at the beginning of the growth season. Lastly, we find substantial inter-annual variability in the regional distribution of sea ice survivability that seems to be strongly linked to sea ice parcel drift between regions of the Arctic Ocean. Moreover, we use these results to discuss the potential contribution of changes in sea ice motion on the observed slowdown in Arctic sea ice extent loss.

Patrick C Taylor↗
NASA NTRS2006

Introduction to MODIS Cloud Products

The Earth's radiative energy balance and hydrological cycle are fundamentally coupled with the distribution and properties of clouds. Therefore, the ability to remotely infer cloud properties and their variation in space and time is crucial for establishing climatologies as a reference for validation of present-day climate models and in assessing future climate change. Remote cloud observations also provide data sets useful for testing and improving cloud model physics, and for assimilation into numerical weather prediction models. The MODerate Resolution Imaging Spectroradiometer (MODIS) imagers on the Terra and Aqua Earth Observing System (EOS) platforms provide the capability for globally retrieving these properties using passive solar reflectance and infrared techniques. In addition to providing measurements similar to those offered on a suite of historical operational weather platforms such as the Advanced Very High Resolution Radiometer (AVHRR), the High-resolution Infrared Radiation Sounder (HIRS), and the Geostationary Operational Environmental Satellite (GOES), MODIS provides additional spectral and/or spatial resolution in key atmospheric bands, along with on-board calibration, to expand the capability for global cloud property retrievals. The core MODIS operational cloud products include cloud top pressure, thermodynamic phase, optical thickness, particle size, and water path, and are derived globally at spatial resolutions of either 1- or 5-km (referred to as Level-2 or pixel-level products). In addition, the MODIS atmosphere team (collectively providing cloud, aerosol, and clear sky products) produces a combined gridded product (referred to as Level-3) aggregated to a 1 equal-angle grid, available for daily, eight-day, and monthly time periods. The wealth of information available from these products provides critical information for climate studies as well as the continuation and improved understanding of existing satellite-based cloud climatologies derived from heritage instruments. This chapter provides an overview of the MODIS Level-2 and -3 operational cloud products.

Baum, Bryan A.↗
NASA NTRSDate not supplied

Cloud-Type Mean Cloud Properties and Associated Cloud Radiative Effects in the Tropical “Chimney” Zones Using 20 Years High-Resolution CERES Satellite Data

It has been well known that convection tends to be more intense over land than over ocean and continental convection generally contains wider cores that are protected from entrainment than their oceanic counterparts. How does this difference in convective intensity impact the rest of cloud types and associated cloud radiative effects (CREs)? We select five regions in the tropical “chimney” zones, with two regions over land (Africa and Amazon) and three regions over ocean (eastern and western Pacific, and Atlantic) to understand the differences in the cloud-type mean cloud properties and associated CREs using 20-years high-resolution (2x2 km2) CERES satellite data. The cloud types are based upon the ISCCP classification using the joint cloud-top pressure and cloud optical depth distribution. These five regions are compared to the entire tropics (25°S to 25°N). First, we compare the frequencies of occurrence of cloud types among the regions. Compared to the entire tropics, the “chimney” zones are cloudier, particularly, for cloud types with moderate cloud optical depths in the middle and upper troposphere. Except for the eastern Pacific region, the low-level cloud types are less abundant over land (due to higher boundary layers) and ocean (due to weaker lower-tropospheric subsidence). Over Africa, upper-level clouds are scarcer, due perhaps to weak production of convective anvils related to the drier atmosphere. Second, liquid water path (LWP) and ice water path (IWP) are compared among the regions. Compared to the entire tropics, mean LWPs for moderate cloud optical depths are higher for oceanic regions but lower for Africa while they are slightly lower in the lower troposphere but slightly higher in the middle troposphere for Amazon. Not surprisingly, mean IWPs are higher over land than over ocean, which is directly related to the difference in convective intensity. That is, the explanations to this result are the relatively weak convective intensity over the eastern Pacific and the abundance of anvil clouds over the western Pacific. Both reduce the mean IWPs. The CRE results are still being analyzed because the totally-clear grids (1° x 1°) of the CERES data product should be removed before taking a regional mean of clear-sky fluxes, which cause the wrong signs of CREs for a few cloud types.

cloud radiative effects↗
NASA NTRS1996

Development of Multiobjective Optimization Techniques for Sonic Boom Minimization

A discrete, semi-analytical sensitivity analysis procedure has been developed for calculating aerodynamic design sensitivities. The sensitivities of the flow variables and the grid coordinates are numerically calculated using direct differentiation of the respective discretized governing equations. The sensitivity analysis techniques are adapted within a parabolized Navier Stokes equations solver. Aerodynamic design sensitivities for high speed wing-body configurations are calculated using the semi-analytical sensitivity analysis procedures. Representative results obtained compare well with those obtained using the finite difference approach and establish the computational efficiency and accuracy of the semi-analytical procedures. Multidisciplinary design optimization procedures have been developed for aerospace applications namely, gas turbine blades and high speed wing-body configurations. In complex applications, the coupled optimization problems are decomposed into sublevels using multilevel decomposition techniques. In cases with multiple objective functions, formal multiobjective formulation such as the Kreisselmeier-Steinhauser function approach and the modified global criteria approach have been used. Nonlinear programming techniques for continuous design variables and a hybrid optimization technique, based on a simulated annealing algorithm, for discrete design variables have been used for solving the optimization problems. The optimization procedure for gas turbine blades improves the aerodynamic and heat transfer characteristics of the blades. The two-dimensional, blade-to-blade aerodynamic analysis is performed using a panel code. The blade heat transfer analysis is performed using an in-house developed finite element procedure. The optimization procedure yields blade shapes with significantly improved velocity and temperature distributions. The multidisciplinary design optimization procedures for high speed wing-body configurations simultaneously improve the aerodynamic, the sonic boom and the structural characteristics of the aircraft. The flow solution is obtained using a comprehensive parabolized Navier Stokes solver. Sonic boom analysis is performed using an extrapolation procedure. The aircraft wing load carrying member is modeled as either an isotropic or a composite box beam. The isotropic box beam is analyzed using thin wall theory. The composite box beam is analyzed using a finite element procedure. The developed optimization procedures yield significant improvements in all the performance criteria and provide interesting design trade-offs. The semi-analytical sensitivity analysis techniques offer significant computational savings and allow the use of comprehensive analysis procedures within design optimization studies.

Chattopadhyay, Aditi↗
NASA NTRS2005

NSTAR Extended Life Test Discharge Chamber Flake Analysis

The Extended Life Test (ELT) of the NASA Solar Electric Propulsion Technology Readiness (NSTAR) ion thruster was concluded after 30,352 hours of operation. The ELT was conducted using the Deep Space 1 (DS1) back-up flight engine, a 30 cm diameter xenon ion thruster. Post-test inspection of the ELT engine revealed numerous contaminant flakes distributed over the bottom of the cylindrical section of the anode within the discharge chamber (DC). Extensive analyses were conducted to determine the source of the particles, which is critical to the understanding of degradation mechanisms of long life ion thruster operation. Analyses included: optical microscopy (OM) and particle length histograms, field emission scanning electron microscopy (FESEM) combined with energy dispersive spectroscopy (EDS), and atomic oxygen plasma exposure tests. Analyses of the particles indicate that the majority of the DC flakes consist of a layered structure, typically with either two or three layers. The flakes comprising two layers were typically found to have a molybdenum-rich (Mo-rich) layer on one side and a carbon-rich (C-rich) layer on the other side. The flakes comprising three layers were found to be sandwich-like structures with Mo-rich exterior layers and a C-rich interior layer. The presence of the C-rich layers indicates that these particles were produced by sputter deposition build-up on a surface external to the discharge chamber from ion sputter erosion of the graphite target in the test chamber. This contaminant layer became thick enough that particles spalled off, and then were electro-statically attracted into the ion thruster interior, where they were coated with Mo from internal sputter erosion of the screen grid and cathode components. Atomic oxygen tests provided evidence that the DC chamber flakes are composed of a significant fraction of carbon. Particle size histograms further indicated that the source of the particles was spalling of carbon flakes from downstream surfaces. Analyses of flakes taken from the downstream surface of the accelerator grid provided additional supportive information. The production of the downstream carbon flakes, and hence the potential problems associated with the flake particles in the ELT ion thruster engine is a facility induced effect and would not occur in the space environment.

deGroh, Kim K.↗
NASA NTRS1999

Computations of Boiling in Microgravity

The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work is limited to very simple models. In this project, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. The proposed work is based on previously funded NASA work that allowed us to develop a two-dimensional numerical method for boiling flows and to demonstrate the ability of the method to simulate film boiling. While numerical simulations of multi-fluid flows have been advanced in a major way during the last five years, or so, similar capability for flows with phase change are still in their infancy. Although the feasibility of the proposed approach has been demonstrated, it has yet to be extended and applied to fully three-dimensional simulations. Here, a fully three-dimensional, parallel, grid adaptive code will be developed. The numerical method will be used to study nucleate boiling in microgravity, with particular emphasis on two aspects of the problem: 1) Examination of the growth of bubbles at a wall nucleation site and the instabilities of rapidly growing bubbles. Particular emphasis will be put on accurately capturing the thin wall layer left behind as a bubble expands along a wall, on computing instabilities on bubble surfaces as bubbles grow, and on quantifying the effects of both these phenomena on heat transfer; and 2) Examination of the effect of shear flow on bubble growth and heat transfer.

Tryggvason, Gretar↗
NASA NTRS2007

Comparisons of Instantaneous TRMM Ground Validation and Satellite Rain Rate Estimates at Different Spatial Scales

This study provides a comprehensive inter-comparison of instantaneous rain estimates from the two rain sensors aboard the TRMM satellite with ground data from thee designated Ground Validation Sites: Kwajalein Atoll, Melbourne, Florida and Houston, Texas. The satellite rain retrievals utilize rain observations collected by the TRMM microwave imager (TMI) and the Precipitation Radar (PR) aboard the TRMM satellite. Three standard instantaneous rain products are the generated from the rain information retrieved from the satellite using the TMI, PR and Combined (COM) rain algorithms. The validation data set used in this study was obtained from instantaneous rain rates inferred from ground radars at each GV site. The first comparison used 0.5(sup 0) x 0.5(sup 0) gridded data obtained from the TRMM 3668 product, and similarly gridded GV data obtained from ground-based radars. The comparisons were made at the same spatial and temporal scales in order to eliminate sampling biases in our comparisons. An additional comparison was made by averaging rain rates for the PR, COM and GV estimates within each TMI footprint (approx. 150 square kilometers). For this analysis, unconditional mean rain rates from PR, COM and GV estimates were calculated within each TMI footprint that was observed within 100 km from the respective GV site (and also observed by the PR). This analysis used all the available matching data from the period 1999-2004, representing a sample size of over 50,000 footprints for each site. In the first analysis our results showed that all of the respective rain rate estimates agree well, with some exceptions. The more salient differences were associated with heavy rain events in which one or more of the algorithms failed to properly retrieve these extreme events. Also, it appears that there is a preferred mode of precipitation for TMI rain rates at or near 2 mm per hour over the ocean. This mode was noted over ocean areas of Melbourne, Florida and Kwajalein, Republic of the Marshall Islands, and is shown to exist in TRMM tropical-global ocean areas as well. Further research by algorithm developers is needed to explain or justify the seemingly errant observed probability distributions.

Wolff, David B.↗
NASA NTRS1994

Applications Performance on NAS Intel Paragon XP/S - 15#

The Numerical Aerodynamic Simulation (NAS) Systems Division received an Intel Touchstone Sigma prototype model Paragon XP/S- 15 in February, 1993. The i860 XP microprocessor with an integrated floating point unit and operating in dual -instruction mode gives peak performance of 75 million floating point operations (NIFLOPS) per second for 64 bit floating point arithmetic. It is used in the Paragon XP/S-15 which has been installed at NAS, NASA Ames Research Center. The NAS Paragon has 208 nodes and its peak performance is 15.6 GFLOPS. Here, we will report on early experience using the Paragon XP/S- 15. We have tested its performance using both kernels and applications of interest to NAS. We have measured the performance of BLAS 1, 2 and 3 both assembly-coded and Fortran coded on NAS Paragon XP/S- 15. Furthermore, we have investigated the performance of a single node one-dimensional FFT, a distributed two-dimensional FFT and a distributed three-dimensional FFT Finally, we measured the performance of NAS Parallel Benchmarks (NPB) on the Paragon and compare it with the performance obtained on other highly parallel machines, such as CM-5, CRAY T3D, IBM SP I, etc. In particular, we investigated the following issues, which can strongly affect the performance of the Paragon: a. Impact of the operating system: Intel currently uses as a default an operating system OSF/1 AD from the Open Software Foundation. The paging of Open Software Foundation (OSF) server at 22 MB to make more memory available for the application degrades the performance. We found that when the limit of 26 NIB per node out of 32 MB available is reached, the application is paged out of main memory using virtual memory. When the application starts paging, the performance is considerably reduced. We found that dynamic memory allocation can help applications performance under certain circumstances. b. Impact of data cache on the i860/XP: We measured the performance of the BLAS both assembly coded and Fortran coded. We found that the measured performance of assembly-coded BLAS is much less than what memory bandwidth limitation would predict. The influence of data cache on different sizes of vectors is also investigated using one-dimensional FFTs. c. Impact of processor layout: There are several different ways processors can be laid out within the two-dimensional grid of processors on the Paragon. We have used the FFT example to investigate performance differences based on processors layout.

Saini, Subhash↗
NASA NTRS2015

A Framework for Orbital Performance Evaluation in Distributed Space Missions for Earth Observation

Distributed Space Missions (DSMs) are gaining momentum in their application to earth science missions owing to their unique ability to increase observation sampling in spatial, spectral and temporal dimensions simultaneously. DSM architectures have a large number of design variables and since they are expected to increase mission flexibility, scalability, evolvability and robustness, their design is a complex problem with many variables and objectives affecting performance. There are very few open-access tools available to explore the tradespace of variables which allow performance assessment and are easy to plug into science goals, and therefore select the most optimal design. This paper presents a software tool developed on the MATLAB engine interfacing with STK, for DSM orbit design and selection. It is capable of generating thousands of homogeneous constellation or formation flight architectures based on pre-defined design variable ranges and sizing those architectures in terms of predefined performance metrics. The metrics can be input into observing system simulation experiments, as available from the science teams, allowing dynamic coupling of science and engineering designs. Design variables include but are not restricted to constellation type, formation flight type, FOV of instrument, altitude and inclination of chief orbits, differential orbital elements, leader satellites, latitudes or regions of interest, planes and satellite numbers. Intermediate performance metrics include angular coverage, number of accesses, revisit coverage, access deterioration over time at every point of the Earth's grid. The orbit design process can be streamlined and variables more bounded along the way, owing to the availability of low fidelity and low complexity models such as corrected HCW equations up to high precision STK models with J2 and drag. The tool can thus help any scientist or program manager select pre-Phase A, Pareto optimal DSM designs for a variety of science goals without having to delve into the details of the engineering design process.

Earth Science↗
NASA NTRS2009

Cold Season QPF: Sensitivities to Snow Parameterizations and Comparisons to NASA CloudSat Observations

As operational numerical weather prediction is performed at increasingly finer spatial resolution, precipitation traditionally represented by sub-grid scale parameterization schemes is now being calculated explicitly through the use of single- or multi-moment, bulk water microphysics schemes. As computational resources grow, the real-time application of these schemes is becoming available to a broader audience, ranging from national meteorological centers to their component forecast offices. A need for improved quantitative precipitation forecasts has been highlighted by the United States Weather Research Program, which advised that gains in forecasting skill will draw upon improved simulations of clouds and cloud microphysical processes. Investments in space-borne remote sensing have produced the NASA A-Train of polar orbiting satellites, specially equipped to observe and catalog cloud properties. The NASA CloudSat instrument, a recent addition to the A-Train and the first 94 GHz radar system operated in space, provides a unique opportunity to compare observed cloud profiles to their modeled counterparts. Comparisons are available through the use of a radiative transfer model (QuickBeam), which simulates 94 GHz radar returns based on the microphysics of cloudy model profiles and the prescribed characteristics of their constituent hydrometeor classes. CloudSat observations of snowfall are presented for a case in the central United States, with comparisons made to precipitating clouds as simulated by the Weather Research and Forecasting Model and the Goddard single-moment microphysics scheme. An additional forecast cycle is performed with a temperature-based parameterization of the snow distribution slope parameter, with comparisons to CloudSat observations provided through the QuickBeam simulator.

Molthan, A. L.↗
NASA NTRS1993

Ionospheric calibration for single frequency altimeter measurements

This report investigates the potential of using Global Positioning System (GPS) data and a model of the ionosphere to supply a measure of the sub-satellite Total Electron Current (TEC) of the required accuracy (10 TECU rms) for the purpose of calibrating single frequency radar altimeter measurements. Since climatological (monthly mean) models are known to be in error by as much as 50 percent, this work focused on the Parameterized Real-Time Ionospheric Specification Model (PRISM) which has the capability to improve model accuracy by ingesting (adjusting to) in situ ionospheric measurements. A set of globally distributed TEC measurements were generated using GPS data and were used as input to improve the accuracy of the PRISM model. The adjusted PRISM TEC values were compared to TOPEX dual frequency TEC measurements (which are considered truth) for a number of TOPEX sub-satellite tracks. The adjusted PRISM values generally compared to the TOPEX measurements within the 10 TECU accuracy requirements when the sub-satellite track passed within 300 to 400 km of the GPS TEC data or when the track passed through a night time ionosphere. However, when the sub-satellite points were greater than 300 to 400 km away from the GPS TEC data or when a local noon ionosphere was sampled, the adjusted PRISM values generally differed by greater than 10 TECU rms with data excursions from the TOPEX TEC measurements of as much as 40 TECU (an 8 cm path delay error at K band). Therefore, it can be concluded from this analysis that an unrealistically large number of GPS stations would be needed to predict sub-satellite TEC at the 10 TECU level in the day time ionosphere using a model such as PRISM. However, a technique currently being studied at the Jet Propulsion Laboratory (JPL) may provide a means of supplying adequate TEC data to meet the 10 TECU ionospheric correction accuracy when using a realistic number of ionospheric stations. This method involves using global GPS TEC data to estimate a global grid of vertical ionospheric TEC as a function of time (i.e. every one half hour) in a sun-fixed longitude frame. Working in a sun-fixed longitude frame, one is not limited by the spatial decorrelation distance of the ionosphere, but instead is limited more by the temporal correlations of the ionosphere in the sun-fixed frame which are a smaller effect. It is the opinion of the authors that using the global sun-fixed TEC grid data, in particular, ingesting it into PRISM, offers the best possibility of meeting the 10 TECU ionospheric correction accuracy requirement, and should be the subject of further study.

Schreiner, William S.↗
NASA NTRS2021

Maximum entropy modeling to identify physical drivers of shallow snowpack heterogeneity using unpiloted aerial system (UAS) lidar

Understanding the spatial variability of the snowpack is valuable for hydrologists and ecologists seeking to predict hydrological processes in a cold region. Snow distribution is a function of interactions among static variables, such as terrain, vegetation, and soil properties, and dynamic meteorological variables, such as solar radiation, wind speed and direction, and soil moisture. However, identifying the dominant physical drivers responsible for spatial patterns of the snowpack, particularly for ephemeral, shallow snowpacks, has been challenging due to the lack of the high-resolution snowpack and physical variables with high vertical accuracy as well as inherent limitations in traditional approaches. This study uses an Unpiloted Aerial System (UAS) lidar-based snow depth and static variables (1-m spatial resolution) to analyze field-scale spatial structures of snow depth and apply the Maximum Entropy (MaxEnt) model to identify primary controls over open terrain and forests at the University of New Hampshire Thompson Farm Research Observatory, New Hampshire, United States. We found that, among nine topographic and soil variables, plant functional type and terrain roughness contribute up to 80% and 76% of relative importance in the MaxEnt framework to predict locations of deeper or shallower snowpacks, respectively, across a mixed temperate forested and field landscape. Soil variables, such as organic matter and saturated hydraulic conductivity, were also important controls (up to 70% and 81%) on snow depth spatial variations for both open and forested landscapes suggesting spatial variations in soil variables under snow can control thermal transfer among soil, snowpack, and surface-atmosphere. This work contributes to improving land surface and snow models by informing parameterization of the sub-grid scale snow depths, down-scaling remotely sensed snow products, and understanding field scale snow states.

Eunsang Cho↗
NASA NTRS2000

The Geoscience Laser Altimeter System (GLAS) for the ICESAT Mission

The Laser In space Technology Experiment, Shuttle Laser Altimeter and the Mars Observer Laser Altimeter have demonstrated accurate measurements of atmospheric backscatter and Surface heights from space. The recent MOLA measurements of the Mars surface have 40 cm vertical resolution and have reduced the global uncertainty in Mars topography from a few km to about 5 m. The Geoscience Laser Altimeter System (GLAS) is a next generation lidar for Earth orbit being developed as part of NASA's Icesat Mission. The GLAS design combines a 10 cm precision surface lidar with a sensitive dual wavelength cloud and aerosol lidar. GLAS will precisely measure the heights of the Earth's polar ice sheets, establish a grid of accurate height profiles of the Earth's land topography, and profile the vertical backscatter of clouds and aerosols on a global scale. GLAS is being developed to fly on a small dedicated spacecraft in a polar orbit with a 590 630 km altitude at inclination of 94 degrees. GLAS is scheduled to launch in the summer 2001 and to operate continuously for a minimum of 3 years with a goal of 5 years. The primary mission for GLAS is to measure the seasonal and annual changes in the heights of the Greenland and Antarctic ice sheets. GLAS will continuously measure the vertical distance from orbit to the Earth's surface with 1064 nm pulses from a ND:YAG laser at a 40 Hz rate. Each 5 nsec wide laser pulse is used to produce a single range measurement, and the laser spots have 66 m diameter and about 170 m center-center spacings. When over land GLAS will profile the heights of the topography and vegetation. The GLAS receiver uses a 1 m diameter telescope and a Si APD detector. The detector signal is sampled by an all digital receiver which records each surface echo waveform with I nsec resolution and a stored echo record lengths of either 200, 400, or 600 samples. Analysis of the echo waveforms within the instrument permits discrimination between cloud and surface echoes. Ground based echo analysis permits precise ranging, determining the roughness or slopes of the surface as well as the vertical distributions of vegetation illuminated by the laser. Accurate knowledge of the laser beam's pointing angle is needed to prevent height biases when over sloped surfaces. For surfaces with 2 deg. slopes, knowledge of pointing angle of the beam's centroid to about 8 urad is needed to achieve 10 cm height accuracy. GLAS uses a stellar reference system (SRS) to determine the pointing angle of each laser firing relative to inertial space. The SRS uses a high precision star camera oriented toward local zenith and a gyroscope to determine the inertial orientation of the SRS optical bench. The far field pattern of each laser is measured pulse relative to the star camera with a laser reference system (LRS). Optically measuring each laser far field pattern relative to the orientation of the star camera and gyroscope permits the precise pointing angle of each laser pulse to be determined. GLAS will also determine the vertical distributions of clouds and aerosols by measuring the vertical profile of laser energy backscattered by the atmosphere at both 1064 and 532 nm. The 1064 nm measurements use the Si APD detector and profile the height and vertical structure of thicker clouds. The measurements at 532 nm use new highly sensitive photon counting, detectors, and measure the height distributions of very thin Clouds and aerosol layers. With averaging these can be used to determine the height of the planetary boundary layer. The instrument design and expected performance will be discussed.

Abshire, James B.↗
NASA NTRS1994

Aeroacoustic computation of gust-blade interaction

To better understand and address the challenges faced in computing the acoustics of flow fields, test problems must be considered. In the present study, the sound radiated by the interaction of a flat plate with an oncoming gust containing a two component, mean velocity is computed. The gust has a uniform mean flow in x with Mach number M(infinity) equal to 0.5. The gust's mean velocity in y is of smaller amplitude and is given by: v = 0.1 sin(pi/8(x/M(sub infinity) - t)). This problem has been posed for an upcoming ICASE/LaRC workshop on benchmark problems in computational aeroacoustics. A plate with a length of 30 units in x is used. The plate is assumed to be infinitesimally thin and is centered at the origin. All variables are made dimensionless using the scales specified. Acoustic quantities are obtained by numerically integrating the linearized Euler equations. Integration is performed on the computational domain -100.0 less than or equal to x less than or equal to 100.0, -100.0 less than or equal to y less than or equal to 100.0, using unit length grid spacing in x and in y. An integration scheme is sought which will provide accurate solution to the small quantities of interest at a minimal computational expense. Results indicate that with the given discretization a scheme of minimal fourth order accuracy might be adequate to approximate the waves within the given flow. Thus, a variation of the MacCormack scheme with fourth order accuracy in space and second order accuracy in time was chosen. A scheme with sixth order accuracy in space has also been implemented and results compared with those of the fourth order accurate scheme. To ensure no mass flux, zero normal velocity is assigned at the plate. This condition will induce a discontinuity in the pressure across the plate location. Values for the perturbation pressure p' along the surface of the plate are obtained using a one-sided, third order Taylor expansion, such that p'(sub y) = O. In accordance with the Kutta condition, perturbation pressure at the trailing edge is assigned to zero. In the far field, radiation boundary conditions have been implemented. The effectiveness of the far field conditions are validated by computing in a larger computational domain and comparing the results. Early in time, sound waves begin to radiate from the plate. A Doppler effect is observed. After the initial transients disappear, the strongest waves leave the trailing edge at an approximate 45 degree angle. The intensity pattern of pressure fluctuations shows five lobes (of increasing magnitude with increasing downstream direction) emerging from the plate. Undesirable short wave contaminants are observed in the computed pressure distribution along the plate surface. For a more accurate solution at small scales, a more refined discretization will be required.

Martin, James E.↗
NASA NTRS2013

SIM_EXPLORE: Software for Directed Exploration of Complex Systems

Physics-based numerical simulation codes are widely used in science and engineering to model complex systems that would be infeasible to study otherwise. While such codes may provide the highest- fidelity representation of system behavior, they are often so slow to run that insight into the system is limited. Trying to understand the effects of inputs on outputs by conducting an exhaustive grid-based sweep over the input parameter space is simply too time-consuming. An alternative approach called "directed exploration" has been developed to harvest information from numerical simulators more efficiently. The basic idea is to employ active learning and supervised machine learning to choose cleverly at each step which simulation trials to run next based on the results of previous trials. SIM_EXPLORE is a new computer program that uses directed exploration to explore efficiently complex systems represented by numerical simulations. The software sequentially identifies and runs simulation trials that it believes will be most informative given the results of previous trials. The results of new trials are incorporated into the software's model of the system behavior. The updated model is then used to pick the next round of new trials. This process, implemented as a closed-loop system wrapped around existing simulation code, provides a means to improve the speed and efficiency with which a set of simulations can yield scientifically useful results. The software focuses on the case in which the feedback from the simulation trials is binary-valued, i.e., the learner is only informed of the success or failure of the simulation trial to produce a desired output. The software offers a number of choices for the supervised learning algorithm (the method used to model the system behavior given the results so far) and a number of choices for the active learning strategy (the method used to choose which new simulation trials to run given the current behavior model). The software also makes use of the LEGION distributed computing framework to leverage the power of a set of compute nodes. The approach has been demonstrated on a planetary science application in which numerical simulations are used to study the formation of asteroid families.

Burl, Michael↗
NASA NTRS2000

Bridging the Gap from Networking Technologies to Applications: Workshop Report

The objective of the Next Generation Internet (NGI) Federal program is threefold, encompassing development of networking technologies, high-performance network testbeds, and revolutionary applications. There have been notable advances in emerging network technologies and several nationwide testbeds have been established, but the integration of emerging technologies into applications is lagging. To help bridge this gap between developers of NGI networking technologies and developers of NGI applications, the NASA Research and Education Network (NREN) project hosted a two-day workshop at NASA Ames Research Center in August 1999. This paper presents a summary of the results of this workshop and also describes some of the challenges NREN is facing while incorporating new technologies into HPCC and other NASA applications. The workshop focused on three technologies - Quality of Service (QoS), advanced multicast, and security-and five major NGI application areas - telemedicine, digital earth, digital video, distributed data-intensive applications, and computational infrastructure applications. Network technology experts, application developers, and NGI testbed representatives came together at the workshop to promote cross-fertilization between the groups. Presentations on the first day, including an overview of the three technologies, application case studies and testbed status reports, laid the foundation for discussions on the second day. The objective of these latter discussions, held within smaller breakout groups, was to establish a coherent picture of the current status of the various pieces of each of the three technologies, to create a roadmap outlining future technology development, and to offer technological guidance to application developers. In this paper we first present a brief overview of the NGI applications that were represented at the workshop, focusing on the identification of technological advances that have successfully been incorporated in each application and technological challenges that remain. Next we present the technology roadmaps that were created at the workshop, summarizing the status of various mechanisms that are currently under development and forecasting when various advances are likely to occur within the next one-to-three-year time span. Then we identify issues that were raised at the workshop that might hinder technology development or that might impede integration into NGI applications. We also report some specific guidelines that were offered at the workshop to enable application developers to integrate and effectively use emerging NGI technology building blocks. Finally, we describe NREN activities to incorporate emerging technologies into NASA applications. These activities include support for other NASA High-Performance Computing and Communications Program areas such as IPG (Information Power Grid), support for NASA science enterprises such as Earth science and Mars program prototyping activities, support for satellite/terrestrial networking applications such as the TransAtlantic and TransPacific demonstrations and the Interplanetary Internet, support for NASA telemedicine applications such as the Virtual Collaborative Clinic, and participation in NGI advanced technology testbed initiatives such as the QBone and the NTON/Supernet. For each activity we highlight the primary technological challenge that is associated with it.

Johnson, Marjory J.↗