Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “high performance and distributed 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 847 records · Page 47

Bayesian model updating with finite element vs surrogate models: Application to a miter gate structural system

Bayesian finite element (FE) model updating using direct model evaluations of large-scale high-fidelity FE models is extremely computationally expensive. Surrogate models can be used as fast emulators of FE models to accelerate the model calibration process. The physics/mechanics-based FE models are still the underpinning behind the surrogate models. Here, this paper evaluates the loss in accuracy and the gain in computational time while performing Bayesian model updating by using surrogate model evaluations compared to using direct FE model evaluations. This evaluation is crucial before entirely relying on surrogate models in model updating for structural health monitoring (SHM) and damage prognosis (DP) purposes. This paper also demonstrates Bayesian updating and surrogate model construction of large-scale high-fidelity FE models of infrastructure systems. In this regard, the miter gate structural system is considered as the testbed structure. Three predominant damage modes (loss of contact between gate and wall, loss of thickness due to corrosion, and loss of tension in the diagonal rods) are considered for model updating purposes. Bayesian model updating is performed using direct FE evaluations by leveraging parallel computing. Two types of surrogates, namely polynomial chaos expansion (PCE) and Gaussian process regression (GPR), are developed for the miter gate. Model updating is performed again using the trained surrogate models, and the updating results are compared with their counterparts obtained using the direct FE evaluation results. The posterior distribution of the FE model parameters obtained using the trained surrogates are sufficiently accurate with respect to the posterior obtained utilizing the direct FE evaluations. In addition, an approximate 4-fold decrease in the computational time was observed when using surrogate model evaluations instead of direct FE evaluations for model updating.

42 ENGINEERING↗

Modeling of the Tau and Muon Neutrino-induced Optical Cherenkov Signals from Upward-moving Extensive Air Showers

We present a detailed modeling and computation methodology to determine the optical Cherenkovsignals produced by upward-moving extensive air showers (EASs) induced byτ-leptons and muons,sourced from the interaction of high-energy astrophysical neutrinos interacting in the Earth. Fol-lowing and extending the physics modeling and Cherenkov signal simulations performed in [1], thisscheme encompasses a new, state-of-the-art computation of the muon neutrino propagation insidethe Earth and the contribution to theτ-lepton muon decay channel. The modeling takes into ac-count all possibleτ-lepton decay and muon energy loss channels that feed the optical Cherenkovemission, produced by both tau and muon initiated EASs. The EAS modeling uses the electronenergy, angular, and lateral distributions in the EAS and their evolution as well as the wavelengthdependence of the Cherenkov emission and its atmospheric attenuation. The results presented hereare focused on the detection capabilities of sub-orbital (balloon-borne) and orbital (satellite) basedinstruments. The latter case was calculated for POEMMAato compare to that presented in [1],specifically including the muon-decay channel ofτ-leptons and the muonic EAS Cherenkov signalfrom muon neutrino interactions in the Earth. By detailing all these individual contributions tothe optical Cherenkov emission and detection, we show how the ensemble that includes muonicchannels provides a large detection capability for space-based, high-energy cosmic neutrino detec-tion. Specifically, we show that for neutrino energies∼<10 PeV, the upward-EAS sensitivity due tomuon neutrino interactions in the Earth begin to dominate over that for tau neutrino interactions,effectively extending the neutrino sensitivity to lower energies.

A L Cummings↗

Framework and Method for Controlling a Robotic System Using a Distributed Computer Network

A robotic system for performing an autonomous task includes a humanoid robot having a plurality of compliant robotic joints, actuators, and other integrated system devices that are controllable in response to control data from various control points, and having sensors for measuring feedback data at the control points. The system includes a multi-level distributed control framework (DCF) for controlling the integrated system components over multiple high-speed communication networks. The DCF has a plurality of first controllers each embedded in a respective one of the integrated system components, e.g., the robotic joints, a second controller coordinating the components via the first controllers, and a third controller for transmitting a signal commanding performance of the autonomous task to the second controller. The DCF virtually centralizes all of the control data and the feedback data in a single location to facilitate control of the robot across the multiple communication networks.

Sanders, Adam M.↗

Addressing GPU memory limitations for Graph Neural Networks in High-Energy Physics applications

Introduction Reconstructing low-level particle tracks in neutrino physics can address some of the most fundamental questions about the universe. However, processing petabytes of raw data using deep learning techniques poses a challenging problem in the field of High Energy Physics (HEP). In the Exa.TrkX Project, an illustrative HEP application, preprocessed simulation data is fed into a state-of-art Graph Neural Network (GNN) model, accelerated by GPUs. However, limited GPU memory often leads to Out-of-Memory (OOM) exceptions during training, due to the large size of models and datasets. This problem is exacerbated when deploying models on High-Performance Computing (HPC) systems designed for large-scale applications. Methods We observe a high workload imbalance issue during GNN model training caused by the irregular sizes of input graph samples in HEP datasets, contributing to OOM exceptions. We aim to scale GNNs on HPC systems, by prioritizing workload balance in graph inputs while maintaining model accuracy. Our paper introduces diverse balancing strategies aimed at decreasing the maximum GPU memory footprint and avoiding the OOM exception, across various datasets. Results Our experiments showcase memory reduction of up to 32.14% compared to the baseline. We also demonstrate the proposed strategies can avoid OOM in application. Additionally, we create a distributed multi-GPU implementation using these samplers to demonstrate the scalability of these techniques on the HEP dataset. Discussion By assessing the performance of these strategies as data loading samplers across multiple datasets, we can gauge their effectiveness in both single-GPU and distributed environments. Our experiments, conducted on datasets of varying sizes and across multiple GPUs, broaden the applicability of our work to various GNN applications that handle input datasets with irregular graph sizes.

Lee, Claire Songhyun↗

Data-driven quantum approximate optimization algorithm for power systems

Abstract Quantum technology provides a ground-breaking methodology to tackle challenging computational issues in power systems. It is especially promising for Distributed Energy Resources (DERs) dominant systems that have been widely developed to promote energy sustainability. In those systems, knowing the maximum sections of power and data delivery is essential for monitoring, operation, and control. However, high computational effort is required. By leveraging quantum resources, Quantum Approximate Optimization Algorithm (QAOA) provides a means to search for these sections efficiently. However, QAOA performance relies heavily on critical parameters, especially for weighted graphs. Here we present a data-driven QAOA, which transfers quasi-optimal parameters between weighted graphs based on the normalized graph density. We verify the strategy with 39,774 expectation value calculations. Without parameter optimization, our data-driven QAOA is comparable with the Goemans-Williamson algorithm. This work advances QAOA and pilots its practical application to power systems in noisy intermediate-scale quantum devices.

Jing, Hang (ORCID:0000000256343489)↗

SU-D-206-01: Employing a Novel Consensus Optimization Strategy to Achieve Iterative Cone Beam CT Reconstruction On a Multi-GPU Platform

Purpose: While compressed sensing-based cone-beam CT (CBCT) iterative reconstruction techniques have demonstrated tremendous capability of reconstructing high-quality images from undersampled noisy data, its long computation time still hinders wide application in routine clinic. The purpose of this study is to develop a reconstruction framework that employs modern consensus optimization techniques to achieve CBCT reconstruction on a multi-GPU platform for improved computational efficiency. Methods: Total projection data were evenly distributed to multiple GPUs. Each GPU performed reconstruction using its own projection data with a conventional total variation regularization approach to ensure image quality. In addition, the solutions from GPUs were subject to a consistency constraint that they should be identical. We solved the optimization problem with all the constraints considered rigorously using an alternating direction method of multipliers (ADMM) algorithm. The reconstruction framework was implemented using OpenCL on a platform with two Nvidia GTX590 GPU cards, each with two GPUs. We studied the performance of our method and demonstrated its advantages through a simulation case with a NCAT phantom and an experimental case with a Catphan phantom. Result: Compared with the CBCT images reconstructed using conventional FDK method with full projection datasets, our proposed method achieved comparable image quality with about one third projection numbers. The computation time on the multi-GPU platform was ∼55 s and ∼ 35 s in the two cases respectively, achieving a speedup factor of ∼ 3.0 compared with single GPU reconstruction. Conclusion: We have developed a consensus ADMM-based CBCT reconstruction method which enabled performing reconstruction on a multi-GPU platform. The achieved efficiency made this method clinically attractive.

60 APPLIED LIFE SCIENCES↗

Multi-kw dc power distribution system study program

The first phase of the Multi-kw dc Power Distribution Technology Program is reported and involves the test and evaluation of a technology breadboard in a specifically designed test facility according to design concepts developed in a previous study on space vehicle electrical power processing, distribution, and control. The static and dynamic performance, fault isolation, reliability, electromagnetic interference characterisitics, and operability factors of high distribution systems were studied in order to gain a technology base for the use of high voltage dc systems in future aerospace vehicles. Detailed technical descriptions are presented and include data for the following: (1) dynamic interactions due to operation of solid state and electromechanical switchgear; (2) multiplexed and computer controlled supervision and checkout methods; (3) pulse width modulator design; and (4) cable design factors.

Berkery, E. A.↗

Computation of scalar far-field patterns of large-aperture antennas

In computer programs used for evaluating the performance of high-gain antennas, efficient numerical methods for calculating the far-field patterns must be used since the majority of computer time and storage requirements may be attributed to this phase of the program. The numerical method most frequently used is the Fast Fourier Transform (FFT), which computes the far field as the Fourier transform of the field distribution in the antenna aperture. A new numerical method that in many applications is superior to the FFT in terms of reducing computer time and storage requirements is described.

Omalley, T. A.↗

Parallel implementation and performance optimization of the configuration-interaction method

The configuration-interaction (CI) method, long a popular approach to describe quantum many-body systems, is cast as a very large sparse matrix eigenpair problem with matrices whose dimension can exceed one billion. Such formulations place high demands on memory capacity and memory bandwidth - - two quantities at a premium today. In this paper, we describe an efficient, scalable implementation, BIGSTICK, which, by factorizing both the basis and the interaction into two levels, can reconstruct the nonzero matrix elements on the fly, reduce the memory requirements by one or two orders of magnitude, and enable researchers to trade reduced resources for increased computational time. We optimize BIGSTICK on two leading HPC platforms - - the Cray XC30 and the IBM Blue Gene/Q. Specifically, we not only develop an empirically-driven load balancing strategy that can evenly distribute the matrix-vector multiplication across 256K threads, we also developed techniques that improve the performance of the Lanczos reorthogonalization. Combined, these optimizations improved performance by 1.3-8× depending on platform and configuration.

97 MATHEMATICS AND COMPUTING↗

A forward-swept wing configuration designed for high maneuverability by use of three-dimensional transonic theory

Supercritical technology has been applied to the design of a forward-swept-wing fighter configuration, and an assessment has been made of the relative performance of forward versus aft sweep. The wing and canard for this forward-swept wing configuration were designed for transonic maneuver by the use of a transonic computational analysis method and a transonic design procedure. The computational method calculates the transonic flow over a canard-wing-fuselage combination so that the strong transonic induced-flow effects of the canard on the wing are taken into account. A model of this configuration was constructed and was tested in the Langley 16-Foot Transonic Tunnel. The transonic theory gave a reasonably good estimate of the wing pressure distributions at transonic maneuver conditions. Comparison of this configuration with an equivalent aft-swept wing configuration showed that, at a Mach number of 0.9 and a lift coefficient of 0.9, the two configurations have essentially the same drag. This forward-swept wing configuration was also found to have very good maneuver performance relative to the Rockwell International HiMAT highly-maneuverable aircraft configuration.

Mann, M. J.↗

Transport of Space Environment Electrons: A Simplified Rapid-Analysis Computational Procedure

A computational procedure for describing transport of electrons in condensed media has been formulated for application to effects and exposures from spectral distributions typical of electrons trapped in planetary magnetic fields. The procedure is based on earlier parameterizations established from numerous electron beam experiments. New parameterizations have been derived that logically extend the domain of application to low molecular weight (high hydrogen content) materials and higher energies (approximately 50 MeV). The production and transport of high energy photons (bremsstrahlung) generated in the electron transport processes have also been modeled using tabulated values of photon production cross sections. A primary purpose for developing the procedure has been to provide a means for rapidly performing numerous repetitive calculations essential for electron radiation exposure assessments for complex space structures. Several favorable comparisons have been made with previous calculations for typical space environment spectra, which have indicated that accuracy has not been substantially compromised at the expense of computational speed.

Nealy, John E.↗

Developing and Testing a New Gas Distribution System for the Fluidized Bed of the Integrated Waste Treatment Unit (IWTU) Steam Reformer - 19219

The Integrated Waste Treatment Unit (IWTU) has been tasked with using steam reforming technology to treat 3.4 million liters of radioactive sodium-bearing waste (SBW) at the INL. Progress in developing steam reforming for SBW treatment has been hampered by numerous processing issues associated with the existing process equipment during simulant testing. The primary issue was high differential temperatures that developed in the fluidized bed of the denitration mineralization reformer (DMR). Based on technical evaluations and computer modeling, it was determined that the primary cause of poor fluidization in the DMR was the existing fluidized gas distributor. This distributor consisted of a combination of a ring header and four removable fluidizing gas rails. The combined ring header/rail system was unable to effectively fluidize the bed at the bottom of the DMR (around and below the ring header). As a result, a decision was made to design and test a new fluidized gas distributor. Conceptual designs focused on three different distributors: a modified design of the current system; a flat plate with bubble caps; and a conical bottom with dual plena. After testing the conceptual prototypes at Hazen Research Inc., the conical bottom distributor was chosen. A 30-day test with simulant SBW feed was then performed with the new distributor. During testing, a temperature differential of less than 2 deg. C was maintained in the DMR fluidized bed. The results supported the conclusion that poor fluidized gas distribution was the primary cause of high temperature differentials in the DMR fluidized bed. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Rapid QSTS Simulations for High-Resolution Comprehensive Assessment of Distributed PV

The rapid increase in penetration of distributed energy resources on the electric power distribution system has created a need for more comprehensive interconnection modeling and impact analysis. Unlike conventional scenario-based studies, quasi-static time-series (QSTS) simulations can realistically model time-dependent voltage controllers and the diversity of potential impacts that can occur at different times of year. However, to accurately model a distribution system with all its controllable devices, a yearlong simulation at 1-second resolution is often required, which could take conventional computers a computational time of 10 to 120 hours when an actual unbalanced distribution feeder is modeled. This computational burden is a clear limitation to the adoption of QSTS simulations in interconnection studies and for determining optimal control solutions for utility operations. The solutions we developed include accurate and computationally efficient QSTS methods that could be implemented in existing open-source and commercial software used by utilities and the development of methods to create high-resolution proxy data sets. This project demonstrated multiple pathways for speeding up the QSTS computation using new and innovative methods for advanced time-series analysis, faster power flow solvers, parallel processing of power flow solutions and circuit reduction. The target performance level for this project was achieved with year-long high-resolution time series solutions run in less than 5 minutes within an acceptable error.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Landscape rules predict optimal superhighways for the first peopling of Sahul

Archaeological data and demographic modelling suggest that the peopling of Sahul required substantial populations, occurred rapidly within a few thousand years and encompassed environments ranging from hyper-arid deserts to temperate uplands and tropical rainforests. How this migration occurred and how humans responded to the physical environments they encountered have, however, remained largely speculative. By constructing a high-resolution digital elevation model for Sahul and coupling it with fine-scale viewshed analysis of landscape prominence, least-cost pedestrian travel modelling and high-performance computing, here we create over 125 billion potential migratory pathways, whereby the most parsimonious routes traversed emerge. Our analysis revealed several major pathways—superhighways—transecting the continent, that we evaluated using archaeological data. These results suggest that the earliest Australian ancestors adopted a set of fundamental rules shaped by physiological capacity, attraction to visually prominent landscape features and freshwater distribution to maximize survival, even without previous experience of the landscapes they encountered.

54 ENVIRONMENTAL SCIENCES↗

Differentiable Earth mover’s distance for data compression at the high-luminosity LHC

Abstract The Earth mover’s distance (EMD) is a useful metric for image recognition and classification, but its usual implementations are not differentiable or too slow to be used as a loss function for training other algorithms via gradient descent. In this paper, we train a convolutional neural network (CNN) to learn a differentiable, fast approximation of the EMD and demonstrate that it can be used as a substitute for computing-intensive EMD implementations. We apply this differentiable approximation in the training of an autoencoder-inspired neural network (encoder NN) for data compression at the high-luminosity LHC at CERN The goal of this encoder NN is to compress the data while preserving the information related to the distribution of energy deposits in particle detectors. We demonstrate that the performance of our encoder NN trained using the differentiable EMD CNN surpasses that of training with loss functions based on mean squared error.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Using photoelectron spectroscopy to measure resonant inelastic X-ray scattering: a computational investigation

Resonant inelastic X-ray scattering (RIXS) has become an important scientific tool. Nonetheless, conventional high-resolution (few hundred meV or less) RIXS measurements, especially in the soft X-ray range, require low-throughput grating spectrometers, which limits measurement accuracy. Here, the performance of a different method for measuring RIXS, i.e. photoelectron spectrometry for analysis of X-rays (PAX), is computationally investigated. This method transforms the X-ray measurement problem of RIXS to an electron measurement problem, enabling use of high-throughput, compact electron spectrometers. X-rays to be measured are incident on a converter material and the energy distribution of the resultant photoelectrons, the PAX spectrum, is measured with an electron spectrometer. A deconvolution algorithm for analysis of such PAX data is proposed. It is shown that the deconvolution algorithm works well on data recorded with ~0.5 eV resolution. Additional simulations show the potential of PAX for estimation of RIXS features with smaller widths. For simulations using the 3d levels of Ag as a converter material, and with 10 5 simulated detected electrons, it is estimated that features with a few hundred meV width can be accurately estimated in a model RIXS spectrum. For simulations using a sharp Fermi edge to encode RIXS spectra, it is estimated that one can accurately distinguish 100 meV FWHM peaks separated by 45 meV with 10 5 simulated detected electrons that were photoemitted from within 0.4 eV of the Fermi level.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Benchmarking of the Updated Resolved Resonance Region Evaluations of Copper

A new resolved resonance region (RRR) evaluation for the two naturally occurring isotopes of copper, {sup 63}Cu and {sup 65}Cu was performed at Oak Ridge National Laboratory (ORNL), and the benchmarking results based on the four highly enriched uranium (HEU), intermediate spectrum, copper-reflected Zeus experiments were presented in June 2014. Since then, two significant updates to the evaluations have been performed: first, the data for the experimental measurement of the capture cross section between 220 keV and 300 keV were renormalized and the resonance parameters were refitted. Second, the angular distribution of elastic scattering initially reported in point-wise form in File 4 of the ENDF file format with high fidelity was significantly reduced in size to allow for better computer storage. Twenty-seven criticality safety benchmarks from the International Criticality Safety Benchmark Evaluation Project (ICSBEP) were calculated using these updated resolved resonance region evaluations for the two isotopes of copper. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Rotating reactor studies

Undesired gravitational effects such as convection or sedimentation in a fluid can sometimes be avoided or decreased by the use of a closed chamber uniformly rotated about a horizontal axis. In a previous study, the spiral orbits of a heavy or buoyant particle in a uniformly rotating fluid were determined. The particles move in circles, and spiral in or out under the combined effects of the centrifugal force and centrifugal buoyancy. A optimization problem for the rotation rate of a cylindrical reactor rotated about its axis and containing distributed particles was formulated and solved. Related studies in several areas are addressed. A computer program based on the analysis was upgraded by correcting some minor errors, adding a sophisticated screen-and-printer graphics capability and other output options, and by improving the automation. The design, performance, and analysis of a series of experiments with monodisperse polystyrene latex microspheres in water were supported to test the theory and its limitations. The theory was amply confirmed at high rotation rates. However, at low rotation rates (1 rpm or less) the assumption of uniform solid-body rotation of the fluid became invalid, and there were increasingly strong secondary motions driven by variations in the mean fluid density due to variations in the particle concentration. In these tests the increase in the mean fluid density due to the particles was of order 0.015 percent. To a first approximation, these flows are driven by the buoyancy in a thin crescent-shaped depleted layer on the descending side of the rotating reactor. This buoyancy distribution is balanced by viscosity near the walls, and by the Coriolis force in the interior. A full analysis is beyond the scope of this study. Secondary flows are likely to be stronger for buoyant particles, which spiral in towards the neutral point near the rotation axis under the influence of their centrifugal buoyancy. This is because the depleted layer is thicker and extends all the way around the reactor.

Roberts, Glyn O.↗