Search NASA⌕ Search

SEARCH · Search NASA

Results for “cluster computing”

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 379 records · Page 21

Phase-Retrieval Uncertainty Estimation and Algorithm Comparison for the JWST-ISIM Test Campaign

Phase retrieval, the process of determining the exitpupil wavefront of an optical instrument from image-plane intensity measurements, is the baseline methodology for characterizing the wavefront for the suite of science instruments (SIs) in the Integrated Science Instrument Module (ISIM) for the James Webb Space Telescope (JWST). JWST is a large, infrared space telescope with a 6.5-meter diameter primary mirror. JWST is currently NASA's flagship mission and will be the premier space observatory of the next decade. ISIM contains four optical benches with nine unique instruments, including redundancies. ISIM was characterized at the Goddard Space Flight Center (GSFC) in Greenbelt, MD in a series of cryogenic vacuum tests using a telescope simulator. During these tests, phase-retrieval algorithms were used to characterize the instruments. The objective of this paper is to describe the Monte-Carlo simulations that were used to establish uncertainties (i.e., error bars) for the wavefronts of the various instruments in ISIM. Multiple retrieval algorithms were used in the analysis of ISIM phase-retrieval focus-sweep data, including an iterativetransform algorithm and a nonlinear optimization algorithm. These algorithms emphasize the recovery of numerous optical parameters, including low-order wavefront composition described by Zernike polynomial terms and high-order wavefront described by a point-by-point map, location of instrument best focus, focal ratio, exit-pupil amplitude, the morphology of any extended object, and optical jitter. The secondary objective of this paper is to report on the relative accuracies of these algorithms for the ISIM instrument tests, and a comparison of their computational complexity and their performance on central and graphical processing unit clusters. From a phase-retrieval perspective, the ISIM test campaign includes a variety of source illumination bandwidths, various image-plane sampling criteria above and below the Nyquist- Shannon critical sampling value, various extended object sizes, and several other impactful effects.

Design Analysis↗

Software Testing and Verification in Climate Model Development

Over the past 30 years most climate models have grown from relatively simple representations of a few atmospheric processes to a complex multi-disciplinary system. Computer infrastructure over that period has gone from punch card mainframes to modem parallel clusters. Model implementations have become complex, brittle, and increasingly difficult to extend and maintain. Existing verification processes for model implementations rely almost exclusively upon some combination of detailed analysis of output from full climate simulations and system-level regression tests. In additional to being quite costly in terms of developer time and computing resources, these testing methodologies are limited in terms of the types of defects that can be detected, isolated and diagnosed. Mitigating these weaknesses of coarse-grained testing with finer-grained "unit" tests has been perceived as cumbersome and counter-productive. In the commercial software sector, recent advances in tools and methodology have led to a renaissance for systematic fine-grained testing. We discuss the availability of analogous tools for scientific software and examine benefits that similar testing methodologies could bring to climate modeling software. We describe the unique challenges faced when testing complex numerical algorithms and suggest techniques to minimize and/or eliminate the difficulties.

Clune, Thomas L.↗

H2O/Ni(100) and NH3/Ni(100) - A computational approach

The adsorption of NH3 and H2O on the Ni(100) is treated using a cluster model. The adsorption is found to have only a small effect on the HXH angle. Rotation about the principal ligand axis requires virtually no energy. Ligand tilts of 15 deg are found to require small amounts of energy, leading to the suggestion that the ESDIAD signal is a result of excited tilting modes, not a change in the HXH angle. For H2O additional bending modes are considered and all are found to be quite flat in energy.

Bauschlicher, C. W., Jr.↗

A distributed version of the NASA Engine Performance Program

Distributed NEPP, a version of the NASA Engine Performance Program, uses the original NEPP code but executes it in a distributed computer environment. Multiple workstations connected by a network increase the program's speed and, more importantly, the complexity of the cases it can handle in a reasonable time. Distributed NEPP uses the public domain software package, called Parallel Virtual Machine, allowing it to execute on clusters of machines containing many different architectures. It includes the capability to link with other computers, allowing them to process NEPP jobs in parallel. This paper discusses the design issues and granularity considerations that entered into programming Distributed NEPP and presents the results of timing runs.

Cours, Jeffrey T.↗

Adiabatic invariants in stellar dynamics, 3: Application to globular cluster evolution

The previous two companion papers demonstrate that slowly varying perturbations may not result in adiabatic cutoffs and provide a formalism for computing the long-term effects of time-dependent perturbations on stellar systems. Here, the theory is implemented in a Fokker-Planck code and a suite of runs illustrating the effects of shock heating on globular cluster evolution are described. Shock heating alone results in considerable mass loss for clusters with R(sub g) less than or approximately 8 kpc: a concentration c = 1.5 cluster with R(sub g) kpc loses up to 95% of its initial mass in 15 Gyr. Only those with concentration c greater than or approximately 1.3 survive disk shocks inside of this radius. Other effects, such as mass loss by stellar evolution, will decrease this survival bound. Loss of the initial halo together with mass segregation leads to mass spectral indices, x, which may be considerably larger than their initial values.

Weinberg, Martin D.↗

Deciphering the Spectra of Flowers to Map Landscape-scale Blooming Dynamics

Like leaves, floral coloration is driven by inherent optical properties, which are determined by pigments, scattering structure, and thickness. However, establishing the relative contribution of these factors to canopy spectral signals is usually limited to in-situ observations. Modeling flowering dynamics (e.g., blooming duration, spatial distribution) at the landscape scale may reveal insights into ecological processes and phenological adaptations to environmental changes. Multitemporal visible to shortwave infrared (VSWIR) imaging spectroscopy observations are especially suited for such efforts. Reflectance in this spectral range is sensitive to major flower pigments, flowering phenology traces, and biophysical differences between flowers and other plant parts. We explored how flowers contribute to spectral signals using a time series of imagery from the Airborne Visible InfraRed Imaging Spectrometer - Next Generation (AVIRIS-NG) collected as part of the SBG High-Frequency Time Series (SHIFT) campaign as a case study. Airborne data were collected weekly during the spring of 2022 across two natural reserves in California. Field spectra were gathered from blooming plots at leaf, flower, and canopy levels at two time points during the campaign. The processed data was used to investigate flowering species' spectro-temporal variation and spatial distribution using Spectral Mixture Residual, Gaussian clustering techniques, and a proposed narrow-band flowering index. Linear spectral unmixing allowed the computation of the weighted contribution of four major high-variance endmembers (leaves, flowers, soil, dark) and low-variance residual signal that comprises subtle spectral features used to track biophysical processes. The reflectance residual was projected on a low principal component basis to characterize flowering clusters' variation and spatial distribution based on the Gaussian mixture model, providing an uncertainty metric to assess the results. Mapping flowering events from modeling spectro-temporal dynamics throughout the season, from pre-blooming to post-flowering stages, allowed us to identify gradient variations in spectral features within the VSWIR spectral range linked to flowering pigments. Time series of the Mixture Residual Blooming Index and the Red-Edge Normalized Difference Vegetation Index revealed specific flowering and greenness phenophases across the two main species (Coreopsis gigantea, Artemisia californica) in the flowering areas. Overall, our approach opens opportunities for future satellite monitoring of floral cycles at broader scales.

Yoseline Angel↗

Evaluating computer-drawn ground-cover maps

Computer-generated character maps from LANDSAT data are compared to aerial photos for test sites in Florida. Report Describes extraction of ground features by two analytical techniques: unsupervised clustering algorithm, called LANDSAT Signature Development Program (LSDP), and interactive algorithm based on multispectral image analyzer. Study concluded that computer classification of digital LANDSAT multispectral data, supplemented with certain ground-cover information, is valuable tool for analysis of renewable resources.

Arvanitis, L. G.↗

Numerical solution of Space Shuttle Orbiter flow field

The supersonic, viscous laminar flow around the Space Shuttle Orbiter forebody has been computed with a parabolized Navier-Stokes code using a generalized coordinate transformation. The initial solution for the nose part of the Orbiter geometry was obtained with a three-dimensional time-dependent Navier-Stokes solver. It was necessary to employ a wind axis oriented coordinate system to obtain the initial solution with the time-dependent code. The generalized PNS technique used in this study allows the solution to be marched from the given initial data surface to any desired surface downstream. A grid point clustering scheme was employed to accurately describe the body shape by clustering points at the wing tip and at the wing body juncture. The computed heat transfer coefficients, pressure coefficients, and shock shapes are compared with the available experimental data for 0 degrees and 30 degrees angle of attack.

Venkatapathy, E.↗

Buckets, Clusters and Dienst

In this paper we describe NCSTRL+, a unified, canonical digital library for scientific and technical information (STI). NCSTRL+ is based on the Networked Computer Science Technical Report Library (NCSTRL), a World Wide Web (WWW) accessible digital library (DL) that provides access to over 80 university departments and laboratories. NCSTRL+ implements two new technologies: cluster functionality and publishing "buckets." We have extended the Dienst protocol, the protocol underlying NCSTRL, to provide the ability to "cluster" independent collections into a logically centralized digital library based upon subject category classification, type of organization, and genres of material. The concept of "buckets" provides a mechanism for publishing and managing logically linked entities with multiple data formats. The NCSTRL+ prototype DL contains the holdings of NCSTRL and the NASA Technical Report Server (NTRS). The prototype demonstrates the feasibility of publishing into a multi-cluster DL, searching across clusters, and storing and presenting buckets of information. We show that the overhead for these additional capabilities is minimal to both the author and the user when compared to the equivalent process within NCSTRL.

Nelson, Michael L.↗

Multiple-fold clustered processor mesh array

The multiple-fold clustered processor mesh array is a triangular organization of clustered processing elements. This multiple-fold array maintains functional equivalence to the nearest neighbor mesh computer with uni-directional interprocessor communications, but with half the number of connection wires. In addition, the connectivity of the multiple-folded organization is superior to the standard square mesh due to the improved connectivity between the clustered processors. One of the primary application areas targeted is High Performance Architectures for image processing.

Pechanek, Gerald G.↗

DHARMA - Discriminant hyperplane abstracting residuals minimization algorithm for separating clusters with fuzzy boundaries

Learning of discriminant hyperplanes in imperfectly supervised or unsupervised training sample sets with unreliably labeled samples along the fuzzy joint boundaries between sample clusters is discussed, with the discriminant hyperplane designed to be a least-squares fit to the unreliably labeled data points. (Samples along the fuzzy boundary jump back and forth from one cluster to the other in recursive cluster stabilization and are considered unreliably labeled.) Minimization of the distances of these unreliably labeled samples from the hyperplanes does not sacrifice the ability to discriminate between classes represented by reliably labeled subsets of samples. An equivalent unconstrained linear inequality problem is formulated and algorithms for its solution are indicated. Landsat earth sensing data were used in confirming the validity and computational feasibility of the approach, which should be useful in deriving discriminant hyperplanes separating clusters with fuzzy boundaries, given supervised training sample sets with unreliably labeled boundary samples.

Dasarathy, B. V.↗

MEDPRAT Treatment Clusters: Improving Representation of Mission Medical Risk

INTRODUCTION The Medical Extensible Dynamic Probabilistic Risk Assessment Tool (MEDPRAT) implements a computational model that aims to quantify spaceflight medical risk by utilizing probabilistic techniques to simulate critical event incidence and outcomes over thousands of simulated mission trials. The goal of MEDPRAT is to characterize mission medical risk and provide insight into medical resource utilization. In order to analyze the medical resource space, treatment must be mapped from each simulated condition, and resources consumed as a result of this treatment must be tracked throughout the course of the mission. A new MEDPRAT feature, ‘treatment clusters’, provide a more sophisticated method of defining the structure and interaction between resources, more closely mimicking the way treatment is carried out clinically. METHODS Treatment clusters expand on the two existing treatment groupings (combination and alternate) adding a new grouping: bundled treatment. Treatment clusters may be combined to any depth, giving users the ability to specify complex treatment trees whose behavior is governed by several user-specified parameters. This approach emphasizes reusability, as treatment clusters, once defined, can be used to create more complex treatment trees or applied to many conditions. By configuring parameters for contribution, efficacy, necessity, primacy, and equivalence, resource relationships and dependencies can be more accurately represented, thereby allowing users to build capabilities with desired treatment properties, for example an intravenous capability for conditions such as anaphylaxis, acute radiation syndrome, etc. MEDPRAT v1.0 remains backward compatible with existing treatment structures, giving users the ability to define new treatment clusters as evidence becomes available, without having to recode their existing treatment databases. In addition to facilitating the representation of more complex treatment options, by pairing treatment clusters with the internal optimization routine, the MEDPRAT set selector, medical resources can be identified as organized in bundles, where appropriate, so that optimized resource sets include groups of highly-dependent resources only when all resources of the group are together. For example, it would be wasteful to include ultrasound gel but not an ultrasound machine, since the gel provides no benefit as a treatment without the ultrasound machine. With treatment clusters, the user may require that both resources are available to provide any benefit as treatment, so that if one resource is optimized out of the set, the other resource will be optimized out as well. RESULTS AND CONCLUSIONS We will report on MEDPRAT treatment clusters used in a bundling study under the IMPACT project of the ExMC element. We will discuss an example of a complex treatment tree. Through the implementation of this feature MEDPRAT enables treatment to be defined and applied in a way that is more representative of the real world, providing more accurate insight into mission medical risk and the medical resource space.

Lawrence Leinweber↗

Implicit finite-difference procedures for the computation of vortex wakes

Implicit finite-difference procedures for the primitive form of the incompressible Navier-Stokes and the compressible Euler equations are used to compute vortex wake flows. The partial differential equations in strong conservation-law form are transformed to cluster grid points in regions with large changes in vorticity. In addition to clustering, fourth-order accurate, spatial difference operators are used to help resolve the flow-field gradients. The use of implicit time-differencing permits large time steps to be taken since temporal variations are typically small. Computational efficiency is achieved by approximate factorization. Both two-dimensional and preliminary three-dimensional calculations are described and qualitatively compared with existing experimental data.

Steger, J. L.↗

Small Al clusters. I - The effect of basis set and correlation on the geometry of small Al clusters

A detailed study is presented of the basis set requirements and effects of correlation on the geometry and structure of small Al(n) clusters n = 2, 4, and 13. An effective core potential (ECP) is developed from the Huzinaga basis which accurately reproduces the all-electron results. It is found that depolarization functions are very important in computing the bond length, and that the bond contraction obtained is about as large for a 13-atom cluster as for the Al(4) rhombus structure. With d functions on each center in Al(13), a bond distance shorter than the bulk is obtained, as expected. The inclusion of correlation in Al(4) is found to be less important than the addition of polarization functions for determining the bond length, but increases substantially the calculated atomization energy. These calibration calculations allow an accurate treatment for Al(n) clusters to be designed.

Bauschlicher, Charles W., Jr.↗

Parallel and Distributed Computational Fluid Dynamics: Experimental Results and Challenges

This paper describes several results of parallel and distributed computing using a large scale production flow solver program. A coarse grained parallelization based on clustering of discretization grids combined with partitioning of large grids for load balancing is presented. An assessment is given of its performance on distributed and distributed-shared memory platforms using large scale scientific problems. An experiment with this solver, adapted to a Wide Area Network execution environment is presented. We also give a comparative performance assessment of computation and communication times on both the tightly and loosely-coupled machines.

Djomehri, Mohammad Jahed↗

Scalable Multiprocessor for High-Speed Computing in Space

A report discusses the continuing development of a scalable multiprocessor computing system for hard real-time applications aboard a spacecraft. "Hard realtime applications" signifies applications, like real-time radar signal processing, in which the data to be processed are generated at "hundreds" of pulses per second, each pulse "requiring" millions of arithmetic operations. In these applications, the digital processors must be tightly integrated with analog instrumentation (e.g., radar equipment), and data input/output must be synchronized with analog instrumentation, controlled to within fractions of a microsecond. The scalable multiprocessor is a cluster of identical commercial-off-the-shelf generic DSP (digital-signal-processing) computers plus generic interface circuits, including analog-to-digital converters, all controlled by software. The processors are computers interconnected by high-speed serial links. Performance can be increased by adding hardware modules and correspondingly modifying the software. Work is distributed among the processors in a parallel or pipeline fashion by means of a flexible master/slave control and timing scheme. Each processor operates under its own local clock; synchronization is achieved by broadcasting master time signals to all the processors, which compute offsets between the master clock and their local clocks.

Lux, James↗

A bound on the rate-distortion function and application to images.

An upper bound on the rate-distortion function for discrete ergodic sources with memory is developed by partitioning the source sample space into a finite number of disjoint subsets and bounding the rates for each subset. The bound depends only on the mean vectors and covariance matrices for the subsets and is easy to compute. It is tighter than the Gaussian bound for sources that exhibit clustering of either the values or covariances of successive source outputs. The bound is evaluated for a certain class of pictorial data using both one-dimensional and two-dimensional blocks of picture elements. Two-dimensional blocks yield a tighter bound than one-dimensional blocks; both result in a significantly tighter bound than the Gaussian bound.

Tasto, M.↗

Utilization of ERTS-1 data in the Houston area

The author has identified the following significant results. Using clustering techniques, several large lakes in Texas have been accurately delineated in computer printout graymaps. It was also found that small bodies of water (one to two acres in size) could be detected by searching for small reflectance values in the infrared data. A graymap printout of a lake described a shore outline that was not consistent with available maps. Field examination revealed that the actual level of the lake was below that for which the map was drawn. The current lake configuration agrees in shape and relative size with the ERTS-1 data printout. Water turbidity causes reflectance changes which are detectable in ERTS-1 band 7 data. A comparison has been made of the Monterey Bay, California area using 1971 aerial color infrared photography and a 1972 ERTS-1 band 7 infrared image. This comparison revealed that some event has occurred to impound a significant amount of water in the area since the infrared photography was taken. Data values in the ERTS-1 infrared image exhibit detectable changes in brightness at inflow points, where high turbidity would be present. Researchers had not expected to detect water turbidity patterns in band 7 (800 to 1100 nanometers).

Erb, R. B.↗