SEARCH · Engineering Papers
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.
The effect of defocusing due to the finite thickness of the crystal, in which parametric scattering occurs, on the spatial resolution of ghost images is considered. The necessary relationships are presented, and methods of levelling this distorting factor and improving the quality of ghost images are proposed. (paper)
We propose a reformulation of the streaming dynamic mode decomposition method that requires maintaining a single orthonormal basis, thereby reducing computational redundancy. The proposed efficient streaming dynamic mode decomposition method results in a constant-factor reduction in computational complexity and memory storage requirements. Numerical experiments on representative canonical dynamical systems show that the enhanced computational efficiency does not compromise the accuracy of the proposed method.
Abstract We discuss the use of Gaussian random fields to estimate the look-elsewhere effect correction. We show that Gaussian random fields can be used to model the null-hypothesis significance maps from a large set of statistical problems commonly encountered in physics, such as template matching and likelihood ratio tests. Some specific examples are searches for dark matter using pixel arrays, searches for astronomical transients, and searches for fast-radio bursts. Gaussian random fields can be sampled efficiently in the frequency domain, and the excursion probability can be fitted with these samples to extend any estimation of the look-elsewhere effect to lower p values. In addition, in cases where the Gaussian random field is stationary and the parameter space is Euclidean, the look-elsewhere effect correction can be computed analytically. We demonstrate these methods using two example template matching problems. Finally, we apply these methods to estimate the trial factor of a $$4^3$$ 4 3 accelerometer array for the detection of dark matter tracks in the Windchime project. When a global significance of $$3\sigma $$ 3 σ is required, the estimated trial factor for such an accelerometer array is $$10^{14}$$ 10 14 for a one-second search, and $$10^{22}$$ 10 22 for a 1-year search.
Joint Nonnegative Matrix Factorization (JointNMF) is a hybrid method for mining information from datasets that contain both feature and connection information. We propose distributed-memory parallelizations of three algorithms for solving the JointNMF problem based on Alternating Nonnegative Least Squares, Projected Gradient Descent, and Projected Gauss-Newton. We extend well-known communication-avoiding algorithms using a single processor grid case to our coupled case on two processor grids. We demonstrate the scalability of the algorithms on up to 960 cores (40 nodes) with 60% parallel efficiency. The more sophisticated Alternating Nonnegative Least Squares (ANLS) and Gauss-Newton variants outperform the first-order gradient descent method in reducing the objective on large-scale problems. We perform a topic modelling task on a large corpus of academic papers that consists of over 37 million paper abstracts and nearly a billion citation relationships, demonstrating the utility and scalability of the methods.
Cross section self-shielding methodologies for TRISO fuel were assessed to provide accurate multigroup cross sections for a high-fidelity reactor physics code so that the code is able to accurately model and simulate advanced reactors with TRISO fuel. Initially, the two existing methodologies (the SCALE method and the Sanchez-Pomraning method) were studied and implemented to MC2-3 for detailed performance tests. Additionally, a new spatial self-shielding method, named the iterative local spatial self-shielding (ILSS) method, for particulate fuels was developed based on the disadvantage factor and implemented to MC2-3 as well. The new method approximately accounts for the effect of randomly distributed particles on the particle shadowing effect using a homogenized compact region surrounding a particle of interest at the center. The self-shielded cross sections of the particle at the center are determined iteratively since the cross sections of the homogenized compact region are calculated using them. For the energy range above 100 keV where the fuel-to-moderator ratio is more important than the random distribution of particles, a single particle unit-cell model is used by preserving the average amount of moderator per fuel particle in the system. The three self-shielding methods implemented in MC2-3 were tested using numerical benchmark problems made based on fuel compact problems of a prismatic-type very high temperature reactor. Test results indicated that the ILSS method produced slightly better results than the SCALE and Sanchez-Pomraning methods, compared to the Serpent-2 Monte Carlo results obtained with 25 independent random particle configurations. The SCALE and Sanchez-Pomraning methods tend to underestimate the heterogeneity effect by 150 and 100 pcm, respectively, while the new ILSS method overestimates the heterogeneity effect by 70 pcm. In future, the new self-shielding method will be extended to perform pebble calculations and compare results with those from the SCALE and Sanchez-Pomraning methods. Furthermore, the new method will be optimized for practical applications to on-the-fly resonance treatment for lattice or whole-core calculations for advanced reactors.
This document provides a technical basis for establishing a Dosimetry-based Area Monitoring Program (DAMP) at Lawrence Livermore National Laboratory (LLNL). A DAMP is part of a comprehensive routine monitoring program and provides information about radiation levels inside and outside of Radiologically-Controlled Areas (RCAs). The routine, hand-held radiation survey program driven by the Health Physics Discipline Action Plan (HP-DAP) helps to ensure radiological conditions within RCAs are well-characterized and routinely monitored; however, such surveys are a snapshot in time, whereas area monitoring dosimeters (AMDs) continuously monitor radiation doses in the areas where they are installed. Together, data from the hand-held radiation surveys and the DAMP provides a comprehensive picture of radiation environment at LLNL. A DAMP helps to verify the effectiveness of established engineered and administrative controls while documenting that radiation doses in RCAs are below that which requires individual monitoring. This document establishes the basis for determining which model of dosimeter to use, the exchange frequency, the occupancy factor, monitoring locations, and a method for evaluating AMD results.
Randomized algorithms have propelled advances in artificial intelligence (AI) and represent a foundational research area in advancing AI for Science. Future advancements in DOE Office of Science priority areas such as climate science, astrophysics, fusion, advanced materials, combustion, and quantum computing all require randomized algorithms for surmounting challenges of complexity, robustness, and scalability. Advances in data collection and numerical simulation have changed the dynamics of scientific research and motivate the need for randomized algorithms. For instance, advances in imaging technologies such as X-ray ptychography, electron microscopy, electron energy loss spectroscopy, or adaptive optics lattice light-sheet microscopy collect hyperspectral imaging and scattering data in terabytes, at breakneck speed enabled by state-of-the-art detectors. The data collection is exceptionally fast compared with its analysis. Likewise, advances in high-performance architectures have made exascale computing a reality and changed the economies of scientific computing in the process. Floating-point operations that create data are essentially free in comparison with data movement. Thus far, most approaches have focused on creating faster hardware. Ironically, this faster hardware has exacerbated the problem by making data still easier to create. Under such an onslaught, scientists often resort to heuristic deterministic sampling schemes (e.g., low-precision arithmetic, sampling every nth element) and sacrifice potentially valuable accuracy. Dramatically better results can be achieved via randomized algorithms, reducing the data size as much as or more than naive deterministic subsampling can achieve, while retaining the high accuracy of computing on the full data set. By randomized algorithms we mean those algorithms that employ some form of randomness in internal algorithmic decisions to accelerate time to solution, increase scalability, or improve reliability. Examples include matrix sketching for solving large-scale least-squares problems (see Figure 1) and stochastic gradient descent for training machine learning models. We are not recommending heuristic methods but rather randomized algorithms that have certificates of correctness and probabilistic guarantees of optimality and near-optimality. Such approaches can be useful beyond acceleration, for example, in understanding how to avoid measure zero worst-case scenarios that plague methods such as QR matrix factorization.
This Core Model Proposal adds fugitive CO2 emissions from coal, crude oil, and natural gas resources to the data system. It also includes an improved method to replace outlier emissions factors for fossil resources. Lastly, it makes some minor changes to the data system and Model Interface headers to fix some inconsistencies in how resource technologies are named in the XMLs.
Low Temperature Detectors have been used to measure embedded radioisotopes in a measurement mode known as Decay Energy Spectroscopy (DES) since 1992. DES microcalorimeter measurements have been used for applications ranging from neutrino mass measurements to metrology to measurements for safeguards and medical nuclides. While the low temperature detectors have extremely high intrinsic energy resolution (several times better than semiconductor detectors), the energy resolution achieved in practice is strongly dependent on factors such as sample preparation method. This review seeks to present the literature consensus on what has been learned by looking at the energy resolution as a function of various choices of detector, absorber, and sample preparation methods.
This paper aims to forecast high-resolution (hourly) aggregated load for a certain region in the medium term (a few days to over a year). One region is defined as some places with similar climate characteristics because the climate influences people's daily lifestyles and hence the electric usage. We decompose the electric usage records into two parts: base load and seasonal load. Considering both temperature and time factors, different deep-learning methods are adopted to characterize them. The first goal of our approach is to predict the peak load which is critical for power system planning. Furthermore, our proposed forecast method can provide the depiction of the hourly load profile to provide customized load curves for high-level real-time applications. The proposed method is tested on real-world historical data collected by CAISO, BPA, and PACW. The experimental results show that trained by three years of data, our method could reduce the prediction error for a one-year lead hourly load below $5\%$ MAPE, and predict the occurrence of the peak load for next year in CAISO with an error within three days. Furthermore, as a byproduct, an interesting observation on the impact of COVID-19 on human life was made and discussed based on these case studies.
This paper aims to forecast high-resolution (hourly) aggregated load for a certain region in the medium term (a few days to over a year). One region is defined as some places with similar climate characteristics because the climate influences people's daily lifestyles and hence the electric usage. We decompose the electric usage records into two parts: base load and seasonal load. Considering both temperature and time factors, different deep-learning methods are adopted to characterize them. The first goal of our approach is to predict the peak load which is critical for power system planning. Furthermore, our proposed forecast method can provide the depiction of the hourly load profile to provide customized load curves for high-level real-time applications. The proposed method is tested on real-world historical data collected by CAISO, BPA, and PACW. The experimental results show that trained by three years of data, our method could reduce the prediction error for a one-year lead hourly load below 5% MAPE, and predict the occurrence of the peak load for next year in CAISO with an error within three days. Furthermore, as a byproduct, an interesting observation on the impact of COVID-19 on human life was made and discussed based on these case studies.
The ASME AG-1 Code specifies minimum rated airflow based on standard cubic feet per minute (SCFM) for axial high efficiency particulate air (HEPA) filters and actual cubic feet per minute (ACFM) for radial HEPA filters with corresponding maximum pressure drop. A quantitative analysis can be performed to demonstrate the variation in filter efficiency and develop a correlation of performance between axial filters subjected to ACFM and SCFM. This study illustrates the variation in filter efficiency of axial HEPA filters subjected to ACFM and SCFM flow rate conditions at varying atmospheric pressures. It also presents the methods used to calculate correction factors to convert filter efficiency and pressure drop from elevated to standard conditions. The Axial Large Scale Test Stand (ALSTS) in the Institute for Clean Energy Technology (ICET) at Mississippi State University has been utilized to control aerosol concentration, flowrate, atmospheric pressure, relative humidity, and temperature. The filter test matrix consists of axial HEPA filters rated at 1,000 and 1,500 CFM with varying filter pack types. Each filter is subjected to varying pressures of 14.70 psi, 13.55 psi, 12.40 psi, and 11.25 psi, temperatures of 70 deg. F, 110 deg. F, and 150 deg. F, and relative humidities of 50% and 80%, resulting in twenty-four unique testing points. The corresponding SCFM flow rates are computed for each point in order to form a comparison. Evaluation of the effects of designating ACFM or SCFM will be based upon the filter efficiencies of generated Dioctyl Phthalate (DOP) particles at a diameter of 0.3 micrometers. The results show filter efficiency decreasing as atmospheric pressure decreases due to increased filter face velocity. A correlation model is developed to convert filter efficiency and pressure drop at elevated conditions to those of standard pressure, relative humidity, and temperature.
We perform first-principle calculations of electron-nucleus scattering on 3 He and 3 H using the Green's function Monte Carlo method and two approaches based on the factorization of the final hadronic state: the spectral-function formalism and the short-time approximation. These three methods are benchmarked among each other and compared to the experimental data for the longitudinal and transverse electromagnetic response functions of 3 He, and the inclusive cross sections of both 3 He and 3 H. Since these three approaches are based on the same description of nuclear dynamics of the initial target state, comparing their results enables a precise quantification of the uncertainties inherent to factorization schemes. At sufficiently large values of the momentum transfer, we find an excellent agreement of the Green's function Monte Carlo calculation with experimental data and with both the spectral-function formalism and the short-time approximation. Here, we also analyze the relevance of relativistic effects, whose inclusion becomes crucial to explain data at high momentum and energy transfer.
We perform first-principle calculations of electron-nucleus scattering on $^3$He and $^3$H using the Green's function Monte Carlo method and two approaches based on the factorization of the final hadronic state: the spectral-function formalism and the short-time approximation. These three methods are benchmarked among each other and compared to the experimental data for the longitudinal and transverse electromagnetic response functions of $^3$He, and the inclusive cross sections of both $^3$He and $^3$H. Since these three approaches are based on the same description of nuclear dynamics of the initial target state, comparing their results enables a precise quantification of the uncertainties inherent to factorization schemes. At sufficiently large values of the momentum transfer, we find an excellent agreement of the Green's function Monte Carlo calculation with experimental data and with both the spectral-function formalism and the short-time approximation. We also analyze the relevance of relativistic effects, whose inclusion becomes crucial to explain data at high momentum and energy transfer.
The standard LU factorization-based solution process for linear systems can be enhanced in speed or accuracy by employing mixed-precision iterative refinement. Most recent work has focused on dense systems. We investigate the potential of mixed-precision iterative refinement to enhance methods for sparse systems based on approximate sparse factorizations. In doing so, we first develop a new error analysis for LU- and GMRES-based iterative refinement under a general model of LU factorization that accounts for the approximation methods typically used by modern sparse solvers, such as low-rank approximations or relaxed pivoting strategies. We then provide a detailed performance analysis of both the execution time and memory consumption of different algorithms, based on a selected set of iterative refinement variants and approximate sparse factorizations. Our performance study uses the multifrontal solver MUMPS, which can exploit block low-rank factorization and static pivoting. We evaluate the performance of the algorithms on large, sparse problems coming from a variety of real-life and industrial applications showing that mixed-precision iterative refinement combined with approximate sparse factorization can lead to considerable reductions of both the time and memory consumption.
Accurate synthetic spectra that rely on large Line-By-Line (LBL)-databases are used in a wide range of applications such as high temperature combustion, atmospheric re-entry, planetary surveillance and laboratory plasmas. Conventionally synthetic spectra are calculated by computing a lineshape for every spectral line in the database and adding those together, which may take multiple hours for large databases. In this paper we propose a new approach for spectral synthesis based on an integral transform: the synthetic spectrum is calculated as the integral over the product of a Voigt profile and a newly proposed three-dimensional “lineshape distribution function”, which is a function of spectral position and Gaussian- & Lorentzian width coordinates. A fast discrete version of this transform based on the Fast Fourier Transform (FFT) is proposed, which improves performance compared to the conventional approach by several orders of magnitude while maintaining accuracy. Strategies that minimize the discretization error are discussed. A Python implementation of the method is compared against state-of-the-art spectral code RADIS, and is since adopted as RADIS's default synthesis method. The synthesis of a benchmark CO2 spectrum consisting of 1.8 M spectral lines and 200k spectral points took only 3.1 s using the proposed method (1011 lines × spectral points/s), a factor ~300 improvement over the state-of-the-art, with the relative improvement generally increasing for higher number of lines and/or number of spectral points. Finally, an experimental GPU-implementation of the method was also benchmarked, which demonstrated another 2~3 orders performance increase, achieving up to 5 ∙ 10 14 lines × spectral points/s.
The accurate and reliable prediction of protein–ligand binding affinities can play a central role in the drug discovery process as well as in personalized medicine. Of considerable importance during lead optimization are the alchemical free energy methods that furnish an estimation of relative binding free energies (RBFE) of similar molecules. Recent advances in these methods have increased their speed, accuracy, and precision. This is evident from the increasing number of retrospective as well as prospective studies employing them. However, such methods still have limited applicability in real-world scenarios due to a number of important yet unresolved issues. Here, we report the findings from a large data set comprising over 500 ligand transformations spanning over 300 ligands binding to a diverse set of 14 different protein targets which furnish statistically robust results on the accuracy, precision, and reproducibility of RBFE calculations. We use ensemble-based methods which are the only way to provide reliable uncertainty quantification given that the underlying molecular dynamics is chaotic. These are implemented using TIES (Thermodynamic Integration with Enhanced Sampling). Results achieve chemical accuracy in all cases. Ensemble simulations also furnish information on the statistical distributions of the free energy calculations which exhibit non-normal behavior. We find that the “enhanced sampling” method known as replica exchange with solute tempering degrades RBFE predictions. We also report definitively on numerous associated alchemical factors including the choice of ligand charge method, flexibility in ligand structure, and the size of the alchemical region including the number of atoms involved in transforming one ligand into another. Our findings provide a key set of recommendations that should be adopted for the reliable application of RBFE methods.
Adsorbate–adsorbate lateral interactions at relevant surface coverages have a significant effect on chemical kinetics, thereby influencing the activity of a heterogeneous catalyst. Coverage-dependent kinetic and thermodynamic parameters therefore must be included in studies of such complex systems to properly predict the turnover frequencies and kinetic trends. Thus, it becomes extremely important to accurately capture the strength of lateral interactions between neighboring species under realistic reaction conditions. In this Perspective, we discuss the various existing computational and experimental methods for determining adspecies coverage and configurational effects. The choice of the tools and methods employed in such studies depends on factors such as time, length scales, computational cost, the presence of solvents, and reaction conditions. The applications of each method and the respective challenges are also discussed here. As a result, we discuss the recent developments and future of the state-of-the-art for inclusion of surface coverage and configuration into a holistic picture for accurate predictions of catalytic behavior.