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 883 records · Page 49

Evolutionary models of halo stars with rotation. I - Evidence for differential rotation with depth in stars

Evolutionary models of metal-poor stars are computed including the effects of rotation, and their properties are compared with observations. The models rotate slowly at the surface, in agreement with the observed upper limits on the rotation velocity at main-sequence turnoff; they also have substantial differential rotation with depth. This differential rotation preserves a sufficient amount of internal angular momentum to explain the rapid rotation of evolved horizontal-branch stars. These results hold for a wide range of angular momentum loss and transport parameter values. Differences and similarities between the surface and internal rotation of solar metallicity and metal-poor models are discussed. Rigidly rotating models are found to be incompatible with the observations once giant branch mass loss is taken into account. Horizontal-branch rotation velocity measurements as a function of color are proposed as a test of the rotation law enforced in convection zones, and their dependence on cluster age and metallicity are discussed.

Pinsonneault, Marc H.↗

Is there a classical Hyades lithium problem?

The distribution of surface lithium abundances in the Hyades cluster, as a function of effective temperature, has resisted a satisfactory theoretical explanation for more than two decades. Trends in the historically recomputed values of stellar opacity have stimulated a systematic investigation of the effects of larger opacities on lithium depletion during the premain-sequence evolution of low-mass stars. It is shown that the calculated depletion depends quite strongly on the opacities employed, and that both surface and interior opacities play an important role. By performing standard evolutionary calculations with parameterized increments in the opacities, the observed Hyades Li and T(eff) relation for G and K dwarfs have been determined in a manner consistent with the other observational constraints. The solution presented herein requires an increase in interior opacities at about 4 million K over older values by approximately less than 37 percent. Such a change is not incompatible with the historical trends in opacity computations and suggests that the resolution of this problem lies solely in the opacities, without recourse to other physical mechanisms such as turbulent diffusion or 'extra mixing'.

Swenson, Fritz J.↗

Analysis and visualization of single-trial event-related potentials

In this study, a linear decomposition technique, independent component analysis (ICA), is applied to single-trial multichannel EEG data from event-related potential (ERP) experiments. Spatial filters derived by ICA blindly separate the input data into a sum of temporally independent and spatially fixed components arising from distinct or overlapping brain or extra-brain sources. Both the data and their decomposition are displayed using a new visualization tool, the "ERP image," that can clearly characterize single-trial variations in the amplitudes and latencies of evoked responses, particularly when sorted by a relevant behavioral or physiological variable. These tools were used to analyze data from a visual selective attention experiment on 28 control subjects plus 22 neurological patients whose EEG records were heavily contaminated with blink and other eye-movement artifacts. Results show that ICA can separate artifactual, stimulus-locked, response-locked, and non-event-related background EEG activities into separate components, a taxonomy not obtained from conventional signal averaging approaches. This method allows: (1) removal of pervasive artifacts of all types from single-trial EEG records, (2) identification and segregation of stimulus- and response-locked EEG components, (3) examination of differences in single-trial responses, and (4) separation of temporally distinct but spatially overlapping EEG oscillatory activities with distinct relationships to task events. The proposed methods also allow the interaction between ERPs and the ongoing EEG to be investigated directly. We studied the between-subject component stability of ICA decomposition of single-trial EEG epochs by clustering components with similar scalp maps and activation power spectra. Components accounting for blinks, eye movements, temporal muscle activity, event-related potentials, and event-modulated alpha activities were largely replicated across subjects. Applying ICA and ERP image visualization to the analysis of sets of single trials from event-related EEG (or MEG) experiments can increase the information available from ERP (or ERF) data. Copyright 2001 Wiley-Liss, Inc.

Non-NASA Center↗

An Innovative Approach to Modeling VIPER Rover Software Life Cycle Cost

NASA’s “Volatiles Investigating Polar Exploration Rover” (VIPER) will be the first robotic mission to prospect for water ice near the south pole of the Moon in late 2023 on a 100-Earth-day mission. The information that the VIPER rover provides will help improve understanding of the composition, distribution, and accessibility of Lunar polar volatiles and will help determine how the Moon’s resources can support future human space exploration. VIPER, however, represents a radical departure from the way that NASA has traditionally developed planetary robotic missions. A key consequence of these differences is that estimating the cost of VIPER’s rover software is challenging and complex.For example, VIPER is being developed using management procedures typically applied to NASA research and technology projects, rather than space flight programs. In addition, key portions of the rover’s software are being designed as ground software to run on mission control computers (rather than on-board the rover as flight software as with prior planetary missions) taking advantage of continuous, interactive data communications between the Moon and Earth and higher performance computing available on the ground. Moreover, the rover’s software is being engineered using Agile software development practices and incorporates a significant amount of open-source, rather than following traditional (spiral, waterfall, etc.) development methods and in-house code. In this paper, we present an innovative process to estimate the life cycle cost of VIPER’s rover software. We first describe how we modeled the architecture and code counts for three software elements: Rover Flight Software (RFSW), Rover Ground Software (RGSW), and Rover Simulation Software (RSIM). We then discuss key challenges and unique aspects of our approach, such as the lack of Lunar rover analogies, the need to integrate and test large open source software, and the strategies developed to account for use of non-space flight management practices and the impact of the COVID-19 pandemic. We conclude with a summary of our results, including cumulative distribution, nearest neighbors and cluster analysis, as well as heuristics used to confirm the reasonableness of the cost estimate.

Utz, Hans↗

Lunar Gateway Charging and Effects on Low-Energy Electron Measurements

The Electron Electrostatic Analyzer (EEA) is part of the HERMES instrument suite on Lunar Gateway and is designed to measure electrons in the energy range from 1 eV to 18 keV. Electrons in this range can be deflected or blocked from entering a detector by electric fields from a charged spacecraft platform, so many heliophysics missions (MMS, Cluster, etc.) use active potential control to minimize stray electric fields which would corrupt low-energy plasma measurements. Gateway, however, does not have such capabilities so the surface charging effects will influence the electron data. Further, some elements of Gateway partially block the instrument’s field of view to space. The effects of the Gateway platform on the low-energy electron population must be understood to maximize the accuracy of the EEA data. To investigate the electrostatic environment near Gateway and its effects on the EEA measurements, we use Nascap-2K to model the expected charging of Gateway under representative solar wind, magnetosheathic, and magnetospheric conditions. Next, we use the potential and electric field outputs from Nascap-2K along with in-house particle tracing codes to investigate how low-energy electrons are perturbed, blocked, or deflected by the spacecraft. We discuss how these perturbations impact the plasma moments computed from the EEA data and possible approaches for correcting the charging effects including application of advanced methods such as neural networks. In addition to electrons from the space environment being influenced by the electrostatic environment around Gateway, photoelectrons generated on Gateway surfaces will also be measured by EEA. We are investigating how knowledge of the potentials of various Gateway surfaces can be extracted from these measured photoelectrons using advanced numerical modeling and analysis. Such information would be valuable for validating charging models and planning mission operations around the platform. Finally, we discuss how the specific workflow developed for Gateway can be extended to a generalized process that is applicable to other future missions. As the space industry prepares for human exploration beyond low-Earth orbit, the need to obtain accurate space weather data from rideshare science payloads is critical. Our process for removing platform effects from the science data serves as a pathfinder for leveraging advanced simulations to maximize science returns and ensure the safety and success of future crewed missions.

Miles Bengtson↗

A Coupled-Cluster Study of the Molecular Structure, Vibrational Spectrum, and Heats of Formation of XONO2 (X=H, F, Cl)

The equilibrium structures, harmonic vibrational frequencies, dipole moments, and IR intensities of nitric acid, fluorine nitrate, and chlorine nitrate have been investigated by using the singles and doubles coupled cluster method that also includes a perturbational estimate of the effects of connected triple excitations, CCSD(T). A standard triple-zeta double-polarized basis set was utilized. The equilibrium geometries and vibrational spectra of HONO2 and ClONO2 are shown to be in excellent agreement with the available experimental data. The ab initio vibrational spectrum of FONO2 is also shown to be in excellent agreement with experiment. Unlike the FOOF and FON molecules, but similar to the cis- and trans-FONO molecules, FONO2 is shown to possess normal bond distances. The bonding in FONO2 is shown to be more similar to that in ClONO2 than that in HONO2, although there are still significant differences, especially in the partial atomic charges as deduced from Mulliken populations. This causes FONO2 to possess almost no dipole moment, which is very different to both HONO2 and ClONO2. By using large atomic natural orbital basis sets, CCSD(T) energies are computed for four isodesmic reactions in order to determine an accurate heat of formation for FONO2. Our best estimate for DELTA H(raised circle) (sub f,298) is 3.1 plus or minus 2.0 kcal/mol, indicating that the F-ONO2 bond energy is 31.3 kcal/mol.

Lee, Timothy J.↗

Bayesian Estimation of Precipitation from Satellite Passive Microwave Observations Using Combined Radar-Radiometer Retrievals

Precipitation estimation from satellite passive microwave radiometer observations is a problem that does not have a unique solution that is insensitive to errors in the input data. Traditionally, to make this problem well posed, a priori information derived from physical models or independent, high-quality observations is incorporated into the solution. In the present study, a database of precipitation profiles and associated brightness temperatures is constructed to serve as a priori information in a passive microwave radiometer algorithm. The precipitation profiles are derived from a Tropical Rainfall Measuring Mission (TRMM) combined radar radiometer algorithm, and the brightness temperatures are TRMM Microwave Imager (TMI) observed. Because the observed brightness temperatures are consistent with those derived from a radiative transfer model embedded in the combined algorithm, the precipitation brightness temperature database is considered to be physically consistent. The database examined here is derived from the analysis of a month-long record of TRMM data that yields more than a million profiles of precipitation and associated brightness temperatures. These profiles are clustered into a tractable number of classes based on the local sea surface temperature, a radiometer-based estimate of the echo-top height (the height beyond which the reflectivity drops below 17 dBZ), and brightness temperature principal components. For each class, the mean precipitation profile, brightness temperature principal components, and probability of occurrence are determined. The precipitation brightness temperature database supports a radiometer-only algorithm that incorporates a Bayesian estimation methodology. In the Bayesian framework, precipitation estimates are weighted averages of the mean precipitation values corresponding to the classes in the database, with the weights being determined according to the similarity between the observed brightness temperature principal components and the brightness temperature principal components of the classes. Because the classes are stratified by the sea surface temperature and the echo-top-height estimator, the number of classes that are considered for retrieval is significantly smaller than the total number of classes, making the algorithm computationally efficient. The radiometer-only algorithm is applied to TMI observations, and precipitation estimates are compared with combined TRMM precipitation radar (PR) TMI reference estimates. The TMI-only algorithm, supported by the empirically derived database, produces estimates that are more consistent with the reference values than the precipitation estimates from the version-6 TRMM facility TMI algorithm. Cloud-resolving model simulations are used to assign a latent heating profile to each precipitation profile in the empirically derived database, making it possible to estimate latent heating using the radiometer-only algorithm. Although the evaluation of latent heating estimates in this study is preliminary, because realistic conditional probability distribution functions are attached to latent heating structures in the algorithm s database, a generally positive impact on latent heating estimation from passive microwave observations is expected.

Grecu, Mircea↗

Self-consistent Quantum Iteratively Sparsified Hamiltonian Algorithm (SQuISH)

Due to coherence time limitations, reducing the resources required to run quantum algorithms and simulate physical systems on a quantum computer is crucial. With regards to Hamiltonian simulation, a significant effort has focused on building efficient algorithms using various factorizations and truncations, typically derived from the Hamiltonian alone. We introduce a new paradigm for improving Hamiltonian simulation and reducing the cost of ground state problems based on ideas recently developed for classical chemistry simulations. The key idea is that one can find efficient ways to reduce resources needed by quantum algorithms by making use of two key pieces of information: the Hamiltonian operator and an approximate ground state wavefunction. We refer to our algorithm as the self-consistent quantum iteratively sparsified Hamiltonian (SQuISH). By performing our scheme iteratively, one can drive SQuISH to create an accurate wavefunction using a truncated, resource-efficient Hamiltonian. By utilizing this more compact Hamiltonian, our algorithm provides an approach to reduce the gate complexity of ground state calculations on quantum hardware. As proof of principle, we implement SQuISH using configuration interaction for small molecules and coupled cluster for larger systems. Through our combination of approaches, we demonstrate how it performs on a range of systems, the largest of which would require more than 200 qubits to run on quantum hardware.

Diana Chamaki↗

Adaptive pattern recognition by mini-max neural networks as a part of an intelligent processor

In this decade and progressing into 21st Century, NASA will have missions including Space Station and the Earth related Planet Sciences. To support these missions, a high degree of sophistication in machine automation and an increasing amount of data processing throughput rate are necessary. Meeting these challenges requires intelligent machines, designed to support the necessary automations in a remote space and hazardous environment. There are two approaches to designing these intelligent machines. One of these is the knowledge-based expert system approach, namely AI. The other is a non-rule approach based on parallel and distributed computing for adaptive fault-tolerances, namely Neural or Natural Intelligence (NI). The union of AI and NI is the solution to the problem stated above. The NI segment of this unit extracts features automatically by applying Cauchy simulated annealing to a mini-max cost energy function. The feature discovered by NI can then be passed to the AI system for future processing, and vice versa. This passing increases reliability, for AI can follow the NI formulated algorithm exactly, and can provide the context knowledge base as the constraints of neurocomputing. The mini-max cost function that solves the unknown feature can furthermore give us a top-down architectural design of neural networks by means of Taylor series expansion of the cost function. A typical mini-max cost function consists of the sample variance of each class in the numerator, and separation of the center of each class in the denominator. Thus, when the total cost energy is minimized, the conflicting goals of intraclass clustering and interclass segregation are achieved simultaneously.

Szu, Harold H.↗

Astrometry With the Hubble Space Telescope: Trigonometric Parallaxes of Planetary Nebula Nuclei NGC 6853, NGC 7293, ABELL 31, and DeHt 5

We present absolute parallaxes and relative proper motions for the central stars of the planetary nebulae NGC 6853 (The Dumbbell), NGC 7293 (The Helix), Abell 31, and DeHt 5. This paper details our reduction and analysis using DeHt 5 as an example. We obtain these planetary nebula nuclei (PNNi) parallaxes with astrometric data from Fine Guidance Sensors FGS 1r and FGS 3, white-light interferometers on the Hubble Space Telescope. Proper motions, spectral classifications and VJHKT2M and DDO51 photometry of the stars comprising the astrometric reference frames provide spectrophotometric estimates of reference star absolute parallaxes. Introducing these into our model as observations with error, we determine absolute parallaxes for each PNN. Weighted averaging with previous independent parallax measurements yields an average parallax precision, sigma (sub pi)/ pi = 5%. Derived distances are: d(sub NGC6853) = 405(exp +28 sub -25) pc, d(sub NGC7293) = 216(exp +14 sub -12) pc, d(sub Abell31) = 621(exp +91 sub -70) pc, and d(sub DeHt5) = 345(exp +19 sub -17) pc. These PNNi distances are all smaller than previously derived from spectroscopic analyses of the central stars. To obtain absolute magnitudes from these distances requires estimates of interstellar extinction. We average extinction measurements culled from the literature, from reddening based on PNNi intrinsic colors derived from model SEDs, and an assumption that each PNN experiences the same rate of extinction as a function of distance as do the reference stars nearest (in angular separation) to each central star. We also apply Lutz-Kelker bias corrections. The absolute magnitudes and effective temperatures permit estimates of PNNi radii through both the Stefan-Boltzmann relation and Eddington fluxes. Comparing absolute magnitudes with post-AGB models provides mass estimates. Masses cluster around 0.57 solar Mass, close to the peak of the white dwarf mass distribution. Adding a few more PNNi with well-determined distances and masses, we compare all the PNNi with cooler white dwarfs of similar mass, and confirm, as expected, that PNNi have larger radii than white dwarfs that have reached their final cooling tracks.

STARS↗

Commercial-Off-The-Shelf Microelectromechanical Systems (MEMS) Flow-Measurement Probes Fabricated And Assembled

As an alternative to conventional tubing instrumentation for measuring airflow, designers and technicians at the NASA Glenn Research Center have been fabricating packaging components and assembling a set of unique probes using commercial-off-the-shelf microelectromechanical systems (MEMS) integrated circuits (computer chips). Using MEMS as an alternative has some compelling advantages over standard measurement devices. Sensor technologies have matured through high-production usage in industries such as automotive and aircraft manufacturers. Currently, MEMS are the choice in applications such as tire pressure monitors, altimeters, pneumatic controls, cable leak detectors, and consumer appliances. Conventional instrumentation uses tubing buried in the model aerodynamic surfaces or wind tunnel walls. The measurements are made when pressure is introduced at the tube opening. The pressure then must travel the tubing for lengths ranging from 20 to hundreds of feet before reaching an electronic signal conditioner. This condition causes a considerable amount of damping and requires measurements to be made only after the test rig has reached steady-state operation. The electronic MEMS pressure sensor is able to take readings continuously under dynamic states in nearly real time. The use of stainless steel tubing for pressure measurements requires many tubes to be cleaned, cut to length, carefully installed, and delicately deburred and spliced for use. A cluster of a few hundred 1/16-in.- (0.0625-in.-) diameter tubes (not uncommon in research testing facilities) can be several inches in diameter and may weigh enough to require two men to handle. Replacing hard tubing with electronic chips can eliminate much of the bulk. Each sensor would fit on the tip of the 1/16-in. tubing with room to spare. The P592 piezoresistive silicon pressure sensor (Lucas NovaSensor, Fremont, CA) was chosen for this project because of its cost, availability, and tolerance to extreme ambient conditions. The chip is 1 mm square by 0.6 mm thick (0.039 by 0.023 in.) with 0.12-mm (0.005-in.) wire connection tabs. Three MEMS chips were used to build the first type of flow-angularity probe. This MEMS probe will be demonstrated as an alternative to a standard tube type "Cobra Probe" now used routinely in wind tunnel and aeronautical hardware applications. The response time and accuracy would allow the probe to be translated on an actuator across a flow field, yielding precision dynamic measurements not possible with conventional instrumentation. The low profile, the minimal power requirement, the rugged construction, and the moderate cost all contribute to making MEMS sensors the enticing choice instrument in future research measurement needs. The MEMS probe efforts are a continuation of work initiated by Brian Willis, without whose foresight and efforts this project would never have been realized. This task was funded through cooperation with the NASA Electronic Parts and Packaging (NEPP) program at the Jet Propulsion Laboratory.

Redding, Chip↗

Using Docker Containers to Extend Reproducibility Architecture for the NASA Earth Exchange (NEX)

NASA Earth Exchange (NEX) is a data, supercomputing and knowledge collaboratory that houses NASA satellite, climate and ancillary data where a focused community can come together to address large-scale challenges in Earth sciences. As NEX has been growing into a petabyte-size platform for analysis, experiments and data production, it has been increasingly important to enable users to easily retrace their steps, identify what datasets were produced by which process chains, and give them ability to readily reproduce their results. This can be a tedious and difficult task even for a small project, but is almost impossible on large processing pipelines. We have developed an initial reproducibility and knowledge capture solution for the NEX, however, if users want to move the code to another system, whether it is their home institution cluster, laptop or the cloud, they have to find, build and install all the required dependencies that would run their code. This can be a very tedious and tricky process and is a big impediment to moving code to data and reproducibility outside the original system. The NEX team has tried to assist users who wanted to move their code into OpenNEX on Amazon cloud by creating custom virtual machines with all the software and dependencies installed, but this, while solving some of the issues, creates a new bottleneck that requires the NEX team to be involved with any new request, updates to virtual machines and general maintenance support. In this presentation, we will describe a solution that integrates NEX and Docker to bridge the gap in code-to-data migration. The core of the solution is saemi-automatic conversion of science codes, tools and services that are already tracked and described in the NEX provenance system, to Docker - an open-source Linux container software. Docker is available on most computer platforms, easy to install and capable of seamlessly creating and/or executing any application packaged in the appropriate format. We believe this is an important step towards seamless process deployment in heterogeneous environments that will enhance community access to NASA data and tools in a scalable way, promote software reuse, and improve reproducibility of scientific results.

earth exchange↗

SIM Planetquest Science and Technology: A Status Report

Optical interferometry will open new vistas for astronomy over the next decade. The Space Interferometry Mission (SIM-PlanetQuest), operating unfettered by the Earth's atmosphere, will offer unprecedented astrometric precision that promises the discovery of Earth-analog extra-solar planets as well as a wealth of important astrophysics. Results from SIM will permit the determination of stellar masses to accuracies of 2% or better for objects ranging from brown dwarfs through main sequence stars to evolved white dwarfs, neutron stars, and black holes. Studies of star clusters will yield age determinations and internal dynamics. Microlensing measurements will present the mass spectrum of the Milky Way internal to the Sun while proper motion surveys will show the Sun's orbital radius and speed. Studies of the Galaxy's halo component and companion dwarf galaxies permit the determination of the Milky Way's mass distribution, including its Dark Matter component and the mass distribution and Dark Matter component of the Local Group. Cosmology benefits from precision (1-2%) determination of distances to Cepheid and RR Lyrae standard candles. The emission mechanism of supermassive black holes will be investigated. Finally, radio and optical celestial reference frames will be tied together by an improvement of two orders of magnitude. Optical interferometers present severe technological challenges. The Jet Propulsion Laboratory, with the support of Lockheed Martin Advanced Technology Center (LM ATC) and Northrop Grumman Space Technology (NGST), has addressed these challenges with a technology development program that is now complete. The requirements for SIM have been satisfied, based on outside peer review, using a series of laboratory tests and appropriate computer simulations: laser metrology systems perform with 10 picometer precision; mechanical vibrations have been controlled to nanometers, demonstrating orders of magnitude disturbance rejection; and knowledge of component positions throughout the whole test assembly has been demonstrated to the required picometer level. Technology transfer to the SIM flight team is now well along.

astrometry↗

Level 1 Processing of MODIS Direct Broadcast Data From Terra

In February 2000, an effort was begun to adapt the Moderate Resolution Imaging Spectroradiometer (MODIS) Level 1 production software to process direct broadcast data. Three Level 1 algorithms have been adapted and packaged for release: Level 1A converts raw (level 0) data into Hierarchical Data Format (HDF), unpacking packets into scans; Geolocation computes geographic information for the data points in the Level 1A; and the Level 1B computes geolocated, calibrated radiances from the Level 1A and Geolocation products. One useful aspect of adapting the production software is the ability to incorporate enhancements contributed by the MODIS Science Team. We have therefore tried to limit changes to the software. However, in order to process the data immediately on receipt, we have taken advantage of a branch in the geolocation software that reads orbit and altitude information from the packets themselves, rather than external ancillary files used in standard production. We have also verified that the algorithms can be run with smaller time increments (2.5 minutes) than the five-minute increments used in production. To make the code easier to build and run, we have simplified directories and build scripts. Also, dependencies on a commercial numerics library have been replaced by public domain software. A version of the adapted code has been released for Silicon Graphics machines running lrix. Perhaps owing to its origin in production, the software is rather CPU-intensive. Consequently, a port to Linux is underway, followed by a version to run on PC clusters, with an eventual goal of running in near-real-time (i.e., process a ten-minute pass in ten minutes).

Lynnes, Christopher↗

Parallel architectures for iterative methods on adaptive, block structured grids

A parallel computer architecture well suited to the solution of partial differential equations in complicated geometries is proposed. Algorithms for partial differential equations contain a great deal of parallelism. But this parallelism can be difficult to exploit, particularly on complex problems. One approach to extraction of this parallelism is the use of special purpose architectures tuned to a given problem class. The architecture proposed here is tuned to boundary value problems on complex domains. An adaptive elliptic algorithm which maps effectively onto the proposed architecture is considered in detail. Two levels of parallelism are exploited by the proposed architecture. First, by making use of the freedom one has in grid generation, one can construct grids which are locally regular, permitting a one to one mapping of grids to systolic style processor arrays, at least over small regions. All local parallelism can be extracted by this approach. Second, though there may be a regular global structure to the grids constructed, there will be parallelism at this level. One approach to finding and exploiting this parallelism is to use an architecture having a number of processor clusters connected by a switching network. The use of such a network creates a highly flexible architecture which automatically configures to the problem being solved.

Gannon, D.↗

From the Cluster Temperature Function to the Mass Function at Low Z

The project consists of three X-ray Multimirror Mission (XMM) observations of the hot, nearby galaxy cluster Triangulum Australis (one pointing to the cluster center and two offsets). The major goal of the project is to obtain an accurate radial gas temperature profile out to large radii and derive the total gravitating mass within the radius of over density 500. A secondary goal is to obtain a detailed two-dimensional gas temperature map of the central region. We have recently received the central pointing (of the two CDROMs received, one contains the central pointing and the other, a zero-exposure dummy observation). The analysis of the data at hand is underway with the aim of deriving the two-dimensional temperature map. At present, we have to limit quantitative analysis to the EPIC-MOS data because of the absence of the necessary calibration products for the Extended Full Frame mode used for EPIC-pn in these observations. We anticipate the arrival of the data from the remaining two offset pointings in the near future, after which we will be able to undertake the total mass derivation. Money from this year's installment of the grant was partially spent on an upgrade of the computer used for the above data analysis, and on travel to a conference that extensively discussed the XMM calibration issues.

Markevitch, Maxim↗

Derived parameters for NGC 6791 from high-metallicity isochrones

We have computed 8, 10, and 12 Gyr isochrones and physically consistent models of zero-age red horizontal branch stars for stellar masses between 0.55 and 1.3 solar mass, all at (fe/H) = +0.15. Comparison to the NGC 6791 BVI photometry of Kaluzny & Udalski (1992) and Montgomery et al. (1994) yields an age of 10.0 +/- 0.5 Gyr at an apparent distance modulus 13.49 less than (m-M)(sub V) less than 13.70. The color offsets required to fit the isochrones, combined with the spectroscopic results of Friel & Janes (1993), imply that the foreground reddening to NGC 6791 lies in the range 0.24 greater than E(B -V) greater than 0.19 with +0.27 less than (Fe/H) less than + 0.44. These results are derived using a technique by which we predict color shifts and apply these to the isochrones to simulate progressively higher metallicities. The zero-age horizontal branch model suggest that the red horizontal branch stars of NGC 6791 have masses approximately less than 0.7 solar mass. The masses are similar to those found for M67 red horizontal branch stars by Tripicco et al. (1993) and for globular cluster red horizontal branch stars, even though the M67 progenitors are approximately 0.2 solar-mass more massive, while the progenitors of globular cluster horizontal branch stars are similarly less massive. This suggests the presence of a mechanism, not strongly dependent on metallicity, which reduces stellar envelopes on the zero-age horizontal branch to a given mass rather than by a given amount.

Dorman, Ben↗

Navier-Stokes computations with finite-rate chemistry for LO2/LH2 rocket engine plume flow studies

Computational fluid dynamics methods have been developed and applied to Space Shuttle Main Engine LO2/LH2 plume flow simulation/analysis of airloading and convective base heating effects on the vehicle at high flight velocities and altitudes. New methods are described which were applied to the simulation of a Return-to-Launch-Site abort where the vehicle would fly briefly at negative angles of attack into its own plume. A simplified two-perfect-gases-mixing approach is used where one gas is the plume and the other is air at 180-deg and 135-deg flight angle of attack. Related research has resulted in real gas multiple-plume interaction methods with finite-rate chemistry described herein which are applied to the same high-altitude-flight conditions of 0 deg angle of attack. Continuing research plans are to study Orbiter wake/plume flows at several Mach numbers and altitudes during ascent and then to merge this model with the Shuttle 'nose-to-tail' aerodynamic and SRB plume models for an overall 'nose-to-plume' capability. These new methods are also applicable to future launch vehicles using clustered-engine LO2/LH2 propulsion.

Dougherty, N. Sam↗