Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tensor completion”

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 109 records · Page 6

Structure and energetics of the elbows in the Au(111) herringbone reconstruction

We study the structure of the threading edge dislocations, or “elbows,” which are an essential component of the well-known herringbone reconstruction of the (111) surface of Au. Previous work had shown that these dislocations can be stabilized by long-range elastic relaxations into the bulk. However, the validity of the harmonic spring model that had been used to estimate the energies of the dislocations is uncertain. To enable a more refined model of the dislocation energetics, we have imaged the atomic structure of these dislocations using scanning tunneling microscopy. We find that the harmonic spring model does not adequately reproduce the observed structure. We are able to reproduce the structure, however, with a two-dimensional Frenkel-Kontorova (FK) model that uses a pairwise Morse potential to describe the interactions between the top layer Au atoms on a rigid substrate. The parameters of the potential were obtained by fitting the energy of uniaxially compressed phases, or “stripes”, computed with density functional theory, as a function of surface Au density. Within this model, the formation of the threading dislocations remains unfavorable. However, the large forces on the substrate atoms near the threading-dislocation cores, render the assumption of a completely rigid substrate questionable. Indeed, if the FK parameters are modified to account for the relaxation of just one more atomic layer, threading dislocations can, in principle, become favorable, even without bulk elastic relaxations. Additional evidence for a small elbow energy is that our computed change in the Au(111) surface stress tensor caused by the (√ 3 × 22) reconstruction is considerably smaller than previous estimates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A self-documenting source-independent data format for computer processing of tensor time series

The UCLA Space Science Group has developed a fixed format intermediate data set called a block data set, which is designed to hold multiple segments of multicomponent sampled data series. The format is sufficiently general so that tensor functions of one or more independent variables can be stored in the form of virtual data. This makes it possible for the unit data records of the block data set to be arrays of a single dependent variable rather than discrete samples. The format is self-documenting with parameter, label and header records completely characterizing the contents of the file. The block data set has been applied to the filing of satellite data (of ATS-6 among others).

Mcpherron, R. L.↗

Spin Decoherence in VOPc@graphene Nanoribbon Complexes

Carbon nanoribbons or nanographene qubit arrays can facilitate quantum-to-quantum transduction between light, charge, and spin, making them an excellent testbed for fundamental science in quantum coherent systems and for the construction of higher-level qubit circuits. Here, in this work, we study spin decoherence due to coupling with a surrounding nuclear spin bath of an electronic molecular spin of a vanadyl phthalocyanine (VOPc) molecule integrated onto an armchair-edged graphene nanoribbon (GNR). Density functional theory (DFT) is used to obtain ground-state atomic configurations. Decay of spin coherence in Hahn echo experiments is then simulated using the cluster correlation expansion method with a spin Hamiltonian involving hyperfine and electric field gradient tensors calculated from DFT. We find that the decoherence time T 2 is anisotropic with respect to magnetic field orientation and determined only by the hydrogen nuclear spins on both VOPc and GNR. Large electron spin echo envelope modulation (ESEEM) due to nitrogen and vanadium nuclear spins is present at specific field ranges and can be completely suppressed by tuning the magnetic field. The relation between these field ranges and the hyperfine interactions is analyzed. The effects of interactions with the nuclear quadrupole moments are also studied, validating the applicability and limitations of the spin Hamiltonian when they are disregarded.

Hamiltonians↗

Proton’s gluon GPDs at large skewness and gravitational form factors from near threshold heavy quarkonium photoproduction

We study the exclusive near threshold photoproduction of heavy quarkonium in the framework of the generalized parton distribution (GPD) factorization, taking the 𝐽/𝜓 production as an example. Because of the threshold kinematics, the Compton-like amplitudes are related to gluon GPDs at large skewness 𝜉, distinct from the common kinematics in asymptotic high energy where the skewness is typically small. We discuss the nature of large-𝜉 expansion of these amplitudes in terms of the moments of gluon GPDs in the large-𝜉 limit. Based on that, we propose several ways to extract the first few moments of the gluon GPDs from these amplitudes, with the leading ones corresponding to the gluonic gravitational or energy-momentum tensor form factors (GFFs). We apply these methods to analyze the recent near threshold 𝐽/𝜓 production measurements by the 𝐽/𝜓 007 experiment and GlueX Collaboration and find that the 𝜉 scaling of the measured differential cross sections is consistent with the asymptotic behavior. However, the current data are not accurate enough yet for a complete determination of the gluonic GFFs, and therefore we consider some prospects for better extractions in the future.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Automated calculation and convergence of defect transport tensors

Defect diffusion is a key process in materials science and catalysis, but as migration mechanisms are often too complex to enumerate a priori, calculation of transport tensors typically have no measure of convergence and require significant end-user intervention. These two bottlenecks prevent high-throughput implementations essential to propagate model-form uncertainty from interatomic interactions to predictive simulations. In order to address these issues, we extend a massively parallel accelerated sampling scheme, autonomously controlled by Bayesian estimators of statewide sampling completeness, to build atomistic kinetic Monte Carlo models on a state-space irreducible under exchange and space group symmetries. Focusing on isolated defects, we derive analytic expressions for drift and diffusion coefficients, providing a convergence metric by calculating the Kullback–Leibler divergence across the ensemble of diffusion processes consistent with the sampling uncertainty. The autonomy and efficacy of the method is demonstrated on surface trimers in tungsten and Hexa-interstitials in magnesium oxide, both of which exhibit complex, correlated migration mechanisms.

36 MATERIALS SCIENCE↗

A Robust Event Diagnostics Platform: Integrating Tensor Analytics and Machine Learning into Real-time Grid Monitoring

The objective of this project is to develop a robust event diagnostics (RED) platform by integrating state-of-the-art tensor analytics and machine learning into real-time grid monitoring. The proposed platform can effectively analyze and discover the information hiding within the provided PMU data for effective real-time grid monitoring. The proposed RED platform provides a set of robust diagnostics tools for grid operation and management, including 1) data quality assessment, 2) data completion, 3) event detection, and 4) robust event classification. All the functionalities of the RED platform can help the operator to make informed decisions and respond in a timely manner. The developed RED platform will serve as an innovative advisory tool to reliably identify key events and discover new insights about the events and grid characteristics in the PMU data, and contribute to the efficient, safe, reliable operation and design of the nation’s electric system.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Microwave Scattering Model for Grass Blade Structures

In this paper, the electromagnetic scattering solution for a grass blade with complex cross-section geometry is considered. It is assumed that the blade cross section is electrically small, but its length is large compared to the incident wavelength. In a recent study it has been shown that the scattering solution for such problems, in the form of a polarizability tensor, can be obtained using the low-frequency approximation in conjunction with the method of moments. In addition, the study shows that the relationship between the polarizability tensor of a dielectric cylinder and its dielectric constant can be approximated by a simple algebraic expression. The results of this study are used to show that this algebraic approximation is valid also for cylinders with cross sections the shape of grass blades, providing that proper values am selected for each of three constants appearing in the expression. These constants are dependent on cylinder shape, and if the relationship between the constants and the three parameters describing a grass blade shape can be determined, an algebraic approximation relating polarizability tensor to blade shape, as well as dielectric constant, can be formed. Since the elements of the polarizability tensor are dependent on only these parameters, this algebraic approximation can replace the cumbersome method of moments model. A conjugate gradient method is then implemented to correctly determine the three constants of the algebraic approximation for each blade shape. A third-order polynomial fit to the data is then determined for each constant, thus providing a complete analytic replacement to the numerical (moment method) scattering model. Comparisons of this approximation to the numerical model show an average error of less than 3%.

Stiles, James M.↗

Proving the 6d Cardy formula and matching global gravitational anomalies

A Cardy formula for 6d superconformal field theories (SCFTs) conjectured by Di Pietro and Komargodski in [1] governs the universal behavior of the supersymmetric partition function on S^1_\beta \times S^5 S β 1 × S 5 in the limit of small \beta β and fixed squashing of the S^5 S 5 . For a general 6d SCFT, we study its 5d effective action, which is dominated by the supersymmetric completions of perturbatively gauge-invariant Chern-Simons terms in the small \beta β limit. Explicitly evaluating these supersymmetric completions gives the precise squashing dependence in the Cardy formula. For SCFTs with a pure Higgs branch (also known as very Higgsable SCFTs), we determine the Chern-Simons levels by explicitly going onto the Higgs branch and integrating out the Kaluza-Klein modes of the 6d fields on S^1_\beta S β 1 . We then discuss tensor branch flows, where an apparent mismatch between the formula in [1] and the free field answer requires an additional contribution from BPS strings. This ``missing contribution’’ is further sharpened by the relation between the fractional part of the Chern-Simons levels and the (mixed) global gravitational anomalies of the 6d SCFT. We also comment on the Cardy formula for 4d \mathcal{N}=2 𝒩 = 2 SCFTs in relation to Higgs branch and Coulomb branch flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Mechanical Erosion Modeling of TPS Materials

This work describes the development of a model that accounts for the additional total surface recession in Thermal Protection Systems materials as a result of mechanical erosion due to high shear conditions during atmospheric entry. A computational solid mechanics capability module was integrated within the Porous material Analysis Toolbox based on OpenFOAM, PATO. The mechanical erosion was modeled in three steps: first, the implemented stress analysis solver module computes the stress and the displacement fields for orthotropic materials using the wall shear stress tensor as a boundary condition; then, regions on the surface where the stress meets the failure criteria are identified; and last, the mesh is moved accordingly to remove the failed material. The outcome is a model of the total recession of the material due to surface chemistry and mechanical erosion. Results will be included in the final paper after verifying and completing the study.

Stress Analysis↗

Mechanical Erosion Modeling of TPS Materials

This work describes the development of a model that accounts for the additional surface recession in Thermal Protection Systems (TPS) materials as a result of mechanical erosion due to high shear conditions during atmospheric entry. A computational solid mechanics module was integrated within the Porous material Analysis Toolbox (PATO) based on OpenFOAM. The mechanical erosion was modeled in three steps: first, the implemented stress analysis solver computes the stress and the displacement fields for orthotropic materials using the wall shear stress tensor as a boundary condition; then, regions on the surface where the stress meets the failure criteria are identified; and last, the mesh is moved accordingly to remove the failed material. The outcome is a model capable to predict the total recession of the material due to surface chemistry and mechanical erosion. Results will be included in the final paper after verifying and completing the study.

Stress Analysis↗

Mechanical Erosion Modeling of TPS Materials

This work describes the development of a model that accounts for the additional surface recession in Thermal Protection Systems (TPS) materials as a result of mechanical erosion due to high shear conditions during atmospheric entry. A computational solid mechanics module was integrated within the Porous material Analysis Toolbox (PATO) based on OpenFOAM. The mechanical erosion was modeled in three steps: first, the implemented stress analysis solver computes the stress and the displacement fields for orthotropic materials using the wall shear stress tensor as a boundary condition; then, regions on the surface where the stress meets the failure criteria are identified; and last, the mesh is moved accordingly to remove the failed material. The outcome is a model capable to predict the total recession of the material due to surface chemistry and mechanical erosion. Results will be included in the final paper after verifying and completing the study.

Stress Analysis↗

Mechanical form factors and densities of nonrelativistic fermions

The hadron physics community has been actively debating the interpretation of so-called mechanical properties of hadrons. Nonrelativistic quantum-mechanical systems like the hydrogen atom have been appealed to in these debates as analogies. Since such appeals are likely to continue, it is important to have Galilei-covariant expressions for matrix elements of the energy-momentum tensor. In this work, I obtain Galilei-covariant breakdowns of such matrix elements into mechanical form factors, with a special focus on spin-half states. I additionally study the spatial densities associated with these form factors, using the pilot wave interpretation to guide their breakdown into contributions from internal structure and from quantum-mechanical effects such as wave packet dispersion. For completeness, I also obtain nonrelativistic Breit frame densities.

form factors↗

One- and two-dimensional higher-point conformal blocks as free-particle wavefunctions in $$ {\textrm{AdS}}_3^{\otimes m} $$

Abstract We establish that all of the one- and two-dimensional global conformal blocks are, up to some choice of prefactor, free-particle wavefunctions in tensor products of AdS 3 or limits thereof. Our first core observation is that the six-point comb-channel conformal blocks correspond to free-particle wavefunctions on an AdS 3 constructed directly in cross-ratio space. This construction generalizes to blocks for a special class of diagrams, which are determined as free-particle wavefunctions in tensor products of AdS 3 . Conformal blocks for all the remaining topologies are obtained as limits of the free wavefunctions mentioned above. Our results show directly that the integrable models associated with all one- and two-dimensional conformal blocks can be seen as limits of free theory, and manifest a relation between AdS and CFT kinematics that lies outside of the standard AdS/CFT dictionary. We complete the discussion by providing explicit Feynman-like rules that can be used to work out blocks for all topologies, as well as a Mathematica notebook that allows simple computation of Casimir equations and series expansions for blocks, by requiring just an OPE diagram as input.

Physics↗

Source Characterization at Low Yields from Physics Experiment One A

The detonation of a well-recorded coupled explosion at the well-characterized Nevada National Security Site, as part of Physics Experiment One, has provided us with an opportunity to compare and contrast methods and results of yield estimation and other event characterization methods. Yield estimates are made using several methods, including local/regional magnitudes, waveform envelopes, acoustic amplitudes, joint seismoacoustic yield, and moment tensors. We include in our analysis a discussion and comparison of applicability, calibrations, assumptions, and uncertainties involved in each. With the exception of magnitude-based methods, all the techniques were more or less successful in recovering the known yield. It is, however, useful to be cognizant of requirements and calibrations, such as station coverage, emplacement conditions, propagation characterization, and uncertainties, to apply the methods to other areas with less complete information. We demonstrate that we can be successful at monitoring low-yield explosions using a variety of methods assuming that the signal-to-noise ratio is adequate at recording stations, and that the material properties are taken into account. Depth estimates are also important, particularly if one allows the possibility of a surface explosion.

Pasyanos, Michael E. [Lawrence Livermore National ↗

Impact of tensor couplings with scalar mixing on covariant energy density functionals

The recent pioneering campaigns conducted by the Lead Radius Experiment (PREX) and the Calcium Radius Experiment (CREX) Collaborations have uncovered major deficiencies in the theoretical description of some fundamental properties of atomic nuclei. Following a recent refinement to the isovector sector of covariant energy density functionals we present here additional improvements to the functional by including both tensor couplings and an isoscalar-isovector mixing term in the scalar sector. Motivated by the distinct surface properties of calcium and lead, we expect that the tensor terms that generate derivative couplings will help break the linear correlation between the neutron skin thickness of these two nuclei. Moreover, the addition of these new terms mitigates most of the problems identified by Reed et al. in describing the properties of both finite nuclei and neutron stars. While significant progress has been made in reconciling the PREX-CREX results without compromising other observables, the final resolution awaits the completion of a proper calibration for this new class of functionals. As a result, we expect that powerful reduced basis methods used recently to create efficient emulators will be essential to accomplish this task.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling and simulation of a turbulent far wake

Work continued on two projects which had been started during previous years. Both projects involve calculations of the subsonic, turbulent far wake of a two-dimensional object at a Reynolds number of 1000 (based on wake momentum thickness). This flow was used as a test case for direct comparison of various turbulence models and a direct numerical simulations (DNS) of this flow were undertaken. In the turbulence model comparison studies, for any particular model tested, a unique self-similar solution was obtained far enough downstream, regardless of inlet conditions. Furthermore, different turbulence models led to different far-wake equilibrium solutions. No turbulence model could correctly predict all features of the turbulent far wake. For example, the spreading rate and turbulent shear stresses were underpredicted by all the standard models (both two-equation and full Reynolds stress models). In cases where a more correct spreading rate was achieved, it was at the expense of the turbulent kinetic energy, which was overpredicted. In general, the Algebraic Dissipation Rate Model of Gatski and Speziale, 1992, when added to any of the standard models, improved the results dramatically. Also, full Reynolds stress closure models did a much better job at predicting the shapes of both the mean and turbulence profiles, but the spreading rate was not significantly improved over that predicted by the simpler two-equation models. There are two main conclusions from these studies: First, in a comparison such as this, it is not enough to compare just one parameter, like the spreading rate. A good prediction for one parameter does not necessarily imply good predictions for all parameters in a flow. Second, since no turbulence model could correctly predict the turbulent far wake, much of the important physics of turbulent free shear flows is apparently lost by the assumptions inherent in today's methods of turbulence modeling; turbulence models must be improved. Direct simulations of this flow were begun last year in order to provide a data base through which some of the deficiencies of the existing turbulence models could be identified. Quantities such as the pressure-strain correlation, turbulent diffusion, and the dissipation rate tensor can be easily calculated from the DNS results, whereas these quantities are nearly impossible to measure experimentally. Improvements to existing turbulence models (and development of new models) require knowledge about flow quantities such as these. During this summer, diagnostics codes were written which will calculate the parameters mentioned above, along with other single-point and multi-point statistics. The DNS calculations are still in progress at the time of this writing. When these calculations are complete, the diagnostics codes will be applied so that the results can aid turbulence modelers. In addition, the results will show whether or not there exists a universal equilibrium turbulent far wake, independent of initial conditions.

Cimbala, John M.↗

Uncertainty Quantification of Geophysical Inversion Using Stochastic Partial Differential Equations (LDRD #218329)

This report summarizes work completed under the Laboratory Directed Research and Development (LDRD) project "Uncertainty Quantification of Geophysical Inversion Using Stochastic Differential Equations." Geophysical inversions often require computationally expensive algorithms to find even one solution, let alone propagating uncertainties through to the solution domain. The primary purpose of this project was to find more computationally efficient means to approximate solution uncertainty in geophysical inversions. We found multiple computationally efficient methods of propagating Earth model uncertainty into uncertainties in solutions of full waveform seismic moment tensor inversions. However, the optimum method of approximating the uncertainty in these seismic source solutions was to use the Karhunen-Love theorem with data misfit residuals. This method was orders of magnitude more computationally efficient than traditional Monte Carlo methods and yielded estimates of uncertainty that closely approximated those of Monte Carlo. We will summarize the various methods we evaluated for estimating uncertainty in seismic source inversions as well as work toward this goal in the realm of 3-D seismic tomographic inversion uncertainty.

58 GEOSCIENCES↗

Loop-string-hadron approach to SU(3) lattice Yang-Mills theory: Hilbert space of a trivalent vertex

The construction of gauge-invariant states of SU(3) lattice gauge theories has garnered new interest in recent years, but implementing them is complicated by the need for SU(3) Clebsch-Gordon coefficients. In the loop-string-hadron (LSH) approach to lattice gauge theories, the elementary excitations are strictly gauge invariant, and constructing the basis requires no knowledge of Clebsch-Gordon coefficients. Originally developed for SU(2), the LSH formulation was recently generalized to SU(3), but limited to one spatial dimension. In this work, we generalize the LSH approach to constructing the basis of SU(3) gauge-invariant states at a trivalent vertex—the essential building block to multidimensional space. A direct generalization from the SU(2) vertex yields a legitimate basis; however, in certain sectors of the Hilbert space, the naive LSH basis vectors so defined suffer from being nonorthogonal. The issues with orthogonality are directly related to the “missing label” or “outer multiplicity” problem associated with SU(3) tensor products and may also be phrased in terms of Littlewood-Richardson coefficients or the need for a “seventh Casimir” operator. The states that are unaffected by the problem are orthonormalized in closed form. For the sectors that are afflicted, we discuss the nonorthogonal bases and their orthogonalization. A few candidates for seventh Casimir operators are readily constructed from the suite of LSH gauge-singlet operators. The diagonalization of a seventh Casimir represents one prescriptive solution toward obtaining a complete orthonormal basis, but a closed-form general solution remains to be found. Published by the American Physical Society 2025

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗