Search NASA⌕ Search

SEARCH · Search NASA

Results for “TENSOR ANALYSIS”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Broadscale Postseismic Gravity Change Following the 2011 Tohoku-Oki Earthquake and Implication for Deformation by Viscoelastic Relaxation and Afterslip

The analysis of GRACE gravity data revealed post-seismic gravity increase by 6 micro-Gal over a 500 km scale within a couple of years after the 2011 Tohoku-Oki earthquake, which is nearly 40-50% of the co-seismic gravity change. It originates mostly from changes in the isotropic component corresponding to the M(sub rr) moment tensor element. The exponential decay with rapid change in a year and gradual change afterward is a characteristic temporal pattern. Both viscoelastic relaxation and afterslip models produce reasonable agreement with the GRACE free-air gravity observation, while their Bouguer gravity patterns and seafloor vertical deformations are distinctly different. The post-seismic gravity variation is best modeled by the bi-viscous relaxation with a transient and steady state viscosity of 10(exp 18) and 10(exp 19) Pa s, respectively, for the asthenosphere. Our calculated higher-resolution viscoelastic relaxation model, underlying the partially ruptured elastic lithosphere, yields the localized post-seismic subsidence above the hypocenter reported from the GPS-acoustic seafloor surveying.

GRACE↗
DOE OSTI2024

Random insights into the complexity of two-dimensional tensor network calculations

Projected entangled pair states (PEPS) offer memory-efficient representations of some quantum many-body states that obey an entanglement area law and are the basis for classical simulations of ground states in two-dimensional (2d) condensed matter systems. However, rigorous results show that exactly computing observables from a 2d PEPS state is generically a computationally hard problem. Yet approximation schemes for computing properties of 2d PEPS are regularly used, and empirically seen to succeed, for a large subclass of (“not too entangled”) condensed matter ground states. Adopting the philosophy of random matrix theory, in this work, we analyze the complexity of approximately contracting a 2d random PEPS by exploiting an analytic mapping to an effective replicated statistical mechanics model that permits a controlled analysis at a large bond dimension. Through this statistical-mechanics lens, we argue that (i) although approximately sampling wave-function amplitudes of random PEPS faces a computational-complexity phase transition above a critical bond dimension, and (ii) one can generically efficiently estimate the norm and correlation functions for any finite bond dimension. Furthermore, these results are supported numerically for various bond-dimension regimes. It is an important open question whether the above results for random PEPS apply more generally also to PEPS representing physically relevant ground states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

The use of digital thread for reconstruction of local fiber orientation in a compression molded pin bracket via deep learning

A deep convolutional neural network (DCNN) was used for microstructure reconstruction using artificial intelligence (MR-AI) by predicting local average fiber orientation distributions (FOD) in a 3D prepreg platelet molded composite (PPMC) pin bracket. To train the MR-AI model, surface strain fields from residual stresses simulated in PPMC plates were used as the input to the DCNN. A training dataset included PPMC plates with various degrees of global fiber alignment, based on the information obtained from high-fidelity flow simulation of a pin bracket. Further, the MR-AI model was then deployed to analyze FOD in the 3D pin bracket by conducting thermo-elastic residual stress analysis. Initially, the MR-AI model was established entirely on the synthetic simulation data. Then, a μCT scan of a physically molded pin bracket was used to create a finite element model that provided data for additional validation of the DCNN model. For the μCT scan finite element pin bracket the MR-AI model predicted the distribution of fiber orientation tensor components with MAE of 0.10 indicating a global prediction error of 10%. For the flow simulated pin bracket, the MR-AI model predicted the distribution of fiber orientation tensor components with a global prediction error of 11%. The MR-AI model showed the ability to predict regions of varying alignment in the base and flange of the pin bracket. The proposed MR-AI methodology allows for rapid prediction of FOD in geometrically complex parts and offers a promising path to detecting unique fiber orientation states in molded components.

42 ENGINEERING↗
DOE OSTI2023

Third-order Hall effect in the surface states of a topological insulator

Time reversal and inversion symmetric materials fail to yield linear and nonlinear responses since they possess net zero Berry curvature. However, higher-order Hall response can be generated in these systems upon constraining the crystalline symmetries. Motivated by the recently discovered third-order Hall (TOH) response mediated by Berry connection polarizability, namely, the variation of the Berry connection with respect to an applied electric field, here, we investigate the existence of such a Hall effect in the surface states of hexagonal warped topological insulators (e.g., Bi 2 Te 3 ) under the application of only electric field. Using the semiclassical Boltzmann formalism, we investigate the effect of tilt and hexagonal warping on the Berry connection polarizability tensor and, consequently, the TOH effect, provided the Dirac cone remains gapless. We find that the magnitude of the response increases significantly with increasing tilt strength and warping, and therefore, they can provide the tunability of this effect. In addition, we also explore the effect of chemical doping on the TOH response in this system. Interestingly, we show, based on the symmetry analysis, that the TOH can be the leading-order response in this system, which can directly be verified in experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

TuckerMPI: A Parallel C++/MPI Software Package for Large-scale Data Compression via the Tucker Tensor Decomposition

With this study, our goal is compression of massive-scale grid-structured data, such as the multi-terabyte output of a high-fidelity computational simulation. For such data sets, we have developed a new software package called TuckerMPI, a parallel C++/MPI software package for compressing distributed data. The approach is based on treating the data as a tensor, i.e., a multidimensional array, and computing its truncated Tucker decomposition, a higher-order analogue to the truncated singular value decomposition of a matrix. The result is a low-rank approximation of the original tensor-structured data. Compression efficiency is achieved by detecting latent global structure within the data, which we contrast to most compression methods that are focused on local structure. In this work, we describe TuckerMPI, our implementation of the truncated Tucker decomposition, including details of the data distribution and in-memory layouts, the parallel and serial implementations of the key kernels, and analysis of the storage, communication, and computational costs. We test the software on 4.5 and 6.7 terabyte data sets distributed across 100 s of nodes (1,000 s of MPI processes), achieving compression ratios between 100 and 200,000×, which equates to 99--99.999% compression (depending on the desired accuracy) in substantially less time than it would take to even read the same dataset from a parallel file system. Moreover, we show that our method also allows for reconstruction of partial or down-sampled data on a single node, without a parallel computer so long as the reconstructed portion is small enough to fit on a single machine, e.g., in the instance of reconstructing/visualizing a single down-sampled time step or computing summary statistics. The code is available at https://gitlab.com/tensors/TuckerMPI.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1971

On the dynamics of a human body model.

Equations of motion for a model of the human body are developed. Basically, the model consists of an elliptical cylinder representing the torso, together with a system of frustrums of elliptical cones representing the limbs. They are connected to the main body and each other by hinges and ball and socket joints. Vector, tensor, and matrix methods provide a systematic organization of the geometry. The equations of motion are developed from the principles of classical mechanics. The solution of these equations then provide the displacement and rotation of the main body when the external forces and relative limb motions are specified. Three simple example motions are studied to illustrate the method. The first is an analysis and comparison of simple lifting on the earth and the moon. The second is an elementary approach to underwater swimming, including both viscous and inertia effects. The third is an analysis of kicking motion and its effect upon a vertically suspended man such as a parachutist.

Huston, R. L.↗
NASA NTRS1982

Prediction of molecular infrared band and line intensities

The general problem of predicting infrared intensities is examined. The problem of predicting line intensities is discussed with attention focussed specifically on the problem of predicting vibrational strengths (or total band intensities). Two methods are considered: predictions based on ab initio quantum mechanical calculations, and predictions based on the transfer of atomic polar tensor parameters from small model molecules to the molecule of interest. Comparisons of predicted spectra with experimental spectra are shown for CH3OH, CH3OCH3, C2H5OH, and HOC2. It is suggested that the most important application of this work may be to the prediction of intensities for free radicals, ions, and reactive intermediates so that the infrared spectrum can be used in the quantitative analysis for these interesting species.

Person, W. B.↗
DOE OSTI2025

The structure of quantum singularities on a Cauchy horizon

Spacetime singularities pose a long-standing puzzle in quantum gravity. Unlike Schwarzschild, a generic family of black holes gives rise to a Cauchy horizon on which, even in the Hartle-Hawking state, quantum observables such as ⟨T μν ⟩ — the expectation value of the stress-energy tensor — can diverge, causing a breakdown of semiclassical gravity. Because they are diagnosed within quantum field theory (QFT) on a smooth background, these singularities may provide a better-controlled version of the spacetime singularity problem, and merit further study. Here, I highlight a mildness puzzle of Cauchy horizon singularities: the ⟨T μν ⟩ singularity is significantly milder than expected from symmetry and dimensional analysis. I address the puzzle in a simple spacetime ${\mathcal{W}}_P$, which arises universally near all black hole Cauchy horizons: the past of a codimension-two spacelike plane in flat spacetime. Specifically, I propose an extremely broad QFT construction in which, roughly speaking, Cauchy horizon singularities originate from operator insertions in the causal complement of the spacetime. The construction reproduces well-known outer horizon singularities (e.g., in the Boulware state), and remarkably, when applied to ${\mathcal{W}}_P$, gives rise to a universal mild singularity structure for robust singularities, ones whose leading singular behavior is state-independent. I make non-trivial predictions for all black hole Cauchy horizon singularities using this, and discuss extending the results beyond robust singularities and the strict near Cauchy horizon limit.

black holes↗
NASA NTRS1973

Analysis of multiple pulse NMR in solids

The general problems associated with the removal of the effects of dipolar broadening from solid-state NMR spectra are analyzed. The effects of finite pulse width and H sub 1 inhomogeneity are shown to have limited the resolution of previous pulse cycles, and a new eight-pulse cycle designed to minimize these problems is discussed. Spectra for F-19 in CaF2 taken with this cycle are presented which show residual linewidth near 10 Hz. The feasibility of measuring proton chemical shift tensors is discussed.

Rhim, W.-K.↗
DOE OSTI2021

Strongly Anisotropic Thermomechanical Response to Shock Wave Loading in Oriented Samples of the Triclinic Molecular Crystal 1,3,5-Triamino-2,4,6-trinitrobenzene

All-atom molecular dynamics (MD) simulations were used to study shock wave loading in oriented single crystals of the highly anisotropic triclinic molecular crystal 1,3,5-triamino-2,4,6-trinitrobenzene (TATB). The crystal structure consists of planar hydrogen-bonded sheets of individually planar TATB molecules that stack into graphitic-like layers. Shocks were studied for seven systematically prepared crystal orientations with limiting cases that correspond to shock propagation exactly perpendicular and exactly parallel to the graphitic-like layers. The simulations were performed for initially defect-free crystals using a reverse-ballistic configuration that generates explicit, supported shocks. Final longitudinal stress components are between ≈8.5 and ≈10.5 GPa for the 1.0 km s –1 impact speed studied. Orientation-dependent properties are reported including shock speeds, stresses, temperatures, compression ratios, and local material strain rates. Spatiotemporal maps of the temperature, stress tensor, material flow, and molecular orientations reveal complicated processes that arise for specific shock directions. Furthermore, the results indicate that TATB shock response is highly sensitive to crystal orientation, with significant qualitative differences for the time evolution of the stress tensor and temperature, elastic/inelastic compression response, defect formation and growth, critical von Mises stress, and strain rates during shock rise that span nearly an order of magnitude. A variety of inelastic deformation mechanisms are identified, ranging from crumpling of graphitic-like layers to dislocation-mediated plasticity to intense shear strain localization. To our knowledge, these are the first systematic MD simulations and analysis of explicit shock wave propagation along nontrivial crystal directions in a triclinic molecular crystal.

36 MATERIALS SCIENCE↗
NASA NTRS1982

Reliability of the finite element method for calculating free edge stresses in composite laminates

The interlaminar normal stress distributions along the interface between the +45 deg and -45 deg plies of a graphite/epoxy laminate, obtained by various investigators, were found to disagree in both magnitude and sign. The reliability of the displacement-formulated finite element method in analyzing the edge-stress problem of a composite laminate is investigated. The history of the edge-stress problem is reviewed, and two well-known elasticity problems, one involving a stress discontinuity and one a singularity, are analyzed. The finite element analysis in these problems yields accurate stress distributions everywhere except in two elements closest to the stress discontinuity or singularity. Stress distributions for a + or -45 deg ply laminate near the singularity were similar to those of the two elasticity problems, demonstrating the methods, accuracy for calculating interlaminar stresses in composite laminates. The disagreement between the numerical methods was attributed to the unsymmetric stress tensor at singularity.

Whitcomb, J. D.↗
DOE OSTI2022

G2Aero (Separable shape tensors for aerodynamic design) [SWR-22-44]

G2Aero is a flexible and practical tool for design and deformation of 2D airfoils and 3D blades using data-driven approaches. G2Aero utilizes the geometry of matrix manifolds—specifically the Grassmannian—to build a novel framework for representing physics-based separable deformations of shapes. G2Aero offers the flexibility to generate perturbations in a customizable way over any portion of the blade. The G2Aero framework utilizes data-driven methods based on a curated database of physically relevant airfoils. Specific tools include: - principal geodesic analysis over normal coordinate neighborhoods of matrix manifolds; - a variety of data-regularized deformations to nominal 2D airfoil shapes; - Riemannian interpolation connecting a sequence of airfoil cross-sections to build 3D blades from 2D data; - consistent perturbations over the span of interpolated 3D blades based on dominant modes from the data-driven analysis.

Doronina, Olga↗
DOE OSTI2020

Normal modes in thermal AdS via the Selberg zeta function

The heat kernel and quasinormal mode methods of computing 1-loop partition functions of spin s s fields on hyperbolic quotient spacetimes \mathbb{H}^{3}/\mathbb{Z} ℍ 3 / ℤ are related via the Selberg zeta function. We extend that analysis to thermal \text{AdS}_{2n+1} AdS 2 n + 1 backgrounds, with quotient structure \mathbb{H}^{2n+1}/\mathbb{Z} ℍ 2 n + 1 / ℤ . Specifically, we demonstrate the zeros of the Selberg function encode the normal mode frequencies of spin fields upon removal of non-square-integrable modes. With this information we construct the 1-loop partition functions for symmetric transverse traceless tensors in terms of the Selberg zeta function and find exact agreement with the heat kernel method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Physics Division Strategic Plan Fiscal Years 2020-2024

The vision of the Physics Division (PHY) at Argonne National Laboratory is to continue enhancing its role as a world-leading institution in basic nuclear physics research and its applications. Key to this vision is for PHY to continue to safely and effectively operate and evolve the capabilities of the Argonne Tandem Linac Accelerator System (ATLAS) facility to best serve its users. ATLAS is the Department of Energy (DOE) accelerator facility for low-energy nuclear physics research. The research carried out by PHY covers many themes in contemporary science but can be distilled into five areas of focus. These themes are interconnected, with continuously evolving synergies between the various groups and facilities in PHY and the broader Laboratory. This five-year strategic plan serves to illustrate our current capabilities and new directions related to these five areas of focus (ATLAS also develops an independent strategic plan): Accelerator research and design. The goal of this theme is to design, fabricate, test, and implement novel accelerator systems, with a focus on high-intensity ion and electron systems. These R&D activities have led to enhancements in the ATLAS accelerator system. Among others, the division’s accelerator systems are in use or planned for use at Fermilab, the Advanced Photon Source, and in a broad variety of applications at the Facility for Rare Isotope Beams (FRIB). Atom trapping and fundamental symmetries. The goal of this theme is to explore and exploit the uses of advanced laser cooling and trapping techniques to manipulate atoms. There are three main areas of focus. First is in the application of the atom trap trace analysis (ATTA) technique for age determination of groundwater and ice by using radio-krypton dating. The next two trapping-based programs involve tests of fundamental symmetries in nature: the cooling and trapping of radium-225 with the aim of determining limits on an observation of its electric dipole moment and precision measurements of the beta decay properties of helium-6 to set limits on the tensor coupling constant. This is complemented by measurements of similar properties in lithium-8 and boron-8. Nuclear astrophysics. The goal of this theme is to enhance our understanding of how elements are created in the universe via explosive nucleosynthesis and how stars evolve. To meet the challenges of this theme, many of the capabilities of ATLAS are being enhanced, including the development of new beams through a new in-flight separator (RAISOR) and the anticipated neutron-generator upgrade of the Californium Rare Isotope Breeder Upgrade facility (nuCARIBU). These are coupled to state-of-the-art instruments such as the Canadian Penning Trap, Helical Orbit Spectrometer, Gammasphere, GRETINA, Multi-Sampling Ionization Chamber, the Fragment Mass Analyzer, and Argonne Gas-filled Fragment Analyzer, and a new low-background experimental area for decay studies. Nuclear structure. The goal of this theme is to understand the structure of nuclei, both stable and radioactive, in terms of single-particle properties, their shapes, and the dynamics governing reactions between them. These include questions such as what are the limits of nuclear stability and what are the properties of super heavy nuclei. As with the nuclear astrophysics theme, the capabilities of ATLAS, guided by the ATLAS user community, are continuously being enhanced to this end. The Division has strategic initiatives to play a leading role in the development of instrumentation and research programs at the FRIB. Quantum chromodynamics (QCD) and hadron physics. The goal of this theme is to lead major research programs focused on revealing the quark and gluon structure of protons, neutrons, nuclei, and short-lived mesons and baryons. These involve major programs at Jefferson Lab, Fermilab, and smaller facilities. They are complemented by theoretical endeavors centered on the question of how hadrons and their properties emerge from QCD. Out of these activities arise strategic initiatives to play a major role in the forthcoming Electron-Ion Collider (EIC).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Demonstration of CFD to support the coupled analysis of a reactor pressure vessel subjected to pressurized thermal shock

The structural components that comprise nuclear reactors and their supporting structures are subjected to harsh operating environments that can challenge their integrity, especially after exposure for extended duration or under accident condition. As one of the most significant components of a reactor, the Reactor Pressure Vessel (RPV) is exposed to an aggressive environment during the operation time (e.g. more than 40 years). Aging degradation mechanisms (e.g. thermo-fatigue) could grow initial defects up to a critical size, increasing the susceptibility to failure in the RPV. The conventional methods are mostly based on simple crack and structure geometries. Very limited studies consider the real conditions of the RPV subjected to a thermal shock due to a Loss of Coolant Accident (LOCA). During a LOCA event, the most severe conditions take place when the emergency core cooling (ECC) water is injected inside the cold legs filled initially with hotter water and/or steam. The rapid cooling of the down-comer and the internal RPV surface followed probably by re-pressurization of the RPV causes large temperature gradients and variation of pressure which induces thermal-mechanical stresses. In order to develop the model for integrity assessment of a reactor pressure vessel (RPV) subjected to pressurized thermal shock (PTS), a multi-physics simulation, which includes the thermo-hydraulic, thermo-mechanical and fracture mechanics analyses is necessary. The prediction of the temperature field is achieved by using computational fluid dynamics (CFD) simulation. In this report, a demonstration CFD standalone simulation is performed to support coupled analysis for Reactor Pressure Vessel (RPV) subjected to Pressurized Thermal Shock (PTS). The study use a simplified computational domain to represents a real RPV. The purpose of the study is to demonstrate the transient temperature response of RPV to ECC injection. The CFD model is built in a robust and efficient way for further coupled calculation. The next steps of this work, including the coupled thermal and tensor mechanics capabilities using Cardinal are expected to be complete by the end of FY21 for the demo problem. After this, into FY22, the capability will be demonstrated for a realistic RPV.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1991

Seals Research at Texas A/M University

The Turbomachinery Laboratory at Texas A&M has been providing experimental data and computational codes for the design seals for many years. The program began with the development of a Halon based seal test rig. This facility provided information about the effective stiffness and damping in whirling seals. The Halon effectively simulated cryogenic fluids. Another test facility was developed (using air as the working fluid) where the stiffness and damping matrices can be determined. This data was used to develop bulk flow models of the seal's effect upon rotating machinery; in conjunction with this research, a bulk flow model for calculation of performance and rotordynamic coefficients of annular pressure seals of arbitrary non-uniform clearance for barotropic fluids such as LH2, LOX, LN2, and CH4 was developed. This program is very efficient (fast) and converges for very large eccentricities. Currently, work is being performed on a bulk flow analysis of the effects of the impeller-shroud interaction upon the stability of pumps. The data was used along with data from other researchers to develop an empirical leakage prediction code for MSFC. Presently, the flow field inside labyrinth and annular seals are being studied in detail. An advanced 3-D Doppler anemometer system is being used to measure the mean velocity and entire Reynolds stress tensor distribution throughout the seals. Concentric and statically eccentric seals were studied; presently, whirling seals are being studied. The data obtained are providing valuable information about the flow phenomena occurring inside the seals, as well as a data base for comparison with numerical predictions and for turbulence model development. A finite difference computer code was developed for solving the Reynolds averaged Navier Stokes equation inside labyrinth seals. A multi-scale k-epsilon turbulence model is currently being evaluated. A new seal geometry was designed and patented using a computer code. A large scale, 2-D seal flow visualization facility is also being developed.

Morrison, Gerald L.↗
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.↗
NASA NTRS1998

Quantity, Revisited: An Object-Oriented Reusable Class

"Quantity", a prototype implementation of an object-oriented class, was developed for two reasons: to help engineers and scientists manipulate the many types of quantities encountered during routine analysis, and to create a reusable software component to for large domain-specific applications. From being used as a stand-alone application to being incorporated into an existing computational mechanics toolkit, "Quantity" appears to be a useful and powerful object. "Quantity" has been designed to maintain the full engineering meaning of values with respect to units and coordinate systems. A value is a scalar, vector, tensor, or matrix, each of which is composed of Value Components, each of which may be an integer, floating point number, fuzzy number, etc., and its associated physical unit. Operations such as coordinate transformation and arithmetic operations are handled by member functions of "Quantity". The prototype has successfully tested such characteristics as maintaining a numeric value, an associated unit, and an annotation. In this paper we further explore the design of "Quantity", with particular attention to coordinate systems.

Funston, Monica Gayle↗