Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data analysis cluster”

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 739 records · Page 41

Predicting Cancer Risk from Ionizing Radiation

The ability to predict cancer risk associated with exposure to low doses of high-LET ionizing radiation (IR) remains a challenge. Epidemiological methods lack the sensitivity and power to provide detailed risk estimates for cancer and ignore individual variance in IR sensitivity. We have hypothesized that DNA repair capacity can be used as a marker to evaluate and differentiate individual radiation sensitivity. More specifically, this work is based on the concept that the combined time-dose dependence of radiation-induced foci (RIF) of p53-binding protein 1 (53BP1) following low-LET exposure contains sufficient information to infer sensitivity to any other LET. Our hypothesis was tested in 15 different mouse strains as well as in primary human immune cells. We first approached individual ionizing radiation sensitivity in a mouse model by culturing primary skin fibroblasts extracted from 76 mice of 15 different genetic backgrounds and exposing them to HZE particles and X-rays. This work is one of the most extensive studies on the kinetics and possible genetic underpinnings of radiation-induced DNA damage and repair. Our results is in agreement with a DNA repair model we previously postulated, where nearby DNA double strand breaks (DSB) in the nucleus are brought together for more efficient repair, leading to RIF clustering. Such mechanism was evidenced by a specific dose and LET dependence of RIF numbers. Briefly, RIF quantification after low-LET X-ray exposure showed an asymptotic saturation for doses between 1 Gy and 4 Gy 4 hours post-irradiation across all 15 strains. The clustering of DSB across all strains also led to more RIF/Gy for lower LET (X-ray and 350 MeV/n Ar) than for higher LET (600 MeV/n Fe) 4 hours post-exposure. Considering the fact that the number of DSB/Gy should be independent of LET, our data suggest there are more DSB in individual RIF as the LET increases. RIF numbers for 24 and 48 hours post-exposure led to the inverse trend, with more remaining RIF/Gy for higher LET (by 600 MeV/n Fe). This result suggests cells have more difficulty resolving RIF from higher LET as they the number DSB/RIF increases. Note that for most conditions, the variance of RIF/Gy was small within individual animals of the same strain and large between strains, suggesting a strong genetics component. Furthermore, we present our preliminary data from an ongoing study on human genetic associations with IR sensitivity. To address the human variability in responses to HZE particle irradiation in a maximally comprehensive manner, we are in the process of collecting and isolating primary blood mononuclear cells from 768 healthy subjects of European descent, 18-75 years of age, 50/50 male/female distribution. We have analyzed 53BP1+ RIF formation as well as oxidative stress and cell death in primary cells from 192 subjects in response to the same HZE particles as used in mice: 600 MeV/n Fe, 350 MeV/n Ar and 350 MeV/n Si, 1.1 and 3 particles/100m2, 4 and 24 hours after irradiation. We will next complete the quantification of HZE particle-induced DNA and cellular damage in the remaining subjects and compare it to their responses to low-LET irradiation. Finally, we will perform GWAS analysis to identify human genomic associations with IR sensitivity and potential targets for biomarker development.

High-LET radiation↗

Phylogenetic and physiological diversity of microorganisms isolated from a deep greenland glacier ice core

We studied a sample from the GISP 2 (Greenland Ice Sheet Project) ice core to determine the diversity and survival of microorganisms trapped in the ice at least 120,000 years ago. Previously, we examined the phylogenetic relationships among 16S ribosomal DNA (rDNA) sequences in a clone library obtained by PCR amplification from genomic DNA extracted from anaerobic enrichments. Here we report the isolation of nearly 800 aerobic organisms that were grouped by morphology and amplified rDNA restriction analysis patterns to select isolates for further study. The phylogenetic analyses of 56 representative rDNA sequences showed that the isolates belonged to four major phylogenetic groups: the high-G+C gram-positives, low-G+C gram-positives, Proteobacteria, and the Cytophaga-Flavobacterium-Bacteroides group. The most abundant and diverse isolates were within the high-G+C gram-positive cluster that had not been represented in the clone library. The Jukes-Cantor evolutionary distance matrix results suggested that at least 7 isolates represent new species within characterized genera and that 49 are different strains of known species. The isolates were further categorized based on the isolation conditions, temperature range for growth, enzyme activity, antibiotic resistance, presence of plasmids, and strain-specific genomic variations. A significant observation with implications for the development of novel and more effective cultivation methods was that preliminary incubation in anaerobic and aerobic liquid prior to plating on agar media greatly increased the recovery of CFU from the ice core sample.

Bacteria/classification/genetics/isolation & purif↗

A Photo Album of Earth Scheduling Landsat 7 Mission Daily Activities

Landsat7 is a member of a new generation of Earth observation satellites. Landsat7 will carry on the mission of the aging Landsat 5 spacecraft by acquiring high resolution, multi-spectral images of the Earth surface for strategic, environmental, commercial, agricultural and civil analysis and research. One of the primary mission goals of Landsat7 is to accumulate and seasonally refresh an archive of global images with full coverage of Earth's landmass, less the central portion of Antarctica. This archive will enable further research into seasonal, annual and long-range trending analysis in such diverse research areas as crop yields, deforestation, population growth, and pollution control, to name just a few. A secondary goal of Landsat7 is to fulfill imaging requests from our international partners in the mission. Landsat7 will transmit raw image data from the spacecraft to 25 ground stations in 20 subscribing countries. Whereas earlier Landsat missions were scheduled manually (as are the majority of current low-orbit satellite missions), the task of manually planning and scheduling Landsat7 mission activities would be overwhelmingly complex when considering the large volume of image requests, the limited resources available, spacecraft instrument limitations, and the limited ground image processing capacity, not to mention avoidance of foul weather systems. The Landsat7 Mission Operation Center (MOC) includes an image scheduler subsystem that is designed to automate the majority of mission planning and scheduling, including selection of the images to be acquired, managing the recording and playback of the images by the spacecraft, scheduling ground station contacts for downlink of images, and generating the spacecraft commands for controlling the imager, recorder, transmitters and antennas. The image scheduler subsystem autonomously generates 90% of the spacecraft commanding with minimal manual intervention. The image scheduler produces a conflict-free schedule for acquiring images of the "best" 250 scenes daily for refreshing the global archive. It then equitably distributes the remaining resources for acquiring up to 430 scenes to satisfy requests by international subscribers. The image scheduler selects candidate scenes based on priority and age of the requests, and predicted cloud cover and sun angle at each scene. It also selects these scenes to avoid instrument constraint violations and maximizes efficiency of resource usage by encouraging acquisition of scenes in clusters. Of particular interest to the mission planners, it produces the resulting schedule in a reasonable time, typically within 15 minutes.

Potter, William↗

Cosmic Infrared Background Fluctuations in Deep Spitzer Infrared Array Camera Images: Data Processing and Analysis

This paper provides a detailed description of the data reduction and analysis procedures that have been employed in our previous studies of spatial fluctuation of the cosmic infrared background (CIB) using deep Spitzer Infrared Array Camera observations. The self-calibration we apply removes a strong instrumental signal from the fluctuations that would otherwise corrupt the results. The procedures and results for masking bright sources and modeling faint sources down to levels set by the instrumental noise are presented. Various tests are performed to demonstrate that the resulting power spectra of these fields are not dominated by instrumental or procedural effects. These tests indicate that the large-scale ([greater, similar]30') fluctuations that remain in the deepest fields are not directly related to the galaxies that are bright enough to be individually detected. We provide the parameterization of these power spectra in terms of separate instrument noise, shot noise, and power-law components. We discuss the relationship between fluctuations measured at different wavelengths and depths, and the relations between constraints on the mean intensity of the CIB and its fluctuation spectrum. Consistent with growing evidence that the [approx]1-5 [mu]m mean intensity of the CIB may not be as far above the integrated emission of resolved galaxies as has been reported in some analyses of DIRBE and IRTS observations, our measurements of spatial fluctuations of the CIB intensity indicate the mean emission from the objects producing the fluctuations is quite low ([greater, similar]1 nW m-2 sr-1 at 3-5 [mu]m), and thus consistent with current [gamma]-ray absorption constraints. The source of the fluctuations may be high-z Population III objects, or a more local component of very low luminosity objects with clustering properties that differ from the resolved galaxies. Finally, we discuss the prospects of the upcoming space-based surveys to directly measure the epochs inhabited by the populations producing these source-subtracted CIB fluctuations, and to isolate the individual fluxes of these populations.

Arendt, Richard↗

Landslides in West Coast Metropolitan Areas: The Role of Extreme Weather Events

Rainfall-induced landslides represent a pervasive issue in areas where extreme rainfall intersects complex terrain. A farsighted management of landslide risk requires assessing how landslide hazard will change in coming decades and thus requires, inter alia, that we understand what rainfall events are most likely to trigger landslides and how global warming will affect the frequency of such weather events. We take advantage of 9 years of landslide occurrence data compiled by collating Google news reports and of a high-resolution satellite-based daily rainfall data to investigate what weather triggers landslide along the West Coast US. We show that, while this landslide compilation cannot provide consistent and widespread monitoring everywhere, it captures enough of the events in the major urban areas that it can be used to identify the relevant relationships between landslides and rainfall events in Puget Sound, the Bay Area, and greater Los Angeles. In all these regions, days that recorded landslides have rainfall distributions that are skewed away from dry and low-rainfall accumulations and towards heavy intensities. However, large daily accumulation is the main driver of enhanced hazard of landslides only in Puget Sound. There, landslide are often clustered in space and time and major events are primarily driven by synoptic scale variability, namely "atmospheric rivers" of high humidity air hitting anywhere along the West Coast, and the interaction of frontal system with the coastal orography. The relationship between landslide occurrences and daily rainfall is less robust in California, where antecedent precipitation (in the case of the Bay area) and the peak intensity of localized downpours at sub-daily time scales (in the case of Los Angeles) are key factors not captured by the same-day accumulations. Accordingly, we suggest that the assessment of future changes in landslide hazard for the entire the West Coast requires consideration of future changes in the occurrence and intensity of atmospheric rivers, in their duration and clustering, and in the occurrence of short-duration (sub-daily) extreme rainfall as well. Major regional landslide events, in which multiple occurrences are recorded in the catalog for the same day, are too rare to allow a statistical characterization of their triggering events, but a case study analysis indicates that a variety of synoptic-scale events can be involved, including not only atmospheric rivers but also broader cold- and warm-front precipitation. That a news-based catalog of landslides is accurate enough to allow the identification of different landslide/ rainfall relationships in the major urban areas along the US West Coast suggests that this technology can potentially be used for other English-language cities and could become an even more powerful tool if expanded to other languages and non-traditional news sources, such as social media.

landslides↗

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗

Characterization and detection of Anopheles vestitipennis and Anopheles punctimacula (Diptera: Culicidae) larval habitats in Belize with field survey and SPOT satellite imagery

Surveys of larval habitats of Anopheles vestitipennis and Anopheles punctimacula were conducted in Belize, Central America. Habitat analysis and classification resulted in delineation of eight habitat types defined by dominant life forms and hydrology. Percent cover of tall dense macrophytes, shrubs, open water, and pH were significantly different between sites with and without An. vestitipennis. For An. punctimacula, percent cover of tall dense macrophytes, trees, detritus, open water, and water depth were significantly different between larvae positive and negative sites. The discriminant function for An. vestitipennis correctly predicted the presence of larvae in 65% of sites and correctly predicted the absence of larvae in 88% of sites. The discriminant function for An. punctimacula correctly predicted 81% of sites for the presence of larvae and 45% for the absence of larvae. Canonical discriminant analysis of the three groups of habitats (An. vestitipennis positive; An. punctimacula positive; all negative) confirmed that while larval habitats of An. punctimacula are clustered in the tree dominated area, larval habitats of An. vestitipennis were found in both tree dominated and tall dense macrophyte dominated environments. The forest larval habitats of An. vestitipennis and An. punctimacula seem to be randomly distributed among different forest types. Both species tend to occur in denser forests with more detritus, shallower water, and slightly higher pH. Classification of dry season (February) SPOT multispectral satellite imagery produced 10 land cover types with the swamp forest and tall dense marsh classes being of particular interest. The accuracy assessment showed that commission errors for the tall, dense marsh and swamp forest appeared to be minor; but omission errors were significant, especially for the swamp forest (perhaps because no swamp forests are flooded in February). This means that where the classification indicates there are An. vestitipennis breeding sites, they probably do exist; but breeding sites in many locations are not identified and could be more abundant than indicated.

NASA Discipline General Space Life Sciences↗

Simulation of Radiation-Induced DNA Damage With the Code RITRACKS

INTRODUCTION DNA damage is one of the most physiologically important effects of ionizing radiation. Clustered DNA damage events, like double-strand breaks (DSBs), have the most notable biological consequences. DNA damage types depend on both the track structure of the radiation and the spatial organization of the DNA. High linear energy transfer (LET) charged nuclei, found in galactic cosmic rays (GCR), are known to produce large numbers of complex DNA damage events. The human genome is packaged into chromatin, which can take on locus-dependent and cell type-dependent spatial conformations that correspond to epigenetic states, such as more open, extended structures in transcriptionally active chromatin. These epigenetic differences can affect DNA break patterns in response to ionizing radiation, potentially creating distinct DNA repair and signaling outcomes across the genome in different cells. MATERIAL AND METHODS The code RITRACKS (Relativistic Ion Tracks), which simulates stochastic radiation track structures and radiation chemistry, was used to model damage on isolated and histone-bound DNA by various types of ions and photons. The changes made to the code to perform radiation-induced DNA damage, and simulation results on single nucleosomes are given in our recent paper. In this work, the DNA building capabilities of RITRACKS have been extended to simulate more complex DNA structures build on the coarse-grain simulation framework meso-WLCsim. This code can sample generic chromatin fiber conformation ensembles based on the geometry of nucleosomes and mechanical properties of DNA. Using RITRACKS, we simulated the fragment length distributions (FLD) of irradiated DNA structures built using the chromatin conformations of WLCsim and obtained results representative of those obtained with Radiation-Induced Correlated Cleavage with sequencing (RICC-Seq) experiments [6]. We have also performed Fe ion and photon irradiations of K562, IMR90, BJ and RPE-1 cells at NSRL to experimentally validate results. Sample processing and data analysis are in progress and any available preliminary results will be discussed. DISCUSSION The recent updates in the code RITRACKS allow the calculation of several quantities such as the DNA damage yield and the FLD. This approach can be used to model epigenetic state-specific chromatin structure parameters to leverage the epigenetic state data available for many human cell types to infer relative DNA damage sensitivity among genomic loci.

I Plante↗

Numerical Propulsion System Simulation: An Overview

The cost of implementing new technology in aerospace propulsion systems is becoming prohibitively expensive and time consuming. One of the main contributors to the high cost and lengthy time is the need to perform many large-scale hardware tests and the inability to integrate all appropriate subsystems early in the design process. The NASA Glenn Research Center is developing the technologies required to enable simulations of full aerospace propulsion systems in sufficient detail to resolve critical design issues early in the design process before hardware is built. This concept, called the Numerical Propulsion System Simulation (NPSS), is focused on the integration of multiple disciplines such as aerodynamics, structures and heat transfer with computing and communication technologies to capture complex physical processes in a timely and cost-effective manner. The vision for NPSS, as illustrated, is to be a "numerical test cell" that enables full engine simulation overnight on cost-effective computing platforms. There are several key elements within NPSS that are required to achieve this capability: 1) clear data interfaces through the development and/or use of data exchange standards, 2) modular and flexible program construction through the use of object-oriented programming, 3) integrated multiple fidelity analysis (zooming) techniques that capture the appropriate physics at the appropriate fidelity for the engine systems, 4) multidisciplinary coupling techniques and finally 5) high performance parallel and distributed computing. The current state of development in these five area focuses on air breathing gas turbine engines and is reported in this paper. However, many of the technologies are generic and can be readily applied to rocket based systems and combined cycles currently being considered for low-cost access-to-space applications. Recent accomplishments include: (1) the development of an industry-standard engine cycle analysis program and plug 'n play architecture, called NPSS Version 1, (2) A full engine simulation that combines a 3D low-pressure subsystem with a 0D high pressure core simulation. This demonstrates the ability to integrate analyses at different levels of detail and to aerodynamically couple components, the fan/booster and low-pressure turbine, through a 3D computational fluid dynamics simulation. (3) Simulation of all of the turbomachinery in a modern turbofan engine on parallel computing platform for rapid and cost-effective execution. This capability can also be used to generate full compressor map, requiring both design and off-design simulation. (4) Three levels of coupling characterize the multidisciplinary analysis under NPSS: loosely coupled, process coupled and tightly coupled. The loosely coupled and process coupled approaches require a common geometry definition to link CAD to analysis tools. The tightly coupled approach is currently validating the use of arbitrary Lagrangian/Eulerian formulation for rotating turbomachinery. The validation includes both centrifugal and axial compression systems. The results of the validation will be reported in the paper. (5) The demonstration of significant computing cost/performance reduction for turbine engine applications using PC clusters. The NPSS Project is supported under the NASA High Performance Computing and Communications Program.

Lytle, John K.↗

Spatial Correlation Function of the Chandra Selected Active Galactic Nuclei

We present the spatial correlation function analysis of non-stellar X-ray point sources in the Chandra Large Area Synoptic X-ray Survey of Lockman Hole Northwest (CLASXS). Our 9 ACIS-I fields cover a contiguous solid angle of 0.4 deg(exp 2) and reach a depth of 3 x 10(exp -15) erg/square cm/s in the 2-8 keV band. We supplement our analysis with data from the Chandra Deep Field North (CDFN). The addition of this field allows better probe of the correlation function at small scales. A total of 233 and 252 sources with spectroscopic information are used in the study of the CLASXS and CDFN fields respectively. We calculate both redshift-space and projected correlation functions in co-moving coordinates, averaged over the redshift range of 0.1 < z < 3.0, for both CLASXS and CDFN fields for a standard cosmology with Omega(sub Lambda) = 0.73,Omega(sub M) = 0.27, and h = 0.71 (H(sub 0) = 100h km/s Mpc(exp -1). The correlation function for the CLASXS field over scales of 3 Mpc< s < 200 Mpc can be modeled as a power-law of the form xi(s) = (S/SO)(exp - gamma), with gamma = 1.6(sup +0.4 sub -0.3) and S(sub o) = 8.0(sup +.14 sub -1.5) Mpc. The redshift-space correlation function for CDFN on scales of 1 Mpc< s < 100 Mpc is found to have a similar correlation length so = 8.55(sup +0.74 sub -0.74) Mpc, but a shallower slope (gamma = 1.3 +/- 0.1). The real-space correlation functions derived from the projected correlation functions, are found to be tau(sub 0 = 8.1(sup +1.2 sub -2.2) Mpc, and gamma = 2.1 +/- 0.5 for the CLASXS field, and tau(sub 0) = 5.8(sup +.1.0 sub -1.5) Mpc, gamma = 1.38(sup +0.12 sub -0.14 for the CDFN field. By comparing the real- and redshift-space correlation functions in the combined CLASXS and CDFN samples, we are able to estimate the redshift distortion parameter Beta = 0.4 +/- 0.2 at an effective redshift z = 0.94. We compare the correlation functions for hard and soft spectra sources in the CLASXS field and find no significant difference between the two groups. We have also found that the correlation between X-ray luminosity and clustering amplitude is weak, which, however, is fully consistent with the expectation using the simplest relations between X-ray luminosity, black hole mass, and dark halo mass. We study the evolution of the AGN clustering by dividing the samples into 4 redshift bins over 0.1 Mpc< z <3.0 Mpc. We find a very mild evolution in the clustering amplitude, which show the same evolution trend found in optically selected quasars in the 2dF survey. We estimate the evolution of the bias, and find that the bias increases rapidly with redshift (b(z = 0.45) = 0.95 +/- 0.15 and b(z = 2.07) = 3.03 +/- 0.83): The typical mass of the dark matter halo derived from the bias estimates show little change with redshift. The average halo mass is found to be log (M(sub halo)/M(sun))approximates 12.1. Subject headings: cosmology: observations - large-scale structure of the universe - x-rays: diffuse background - galaxies: nuclei

Yang, Y.↗

The Star Formation History of Orion and its Environs

During the period of performance, we have obtained the following observations and carried out the analysis of the Orion associations itemized below. Optical Photometry: The Quest camera that had been damaged by lightning on Sept. 2001 was reconfigured to use 10 chips and only 3 filters instead of the original 4 but retaining its spatial coverage of 2.3deg wide in declination. During the Winter and Fall 2002 observing campaigns, we obtained multifilter photometry using the Quest camera on the 1mt Schmidt of the Venezuelan National Observatory, to continue the survey of the Orion clouds. In the Winter 2002 season, we obtained photometry for 10 scans at DEC=+5, 7 scans at DEC=-5. Several additional scans at DEC=-3 were obtained during the Fall 2002 observing campaign. Slit spectroscopy of bright candidates: During 2002 spectra for some 235 candidate PMS stars were obtained with the FAST spectrograph at the SAO 1.5m telescope in Mt. Hopkins. These objects, all brighter than V=16, are located in two strips centered at DEC=+l and -5, and spanning from RA=5h-6h. The spectra are being analyzed. Multifiber spectroscopy: Spectroscopy for our fainter (V less than or equal to 16) candidates has not kept pace with our spectroscopy of brighter candidates, because of the delay in the commissioning of the Hectospec and Hectoechelle multi-object spectrographs at the new 6.5m MMT. We obtained some data with the Hydra spectrograph on the WIYN telescope: Feb 1-4, 2002. Time granted to the project: "A Large Scale Survey of the Low-Mass Young Population in the Orion OB1 Association" (P.I.: C. Briceno). Three half-nights (we used the first half of the night, the second half of the night was used with the Mini-Mosaic imager by another team). Only the first half night was clear. We observed 3 fields centered in a 2.3 deg wide strip centered at DEC=-1 and spanning RA=Sh to 6h. A total of 250 spectra were obtained for the same number of candidate PMS stars. Nov. 13-15, 2002. Time granted to the project "Cluster survey of protoplanetary disk evolution", (P.I.: A. Sicilia-Aguilar). Clear weather throughout. We obtained spectra for some 200 candidate PMS stars in 3 fields located in a 2.3 deg wide strip centered at DEC=+l and spanning from RA=5h20m to 5h35m. Data for both observing runs is partially reduced and is being analyzed. U photometry: During Nov.29-Dec.4, 2002 we were granted time with the 4-shooter CCD Mosaic Camera at the SAO 1.2m telescope, to obtain U-band photometry of a subset of the newly identified T Tauri stars in the strip centered at DEC=-1. This sample is composed of strong H(alpha) emitting PMS stars (Classical T Tauri stars) located mostly in the Orion OB lb association, around the Orion Belt area. Our data will be combined with our calibration of the U-band excess (Gullbring et al. 1998) to derive mass accretion rates. Because of bad weather we observed only 20 stars during 1.5 clear nights. Near IR photometry: During Dec. 14-19, 2002 we used the IR Camera on the SAO 1.2m telescope, to obtain L-band photometry of a set of 17 the newly identified Orion OB 1a and 1b stars, in order to look for IR excess emission coming from the hotter inner parts of circumstellar disks. Again, mediocre weather prevented a larger set to be observed. These stars were also observed with the OSCIR mid-IR camera on Gemini North during later 2001. Our goal is to construct SEDs from the visual out to the mid-IR (by combining our optical/IR data with 2MASS JHK magnitudes) in order to explore disk dissipation at the critical ages of 5-10 Myr spanned by these stars.

Calvet, Nuria↗

XMM-Newton and Chandra Observations of the Galaxy Group NGC 5044: Evidence for Limited Multiphase Hot Gas - 1

Using new XMM and Chandra observations, we present an analysis of the temperature structure of the hot gas within a radius of 100 kpc of the bright nearby galaxy group NGC 5044. A spectral deprojection analysis of data extracted from circular annuli reveals that a two-temperature model (2T) of the hot gas is favored over single-phase or cooling flow (M = 4.5 +/- 0.2 solar mass/yr) models within the central approx.30 kpc. Alternatively, the data can be fitted equally well if the temperature within each spherical shell varies continuously from approx.T(sub h) to T(sub c) approx. T(sub h)/2, but no lower. The high spatial resolution of the Chandra data allows us to determine that the temperature excursion T(sub h) approaches T(sub c) required in each shell exceeds the temperature range between the boundaries of the same shell in the best-fitting single-phase model. This is strong evidence for a multiphase gas having a limited temperature range. We do not find any evidence that azimuthal temperature variations within each annulus on the sky can account for the range in temperatures within each shell. We provide a detailed investigation of the systematic errors on the derived spectral models considering the effects of calibration, plasma codes, bandwidth, variable NH, and background rate. We find that the RGS gratings and the EPIC and ACIS CCDs give fully consistent results when the same models are fitted over the same energy ranges for each instrument. The cooler component of the 2T model has a temperature (T(sub c) approx. 0.7 keV) similar to the kinetic temperature of the stars. The hot phase has a temperature (T(sub h) approx. 1.4 keV) characteristic of the virial temperature of the solar mass halo expected in the NGC 5044 group. However, in view of the morphological disturbances and X-ray holes visible in the Chandra image within R approx. equals 10 kpc, bubbles of gas heated to approx.T(sub h) in this region may be formed by intermittent AGN feedback. Some additional heating at larger radii may be associated with the evolution of the cold front near R approx. 50 kpc, as suggested by the sharp edge in the EPIC images. Subject headings: cooling flows - galaxies: formation - galaxies: halos - intergalactic medium - X-rays: galaxies: clusters

Buote, David A.↗

NASA Taxonomies for Searching Problem Reports and FMEAs

Many types of hazard and risk analyses are used during the life cycle of complex systems, including Failure Modes and Effects Analysis (FMEA), Hazard Analysis, Fault Tree and Event Tree Analysis, Probabilistic Risk Assessment, Reliability Analysis and analysis of Problem Reporting and Corrective Action (PRACA) databases. The success of these methods depends on the availability of input data and the analysts knowledge. Standard nomenclature can increase the reusability of hazard, risk and problem data. When nomenclature in the source texts is not standard, taxonomies with mapping words (sets of rough synonyms) can be combined with semantic search to identify items and tag them with metadata based on a rich standard nomenclature. Semantic search uses word meanings in the context of parsed phrases to find matches. The NASA taxonomies provide the word meanings. Spacecraft taxonomies and ontologies (generalization hierarchies with attributes and relationships, based on terms meanings) are being developed for types of subsystems, functions, entities, hazards and failures. The ontologies are broad and general, covering hardware, software and human systems. Semantic search of Space Station texts was used to validate and extend the taxonomies. The taxonomies have also been used to extract system connectivity (interaction) models and functions from requirements text. Now the Reconciler semantic search tool and the taxonomies are being applied to improve search in the Space Shuttle PRACA database, to discover recurring patterns of failure. Usual methods of string search and keyword search fall short because the entries are terse and have numerous shortcuts (irregular abbreviations, nonstandard acronyms, cryptic codes) and modifier words cannot be used in sentence context to refine the search. The limited and fixed FMEA categories associated with the entries do not make the fine distinctions needed in the search. The approach assigns PRACA report titles to problem classes in the taxonomy. Each ontology class includes mapping words - near-synonyms naming different manifestations of that problem class. The mapping words for Problems, Entities and Functions are converted to a canonical form plus any of a small set of modifier words (e.g. non-uniformity NOT + UNIFORM.) The report titles are parsed as sentences if possible, or treated as a flat sequence of word tokens if parsing fails. When canonical forms in the title match mapping words, the PRACA entry is associated with the corresponding Problem, Entity or Function in the ontology. The user can search for types of failures associated with types of equipment, clustering by type of problem (e.g., all bearings found with problems of being uneven: rough, irregular, gritty ). The results could also be used for tagging PRACA report entries with rich metadata. This approach could also be applied to searching and tagging failure modes, failure effects and mitigations in FMEAs. In the pilot work, parsing 52K+ truncated titles (the test cases that were available), has resulted in identification of both a type of equipment and type of problem in about 75% of the cases. The results are displayed in a manner analogous to Google search results. The effort has also led to the enrichment of the taxonomy, adding some new categories and many new mapping words. Further work would make enhancements that have been identified for improving the clustering and further reducing the false alarm rate. (In searching for recurring problems, good clustering is more important than reducing false alarms). Searching complete PRACA reports should lead to immediate improvement.

Malin, Jane T.↗

The Star Formation History of Orion and its Environs

During this period of performance, we have obtained the following observations and carried out the analysis of the Orion associations itemized below. 1. Quest Optical Photometry: The photometric survey carries out with the Quest camera is finished. The strips at declinations +1 degree and -1 degree have already been processed. Photometry for a total of around 400,000 stars was obtained in these two strips with the Quest camera pipeline. Around 24,000 variables were picked out of this total by our variability software. Of these, around 2,500 stars fall above the main sequence and so were picked as candidates for spectroscopic follow-up. 2. Slit spectroscopy of bright candidates: Spectra for some 800 candidate PMS stars were obtained with the FAST spectrograph at the SAO 1.5m telescope in Mt. Hopkins. The spectra are being analyzed; 300 stars have been confirmed as young. 3. Multifiber spectroscopy: The first test of the multifiber spectrograph Hectoechelle were carried out in December 2003. One field of the Orion Nebula Cluster was observed with Hectochelle at three wavelength settings. A total of 157 spectra were obtained. Of these, 63 stars have been confirmed as Classical T Tauri stars, and 36 additional stars need further follow up. A paper is in preparation. 4. UBVRI photometry: We were granted time with the 4-shooter CCD Mosaic Camera at the SAO 1.2m telescope, to obtain UBVRI photometry of a subset of 53 newly identified T Tauri stars in the strips centered at DEC=-1 and +l. This sample is composed of strong Halpha emitting PMS stars (Classical T Tauri stars) located mostly in the Orion OB l b association, around the Orion Belt area. We have estimated mass accretion rates for 22 for these stars using the U photometry and the calibration of Gullbring et al. (1998), and found it to be similar to that of young stars in associations of similar age. 5. Near and mid-IR photometry: During the winter of 2003, we used the IR Camera on the SAO 1.2m telescope, to obtain L- band photometry of a set of the newly identified Orion OB l a and l b stars, in order to look for IR excess emission coming from the hotter inner parts of circumstellar disks. These data is being combined with the GEMINI/OSCIR 10 micron photometry obtained during 2002 and the mass accretion rates to determine the properties of disks in Ori lb. A paper is presently on preparation.

Calvet, Nuria↗

Atacama Cosmology Telescope: Combined kinematic and thermal Sunyaev-Zel’dovich measurements from BOSS CMASS and LOWZ halos

The scattering of cosmic microwave background (CMB) photons off the free-electron gas in galaxies and clusters leaves detectable imprints on high resolution CMB maps: the thermal and kinematic Sunyaev-Zel’dovich effects (tSZ and kSZ respectively). We use combined microwave maps from the Atacama Cosmology Telescope DR5 and Planck in combination with the CMASS (mean redshifthzi¼0.55and host halo masshMviri¼3×1013M⊙) and LOWZ (hzi¼0.31,hMviri¼5×1013M⊙) galaxy catalogs from the Baryon Oscillation Spectroscopic Survey (BOSS DR10 and DR12), to study the gas associated with these galaxy groups. Using individual reconstructed velocities, we perform a stacking analysis and reject the no-kSZ hypothes is at 6.5σ, the highest significance to date. This directly translates into a measurement of the electron number density profile, and thus of the gas density profile. Despite the limited signal to noise, the measurement shows at high significance that the gas density profile is more extended than the dark matter density profile, for any reasonable baryon abundance (formally>90σfor the cosmic baryon abundance). We simultaneously measure the tSZ signal, i.e., the electron thermal pressure profile of the same CMASS objects, and reject theno-tSZ hypothesis at10σ. We combine tSZ and kSZ measurements to estimate the electron temperature to20% precision in several aperture bins, and find it comparable to the virial temperature. In a companion paper, we analyze these measurements to constrain the gas thermodynamics and the properties of feedback inside galaxy groups. We present the corresponding LOWZ measurements in this paper, ruling out a null kSZ (tSZ)signal at 2.9ð13.9Þσ, and leave their interpretation to future work. This paper and the companion paper demonstrate that current CMB experiments can detect and resolve gas profiles in low mass halos and at high redshifts, which are the most sensitive to feedback in galaxy formation and the most difficult to measure any other way. They will be a crucial input to cosmological hydrodynamical simulations, thus improving our understanding of galaxy formation. These precise gas profiles are already sufficient to reduce the main limiting theoretical systematic in galaxy-galaxy lensing: baryonic uncertainties. Future such measurements will thus unleash the statistical power of weak lensing from the Rubin, Euclid and Roman observatories. Our stacking software Thumb Stack is publicly available and directly applicable to future Simons Observatory andCMB-S4 data.

Emmanuel Schaan↗

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively. The Mach Number distribution shows small yet distinct differences between the two cases such as locations of shocks/shock reflections and a slightly different impingement point on the wall of the diffuser from the expansion at the exit of the nozzle. Similarly the temperature distribution indicates different flow recirculation patterns in the test cell. Both cases capture all the essential flow phenomena such as the shock-boundary layer interaction, plume expansion, expansion of the first stage ejectors, mixing between the engine plume and the first stage ejector flow and pressurization due to the first stage ejectors. The final paper will discuss thermal loads on the walls of the diffuser and cooling mechanisms investigated.

Allgood, Daniel C.↗

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively.

Allgood, Daniel C.↗

Using XR for Improving Scientific Discovery With Numerical Weather Models

Earth science (ES) digital twins will help us understand the complex interactions and interrelationships that make up our Earth system and the impacts of earth science phenomena on it. Our work addresses two underdeveloped areas in current ES digital twin work: improving the understanding and interaction with ES model outputs by using Virtual and Mixed Reality (XR) tools and improving the non-intuitive mapping of continuous ES natural phenomena to gridded reference frames in current numerical models. Traditionally, scientists working on ES view and analyze the results of calculated or measured observables with static 1-dimensional (1D), 2D or 3D plots displayed on flat computer screens or paper. Using such limited mediums, it can be very difficult to identify, track and understand the evolution of key features due to poor viewing angles and the nature of flat computer screens. In addition, numerical models, such as the NASA Goddard Earth Observing System (GEOS) ES model, are almost exclusively formulated, visualized and analyzed in an Eulerian reference frame with fixed grid points in space and time. However, ES phenomena such as convective clouds, hurricanes and wildfire smoke plumes are visualized and analyzed in a Lagrangian reference frame: therefore it is often difficult and unnatural to understand these phenomena in relation to each other, visualized either in an Eulerian or Lagrangian context. In 3D visualizations, data generally takes one of three forms: gridded (e.g., voxelized) data, where space is divided into regions; point clouds, where data is represented as a set of points; and meshes, where objects are rendered as surfaces composed of small polygons (usually triangles). A gridded, Eulerian reference frame has been the default representation for the 2D visual analysis of atmospheric data in part because the numerical methods used to generate atmospheric model data in the first place use a gridded approach, with equations defining the relationships between the physical variables in each of a grid's cells across successive timesteps. In our work, we are particularly interested in data from GEOS. Another reason why gridded representations tend to be used for visualizing data from such models is because trajectories are difficult to interpret from representations on 2D surfaces, due to line-of-sight ambiguity. Instead of a fixed grid from GEOS, we embed a trajectory model to simulate particles' movement throughout a GEOS run. We then ingest these particle trajectories as animated point clouds with a NASA open source XR toolkit, the Mixed Reality Exploration Toolkit (MRET), and merge GEOS data with ES phenomena data onto one combined visualization that the user can intuitively interact with. Efficient rendering of arbitrarily large point clouds is an ongoing challenge being addressed by the computer science community, with the GPU-based optimizations and efficient GPU memory utilization a common theme of recent advances, especially for XR, where sustained high frame rate is mandatory to save the user from suffering due to simulation sickness. In this work, we describe and evaluate our progress in choosing and implementing appropriate methods for rendering arbitrarily large point clouds within MRET for XR. While tracking the XR headset enables the immersion of a user within a 3D scene of a data visualization, tracking of XR handheld controllers or user’s hands enables us to implement intuitive user interactions with the visualized datasets. Conventional tools require a user working with an ES visualization to conduct many interactions to commit their intended selections or manipulations with a visualized dataset; for example to specify a set of points in 3D space. Doing so in a 2D flat screen interface has traditionally required specifying a set of points in three distinct 2D coordinate systems (XY, XZ, and YZ), which is cumbersome. In other scientific domains, it has been shown that specifying or selecting a location or volume in XR using handheld controllers or tracked hands allows for greater speed and accuracy. We anticipate the same will hold true for atmospheric data, and we will share initial results of measuring the utility of such an interface. Notably, as the data being visualized is generated by GEOS as a prediction based on initial conditions, an intended application of our tool is to serve as part of an iterative feedback loop. Through XR, a scientist will review and manipulate a GEOS model run, modifying the conditions as needed to do subsequent runs of GEOS. Thereby, XR-based improvements to speed and accuracy of 3D tagging of points minimizes the effort required by both the scientist and the computer cluster conducting the necessary calculations.

Thomas Grubb↗