Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computer Graphics”

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 685 records · Page 38

QMCPACK: Advances in the development, efficiency, and application of auxiliary field and real-space variational and diffusion quantum Monte Carlo

We review recent advances in the capabilities of the open source ab initio Quantum Monte Carlo (QMC) package QMCPACK and the workflow tool Nexus used for greater efficiency and reproducibility. The auxiliary field QMC (AFQMC) implementation has been greatly expanded to include k-point symmetries, tensor-hypercontraction, and accelerated graphical processing unit (GPU) support. These scaling and memory reductions greatly increase the number of orbitals that can practically be included in AFQMC calculations, increasing the accuracy. Advances in real space methods include techniques for accurate computation of bandgaps and for systematically improving the nodal surface of ground state wavefunctions. Results of these calculations can be used to validate application of more approximate electronic structure methods, including GW and density functional based techniques. To provide an improved foundation for these calculations, we utilize a new set of correlation-consistent effective core potentials (pseudopotentials) that are more accurate than previous sets; these can also be applied in quantum-chemical and other many-body applications, not only QMC. These advances increase the efficiency, accuracy, and range of properties that can be studied in both molecules and materials with QMC and QMCPACK.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Hybridized Discontinuous Galerkin Methods for Computational Fluid Dynamics

Hybridizable Discontinuous Galerkin (HDG) methods hold promise for any applications with significant advection character, including thermal hydraulics in light water reactors and advanced reactor concepts and fluid models of plasmas in magnetic confinement fusion. Its features include natural upwinding, local element conservation, and extensibility to arbitrarily high order accuracy. In the last fiscal year we have implemented HDG in the Multiphysics Object-Oriented Simulation Environment (MOOSE). We developed a first-of-its-kind automatic static condensation system in MOOSE’s underlying finite element library libMesh which can condense out arbitrarily many internal variables. Finally, we developed the first preconditioner for HDG discretizations of the Navier-Stokes equations which shows robust performance across a wide range of problem sizes and Reynolds numbers. This preconditioner yields solution times that are equivalent to the fastest developed for industry standard finite volume methods. Moreover, the arbitrarily high-order nature of HDG makes it a prime candidate for acceleration via graphical processing units (GPUs). We believe these developments will hold significant importance for future DOE Nuclear Energy (NE) and Fusion Energy Science (FES) programs.

97 MATHEMATICS AND COMPUTING↗

pyRMG: A framework for high-throughput, large-cell DFT calculations on supercomputers

Exascale computing delivers the raw power to simulate ever larger and more chemically realistic systems, but realizing this potential requires codes that can efficiently use thousands of processors. Our real-space multigrid (RMG) density functional theory (DFT) code’s grid-decomposition approach scales nearly linearly with the number of graphics processing units (GPUs), even for simulations exceeding thousands of atoms. This scalability makes RMG a compelling tool for high-throughput DFT studies of materials that would otherwise be bottlenecked in other codes (for example, by global fast Fourier transforms in plane-wave DFT). However, the limited workflow infrastructure for RMG has thus far constrained its adoption to a small user community. In this work, we present pyRMG, a Python package designed to streamline the setup and execution of RMG DFT calculations. Built on the pymatgen and ASE (Atomic Simulation Environment) computational materials science Python packages, pyRMG automates input generation and convergence checking, and it integrates with modern job schedulers (e.g., Flux) on leadership-class platforms such as Frontier and Perlmutter. Here, we demonstrate pyRMG for a high-throughput study of strain effects in 2D 2L-Bi 2 Se 3 /2L-NbSe 2 heterostructures, which offers chemical insights into this system and shows that RMG-based workflows can converge with limited user intervention.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nuclear, Land marker for a Waste Isolation Site: International Architecture Competition - 18240

The design and construction of a nuclear waste marker must include many architectural considerations. Very little research has been conducted by architects to address this issue to date. Historically, the task of designating a nuclear waste site has been a discussion predominantly contained within the scientific and engineering communities. In an effort to engage architects within this dialogue, arch out loud hosted the NUCLEAR competition as an international, architectural investigation into the design of a nuclear waste site marker. The NUCLEAR competition asks participants to design a marker to designate the site of the Waste Isolation Pilot Plant (WIPP) located near Carlsbad, New Mexico. The marker will be meant to be a means of passive institutional controls (PICs) to deter inadvertent human intrusion upon the site and should last 10,000 years. The NUCLEAR competition defined the purpose of the nuclear waste site marker as a means to communicate the danger of inadvertent intrusion from future generations. The competition requires that entries must be both legible in communication and durable in their construction. The competition brief challenged participants to consider how customs and languages will evolve over time and to predict what motivations future generations might have for digging at the site. Other background information was provided about the site such as how waste was stored and how the unique geology of the site seals the waste. The 1993 Sandia National Laboratories report titled Expert Judgment on Markers to Deter Inadvertent Human Intrusion into the Waste Isolation Pilot Plant [1] was referenced along with the 1992 WIPP Land Withdrawal Act [2], interviews with WIPP scientists [3] involved with the WIPP PICs program and the Nuclear Energy Agency Records Knowledge and Memory working group. Competition participants were given digital 3-D Computer-Aided Design models for reference that included the underground storage portion of the site, existing buildings and the 16-square-mile site boundary identified in the 1992 WIPP Land Withdrawal Act [2]. This paper will: - Provide a detailed description of the submitted proposals including graphic analysis of the visualizations, diagrams, technical drawings, jury comments and statements. - Catalogue and identify the consistent themes within all submissions with possible categories including: 1) overall architectural form and language, 2) integration of written language and symbols within architecture, 3) use of technology and robotics, 4) construction methods and materials. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Analysis of Plane-Parallel Electron Beam Propagation in Different Media by Numerical Simulation Methods

Simulation by the Monte Carlo method is widely used to calculate the character of ionizing radiation interaction with substance. A wide variety of programs based on the given method allows users to choose the most suitable package for solving computational problems. In turn, it is important to know exactly restrictions of numerical systems to avoid gross errors. Results of estimation of the feasibility of application of the program PCLab (Computer Laboratory, version 9.9) for numerical simulation of the electron energy distribution absorbed in beryllium, aluminum, gold, and water for industrial, research, and clinical beams are presented. The data obtained using programs ITS and Geant4 being the most popular software packages for solving the given problems and the program PCLab are presented in the graphic form. A comparison and an analysis of the results obtained demonstrate the feasibility of application of the program PCLab for simulation of the absorbed energy distribution and dose of electrons in various materials for energies in the range 1–20 MeV.

61 RADIATION PROTECTION AND DOSIMETRY↗

Mixed-Precision S/DGEMM Using the TF32 and TF64 Frameworks on Low-Precision AI Tensor Cores

Using NVIDIA graphics processing units (GPUs) equipped with Tensor Cores has enabled the significant acceleration of general matrix multiplication (GEMM) for applications in machine learning (ML) and artificial intelligence (AI) and in high-performance computing (HPC) generally. The use of such power-efficient, specialized accelerators can provide a performance increase between 8 × and 20 ×, albeit with a loss in precision. However, a high level of precision is required in many large scientific and HPC applications, and computing in single or double precision is still necessary for many of these applications to maintain accuracy. Fortunately, mixed-precision methods can be employed to maintain a higher level of numerical precision while also taking advantage of the performance increases from computing with lower-precision AI cores. With this in mind, we extend the state of the art by using NVIDIA’s new TF32 framework. This new framework not only burdens some constraints of the previous frameworks, such as costly 32 16-bit castings but also provides an equivalent precision and performance by using a much simpler approach. We also propose a new framework called TF64 that attempts double-precision arithmetic with low-precision Tensor Cores. Although this framework does not exist yet, we validated the correctness of this idea and achieved an equivalent of 64-bit precision on 32-bit hardware.

Valero Lara, Pedro↗

HipGISAXS

Implementation of a flexible Grazing Incidence Small-Angle Scattering (GISAXS). Simulation code based on the Distorted Wave Born Approximation (DWBA) theory that effectively utilizes the parallel processing power provided by graphic processors and multicore processors. This constitutes a handy tool for experimentalists facing a massive flux of data, allowing them to accurately simulate the GISAXS process and analyze the produced data. The software computes the diffraction image for any given super position of custom shapes or morphologies in a user-defined region of the reciprocal space for all possible grazing incidence angles and in-plane sample rotations. This flexibility then allows to easily tackle a wide range of possible sample structures such as nanoparticles o top of or embedded in a substrate or a multilayered structure. In cases where the sample displays regions of significant refractive index contrast, an algorithm has been implemented to perform a slicing of the sample along the vertical direction and compute the averaging refractive index profile to be used as the reference geometry of the unpertubed system. Preliminary tests show good agreement with experimental data for a variety of commonly encountered nanostructures.

Xiaoye, Sherry Li↗

The MARS Code System User’s Guide Version 15(2016)

This paper is a user’s guide to the current version of the MARS15 Monte Carlo code. MARS allows inclusive and, in many cases, exclusive simulations of three-dimensional hadronic and electromagnetic cascades and modeling of heavy ion, muon and low energy neutron-photon transport in accelerator, detector, spacecraft and shielding components. Particles can have energies from a fraction of an electronvolt up to about 100 TeV. The code has undergone substantial improvements since the last documented version MARS13(95) and since previous releases of MARS15; all these as well as other specific features of the MARS code system are explained in detail. Descriptions of general input and output with commentary and recommendations are given. Examples are given for running the program with distributed sources, complex compounds, arbitrary geometries, and magnetic fields. A built-in graphical-user interface allows the visualization and detailed analysis of the geometry description. Enhanced options are described for the MARS code being coupled with: the DPMJET hadron event generator, the MCNP code for low-energy neutron transport, the ANSYS system for thermal and stress analysis, physics analysis and graphics packages, the MAD accelerator lattice description via a universal MAD-MARS beam line builder, the STRUCT program for tracking particles in accelerator lattices with beam loss recording, the DETRA code for nuclide decay and transmutation calculations, the EGS5 code for precise modeling of low-energy electromagnetic showers, and the ROOT system for description and visualization of arbitrary geometries. Use of the code in a multistage mode, coupled with event generators (DPMJET), with the STRUCT program for tracking particles in accelerator lattices with beam loss recording, and with physics analysis and graphics packages is demonstrated with typical input and output examples.

43 PARTICLE ACCELERATORS↗

Fluctuating Hydrodynamics Methods for Dynamic Coarse-Grained Implicit-Solvent Simulations in LAMMPS

We present a software package integrated with the molecular dynamics software LAMMPS for fluctuating hydrodynamics simulations of fluid-structure interactions subject to thermal fluctuations. The package is motivated to provide dynamic thermostats to extend implicit-solvent coarse-grained (IS-CG) models by incorporating kinetic contributions from the solvent to facilitate their use in a wider range of applications. To capture the thermal and hydrodynamic contributions of the solvent to dynamics, we introduce momentum conserving thermostats and computational methods based on fluctuating hydrodynamics and the stochastic Eulerian Lagrangian method (SELM). SELM couples the coarse-grained microstructure degrees of freedom to continuum stochastic fields to capture both the relaxation of hydrodynamic modes and thermal fluctuations. Features of the SELM software include (i) numerical time-step integrators for SELM fluctuating hydrodynamics in inertial and quasi-steady regimes, (ii) Lees--Edwards-style methods for imposing shear, (iii) a Java-based graphical user interface (GUI) for setting up models and simulations, (iv) standardized XML formats for parametrization and data output, and (v) standardized formats VTK for continuum fields and microstructures. The SELM software package facilitates for pre-established models in LAMMPS easy adoption of the SELM fluctuating hydrodynamics thermostats. We offer here an overview of the SELM software package, computational methods, and applications.

97 MATHEMATICS AND COMPUTING↗

SCALE Code System 6.2.1

The SCALE Code System is a widely-used modeling and simulation suite for nuclear safety analysis and design that is developed, maintained, tested, and managed by the Reactor and Nuclear Systems Division (RNSD) of Oak Ridge National Laboratory (ORNL). SCALE provides a comprehensive, verified and validated, user-friendly tool set for criticality safety, reactor and lattice physics, radiation shielding, spent fuel and radioactive source term characterization, and sensitivity and uncertainty analysis. Since 1980, regulators, licensees, and research institutions around the world have used SCALE for safety analysis and design. SCALE provides an integrated framework with dozens of computational modules including three deterministic and three Monte Carlo radiation transport solvers that are selected based on the desired solution strategy. SCALE includes current nuclear data libraries and problem-dependent processing tools for continuous-energy (CE) and multigroup (MG) neutronics and coupled neutron-gamma calculations, as well as activation, depletion, and decay calculations. SCALE includes unique capabilities for automated variance reduction for shielding calculations, as well as sensitivity and uncertainty analysis. SCALE’s graphical user interfaces assist with accurate system modeling, visualization of nuclear data, and convenient access to desired results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Enabling real-time adaptation of machine learning models at x-ray Free Electron Laser facilities with high-speed training optimized computational hardware

The emergence of novel computational hardware is enabling a new paradigm for rapid machine learning model training. For the Department of Energy’s major research facilities, this developing technology will enable a highly adaptive approach to experimental sciences. In this manuscript we present the per-epoch and end-to-end training times for an example of a streaming diagnostic that is planned for the upcoming high-repetition rate x-ray Free Electron Laser, the Linac Coherent Light Source-II. We explore the parameter space of batch size and data parallel training across multiple Graphics Processing Units and Reconfigurable Dataflow Units. We show the landscape of training times with a goal of full model retraining in under 15 min. Although a full from scratch retraining of a model may not be required in all cases, we nevertheless present an example of the application of emerging computational hardware for adapting machine learning models to changing environments in real-time, during streaming data acquisition, at the rates expected for the data fire hoses of accelerator-based user facilities.

97 MATHEMATICS AND COMPUTING↗

Exploring the interoperability of remote GPGPU virtualization using rCUDA and directive-based programming models

Directive-based programming models, such as OpenMP, OpenACC, and OmpSs, enable users to accelerate applications by using coprocessors with little effort. These devices offer significant computing power, but their use can introduce two problems: an increase in the total cost of ownership and their underutilization because not all codes match their architecture. Remote accelerator virtualization frameworks address those problems. In particular, rCUDA provides transparent access to any graphic processor unit installed in a cluster, reducing the number of accelerators and increasing their utilization ratio. Joining these two technologies, directive-based programming models and rCUDA, is thus highly appealing. In this work, we study the integration of OmpSs and OpenACC with rCUDA, describing and analyzing several applications over three different hardware configurations that include two InfiniBand interconnections and three NVIDIA accelerators. Our evaluation reveals favorable performance results, showing low overhead and similar scaling factors when using remote accelerators instead of local devices.

97 MATHEMATICS AND COMPUTING↗

Implications of doping on microstructure, processing, and thermoelectric performance: The case of PbSe

Abstract In this work, we highlight the often-overlooked effects of doping on the microstructure and performance of bulk thermoelectric materials to offer a broader perspective on how dopants interact with their parent material. Using PbSe doped with Na, Ag, and K as a model material system, we combine original computational, experimental, and microscopy data with established trends in material behavior, to provide an in-depth discussion of the relationship between dopants, processing, and microstructure, and their effects on thermoelectric efficiency and thermal stability. Notable observations include differences in the microstructure and mass loss of thermally treated samples of Na- and Ag-doped PbSe, as well as findings that Na and K cations exist predominantly as substitutional point defects while Ag also occupies interstitial sites and exhibits lower solubility. We discuss how these differences in point defect populations are known to affect a dopants’ ability to alter carrier concentration and how they may affect the mechanical properties of PbSe during processing. Graphic Abstract

Grovogui, Jann A. (ORCID:0000000200481049)↗

Surface Triangulation for CSG in Mercury

Visualization routines for rendering complicated geometries are very useful for engineers and scientists who are trying to build 3D prototypes of their designs. A common way to rapidly add interesting features to a 3D model is through the use of a concept called Constructive Solid Geometry. CSG uses compositions of the boolean set operations to manipulate basic geometric primitives to form more complicated objects. The most common boolean operations employed are union, intersection, and subtraction. Most computer-aided design software packages contain some sort of ability visualize CSG. The typical workflow for the user is as follows: The user specifies the individual primitive components, the user arbitrarily combines each of these primitives with boolean operations, the software generates a CSG tree structure which normally stores these solids implicitly with their defining equation, the tree is traversed and a general algorithm is applied to render the appropriate geometry onto the screen. Algorithms for visualizing CSG have been extensively developed for over a decade. Points sampled from the implicit solids are typically used as input by variations of algorithms like marching cubes and point-cloud surface reconstruction. Here, we explain a surface triangulation method from the graphics community that is being used for surface visualization in the framework of a Monte-Carlo neutron transport code called Mercury.

97 MATHEMATICS AND COMPUTING↗

Analysis of the performance of a hybrid CPU/GPU 1D2D coupled model for real flood cases

Coupled 1D2D models emerged as an efficient solution for a two-dimensional (2D) representation of the floodplain combined with a fast one-dimensional (1D) schematization of the main channel. At the same time, high-performance computing (HPC) has appeared as an efficient tool for model acceleration. In this work, a previously validated 1D2D Central Processing Unit (CPU) model is combined with an HPC technique for fast and accurate flood simulation. Due to the speed of 1D schemes, a hybrid CPU/GPU model that runs the 1D main channel on CPU and accelerates the 2D floodplain with a Graphics Processing Unit (GPU) is presented. Since the data transfer between sub-domains and devices (CPU/GPU) may be the main potential drawback of this architecture, the test cases are selected to carry out a careful time analysis. Here, the results reveal the speed-up dependency on the 2D mesh, the event to be solved and the 1D discretization of the main channel. Additionally, special attention must be paid to the time step size computation shared between sub-models. In spite of the use of a hybrid CPU/GPU implementation, high speed-ups are accomplished in some cases.

54 ENVIRONMENTAL SCIENCES↗

SU-F-T-242: A Method for Collision Avoidance in External Beam Radiation Therapy

Purpose: We proposed a method for collision avoidance (CA) in external beam radiation therapy (EBRT). The method encompasses the analysis of all positions of the moving components of the beam delivery system such as the treatment table and gantry, including patient specific information obtained from the CT images. This method eliminates the need for time-consuming dry-runs prior to the actual treatments. Methods: The QA procedure for EBRT requires that the collision should be checked prior to treatment. We developed a system capable of a rigorous computer simulation of all moving components including positions of the couch and gantry during the delivery, position of the patients, and imaging equipment. By running this treatment simulation it is possible to quantify and graphically represent all positions and corresponding trajectories of all points of the moving parts during the treatment delivery. The development of the workflow for implementation of the CA includes several steps: a) derivation of combined dynamic equation of motion of the EBRT delivery systems, b) developing the simulation model capable of drawing the motion trajectories of the specific points, c) developing the interface between the model and the treatment plan parameters such as couch and gantry parameters for each field. Results: The patient CT images were registered to the treatment couch so the patient dimensions were included into the simulation. The treatment field parameters were structured in the xml-file which was used as the input into the dynamic equations. The trajectories of the moving components were plotted on the same graph using the dynamic equations. If the trajectories intersect that was the signal that collision exists. Conclusion: This CA method was proved to be effective in the simulation of treatment delivery. The proper implementation of this system can potentially improve the QA program and increase the efficacy in the clinical setup.

60 APPLIED LIFE SCIENCES↗

Optimization Studies of Radiation Shielding for the PIP-II Project at Fermilab

The PIP-II project at Fermilab, which includes an 800-MeV superconducting LINAC, demands rigorous radiation shielding optimization to meet safety requirements. We updated the MARS geometry model to reflect new magnet and collimator designs and introduced high-resolution detector planes to better capture radiation field distributions. To overcome the significant computational demands, we implemented a well-known branching technique that drastically reduced simulation runtimes while maintaining statistical integrity. This was achieved through particle splitting and the application of Russian Roulette techniques. Additionally, new graphical tools were created to streamline data visualization and MARS code usability.

Makovec, Alajos [Fermilab] (ORCID:0000000286157492↗

Providing a Graphical Tool for Modeling Reactor Cores

Complex, coupled physical systems that occur in various domain sciences present several serious challenges in implementing an efficient and seamless simulation workflow. These challenges impact the productivity of scientific research and the overall time it takes to verify simulation results. To overcome these computational modeling barriers, specifically in nuclear reactor analysis, one must rely on flexible and scalable tools that simplify modeling complexity.

97 MATHEMATICS AND COMPUTING↗