Search NASA⌕ Search

SEARCH · Search NASA

Results for “COUNTING RATE COMPUTER”

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 145 records · Page 8

Correlation of the cosmic-ray intensity with solar-terrestrial parameters

Results are presented of numerical simulations of cosmic ray drift due to perturbations in the interplanetary medium affecting solar-terrestrial parameters. Computations were performed for a model spherical heliosphere of radius 25 AU, with an interplanetary magnetic field in the shape of an Archimedean spiral with magnitude of 5 gamma at 1 AU, corresponding to a wind velocity of 400 km/sec, and diffusion coefficients corresponding to the assumption that the parallel mean free path is somewhat larger han the cyclotron radius at earth orbit, and its spatial dependence is proportional to the gyroradius. It is found that perturbations in the heliospheric equatorial regions affect galactic cosmic rays more effectively when the Northern Hemisphere heliospheric magnetic field is directed inward rather than outward, with perturbations in the parallel drift coefficient more effective than perturbations in solar wind velocity or perpendicular drift velocity. Results are consistent with the observed variation of the correlation between the geomagnetic aa index and the neutron monitor counting rate on Mt. Washington with the solar cycle.

Jokipii, J. R.↗

Dynamic Measurement of Temperature, Velocity, and Density in Hot Jets Using Rayleigh Scattering

A molecular Rayleigh scattering technique is utilized to measure gas temperature, velocity, and density in unseeded gas flows at sampling rates up to 10 kHz, providing fluctuation information up to 5 kHz based on the Nyquist theorem. A high-power continuous-wave laser beam is focused at a point in an air flow field and Rayleigh scattered light is collected and fiber-optically transmitted to a Fabry-Perot interferometer for spectral analysis. Photomultiplier tubes operated in the photon counting mode allow high-frequency sampling of the total signal level and the circular interference pattern to provide dynamic density, temperature, and velocity measurements. Mean and root mean square velocity, temperature, and density, as well as power spectral density calculations, are presented for measurements in a hydrogen-combustor heated jet facility with a 50.8-mm diameter nozzle at NASA John H. Glenn Research Center at Lewis Field. The Rayleigh measurements are compared with particle image velocimetry data and computational fluid dynamics predictions. This technique is aimed at aeronautics research related to identifying noise sources in free jets, as well as applications in supersonic and hypersonic flows where measurement of flow properties, including mass flux, is required in the presence of shocks and ionization occurrence.

Mielke, Amy F.↗

Gamma and Neutron Measurement and Modeling of Irradiated TRISO Fuel

Given the unique characteristics of the PBR fuel cycle, both gamma and neutron measurements are expected to play important roles in performing and maintaining nuclear material control and accounting for spent pebbles to safeguard the fuel cycle. Given the lack of irradiated pebbles in the US, a variety of irradiated TRISO fuel samples with wide ranges of burnups and cooling times available at ORNL were used in this work. A large number of gamma and neutron measurements have been performed on these samples to collect data to test the various detectors and to benchmark the computer models to simulate the depletion and decay of the fuel and the measurements themselves. Two neutron detectors, including a custom-made detector and the Very High-Performance Neutron Multiplicity Counting, were used to measure the neutrons emitted by these TRISO samples. Three gamma spectrometry detectors, including an HPGe and the M400 CZT detector, were used to measure gamma-ray emissions from these samples. The M400 was recently adopted by the IAEA for fresh uranium measurements, but it was tested for spent fuel measurements prior to this project. Detailed MCNP models were developed to simulate these neutron and gamma measurements. Some GADRAS models were also developed to cross check the MCNP models for the gamma measurements. It was found challenging to perform neutron measurements in the hot cell due to the high background counts. Close agreements were observed between the simulated and measured neutron count rates in both detectors’ measurements of californium calibration sources. Both the HPGe and M400 detectors were able to measure the 604 and 662 keV peaks from these samples, which are the two most important peaks used to infer fuel burnup. Although the M400 detector did not have nearly good energy resolution and did not detect some of the minor peaks as the HPGe detector, it was found to be capable of handling significantly higher dose rates than HPGe. Given the complexities in the TRISO samples (e.g., different samples sizes) and uncertainties in the alignments between the detector and the TRISO fuel inside the containers, large scatters were found between the peak area rates and the samples’ burnups. However, the 604/662 peak ratios were found to trend well with the samples’ burnups among most samples in both measured and simulated results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Are Stronger North-Atlantic Southwesterlies the Forcing to the Late-Winter Warming in Europe?

We examine a possible mechanism leading to late-winter warming and thus to an early spring in Europe. From the NCEP Reanalysis, we extract for the years 1948-1999 ocean-surface winds over the eastern North Atlantic, and air temperatures at the surface, T(sub s), and at the 500 mb level, T(sub 500), in late-winter and spring. T(sub s) is extracted at six European locations, all at 50.5 N, ranging in longitude from 1.9 E (northeastern France) to 26.2 E (Ukraine). To quantify the advection of maritime air into Europe, we evaluate for 3-pentad groups the Index I(sub na) of the southwesterlies at 45 N; 20 W: I(sub na) is the average wind speed at this point if the direction is from the quadrant 180-270 deg (when the direction is different, the contribution counts as zero). In late winter correlations C(sub it) between the Index I(sub na) and the temperature T(sub s) are substantial, up to the 0.6 level, in western Europe (but weaker correlations for Poland and Ukraine). C(sub it) drops sharply by mid-March, taking occasionally negative values subsequently. This drop in C(sub it) indicates that maritime air advection is no longer associated closely with the surface-air warming, the role of immolation becomes important, and thus the drop in C(sub it) marks the arrival of spring. Correlations C(sub i delta) between I(sub na) and our lapse-rate parameter delta, the difference between T(sub s) and T(sub 500), indicate that the flow of warm maritime-air from the North Atlantic into this 'corridor' at 50.5 N is predominantly at low tropospheric level. By computing the best linear fit to I(sub na) and T(sub s), the trends for the period 1948-1999 are evaluated. The trends are appreciable in the second half of February and the first half of March. Our 3-pentad analysis points to the interval from mid-February to mid-March as the end-of-winter period in which the southwesterlies over the eastern North Atlantic become stronger and the surface-air temperatures in Europe rise markedly, the lapse rate becomes steeper, and concurrently the longitudinal temperature gradient between the Somme (France) and the Oder (Germany/Poland border) is reduced by 0.8 C, that is, by 20% of its 1948 value. Our thesis, that the observed late-winter warming and the corollary advancement of spring in Europe resulted at least in part from stronger southwesterlies over the North Atlantic, merits further investigations.

Otterman, J.↗

Harvesting the low-hanging fruit of high energy savings -- Virtual Occupancy using Wi-Fi Data

Approximately 20% of primary energy consumed in the U.S. is attributed to HVAC use. Ideally, HVAC operation would be driven by actual building occupancy, but lack of reliable occupancy information often results in the use of conservative static schedules. This disparity is even more pronounced in a college campus, where the function of each space differs by building (classrooms, offices, libraries) and the class schedules change frequently -- every semester, day of week, and hour. While several research papers propose the use of counts of the Wi-Fi connections (e.g., phones, computers) as a proxy for occupancy, few real-world implementations exist. This paper describes the development and deployment of an open-source Wi-Fi-to-Occupancy software library in 65 buildings of a college campus, and the planned integration with the building energy management and control system at the building scale. Over a year of Wi-Fi data was gathered into distinct academic periods, including fall and spring semester, academic breaks, and summer sessions. Patterns such as students moving between classrooms, closing laptops before exams, etc., can be visualized from the data. Approximating occupancy from Wi-Fi data presents challenges which we address in this project -- for example, identifying static devices, or estimating the ratio of devices per person. Utilizing real-time occupancy data to inform optimal HVAC schedules and ventilation rates creates the potential to identify and reduce energy waste. Other potential applications include forecasting occupancy, and using Wi-Fi data to predict peak demands. Finally, the paper discusses how to easily scale these tools to other buildings.

Pritoni, Marco↗

Negative Muon Count for the FNAL Muon g-2 Experiment

The primary goal of the FNAL Muon g-2 Experiment, located in Batavia, IL, USA, is to measure the anomalous magnetic dipole moment, a μ = (g-2)/2, of a muon to a precision of 0.14 ppm. Initial results from the experiment confirm a discrepancy between the experimental value of the muon anomaly compared with the theoretical value first discovered by the Brookhaven Muon g-2 Experiment in 2001. This discrepancy hints at new physics. However, running the experiment with negative muons is essential in order to verify the experimental value of the muon anomaly. As neither experiments nor simulation studies pertaining to the muon anomaly have been performed with negative muons, the goal of this work is to estimate the rates of negative muons arriving at the FNAL Muon g-2 experiment so that scientists can estimate the time required to run the experiment with negative muons. In order to do this, we performed several simulations of the Muon Campus, which delivers muons to the FNAL Muon g-2 Experiment, using the G4Beamline v. 3.06 at the National Energy Research Computing Center. We performed separate simulations for different sections of the muon campus, and simulations were repeated twice, once for negative muons and once for positive muons. An initial 10 9 protons were shot on target for both the negative and positive muon simulations. We found that the rate of negative muons arriving at the experiment was 0.549 times the rate of positive muons arriving at the experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

An Evaluation and Recommendations for Further CFD Research Based on the NASA Common Research Model (CRM) Analysis from the AIAA Drag Prediction Workshop (DPW) Series

The six Drag Prediction Workshops (DPW) have drawn a total of 125 participants submitting 225 data entries. Over half of these occurred in the last three workshops dealing with the NASA Common Research Model (CRM). Participants have ranged from North America, Europe, Asia, and South America, representing Government, Industry, Academia, and Commercial/Vendors. The DPW Series has generated a tremendous amount of CFD data, freely available to the public domain. A continuous improvement of the results as measured by: the reduction of the spread of drag predictions for the specified “nominal” cruise point, accuracy of drag increment, and the angle-of-attack sweep behavior, has been noted. Although it has been difficult separating out the parts and pieces that contribute to errors in CFD drag predictions, part of the improvements over the years can be attributed to grid resolution. The scatter in predicted drag in the continuum has been reduced dramatically from plus or minus 50 counts in DPW-I to around plus or minus 5 counts for DPW-V and - VI. During the six workshops, the grid size has grown dramatically. The average size of the medium wing-body meshes in DPW-I through DPW-IV have been 3.2, 5.4, 7.8 and 10.9 million, respectively. This represents a growth rate of approximately 17 percent per year during the eight years between DPW-I and DPW-IV. For DPW-VI, this number has grown from 25 to 50 million points for the various families of grids available. The finest level grids have increased steadily, from just over 3 million grid points in DPW-I to 225 million for the WBNP (Wing-Body-Nacelle-Pylon) in DPW-VI. While progress has been made, persistent problems with grid generation and separated flow prediction continue. The experiences and recommendations presented here are consistent with those espoused in the “CFD Vision 2030 Study: A Path to Revolutionary Computational Aerosciences,” NASA/CR-2014-218178.

Tinoco, Edward N.↗

Rapid neutron and gamma-ray source localization using machine learning

Rapid localization of radiation sources is critical for applications including nuclear emergency response, safeguards, and security. However, conventional imaging systems such as neutron scatter cameras and Compton cameras depend on rare coincidence events, which often result in long acquisition times. In this work, we address the challenge of rapid source localization by developing a machine learning approach to predict the direction of a single radiation source using only count rates from an array of neutron and gamma-ray detectors. The proposed model is a fully connected neural network (FCNN) trained using Monte Carlo simulation data from a 252 Cf source. The model hyperparameters are optimized with a small set of routine 252 Cf measurements. We benchmarked the performance of the trained and optimized machine learning model using additional 252 Cf , 137 Cs , and PuBe measurements under laboratory conditions with varying source-detector configurations. For these measurements, the machine learning model achieved a mean localization error smaller than 30° with 3 x 10 3 system counts, corresponding to 8 s measurement time for the imaging system used in this work. In this low-statistics regime, the method outperformed traditional scatter-based imaging by more than 75% in localization accuracy for the evaluated measurement configurations. These results demonstrate that a machine learning-based approach can significantly reduce the time required for accurate single-source localization, providing a robust and computationally efficient alternative to traditional imaging systems in time-critical nuclear security and emergency response scenarios.

Gamma-ray imaging↗

The Apollo Sample Suite: 50 Years of Solar System Insight

The Apollo program was undoubtable a crowning achievement in human history. In addition to the obvious cultural significance, scientific results from the Apollo program had a lasting impression on a range of scientific fields, none more so that the effect the samples had on the fields of geology and cosmochemistry. The six Apollo missions collected 382 kg of rock, regolith, and core samples from geologically diverse locations on the Moon. In the nearly 50 years since the first samples were returned, there have been over 3000 different requests for samples, each yielding insights into fields as disparate as biology, medicine, astronomy, engineering, material science, and of course geology. Early studies of the Apollo samples revealed primary insights into the origin and evolution of the Moon, and of the Earth-Moon system, but the results also had implications for bodies throughout the solar system, e.g., defining crater counting rates. Over the decades, continued study of the Apollo samples by new generations of scientists using new instruments have continued to yield significant new discoveries, including the presence of endogenous water in the Moon and the possible presence of a lunar cataclysm, that in turn has contributed to new models of solar system formation and evolution. The Apollo samples have often been used as a proxy for studying other bodies like Mercury or asteroids. The Apollo samples have also directly contributed to the interpretation of remotely sensed data sets, including their use as ground truth for both Clementine and Lunar Prospector global geochemical maps. Despite the Apollo samples being a static collection, recent efforts will ensure that investigators continue to have access to new samples. For example, there was a recent solicitation for study of previously unopened Apollo samples in vacuum-sealed containers, as well as new access to samples stored frozen or in a He atmosphere. Similarly, the use of X-ray computed tomography as part of the curation process is identifying new clasts within polymict breccias that are available for study. Finally, the MoonDB project is putting all previously published lunar geochemical analyses into a searchable database, which should facilitate new investigations.

Zeigler, Ryan↗

High resolution time interval counter

In recent years, we have developed two types of high resolution, multi-channel time interval counters. In the NIST two-way time transfer MODEM application, the counter is designed for operating primarily in the interrupt-driven mode, with 3 start channels and 3 stop channels. The intended start and stop signals are 1 PPS, although other frequencies can also be applied to start and stop the count. The time interval counters used in the NIST Frequency Measurement and Analysis System are implemented with 7 start channels and 7 stop channels. Four of the 7 start channels are devoted to the frequencies of 1 MHz, 5 MHz or 10 MHz, while triggering signals to all other start and stop channels can range from 1 PPS to 100 kHz. Time interval interpolation plays a key role in achieving the high resolution time interval measurements for both counters. With a 10 MHz time base, both counters demonstrate a single-shot resolution of better than 40 ps, and a stability of better than 5 x 10(exp -12) (sigma(sub chi)(tau)) after self test of 1000 seconds). The maximum rate of time interval measurements (with no dead time) is 1.0 kHz for the counter used in the MODEM application and is 2.0 kHz for the counter used in the Frequency Measurement and Analysis System. The counters are implemented as plug-in units for an AT-compatible personal computer. This configuration provides an efficient way of using a computer not only to control and operate the counters, but also to store and process measured data.

Zhang, Victor S.↗

Single-photon detection using high-temperature cuprate superconductors

The detection of individual quanta of light is important for quantum computation, fluorescence lifetime imaging, single-molecule detection, remote sensing, correlation spectroscopy, and more. Thanks to their broadband operation, high detection efficiency, exceptional signal-to-noise ratio, and fast recovery times, superconducting nanowire single-photon detectors (SNSPDs) have become a critical component in these applications. The operation of SNSPDs based on conventional superconductors, which have a low critical temperature (T c ), requires costly and bulky cryocoolers. This motivated exploration of other superconducting materials with higher T c that would enable single-photon detection at elevated temperatures, yet this task has proven exceedingly difficult. Here we show that with proper processing, high-T c cuprate superconductors can meet this challenge. We fabricated superconducting nanowires (SNWs) out of thin flakes of Bi 2 Sr 2 CaCu2O 8+δ and La 1.55 Sr 0.45 CuO 4 /La 2 CuO 4 (LSCO-LCO) bilayer films and demonstrated their single-photon response up to 25 and 8 K, respectively. The single-photon operation is revealed through the linear scaling of the photon count rate (PCR) on the radiation power. Both of our cuprate-based SNSPDs exhibited single-photon sensitivity at the technologically-important 1.5 μm telecommunications wavelength. Our work expands the family of superconducting materials for SNSPD technology, opens the prospects of raising the temperature ceiling, and raises important questions about the underlying mechanisms of single-photon detection by unconventional superconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Evaluation of gamma-ray transmission through rectangular collimator slits for application in nuclear fuel spectrometry

Gamma-ray spectrometry is widely applied in several science fields, and in particular in non-destructive gamma scanning and gamma emission tomography of irradiated nuclear fuel. Usually, a collimator is used in the experimental setup, to selectively interrogate a region of interest in the fuel. For the optimization of instrument design, as well as for planning measurement campaigns, predictive models for the transmitted gamma-ray intensity through the collimator are needed. Commonly, Monte Carlo Radiation Transport tools are used for accurate prediction of gamma-ray transport, however, the long computation time requirements when used in low-efficiency experimental setups present challenges. Here, the full-energy peak intensity transmitted through a rectangular collimator slit was examined. A uniform planar surface source emitting isotropically was considered, and the rate of photons reaching an ideal counter plane on the opposite side of the collimator was evaluated by analytical integration. To find a closed-form primitive function, some idealizations were required, and thereby parametric models were obtained for the optical field of view, dependent on slit dimensions (length, height and width) and source-to-collimator distance. For contributions from outside the optical field of view, where a closed-form expression cannot be found, fast numerical integral methods were instead used. The results were validated using the Monte Carlo code MCNP6 and show an agreement within three percent for the numerical method. For the analytical method, deviations up to tens of percent were obtained, which is deemed to still be sufficient for instrument design and measurement planning, where often the order of magnitude of the count rate is not a priori known. The method is planned for use in iterative optimization routines in the design of Gamma Emission Tomography devices, as well as for the prediction of gamma spectra obtained in the planning of fuel inspections. An application of the proposed method was demonstrated in spectrum prediction for a short cooling-time fuel rod test from the Halden reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Role of Temperature, Humidity and Rainfall on Influenza Transmission in Guatemala, El Salvador and Panama

Worldwide, seasonal influenza causes about 500,000 deaths and 5 million severe illnesses per year. The environmental drivers of influenza transmission are poorly understood especially in the tropics. We aimed to identify meteorological factors for influenza transmission in tropical Central America. We gathered laboratory-confirmed influenza case-counts by week from Guatemala City, San Salvador Department (El Salvador) and Panama Province from 2006 to 2010. The average total cases per year were: 390 (Guatemala), 99 (San Salvador) and 129 (Panama). Meteorological factors including daily air temperature, rainfall, relative and absolute humidity (RH, AH) were obtained from ground stations, NASA satellites and land models. For these factors, we computed weekly averages and their deviation from the 5-yr means. We assessed the relationship between the number of influenza case-counts and the meteorological factors, including effects lagged by 1 to 4 weeks, using Poisson regression for each site. Our results showed influenza in San Salvador would increase by 1 case within a week of every 1 day with RH>75% (Relative Risk (RR)= 1.32, p=.001) and every 1C increase in minimum temperature (RR=1.29, p=.007) but it would decrease by 1 case for every 1mm-above mean weekly rainfall (RR=0.93,p<.001) (model pseudo-R2=0.55). Within 2 weeks, influenza in Panama was increased by 1 case for every 1% increase in RH (RR=1.04, p=.003), and it was increased by 2 cases for every 1C increase of minimum temperature (RR=2.01, p<.001) (model pseudo-R2=0.4). Influenza counts in Guatemala had 1 case increase for every 1C increase in minimum temperature in the previous week (RR=1.21, p<.001), and for every 1mm/day-above normal increase of rainfall rate (RR=1.03, p=.03) (model pseudo-R2=0.54). Our findings that cases increase with temperature and humidity differ from some temperate-zone studies. But they indicate that climate parameters such as humidity and temperature could be predictive of influenza activity and should be incorporated into country-specific influenza transmission models

Soebiyanto, Radina P.↗

A shared-mobility-based framework for evacuation planning and operations under forecast uncertainty

To meet evacuation needs from carless populations who need personalized assistance to evacuate safely, in this article we propose a ridesharing-based evacuation program that recruits volunteer drivers before a disaster strikes, and then matches volunteer drivers with evacuees once demand is realized. Here we optimize resource planning and evacuation operations under uncertain spatiotemporal demand, and construct a two-stage stochastic mixed-integer program to ensure high demand fulfillment rates. We consider three formulations to improve the number of evacuees served, by minimizing an expected penalty cost, imposing a probabilistic constraint, and enforcing a constraint on the conditional value at risk of the total number of unserved evacuees, respectively. We discuss the benefits and disadvantages of the different risk measures used in the three formulations, given certain carless population sizes and the variety of evacuation modes available. We also develop a heuristic approach to provide quick, dynamic and conservative solutions. We demonstrate the performance of our approaches using five different networks of varying sizes based on regions of Charleston County, South Carolina, an area that experienced a mandatory evacuation order during Hurricane Florence, and utilize real demographic data and hourly traffic count data to estimate the demand distribution.

97 MATHEMATICS AND COMPUTING↗

Toward an AI-Powered Software Pipeline for Real-Time Tracking and Analysis of Wildfire and Smoke

Real-time tracking of wildfires and smoke is crucial for effective response, minimizing damage, protecting lives, and efficiently managing resources during fire emergencies. We develop a web-based AI-powered pipeline that detects wildfires in aerial video and estimates deployment-relevant behavior metrics, including cumulative burned area, burned-area growth rate, fire spread direction, and smoke dispersion. The system combines a YOLO-based detector with YCbCr-based fire segmentation, HSV-based smoke segmentation, Farneback optical flow, and centroid-based spatiotemporal tracking. Using ground sampling distance (GSD), pixel-level fire masks are converted to physical burned-area measurements by correlating fire pixel counts with camera altitude and tilt angle. We benchmark YOLO variants and non-YOLO baselines (GoogLeNet, CNN, DBN, Autoencoder, U-Net, and AlexNet) on the IEEE FLAME dataset and a newly created aerial frame dataset, Wildfire-DB. Cross-dataset evaluation uses a strict threshold-transfer protocol: decision thresholds are selected on FLAME validation and transferred unchanged to Wildfire-DB to quantify generalization under domain shift. YOLOv6 achieves the strongest cross-dataset frame-level fire detection on Wildfire-DB (ROC-AUC 0.8200, PR-AUC 0.8044, and transferred-threshold F1 0.7596). For tracking-oriented deployment requiring oriented localization, YOLO11-OBB provides the most reliable cross-dataset behavior among OBB-capable models while remaining computationally feasible. To analyze the feasibility of UAV deployment, we further measure inference efficiency using synchronized GPU and CPU power logs on a fixed workload of 1569 frames. YOLO-family models process the video in 5.73–12.47 seconds with net energy of 1247.28–1775.39 J, substantially lower latency and energy than heavier classification and reconstruction baselines. Overall, model optimality depends on operational objectives: YOLOv6 is best for cross-dataset detection robustness, whereas YOL...

Color segmentation↗

Hydrogen Plus Other Alternative Fuels Risk Assessment Models (HyRAM+) Version 5.1 Technical Reference Manual

The HyRAM+ software toolkit provides a basis for conducting quantitative risk assessment and consequence modeling for hydrogen, natural gas, and autogas systems. HyRAM+ is designed to facilitate the use of state-of-the-art models to conduct robust, repeatable assessments of safety, hazards, and risk. HyRAM+ integrates deterministic and probabilistic models for quantifying leak sizes and rates, predicting physical effects, characterizing hazards (thermal effects from jet fires, overpressure effects from delayed ignition), and assessing impacts on people. HyRAM+ is developed at Sandia National Laboratories to support the development and revision of national and international codes and standards, and to provide developed models in a publicly-accessible toolkit usable by all stakeholders. This document provides a description of the methodology and models contained in HyRAM+ version 5.1. The most significant changes for HyRAM+ version 5.1 from HyRAM+ version 5.0 are updated default leak frequency values for propane, new default component counts for different fuel types, and an improved fuel specification view in the graphical user interface.

02 PETROLEUM↗

Performance of the OVERFLOW-MLP and LAURA-MLP CFD Codes on the NASA Ames 512 CPU Origin System

The shared memory Multi-Level Parallelism (MLP) technique, developed last year at NASA Ames has been very successful in dramatically improving the performance of important NASA CFD codes. This new and very simple parallel programming technique was first inserted into the OVERFLOW production CFD code in FY 1998. The OVERFLOW-MLP code's parallel performance scaled linearly to 256 CPUs on the NASA Ames 256 CPU Origin 2000 system (steger). Overall performance exceeded 20.1 GFLOP/s, or about 4.5x the performance of a dedicated 16 CPU C90 system. All of this was achieved without any major modification to the original vector based code. The OVERFLOW-MLP code is now in production on the inhouse Origin systems as well as being used offsite at commercial aerospace companies. Partially as a result of this work, NASA Ames has purchased a new 512 CPU Origin 2000 system to further test the limits of parallel performance for NASA codes of interest. This paper presents the performance obtained from the latest optimization efforts on this machine for the LAURA-MLP and OVERFLOW-MLP codes. The Langley Aerothermodynamics Upwind Relaxation Algorithm (LAURA) code is a key simulation tool in the development of the next generation shuttle, interplanetary reentry vehicles, and nearly all "X" plane development. This code sustains about 4-5 GFLOP/s on a dedicated 16 CPU C90. At this rate, expected workloads would require over 100 C90 CPU years of computing over the next few calendar years. It is not feasible to expect that this would be affordable or available to the user community. Dramatic performance gains on cheaper systems are needed. This code is expected to be perhaps the largest consumer of NASA Ames compute cycles per run in the coming year.The OVERFLOW CFD code is extensively used in the government and commercial aerospace communities to evaluate new aircraft designs. It is one of the largest consumers of NASA supercomputing cycles and large simulations of highly resolved full aircraft are routinely undertaken. Typical large problems might require 100s of Cray C90 CPU hours to complete. The dramatic performance gains with the 256 CPU steger system are exciting. Obtaining results in hours instead of months is revolutionizing the way in which aircraft manufacturers are looking at future aircraft simulation work. Figure 2 below is a current state of the art plot of OVERFLOW-MLP performance on the 512 CPU Lomax system. As can be seen, the chart indicates that OVERFLOW-MLP continues to scale linearly with CPU count up to 512 CPUs on a large 35 million point full aircraft RANS simulation. At this point performance is such that a fully converged simulation of 2500 time steps is completed in less than 2 hours of elapsed time. Further work over the next few weeks will improve the performance of this code even further.The LAURA code has been converted to the MLP format as well. This code is currently being optimized for the 512 CPU system. Performance statistics indicate that the goal of 100 GFLOP/s will be achieved by year's end. This amounts to 20x the 16 CPU C90 result and strongly demonstrates the viability of the new parallel systems rapidly solving very large simulations in a production environment.

Taft, James R.↗

Statistical Validation of Multiple Related Data Sets—Case Study Using Interstellar Boundary Explorer Satellite Data

Abstract Space scientists often face the question of whether data collected by different instruments are measurements of the same source population. This paper proposes a statistical validation method for evaluating the agreement between such related data sets. It offers a detailed case study focused on validating a new data set from the Interstellar Boundary Explorer (IBEX) mission, which serves as a practical how-to guide for similar analyses. Since 2008, the IBEX satellite has been gathering data on heliospheric energetic neutral atoms (ENAs) while being exposed to various sources of background noise, such as cosmic rays and solar energetic particles. The IBEX mission initially released only a qualified triple-coincidence (qABC) data product, which was designed to provide observations of ENAs free of background contamination. Further measurements revealed that the qABC data were in fact susceptible to contamination, having relatively low ENA counts and high background rates. To mitigate this issue, the mission team recently considered releasing a certain qualified double-coincidence (qBC) data product, which has roughly twice the detection rate of the qABC data product. This paper presents a simulation-based validation of the new qBC data product against the already-released qABC data product. The results show that the qBCs can plausibly be said to be measuring the same source population as the qABCs up to an average absolute deviation of 3.6%. Visual diagnostics provide additional confirmation of source rate coherence across data products. The framework introduced here is general and can be applied to other validation problems both within and outside the field of space physics.

79 ASTRONOMY AND ASTROPHYSICS↗