SEARCH · Engineering Papers

Results for “tensor cores”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS1969

ELAS: A general-purpose computer program for the equilibrium problems of linear structures. Volume 2: Documentation of the program

A general purpose digital computer program for the in-core solution of linear equilibrium problems of structural mechanics is documented. The program requires minimum input for the description of the problem. The solution is obtained by means of the displacement method and the finite element technique. Almost any geometry and structure may be handled because of the availability of linear, triangular, quadrilateral, tetrahedral, hexahedral, conical, triangular torus, and quadrilateral torus elements. The assumption of piecewise linear deflection distribution insures monotonic convergence of the deflections from the stiffer side with decreasing mesh size. The stresses are provided by the best-fit strain tensors in the least squares at the mesh points where the deflections are given. The selection of local coordinate systems whenever necessary is automatic. The core memory is used by means of dynamic memory allocation, an optional mesh-point relabelling scheme and imposition of the boundary conditions during the assembly time.

Utku, S.↗

DOE OSTI2021

Calculating the hyperfine tensors for group-IV impurity-vacancy centers in diamond using hybrid density functional theory

The hyperfine interaction is an important probe for understanding the structure and symmetry of defects in a semiconductor. Density functional theory has shown that it can provide useful first-principles predictions for both the hyperfine tensor and the hyperfine constants that arise from it. Recently there has been great interest in using group-IV impurity-vacancy color centers XV – (where X = Si, Ge, Sn, or Pb and V is a carbon vacancy) for important applications in quantum computing and quantum information science. In this paper we have calculated the hyperfine tensors for these XV – color centers using the HSE06 screened Hartree-Fock hybrid exchange-correlation functional with the inclusion of core electron spin polarization. We have compared our results to calculations which only use the PBE exchange-correlation functional without the inclusion of core electron spin polarization and we have found that our results are in very good agreement with available experimental results. Lastly, we have theoretically shown that these XV – color centers exhibit a Jahn-Teller distortion which explains the observed anisotropic distribution of the hyperfine constants among the neighboring 13 C nuclear spins.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Density matrix renormalization group description of the island of inversion isotopes F 28 – 33

Recent experiments have confirmed that the neutron-rich isotopes 28,29 F belong to the so-called island of inversion (IOI), a region of the nuclear chart around Z = 10 and N = 20 where nuclear structure deviates from the standard shell model predictions due to deformation and continuum effects. However, while the general principles leading to the IOI are relatively well understood, the details of the low-lying structure of the exotic fluorine isotopes 28–33 F are basically unknown. In this study, we perform large-scale shell model calculations including continuum states to investigate the properties of the neutron-rich isotopes 25–33 F, from a core of 24 O and using an effective two-body interaction with a small number of adjustable parameters in the central and tensor channels. We develop two models adjusted on experimentally confirmed states in 25,26 O and 25–27 F based on different assumptions concerning the positions of the neutron 0d 3/2 and 1p 3/2 shells, and solve the many-body problem using the density matrix renormalization group (DMRG) method for open quantum systems in an sd–fp model space. We obtain the first detailed spectroscopy of 25–33F in the continuum and show how the interplay between continuum effects and deformation explains the recent data on 28,29 F. Several deformed one- and two-neutron halo states are predicted in 29,31 F, and we provide some information about the possible structure of the heaviest fluorine isotopes. We also suggest several experimental studies of interest to constraint models and test the present predictions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Improved Tensor Current Limit from 8 B 𝛽 Decay Including New Recoil-Order Calculations

A precision measurement of the 𝛽 + decay of 8 B was performed using the Beta-decay Paul Trap to determine the 𝛽−𝜈 angular correlation coefficient 𝑎 𝛽⁢𝜈 . The experimental results were combined with new ab initio symmetry-adapted no-core shell-model calculations to yield the second-most precise measurement from Gamow-Teller decays, 𝑎 𝛽⁢𝜈 = −0.3345 ± 0.001⁢9 stat ± 0.002⁢1 syst . This value agrees with the standard model value of −1/3 and improves uncertainties in 8 B by nearly a factor of 2. By combining results from 8 B and 8 Li , a tight limit on tensor current coupling to right-handed neutrinos was obtained. A recent global evaluation of all other precision 𝛽 decay studies suggested a nonzero value for right-handed neutrino coupling in contradiction with the standard model at just above 3⁢𝜎. Finally, the present results are of comparable sensitivity and do not support this finding.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1992

Near-wall reconstruction of higher order moments and length scales using the POD

An analysis of the near-wall behavior of the proper orthogonal decomposition (POD) eigenfunctions derived from direct numerical simulation (DNS) of channel flow is performed. Consistent with previous studies, a low order multi-mode reconstruction of the kinetic energy and Reynolds shear stress suffices. A similar reconstruction of the isotropic dissipation rate is shown to be insufficient, however. An analysis is performed of the multi-mode composition of the dissipation rate in the near-wall region, and it is shown that a significant number of higher-order modes are required to achieve the correct asymptotic consistency in the near-wall region. In an attempt to avoid this problem, a length scale definition is proposed in terms of an integration of the correlation tensor which factors in the presence of the wall. The wall is accounted for by only integrating out to 2y(+) and not over the entire domain. Viscous and inviscid estimates for the dissipation were used in the near-wall and core regions respectively, in conjunction with this length scale representation to obtain an estimate of the dissipation throughout the domain. The resulting dissipation exhibits the proper behavior near the wall and in the inertial layer. A 1 POD mode estimate of the length scale is computed and found to agree quite well with the length scale obtained when the entire correlation tensor is used.

Glauser, Mark N.↗

DOE OSTI2023

Using the Shallow Strain Tensor to Characterize Deep Geologic Reservoirs

Abstract Storing and recovering water, carbon, and heat from geologic reservoirs is central to managing resources in a changing climate. We tested the hypothesis that the strain tensor caused by injecting or producing fluids can be measured at shallow depths and interpreted to advance understanding of underlying deep aquifers or reservoirs. Geodetic‐grade strainmeters were deployed at 30 m depth overlying the Bartlesville Formation, a 500‐m‐deep sandstone near Tulsa, OK. The strainmeters are 220 m east of injection well 9A completed in a permeable lens at the base of the Bartlesville Formation. Water was injected into well 9A at approximately 1.0 L/s during four tests that ranged in duration from a few hours to a few weeks. The horizontal strain increased (tension) and the circumferential strain was a few times larger than the radial strain. The vertical strain decreased (compression) during injection. Strain rates were approximately 100 nε/day during the first few hours, but the rates decreased and were approximately 10 nε/day during most of the tests. Four independent methods of poroelastic simulation and inversion predict reservoir properties and geometries that are similar to each other and consistent with independent information about the reservoir. All strain interpretations predict that a boundary to the permeable lens occurs beneath the vicinity of the strainmeters, which is consistent with core data from the site. The boundary of the permeable lens is located by matching the vertical, radial and circumferential strains, which demonstrates the value of measuring the strain tensor.

Murdoch, Lawrence C.↗

DOE OSTI2021

A Non-Invasive Approach for Elucidating the Spatial Distribution of In Situ Stress in Deep Subsurface Geologic Formations Considered for CO 2 Storage (Task 5 Report Field Scale Stress Modeling)

This report describes research accomplishments achieved with funding provided through Department of Energy (DOE) Contract DE-FE0031686, for the project “A Non-Invasive Approach for Elucidating the Spatial Distribution of In Situ Stress in Deep Subsurface Geologic Formations Considered for CO 2 Storage.” The purpose of this DOE funding is to develop non-invasive methods to obtain important information about subsurface stresses. The overall research goal of the project was to develop and improve methods for determining the spatial distribution of stresses, including magnitude and orientation of the three principal stress components in the subsurface, based on new methods that extract stress information from seismic data combined with well measurements, extended with well tests and logs, and unified in a numerical model that permits computation of the full stress tensor throughout the domains under investigation.

58 GEOSCIENCES↗

DOE OSTI2024

Polarization near dislocation cores in SrTiO 3 single crystals: The role of flexoelectricity

Spontaneous polarization as large as ~28 μC/cm 2 was recently observed around the dislocation cores in non-polar SrTiO 3 bulk crystals, and its origin was attributed to the flexoelectric effect, i.e., polarization induced by strain gradients. However, the roles of flexoelectricity, relative to other electromechanical contributions, and the nature of dislocations, i.e., edge vs screw dislocations in the induced polarization, are not well understood. In this work, we study the role of flexoelectricity in inducing polarization around three types of dislocation cores in SrTiO 3 : b= a(100) edge dislocation, b= a(110) edge dislocation, and b= a(010) screw dislocation, where b is the Burgers vector. For the edge dislocations, polarization can be induced by electrostriction alone, while flexoelectricity is essential for stabilizing the symmetric polarization pattern. The shear component of the flexoelectric tensor has a dominant effect on the magnitude and spatial distribution of the flexoelectric polarization. In contrast, no polarization is induced around the b= a(010) screw dislocation through either electrostriction or flexoelectricity. Finally, our findings provide an in-depth understanding of the role of flexoelectricity in inducing polarization around dislocation cores and offer insights into the defect engineering of dielectric/ferroelectric materials.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Relativistic Density Functional NMR Tensors Analyzed with Spin–free Localized Molecular Orbitals

The implementation of fast relativistic methods based on density functional theory, in conjunction with localized molecular orbital (LMO) based analysis, allows straightforward interpretations of NMR parameters in terms of contributions from core shells, lone pairs, and bonds, for compounds containing elements from across the periodic table. We present a conceptual review of a frequently used LMO analysis of NMR parameters calculated in the presence of spin-orbit interactions and other relativistic effects. An accompanying example focuses on the 15 N shielding in a heavy metal complex.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Graph Contractions for Calculating Correlation Functions in Lattice QCD

Computing correlation functions for many-particle systems in Lattice QCD is vital to extract nuclear physics observables like the energy spectrum of hadrons such as protons. However, this type of calculation has long been considered to be very challenging and computing-resource intensive because of the complex nature of a hadron composed of quarks with many degrees of freedom. In particular, a correlation function can be calculated through a sum of all possible pairs of quark contractions, each of which is a batched tensor contraction, dictated by Wick's theorem. Because the number of terms of this sum can be very large for any hadronic system of interest, fast evaluation of the sum faces several challenges: an extremely large number of contractions, a huge memory footprint at runtime, and the speed of tensor contractions. In this paper, we present a Lattice QCD analysis software suite, Redstar, which addresses these challenges by utilizing novel algorithmic and software engineering methods targeting modern computing platforms such as many-core CPUs and GPUs. In particular, Redstar represents every term in the sum of a correlation function by a graph, applies efficient graph algorithms to reduce the number of contractions to lower the cost of computations, and minimizes the total memory footprint. Moreover, Redstar carries out the contractions on either CPUs or GPUs utilizing an internal and highly efficient Hadron contraction library. Specifically, we illustrate some important algorithmic optimizations of Redstar, show various key design features of Hadron library, and present the speedup values due to the optimizations along with performance figures for calculating six correlations functions on four computing platforms.

Chen, Jie↗

DOE OSTI2023

Graph Contractions for Calculating Correlation Functions in Lattice QCD

Chen, Jie↗

DOE OSTI2024

NEML2: A High Performance Library for Constitutive Modeling

NEML2, the New Engineering Material model Library, version 2, is an offshoot of NEML, an earlier material modeling code developed at Argonne National Laboratory. NEML2 extends the key philosophy of its predecessor, i.e., material models are flexible, modular, and can be built from smaller blocks. It also provides modern features that do not exist in the framework of its predecessor such as material model vectorization, automatic differentiation, device-portable just-in-time compilation, operator fusion, lazy tensor evaluation, etc. Moreover, NEML2 can seamlessly integrate with the popular machine learning package PyTorch to take advantage of modern and fast-growing machine learning techniques. In this fiscal year, the development of core library features and capabilities are complete. The purpose of this report is not to serve as a verbatim copy of the software API reference (which is available online at https://reverendbedford.github.io/neml2/). Instead, this report documents the motivation, implementation, design choices, and usage of each core capability as well as their applications in solving practical engineering problems. This report is compiled based on the NEML2 major release 2.0.0.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Violation of energy conditions and the gravitational radius of the proton

The energy-momentum tensor (EMT) of the proton encodes fundamental information about its mass, pressure, and shear distributions. Using recent lattice QCD data for the gravitational form factors, we show that the Breit-frame Wigner EMT may be of Hawking-Ellis type IV in the proton’s core. Such EMT violates all pointwise energy conditions and lacks a causal rest frame so that the usual mechanical picture fails at short distances. We define the —a new hadronic observable—marking the scale where the EMT becomes ordinary (type I) and the classical interpretation is restored. We also derive from the averaged null energy condition nonperturbative, model-independent quantum field theory constraints on gravitational form factors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Accelerated Constrained Sparse Tensor Factorization on Massively Parallel Architectures

This study presents the first constrained sparse tensor factorization (cSTF) framework that optimizes and fully offloads computation to massively parallel GPU architectures, and the first performance characterization of cSTF on GPU architectures. In contrast to prior work on tensor factorization, where the matricized tensor times Khatri-Rao product (MTTKRP) is the primary performance bottleneck, our systematic analysis of the cSTF algorithm on GPUs reveals that adding constraints creates an additional bottleneck in the update operation for many real-world sparse tensors. While executing the update operation on the GPU brings significant speedup over its CPU counterpart, it remains a significant bottleneck. To further accelerate the update operation, we propose cuADMM, a new update algorithm that leverages algorithmic and code optimization strategies to minimize both computation and data movement on GPUs. As a result, our framework delivers significantly improved performance compared to prior state-of-the-art. On 10 real-world sparse tensors, our framework achieves geometric mean speedup of 5.1 × (max 41.59 ×) and 7.01 × (max 58.05 ×) on the NIVIDA A100 and H100 GPUs, respectively, over the state-of-the-art SPLATT library running on a 26-core Intel Ice Lake Xeon CPU.

Soh, Yongseok↗

NASA NTRS2003

Spectral Characteristics of Wake Vortex Sound During Roll-Up

This report presents an analysis of the sound spectra generated by a trailing aircraft vortex during its rolling-up process. The study demonstrates that a rolling-up vortex could produce low frequency (less than 100 Hz) sound with very high intensity (60 dB above threshold of human hearing) at a distance of 200 ft from the vortex core. The spectrum then drops o rapidly thereafter. A rigorous analytical approach has been adopted in this report to derive the spectrum of vortex sound. First, the sound pressure was solved from an alternative treatment of the Lighthill s acoustic analogy approach [1]. After the application of Green s function for free space, a tensor analysis was applied to permit the removal of the source term singularity of the wave equation in the far field. Consequently, the sound pressure is expressed in terms of the retarded time that indicates the time history and spacial distribution of the sound source. The Fourier transformation is then applied to the sound pressure to compute its spectrum. As a result, the Fourier transformation greatly simplifies the expression of the vortex sound pressure involving the retarded time, so that the numerical computation is applicable with ease for axisymmetric line vortices during the rolling-up process. The vortex model assumes that the vortex circulation is proportional to the time and the core radius is a constant. In addition, the velocity profile is assumed to be self-similar along the aircraft flight path, so that a benchmark vortex velocity profile can be devised to obtain a closed form solution, which is then used to validate the numerical calculations for other more realistic vortex profiles for which no closed form solutions are available. The study suggests that acoustic sensors operating at low frequency band could be profitably deployed for detecting the vortex sound during the rolling-up process.

Booth, Earl R., Jr.↗

NASA NTRS1992

Cold dissipationless collapse of spherical systems - Sensitivity to the initial density law

The collapse of cold, initially spherical systems with varying degrees of central condensation is investigated. The way in which the final shape of a collapsing system depends on the initial density law is examined. For an initial stellar number density rho varies as r exp -n, where n is in the range 0-2.5, the final, nearly prolate shape is given by a/c is approximately equal to 1.28(1 + 0.16 n), where a/c is the ratio of long to short axes of the inertia ellipsoid computed from the moment of inertia tensor of the most tightly bound 80 percent of the mass. The properties associated with the final states in the present computations are also studied. The collapsing systems develop an anisotropic halo dominated by radial orbits surrounding an isotropic core as predicted by Burkert (1990).

Cannizzo, John K.↗

DOE OSTI2022

A multiphysics model of the versatile test reactor based on the MOOSE framework

The traditional modeling approach for sodium fast reactor cores relies on separate physics models, where the fuel performance, thermal–hydraulics, and neutronics calculations required to predict the core physics characteristics for nominal conditions are decoupled by relying on user-imposed boundary conditions. Here, this paper aims at evaluating the impact of multiphysics simulations for predicting the core characteristics of the Versatile Test Reactor, which is being designed as a 300-MWt sodium-cooled fast reactor. The purpose of the Versatile Test Reactor is to accelerate the testing of advanced nuclear materials in the United States. The proposed multiphysics model relies on the Griffin reactor physics code, the SAM thermal–hydraulic system code, the BISON fuel performance code, as well as generic Multiphysics Object-Oriented Simulation Environment capabilities implemented in the open-source tensor mechanics module. For k eff calculations, the introduction of a tight coupling between the neutronics, thermo-mechanical and thermal–hydraulics models induces a change of around 543 pcm in the eigenvalue, compared to the traditional standalone neutronics calculation where approximate temperature profiles are used. The multiphysics model is then employed for quantifying the impact of the thermal conductivity uncertainties on some of the key figures of merit, such as the fuel centerline temperature, assembly powers, and keff for nominal core conditions. As anticipated, uncertainties on fuel thermal conductivity mostly impact the fuel centerline temperature, and to a lesser extend the k eff .

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Slowly rotating gravastars

We solve Einstein’s equations for slowly rotating gravitational condensate stars (gravastars) up to second order in the rotation by expanding about the spherically symmetric gravastar with de Sitter interior and Schwarzschild exterior matched at their common horizon. Requiring that the perturbations are finite on the null surface reduces the exterior geometry to that of a Kerr black hole, implying that a slowly rotating gravastar cannot be distinguished from a Kerr black hole by any measurement or observation restricted to the macroscopic spacetime exterior to the horizon. Here we determine the interior solution, the surface stress tensor, and the Komar mass and angular momentum localized on the slowly rotating horizon surface. With the interior equation of state fixed at p = –ρ, finite junction conditions on the null horizon surface necessarily lead to an interior solution with a singular core, where the perturbative expansion breaks down. Comparison to other models and implications for more rapidly rotating gravastars are briefly discussed.

79 ASTRONOMY AND ASTROPHYSICS↗