Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distributed computation”

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,675 records · Page 93

UPC++ Programmer’s Guide, v1.0-2018.3.0

UPC++ is a C++11 library that provides Partitioned Global Address Space (PGAS) programming. It is designed for writing parallel programs that run efficiently and scale well on distributed-memory parallel computers. The PGAS model is single program, multiple-data (SPMD), with each separate thread of execution (referred to as a rank, a term borrowed from MPI) having access to local memory as it would in C++. However, PGAS also provides access to a global address space, which is allocated in shared segments that are distributed over the ranks. UPC++ provides numerous methods for accessing and using global memory. In UPC++, all operations that access remote memory are explicit, which encourages programmers to be aware of the cost of communication and data movement. Moreover, all remote-memory access operations are by default asynchronous, to enable programmers to write code that scales well even on hundreds of thousands of cores.

97 MATHEMATICS AND COMPUTING↗

Model Reduction for Fast Assessment of Grid Impact of High Penetration PV

A key challenge in properly analyzing the impact of solar photovoltaics (PV) during planning, interconnection and operation in the distribution grid is computational complexity. Such studies require the use of Quasi-Steady-Time-Series (QSTS) power flow to capture the complex interaction of PV with voltage regulating equipment and protection devices and a large number of scenario studies are needed to capture location specific constraints and uncertainties. In addition, a typical distribution feeder can have thousands of nodes and there can be thousands of such feeders in a given utility system. To solve this problem, a new and systematic method for distribution network reduction is proposed. The objective is to reduce the size of distribution systems by more than 80% while preserving the accuracy of power flow results. To this end, an extended Ward-equivalent method has been proposed to create an equivalent circuit. Moreover, load and PV aggregation based on the backward sweep method has been developed to further simplify the feeder. The proposed method is tested with a modified IEEE 123-node distribution test system. The simulation results demonstrate the high efficiency and accuracy of the proposed approach for distribution network reduction.

reduced order systems, load flow, load modeling, v↗

Evaluation of a strain-gage load calibration on a low-aspect-ratio wing structure at elevated temperature

The environmental aspect of elevated temperature and its relationship to the science of strain gage calibrations of aircraft structures are addressed. A section of a wing designed for a high-speed aircraft structure was used to study this problem. This structure was instrumented with strain gages calibrated at both elevated and room temperatures. Load equations derived from a high-temperature load calibration were compared with equations derived from an identical load calibration at room temperature. The implications of the high temperature load calibration were studied from the viewpoint of applicability and necessity. Load equations derived from the room temperature load calibration resulted in generally lower equation standard errors than equations derived from the elevated temperature load calibration. A distributed load was applied to the structure at elevated temperature and strain gage outputs were measured. This applied load was then calculated using equations derived from both the room temperature and elevated temperature calibration data. It was found that no significant differences between the two equation systems existed in terms of computing this applied distributed load, as long as the thermal shifts resulting from thermal stresses could be identified. This identification requires a heating of the structure. Therefore, it is concluded that for this structure, a high temperature load calibration is not required. However, a heating of the structure is required to determine thermal shifts.

Reardon, Lawrence F.↗

High-Performance Monitoring Architecture for Large-Scale Distributed Systems Using Event Filtering

Monitoring is an essential process to observe and improve the reliability and the performance of large-scale distributed (LSD) systems. In an LSD environment, a large number of events is generated by the system components during its execution or interaction with external objects (e.g. users or processes). Monitoring such events is necessary for observing the run-time behavior of LSD systems and providing status information required for debugging, tuning and managing such applications. However, correlated events are generated concurrently and could be distributed in various locations in the applications environment which complicates the management decisions process and thereby makes monitoring LSD systems an intricate task. We propose a scalable high-performance monitoring architecture for LSD systems to detect and classify interesting local and global events and disseminate the monitoring information to the corresponding end- points management applications such as debugging and reactive control tools to improve the application performance and reliability. A large volume of events may be generated due to the extensive demands of the monitoring applications and the high interaction of LSD systems. The monitoring architecture employs a high-performance event filtering mechanism to efficiently process the large volume of event traffic generated by LSD systems and minimize the intrusiveness of the monitoring process by reducing the event traffic flow in the system and distributing the monitoring computation. Our architecture also supports dynamic and flexible reconfiguration of the monitoring mechanism via its Instrumentation and subscription components. As a case study, we show how our monitoring architecture can be utilized to improve the reliability and the performance of the Interactive Remote Instruction (IRI) system which is a large-scale distributed system for collaborative distance learning. The filtering mechanism represents an Intrinsic component integrated with the monitoring architecture to reduce the volume of event traffic flow in the system, and thereby reduce the intrusiveness of the monitoring process. We are developing an event filtering architecture to efficiently process the large volume of event traffic generated by LSD systems (such as distributed interactive applications). This filtering architecture is used to monitor collaborative distance learning application for obtaining debugging and feedback information. Our architecture supports the dynamic (re)configuration and optimization of event filters in large-scale distributed systems. Our work represents a major contribution by (1) survey and evaluating existing event filtering mechanisms In supporting monitoring LSD systems and (2) devising an integrated scalable high- performance architecture of event filtering that spans several kev application domains, presenting techniques to improve the functionality, performance and scalability. This paper describes the primary characteristics and challenges of developing high-performance event filtering for monitoring LSD systems. We survey existing event filtering mechanisms and explain key characteristics for each technique. In addition, we discuss limitations with existing event filtering mechanisms and outline how our architecture will improve key aspects of event filtering.

Maly, K.↗

Bulk properties and velocity distributions of water group ions at Comet Halley - Giotto measurements

In the region upstream of Comet Halley, pickup heavy ions of cometary origin were directly observed by the implanted ion spectrometer on Giotto. Diffusion of this population in pitch angle and in energy, during the approach to the comet and on the outbound leg is discussed. The two data sets are compared and qualitative ideas on scattering timescales are inferred. In addition the bulk parameters of these distributions have been computed and a comparison of the observed speed in the solar wind frame and the observed density with expectations is presented. Pitch angle scattering occurs more slowly than expected with filled shells appearing at 2,500,000 km, and significant energy diffusion does not occur until the bow shock region. Also the shell distributions downstream of the shock flow at the bispherical bulk speed (related to the Alfven speed) along the magnetic field with respect to the solar wind in accordance with conservation of energy between the pickup ions and the wave turbulence.

Coates, A. J.↗

VERA and VERA-EDU 3.5 Release Notes

The Virtual Environment for Reactor Applications components included in this distribution include selected computational tools and supporting infrastructure that solve neutronics, thermal-hydraulics, fuel performance, and coupled neutronics-thermal hydraulics problems. The infrastructure components provide a simplified common user input capability and provide for the physics integration with data transfer and coupled-physics iterative solution algorithms. Neutronics analysis can be performed for 2D lattices, 2D core and 3D core problems for pressurized water reactor geometries that can be used to calculate criticality and fission rate distributions by pin for input fuel compositions. MPACT uses the Method of Characteristics transport approach for 2D problems. For 3D problems, MPACT uses the 2D/1D method which uses 2D MOC in a radial plane and diffusion or SPn in the axial direction. MPACT includes integrated cross section capabilities that provide problem-specific cross sections generated using the subgroup methodology. The code can be executed both 2D and 3D problems in parallel to reduce overall run time. A thermal-hydraulics capability is provided with CTF (an updated version of COBRA-TF) that allows thermal-hydraulics analyses for single and multiple assemblies using the simplified VERA common input. This distribution also includes coupled neutronics/thermal-hydraulics capabilities to allow calculations using MPACT coupled with CTF. The VERA fuel rod performance component BISON calculates, on a 2D or 3D basis, fuel rod temperature, fuel rod internal pressure, free gas volume, clad integrity and fuel rod waterside diameter. These capabilities allow simulation of power cycling, fuel conditioning and deconditioning, high burnup performance, power uprate scoping studies, and accident performance. Input/Output capabilities include the VERA Common Input (VERAIn) script which converts the ASCII common input file to the intermediate XML used to drive all of the physics codes in the VERA Core Simulator (VERA-CS). VERA component codes either input the VERA XML format directly, or provide a preprocessor which can convert the XML into native input. VERAView is an interactive graphical interface for the visualization and engineering analyses of output data from VERA. The python-based software is easy to install and intuitive to use, and provides instantaneous 2D and 3D images, 1D plots, and alpha-numeric data from VERA multi-physics simulations. Testing within CASL has focused primarily on Westinghouse four-loop reactor geometries and conditions with example problems included in the distribution.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Structure-Properties Correlation in Si Nanoparticles by Total Scattering and Computer Simulations

High-energy synchrotron X-ray diffraction coupled to atomic pair distribution function analysis and computer simulations is used to determine the atomic-scale structure of silicon (Si) nanoparticles obtained by two different synthetic routes. Results show that Si nanoparticles may have significant structural differences depending on the synthesis route and surface chemistry. In this case, one method produced Si nanoparticles that are highly crystalline but surface oxidized, whereas a different method yields organic ligand-passivated nanoparticles without surface oxide but that are structurally distorted at the atomic scale. Particular structural features of the oxide-free Si nanoparticles such as average first coordination numbers, length of structural coherence, and degree of local distortions are compared to their optical properties such as photoluminescence emission energy, quantum yield, and Raman spectra. A clear structure–properties correlation is observed indicating that the former may need to be taken into account when considering the latter.

Petkov, Valeri↗

SAMICS marketing and distribution model

A SAMICS (Solar Array Manufacturing Industry Costing Standards) was formulated as a computer simulation model. Given a proper description of the manufacturing technology as input, this model computes the manufacturing price of solar arrays for a broad range of production levels. This report presents a model for computing these marketing and distribution costs, the end point of the model being the loading dock of the final manufacturer.

Source record↗

The Martian climate and energy balance models with CO2/H2O atmospheres

The analysis begins with a seasonal energy balance model (EBM) for Mars. This is used to compute surface temperature versus x = sin(latitude) and time over the seasonal cycle. The core model also computes the evolving boundaries of the CO2 icecaps, net sublimational/condensation rates, and the resulting seasonal pressure wave. Model results are compared with surface temperature and pressure history data at Viking lander sites, indicating fairly good agreement when meridional heat transport is represented by a thermal diffusion coefficient D approx. 0.015 W/sq. m/K. Condensational wind distributions are also computed. An analytic model of Martian wind circulation is then proposed, as an extension of the EMB, which incorporates vertical wind profiles containing an x-dependent function evaluated by substitution in the equation defining the diffusion coefficient. This leads to a parameterization of D(x) and of the meridional circulation which recovers the high surface winds predicted by dynamic Mars atmosphere models (approx. 10 m/sec). Peak diffusion coefficients, D approx. 0.6 w/sq m/K, are found over strong Hadley zones - some 40 times larger than those of high-latitude baroclinic eddies. When the wind parameterization is used to find streamline patterns over Martian seasons, the resulting picture shows overturning hemispheric Hadley cells crossing the equator during solstices, and attaining peak intensities during the south summer dust storm season, while condensational winds are most important near the polar caps.

Hoffert, M. I.↗

On the role of topological complexity in spontaneous development of current sheets

The computations presented in this work aim to asses the importance of field line interlacing on spontaneous development of current sheets. From Parker's magnetostatic theorem, such development of current sheets is inevitable in a topologically complex magnetofluid, with infinite electrical conductivity, at equilibrium. Relevant initial value problems are constructed by superposition of two untwisted component fields, each component field being represented by a pair of global magnetic flux surface. The intensity of field line interlacing is then specified by the relative amplitude of the two superposed fields. The computations are performed by varying this relative amplitude. Also to have a direct visualization of current sheet formation, we follow the evolution of flux surfaces instead of the vector magnetic field. An important finding of this paper is in the demonstration that initial field lines having intense interlacing tend to develop current sheets which are distributed throughout the computational domain with no preference for topologically favorable sites like magnetic nulls or field reversal layers. The onsets of these current sheets are attributed to favorable contortions of magnetic flux surfaces where two oppositely directed parts of the same field line or different field lines come to close proximity. However, for less intensely interlaced field lines, the simulations indicate development of current sheets at sites only where the magnetic topology is favorable. These current sheets originate as two sets of anti-parallel complimentary field lines press onto each other.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Brief Summary on Formalizing Parallel Tensor Distributions, Redistributions, and Algorithm Derivations

Large-scale datasets in computational chemistry typically require distributed-memory parallel methods to perform a special operation known as tensor contraction. Tensors are multidimensional arrays, and a tensor contraction is akin to matrix multiplication with special types of permutations. Creating an efficient algorithm and optimized implementation in this domain is complex, tedious, and error-prone. To address this, we develop a notation to express data distributions so that we can apply use automated methods to find optimized implementations for tensor contractions. We consider the spin-adapted coupled cluster singles and doubles method from computational chemistry and use our methodology to produce an efficient implementation. Experiments per- formed on the IBM Blue Gene/Q and Cray XC30 demonstrate impact both improved performance and reduced memory consumption.

97 MATHEMATICS AND COMPUTING↗

NEUTRINO TRANSFER IN THREE DIMENSIONS FOR CORE-COLLAPSE SUPERNOVAE. I. STATIC CONFIGURATIONS

We develop a numerical code to calculate the neutrino transfer with multi-energy and multi-angle in three dimensions (3D) for the study of core-collapse supernovae. The numerical code solves the Boltzmann equations for neutrino distributions by the discrete-ordinate (S{sub n} ) method with a fully implicit differencing for time advance. The Boltzmann equations are formulated in the inertial frame with collision terms being evaluated to the zeroth order of v/c. A basic set of neutrino reactions for three neutrino species is implemented together with a realistic equation of state of dense matter. The pair process is included approximately in order to keep the system linear. We present numerical results for a set of test problems to demonstrate the ability of the code. The numerical treatments of advection and collision terms are validated first in the diffusion and free-streaming limits. Then we compute steady neutrino distributions for a background extracted from a spherically symmetric, general relativistic simulation of a 15 M{sub Sun} star and compare them with the results in the latter computation. We also demonstrate multi-dimensional capabilities of the 3D code solving neutrino transfers for artificially deformed supernova cores in 2D and 3D. Formal solutions along neutrino paths are utilized as exact solutions. We plan to apply this code to the 3D neutrino-radiation hydrodynamics simulations of supernovae. This is the first article in a series of reports on the development.

79 ASTRONOMY AND ASTROPHYSICS↗

A novel framework for hosting capacity analysis with spatio-temporal probabilistic voltage sensitivity analysis

Smart grids are envisioned to accommodate high penetration of distributed photovoltaic (PV) generation, which may cause adverse grid impacts in terms of voltage violations. Therefore, PV Hosting capacity is being used as a planning tool to determine the maximum PV installation capacity that causes the first voltage violation and above which would require infrastructure upgrades. Additionally, traditional methods of Hosting capacity analysis are scenario based and computationally complex as they rely on iterative load flow algorithms that require investigating a large number of scenarios for accurate assessment of PV impacts. Therefore, this paper presents a computationally efficient analytical approach to compute the probability distribution of voltage change due to random behavior of randomly located multiple distributed PVs. The proposed approach is based on Spatio-temporal probabilistic voltage sensitivity analysis that exploits both spatial and temporal uncertainties associated with PV injections. Thereafter, the derived distribution is used to quantify voltage violations for various PV penetration levels and subsequently determine the hosting capacity of the system without the need to examine large number of scenarios. Results of the proposed framework are validated via conventional load flow based simulation approach on the IEEE 37 and IEEE 123 node test systems.

42 ENGINEERING↗

Parallel automated adaptive procedures for unstructured meshes

Consideration is given to the techniques required to support adaptive analysis of automatically generated unstructured meshes on distributed memory MIMD parallel computers. The key areas of new development are focused on the support of effective parallel computations when the structure of the numerical discretization, the mesh, is evolving, and in fact constructed, during the computation. All the procedures presented operate in parallel on already distributed mesh information. Starting from a mesh definition in terms of a topological hierarchy, techniques to support the distribution, redistribution and communication among the mesh entities over the processors is given, and algorithms to dynamically balance processor workload based on the migration of mesh entities are given. A procedure to automatically generate meshes in parallel, starting from CAD geometric models, is given. Parallel procedures to enrich the mesh through local mesh modifications are also given. Finally, the combination of these techniques to produce a parallel automated finite element analysis procedure for rotorcraft aerodynamics calculations is discussed and demonstrated.

Shephard, M. S.↗

Surfer: An Extensible Pull-Based Framework for Resource Selection and Ranking

Grid computing aims to connect large numbers of geographically and organizationally distributed resources to increase computational power; resource utilization, and resource accessibility. In order to effectively utilize grids, users need to be connected to the best available resources at any given time. As grids are in constant flux, users cannot be expected to keep up with the configuration and status of the grid, thus they must be provided with automatic resource brokering for selecting and ranking resources meeting constraints and preferences they specify. This paper presents a new OGSI-compliant resource selection and ranking framework called Surfer that has been implemented as part of NASA's Information Power Grid (IPG) project. Surfer is highly extensible and may be integrated into any grid environment by adding information providers knowledgeable about that environment.

Zolano, Paul Z.↗