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 955 records · Page 53

SharP Hash: A High-Performing Distributed Hash for Extreme-Scale Systems

A high-performing distributed hash is critical for achieving performance in many applications and system software using extreme-scale systems. It is also a central part of many Big-Data frameworks including Memcached, file systems, and job schedulers. However, there is a lack of high-performing distributed hash implementations. In this work, we propose, design, and implement, SharP Hash, a high-performing, RDMA-based distributed hash for extreme-scale systems. SharP Hash's high performance is obtained through the use of high-performing networks and one-sided semantics. We perform an evaluation of SharP Hash and demonstrate its performance characteristics with a synthetic micro-benchmark and implementation of a Key Value (KV) store, Memcached.

97 MATHEMATICS AND COMPUTING↗

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.↗

A Simple Approach to the Technical Aspects of Radiosurgery Treatments

An approach to radiosurgery treatment that can be readily adopted in most radiotherapy centers with linear accelerators is presented. In our institution, a Leksell-type of neurosurgical frame, a computed tomography scanner, locally fabricated cones, and 6 MV X-ray beams are used to perform radiosurgery treatments. Collimated arcs with dose distributions that conform to the shape of the lesion in the transverse and the sagittal planes are used. It is argued that the uncertainties in the localization of the isocenter within a lesion and the specifications of the size of the target volume do not justify high precision mechanical devices for most radiosurgery treatments.

61 RADIATION PROTECTION AND DOSIMETRY↗

Data augmentation for disruption prediction via robust surrogate models

The goal of this work is to generate large statistically representative data sets to train machine learning models for disruption prediction provided by data from few existing discharges. Such a comprehensive training database is important to achieve satisfying and reliable prediction results in artificial neural network classifiers. Here, we aim for a robust augmentation of the training database for multivariate time series data using Student t process regression. We apply Student t process regression in a state space formulation via Bayesian filtering to tackle challenges imposed by outliers and noise in the training data set and to reduce the computational complexity. Thus, the method can also be used if the time resolution is high. We use an uncorrelated model for each dimension and impose correlations afterwards via colouring transformations. We demonstrate the efficacy of our approach on plasma diagnostics data of three different disruption classes from the DIII-D tokamak. To evaluate if the distribution of the generated data is similar to the training data, we additionally perform statistical analyses using methods from time series analysis, descriptive statistics and classic machine learning clustering algorithms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

PATCHY BLAZAR HEATING: DIVERSIFYING THE THERMAL HISTORY OF THE INTERGALACTIC MEDIUM

TeV-blazars potentially heat the intergalactic medium (IGM) as their gamma rays interact with photons of the extragalactic background light to produce electron–positron pairs, which lose their kinetic energy to the surrounding medium through plasma instabilities. This results in a heating mechanism that is only weakly sensitive to the local density, and therefore approximately spatially uniform, naturally producing an inverted temperature–density relation in underdense regions. In this paper we go beyond the approximation of uniform heating and quantify the heating rate fluctuations due to the clustered distribution of blazars and how this impacts the thermal history of the IGM. We analytically compute a filtering function that relates the heating rate fluctuations to the underlying dark matter density field. We implement it in the cosmological code GADGET-3 and perform large-scale simulations to determine the impact of inhomogeneous heating. We show that because of blazar clustering, blazar heating is inhomogeneous for z ≳ 2. At high redshift, the temperature–density relation shows an important scatter and presents a low temperature envelope of unheated regions, in particular at low densities and within voids. However, the median temperature of the IGM is close to that in the uniform case, albeit slightly lower at low redshift. We find that blazar heating is more complex than initially assumed and that the temperature–density relation is not unique. Our analytic model for the heating rate fluctuations couples well with large-scale simulations and provides a cost-effective alternative to subgrid models.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling of irradiation hardening of iron after low–dose and low–temperature neutron irradiation

Irradiation hardening is a prominent low-temperature degradation phenomena in materials, and is characterized both by an irradiation-induced increase in yield strength along with the loss of ductility. In this paper, a reaction–diffusion cluster dynamics model is used to predict the distribution of vacancy and interstitial clusters in iron following low-temperature (<373 K) and low-dose (<0.1 dpa) neutron irradiation. The predicted microstructure evolutions of high-purity iron samples are compared to published experimental data (positron annihilation spectroscopy and transmission electron microscopy) and show good agreement for neutron irradiation in this regime. The defect cluster distributions are then coupled to a dispersed barrier hardening model that assumes a strength factor, α, which varies with cluster type and size to compute the yield strength increase; the results of which agree reasonably well with tensile tests performed in previous studies. Furthermore, the modeling results presented here compare quite well to the experimental observations in the low-dose regime, and provide insight into the underlying microstructure–property relationships and the need for spatially dependent modeling to accurately predict the saturation behavior of yield strength changes observed experimentally at higher dose levels.

36 MATERIALS SCIENCE↗

Fed-DeepONet: Stochastic Gradient-Based Federated Training of Deep Operator Networks

The Deep Operator Network (DeepONet) framework is a different class of neural network architecture that one trains to learn nonlinear operators, i.e., mappings between infinite-dimensional spaces. Traditionally, DeepONets are trained using a centralized strategy that requires transferring the training data to a centralized location. Such a strategy, however, limits our ability to secure data privacy or use high-performance distributed/parallel computing platforms. To alleviate such limitations, in this paper, we study the federated training of DeepONets for the first time. That is, we develop a framework, which we refer to as Fed-DeepONet, that allows multiple clients to train DeepONets collaboratively under the coordination of a centralized server. To achieve Fed-DeepONets, we propose an efficient stochastic gradient-based algorithm that enables the distributed optimization of the DeepONet parameters by averaging first-order estimates of the DeepONet loss gradient. Then, to accelerate the training convergence of Fed-DeepONets, we propose a moment-enhanced (i.e., adaptive) stochastic gradient-based strategy. Finally, we verify the performance of Fed-DeepONet by learning, for different configurations of the number of clients and fractions of available clients, (i) the solution operator of a gravity pendulum and (ii) the dynamic response of a parametric library of pendulums.

Moya, Christian↗

Asymptotic/numerical analysis of supersonic propeller noise

An asymptotic analysis based on the Mach surface structure of the field of a supersonic helical source distribution is applied to predict thickness and loading noise radiated by high speed propeller blades. The theory utilizes an integral representation of the Ffowcs-Williams Hawkings equation in a fully linearized form. The asymptotic results are used for chordwise strips of the blade, while required spanwise integrations are performed numerically. The form of the analysis enables predicted waveforms to be interpreted in terms of Mach surface propagation. A computer code developed to implement the theory is described and found to yield results in close agreement with more exact computations.

Myers, M. K.↗

Thermodynamics, Kinetics and Structural Evolution of ε-LiVOPO 4 over Multiple Lithium Intercalation

ε -Li x VOPO 4 is a promising multi-electron cathode for rechargeable lithium-ion bat- teries that has an extremely high theoretical capacity of 318 mAh/g. In this w ork, we demonstrate the stable cycling of more than one Li in solid-state-syn thesized ε - LiVOPO 4 over more than 20 cycles for the first time. Using a combination of densit y functional theory (DFT) calculations, X-ray pair distribution funct ion (PDF) analy- sis and X-ray Absorption Near Edge Structure (XANES) measurements, we pre sent a comprehensive analysis of the thermodynamics, kinetics and stru ctural evolution of ε -Li x VOPO 4 over the entire lithiation range. We identify two intermediate ph ases at x = 1 . 5 and 1.75 in the low-voltage regime using DFT calculations, and the comput ed and electrochemically measured voltage profiles are in excellent agree ment. Operando PDF techniques show a reversible hysteretic change in the short ( < 2° A) V-O bond lengths coupled with an irreversible extension of the long V-O bond ( > 2.4 °A) dur- ing low-voltage cycling. These observations are confirmed with EXAFS spe ctra. We identify hydrogen intercalation from the electrolyte decomposition as a likely explana- tion for the ~ 2 . 4°A V-O bond and its irreversible extension. In terms of electronic conductivity, we find Li x VOPO 4 to be a large band gap insulator across the entire lithiation range, and calculated small polaron migration barriers are similar to those of the olivine LiMPO 4 cathodes. Finally, we demonstrate that ε -LiVOPO 4 is likely to be a one-dimensional diffuser using climbing-image nudged elastic ban d calculations. These results highlight the importance of nano-sizing and carbon coating in achieving good electrochemical performance in this material.

Lin, Yuh-Chieh↗

A Navier-Strokes Chimera Code on the Connection Machine CM-5: Design and Performance

We have implemented a three-dimensional compressible Navier-Stokes code on the Connection Machine CM-5. The code is set up for implicit time-stepping on single or multiple structured grids. For multiple grids and geometrically complex problems, we follow the 'chimera' approach, where flow data on one zone is interpolated onto another in the region of overlap. We will describe our design philosophy and give some timing results for the current code. A parallel machine like the CM-5 is well-suited for finite-difference methods on structured grids. The regular pattern of connections of a structured mesh maps well onto the architecture of the machine. So the first design choice, finite differences on a structured mesh, is natural. We use centered differences in space, with added artificial dissipation terms. When numerically solving the Navier-Stokes equations, there are liable to be some mesh cells near a solid body that are small in at least one direction. This mesh cell geometry can impose a very severe CFL (Courant-Friedrichs-Lewy) condition on the time step for explicit time-stepping methods. Thus, though explicit time-stepping is well-suited to the architecture of the machine, we have adopted implicit time-stepping. We have further taken the approximate factorization approach. This creates the need to solve large banded linear systems and creates the first possible barrier to an efficient algorithm. To overcome this first possible barrier we have considered two options. The first is just to solve the banded linear systems with data spread over the whole machine, using whatever fast method is available. This option is adequate for solving scalar tridiagonal systems, but for scalar pentadiagonal or block tridiagonal systems it is somewhat slower than desired. The second option is to 'transpose' the flow and geometry variables as part of the time-stepping process: Start with x-lines of data in-processor. Form explicit terms in x, then transpose so y-lines of data are in-processor. Form explicit terms in y, then transpose so z-lines are in processor. Form explicit terms in z, then solve linear systems in the z-direction. Transpose to the y-direction, then solve linear systems in the y-direction. Finally transpose to the x direction and solve linear systems in the x-direction. This strategy avoids inter-processor communication when differencing and solving linear systems, but requires a large amount of communication when doing the transposes. The transpose method is more efficient than the non-transpose strategy when dealing with scalar pentadiagonal or block tridiagonal systems. For handling geometrically complex problems the chimera strategy was adopted. For multiple zone cases we compute on each zone sequentially (using the whole parallel machine), then send the chimera interpolation data to a distributed data structure (array) laid out over the whole machine. This information transfer implies an irregular communication pattern, and is the second possible barrier to an efficient algorithm. We have implemented these ideas on the CM-5 using CMF (Connection Machine Fortran), a data parallel language which combines elements of Fortran 90 and certain extensions, and which bears a strong similarity to High Performance Fortran. We make use of the Connection Machine Scientific Software Library (CMSSL) for the linear solver and array transpose operations.

Jespersen, Dennis C.↗

Evaluating the effect of internal aperture variability on transport in kilometer scale discrete fracture networks

The apertures of natural fractures in fractured rock are highly heterogeneous. However, in-fracture aperture variability is often neglected in flow and transport modeling and individual fractures are assumed to have uniform aperture distribution. The relative importance of in-fracture variability in flow and transport modeling within kilometer18 scale field–scale fracture networks has been under a matter of debate for a long time because the flow in each single fracture is controlled not only by in-fracture variability but also by boundary conditions. Computational limitations have previously prohibited researchers from investigating the relative importance of in-fracture variability in flow and transport modeling within large-scale fracture networks. We address this question by incorporating internal heterogeneity of individual fractures into 23 flow simulations within kilometer scale three-dimensional fracture networks, where fracture intensity, P32 (ratio between total fracture area and domain volume) is between 0.027 and 0.031 [1/m]. A recently developed discrete fracture network (DFN) simulation capability, dfnWorks, is used to generate DFNs that include in-fracture aperture variability represented by a stationary log-normal stochastic field with various correlation lengths and variances. The Lagrangian transport parameters, non-reacting travel time and cumulative retention, are calculated along particles streamlines. It is observed that due to local flow channeling early particle travel times are more sensitive to in-fracture variability than the tails of travel time distributions, where no significant effect of the in-fracture transmissivity variations and spatial correlation length is observed.

54 ENVIRONMENTAL SCIENCES↗

Overall Traveling-Wave-Tube Efficiency Improved By Optimized Multistage Depressed Collector Design

Depressed Collector Design The microwave traveling wave tube (TWT) is used widely for space communications and high-power airborne transmitting sources. One of the most important features in designing a TWT is overall efficiency. Yet, overall TWT efficiency is strongly dependent on the efficiency of the electron beam collector, particularly for high values of collector efficiency. For these reasons, the NASA Glenn Research Center developed an optimization algorithm based on simulated annealing to quickly design highly efficient multistage depressed collectors (MDC's). Simulated annealing is a strategy for solving highly nonlinear combinatorial optimization problems. Its major advantage over other methods is its ability to avoid becoming trapped in local minima. Simulated annealing is based on an analogy to statistical thermodynamics, specifically the physical process of annealing: heating a material to a temperature that permits many atomic rearrangements and then cooling it carefully and slowly, until it freezes into a strong, minimum-energy crystalline structure. This minimum energy crystal corresponds to the optimal solution of a mathematical optimization problem. The TWT used as a baseline for optimization was the 32-GHz, 10-W, helical TWT developed for the Cassini mission to Saturn. The method of collector analysis and design used was a 2-1/2-dimensional computational procedure that employs two types of codes, a large signal analysis code and an electron trajectory code. The large signal analysis code produces the spatial, energetic, and temporal distributions of the spent beam entering the MDC. An electron trajectory code uses the resultant data to perform the actual collector analysis. The MDC was optimized for maximum MDC efficiency and minimum final kinetic energy of all collected electrons (to reduce heat transfer). The preceding figure shows the geometric and electrical configuration of an optimized collector with an efficiency of 93.8 percent. The results show the improvement in collector efficiency from 89.7 to 93.8 percent, resulting in an increase of three overall efficiency points. In addition, the time to design a highly efficient MDC was reduced from a month to a few days. All work was done in-house at Glenn for the High Rate Data Delivery Program. Future plans include optimizing the MDC and TWT interaction circuit in tandem to further improve overall TWT efficiency.

Vaden, Karl R.↗

High-resolution wall temperature measurements with distributed fiber optic sensors

Conventional approaches of surface temperature measurement are often unreliable, and other methods may not be practical in closed, high-temperature spaces. In this study, surface temperature measurements of the cooling panel of a natural circulation facility were performed with fiber optic distributed temperature sensors (DTS). High resolution surface temperature profiles of the cooling panel were acquired under two steady experimental conditions, and compared with measurements performed with nearby thermocouples. The total measurement error of the DTS results presented has been estimated to be ±4.4 °C. The findings reflected herein show that these DTS can reasonably measure surface temperature when installed with the methods developed in this study. Through the analysis of the DTS data collected, the interesting behavior of the cooling panel was revealed, providing new insights on the system's behavior that may be beneficial to the design and optimization of such heat transfer devices. These sensors can also provide the measurement density necessary for validation efforts with highly-resolved prediction methods such as finite element analysis (FEA) or computational fluid dynamics (CFD).

47 OTHER INSTRUMENTATION↗

Design consideration in constructing high performance embedded Knowledge-Based Systems (KBS)

As the hardware trends for artificial intelligence (AI) involve more and more complexity, the process of optimizing the computer system design for a particular problem will also increase in complexity. Space applications of knowledge based systems (KBS) will often require an ability to perform both numerically intensive vector computations and real time symbolic computations. Although parallel machines can theoretically achieve the speeds necessary for most of these problems, if the application itself is not highly parallel, the machine's power cannot be utilized. A scheme is presented which will provide the computer systems engineer with a tool for analyzing machines with various configurations of array, symbolic, scaler, and multiprocessors. High speed networks and interconnections make customized, distributed, intelligent systems feasible for the application of AI in space. The method presented can be used to optimize such AI system configurations and to make comparisons between existing computer systems. It is an open question whether or not, for a given mission requirement, a suitable computer system design can be constructed for any amount of money.

Dalton, Shelly D.↗

Excited State Chemistry in the Free Stream of the NASA IHF Arc Jet Facility Observed by Emission Spectroscopy

Spectroscopic measurements of non-equilibrium emission were made in the free stream of the 60 megawatts Interaction Heating Facility at NASA Ames Research Center. In the visible near infrared wavelength region, the most prominent emission was from molecular N2, and in the ultra violet region, the spectra were dominated by emission from molecular NO. The only atomic lines observed were those of copper (an erosion product of the electrodes). The bands of the 1st Positive system of N2 (if B is true then A is true) differed significantly from spectra computed spectra assuming only thermal excitation, suggesting overpopulation of the high vibrational states of the B state of N2. Populations of these high vibrational levels (peaking at v (sub upper) equals 13) of the N2 B state were determined by scaling simulated spectra; calculations were performed for each upper vibrational state separately. The experimental-theoretical procedure was repeated for several radial positions away from the nozzle axis to obtain spatial distributions of the upper state populations; rotational symmetry of the flow was assumed in simulations. The overpopulation of the high vibrational levels has been interpreted as the effect of inverse pre-dissociation of neutral atoms in the N2 A state, which populates the N2 B state through a level crossing process at v (sub upper) is greater than 10.

arc-jet testing↗

An algorithm for resolving intragranular orientation fields using coupled far-field and near-field high energy $\mathrm{X}$-ray diffraction microscopy

Here, we present a novel algorithm for reconstructing spatial intragranular lattice orientation fields using both far-field and near-field high energy X-ray diffraction microscopy (HEDM) measurements. An established far-field indexing algorithm is modified to include lattice orientation distribution information (grain orientation envelopes) in addition to average grain orientations. The near-field data reduction algorithm utilizes this enriched far-field orientation data as a seed for reconstructing intragranular spatial maps of lattice orientation from the diffraction images. The primary benefit of the new algorithm is a significant decrease in the number of trial calculations that must be performed in the reconstruction process compared to existing methodologies while maintaining scalability. The resulting gains in efficiency facilitate the use of relatively modest computational resources and improve throughput at the point of measurement. We provide two example applications: volumetric orientation field reconstructions for a Ti-Al alloy both before and after the application of 3% uniaxial strain. The results showcase the efficiency of the new method and the ability to resolve subtle changes in microstructure, which are associated with incipient plastic deformation.

36 MATERIALS SCIENCE↗

Accelerated coupled Monte Carlo-Thermal hydraulic calculations using a hybrid GTF-diffusion-based prediction block: first results

Accurate predictions of spatial power and temperature distributions require the coupling of a neutron transport solver with a thermal-hydraulic (TH) feedback. Nowadays, Monte Carlo (MC) codes are widely coupled to TH solvers, typically via a Picard iteration (PI) method, due to the higher fidelity that such frameworks can produce. To speed up a PI, a prediction step can produce an improved initial guess for a source distribution and feed it to the MC code. Recent work investigated a prediction step that uses generalized transfer functions (GTFs) to predict the macroscopic cross sections' variations following a perturbation in TH properties, such as coolant density. The previous method also relied on first order perturbation (FOP) theory to predict perturbed power profiles, rather than using an expensive MC iterate. The implemented FOP method relied on generating a fission matrix from which the forward and adjoint Eigenmodes were extracted and later used to for power calculations. The generation of the fission matrix can introduce a significant computational overhead, therefore undermining the performance of the proposed hybrid technique when applied to high-dimensional problems, e.g., full core calculations. This work attempts to improve the GTF-FOP prediction step by replacing the FOP solver with a nodal diffusion solver, thus eliminating the need to calculate a fission matrix. The GTF-diffusion step was tested for various moderator density perturbations. In each case, the predicted power distribution showed good agreement with the reference case. The latter is attributed to the generally good prediction of most spatially distributed macroscopic cross sections, except the transport cross section, which will become the focus of future work. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Impact of Time Dependent Reactor and Sensor Physics on Core Power Synthesis

Online synthesis of the power distribution is critical in the operation and control of nuclear power reactors to ensure that the core is operating within safety margins, and to provide essential knowledge associated with the burnup of the fuel. In light water reactors (LWRs), power synthesis is achieved by using some a priori knowledge of the state of the reactor core and updating based on the signals coming from in-core sensors—namely, self-powered neutron detectors (SPNDs). This report aims to study the effects of fuel burnup and sensor degradation on the ability to accurately synthesize the power distribution in a LWR. Several modeling tools were used to simulate power synthesis based on the responses of SPNDs, with emitters made out of Rh or V. A representative pressurized water reactor low-enriched uranium (LEU) core was modeled using the Polaris/Purdue Advanced Reactor Core Simulator (PARCS) approach. The Monte Carlo N-Particle Transport 6 (MCNP6) code was used, as well, to calculate response functions between different segments of fuel to individual SPNDs; this is a crucial parameter for power synthesis. The Oak Ridge Isotope GENeration (ORIGEN) package in the Standardized Computer Analyses for Licensing Evaluation (SCALE) code was used to model the time-dependent isotopic transmutation in the SPND emitters. All these data were fed into a custom code that enacted the point-based iterative (PBI) method to simulate power synthesis. Developmental work was also performed on high-fidelity SPND models in the GEometry ANd Tracking 4 (Geant4) code, which enables higher-accuracy modeling of the current responses from SPNDs. In this work, five sets of time-dependent power synthesis test cases were conducted. In these test cases, systematic changes in the input conditions enabled an analysis of the effect of (1) slightly inaccurate a priori power distribution assumptions with respect to fuel burnup, (2) highly inaccurate a priori assumptions with respect to fuel burnup (such that burnup is not included in the a priori assumed distribution), and (3) differences between Rh and V SPNDs in terms of downstream consequences of the transmutation in the emitters. The authors discovered that one may permissibly have slightly inaccurate a priori assumptions of the fuel burnup (such that the level of burnup may be slightly under- or over-approximated by the accumulated burnup in approximately 9.3 full power days), but to not account for burnup at all in the a priori assumption leads to severe levels of error, approaching 25% at maximum. The authors also discovered that V SPNDs are extraordinarily robust in the low-enriched uranium fuel cycle considered in this modeling work, whereas Rh SPNDs undergo significant transmutation that can result in large errors in the synthesized power distribution.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗