Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tensor”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

MOCMC: Method of Characteristics Moment Closure, a Numerical Method for Covariant Radiation Magnetohydrodynamics

In this work, we present a conservative numerical method for radiation magnetohydrodynamics with frequency-dependent full transport in stationary spacetimes. This method is stable and accurate for both large and small optical depths and radiation pressures. The radiation stress–energy tensor is evolved in flux-conservative form, and closed with a swarm of samples that each transport a multigroup representation of the invariant specific intensity along a null geodesic. In each zone, the enclosed samples are used to efficiently construct a Delaunay triangulation of the unit sphere in the comoving frame, which in turn is used to calculate the Eddington tensor, average source terms, and adaptively refine the sample swarm. Furthermore, radiation four-forces are evaluated in the moment sector in a semi-implicit fashion. The radiative transfer equation is solved in invariant form deterministically for each sample. Since each sample carries a discrete representation of the full spectrum, the cost of evaluating the transport operator is independent of the number of frequency groups, representing a significant reduction of algorithmic complexity for transport in frequency-dependent problems. The major approximation we make in this work is performing scattering in an angle-averaged way. Local adaptivity in samples also makes this scheme more amenable to nonuniform meshes than a traditional Monte Carlo method. We describe the method and present results on a suite of test problems. We find that Method of Characteristics Moment Closure converges at least as ~N -1 , rather than the canonical Monte Carlo N -1/2 , where N is the number of samples per zone.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

CMB-S4: Iterative Internal Delensing and r Constraints

Abstract The tightest constraints on the tensor-to-scalar ratio r can only be obtained after removing a substantial fraction of the lensing B -mode sample variance. The planned Cosmic Microwave Background (CMB)-S4 experiment ( cmb-s4.org ) will remove the lensing B -mode signal internally by reconstructing the gravitational lenses from high-resolution observations. We document here a first lensing reconstruction pipeline able to achieve this optimally for arbitrary sky coverage. We make it part of a map-based framework to test CMB-S4 delensing performance and its constraining power on r , including inhomogeneous noise and two non-Gaussian Galactic polarized foreground models. The framework performs component separation of the high-resolution maps, followed by the construction of lensing B -mode templates, which are then included in a parametric small-aperture map cross-spectra-based likelihood for r . We find that the lensing reconstruction and framework achieve the expected performance, compatible with the target σ ( r ) ≃ 5 · 10 −4 in the absence of a tensor signal, after an effective removal of 92%–93% of the lensing B -mode variance, depending on the simulation set. The code for the lensing reconstruction can also be used for cross-correlation studies with large-scale structures, lensing spectrum reconstruction, cluster lensing, or other CMB lensing-related purposes. As part of our tests, we also demonstrate the joint optimal reconstruction of the lensing potential with the lensing curl potential mode at second order in the density fluctuations.

Astronomy & Astrophysics↗
DOE OSTI2026

Velocity-space Origins of the Pressure–Strain Interaction in Multipopulation Distributions and Its Application to Magnetic Reconnection

A forefront research question is how energy evolves in weakly collisional plasmas for which departures from local thermodynamic equilibrium (LTE) are significant. The standard approach is studying the terms in the non-LTE energy evolution equation derived by taking the second moment of the Boltzmann equation, but the resultant fluid metrics do not retain information about which particles at which velocities drive energy evolution. A widely studied channel for internal energy density evolution is the pressure–strain interaction. Here, we employ the kinetic pressure–strain, a phase-space diagnostic whose velocity-space integral recovers the pressure–strain interaction to disambiguate the contributions to the pressure–strain interaction from disparate particle populations in composite phase-space densities. We develop phase-space analogs of the pressure–strain interaction decompositions to provide the phase-space origins of normal versus sheared flow. We introduce the “kinetic strain-rate” tensor, the phase-space analog of the strain-rate tensor, which we argue is needed to interpret the phase-space origins of the pressure–strain interaction. To demonstrate the utility of these quantities, we investigate them for composite electron distributions near the electron diffusion region in two-dimensional particle-in-cell simulations of antiparallel symmetric magnetic reconnection. We find that the phase-space-based diagnostics isolate the roles of distinct populations. These results contribute to a growing body of work providing new methods for quantifying phase-space energy evolution for a broad array of processes, from magnetic reconnection to collisionless shocks and turbulence, opening new pathways for answering longstanding problems of particle energization in weakly collisional plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Short-range correlations in nuclei

Atomic nuclei are held together by the strong nuclear force acting between protons and neutrons (nucleons). While the long range, averaged part of this force is well described by the nuclear shell model, the short-range and tensor components create a fascinating substructure: pairs of nucleons that momentarily approach each other very closely, acquiring large relative momenta. These short-range correlated (SRC) pairs account for roughly 20% of all nucleons in any nucleus and almost all of the high momentum nucleons. This chapter provides an introduction to SRC pairs: their origin in the nucleon-nucleon tensor force, the experimental methods used to study them, principally deep inelastic and quasielastic electron and proton scattering, and the comprehensive picture that has emerged over the past three decades.

FOS: Physical sciences↗
DOE OSTI2026

Anti-Ultralocality and Plateau Models of Inflation

Anti-ultralocality refers to the growth of spatial gradient terms relative to velocity terms in the coupled Einstein--scalar field equations. It is a characteristic feature of decelerated expansion before the onset of inflation. Previous numerical relativity studies have shown that anti-ultralocality prevents the onset of inflation in models with power-law inflaton potentials. In this paper, we show that models with plateau-shaped inflaton potentials, which are considered to be the simplest way to generate a tensor-to-scalar ratio below current observational upper limits, are especially vulnerable to anti-ultralocality effects. The reasons are the flatness of the plateau and the energy density gap of $\sim 10$ orders of magnitude between the Planck density and the plateau potential energy. To study the problem, we develop a protocol for assessing the viability of inflationary models in general, and we apply it to a plateau potential using a previously validated numerical relativity code. We find that, starting from generic initial conditions, the growth of gradient terms in the Einstein equations relative to non-gradient terms either prevents inflation from lasting for enough $e$-folds or triggers a phase of quantum runaway. We show that the fine-tuning of initial conditions necessary to avoid these issues becomes more severe as the energy scale of inflation is made smaller, disfavoring common approaches for reducing the tensor-to-scalar ratio.

FOS: Physical sciences↗
DOE OSTI2023

Semi-classical Kinetic Theory for Massive Spin-half Fermions with Leading-order Spin Effects

We consider the quantum kinetic-theory description for interacting massive spin-half fermions using the Wigner function formalism. We derive a general kinetic theory description assuming that the spin effects appear at the classical and quantum level. To track the effect of such different contributions we use the semi-classical expansion method to obtain the generalized dynamical equations including spin, analogous to classical Boltzmann equation. This approach can be used to obtain a collision kernel involving local as well as non-local collisions among the microscopic constituent of the system and eventually, a framework of spin hydrodynamics ensuring the conservation of the energy-momentum tensor and total angular momentum tensor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Scattering of lower hybrid radio frequency waves by cylindrical turbulent structures in the plasma edge in tokamaks

Lower hybrid (LH) radio frequency (RF) waves are used in fusion tokamak devices to generate non-inductively toroidal currents. LH waves are effective in imparting toroidal momentum to electrons in the core of the confined plasma. Since the LH waves are generated by wave guide structures near the wall of a tokamak, the waves have to propagate through a turbulent plasma in the edge region before coupling to the core. It is especially important to quantify the effect of this plasma region on the propagation characteristics of the LH waves, since fusion reactors like ITER will have an extended edge region. Structures in the plasma edge like filaments and blobs which have highly varying density fluctuations compared to the background density, span radial spatial scales that are comparable to the LH wavelength. We study the scattering of the LH waves by the filamentary structures using Maxwell's equations, in which the plasma permittivity is given by the cold plasma dispersion tensor. The collisional absorption of LH waves in the edge region is included by a modification to the elements of the dispersion tensor. The filaments are assumed to be cylindrical with the axis predominantly aligned along the direction of the toroidal magnetic field. Our studies are both analytical and numerical [1,2] and show that these structures can lead to reflection, refraction, diffraction, and side-scattering of both an incident LH plane wave and a Gaussian beam. We will present a variety of different density variations and collisional absorption rates in plasmas with filamentary structures of varying sizes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

The Sylvester Graphical Lasso (SyGlasso)

This paper introduces the Sylvester graphical lasso (SyGlasso) that captures multiway dependencies present in tensor-valued data. The model is based on the Sylvester equation that defines a generative model. The proposed model complements the tensor graphical lasso that imposes a Kronecker sum model for the inverse covariance matrix by providing an alternative Kronecker sum model that is generative and interpretable. A nodewise regression approach is adopted for estimating the conditional independence relationships among variables. The statistical convergence of the method is established, and empirical studies are provided to demonstrate the recovery of meaningful conditional dependency graphs. Here we apply the SyGlasso to an electroencephalography (EEG) study to compare the brain connectivity of alcoholic and nonalcoholic subjects. Here, we demonstrate that our model can simultaneously estimate both the brain connectivity and its temporal dependencies.

Wang, Yu↗
DOE OSTI2020

Scalable Heterogeneous Execution of a Coupled-Cluster Model with Perturbative Triples

The CCSD(T) coupled-cluster model with perturbative triples is considered a gold standard for computational modeling of the correlated behavior of electrons in molecular systems. A fundamental constraint is the relatively small global-memory capacity in GPUs compared to the main-memory capacity on host nodes, necessitating relatively smaller tile sizes for high-dimensional tensor contractions in NWChem's GPU-accelerated implementation of the CCSD(T) method. A coordinated redesign is described to address this limitation and associated data movement overheads, including a novel fused GPU kernel for a set of tensor contractions, along with inter-node communication optimization and data caching. The new implementation of GPU-accelerated CCSD(T) improves overall performance by 3.4x. Finally, we discuss the trade-offs in using this fused algorithm on current and future supercomputing platforms.

Kim, Jinsung↗
DOE OSTI2020

Theoretical and Phenomenological Studies of the Nucleon Structure in High Energy Reactions

The understanding of the internal structure of the proton and other strongly interacting particles is at the forefront of modern nuclear physics research. Generalized Parton Distribution Functions (GPDs) are a powerful tool to advance the understanding of hadron structure. In addition to the information about the one- dimensional collinear momentum distributions of partons (quarks, anti-quarks, and gluons) known from studies of high energy deep-inelastic reactions, GPDs also carry information on the distribution of partons in the transverse plane, and allow us in this way to access the three-dimensional structure of the nucleon. GPDs can be studied in hard-exclusive reactions and contain also information on the energy-momentum tensor form factors which will allow us to gain insights on quantities like pressure or angular momentum distribution inside the nucleon. The goal of this thesis is to deepen our understanding of the three-dimensional structure of the nucleon. We investigate energy-momentum tensor form factors and densities, and all leading-twist GPDs in the bag model. This quark model provides a consistent theoretical framework to investigate many general concepts that have recently attracted interest, and allows one to study insightful limits like the large-Nc limit, heavy-quark limit or the non-relativistic limit. Another important aspect of this thesis is the discussion of the monopole and quadrupole contributions to the angular momentum density. Finally, the description of pseudoscalar meson production in exclusive processes in terms of GPDs is discussed in the Goloskokov-Kroll model and implemented in the PARTONS framework, a software development project which will provide direct support for experiments at the Jefferson National Lab and the future Electron-Ion Collider.

Tezgin, Kemal↗
DOE OSTI2024

Probing the Isospin Composition of Short-Range Correlated Pairs at Jefferson Lab Hall B

Nucleons in short-range correlated (SRC) pairs, due to their close proximity and high relative momentum, can provide insight into the short-range part of the strong nuclear interaction. In particular, the prevalence of np pairs is due to the dominance of a tensor term for correlated nucleons with momenta of approximately 400?600 MeV/c. This dissertation comprises two studies advancing the community?s understanding of the isospin composition of SRC pairs. First, I performed a study of proton and neutron knockout from initially low-momentum and high-momentum states in 3He. Previous work has shown that protons are disproportionately represented in high-momentum states in neutron-rich nuclei. I demonstrate that spectral functions for the proton-rich nucleus 3He predict, in agreement with data, that neutrons are disproportionately represented in high-momentum states, but that 3He does not display the same strong prevalence of np pairs that is observed in larger nuclei. Second, Generalized Contact Formalism (GCF), a well-supported theory for predicting SRC behavior, predicts the transition from an isospin-dependent, tensor-dominant interaction at intermediate distances to a scalar-dominant, isospin-independent interaction at very short distances. This dissertation uses data from the CLAS12 Nuclear Targets Experiments in Hall B at Jefferson Lab to measure the relative abundances of pp and pn pairs for increasing relative momentum and decreasing separation. I provide an independent confirmation of the previously-observed increase in pp pairs at increasing momentum of the struck nucleon. I also contribute to the application of the new CLAS12 Central Neutron Detector by precisely measuring the neutron detection efficiency and developing a machine learning model for rejecting charged particle background.

Seroka, Erin↗
DOE OSTI2023

Utah FORGE 3-2514: A Strain Sensing Array to Characterize Deformation at the FORGE Site - Workshop Presentation

This is a presentation on the Strain Sensing Array to Characterize Deformation at the FORGE Site project by Clemson University, presented by Lawrence Murdoch. The project's objective was to evaluate the feasibility of measuring and interpreting tensor strain data to improve the performance of EGS. This presentation was featured in the Utah FORGE R&D Annual Workshop on September 8, 2023. The workshop provided a valuable opportunity to explore the progress made in each of the 17 Research and Development projects funded under Solicitation 2020-1 which aim to enhance our understanding of the crucial factors influencing the development of Enhanced Geothermal Systems (EGS) reservoirs and resources.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Direct interpolative construction of the discrete Fourier transform as a matrix product operator

The quantum Fourier transform (QFT), which can be viewed as a reindexing of the discrete Fourier transform (DFT), has been shown to be compressible as a low-rank matrix product operator (MPO) or quantized tensor train (QTT) operator. However, the original proof of this fact does not furnish a construction of the MPO with a guaranteed error bound. Meanwhile, the existing practical construction of this MPO, based on the compression of a quantum circuit, is not as efficient as possible. We present a simple closed-form construction of the QFT MPO using the interpolative decomposition, with guaranteed near-optimal compression error for a given rank. This construction can speed up the application of the QFT and the DFT, respectively, in quantum circuit simulations and QTT applications. We also connect our interpolative construction to the approximate quantum Fourier transform (AQFT) by demonstrating that the AQFT can be viewed as an MPO constructed using a different interpolation scheme.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Utah FORGE 3-2514: Strain Sensing Array for Deformation Characterization - 2024 Annual Workshop Presentation

This is a presentation on the Strain Sensing Array to characterize deformation at the FORGE site project by Clemson University, presented by Lawrence Murdoch. The project's objective was to evaluate the feasibility of measuring and interpreting tensor strain data to improve the performance of EGS. This presentation was featured in the Utah FORGE R&D Annual Workshop on August 23, 2024.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

A Low-Rank QTT-based Finite Element Method for Elasticity Problems

We present an efficient and robust numerical algorithm for solving the linear elasticity problem that combines the Quantized Tensor Train format and a domain partitioning strategy. This approach makes it possible to solve the linear elasticity problem on a computational domain that is more general than a square. By integrating Z-ordering and subdomain concatenation, our method substantially decreases memory usage and achieves a notable reduction in rank compared to established Finite Element implementations like the FEniCS platform. This efficiency is maintained while still guaranteeing exponential convergence with respect to the number of degrees of freedom. This performance gain, however, requires a fundamental rethinking of how core finite element operations are implemented. This includes changes to mesh discretization, node and degree of freedom ordering, stiffness matrix and internal nodal force assembly, and the execution of algebraic matrix-vector operations. In this work, we discuss all these aspects in detail and assess the method’s performance in the numerical approximation of three representative test cases.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Comparison of DeePMD, MTP, GAP, ACE and MACE Machine‐Learned Potentials for Radiation‐Damage Simulations: A User Perspective

Accurate and efficient interatomic potentials are essential for molecular dynamics (MD) simulations of radiation damage, gas diffusion, and phase stability in complex ceramics such as LiAlO 2 , especially under extreme conditions relevant to tritium production. Here, we evaluate the performance of six machine-learned interatomic potentials (MLIPs), moment tensor potential (MTP), Gaussian approximation potential, deep potential (DeePMD), atomic cluster expansion (ACE), message-passing ACE (multilayer atomic cluster expansion (MACE) pretrained) and MACE (trained from-scratch), all trained on the same density functional theory dataset with inclusion of tritium. The MLIPs are benchmarked against traditional Buckingham and ReaxFF potentials in terms of energy accuracy, density predictions, thermal equilibration behavior, threshold displacement energy (E d ), tritium diffusivity, and computational cost. Among the models, MTP shows the best overall balance between efficiency and accuracy, with low force and energy errors and realistic E d values for Li and Al. The ACE and MACE (pretrained and trained from scratch) models exhibit high E d (>200 eV) and unphysical pair interactions. DeePMD underestimates Ed due to overly repulsive behavior even at equilibrium distances. All models over-estimate tritium diffusion but the pretrained MACE model behaves well during tritium-diffusion simulations up to 500 K, maintaining diffusivities in the physically consistent 10 −11 m 2 /s range. Finally, we quantify the computational cost of each potential in large-scale atomic/molecular massively parallel simulator, finding that only MTP is more efficient than traditional empirical potentials, while others are significantly more expensive. These findings explain the trade-offs between accuracy and computational cost in MLIP development and provide essential guidance for use in high-throughput radiation damage and gas diffusion simulations in nuclear ceramics.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Nanoscale lattice strains in self-ion implanted tungsten

Developing a comprehensive understanding of the modification of material properties by neutron irradiation is important for the design of future fission and fusion power reactors. Self-ion implantation is commonly used to mimic neutron irradiation damage, however an interesting question concerns the effect of ion energy on the resulting damage structures. The reduction in the thickness of the implanted layer as the implantation energy is reduced results in the significant quandary: Does one attempt to match the primary knock-on atom energy produced during neutron irradiation or implant at a much higher energy, such that a thicker damage layer is produced? In this work, we address this question by measuring the full strain tensor for two ion implantation energies, 2 MeV and 20 MeV in self-ion implanted tungsten, a critical material for the first wall and divertor of fusion reactors. A comparison of 2 MeV and 20 MeV implanted samples is shown to result in similar lattice swelling. Multi-reflection Bragg coherent diffractive imaging (MBCDI) shows that implantation induced strain is in fact heterogeneous at the nanoscale, suggesting that there is a non-uniform distribution of defects, an observation that is not fully captured by micro-beam Laue diffraction. At the surface, MBCDI and high-resolution electron back-scattered diffraction (HR-EBSD) strain measurements agree quite well in terms of this clustering/non-uniformity of the strain distribution. However, MBCDI reveals that the heterogeneity at greater depths in the sample is much larger than at the surface. This combination of techniques provides a powerful method for detailed investigation of the microstructural damage caused by ion bombardment, and more generally of strain related phenomena in micro-volumes that are inaccessible via any other technique.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Optimized Lie–Trotter–Suzuki decompositions for two and three non-commuting terms

Lie–Trotter–Suzuki decompositions are an efficient way to approximate operator exponentials exp ( t H ) when H is a sum of n (non-commuting) terms which, individually, can be exponentiated easily. They are employed in time-evolution algorithms for tensor network states, digital quantum simulation protocols, path integral methods like quantum Monte Carlo, and splitting methods for symplectic integrators in classical Hamiltonian systems. Here, we provide optimized decompositions up to order t 6 . The leading error term is expanded in nested commutators (Hall bases) and we minimize the 1-norm of the coefficients. For n = 2 terms, several of the optima we find are close to those in McLachlan (1995). Generally, our results substantially improve over unoptimized decompositions by Forest, Ruth, Yoshida, and Suzuki. We explain why these decompositions are sufficient to efficiently simulate any one- or two-dimensional lattice model with finite-range interactions. This follows by solving a partitioning problem for the interaction graph.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗