Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “code optimization”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

GPU accelerated flow solver for direct numerical simulation of turbulent flows

Graphical processing units (GPUs), characterized by significant computing performance, are nowadays very appealing for the solution of computationally demanding tasks in a wide variety of scientific applications. However, to run on GPUs, existing codes need to be ported and optimized, a procedure which is not yet standardized and may require non trivial efforts, even to high-performance computing specialists. In the present paper we accurately describe the porting to CUDA (Compute Unified Device Architecture) of a finite-difference compressible Navier–Stokes solver, suitable for direct numerical simulation (DNS) of turbulent flows. Porting and validation processes are illustrated in detail, with emphasis on computational strategies and techniques that can be applied to overcome typical bottlenecks arising from the porting of common computational fluid dynamics solvers. We demonstrate that a careful optimization work is crucial to get the highest performance from GPU accelerators. The results show that the overall speedup of one NVIDIA Tesla S2070 GPU is approximately 22 compared with one AMD Opteron 2352 Barcelona chip and 11 compared with one Intel Xeon X5650 Westmere core. The potential of GPU devices in the simulation of unsteady three-dimensional turbulent flows is proved by performing a DNS of a spatially evolving compressible mixing layer.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Electron-optical systems for planar gyrotrons

The methodology of designing an electron-optical system (EOS) that forms sheet helical electron beams (HEBs) for high-power gyrotrons is developed. As an example, we consider the EOS for a 140-GHz gyrotron operated at the first harmonic of the cyclotron frequency with an accelerating voltage of 50 kV, a beam current of 30 A, and a magnetic field compression of 36. A planar geometry of the magnetron-injection gun (MIG) is suggested. The adiabatic theory of MIGs modified for the planar geometry of EOS is used for preliminary estimations of MIG parameters. Numerical simulation of the HEB properties based on the CST STUDIO SUITE 3D code is performed to find the optimal configuration of a planar MIG. The accuracy of the calculated data is discussed. The main factors that affect the HEB quality are considered. It is shown that a sheet HEB with a pitch-factor of 1.3 and velocity spread not exceeding 25%–30% can be formed; this is quite acceptable for high-efficiency operation of modern gyrotrons. Calculation of the beam-wave interaction with the obtained HEB parameters proved that a high output power with a sufficiently good efficiency of about 20% can be reached. Simulations show the feasibility of the experimental implementation of a novel planar EOS and its use in short-wave planar gyrotrons. The developed technique can be used for the study and optimization of planar gyrotrons of different frequency bands and power levels.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Utilization of Hydroturbines in Wastewater Treatment Plants

Wastewater treatment plants (WWTPs) are a significant energy consumer, yet there are several opportunities for implementing on-site power generation systems. Within the treatment process, the high flow rate of effluent is produced and discharged to a nearby water body by gravity. Thus, hydroturbines can be utilized to generate power in such an application due to a difference in elevation and high flow rate. This paper presents a case study of introducing a hydroturbine in a WWTP in Wisconsin and evaluating the power output in addition to determining the energy savings. The WWTP considered in this study has an effluent flow rate of 190 MGD (million gallons per day) and elevation difference of 3 m (10 ft) between the final stage of treatment and the discharge point. Based on the parameters above; hubless rim-drive Kaplan type hydroturbine (RDT) is the optimal choice to be used in such an application. The RDT is designed and optimized by using in-house code. A computational fluid dynamics (CFD) software is applied to evaluate the performance of the proposed model, and the system is simulated through homer software to validate the results generated by the CFD. The expected savings is estimated to be 1564 MWh/yr.

Energy & Fuels↗
DOE OSTI2016

Measurement of the production cross section for charmed baryons in proton–nucleus interactions at 70 GeV

The results of an analysis of data from the SERP-E-184 experiment devoted to studying mechanisms of the production of charmed particles in proton–nucleus interactions at 70 GeV and their decays are presented. The data in question were obtained upon irradiating the SVD-2 active target consisting of carbon, silicon, and lead plates with a beam of 70-GeV protons. A detailed simulation on the basis of the FRITIOF7.02 and GEANT3.21 code packages made it possible to optimize event-selection criteria and to calculate the detection efficiency for Λ{sub c}{sup +} baryons. After selecting a signal from the threebody decay of a Λ{sub c}{sup +} baryon, the inclusive cross section for its production at near-threshold energies, its lifetime, and the parameter of the A dependence of the cross section were found. The Λ{sub c}{sup +} -baryon yields are tabulated along with data from other experiments and theoretical predictions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2014

Trace-Driven Memory Access Pattern Recognition in Computational Kernels. In: WOSC '14 Proceedings of the Second Workshop on Optimizing Stencil Computations

Classifying memory access patterns is paramount to the selection of the right set of optimizations and determination of the parallelization strategy. Static analyses suffer from ambiguities present in source code, which modern compilation techniques, such as profile-guided optimization, alleviate by observing runtime behavior and feeding back into the compilation flow. This paper discusses a dynamic analysis technique for recognizing memory access patterns, with application to the stencils domain, and presents our design and C++ implementation using the memory-tracing tool Gleipnir. Finally, we evaluate and discuss the performance and matching capability of our classifiers in the context of the Polybench scientific benchmark suite, which includes both stencil and matrix computations.

Park, Eunjung↗
DOE OSTI2022

SODA Synthesizer: an Open-source, Multi-level, Modular, Extensible Compiler from High-level Frameworks to Silicon

The SODA Synthesizer is an open-source modular, end-to-end hardware compiler framework. The SODA frontend, developed in MLIR, performs system-level design, code partitioning, and high-level optimizations to prepare the specifications for the hardware synthesis. The backend is based on a state-of-the-art high-level synthesis tool, and generates the final hardware design. The backend can interface with logic synthesis tools for field programmable gate arrays or with commercial and open-source logic synthesis tools for application-specific integrated circuits. We discuss the opportunities and challenges in integrating with commercial and open-source tools both at the frontend and backend, and the unique opportunities that an open-source hardware design ecosystem provides.

Bohm Agostini, Nicolas↗
DOE OSTI2024

High-Throughput Computing: Case Study of Medical Image Processing Applications

HPC is designed for large-scale simulations using monolithic codes of tightly coupled processes highly optimized to deliver decreased time to solution. Medical image processing is not a traditional field of HPC. Similar to AI applications, medical image processing parses large datasets, typically multiple times, to support a variety of studies for classification, diagnosis or monitoring purposes. The convergence of AI, HPC and Big Data encouraged more fields using image processing to transition to HPC. However, not all applications benefit from the same optimizations. In this paper we focus on high throughput medical image processing applications that analyze a huge dataset of small MRI images and that require HPC systems to decrease the time of parsing the entire dataset and not individual MRIs. We show in this research the performance of running SLANT, an image processing application for a whole brain segmentation, on large-scale systems and highlight performance limitations. We present optimizations prioritizing throughput that exhibit a 3.5x speed-up on the Summit Supercomputer that can be used as a baseline for building a high-throughput execution framework for other HPC systems.

Predescu, Maria↗
DOE OSTI2026

torch-einshard v1.0

torch-einshard is a Python library for describing local and distributed PyTorch tensor computations with compact, einsum-like notation. Its expressions name logical axes, specify how they are sharded across a PyTorch DeviceMesh, and represent partial reductions. The library automatically performs contractions, permutations, reshaping, splitting, gathering, reduction, reduce-scatter, and repartitioning while preserving autograd. Additional features include sharding-aware FFTs, tensor rolls, halo exchange, sliding windows, 1D–3D convolutions, uneven-shard handling, parameter initialization and gradient management, and cost-based execution planning. It is designed for scientific machine learning and large-model workloads, including tensor-, sequence-, and spatial-parallel MLPs, attention, convolutions, and spectral operations. Compared with manually combining torch.einsum and distributed collectives, torch-einshard expresses both the mathematical operation and data placement in one readable formula. This reduces boilerplate and synchronization errors, keeps forward and backward communication consistent, and allows the library to select optimized collective strategies without changing model code.

Morozov, Dmitriy [Lawrence Berkeley National Labor↗
DOE OSTI2011

Monte Carlo Shielding Analysis Capabilities with MAVRIC

Monte Carlo shielding analysis capabilities in SCALE 6 are centered on the CADIS methodology Consistent Adjoint Driven Importance Sampling. CADIS is used to create an importance map for space/energy weight windows as well as a biased source distribution. New to SCALE 6 are the Monaco functional module, a multi-group fixed-source Monte Carlo transport code, and the MAVRIC sequence (Monaco with Automated Variance Reduction Using Importance Calculations). MAVRIC uses the Denovo code (also new to SCALE 6) to compute coarse-mesh discrete ordinates solutions which are used by CADIS to form an importance map and biased source distribution for the Monaco Monte Carlo code. MAVRIC allows the user to optimize the Monaco calculation for a specify tally using the CADIS method with little extra input compared to a standard Monte Carlo calculation. When computing several tallies at once or a mesh tally over a large volume of space, an extension of the CADIS method called FW-CADIS can be used to help the Monte Carlo simulation spread particles over phase space to get more uniform relative uncertainties.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2022

Optimization of the deep neural network parameters for generating homogenized fuel assembly data for nodal codes

Homogenized fuel assembly (FA) data is a typical input data for nodal codes. Generating that data, however, could be time-consuming. One of promising ways to mitigate the computational burden of generating macroscopic cross-sections is to use trained artificial neural network (ANN) models for predicting nuclear data. However, there is a challenge to make the model support variable FA geometry. In this work, two most common types of FA were combined in one ANN model. Since there could be multiple ways of converting 2-dimensional FA data into 1-dimensional input vector for ANN, three different approaches of data flattening were evaluated. The input parameters included each fuel pin enrichment, fuel temperature, moderator temperature and boron concentration. The output parameters were 2-group macroscopic cross-sections (XS) and pin power distribution (HFF). A fully connected deep neural network (DNN) model was trained and tested using pre-generated data obtained with lattice physics code STREAM. The results of this study showed no statistically significant difference in the accuracy of XS and HFF generation for all 3 tested input vector orders. This means that fully connected DNN for XS generation demonstrated input sequence invariance. Results of comparing predicted XS data with reference solutions were found sufficiently close considering the reduction of computation time offered by ANN. Mean relative difference (MRD) for all output XS parameters was found below 0.7%, while HFF MRD was found higher compared to XS values, in some cases slightly exceeding 1%, mostly near guide tube locations. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

LDRD Final Report: Advanced Hohlraum Concepts

Indirect drive inertial confinement fusion (ICF) experiments to date have mostly used cylindrical, laser-heated, gas-filled hohlraums to produce the radiation drive needed to symmetrically implode DT-filled fusion capsules. These hohlraums have generally been unable to produce a symmetric radiation drive through the end of the desired drive pulse, and are plagued with complications due to laser-plasma interactions (LPI) that have made it difficult to predict their performance. In this project we developed several alternate hohlraum concepts. These new hohlraums utilize different hohlraum geometries, radiation shields, and foam materials in an attempt to improve performance relative to cylindrical hohlraums. Each alternate design was optimized using radiation hydrodynamic (RH) design codes to implode a reference DT capsule with a high-density carbon (HDC) ablator. The laser power and energy required to produce the desired time-dependent radiation drive, and the resulting time-dependent radiation symmetry for each new concept were compared to the results for a reference cylindrical hohlraum. Since several of the new designs needed extra laser entrance holes (LEHs), techniques to keep small LEHs open longer, including high-Z foam liners and low-Z wires at the LEH axis, were investigated numerically. Supporting experiments and target fabrication efforts were also done as part of this project. On the Janus laser facility plastic tubes open at one end (halfraums) and filled with SiO 2 or Ta 2 O 5 foam were heated with a single 2w laser. Laser propagation and backscatter were measured. Generally the measured propagation was slower than calculated, and the measured laser backscatter was less than calculated. A comparable, scaled up experiment was designed for the NIF facility and four targets were built. Since low density gold foam was identified as a desirable material for lining the LEH and the hohlraum wall, a technique was developed to produce 550 mg/cc gold foam, and a sample of this material was successfully manufactured.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Developing a Recovery Column for Mo-99 from Highly Concentrated Uranyl Nitrate Solutions

The feasibility of a highly concentrated uranyl nitrate target for 99 Mo production and recovery was investigated with respect to three separate resins. Two resins were inorganic-based – titania and zirconia – and the other was organic-based utilizing Aliquot® 336 extraction chromatography. Results show that pH has a significant effect on select fission products such as Zr whereas Mo is reasonably well retained between pH 0 and 2. Temperature also appreciably affected the retention of Mo. Using a normal-phase titania resin (Sachtapore®), the sorption of Mo was investigated by varying acid and uranium concentrations. The stability of the resin was also studied as a function of acid concentration. Results from the batch studies were then input into the in-house Versatile Reaction Separation (VERSE) code in order to design and optimize a recovery column for 99 Mo. The experimental and model results show that >99.9% of the 99 Mo can be recovered in the presence of saturated uranyl nitrate solutions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Streamlining the Coupling of BISON and Dakota Through the NEAMS Workbench

Metallic nuclear fuels for use in advanced reactors are an active area of research and development. Robust, accurate metallic fuel performance models are necessary for the design, analysis, and licensing of such reactors. However, metallic fuel performance models require additional development; they are not as mature as uranium dioxide fuel performance models. To support further metallic fuel development, Oak Ridge National Laboratory and the University of Florida have streamlined the coupling of the BISON fuel performance code with Design Analysis Kit for Optimization and Terascale Applications (Dakota) statistical analysis tool through the Nuclear Energy Advanced Modeling and Simulation (NEAMS) Workbench. This work included performing three different sensitivity analyses on metallic nuclear fuel models in BISON. The analyses examined were a general model of the IFR-1 experiment, the X430 experiment T654 pin, and the X430 experiment T651 pin. The results suggest that BISON and Dakota can be integrated through NEAMS Workbench to perform sensitivity and uncertainty analyses and visualize the results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

H-AMR: A New GPU-accelerated GRMHD Code for Exascale Computing with 3D Adaptive Mesh Refinement and Local Adaptive Time Stepping

General relativistic magnetohydrodynamic (GRMHD) simulations have revolutionized our understanding of black hole accretion. Here, we present a GPU-accelerated GRMHD code H-AMR with multifaceted optimizations that, collectively, accelerate computation by 2–5 orders of magnitude for a wide range of applications. First, it introduces a spherical grid with 3D adaptive mesh refinement that operates in each of the three dimensions independently. This allows us to circumvent the Courant condition near the polar singularity, which otherwise cripples high-resolution computational performance. Second, we demonstrate that local adaptive time stepping on a logarithmic spherical-polar grid accelerates computation by a factor of ≲10 compared to traditional hierarchical time-stepping approaches. Jointly, these unique features lead to an effective speed of ~10 9 zone cycles per second per node on 5400 NVIDIA V100 GPUs (i.e., 900 nodes of the OLCF Summit supercomputer). We illustrate H-AMR's computational performance by presenting the first GRMHD simulation of a tilted thin accretion disk threaded by a toroidal magnetic field around a rapidly spinning black hole. With an effective resolution of 13,440 × 4608 × 8092 cells and a total of ≲22 billion cells and ~0.65 × 10 8 time steps, it is among the largest astrophysical simulations ever performed. We find that frame dragging by the black hole tears up the disk into two independently precessing subdisks. The innermost subdisk rotation axis intermittently aligns with the black hole spin, demonstrating for the first time that such long-sought alignment is possible in the absence of large-scale poloidal magnetic fields.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

A Spectral Approach for the Design of Experiments: Design, Analysis and Algorithms

This paper proposes a new approach to construct high quality space-filling sample designs. First, we propose a novel technique to quantify the space-filling property and optimally trade-off uniformity and randomness in sample designs in arbitrary dimensions. Second, we connect the proposed metric (defined in the spatial domain) to the quality metric of the design performance (defined in the spectral domain). This connection serves as an analytic framework for evaluating the qualitative properties of space-filling designs in general. Using the theoretical insights provided by this spatial-spectral analysis, we derive the notion of optimal space-filling designs, which we refer to as space-filling spectral designs. Third, we propose an efficient estimator to evaluate the space-filling properties of sample designs in arbitrary dimensions and use it to develop an optimization framework for generating high quality space-filling designs. Finally, we carry out a detailed performance comparison on two different applications in varying dimensions: a) image reconstruction and b) surrogate modeling for several benchmark optimization functions and a physics simulation code for inertial confinement fusion (ICF). Our results clearly evidence the superiority of the proposed space-filling designs over existing approaches, particularly in high dimensions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Pre-Commercial Demonstration of Direct Non-Iterative State Estimator (DNSE+)

The objective of this DOE co-funded project is to demonstrate functionality and performance of a production-grade Direct Non-iterative State Estimator (DNSE) through its integration with NYPA’s Energy Management System (EMS) and an enhanced Real Time Dynamics Monitoring System (RTDMS) synchrophasor platform from Electric Power Group (EPG). A production-grade DNSE aims to overcome a major obstacle to the operational use of synchrophasor data in Synchro-Phasor Management Systems (SPMS) by providing SPMS applications with a consistent and a complete synchrophasor data foundation in the same way that a traditional EMS State Estimator (SE) provides to EMS applications. Specifically, the DNSE uses synchrophasor measurements and Supervisory Control and Data Acquisition (SCADA) measurements, together with the complete power system model, to simultaneously obtain the complete state of the entire utility’s operating model at rates that are close to the synchrophasor data rates. The full system coverage property is what distinguishes the DNSE from the linear synchrophasor SEs (which cover only parts of the operating model that are visible through synchrophasor measurements), and DNSE differs from traditional EMS SE in that the solution is obtained without iterations, without initial values (hence the “direct” in the name), and at a much higher execution rate. This presentation provides an overview of the project objectives, a brief description of the main tasks including the development work to productize and optimize DNSE from its current MATLAB code base, RTDMS enhancement for interfacing with the DNSE, factory acceptance testing, field installation, testing and demonstration, as well as the current status of the project execution.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2015

Wilson Dslash Kernel From Lattice QCD Optimization

Lattice Quantum Chromodynamics (LQCD) is a numerical technique used for calculations in Theoretical Nuclear and High Energy Physics. LQCD is traditionally one of the first applications ported to many new high performance computing architectures and indeed LQCD practitioners have been known to design and build custom LQCD computers. Lattice QCD kernels are frequently used as benchmarks (e.g. 168.wupwise in the SPEC suite) and are generally well understood, and as such are ideal to illustrate several optimization techniques. In this chapter we will detail our work in optimizing the Wilson-Dslash kernels for Intel Xeon Phi, however, as we will show the technique gives excellent performance on regular Xeon Architecture as well.

Lattice QCD, Wilson-Dslash, Parallel Computing, Op↗