Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 721 records · Page 40

Application of High-Dimensional Fuzzy K-Means Cluster Analysis to CALIOP/CALIPSO Version 4.1 Cloud-Aerosol Discrimination

This study applies fuzzy k-means (FKM) cluster analyses to a subset of the parameters reported in the CALIPSO lidar level 2 data products in order to classify the layers detected as either clouds or aerosols. The results obtained are used to assess the reliability of the cloud–aerosol discrimination (CAD) scores reported in the version 4.1 release of the CALIPSO data products. FKM is an unsupervised learning algorithm, whereas the CALIPSO operational CAD algorithm (COCA) takes a highly supervised approach. Despite these substantial computational and architectural differences, our statistical analyses show that the FKM classifications agree with the COCA classifications for more than 94 % of the cases in the troposphere. This high degree of similarity is achieved because the lidar-measured signatures of the majority of the clouds and the aerosols are naturally distinct, and hence objective methods can independently and effectively separate the two classes in most cases. Classification differences most often occur in complex scenes (e.g., evaporating water cloud filaments embedded in dense aerosol) or when observing diffuse features that occur only intermittently (e.g., volcanic ash in the tropical tropopause layer). The two methods examined in this study establish overall classification correctness boundaries due to their differing algorithm uncertainties. In addition to comparing the outputs from the two algorithms, analysis of sampling, data training, performance measurements, fuzzy linear discriminants, defuzzification, error propagation, and key parameters in feature type discrimination with the FKM method are further discussed in order to better understand the utility and limits of the application of clustering algorithms to space lidar measurements. In general, we find that both FKM and COCA classification uncertainties are only minimally affected by noise in the CALIPSO measurements, though both algorithms can be challenged by especially complex scenes containing mixtures of discrete layer types. Our analysis results show that attenuated backscatter and color ratio are the driving factors that separate water clouds from aerosols; backscatter intensity, depolarization, and mid-layer altitude are most useful in discriminating between aerosols and ice clouds; and the joint distribution of backscatter intensity and depolarization ratio is critically important for distinguishing ice clouds from water clouds.

Zeng, Shan↗

Survey of some recent stellar-evolution calculations

The stellar-evolution program developed by Sweigart (1973) was used to compute main sequence turnoff, red giant sequences, and horizontal branch sequences for many values of helium abundance and heavy-element abundance. The immediate aim has been to produce a self-consistent set of evolutionary sequences for wide ranges in the composition and mass. The long-term aim involves application to cluster HR diagrams and the integrated properties of galaxies. Turnoff sequences, red giant sequences, and horizontal branch sequences are discussed.

Sweigart, A. V.↗

Flow Analysis of Space Shuttle Feed Line 17-inch Disconnect Valve

A steady incompressible three-dimensional viscous flow analysis has been conducted for the Space Shuttle External Tank/Orbiter propellant feed line disconnect flapper valves with upstream elbows. The full Navier-Stokes code, INS3D, is modified to handle interior obstacles. Grids are generated by SVTGD3D code. Two dimensional initial grids in the flow cross section with and without the flappers are improved by elliptic smoothing to provide better orthogonality, clustering and smoothness to the three dimensional grid. The flow solver is tested for stability and convergence in the presence of interior flappers. An under-relaxation scheme has been incorporated to improve the solution stability. Important flow characteristics such as secondary flows, recirculation, vortex and wake regions, and separated flows are observed. Computed values for forces, moments, and pressure drop are in satisfactory agreement with water flow test data covering a maximum tube Reynolds number of 3.5 x 10(exp 6). The results will serve as a guide to improved design and enhanced testing of the disconnect.

Kandula, Max↗

When Gravity Fails: Local Search Topology

Local search algorithms for combinatorial search problems frequently encounter a sequence of states in which it is impossible to improve the value of the objective function; moves through these regions, called {\em plateau moves), dominate the time spent in local search. We analyze and characterize {\em plateaus) for three different classes of randomly generated Boolean Satisfiability problems. We identify several interesting features of plateaus that impact the performance of local search algorithms. We show that local minima tend to be small but occasionally may be very large. We also show that local minima can be escaped without unsatisfying a large number of clauses, but that systematically searching for an escape route may be computationally expensive if the local minimum is large. We show that plateaus with exits, called benches, tend to be much larger than minima, and that some benches have very few exit states which local search can use to escape. We show that the solutions (i.e. global minima) of randomly generated problem instances form clusters, which behave similarly to local minima. We revisit several enhancements of local search algorithms and explain their performance in light of our results. Finally we discuss strategies for creating the next generation of local search algorithms.

Frank, Jeremy↗

Time Dependence of Collision Probabilities During Satellite Conjunctions

The NASA Conjunction Assessment Risk Analysis (CARA) team has recently implemented updated software to calculate the probability of collision (P (sub c)) for Earth-orbiting satellites. The algorithm can employ complex dynamical models for orbital motion, and account for the effects of non-linear trajectories as well as both position and velocity uncertainties. This “3D P (sub c)” method entails computing a 3-dimensional numerical integral for each estimated probability. Our analysis indicates that the 3D method provides several new insights over the traditional “2D P (sub c)” method, even when approximating the orbital motion using the relatively simple Keplerian two-body dynamical model. First, the formulation provides the means to estimate variations in the time derivative of the collision probability, or the probability rate, R (sub c). For close-proximity satellites, such as those orbiting in formations or clusters, R (sub c) variations can show multiple peaks that repeat or blend with one another, providing insight into the ongoing temporal distribution of risk. For single, isolated conjunctions, R (sub c) analysis provides the means to identify and bound the times of peak collision risk. Additionally, analysis of multiple actual archived conjunctions demonstrates that the commonly used “2D P (sub c)” approximation can occasionally provide inaccurate estimates. These include cases in which the 2D method yields negligibly small probabilities (e.g., P (sub c)) is greater than 10 (sup -10)), but the 3D estimates are sufficiently large to prompt increased monitoring or collision mitigation (e.g., P (sub c) is greater than or equal to 10 (sup -5)). Finally, the archive analysis indicates that a relatively efficient calculation can be used to identify which conjunctions will have negligibly small probabilities. This small-P (sub c) screening test can significantly speed the overall risk analysis computation for large numbers of conjunctions.

Hall, Doyle T.↗

Hubble Space Telescope Astrometry of the Procyon System

The nearby star Procyon is a visual binary containing the F5 IV-V subgiant Procyon A, orbited in a 40.84-year period by the faint DQZ white dwarf (WD) Procyon B. Using images obtained over two decades with the Hubble Space Telescope, and historical measurements back to the 19th century, we have determined precise orbital elements. Combined with measurements of the parallax and the motion of the A component, these elements yield dynamical masses of 1.478 plus or minus 0.012M and 0.592 plus or minus 0.006M for A and B, respectively. The mass of Procyon A agrees well with theoretical predictions based on asteroseismology and its temperature and luminosity. Use of a standard core-overshoot model agrees best for a surprisingly high amount of core overshoot. Under these modeling assumptions, Procyon A's age is approximately 2.7 Gyr. Procyon B's location in the H-R diagram is in excellent agreement with theoretical cooling tracks for WDs of its dynamical mass. Its position in the mass-radius plane is also consistent with theory, assuming a carbon-oxygen core and a helium-dominated atmosphere. Its progenitor's mass was 1.9-2.2M, depending on its amount of core overshoot. Several astrophysical puzzles remain. In the progenitor system, the stars at periastron were separated by only approximately AU, which might have led to tidal interactions and even mass transfer; yet there is no direct evidence that these have occurred. Moreover the orbital eccentricity has remained high (approximately 0.40). The mass of Procyon B is somewhat lower than anticipated from the initial-to-final-mass relation seen in open clusters. The presence of heavy elements in its atmosphere requires ongoing accretion, but the place of origin is uncertain.

PROCYON↗

Verification and Validation of Elastodynamic Simulation Software for Aerospace Research

Physics-based simulation of nondestructive evaluation (NDE) inspection can help to advance the inspectability and reliability of mechanical systems. However, NDE simulations applicable to non-idealized mechanical components often require large compute domains and long run times. This has prompted development of custom NDE simulation software tailored to high performance computing (HPC) hardware. Verification and validation (V&V) is an integral part of developing this software to ensure implementations are robust and applicable to inspection problems, producing tools and simulations suitable for computational NDE research. This presentation addresses factors common to V&V of several elastodynamic simulation codes applicable to ultrasonic NDE. Examples are drawn from in-house simulation software at NASA Langley Research Center, ranging from ensuring reliability in a 1D heterogeneous media wave equation solver to the V&V needs of 3D cluster-parallel elastodynamic software. Factors specific to a research environment are addressed, where individual simulation results can be as relevant as the software product itself. Distinct facets of V&V are discussed including testing to establish software reliability, employing systematic approaches for consistency with fundamental conservation laws, establishing the numerical stability of algorithms, and demonstrating concurrence with empirical data. This talk also addresses V&V practices for small groups of researchers. This includes establishing resources (e.g. time and personnel) for V&V during project planning to mitigate and control the risk of setbacks. Similarly, we identify ways for individual researchers to use V&V during simulation software development itself to both speed up the development process and reduce incurred technical debt.

NDE↗

Lens Model and Time Delay Predictions for the Sextuply Lensed Quasar SDSS J2222+2745*

SDSS J2222+2745 is a galaxy cluster at z = 0.49, strongly lensing a quasar at z = 2.805 into six widely separated images. In recent Hubble Space Telescope imaging of the field, we identify additional multiply lensed galaxies and confirm the sixth quasar image that was identified by Dahle et al. We used the Gemini-North telescope to measure a spectroscopic redshift of z = 4.56 of one of the lensed galaxies. These data are used to refine the lens model of SDSS J2222+2745, compute the time delay and magnifications of the lensed quasar images, and reconstruct the source image of the quasar host and a lensed galaxy at z = 2.3. This galaxy also appears in absorption in our Gemini spectra of the lensed quasar, at a projected distance of 34 kpc. Our model is in agreement with the recent time delay measurements of Dahle et al., who found T(sub AB) = 47.7 +/- 6.0 days and T(sub AC) = 722 +/- 24 days. We use the observed time delays to further constrain the model, and find that the model-predicted time delays of the three faint images of the quasar are T(sub AD) = 502+/- 68 days, T( sub AE) = 611 +/- 75 days, and T(sub AF) = 415 +/- 72 days. We have initiated a follow-up campaign to measure these time delays with Gemini North. Finally, we present initial results from an X-ray monitoring program with Swift, indicating the presence of hard X-ray emission from the lensed quasar, as well as extended X-ray emission from the cluster itself, which is consistent with the lensing mass measurement and the cluster velocity dispersion.

Sharon, Keren↗

Cluster after 20 Years of Operations: Science Highlights and Technical Challenges

The Cluster mission was the first constellation using four identical spacecraft to study Sun-Earth connection plasma processes. Using four spacecraft in a tetrahedron shape, it could measure, for the first time, 3D quantities such as electrical currents, plasma gradients or divergence of the electron pressure tensor and 3D structures such as boundaries, surface waves or vortices. Launched in pairs in July and August 2000, on two Soyuz rockets from Baikonur, the four spacecraft have been collecting data continuously for more than 20 years. The mission faced many challenges during the years of operations as some spacecraft subsystems had a lifetime of a few years beyond the initial two-year mission. The major one was to operate without functioning batteries and to successfully pass short and long eclipses, up to 3 h long, without damaging the on-board computers and transmitters and without freezing the fuel. More than 1,000 eclipses have been successfully passed since 2010 using a specially made procedure which switches off the complete spacecraft before entering into eclipse and switches it on when the Sun is again illuminating the solar panels. During 20 years, many discoveries and science results have been published in more than 2,700 scientific papers. A few highlights are presented here, focusing on how varying the spacecraft separation was essential to achieve the science goals of the mission. The Cluster Science Data System and the Cluster archive allows public access to all science data as well as spacecraft ancillary data.

C P Escoubet↗

Computational Studies of the Interaction of H/H2 with Diamond and Silicon Surfaces

The interaction of hydrogen atoms and molecules with diamond and silicon surfaces is important in several important applications. Two areas that we are interested are: 1) tribology (molecular level friction) and 2) the role of H atoms in silicon chemical vapor deposition (CVD). In the tribology area, H atoms can be used to tie off dangling bonds, which otherwise form bonds between adjacent surfaces, and lead to resistance to sliding the surfaces by each other. Processes which are important in understanding molecular level friction include barriers to addition of H/H2 to the surface and barriers to migration of H atoms on the surface. In the silicon CVD area, we have studied the process of H2 elimination from the 100 surface of silicon. Cluster models for the dime surfaces of diamond are presented. The unrelaxed 100 surface has carbene like surface carbon atoms; however, for the relaxed surface these dimerize to give rows of surface dimers and there is a significant amount of p bonding between the radical orbitals of the dimer. The 110 surface has zig-zag rows of carbon atoms with a dangling bond on each carbon atom. These dangling bonds are hybridized away from each other and thus interact less strongly than for the 100 surface. Finally, the 111 surface has surface C atoms arranged in a triangular pattern and the surface dangling bonds are well separated from each other (second nearest neighbor distance) leading to almost no interaction between adjacent dangling bonds. These qualitative features may be quantified by computing the overlap of adjacent dangling bonds in a GVB(pp) calculation. The overlaps are 0.462, 0.292, and 0.016 for the diamond 100, 110, and 111 surfaces, respectively.

Walch, Stephen P.↗

Facilitating Navigation Through Large Archives

Automated Visual Access (AVA) is a computer program that effectively makes a large collection of information visible in a manner that enables a user to quickly and efficiently locate information resources, with minimal need for conventional keyword searches and perusal of complex hierarchical directory systems. AVA includes three key components: (1) a taxonomy that comprises a collection of words and phrases, clustered according to meaning, that are used to classify information resources; (2) a statistical indexing and scoring engine; and (3) a component that generates a graphical user interface that uses the scoring data to generate a visual map of resources and topics. The top level of an AVA display is a pictorial representation of an information archive. The user enters the depicted archive by either clicking on a depiction of subject area cluster, selecting a topic from a list, or entering a query into a text box. The resulting display enables the user to view candidate information entities at various levels of detail. Resources are grouped spatially by topic with greatest generality at the top layer and increasing detail with depth. The user can zoom in or out of specific sites or into greater or lesser content detail.

Shelton, Robert O.↗

Discovery of Activities via Statistical Clustering of Fixation Patterns

Human behavior often consists of a series of distinct activities, each characterized by a unique pattern of interaction with the visual environment. This is true even in a restricted domain, such as a piloting an aircraft, where activities with distinct visual signatures might be things like communicating, navigating, and monitoring. We propose a novel analysis method for gaze-tracking data, to perform blind discovery of these hypothetical activities. The method is in some respects similar to recurrence analysis, but here we compare not individual fixations, but groups of fixations aggregated over a fixed time interval. The duration of this interval is a parameter that we will refer to as delta. We assume that the environment has been divided into a set of N different areas-of-interest (AOIs). For a given interval of time of duration delta, we compute the proportion of time spent fixating each AOI, resulting in an N-dimensional vector. These proportions can be converted to integer counts by multiplying by delta divided by the average fixation duration (another parameter that we fix at 280 milliseconds). We compare different intervals by computing the chi-square statistic. The p-value associated with the statistic is the likelihood of observing the data under the hypothesis that the data in the two intervals were generated by a single process with a single set of probabilities governing the fixation of each AOI. The method has been applied to approximately 100 hours of eye movement data collected from pilots in a high-fidelity B747 flight simulator, and the results have been compared to synthetic data in which the each activity is represented as first-order Markov process with random probabilities assigned to the AOIs. Randomly-generated synthetic activities can require thousands of fixations to be discriminated with statistical significance, while the human data can be clustered using averaging windows of some 10's of seconds, suggesting that the actual activities are much more narrowly focused than random Markov models.

activity analysis↗

Shuttle onboard IMU alignment methods

The current approach to the shuttle IMU alignment is based solely on the Apollo Deterministic Method. This method is simple, fast, reliable and provides an accurate estimate for the present cluster to mean of 1,950 transformation matrix. If four or more star sightings are available, the application of least squares analysis can be utilized. The least squares method offers the next level of sophistication to the IMU alignment solution. The least squares method studied shows that a more accurate estimate for the misalignment angles is computed, and the IMU drift rates are a free by-product of the analysis. Core storage requirements are considerably more; estimated 20 to 30 times the core required for the Apollo Deterministic Method. The least squares method offers an intermediate solution utilizing as much data that is available without a complete statistical analysis as in Kalman filtering.

Henderson, D. M.↗

Multiple-Bit Errors Caused By Single Ions

Report describes experimental and computer-simulation study of multiple-bit errors caused by impingement of single energetic ions on 256-Kb dynamic random-access memory (DRAM) integrated circuit. Studies illustrate effects of different mechanisms for transport of charge from ion tracks to various elements of integrated circuits. Shows multiple-bit errors occur in two different types of clusters about ion tracks causing them.

Zoutendyk, John A.↗

High order parallel numerical schemes for solving incompressible flows

The use of parallel computers for numerically solving flow fields has gained much importance in recent years. This paper introduces a new high order numerical scheme for computational fluid dynamics (CFD) specifically designed for parallel computational environments. A distributed MIMD system gives the flexibility of treating different elements of the governing equations with totally different numerical schemes in different regions of the flow field. The parallel decomposition of the governing operator to be solved is the primary parallel split. The primary parallel split was studied using a hypercube like architecture having clusters of shared memory processors at each node. The approach is demonstrated using examples of simple steady state incompressible flows. Future studies should investigate the secondary split because, depending on the numerical scheme that each of the processors applies and the nature of the flow in the specific subdomain, it may be possible for a processor to seek better, or higher order, schemes for its particular subcase.

Lin, Avi↗

Sulfur in Cometary Dust

The computer-intensive project consisted of the analysis and synthesis of existing data on composition of comet Halley dust particles. The main objective was to obtain a complete inventory of sulfur containing compounds in the comet Halley dust by building upon the existing classification of organic and inorganic compounds and applying a variety of statistical techniques for cluster and cross-correlational analyses. A student hired for this project wrote and tested the software to perform cluster analysis. The following tasks were carried out: (1) selecting the data from existing database for the proposed project; (2) finding access to a standard library of statistical routines for cluster analysis; (3) reformatting the data as necessary for input into the library routines; (4) performing cluster analysis and constructing hierarchical cluster trees using three methods to define the proximity of clusters; (5) presenting the output results in different formats to facilitate the interpretation of the obtained cluster trees; (6) selecting groups of data points common for all three trees as stable clusters. We have also considered the chemistry of sulfur in inorganic compounds.

Fomenkova, M. N.↗

Hurricane Forecasting with the High-resolution NASA Finite-volume General Circulation Model

A high-resolution finite-volume General Circulation Model (fvGCM), resulting from a development effort of more than ten years, is now being run operationally at the NASA Goddard Space Flight Center and Ames Research Center. The model is based on a finite-volume dynamical core with terrain-following Lagrangian control-volume discretization and performs efficiently on massive parallel architectures. The computational efficiency allows simulations at a resolution of a quarter of a degree, which is double the resolution currently adopted by most global models in operational weather centers. Such fine global resolution brings us closer to overcoming a fundamental barrier in global atmospheric modeling for both weather and climate, because tropical cyclones and even tropical convective clusters can be more realistically represented. In this work, preliminary results of the fvGCM are shown. Fifteen simulations of four Atlantic tropical cyclones in 2002 and 2004 are chosen because of strong and varied difficulties presented to numerical weather forecasting. It is shown that the fvGCM, run at the resolution of a quarter of a degree, can produce very good forecasts of these tropical systems, adequately resolving problems like erratic track, abrupt recurvature, intense extratropical transition, multiple landfall and reintensification, and interaction among vortices.

Atlas, R.↗

Molecular Dynamics Investigation of the Structural and Mechanical Properties of Off-Stoichiometric Epoxy Resins

Molecular dynamic (MD) simulations were performed to compute the mechanical properties of off-stoichiometric epoxy resins as a function of hardener/epoxy mixture ratio (r). Properties were characterized in relation to their microscopic structures. Such resins have been used recently for adhesive-free bonding of large-scale composite structures using the co-curing-ply method. In this process, two partially precured composite panels with hardener-poor (HP) off-stoichiometric resins are coupled with ply(ies) of complementary hardener-rich (HR) formulations and then cured simultaneously. This bonding process has the potential to produce reliable and certifiable composite joints without the need for additional fasteners, which are often required for many conventional bonding methods because even small amounts of contamination can cause a weak bond. The reflow and mixing of the HP/HR resin in this bonding process result in a joint with no discernable interface that should not be susceptible to surface contamination. However, incomplete mixing of the two offset resins may result in chemical heterogeneity of the cured polymeric joint. Thus, different r values may be obtained across the joint. Classical MD simulations were performed to compute the Young’s modulus of polymers with different r values and correlate their properties to network structures. High stiffness was associated with molecular packing due to chemical crosslinking, leading to a single network structure. Moreover, the networks became denser as the ratio approached the stoichiometric value r = 1. Thus, the r = 1 systems were single clusters, with high stiffness, high molecular weight, and a high degree of crosslinking. Structural properties such as radius of gyration and mean square displacement were determined to investigate the variation in the stiffness with respect to r. This MD simulation study was validated with experimental measurements.

Chang Woon Jang↗