Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data parallelism”

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,567 records · Page 87

Parallelized Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. Model inversion of the thermography data requires a fast method for performing the forward problem. Viable numerical methods for the thermal response forward problem are finite element, finite difference and the quadrupole method. Normally both the finite element and finite difference methods solve for the thermal response in the time domain which limits one’s ability to increase the speed of the simulation by parallelization. In contrast, the quadrupole method solves for the Laplace transform of the thermal response. One of the features of the Laplace transform methodology is the solution at any discrete time is independent of the solution at all other times. Therefore, it is easy to separate into a set of independent calculations with each of the times of interest being performed in parallel. Additionally, the numeric inversion of the Laplace transform typically involves numerically solving for the Laplace transform at multiple Laplace frequencies. Each of those solutions are also independent of solutions at other frequencies and can be calculated in parallel. By parallelization of this method, it is possible to perform the simulations of three-dimensional configurations in seconds. When the input stimulus for thermal response is a delta function heat flux (a reasonable approximation for flash heating), the thermal response is smooth. For this case, it is possible to accurately estimate the thermal response at any time within a given time interval from a set of simulations separated by exponentially increasing time steps. From these simulations, it is possible to accurately interpolate to find the response at intermediate times by a spline interpolation of the logarithm of time versus logarithm of temperature. The thermal response with exponential time stepping is shown to produce values for the thermal response which are within 1% of values within the time interval. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are obtained from other inspection methods.

Thermography↗

Parallelized Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. Model inversion of the thermography data requires a fast method for performing the forward problem. Viable numerical methods for the thermal response forward problem are finite element, finite difference and the quadrupole method. Normally both the finite element and finite difference methods solve for the thermal response in the time domain which limits one’s ability to increase the speed of the simulation by parallelization. In contrast, the quadrupole method solves for the Laplace transform of the thermal response. One of the features of the Laplace transform methodology is the solution at any discrete time is independent of the solution at all other times. Therefore, it is easy to separate into a set of independent calculations with each of the times of interest being performed in parallel. Additionally, the numeric inversion of the Laplace transform typically involves numerically solving for the Laplace transform at multiple Laplace frequencies. Each of those solutions are also independent of solutions at other frequencies and can be calculated in parallel. By parallelization of this method, it is possible to perform the simulations of three-dimensional configurations in seconds. When the input stimulus for thermal response is a delta function heat flux (a reasonable approximation for flash heating), the thermal response is smooth. For this case, it is possible to accurately estimate the thermal response at any time within a given time interval from a set of simulations separated by exponentially increasing time steps. From these simulations, it is possible to accurately interpolate to find the response at intermediate times by a spline interpolation of the logarithm of time versus logarithm of temperature. The thermal response with exponential time stepping is shown to produce values for the thermal response which are within 1% of values within the time interval. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are obtained from other inspection methods.

Thermography↗

First Lunar Wake Passage of ARTEMIS: Discrimination of Wake Effects and Solar Wind Fluctuations by 3D Hybrid Simulations

The spacecraft P1 of the new ARTEMIS (Acceleration, Reconnection, Turbulence, and Electrodynamics of the Moon's Interaction with the Sun) mission passed the lunar wake for the first time on February 13, 2010. We present magnetic field and plasma data of this event and results of 3D hybrid simulations. As the solar wind magnetic field was highly dynamic during the passage, a simulation with stationary solar wind input cannot distinguish whether distortions were caused by these solar wind variations or by the lunar wake; therefore, a dynamic real-time simulation of the flyby has been performed. The input values of this simulation are taken from NASA OMNI data and adapted to the P1 data, resulting in a good agreement between simulation and measurements. Combined with the stationary simulation showing non-transient lunar wake structures, a separation of solar wind and wake effects is achieved. An anisotropy in the magnitude of the plasma bulk flow velocity caused by a non-vanishing magnetic field component parallel to the solar wind flow and perturbations created by counterstreaming ions in the lunar wake are observed in data and simulations. The simulations help to interpret the data granting us the opportunity to examine the entire lunar plasma environment and, thus, extending the possibilities of measurements alone: A comparison of a simulation cross section to theoretical predictions of MHD wave propagation shows that all three basic MHD modes are present in the lunar wake and that their expansion governs the lunar wake refilling process.

Wiehle, S.↗

Estimating Lossy Compressibility of Scientific Data Using Deep Neural Networks

Simulation based scientific applications generate increasingly large amounts of data on high-performance computing (HPC) systems. To allow data to be stored and analyzed efficiently, data compression is often utilized to reduce the volume and velocity of data. However, a question often raised by domain scientists is the level of compression that can be expected so that they can make more informed decisions, balancing between accuracy and performance. In this letter, we propose a deep neural network based approach for estimating the compressibility of scientific data. To train the neural network, we build both general features as well as compressor-specific features so that the characteristics of both data and lossy compressors are captured in training. Our approach is demonstrated to outperform a prior analytical model as well as a sampling based approach in the case of a biased estimation, i.e., for SZ. However, for the unbiased estimation (i.e., ZFP), the sampling based approach yields the best accuracy, despite the high overhead involved in sampling the target dataset.

97 MATHEMATICS AND COMPUTING↗

GPU Acceleration of Large-Scale Full-Frequency GW Calculations

Many-body perturbation theory is a powerful method to simulate electronic excitations in molecules and materials starting from the output of density functional theory calculations. By implementing the theory efficiently so as to run at scale on the latest leadership high-performance computing systems it is possible to extend the scope of GW calculations. Here, we present a GPU acceleration study of the full-frequency GW method as implemented in the WEST code. Excellent performance is achieved through the use of (i) optimized GPU libraries, e.g., cuFFT and cuBLAS, (ii) a hierarchical parallelization strategy that minimizes CPU-CPU, CPU-GPU, and GPU-GPU data transfer operations, (iii) nonblocking MPI communications that overlap with GPU computations, and (iv) mixed precision in selected portions of the code. A series of performance benchmarks has been carried out on leadership high-performance computing systems, showing a substantial speedup of the GPU-accelerated version of WEST with respect to its CPU version. Good strong and weak scaling is demonstrated using up to 25 920 GPUs. Finally, we showcase the capability of the GPU version of WEST for large-scale, full-frequency GW calculations of realistic systems, e.g., a nanostructure, an interface, and a defect, comprising up to 10 368 valence electrons.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A planar-sheet nongraphitic zero-bandgap sp 2 carbon phase made by the low-temperature reaction of γ-graphyne

The highest sheet symmetry form of graphyne, with one triple bond between each neighboring hexagon in graphene, irreversibly transforms exothermically at ambient pressure and low temperatures into a nongraphitic, planar-sheet, zero-bandgap phase consisting of intrasheet-bonded sp 2 carbons. The synthesis of this sp 2 carbon phase is demonstrated, and other carbon phases are described for possible future synthesis from graphyne without breaking graphyne bonds. While measurements and theory indicate that the reacting graphyne becomes nonplanar because of sheet wrinkling produced by dimensional mismatch between reacted and nonreacted sheet regions, sheet planarity is regained when the reaction is complete. Although the observed elimination of triple bonds to make parallel planar sp 2 carbon sheets likely requires ordered transformation within each sheet, diffraction data for reacted multisheet stacks indicate that the relative lateral positions of neighboring sheets are disordered, as predicted, since no crystalline diffraction peak (other than for the intersheet spacing) is observed.

Aliev, Ali E. (ORCID:0000000300349707)↗

Phyllosphere exudates select for distinct microbiome members in sorghum epicuticular wax and aerial root mucilage

Phyllosphere exudates create specialized microhabitats that shape microbial community diversity. We explored the microbiome associated with two sorghum phyllosphere exudates, the epicuticular wax and aerial root mucilage. We assessed the microbiome associated with the wax from sorghum plants over two growth stages, and the root mucilage additionally from nitrogen-fertilized and non-fertilized plants. In parallel, we isolated and characterized hundreds of bacteria from wax and mucilage, and integrated data from cultivation-independent and cultivation-dependent approaches to gain insights into exudate diversity and bacterial phenotypes. We found that Sphingomonadaceae and Rhizobiaceae families were the major taxa in the wax regardless of water availability and plant developmental stage to plants. The cultivation-independent mucilage-associated bacterial microbiome contained Erwiniaceae, Flavobacteriaceae, Rhizobiaceae, Pseudomonadaceae, Sphingomonadaceae, and its structure was strongly influenced by sorghum development but only modestly influenced by fertilization. In contrast, the fungal community structure of mucilage was strongly affected by the year of sampling but not by fertilization or plant developmental stage, suggesting a decoupling of fungal-bacterial dynamics in the mucilage. Our bacterial isolate collection from wax and mucilage had several isolates that matched 100% to detected amplicon sequence variants, and were enriched on media that selected for phenotypes including phosphate solubilization, putative diazotrophy, resistance to desiccation, capability to grow on methanol as a carbon source, and ability to grow in the presence of linalool and β-caryophyllene (terpenes in sorghum wax). This work expands our understanding of the microbiome of phyllosphere exudates and supports our long-term goal to translate microbiome research to support sorghum cultivation.

59 BASIC BIOLOGICAL SCIENCES↗

AMReX and pyAMReX: Looking beyond the exascale computing project

AMReX is a software framework for the development of block-structured mesh applications with adaptive mesh refinement (AMR). AMReX was initially developed and supported by the AMReX Co-Design Center as part of the U.S. DOE Exascale Computing Project (ECP), and is continuing to grow post-ECP. In addition to adding new functionality and performance improvements to the core AMReX framework, we have also developed a Python binding, pyAMReX, that provides a bridge between AMReX-based application codes and the data science ecosystem. pyAMReX provides zero-copy application GPU data access for AI/ML, in situ analysis and application coupling, and enables rapid, massively parallel prototyping. In this paper we review the overall functionality of AMReX and pyAMReX, focusing on new developments, new functionality, and optimizations of key operations. We also summarize capabilities of ECP projects that used AMReX and provide an overview of new, non-ECP applications.

Myers, Andrew↗

A view into FleCSI [Slides]

FleCSI is largely based around sharing of memory spaces. FleCSI facilitates parallel execution and sharing of memory across memory spaces. Execute() transforms sender-side data structures into receiver-side data structures.

97 MATHEMATICS AND COMPUTING↗

Interagency Agreement No. DE-SC0006988 (Final Technical Report)

The overarching goal of this project was to advance understanding of deep convection updraft microphysics—the source of long-lived stratiform ice—by combining detailed mining of in situ and remote sensing data from the MC3E field campaign with detailed 3D simulations. The project began with parallel work, first strictly on the remote-sensing observation side via dedicated analysis of polarimetric radar signatures, which is a relatively new area, on the one hand. On the other hand, a more traditional but well-formulated preparation and comparison of detailed aerosol-aware simulations with in situ measurements was prepared. In the final phase of this work, these parallel elements were brought together into an integrated analysis of in situ and remote-sensing observations with model results. The project also supported team member involvement in collaborative activities.

54 ENVIRONMENTAL SCIENCES↗

Progress Report on Performance of A709 and G91 Steels in Sodium

Specimens in six different processing and heat treatment conditions of A709 H58776 were exposed to sodium at 550, 600, and 650°C, respectively for various exposure times. G91 base metal was tested in sodium at 550, 600, and 650°C, and G91 weldment at 550 and 600°C, respectively. Parallel thermal aging experiments were conducted on G91 and A709 to obtain thermal aging data for comparison with sodium exposure results to separate the thermal and sodium effects. The corrosion data obtained on A709 steel continue to show low corrosion rates, which indicate good compatibility of A709 with sodium when oxygen content is controlled. G91 has also shown acceptable corrosion rates over the temperature range of 550-650°C in sodium. Thermal aging or sodium exposures reduced the tensile strength, uniform elongation, and total elongation for A709 ESR, AOD specimens. These specimens consistently showed higher yield stress and tensile strength and lower uniform and total elongations after sodium exposure than after thermal aging under comparable testing conditions. It implies that there was an additional effect of sodium exposure that gave rise to an increase in tensile strength of A709. Dynamic strain aginginduced flow serrations were completely removed by either thermal aging or sodium exposure at 650°C, which imply that carbon or nitrogen in solution was removed from the solution and formed precipitates during aging or sodium exposures. The tensile data ruled out the possibility of decarburization under the sodium exposure conditions. The HOMO specimens behaved somewhat differently from the ESR and AOD specimens in sodium. While thermal aging and sodium exposures at 550 and 600°C decreased the yield stress and the ultimate tensile strength of G91, there was virtually no additional effect resulting from sodium exposures at these two temperatures. In contrast, sodium exposures at 650°C had a drastic effect on the yield stress and the ultimate tensile strength of G91. It is suggested that G91 experienced decarburization in the 650°C sodium environments. Carbon concentrations in sodium in the SMT-1 and SMT-2 loops were determined by a foil equilibration method. The estimated carbon concentration was in the range of 0.8-1.2 ppm in the SMT-1 loop and 0.3-0.7 ppm in the SMT-2 loop. The carbon activity in sodium in the SMT-2 loop was estimated to be 0.03-0.08 at 600°C, and 0.08-0.2 at 550°C. The carbon activity in sodium at 650°C in the SMT-1 loop was 0.04-0.07. Equilibrium simulation of the carburization – decarburization processes was conducted for A709 and G91 steels exposed in sodium environments at temperatures of 550-700°C. The carbon activity-concentration relationship for G91 was re-evaluated by considering four phases in G91, i.e. bcc ferrite, M 23 C 6 , NbC and VC carbides. It was found that the carburization-decarburization process in G91 steel was dictated by M 23 C 6 carbides at high carbon activities, while NbC carbides dominated the process at low carbon activities. Formation of NbC increases the decarburization resistance of G91 steel. It remains to be understood whether the decarburization resistance provided by MC carbides in G91 can be maintained during long-term operations of SFRs. The carbon concentration-activity relationship for A709 was calculated based on the equilibrium of fcc-austenite and M 23 C 6 carbide phase. The calculations showed that A709 decarburizes at 650°C in the SMT-1 loop environments, which is different from the experimental findings. The carbon concentration-activity relationship in A709 was further evaluated by considering four phases, including fcc-austenite, M 23 C 6 , NbC, and TiC carbide phases. TiC is the most stable carbides among the three carbide phases. The carburization-decarburization process in A709 is dictated by M 23 C 6 carbides at high carbon activities and by TiC at low carbon activities. Formation of TiC increases the decarburization resistance of A709 steel. Because of the complexity of the precipitation process in A709 during sodium exposures, detailed characterization of precipitates in sodium-exposed specimens and improved thermodynamic models are needed to understand the carburization-decarburization behavior of A709. Kinetic analysis of carbon transfer will be investigated for G91 and A709 in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

2011.5 Revision of the Evaluated Nuclear Data Library: ENDL2011.5 (legacy and GNDS), and ENDL2011.5-direct

LLNL’s Nuclear Data and Theory Group has created the 2011.5 revised release of the Evaluated Nuclear Data Library (ENDL2011.5). ENDL2011.5 is designed to support LLNL’s current and future nuclear data needs and will be employed in nuclear reactor, nuclear security and stockpile stewardship simulations with ASC codes. This database is currently the most complete nuclear database for Monte Carlo and deterministic transport of neutrons and charged particles. This library was assembled with strong support from the ASC PEM and Attribution programs, leveraged with support from Campaign 4 and the DOE/Office of Science’s US Nuclear Data Program. This document lists the revisions made in ENDL2011.5 compared with the data existing in the original ENDL2011.2 and ENDL2011.3 releases. These changes are made in parallel with some similar revisions for ENDL2009.5. We also describe ENDL2011.5-direct that directly uses ENDF sources wherever possible.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

PyFLEXTRKR: a flexible feature tracking Python software for convective cloud analysis

Abstract. This paper describes the new open-source framework PyFLEXTRKR (Python FLEXible object TRacKeR), a flexible atmospheric feature tracking software package with specific capabilities to track convective clouds from a variety of observations and model simulations. This software can track any atmospheric 2D objects and handle merging and splitting explicitly. The package has a collection of multi-object identification algorithms, scalable parallelization options, and has been optimized for large datasets including global high-resolution data. We demonstrate applications of PyFLEXTRKR on tracking individual deep convective cells and mesoscale convective systems from observations and model simulations ranging from large-eddy resolving (∼100s m) to mesoscale (∼10s km) resolutions. Visualization, post-processing, and statistical analysis tools are included in the package. New Lagrangian analyses of convective clouds produced by PyFLEXTRKR applicable to a wide range of datasets and scales facilitate advanced model evaluation and development efforts as well as scientific discovery.

54 ENVIRONMENTAL SCIENCES↗

Heat transfer characteristics of porous media

An investigation was conducted regarding the relative effects of conduction and convection in a saturated porous medium. A method reported by Singh et al. (1973) is used to determine the effective thermal conductivity of the saturated porous material. Heat transfer measurements are conducted under conditions of forced convection of the saturated liquid parallel and countercurrent to the flow of heat. The results are compared with the data obtained with the aid of an analytical model.

Singh, B. S.↗

Application of visible linear array technology to earth observation sensors

The present paper identifies the systems engineering aspects of applying solid-state technology to earth observations applications being traditionally performed by point (or multiple-point) detector line scanned mechanisms. It is shown that the translation from a basically serial data flow point-detector mechanically-scanned sensor to a solid state highly parallel linear-array pushbroom sensor results in minimizing mechanical complexity and maximizing electronics complexity, with increased demands upon optical performance in some applications. Technical aspects relevant to highly parallel photodiode linear-array pushbroom applications are discussed. Examples of systems engineering applications are provided.

Noll, R. E.↗

Infrared Imagery of Shuttle (IRIS). Task 2, summary report

End-to-end tests of a 16 element indium antimonide sensor array and 10 channels of associated electronic signal processing were completed. Quantitative data were gathered on system responsivity, frequency response, noise, stray capacitance effects, and sensor paralleling. These tests verify that the temperature accuracies, predicted in the Task 1 study, can be obtained with a very carefully designed electro-optical flight system. Pre-flight and inflight calibration of a high quality are mandatory to obtain these accuracies. Also, optical crosstalk in the array-dewar assembly must be carefully eliminated by its design. Tests of the scaled up tracking system reticle also demonstrate that the predicted tracking system accuracies can be met in the flight system. In addition, improvements in the reticle pattern and electronics are possible, which will reduce the complexity of the flight system and increase tracking accuracy.

Chocol, C. J.↗

Measured effect of wind generation on the fuel consumption of an isolated diesel power system

The Block Island Power Company (BIPCO), on Block Island, Rhode Island, operates an isolated electric power system consisting of diesel generation and an experimental wind turbine. The 150-kW wind turbine, designated MOD-OA by the U.S. Department of Energy is typically operated in parallel with two diesel generators to serve an average winter load of 350 kW. Wind generation serves up to 60 percent of the system demand depending on wind speed and total system load. Results of diesel fuel consumption measurements are given for the diesel units operated in parallel with the wind turbine and again without the wind turbine. The fuel consumption data are used to calculate the amount of fuel displaced by wind energy. Results indicate that the wind turbine displaced 25,700 lbs. of the diesel fuel during the test period, representing a calculated reduction in fuel consumption of 6.7 percent while generating 11 percent of the total electric energy. The amount of displaced fuel depends on operating conditions and system load. It is also shown that diesel engine throttle activity resulting from wind gusts which rapidly change the wind turbine output do not significantly influence fuel consumption.

Stiller, P. H.↗

Efficient implementation of real-time programs under the VAX/VMS operating system

Techniques for writing efficient real-time programs under the VAX/VMS oprating system are presented. Basic operations are presented for executing at real-time priority and for avoiding needlless processing delays. A highly efficient technique for accessing physical devices by mapping to the input/output space and accessing the device registrs directly is described. To illustrate the application of the technique, examples are included of different uses of the technique on three devices in the Langley Avionics Integration Research Lab (AIRLAB): the KW11-K dual programmable real-time clock, the Parallel Communications Link (PCL11-B) communication system, and the Datacom Synchronization Network. Timing data are included to demonstrate the performance improvements realized with these applications of the technique.

Johnson, S. C.↗