Search NASA⌕ Search

SEARCH · Search NASA

Results for “data intensive computing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 775 records · Page 43

Research software engineering: Professionalization, roles, and identity

This report summarizes the results of a study of research software professionals based on 17 interviews conducted at U.S. national laboratories and academic institutions. The study focuses on understanding the emerging professional role of research software engineer (RSE) and examining how the history of software engineering might be relevant to the emergence of research software engineering as a professional movement. The report first uses interview data to identify several models of professional identity for research software professionals. These include dominance models in which a person strongly identifies with a single professional identity (such as RSE), and multiplicity models in which an individual sees an identity like RSE as just one facet of their overall professional role. Next, it identifies three types of work situations for software professionals at national laboratories, covering those who work in dedicated RSE organizations, those who play an integral role in large multiphysics code projects, and those who work in a more isolated capacity to provide software development support for less computationally intensive research projects. It then reviews some key ideas from the sociology of professions in the context of the software engineering profession, suggesting parallels between the current RSE professional movement and aspects of the history of software engineering. It concludes by suggesting some potential practical implications of these findings for the RSE movement and professional organizations.

97 MATHEMATICS AND COMPUTING↗

CFD Benchmark of Pressure Drop in a 61-Pin Wire-Wrapped Assembly with Blocked Channels Using NekRS

Thermal hydraulic behaviors of wire-wrapped rod bundles continue to be a subject of intense research. One of the leading next-generation designs, the sodium fast reactor, relies on a hexagonal assembly of wire-wrapped fuel pins. The issue of potential channel blockage has been raised as a safety concern due to the tightly packed arrangement of the fuel pins. This has led to several recent experimental and computational studies working to quantify the potential impact on the fluid flow and heat transfer behaviors of such blockages. The objective of the present study is to benchmark the high-fidelity NekRS CFD solver in predicting pressure drop for large blockages against available experimental data. A 61-pin wire-wrapped fuel assembly with two flow blockage configurations has been simulated and investigated at various low to moderate Reynolds numbers. The NekRS solver has been shown to yield exponentially decreasing spatial discretization errors with increasing polynomial order. All simulated results agreed well with measured data, which indicates that the overall methodology is adequate and consistent. The results of this benchmark study demonstrate the accuracy of NekRS for sodium fast reactor hydrodynamic simulations, increasing the confidence in its use for design, licensing, and analysis activities.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Challenges and opportunities in understanding microbial communities with metagenome assembly (accompanied by IPython Notebook tutorial)

Metagenomic investigations hold great promise for informing the genetics, physiology, and ecology of environmental microorganisms. Current challenges for metagenomic analysis are related to our ability to connect the dots between sequencing reads, their population of origin, and their encoding functions. Assembly-based methods reduce dataset size by extending overlapping reads into larger contiguous sequences (contigs), providing contextual information for genetic sequences that does not rely on existing references. These methods, however, tend to be computationally intensive and are again challenged by sequencing errors as well as by genomic repeats. While numerous tools have been developed based on these methodological concepts, they present confounding choices and training requirements to metagenomic investigators. To help with accessibility to assembly tools, this review also includes an IPython Notebook metagenomic assembly tutorial. This tutorial has instructions for execution any operating system using Amazon Elastic Cloud Compute and guides users through downloading, assembly, and mapping reads to contigs of a mock microbiome metagenome. Despite its challenges, metagenomic analysis has already revealed novel insights into many environments on Earth. As software, training, and data continue to emerge, metagenomic data access and its discoveries will to grow.

54 ENVIRONMENTAL SCIENCES↗

Nucleosynthesis for Science and Security: Preparing for a Nuclear and Chemical Sciences Program at FRIB

Nuclear reactions are vital to nuclear science and its many applications: they are the mechanism for forging elements in the cosmos, the source of energy for both stars and nuclear weapons, the principal probes used to study nuclear properties, and the primary diagnostic signatures for stockpile stewardship and nuclear forensics. Yet, most reaction experiments have been confined to the small fraction of all nuclei that are stable, and a theoretical foundation for reactions with a true predictive capability is lacking. The Facility for Rare Isotope Beams (FRIB), the flagship U.S. facility in low-energy nuclear science, will offer the opportunity to rectify the long-standing gap in understanding reactions on radioactive nuclei. The work performed in this LDRD project brought together the experimental nuclear physics, radiochemistry, nuclear theory, and high-performance computing expertise of the Nuclear and Chemical Sciences Division to develop new techniques that will be needed to exploit FRIB beams to better understand nuclear reactions. Novel isotope-harvesting techniques were developed to demonstrate that samples could be collected at the National Superconducting Cyclotron Laboratory (NSCL) at Michigan State University (the facility which FRIB builds upon) for direct measurements of long-lived radioisotopes. For isotopes too short lived for direct measurements, indirect approaches were developed that would ultimately be used with high-intensity FRIB radioactive beams. Finally, to better interpret experimental results and reliably extend predictions beyond the available data, a foundation was established for a fully microscopic and unified description of the structure, interactions, and dynamics of two reacting nuclei that moves beyond empirical models. This work has positioned the team for a leadership role at FRIB, enabling LLNL to maintain a world-class nuclear science program and to compete for future talent and funding, both of which are vital for continuing support of the laboratory’s stockpile stewardship mission.

42 ENGINEERING↗

Desmearing Bonse–Hart USANS data using Bayesian Gaussian process regression

Ultra-small-angle neutron scattering (USANS) enables access to micrometer-scale structures but is intrinsically affected by strong, anisotropic resolution smearing arising from slit-geometry optics. As a result, recovery of the intrinsic scattering intensity constitutes an ill-posed inverse problem, and commonly used iterative desmearing methods lack rigorous uncertainty quantification. We present a Bayesian desmearing framework for slit-geometry USANS based on Gaussian process regression. In this approach, the scattering intensity is modeled as a smooth random function, and the instrumental point spread function is incorporated explicitly as a forward operator. The resulting formulation yields a closed-form maximum a posteriori solution with well-defined credibility intervals. Computational benchmarks and experimental validation using combined USANS and small-angle neutron scattering (SANS) measurements demonstrate that the framework enables stable desmearing, suppresses experimental noise, and preserves physically meaningful structural features under realistic conditions.

Tung, Chi-Huan [Oak Ridge National Laboratory (ORN↗

Mechanism-based constitutive modeling of ZEK100 magnesium alloy with crystal plasticity and in-situ HEXRD experiment

The constitutive behavior of a hexagonal close-packed (HCP) polycrystalline ZEK100 magnesium alloy was investigated here using combined high energy X-ray diffraction (HEXRD) from a synchrotron source and crystal plasticity modeling approach. The in-situ tensile test data coupled with the HEXRD enabled the tracking of the lattice strain evolution during deformation. The microscopic behavior represented by lattice strain and the macroscopic behavior represented by stress-strain curves were then used together as objective function to estimate the critical resolved shear stress (CRSS) and hardening parameters of available slip and deformation twin systems in the ZEK100 alloy. An enhanced predominant twinning reorientation (ePTR) scheme was proposed in the current work, and the ePTR parameters were determined for the first time by the use of basal plane peak intensity along loading direction measured from HEXRD. Two crystal plasticity models, the computationally efficient elastic-plastic self-consistent (EPSC) and crystal plasticity finite element (CPFE) models, were developed incorporating the deformation twinning for the HCP-structured metals. The determined constitutive parameters were further validated by comparing the predicted deformation texture with the measured one. The work provides a useful and computationally-efficient modeling scheme to understand the slip/twin induced deformation behaviors of the ZEK100 alloy in micro- and macro-scales.

36 MATERIALS SCIENCE↗

BEYONDPLANCK XIII. Intensity foreground sampling, degeneracies, and priors

We present the intensity foreground algorithms and model employed within the BEYONDPLANCK analysis framework. The BEYONDPLANCK analysis is aimed at integrating component separation and instrumental parameter sampling within a global framework, leading to complete end-to-end error propagation in the Planck Low Frequency Instrument (LFI) data analysis. Given the scope of the BEYONDPLANCK analysis, a limited set of data is included in the component separation process, leading to foreground parameter degeneracies. In order to properly constrain the Galactic foreground parameters, we improve upon the previous Commander component separation implementation by adding a suite of algorithmic techniques. These algorithms are designed to improve the stability and computational efficiency for weakly constrained posterior distributions. These are: (1) joint foreground spectral parameter and amplitude sampling, building on ideas from MIRAMARE; (2) component-based monopole determination; (3) joint spectral parameter and monopole sampling; and (4) application of informative spatial priors for component amplitude maps. We find that the only spectral parameter with a significant signal-to-noise ratio using the current BEYONDPLANCK data set is the peak frequency of the anomalous microwave emission component, for which we find $ν$ p = 25.3 ± 0.5 GHz; all others must be constrained through external priors. Future works will be aimed at integrating many more data sets into this analysis, both map and time-ordered based, thereby gradually eliminating the currently observed degeneracies in a controlled manner with respect to both instrumental systematic effects and astrophysical degeneracies. When this happens, the simple LFI-oriented data model employed in the current work will need to be generalized to account for both a richer astrophysical model and additional instrumental effects. This work will be organized within the Open Science-based COSMOGLOBE community effort.

79 ASTRONOMY AND ASTROPHYSICS↗

Rapid 3D nanoscale coherent imaging via physics-aware deep learning

Phase retrieval, the problem of recovering lost phase information from measured intensity alone, is an inverse problem that is widely faced in various imaging modalities ranging from astronomy to nanoscale imaging. The current process of phase recovery is iterative in nature. As a result, the image formation is time consuming and computationally expensive, precluding real-time imaging. Here, we use 3D nanoscale X-ray imaging as a representative example to develop a deep learning model to address this phase retrieval problem. We introduce 3D-CDI-NN, a deep convolutional neural network and differential programing framework trained to predict 3D structure and strain, solely from input 3D X-ray coherent scattering data. Our networks are designed to be “physics-aware” in multiple aspects; in that the physics of the X-ray scattering process is explicitly enforced in the training of the network, and the training data are drawn from atomistic simulations that are representative of the physics of the material. We further refine the neural network prediction through a physics-based optimization procedure to enable maximum accuracy at lowest computational cost. 3D-CDI-NN can invert a 3D coherent diffraction pattern to real-space structure and strain hundreds of times faster than traditional iterative phase retrieval methods. Our integrated machine learning and differential programing solution to the phase retrieval problem is broadly applicable across inverse problems in other application areas.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Effect of alloying on intrinsic ductility in WTaCrV high entropy alloys

Tungsten (W) exhibits desirable properties for extreme applications, such as the divertor in magnetic fusion reactors, but its practicality remains limited due to poor formability and insufficient irradiation resistance. In this work, we study the intrinsic ductility of body-centered cubic WTaCrV based high entropy alloys (HEAs), which are known to exhibit excellent irradiation resistance. The ductility evaluations are carried out using a criterion based on the competition between the critical stress intensity factors for emission (K Ie ) and cleavage (K Ic ) in the {110} slip planes and {110} crack planes, which are evaluated within the linear elastic fracture mechanics framework and computed using density functional theory. The results suggest that increasing the alloying concentrations of V and reducing the concentrations of W can significantly improve the ductility in these HEAs. The elastic anisotropy for these HEAs is analyzed using the Zener anisotropy ratio and its correlation with the concentration of W in the alloys is studied. Results indicate that these alloys tend to be fairly isotropic independently from the concentration of W in them. The computed data for the elastic constants of these HEAs is also compared against available experimental data. The results are in good agreement, validating the robustness and accuracy of the computational methods. Multiple phenomenological ductility metrics were also computed and analyzed against the analytical model. Some metrics, mainly the surrogate D parameter, show a good correlation with the Rice model. The potential of these empirical metrics to serve as surrogate screening models for optimizing the compositional space is also discussed.

Anisotropy↗

A Modified Model Dielectric Function for Analyzing Optical Spectra of InGaN Nanofilms on Sapphire Substrates

Due to a lower InN bandgap energy Eg~0.7 eV, InxGa1−xN/Sapphire epifilms are considered valuable in the development of low-dimensional heterostructure-based photonic devices. Adjusting the composition x and thickness d in epitaxially grown films has offered many possibilities of light emission across a wide spectral range, from ultraviolet through visible into near-infrared regions. Optical properties have played important roles in making semiconductor materials useful in electro-optic applications. Despite the efforts to grow InxGa1−xN/Sapphire samples, no x- and d-dependent optical studies exist for ultrathin films. Many researchers have used computationally intensive methods to study the electronic band structures Ejk→, and subsequently derive optical properties. By including inter-band transitions at critical points from Ejk→, we have developed a semiempirical approach to comprehend the optical characteristics of InN, GaN and InxGa1−xN. Refractive indices of InxGa1−xN and sapphire substrate are meticulously integrated into a transfer matrix method to simulate d- and x-dependent reflectivity RE and transmission TE spectra of nanostructured InxGa1−xN/Sapphire epifilms. Analyses of RE and TE have offered accurate x-dependent shifts of energy gaps for InxGa1−xN (x = 0.5, 0.7) in excellent agreement with the experimental data.

Talwar, Devki N. (ORCID:0000000202723195)↗

Small-angle X-ray scattering: a bridge between RNA secondary structures and three-dimensional topological structures

Whereas the structures of small to medium-sized well folded RNA molecules often can be determined by either X-ray crystallography or NMR spectroscopy, obtaining structural information for large RNAs using experimental, computational, or combined approaches remains a major interest and challenge. RNA is very sensitive to small-angle X-ray scattering (SAXS) due to high electron density along phosphate-sugar backbones, whose scattering contribution dominates SAXS intensity. For this reason, SAXS is particularly useful in obtaining global RNA structural information that outlines backbone topologies and, therefore, molecular envelopes. Such information is extremely valuable in bridging the gap between the secondary structures and three-dimensional topological structures of RNAmolecules, particularly those that have proven difficult to study using other structuredetermination methods. Here we review published results of RNA topological structures derived from SAXS data or in combination with other experimental data, as well as details on RNA sample preparation for SAXS experiments.

59 BASIC BIOLOGICAL SCIENCES↗

Angular power spectrum of the diffuse gamma-ray emission as measured by the Fermi Large Area Telescope and constraints on its dark matter interpretation

The isotropic gamma-ray background arises from the contribution of unresolved sources, including members of confirmed source classes and proposed gamma-ray emitters such as the radiation induced by dark matter annihilation and decay. Clues about the properties of the contributing sources are imprinted in the anisotropy characteristics of the gamma-ray background. We use 81 months of Pass 7 Reprocessed data from the Fermi Large Area Telescope to perform a measurement of the anisotropy angular power spectrum of the gamma-ray background. Here, we analyze energies between 0.5 and 500 GeV, extending the range considered in the previous measurement based on 22 months of data. We also compute, for the first time, the cross-correlation angular power spectrum between different energy bins. The derived angular spectra are compatible with being Poissonian, i.e. constant in multipole. Furthermore, the energy dependence of the anisotropy suggests that the signal is due to two populations of sources, contributing, respectively, below and above ~ 2 GeV . Finally, using data from state-of-the-art numerical simulations to model the dark matter distribution, we constrain the contribution from dark matter annihilation and decay in Galactic and extra-Galactic structures to the measured anisotropy. These constraints are competitive with those that can be derived from the average intensity of the isotropic gamma-ray background.

79 ASTRONOMY AND ASTROPHYSICS↗

Experimental study of multiple-orientation muon tomography with image optimization in sparse data environments

Due to the high penetrating power of cosmic-ray muons, they can be used to probe very thick and dense objects. As muons are charged particles, they can be tracked by ionization detectors, determining the position and direction of the muons. With detectors on either side of an object to measure particle direction change, scattering information within the object can be found. This can be used to produce a scattering-intensity image within the object related to density and atomic number. Such imaging is typically performed with a single detector-object orientation, taking advantage of the more intense downward flux of muons, producing planar imaging with some depth-of-field information in the third dimension. Several simulation studies were published with multiorientation tomography, which can form a three-dimensional representation faster than a single-orientation view. In this study, experimental muon-scatter-based tomography was performed using a concrete filled steel drum with several different metal wedges inside, with the drum between detector planes. Data were collected from different detector-object orientations by rotating the steel drum. The data collected from each orientation were combined using two different tomographic methods. A traditional inverse Radon transform approach used for computed tomography and a combination of multiple depth-of-field reconstructions were applied to the data. As cosmic-ray muon flux imaging is rate limited, the imaging techniques were compared for sparse data. Using the combined depth-of-field reconstruction technique, fewer detector-object orientations were needed to reconstruct images that could be used to differentiate the metal wedges.

47 OTHER INSTRUMENTATION↗

Feedstock Carbon Intensity Calculator (FD-CIC) Users’ Manual and Technical Documentation

A transparent and easy-to-use tool for feedstock-specific, farm-level CI calculation of feedstocks is helpful. With the ARPA-E support, we have developed a tool–the Feedstock Carbon Intensity Calculator (FD-CIC). The first version of FD-CIC was released with the GREET® model in 2020 (Wang et al., 2020) so that corn feedstock producers could use this publicly available tool (https://greet.es.anl.gov/tool_fd_cic) to quantify corn grain CIs with farm-level input data and management practices. In the 2021 version, we expand the tool’s capabilities by including additional feedstocks such as soybeans, sorghum, and rice. Similar to corn, it calculates the farm-level CI for these feedstocks by allowing user-defined farm-level farming inputs and incorporating the GHG emission intensities of these inputs from GREET (in particular, GREET1, the fuel cycle model of GREET).

09 BIOMASS FUELS↗

Skewering the silos: using Brick to enable portable analytics, modeling and controls in buildings

Nearly all large commercial buildings have heating, ventilation and air conditioning (HVAC) systems, lighting systems, safety and other systems controlled by a computer—a dedicated server with a building energy management system (BMS). However, these BMSs are proprietary with each building’s assets (that is, fans, valves, pumps, and their setpoints) named and coded uniquely by the BMS vendor or engineer; building analytics and control algorithms are written specific to the assets and the building. Thus, any control updates or analytics to improve building performance—especially critical to reduce greenhouse emissions or improve load flexibility—are labor intensive and costly. The Brick schema was developed so the same analysis or control algorithms can work on a variety of buildings if each is digitally represented in a Brick data model. The goal of this project was to further the development of Brick to extend it beyond an academic project with demonstrated success in a small field study, to a practical choice for industrial and commercial stakeholders seeking to realize value from building data. To do this, we executed four objectives: (1) expand the Brick schema including its modeling capabilities and vocabulary, (2) develop tools for integrating Brick with existing digital technologies and representations in buildings, (3) develop an open-source analytics platform to facilitate use of Brick in delivering data value, and (4) demonstrate Brick-driven analytics and controls in real settings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Suite of Comprehensive Rapid Analysis Tools for Estimation (SOCRATES) - 19076

Long-term monitoring and remedy implementation require ongoing assessment and analytics, much of which relies on intensive and time-consuming data integration and engineering, data manipulation, and 'big-data' analytics. Many of these analyses, such as groundwater flow direction, contaminant rate attenuation, and non-parametric statistics, are spatial-temporal in nature and involve coalescing disparate data sets from a variety of sources with measurement frequencies from milliseconds to years, as well as total data volumes often exceeding the petabyte scale. Web-based analytical tools are advantageous for ongoing real-time or near-real-time ingestion of data from various sources for general data access and visualization, and also allow for complex, data rich calculations by leveraging expansive networks and virtually unlimited compute power. A web-based framework offers a natural architecture for such endeavors, as access to compute nodes through auto-scaling and dynamic frameworks enables speedy data manipulation and unlimited access to any data set that is made available via an http request. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Software and Hardware Development for the Mu2e Trigger and Data Acquisition System

The Mu2e experiment at Fermilab will be 10,000 times more sensitive than previous experiments looking for muon-to-electron conversion. This precise and complex experimental apparatus will produce 200 million billion muons per year. The accelerator complex repurposes elements of the infrastructure that produced high energy proton and anti-protons beams for the Tevatron experiments to produce the high-intensity muon beams necessary for Mu2e and the Muon (g-2) experiments.

43 PARTICLE ACCELERATORS↗

The cross-correlation between 21 cm intensity mapping maps and the Lyα forest in the post-reionization era

We investigate the cross-correlation signal between 21cm intensity mapping maps and the Lyα forest in the fully non-linear regime using state-of-the-art hydrodynamic simulations. The cross-correlation signal between the Lyα forest and 21cm maps can provide a coherent and comprehensive picture of the neutral hydrogen (HI) content of our Universe in the post-reionization era, probing both its mass content and volume distribution. We compute the auto-power spectra of both fields together with their cross-power spectrum at z = 2.4 and find that on large scales the fields are completely anti-correlated. This anti-correlation arises because regions with high (low) 21cm emission, such as those with a large (low) concentration of damped Lyα systems, will show up as regions with low (high) transmitted flux. We find that on scales smaller than k ≅ 0.2 h Mpc{sup −1} the cross-correlation coefficient departs from −1, at a scale where non-linearities show up. We use the anisotropy of the power spectra in redshift-space to determine the values of the bias and of the redshift-space distortion parameters of both fields. We find that the errors on the value of the cosmological and astrophysical parameters could decrease by 30% when adding data from the cross-power spectrum, in a conservative analysis. Our results point out that linear theory is capable of reproducing the shape and amplitude of the cross-power up to rather non-linear scales. Finally, we find that the 21cm-Lyα cross-power spectrum can be detected by combining data from a BOSS-like survey together with 21cm intensity mapping observations by SKA1-MID with a S/N ratio higher than 3 in k element of [0.06,1] h Mpc{sup −1}. We emphasize that while the shape and amplitude of the 21cm auto-power spectrum can be severely affected by residual foreground contamination, cross-power spectra will be less sensitive to that and therefore can be used to identify systematics in the 21cm maps.

79 ASTRONOMY AND ASTROPHYSICS↗