Search NASA⌕ Search

SEARCH · Search NASA

Results for “Cluster 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 1,711 records · Page 95

Π4U: A high performance computing framework for Bayesian uncertainty quantification of complex models

We present Π4U,{sup 1} an extensible framework, for non-intrusive Bayesian Uncertainty Quantification and Propagation (UQ+P) of complex and computationally demanding physical models, that can exploit massively parallel computer architectures. The framework incorporates Laplace asymptotic approximations as well as stochastic algorithms, along with distributed numerical differentiation and task-based parallelism for heterogeneous clusters. Sampling is based on the Transitional Markov Chain Monte Carlo (TMCMC) algorithm and its variants. The optimization tasks associated with the asymptotic approximations are treated via the Covariance Matrix Adaptation Evolution Strategy (CMA-ES). A modified subset simulation method is used for posterior reliability measurements of rare events. The framework accommodates scheduling of multiple physical model evaluations based on an adaptive load balancing library and shows excellent scalability. In addition to the software framework, we also provide guidelines as to the applicability and efficiency of Bayesian tools when applied to computationally demanding physical models. Theoretical and computational developments are demonstrated with applications drawn from molecular dynamics, structural dynamics and granular flow.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Multimodal biological analysis using NLP and expression profile

The goal of this project is to gather biological data from different sources and use computational analysis to evaluate it together. Two data sources were used: microarray gene expression data for Arabidopsis thaliana, and gene co-occurrences in scientific literature extracted from bioRxiv using natural language processing (NLP). For analysis, the microarray data was normalized, its dimensionality was reduced using principal component analysis (PCA), and it was grouped into different numbers of clusters using K-means clusters. Then these expression clusters were compared to the co-occurrence pairs in the NLP data, to evaluate the quality of the NLP extractions. This evaluation was done using entropy analysis on the combined data, compared to the maximum entropy in the clustering alone. As a result, the evaluation of the NLP data shows that the results do correspond to the clusters from the microarray data, and may be used for further analysis.

97 MATHEMATICS AND COMPUTING↗

Methodology for Clustering High-Resolution Spatiotemporal Solar Resource Data

In this report, we introduce a methodology to achieve multiple levels of spatial resolution reduction of solar resource data, with minimal impact on data variability, for use in energy systems modeling. The selection of an appropriate clustering algorithm, parameter selection including cluster size, methods of temporal data segmentation, and methods of cluster evaluation are explored in the context of a repeatable process. In describing this process, we illustrate the steps in creating a reduced resolution, but still viable, dataset to support energy systems modeling, e.g. capacity expansion or production cost modeling. This process is demonstrated through the use of a solar resource dataset; however, the methods are applicable to other resource data represented through spatiotemporal grids, including wind data. In addition to energy modeling, the techniques demonstrated in this paper can be used in a novel top-down approach to assess renewable resources within many other contexts that leverage variability in resource data but require reduction in spatial resolution to accommodate modeling or computing constraints.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Online Voltage Event Detection Using Synchrophasor Data with Structured Sparsity-Inducing Norms

This paper develops an accurate and computationally efficient data-driven framework to detect voltage events from PMU data streams. It develops an innovative Proximal Bilateral Random Projection (PBRP) algorithm to quickly decompose the PMU data matrix into a low-rank matrix, a row-sparse event-pattern matrix and a noise matrix. Here, the row-sparse pattern matrix significantly distinguishes events from normal behavior. These matrices are then fed into a clustering algorithm to separate voltage events from normal operating conditions. Large-scale numerical study results on real-world PMU data show that the proposed algorithm is computationally more efficient and achieves higher F scores than state-of-the-art benchmarks.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Scalable semidefinite programming approach to variational embedding for quantum many-body problems

In quantum embedding theories, a quantum many-body system is divided into localized clusters of sites which are treated with an accurate ‘high-level’ theory and glued together self-consistently by a less accurate ‘low-level’ theory at the global scale. The recently introduced variational embedding approach for quantum many-body problems combines the insights of semidefinite relaxation and quantum embedding theory to provide a lower bound on the ground-state energy that improves as the cluster size is increased. The variational embedding method is formulated as a semidefinite program (SDP), which can suffer from poor computational scaling when treated with black-box solvers. Here, we exploit the interpretation of this SDP as an embedding method to develop an algorithm which alternates parallelizable local updates of the high-level quantities with updates that enforce the low-level global constraints. Moreover, we show how translation invariance in lattice systems can be exploited to reduce the complexity of projecting a key matrix to the positive semidefinite cone.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Theory meets experiment: Insights into structure and magnetic properties of Fe 1-x Ni x B alloy

Here, we studied the structural and magnetic properties of the solid solution Fe 1-x Ni x B through theoretical and experimental approaches. Powder X-ray diffraction, X-ray Pair Distribution Function analysis, and energy dispersive X-ray spectroscopy reveal that the Fe 1-x Ni x B solid solution crystallizes in the β-FeB structure type up to x = 0.6–0.7 and exhibits anisotropic unit cell volume contraction with increasing Ni concentration. Magnetic measurements showed a transition from ferromagnetism to paramagnetism around x = 0.7. For x = 0.5, the low (< 0.3 μ B ) magnetic moments suggest itinerant magnetism despite the relatively high Curie temperature (up to 225 K). Theoretical calculations indicated different types of magnetic orderings depending on the Fe/Ni atomic order, with the antiferromagnetic state being stable for ordered FeNiB 2 , whereas the ground state is ferromagnetic for the disordered alloy. Calculations also predicted the coexistence of low- and high-spin states in Fe atoms around the composition with x = 0.5, in line with the experimental evidence from 57 Fe Mössbauer spectroscopy. The two magnetically distinct Fe sites for x = 0.3, 0.4, and 0.5 observed by 57 Fe Mössbauer spectroscopy can also be interpreted as two magnetically different regions or clusters. The formation of these clusters could affect the critical behavior near a quantum magnetic transition based on a potential ferromagnetic quantum critical point identified computationally and experimentally near x = 0.64. This work highlights the complex interplay between structure and magnetism in Fe 1-x Ni x B alloys, suggesting areas for future research on quantum critical behavior.

Bhaskar, Gourab [Iowa State Univ., Ames, IA (Unite↗

Quantum Simulations of Fermionic Hamiltonians with Efficient Encoding and Ansatz Schemes

Here, we propose a computational protocol for quantum simulations of fermionic Hamiltonians on a quantum computer, enabling calculations on spin defect systems which were previously not feasible using conventional encodings and a unitary coupled-cluster ansatz of variational quantum eigensolvers. We combine a qubit-efficient encoding scheme mapping Slater determinants onto qubits with a modified qubitcoupled cluster ansatz and noise-mitigation techniques. Our strategy leads to a substantial improvement in the scaling of circuit gate counts and in the number of required qubits, and to a decrease in the number of required variational parameters, thus increasing the resilience to noise. We present results for spin defects of interest for quantum technologies, going beyond minimum models for the negatively charged nitrogen vacancy center in diamonds and the double vacancy in 4H silicon carbide (4H-SiC) and tackling a defect as complex as negatively charged silicon vacancy in 4H-SiC for the first time.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Correlated pair ansatz with a binary tree structure

We develop an efficient algorithm to implement the recently introduced binary tree state (BTS) ansatz on a classical computer. BTS allows a simple approximation to permanents arising from the computationally intractable antisymmetric product of interacting geminals and respects size-consistency. We show how to compute BTS overlap and reduced density matrices efficiently. We also explore two routes for developing correlated BTS approaches: Jastrow coupled cluster on BTS and linear combinations of BT states. Here, the resulting methods show great promise in benchmark applications to the reduced Bardeen–Cooper–Schrieffer Hamiltonian and the one-dimensional XXZ Heisenberg Hamiltonian.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Closed-form expressions for unitaries of spin-adapted fermionic operators

One of the open challenges in quantum computing simulations of problems of chemical interest is the proper enforcement of spin symmetry. Efficient quantum circuits implementing unitaries generated by spin-adapted operators remain elusive, while naïve Trotterization schemes break spin symmetry. Here, in this work, we analyze the mathematical structure of spin-adapted operators and derive closed-form expressions for unitaries generated by singlet spin-adapted generalised single and double excitations. These results represent significant progress toward the economical enforcement of spin symmetry in quantum simulations.

fermionic algebra↗

Mechanistic Multiscale Uncertainty Propagation in Support of Accelerated Fuel Qualification

Taking a nuclear fuel concept through the research, development, and qualification stages has historically taken on the order of 20 to 25 years because of extensive irradiation tests required for a variety of conditions. The concept of accelerated fuel qualification (AFQ) has been proposed to increase the innovation pace for nuclear fuels. The goal of AFQ is not to replace the traditional qualification approach but rather to reduce the total number of experiments required to ensure approval from the regulatory authority. Of the many AFQ approaches being explored, advanced modeling—and, in particular, mechanistic modeling—is in a uniquely cross-cutting position to reduce the number of required integral tests through the inclusion of separate-effects testing, while helping to extrapolate reactor performance during rare events. We make the case that propagation of uncertainty through various computational length scales helps contextualize mechanistic modeling. We will utilize UO 2 fission gas diffusion predictions from the atomistically informed cluster dynamics code Centipede to inform fuel performance rodlet simulations using the BISON finite element code as the metric for showing how multiscale mechanistic uncertainty quantification can help reduce uncertainty in fuel performance. In conclusion, by quantifying uncertainty and its reduction through multiscale modeling, the qualification process may be accelerated through the reduction of costly irradiation experiments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The impact of anisotropic redshift distributions on angular clustering

A leading way to constrain physical theories from cosmological observations is to test their predictions for the angular clustering statistics of matter tracers, a technique that is set to become ever more central with the next generation of large imaging surveys. Interpretation of this clustering requires knowledge of the projection kernel, or the redshift distribution of the sources, and the typical assumption is an isotropic redshift distribution for the objects. However, variations in the kernel are expected across the survey footprint due to photometric variations and residual observational systematic effects. Here, we develop the formalism for anisotropic projection and present several limiting cases that elucidate the key aspects. We quantify the impact of anisotropies in the redshift distribution on a general class of angular two-point statistics. In particular, we identify a mode-coupling effect that can add power to auto-correlations, including galaxy clustering and cosmic shear, and remove it from certain cross-correlations. If the projection anisotropy is primarily at large scales, the mode-coupling depends upon its variance as a function of redshift; furthermore, it is often of similar shape to the signal. In contrast, the cross-correlation of a field whose selection function is anisotropic with another one featuring no such variations — such as CMB lensing — is immune to these effects. We discuss explicitly several special cases of the general formalism including galaxy clustering, galaxy-galaxy lensing, cosmic shear and cross-correlations with CMB lensing, and publicly release a code to compute the biases.

79 ASTRONOMY AND ASTROPHYSICS↗

An Optimal Parameterization Framework for Infrasonic Tomography of the Stratospheric Winds Using Non-Local Sources

A method is developed to apply acoustic tomography methods to a localized network of infrasound arrays with intention of monitoring the atmosphere state in the region around the network using non-local sources without requiring knowledge of the precise source location or non-local atmosphere state. Closely spaced arrays provide a means to estimate phase velocities of signals that can provide limiting bounds on certain characteristics of the atmosphere. Larger spacing between such clusters provide a means to estimate celerity from propagation times along multiple unique stratospherically or thermospherically ducted propagation paths and compute more precise estimates of the atmosphere state. In order to avoid the commonly encountered complex, multimodal distributions for parametric atmosphere descriptions and to maximize the computational efficiency of the method, an optimal parametrization framework is constructed. This framework identifies the ideal combination of parameters for tomography studies in specific regions of the atmosphere and statistical model selection analysis shows that high quality corrections to the middle atmosphere winds can be obtained using as few as three parameters. Lastly, comparison of the resulting estimates for synthetic data sets shows qualitative agreement between the middle atmosphere winds and those estimated from infrasonic traveltime observations.

58 GEOSCIENCES↗

Emulating ab initio computations of infinite nucleonic matter

We construct efficient emulators for the computation of the infinite nuclear matter equation of state. These emulators are based on the subspace-projected coupled-cluster method for which we here develop a new algorithm called small-batch voting to eliminate spurious states that might appear when emulating quantum many-body methods based on a non-Hermitian Hamiltonian. The efficiency and accuracy of these emulators facilitate a rigorous statistical analysis within which we explore nuclear matter predictions for > 10 6 different parametrizations of a chiral interaction model with explicit Δ -isobars at next-to-next-to leading order. Constrained by nucleon-nucleon scattering phase shifts and bound-state observables of light nuclei up to He 4 , we use history matching to identify nonimplausible domains for the low-energy coupling constants of the chiral interaction. Within these domains we perform a Bayesian analysis using sampling and importance resampling with different likelihood calibrations and study correlations between interaction parameters, calibration observables in light nuclei, and nuclear matter saturation properties. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A large-scale study of soft-errors on GPUs in the field

Parallelism provided by the GPU architecture has enabled domain scientists to simulate physical phenomena at a much faster rate and finer granularity than what was previously possible by CPU-based large-scale clusters. Architecture researchers have been investigating reliability characteristics of GPUs and innovating techniques to increase the reliability of these emerging computing devices. Such efforts are often guided by technology projections and simplistic scientific kernels, and performed using architectural simulators and modeling tools. Lack of large-scale field data impedes the effectiveness of such efforts. This study attempts to bridge this gap by presenting a large-scale field data analysis of GPU reliability. We characterize and quantify different kinds of soft-errors on the Titan supercomputer's GPU nodes. Our study uncovers several interesting and previously unknown insights about the characteristics and impact of soft-errors.

Nie, Bin↗

Asymmetrical Magnet Shape Optimization Based on S-C Mapping for Torque Profile Mitigation in Unidirectional Application of SPMS Machine

Surface-mounted Permanent Magnet Synchronous (SPMS) machines are widely used in applications that require low torque ripple, reduced vibration and acoustic noise. For those applications of SPMSM in which only unidirectional operation is needed, asymmetric design optimization show to have significant advantages over symmetric designs. In this paper, an analytical method based on Schwarz -Christoffel mapping is presented for asymmetrical magnet shape optimization in SPMS machine. In this method, permanent magnet is divided into finite cells so that each cell can be assigned with either air (OFF) or magnet (ON). The Optimization problem is solved using Genetic Algorithms (GAs). In order to remove non-manufacturable shapes, the cluster of materials and cleaning method is used and applied into magnet cells. The analytical method can significantly reduce the computation time used by optimization process as compared with methods based on Finite Element Analysis (FEA). This method is applied to a 12-pole SPMS machine. Validation of the S-C mapping method and final results are carried out using 2D and 3D finite element analysis. Here, the results obtained demonstrate that asymmetrical magnet shape improves cogging torque and torque ripple performance of SPMS machine by a significant margin.

42 ENGINEERING↗

The Attraction Indian Buffet Distribution

We propose the attraction Indian buffet distribution (AIBD), a distribution for binary feature matrices influenced by pairwise similarity information. Binary feature matrices are used in Bayesian models to uncover latent variables (i.e., features) that explain observed data. The Indian buffet process (IBP) is a popular exchangeable prior distribution for latent feature matrices. In the presence of additional information, however, the exchangeability assumption is not reasonable or desirable. The AIBD can incorporate pairwise similarity information, yet it preserves many properties of the IBP, including the distribution of the total number of features. Thus, much of the interpretation and intuition that one has for the IBP directly carries over to the AIBD. A temperature parameter controls the degree to which the similarity information affects feature-sharing between observations. Unlike other nonexchangeable distributions for feature allocations, the probability mass function of the AIBD has a tractable normalizing constant, making posterior inference on hyperparameters straight-forward using standard MCMC methods. A novel posterior sampling algorithm is proposed for the IBP and the AIBD. We demonstrate the feasibility of the AIBD as a prior distribution in feature allocation models and compare the performance of competing methods in simulations and an application.

97 MATHEMATICS AND COMPUTING↗

Critical Configuration and Physics Measurements for Assemblies of U(93.15)O 2 Fuel Rods (1.506-cm Pitch)

A series of critical experiments were completed from 1962–1965 at Oak Ridge National Laboratory’s (ORNL’s) Critical Experiments Facility (CEF) in support of the Medium-Power Reactor Experiments (MPRE) program. In the late 1950s, efforts were made to study “power plants for the production of electrical power in space vehicles.”(a) The MPRE program was a part of those efforts and studied the feasibility of a stainless-steel system, boiling potassium 1 MW(t), or about 140 kW(e), reactor. The program was carried out in [fiscal years] 1964, 1965, and 1966. A summary of the program’s effort was compiled in 1967.a The delayed critical experiments were a mockup of a small, potassium-cooled space power reactor for validation of reactor calculations and reactor physics methods. Initial experiments, performed in November and December of 1962, consisted of a core of unmoderated stainless-steel tubes, each containing 26 UO 2 fuel pellets, surrounded by a graphite reflector. Measurements were performed to determine critical reflector arrangements, relative fission-rate distributions, and cadmium ratio distributions. Subsequent experiments used beryllium reflectors and also measured the reactivity for various materials placed in the core. “The [assemblies were built] on [a] vertical assembly machine so that the movable part was the core and bottom reflector” (see Reference 1). The experiment studied in this evaluation was the second of the series and had the fuel rods in a 1.506-cm-triangular pitch. One critical configuration was found (see Reference 3). Once the critical configuration had been achieved, various measurements of reactivity, relative axial and radial activation rates of 235U,bc and cadmium ratios were performed. The cadmium ratio, reactivity, and activation rate measurements performed on the critical configuration are described in Sections 1.3, 1.4, and 1.7, respectively.

1.506-cm pitch↗

Critical Configuration and Physics Measurements for Beryllium Reflected Assemblies of U(93.15)O 2 Fuel Rods (1.506-cm Pitch)

Cadmium ratios were measured with enriched uranium metal foils at various locations in the assembly with the fuel tube at the 1.506-cm spacing. They are described in the following subsections. The experiment configuration was the same as the first critical configuration described in HEU-COMP-FAST-004 (Case 1). The experimenter placed 0.75-cm-diameter × 0.010-cm-thick 93.15%-235U-enriched uranium metal foils with and without 0.051-cm-thick cadmium covers at various locations in the core and top reflector. One part of the cadmium cover was cupshape and contained the uranium foil. The other part was a lid that fit over the exposed side of the foil when it was in the cup shaped section of the cover. As can be seen in the logbook, two runs were required to obtain all the measurements necessary for the cadmium ratio. The bare foil measurements within the top reflector were run first as part of the axial foil activation measurements. The results of this run are used for both the axial activation results and the cadmium ratios. Cadmium covered foils were then placed at the same location through the top reflector in a different run. Three pairs of bare and cadmium covered foils were also placed through the core tank. One pair was placed at the axial center of a fuel tube 11.35 cm from the center of the core. Two pairs of foils were placed on top of fuel tubes 3.02 and 12.06 cm from the center of the core. The activation of the uranium metal foils was measured after removal from the assembly using two lead shielded NaI scintillation detectors as follows. The NaI scintillators were carefully matched and had detection efficiencies for counting delayed-fission-product gamma rays with energies above 250 KeV within 5%. In all foil activation measurements, one foil at a specific location was used as a normalizing foil to remove the effects of the decay of fission products during the counting measurements with the NaI detectors. The normalization foil was placed on one NaI scintillator and the other foil on the other NaI detector and the activities measured simultaneously. The activation of a particular foil was compared to that of the normalization foil by dividing the count rate for each foil by that of the normalization foil. To correct for the differing efficiencies of the two NaI detectors, the normalization foil was counted in Detector 1 simultaneously with the foil at position x in Detector 2, and then the normalization foil was counted simultaneously in Detector 2 with the foil from position x in Counter 1. The activity of the foil from position x was divided by the activity of the normalization foil counted simultaneously. This resulted in obtaining two values of the ratio that were then averaged. This procedure essentially removed the effect of the differing efficiencies of the two NaI detectors. Differing efficiencies of 10% resulted in errors in the ratios measured to less than 1%. The background counting rates obatined with the foils used for the measurements on the NaI detectors before their irradiation measurement were subtracted from all count rates. The results of the cadmium ratio measurements are given in Table 1.3-1 and Figure 1.3-1. “No correction has been made for self shielding in the foils” (Reference 3).

1.506-cm pitch↗