SEARCH · Search NASA
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.
This dataset contains the Recorder trace files of 91 built-in tests from three widely-used I/O libraries. It was utilized in our IPDPS'25 paper, "VerifyIO: Verifying Adherence to Parallel I/O Consistency Semantics". The dataset can be used to reproduce the results presented in the paper and to conduct additional analyses
Cellular components are non-randomly arranged with respect to the shape and polarity of the whole cell. Patterning within cells can extend down to the level of individual proteins and mRNA. But how much of the proteome is actually localized with respect to cell polarity axes? Proteomics combined with cellular fractionation has shown that most proteins localize to one or more organelles but does not tell us how many proteins have a polarized localization with respect to the large-scale polarity axes of the intact cell. Genome-wide localization studies in yeast found that only a few percent of proteins have a localized position relative to the cell polarity axis defined by sites of polarized cell growth. Here, we describe an approach for analyzing protein distribution within a cell with a visibly obvious global patterning - the giant ciliate Stentor coeruleus. Ciliates, including Stentor, have highly polarized cell shapes with visible surface patterning. A Stentor cell is roughly 2mmlong, allowing a ‘‘proteomic dissection’’ in which microsurgery is used to separate cellular fragments along the anterior-posterior axis, followed by comparative proteomic analysis. In our analysis, 25% of the proteome, including signaling proteins, centrin/SFI proteins, and GAS2 orthologs, shows a polarized location along the cell’s anterior-posterior axis. We conclude that a large proportion of all proteins are polarized with respect to global cell polarity axes and that proteomic dissection provides a simple and effective approach for spatial proteomics.
There exist many methods for sensitivity analysis readily available to the practitioner. While each seeks to help the modeler answer the same general question – How do sources of uncertainty or changes in the model inputs relate to uncertainty in the output? – different methods are associated with different assumptions, constraints, and required resources, leading to conclusions that may vary in interpretability and level of detail. Thus, it is crucial that the practitioner selects the desired sensitivity analysis method judiciously, making sure to match the selected approach to the specifics of their problem and to their desired objectives. In this chapter, we provide a practical overview of a collection of widely used, widely available sensitivity analysis methods. We focus on global sensitivity approaches, which seek to characterize how uncertainty in the model output may be allocated to sources of uncertainty in model inputs across the entire input space. Generally, this will require the practitioner to specify a probability distribution over the input space. On the other hand, methods for local sensitivity analysis do not require this specification but they have more limited utility, providing insight into sources of uncertainty associated only with a particular, specified location in the input space. Our hope is that this chapter may serve as a decision-making tool for practitioners, helping to guide the selection of a sensitivity analysis approach that will best fit their needs. To support this goal, we have selected a suite of approaches to cover, which, while not exhaustive, we believe provides a flexible and robust sensitivity analysis toolkit. All methods included are widely used and available in standard software packages.
Previous studies have shown that the two-layer model more accurately predicts hydrogen dispersion than the conventional notional nozzle models without significantly increasing the computational expense. However, the model was only validated for predicting the concentration distribution and has not been adequately validated for predicting the velocity distributions. In the present study, particle imaging velocimetry (PIV) was used to measure the velocity field of an underexpanded hydrogen jet released at 10 bar from a 1.5 mm diameter orifice. The two-layer model was the used to calculate the inlet conditions for a two-dimensional axisymmetric CFD model to simulate the hydrogen jet downstream of the Mach disk. The predicted velocity spreading and centerline decay rates agreed well with the PIV measurements. The predicted concentration distribution was consistent with data from previous planar Rayleigh scattering measurements used to verify the concentration distribution predictions in an earlier study. The jet spreading was also simulated using several widely used notional nozzle models combined with the integral plume model for comparison. These results show that the velocity and concentration distributions are both better predicted by the two-layer model than the notional nozzle models to complement previous studies verifying only the predicted concentration profiles. Thus, this study shows that the two-layer model can accurately predict the jet velocity distributions as well as the concentration distributions as verified earlier. Though more validation studies are needed to improve confidence in the model and increase the range of validity, the present work indicates that the two-layer model is a promising tool for fast, accurate predictions of the flow fields of underexpanded hydrogen jets.
We discuss how Information and communication technologies (ICT) have emerged as one of the leading technologies to reduce global emissions, particularly in the mobility sector. Automotive electronics, such as sensors/actuators and microcontrollers (commonly known as electronic control units (ECU) which control one or more of the electrical systems or subsystems in a vehicle), play a key role in ICT. Sensors/actuators are key in electronic ICT devices, starting with the data collection and data communication with the internet. The latest two big trends of electrification and automation in vehicles, are projected to increase the use of worldwide automotive sensors from 7.5 billion units in 2017 to 11.0 billion units by the year 2024. A representative state-of-the art automotive sensor system, (i.e., an ultrasonic backup system), including the ECUs controlling the vehicle electrical systems/subsystems, has been considered to estimate the energy footprint in terms of manufacturing and operational energy of global automotive sensors use. A widely used life cycle energy assessment method (i.e., cumulative energy demand) was used as both direct and indirect (including the extraction, manufacturing, and disposal of the raw and auxiliary materials) energy use can be considered for the energy footprint estimation. The embodied manufacturing energy impacts of the system was estimated to be 559 MJ/system, compared to the 417 MJ/system for lifetime system power and additional gasoline use. The share of purchased energy to the embodied energy where the upstream energy isn't included in the former case, is less than 10% and ~ 85% for the component manufacturing and vehicle operation energy uses, respectively. As the purpose of this ultrasonic backup system is to prevent rear crashes, an estimated 1.0 MJ/system is avoided from reduced lifetime vehicular repairs (from an estimated 11% chance of requiring a rear bumper replacement). While all of this is small compared to the overall automotive manufacturing and use energy, the 11 billion automotive sensors expected to be produced in 2024 could require 1540 PJ for manufacturing and those sensors would require an additional 780–1150 PJ for lifetime energy use.
Graphical networks are useful, widely-used modeling approaches to represent complex biological processes with biological measurements generated by platforms such as mass spectrometry. Bayesian analyses of graphical networks for omics data have several advantages over their frequentist counterparts, such as the inclusion of prior knowledge in the estimation of models. However, Bayesian approaches to date have only been feasible for data with a couple hundred biomolecules due to prohibitive computational time, but omics data often contains tens of thousands of biomolecules. Here, we present and illustrate a more computationally efficient approach named BPlane (Bayesian PseudoLikelihood-based Algorithm for Network Estimation) to extend Bayesian modeling capabilities for larger-sized datasets, such as most untargeted proteomics data. Via simulation, we demonstrate that BPlane produces substantial computational savings over a current state-of-the-art Bayesian algorithm while maintaining competitive edge detection accuracy. On a SARS-CoV2 proteomics data with 7000 proteins, the competing algorithm takes three times as long to complete the first iteration as BPlane takes to converge after over 100 iterations.
Understanding the role of ferroelectric polarization in modulating the electronic and structural properties of crystals is critical for advancing these materials for overcoming various technological and scientific challenges. However, due to difficulties in performing experimental methods with the required resolution, or in interpreting the results of methods therein, the nanoscale morphology and response of these surfaces to external electric fields has not been properly elaborated. Here, in this work, we investigate the effect of ferroelectric polarization and local distortions in a BaTiO 3 perovskite, using two widely used computational approaches which treat the many-body nature of X-ray excitations using different philosophies, namely the many-body, delta-self-consistent-field determinant (mb-ΔSCF) and the Bethe–Salpeter equation (BSE) approaches. We show that in agreement with our experiments, both approaches consistently predict higher excitations of the main peak in the O–K edge for the surface with upward polarization. However, the mb-ΔSCF approach mostly fails to capture the L 2,3 separations at the Ti–L edge, due to the absence of spin–orbit coupling in Kohn–Sham density functional theory (KS-DFT) at the generalized gradient approximation level. On the other hand, and most promising, we show that application of the GW/BSE approach successfully reproduces the experimental XAS, both the relative peak intensities as well as the L 2,3 separations at the Ti–L edges upon ferroelectric switching. Thus simulated XAS is shown to be a powerful method for capturing the nanoscale structure of complex materials, and we underscore the need for many-body perturbation approaches, with explicit consideration of core-hole and multiplet effects, for capturing the essential physics in these systems.
Context. Satellite galaxies experience multiple physical processes when interacting with their host halos, often leading to the quenching of star formation. In the Local Group, satellite quenching has been shown to be highly efficient, affecting nearly all satellites except the most massive ones. While recent surveys study Milky Way-analogs to assess how representative our Local Group is, the dominant physical mechanisms behind satellite quenching in Milky Way-mass halos remain under debate. Aims. We analyze satellite quenching within the same Milky Way-mass halo simulated using various widely used astrophysical codes, each using different hydrodynamic methods and implementing different supernovae feedback recipes. The goal is to determine whether quenched fractions, quenching timescales, and the dominant quenching mechanisms are consistent across codes or if they show sensitivity to the specific hydrodynamic method and supernovae feedback physics employed. Methods. We used a subset of high-resolution cosmological zoom-in simulations of a Milky Way-mass halo from the multiple-code AGORA CosmoRun suite. Our analysis focuses on comparing satellite quenching across the different models and against observational data. We also analyzed the dominant mechanisms driving satellite quenching in each model. Results. We find that the quenched fraction is consistent with the latest SAGA Survey results within its 1σ host-to-host scatter across all the models. Regarding quenching timescales, all the models reproduce the trend observed in the ELVES survey, Local Group observations, and previous simulations: The less massive the satellite, the shorter its quenching timescale. All of our models converge on the dominant quenching mechanisms: Strangulation halts cold gas accretion in all satellites, while ram pressure stripping is the predominant mechanism for gas removal, and it is particularly effective in satellites with M * <10 8 M ⊙ . Nevertheless, the efficiency of the stripping mechanisms differs among the codes, showing a strong sensitivity to the different supernovae feedback implementations and/or hydrodynamic methods employed.
The spin cutoff parameter for 59 Ni has been studied from different types of experimental data, including neutron angular distributions from the 56 Fe(α,n) 59 Ni reaction, spin of discrete levels, the level density from the proton evaporation spectrum of the 54 Fe( 6 Li,p) 59 Ni reaction and neutron resonance spacing. Experimental data points were compared with calculations using models widely used in literature. It was found that the available empirical models overestimate data points in the energy region below the neutron separation energy, while microscopic calculations which take into account pairing correlations within Hartree-Fock +Bardeen-Cooper-Schrieffer approach are consistent with data. As a result, it confirmed earlier findings that pairing correlations play an important role and need to be taken into account when the spin cutoff parameter is calculated below the neutron separation energy.
A novel high-temperature particle solar receiver is developed using a light trapping planar cavity configuration. As particles fall through the cavity, the concentrated solar radiation warms the boundaries of the receiver and in turn heats the particles. Particles flow through the system, forming a fluidized bed at the lower section, leaving the system from the bottom at a constant flowrate. Air is introduced to the system as the fluidizing medium to improve particle heat transfer and mixing. A laboratory scale cavity receiver is built by collaborators at the Colorado School of Mines and their data are used for model validation. In this experimental setup, near IR quartz lamp is used to provide flux to the vertical wall of the heat exchanger. The system is modeled using the discrete element method and a continuum two-fluid method. The computational model matches the experimental system size and the particle size distribution is assumed monodisperse. A new continuum conduction model that accounts for the effects of solid concentration is implemented, and the heat flux boundary condition matches the experimental setup. Radiative heat transfer is estimated using a widely used correlation during the post-processing step to determine an overall heat transfer coefficient. The model is validated against testing data and achieves less than 30% discrepancy and a heat transfer coefficient greater than 1000 W/m2 K.
A novel high temperature particle solar receiver is developed by using a light trapping planar cavity configuration. As particles fall through the cavity, the concentrated solar radiation warms the boundaries of the receiver and in turn heats the particles. Particles flow through the system, forming a packed bed at the lower end, leaving the system from the bottom at a constant flow rate. Air is introduced to the system as the fluidizing medium to improve particle heat transfer and mixing. A laboratory scale cavity receiver is built and a near IR quartz lamp is used to provide flux to the vertical wall of the heat exchanger. The system is modeled using a continuum two-fluid method. The computational model matches the experimental system size and the particle size distribution is assumed monodisperse. A conduction model that accounts for the effects of solid concentration is implemented, and the heat flux boundary condition matches the experimental setup. Radiative heat transfer is estimated using a widely used correlation during the post-processing step to determine an overall heat transfer coefficient. The model is validated against testing data and achieves less than 30% discrepancy and a heat transfer coefficient greater than 1000 W/m2K.
Manually tuning physics-based animation parameters to explore a simulation outcome space or achieve desired motion outcomes can be notoriously tedious. This problem has motivated many sophisticated and specialized optimization-based methods for fine-grained (keyframe) control, each of which are typically limited to specific animation phenomena, usually complicated, and, unfortunately, not widely used. In this paper, we propose Unified Many-Worlds Browsing (UMWB), a practical method for sample-level control and exploration of physics-based animations. Our approach supports browsing of large simulation ensembles of arbitrary animation phenomena by using a unified volumetric WORLDPACK representation based on spatiotemporally compressed voxel data associated with geometric occupancy and other low-fidelity animation state. Beyond memory reduction, the WORLDPACK representation also enables unified query support for interactive browsing: it provides fast evaluation of approximate spatiotemporal queries, such as occupancy tests that find ensemble samples ("worlds") where material is either IN or NOT IN a user-specified spacetime region. WORLDPACKS also support real-time hardware-accelerated voxel rendering by exploiting the spatially hierarchical and temporal RLE raster data structure. Our UMWB implementation supports interactive browsing (and offline refinement) of ensembles containing thousands of simulation samples, and fast spatiotemporal queries and ranking. We show UMWB results using a wide variety of physics-based animation phenomena---not just JELL-O ® .
Woody biomass will be an essential feedstock for a large-scale cellulosic biofuel industry. The life cycle carbon accounting for the production of biofuels from woody feedstocks is complex and has engendered significant controversy. Issues such as below ground carbon, carbon debt, varied regional forest practices, multiple parallel forest product lines, and development of realistic counterfactual scenarios all contribute to the complexity. DOE funded CORRIM to develop comprehensive and definitive lifecycle inventories and assessments on the production of fuels from woody feedstocks. CORRIM brought together expertise in forest practices, short rotation woody crops, bioconversion of woody biomass, process modeling, and life cycle assessment to successfully accomplish this goal. The CORRIM team has developed data on forest productivity, fuel usage, and fuel production for six regionally specific forest systems. These six forest systems include three current commercial systems; southern pine plantations, Douglas-fir plantations, naturally regenerated Northeastern (NE) spruce/fir, and three ‘short rotation woody crops’, poplar, eucalyptus and willow, which are at different stages of demonstration in the US. The fuel production systems include cellulosic ethanol and bio-oil based hydrocarbons. The project was successfully reviewed at the BETO Program review in March 2015, 2017, and 2019 and has resulted currently in 19 publications and reports and 35 presentations at both national and international conferences, with more in the pipeline. See publications and presentations for links to each document. Finally, and most importantly, CORRIM has gone beyond the scope of the original DOE proposal to work closely with GREET at Argonne National Laboratory (ANL) to incorporate all the life cycle data and scenario models into their modeling system. GREET is the most widely used and definitive information source for evaluating lifecycle carbon emissions for fuels. The incorporation of CORRIM data from this project guarantees the results of the research will be extensively used and widely disseminated. Technical process improvements and policy relevant accomplishments are detailed in the relevant programmatic sections in the full report including citations therein. Highlights are summarized here for easy reference.
This study was motivated by earlier observations. It is a systematic examination of the adequacy of reporting of information (metadata) necessary to understand x-ray photoelectron spectroscopy (XPS) data collection and data analysis in the scientific literature. The information for this study was obtained from papers published in three high-quality journals over a six-month period in 2019 and throughout 2021. Each paper was evaluated to determine whether the authors had reported (percentages of the papers properly providing the information are given in parentheses) the spectrometer (66%), fitting software (15%), x-ray source (40%), pass energy (10%), spot size (5%), synthetic peak shapes in fits (10%), backgrounds in fits (10%), whether the XPS data are shown in the main body of the paper or in the supporting information (or both), and whether fitted or unfitted spectra were shown (80% of published spectra are fit). The Shirley background is the most widely used background in XPS peak fitting. The Al Kα source is the most widely used x-ray source for XPS data collection. CASAXPS is the most widely used fitting program for XPS data analysis. Further, there is good agreement between the results gathered during the two years of our survey. There are some hints the situation may be improving. This study also provides a list of the information/parameters that should be reported when XPS is performed.
Multibeam sonars are widely used for environmental monitoring of fauna at marine renewable energy sites. However, they can rapidly accrue vast volumes of data, which poses a challenge for data processing. Here, using data from a deployment in a tidal channel with peak currents of 1–2 m s –1 , we demonstrate the data-reduction benefits of real-time automatic classification of targets detected and tracked in multibeam sonar data. First, we evaluate classification capabilities for three machine learning algorithms: random forests, support vector machines, and k-nearest neighbors. For each algorithm, a hill-climbing search optimizes a set of hand-engineered attributes that describe tracked targets. Here, the random forest algorithm is found to be most effective—in postprocessing, discriminating between biological and nonbiological targets with a recall rate of 0.97 and a precision of 0.60. In addition, 89% of biological targets are correctly classified as either seals, diving birds, fish schools, or small targets. Model dependence on the volume of training data is evaluated. Second, a real-time implementation of the model is shown to distinguish between biological targets and nonbiological targets with nearly the same performance as in postprocessing. From this, we make general recommendations for implementing real-time classification of biological targets in multibeam sonar data and the transferability of trained models.
Genome-wide association studies (GWAS) at biobank scale are computationally intensive, especially for admixed populations requiring robust statistical models. SAIGE is a widely used method for generalized linear mixed-model GWAS but is limited by its CPU-based implementation, making phenome-wide association studies impractical for many research groups. We developed SAIGE-GPU, a GPU-accelerated version of SAIGE that replaces CPU-intensive matrix operations with GPU-optimized kernels. The core innovation is distributing genetic relationship matrix calculations across GPUs and communication layers. Applied to 2068 phenotypes from 635 969 participants in the Million Veteran Program, including diverse and admixed populations, SAIGE-GPU achieved a 5-fold speedup in mixed model fitting on supercomputing infrastructure and cloud platforms. We further optimized the variant association testing step through multi-core and multi-trait parallelization. Deployed on Google Cloud Platform and Azure, the method provided substantial cost and time savings. Source code and binaries are available for download at https://github.com/saigegit/SAIGE/tree/SAIGE-GPU-1.3.3. A code snapshot is archived at Zenodo for reproducibility (DOI: [10.5281/zenodo.17642591]). SAIGE-GPU is available in a containerized format for use across HPC and cloud environments and is implemented in R/C++ and runs on Linux systems.
Preconditioning is the most widely used and effective way for treating ill-conditioned linear systems in the context of classical iterative linear system solvers. We introduce a quantum primitive called fast inversion, which can be used as a preconditioner for solving quantum linear systems. The key idea of fast inversion is to directly block encode a matrix inverse through a quantum circuit implementing the inversion of eigenvalues via classical arithmetics. We demonstrate the application of preconditioned linear system solvers for computing single-particle Green's functions of quantum many-body systems, which are widely used in quantum physics, chemistry, and materials science. We analyze the complexities in three scenarios: the Hubbard model, the quantum many-body Hamiltonian in the plane-wave-dual basis, and the Schwinger model. We also provide a method for performing Green's function calculation in second quantization within a fixed-particle manifold and note that this approach may be valuable for simulation more broadly. Aside from solving linear systems, fast inversion also allows us to develop fast algorithms for computing matrix functions, such as the efficient preparation of Gibbs states. Furthermore, we introduce two efficient approaches for such a task, based on the contour-integral formulation and the inverse transform, respectively.
We present a new method to measure the vertical aerosol optical depth (VAOD) during clear nights using a wide-field imager—a CCD camera with a photographic lens on an equatorial mount. A series of 30 s exposures taken at different altitudes above the horizon can be used to measure the VAOD with a precision better than 0.008 optical depths within a few minutes. Such a measurement does not produce any light and is thus suitable for use at sites where other astronomical instruments are located. The precision of the VAOD measurement depends on the laboratory calibration of spectral properties of the system and of the response of the camera electronics to varying illumination levels, as well as careful considerations of the details of stellar photometry and modeling of the dependence of the measured stellar flux, the star color, and the position within the field of view. The results obtained with robotic setups at the future sites of the Cerenkov Telescope Array show good internal consistency and agreement with simultaneous measurements from a Sun/Moon photometer located at the same sites.