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 613 records · Page 34

High-Resolution Simulations of Geological CO 2 Injection: Application to the SPE11 Benchmark

Geological carbon sequestration (GCS) will play a critical role in decarbonization and in facilitating the transition to clean energy systems. Because CO 2 is highly mobile, ensuring its safe and permanent injection into subsurface geological formations involves monitoring over larger spatial domains and longer time periods than is typical for hydrocarbon reservoirs. This can benefit from simulation tools capable of modeling key CO 2 trapping mechanisms, particularly those optimized for speed and scalability on high-performance computing systems. Using isothermal versions of the SPE11B and SPE11C benchmark cases, we conduct a mesh refinement study simulating CO 2 injection into kilometer-scale rock formations at centimeter resolution with the GEOS open-source simulation framework. We focus on how mesh refinement improves the accuracy of convective mixing in both 2D and 3D simulations. The computational costs associated with achieving a converged solution highlight the need for predictive upscaling techniques. A systematic performance scaling analysis—including both central processing unit (CPU) and graphics processing unit (GPU) architectures—complements the “Results” section.

Geosciences↗

Graph 500 in OpenSHMEM, In: OpenSHMEM 2014: OpenSHMEM and Related Technologies. Experiences, Implementations, and Technologies

This document describes the effort to implement the Graph 500 benchmark using OpenSHMEM based on the MPI-2 one-side version. The Graph 500 benchmark performs a breadth-first search in parallel on a large randomly generated undirected graph and can be implemented using basic MPI-1 and MPI-2 one-sided communication. Graph 500 requires atomic bit-wise operations on unsigned long integers but neither atomic bit-wise operations nor OpenSHMEM for unsigned long are available in OpenSHEM. Such needed bit-wise atomic operations and support for unsigned long are implemented using atomic condition swap (CSWAP) on signed long integers. Preliminary results on comparing the OpenSHMEM and MPI-2 one-sided implementations on a Silicon Graphics Incorporated (SGI) cluster and the Cray XK7 are presented.

Computer Science↗

Energy of Motion of Internal and Surface Waves in a Two-Layer Hydrodynamic System

We study the energy of wave motion in a two-layer “layer with solid bottom–layer with free surface” hydrodynamic system. The dependences of the energies of motion of internal and surface waves on the geometric and physical parameters of the system are analyzed. The contribution of the second approximation to the total energy of the system is estimated. The reliability of the obtained results is verified. The results are graphically illustrated.

42 ENGINEERING↗

Sparsity-Independent Lyapunov Exponent in the Sachdev-Ye-Kitaev Model

The saturation of a recently proposed universal bound on the Lyapunov exponent has been conjectured to signal the existence of a gravity dual. This saturation occurs in the low-temperature limit of the dense Sachdev-Ye-Kitaev (SYK) model, N Majorana fermions with q body ( q > 2 ) infinite-range interactions. We calculate certain out-of-time-order correlators (OTOCs) for N ≤ 64 fermions for a highly sparse SYK model and find no significant dependence of the Lyapunov exponent on sparsity up to near the percolation limit where the Hamiltonian breaks up into blocks. This provides strong support to the saturation of the Lyapunov exponent in the low-temperature limit of the sparse SYK. A key ingredient to reaching N = 64 is the development of a novel quantum spin model simulation library that implements highly optimized matrix-free Krylov subspace methods on graphical processing units. This leads to a significantly lower simulation time as well as vastly reduced memory usage over previous approaches, while using modest computational resources. Strong sparsity-driven statistical fluctuations require both the use of a much larger number of disorder realizations with respect to the dense limit and a careful finite size scaling analysis. The saturation of the bound in the sparse SYK points to the existence of a gravity analog that would enlarge substantially the number of field theories with this feature. Published by the American Physical Society 2024

Physics↗

RUScal : Software for the analysis of resonant ultrasound spectroscopy measurements

Resonant ultrasound spectroscopy is used to nondestructively measure the elastic resonances of small solids to elucidate the material's elastic properties or other qualities like size, shape, or composition. In this work, we introduce the software RUScal for the purpose of determining elastic properties by analyzing the eigenfrequencies of solid specimens with common shapes, such as rectangular parallelepipeds, cylinders (solid and hollow tube), ellipsoids, and octahedrons, as well as irregularly shaped ellipsoids that can be described analytically. All symmetry classes are supported, from isotropic to triclinic, along with the option to add or remove up to three orthogonal mirror planes as well as the ability to reorient the crystal axes with respect the sample edges via Euler angles. Additional features include tools to help find initial sets of elastic constants, including grid exploration and Monte Carlo methods, a tool to analyze frequencies as a function of sample length or crystal orientation, an error analysis tool to assess fit quality, and formatting of the input and output files for batch fitting, e.g., as a function of temperature. This software was validated with published resonant ultrasound spectroscopy data for various materials, shapes, and symmetries with noted improvements in calculation time compared to finite element methods.

47 OTHER INSTRUMENTATION↗

Explainable and trustworthy artificial intelligence for correctable modeling in chemical sciences

Data science has primarily focused on big data, but for many physics, chemistry, and engineering applications, data are often small, correlated and, thus, low dimensional, and sourced from both computations and experiments with various levels of noise. Typical statistics and machine learning methods do not work for these cases. Expert knowledge is essential, but a systematic framework for incorporating it into physics-based models under uncertainty is lacking. Here, we develop a mathematical and computational framework for probabilistic artificial intelligence (AI)–based predictive modeling combining data, expert knowledge, multiscale models, and information theory through uncertainty quantification and probabilistic graphical models (PGMs). We apply PGMs to chemistry specifically and develop predictive guarantees for PGMs generally. Our proposed framework, combining AI and uncertainty quantification, provides explainable results leading to correctable and, eventually, trustworthy models. The proposed framework is demonstrated on a microkinetic model of the oxygen reduction reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computationally Efficient Multiscale Neural Networks Applied to Fluid Flow in Complex 3D Porous Media

Abstract The permeability of complex porous materials is of interest to many engineering disciplines. This quantity can be obtained via direct flow simulation, which provides the most accurate results, but is very computationally expensive. In particular, the simulation convergence time scales poorly as the simulation domains become less porous or more heterogeneous. Semi-analytical models that rely on averaged structural properties (i.e., porosity and tortuosity) have been proposed, but these features only partly summarize the domain, resulting in limited applicability. On the other hand, data-driven machine learning approaches have shown great promise for building more general models by virtue of accounting for the spatial arrangement of the domains’ solid boundaries. However, prior approaches building on the convolutional neural network (ConvNet) literature concerning 2D image recognition problems do not scale well to the large 3D domains required to obtain a representative elementary volume (REV). As such, most prior work focused on homogeneous samples, where a small REV entails that the global nature of fluid flow could be mostly neglected, and accordingly, the memory bottleneck of addressing 3D domains with ConvNets was side-stepped. Therefore, important geometries such as fractures and vuggy domains could not be modeled properly. In this work, we address this limitation with a general multiscale deep learning model that is able to learn from porous media simulation data. By using a coupled set of neural networks that view the domain on different scales, we enable the evaluation of large ( $$>512^3$$ > 512 3 ) images in approximately one second on a single graphics processing unit. This model architecture opens up the possibility of modeling domain sizes that would not be feasible using traditional direct simulation tools on a desktop computer. We validate our method with a laminar fluid flow case using vuggy samples and fractures. As a result of viewing the entire domain at once, our model is able to perform accurate prediction on domains exhibiting a large degree of heterogeneity. We expect the methodology to be applicable to many other transport problems where complex geometries play a central role.

36 MATERIALS SCIENCE↗

Final report for the DOE Early Career Award #DE-SC0003912

This DoE supported early career project was aimed at developing computational models, theory and simulation methods that would be then be used to predict assembly and morphology in polymer nanocomposites. In particular, the focus was on composites in active layers of devices, containing conducting polymers that act as electron donors and nanoscale additives that act as electron acceptors. During the course this work, we developed the first of its kind molecular models to represent conducting polymers enabling simulations at the experimentally relevant length and time scales. By comparison with experimentally observed morphologies we validated these models. Furthermore, using these models and molecular dynamics simulations on graphical processing units (GPUs) we predicted the molecular level design features in polymers and additive that lead to morphologies with optimal features for charge carrier behavior in solar cells. Additionally, we also predicted computationally new design rules for better dispersion of additives in polymers that have been confirmed through experiments. Achieving dispersion in polymer nanocomposites is valuable to achieve controlled macroscopic properties of the composite. The results obtained during the course of this DOE funded project enables optimal design of higher efficiency organic electronic and photovoltaic devices and improve every day life with engineering of these higher efficiency devices.

14 SOLAR ENERGY↗

TAMM: Tensor algebra for many-body methods

Tensor algebra operations such as contractions in computational chemistry consume a significant fraction of the computing time on large-scale computing platforms. The widespread use of tensor contractions between large multi-dimensional tensors in describing electronic structure theory has motivated the development of multiple tensor algebra frameworks targeting heterogeneous computing platforms. In this paper, we present Tensor Algebra for Many-body Methods (TAMM), a framework for productive and performance-portable development of scalable computational chemistry methods. TAMM decouples the specification of the computation from the execution of these operations on available high-performance computing systems. With this design choice, the scientific application developers (domain scientists) can focus on the algorithmic requirements using the tensor algebra interface provided by TAMM, whereas high-performance computing developers can direct their attention to various optimizations on the underlying constructs, such as efficient data distribution, optimized scheduling algorithms, and efficient use of intra-node resources (e.g., graphics processing units). The modular structure of TAMM allows it to support different hardware architectures and incorporate new algorithmic advances. We describe the TAMM framework and our approach to the sustainable development of scalable ground- and excited-state electronic structure methods. We present case studies highlighting the ease of use, including the performance and productivity gains compared to other frameworks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Combined experimental and theoretical assessments of the lattice dynamics and optoelectronics of TaON and Ta{sub 3}N{sub 5}

Presented herein is a detailed discussion of the properties of the lattice dynamic and optoelectronic properties of tantalum(V) oxynitride (TaON) and tantalum(V) nitride (Ta{sub 3}N{sub 5}), from experimental and theoretical standpoint. The active Raman and infra red (IR) frequencies of TaON and Ta{sub 3}N{sub 5} were measured using confocal Raman and Fourier Transform Infrared spectroscopies (FTIR) and calculated using the linear response method within the density functional perturbation theory (DFPT). The detailed study leads to an exhaustive description of the spectra, including the symmetry of the vibrational modes. Electronic structures of these materials were computed using DFT within the range-separated hybrid HSE06 exchange–correlation formalism. Electronic and ionic contributions to the dielectric constant tensors of these materials were obtained from DFPT within the linear response method using the PBE functional. Furthermore, effective mass of photogenerated holes and electrons at the band edges of these compounds were computed from the electronic band structure obtained at the DFT/HSE06 level of theory. The results suggest that anisotropic nature in TaON and Ta{sub 3}N{sub 5} is present in terms of dielectric constant and effective masses. - Graphical abstract: Detailed investigation has been conducted from combined experimental and theoretical approaches on Raman and IR spectroscopies, electronic structures, dielectric constants and effective masses of TaON and Ta{sub 3}N{sub 5}. - Highlights: • Crystal structures of TaON and Ta{sub 3}N{sub 5} are discussed based on XRD and DFT calculation. • Raman and IR spectra of TaON and Ta{sub 3}N{sub 5} materials are measured and computed by DFPT method. • Optoelectronic properties of TaON and Ta{sub 3}N{sub 5} are discussed. • Dielectric constant and effective masses of TaON and Ta{sub 3}N{sub 5} are calculated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Finite-temperature many-body perturbation theory for anharmonic vibrations: Recursions, algebraic reduction, second-quantized reduction, diagrammatic rules, linked-diagram theorem, finite-temperature self-consistent field, and general-order algorithm

A unified theory is presented for finite-temperature many-body perturbation expansions of the anharmonic vibrational contributions to thermodynamic functions, i.e., the free energy, internal energy, and entropy. The theory is diagrammatically size-consistent at any order, as ensured by the linked-diagram theorem proved in this study, and, thus, applicable to molecular gases and solids on an equal footing. It is also a basis-set-free formalism, just like its underlying Bose–Einstein theory, capable of summing anharmonic effects over an infinite number of states analytically. It is formulated by the Rayleigh–Schrödinger-style recursions, generating sum-over-states formulas for the perturbation series, which unambiguously converges at the finite-temperature vibrational full-configuration-interaction limits. Two strategies are introduced to reduce these sum-over-states formulas into compact sum-over-modes analytical formulas. One is a purely algebraic method that factorizes each many-mode thermal average into a product of one-mode thermal averages, which are then evaluated by the thermal Born–Huang rules. Canonical forms of these rules are proposed, dramatically expediting the reduction process. The other is finite-temperature normal-ordered second quantization, which is fully developed in this study, including a proof of thermal Wick’s theorem and the derivation of a normal-ordered vibrational Hamiltonian at finite temperature. The latter naturally defines a finite-temperature extension of size-extensive vibrational self-consistent field theory. These reduced formulas can be represented graphically as Feynman diagrams with resolvent lines, which include anomalous and renormalization diagrams. Two order-by-order and one general-order algorithms of computing these perturbation corrections are implemented and applied up to the eighth order. The results show no signs of Kohn–Luttinger-type nonconvergence.

74 ATOMIC AND MOLECULAR PHYSICS↗

ART Front End v1

Enhances accessibility of the Automated Recommendation Tool (ART), by adding on a graphical user interface, user authentication, and by exposing a REST API. Enables authenticated remote users and applications to interface with ART and to leverage remote computing resources to perform runs.

Forrer, Mark↗

Levelized Cost of Avoided CO2 tool (LCAC tool) v0.9

The levelized cost of avoided CO2 (LCAC) tool facilitates techno-economic calculations for various decarbonization technologies. It comprises two components to accomplish this: an input workbook based on MS Excel and a web-based output tool. The workbook acts as the input sheet, using the user-provided decarbonization measures and their impacts on lifetime energy consumption, energy costs, and CO2 emissions to compute LCAC (expressed in $/metric ton of avoided CO2). The workbook is uploaded to a web-based tool to generate an LCAC curve. This graphical representation illustrates the economic comparison of all decarbonization measures for purposes such as implementation prioritization. We are unaware of any publicly available tools with similar capabilities.

Ali Qamar, Muhammad↗

GPU Lifetimes on Titan Supercomputer: Survival Analysis and Reliability

George Ostrouchov, Don Maxwell, Rizwan Ashraf, Mallikarjun Shankar, and James Rogers. 2020. GPU Lifetimes on Titan Supercomputer: Survival Analysis and Reliability. In Proceedings of the International Conference for High Performance Computing, Networking, Storage and Analysis (SC '20). Association for Computing Machinery, New York, NY, USA. Data and code for SC20 paper about Titan GPU reliability analysis: https://github.com/olcf/TitanGPULife. Includes R code to generate graphics for paper and additional analyses. See code/README for instructions. Includes original Titan GPU reliability data on over 100,000 collective hours of operation: data/titan.gpu.history.txt - history data, data/titan.service.txt - service nodes for exclusion. Includes output data files produced by code/TitanGPUmodel.Rmd: data/gc_full.csv - cleaned up data (see paper and R code); data/gc_summary_loc.csv - one record per GPU (variables: SN, time, nlife, nloc, last, col, row, cage, slot, node, max_loc_events, time_max_loc, dbe, dbe_loc, otb, otb_loc, out, batch, days, years, dead, dead_otb, dead_dbe) (see paper and R code). Includes .Rmd analysis document as TitanGPUmode.html. Includes Python code to process data/gc_full.csv into graphics from time-between-failure analyses: See code/tbf-analyses/README for instructions.

42 ENGINEERING↗

GPU Acceleration of Volume Fraction and Centroid Computation from General Shapes on Unstructured Mesh

One important requirement for coupling multi-physics models is a means to remap fields between meshes in an accurate and conservative manner. Since the remapping algorithm is embarrassingly parallel, this motivates the use of Graphical Processing Units (GPU) for such routine. However, since our remapping code, portage, is fairly complex, we explore the use of GPUs through a simpler problem of generating material’s volume fraction and centroid inside each cell of the mesh. We use the same programming structure that is used in portage, namely, on-node parallelism via NVIDIA’s thrust library. Using a similar programming structure enables us to draw lessons from this study that can be applied to the more complex problem. We have found that while the speed-up is pronounced (10-100 times), significant effort is required to re-factor the code so that thrust can be used with a CUDA backend. Moreover, thrust takes away the fine grain control that vanilla CUDA programming gives us.

97 MATHEMATICS AND COMPUTING↗

Hydrothermal reactions: From the synthesis of ligand to new lanthanide 3D-coordination polymers

The organic ligand 2,5-piperazinedione-1,4-diacetic acid (H{sub 2}PDA) was synthesized under hydrothermal conditions starting from the iminodiacetic acid and catalyzed by oxalic acid. The X-ray powder diffraction data indicates that the compound crystallizes in the P2{sub 1}/c monoclinic system as reported in the literature. The ligand was also characterized by elemental analysis, magnetic nuclear resonance, infrared spectroscopy and thermogravimetric analysis. Two new coordination networks based on lanthanide ions were obtained with this ligand using hydrothermal reaction. In addition to single-crystal X-ray diffraction, the compounds were characterized by infrared spectroscopy, thermogravimetric analysis, scanning electron microscopy and elemental analysis. Single-crystal XRD showed that the compounds are isostructural, crystallizing in P2{sub 1}/n monoclinic system with chemical formula [Ln(PDA){sub 1.5}(H{sub 2}O)](H{sub 2}O){sub 3} (Ln=Gd{sup 3+}(1) and Eu{sup 3+}(2)).The luminescence properties of both compounds were studied. In the compound (1), a broad emission band was observed at 479 nm, redshifted by 70 nm in comparison of the free ligand. In (2), the typical f–f transition was observed with a maximum peak at 618 nm, related with the red emission of the europium ions. Computational methods were performed to simulate the crystal structure of (2). The theoretical calculations of the intensity parameters are in good agreement with the experimental values. - Graphical abstract: Scheme of obtaining the ligand 2,5-piperazinedione-1,4-diacetic acid (H{sub 2}PDA) and two new isostructural 3D-coordination polymers [Ln(PDA){sub 1.5}(H{sub 2}O)](H{sub 2}O){sub 3} (Ln=Gd{sup 3+} and Eu{sup 3+}) by hydrothermal synthesis. Display Omitted - Highlights: • The ligand 2,5-piperazinedione-1,4-diacetic acid was synthetized using the hydrothermic method and characterized. • Two new 3D-coordination polymers with this ligand containing Gd{sup 3+} and Eu{sup 3+} ions were also obtained. • The compounds are isostructural and the typical luminescent properties were observed in both structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Complex concentrated alloys and chemical short-range ordering

Abstract Complex concentrated alloys (CCAs) are materials comprising three or more elements in similar proportions and possessing structural but no chemical long-range order. Fascination with CCAs has grown over the last 20 years and to date, CCAs have opened a new materials design paradigm and horizon for discovery of materials to meet the demands of applications in aggressive environments. Understanding the fundamental mechanisms controlling their response, however, is challenging due to the chemical and structural variations that wildly fluctuate over fine atomic and nanoscales. This issue focuses on the experimental, computational, and theoretical investigations that aim to uncover phenomena and processes determining the structure, kinetics, mechanics, or deformation mechanisms in CCAs at the atomic scale. At the atomic scale at which they operate, chemical short-range ordering can be influential. This issue further addresses the capabilities, as well as the debatable need, to characterize, predict, and relate short-range ordering to material performance. Collectively, the articles in this issue highlight the insights, understanding, and experimental and computational tools that attempt to create property-tunable CCAs “from the atom up” by treating short-range ordering and engineering atomic-scale mechanisms. Graphical abstract

36 MATERIALS SCIENCE↗