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.

SourceResearch connectionFrom yearThrough year

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

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

NASA NTRS1978

Characterization of oriented cracks with differential strain analysis

Linear strain of a rock sample as a function of hydrostatic pressure can be measured with a precision of 2 millionths. Such high-precision data for three orthogonal directions allow calculation of the distribution function for the porosity due to cracks' closing completely at a given pressure. Such data for at least six independent directions yield the zero-pressure strain tensor due to cracks' closing completely at a given pressure. The principal values and axes of this tensor distribution function provide information about the orientation of cracks as a function of closure pressure. The mathematical basis for the technique is developed, and the technique is illustrated with differential strain data for two samples, the Westerly (Rhode Island) granite and the Twin Sisters (Washington) dunite. Strain-tensor calculations reveal that each of these samples has a different type of anisotropic crack distribution.

Siegfried, R.↗

DOE OSTI2022

Automatic Calibration of a Geomechanical Model from Sparse Data for Estimating Stress in Deep Geological Formations

Summary In this study, we demonstrate geomechanical modeling with fully automatic parameter calibration to estimate the full geomechanical stress fields of a prospective US carbon dioxide (CO2) storage site, based on sparse measurement data. The goal is to compute full stress tensor field estimates (principal stresses and orientations) that are maximally compatible with observations within the constraints of the model assumptions, thereby extending pointwise, incomplete partial stress measurement to a simulated full formation stress field, as well as a rough assessment of the associated error. We use the Perch site, located in Otsego County, Michigan, USA, as our case study. The input data consist of partial stress tensor information inferred from in-situ borehole tests, geophysical well logs, and processing of seismic data. A static earth model (SEM) of the site was developed, and geomechanical simulation functionality of the open-source MATLAB Reservoir Simulation Toolbox (MRST) was used to model the stress field. Adjoint-based nonlinear optimization was used to adjust boundary conditions and material properties to calibrate simulated results of observations. Results were interpreted through a Bayesian framework. The focus of this paper is to demonstrate how the fully automatic calibration procedure works and discuss the results obtained; it does not attempt a detailed analysis of the stress field in the context of the proposed CO2 storage initiatives. Our work is part of a larger effort to noninvasively determine in-situ stresses in deep formations considered for CO2 storage. Guided by previously published research on geomechanical model calibration, our work presents a novel calibration approach supporting a potentially large number of linear or nonlinear calibration parameters to produce results optimally agreeing with available measurements and thus extend partial pointwise estimates to full tensor fields compatible with the physics of the site.

Engineering↗

DOE OSTI2024

Nucleon Isovector Axial Form Factors

We present results for the isovector axial vector form factors obtained using thirteen 2 + 1 + 1-flavor highly improved staggered quark (HISQ) ensembles generated by the MILC collaboration. The calculation of nucleon two- and three-point correlation functions has been done using Wilson-clover fermions. In the analysis of these data, we quantify the sensitivity of the results to strategies used for removing excited state contamination and invoke the partially conserved axial current relation between the form factors to choose between them. Our data driven analysis includes removing contributions from multihadron $Nπ$ states that make significant contributions. Our final results are $g_A$ = 1.292(53) stat (24) sys for the axial charge; $g_S$ = 1.085(50) stat (103) sys and $g_T$ = 0.991(21) stat (10) sys for the scalar and tensor charges; $\langle{r^2_A}\rangle$ = 0.439(56) sta t(34) sys fm 2 for the mean squared axial charge radius, $g^*_P$ = 9.03(47) stat (42) sys for the induced pseudoscalar charge; and $g_{πNN}$ = 14.14(81) stat (85) sys for the pion-nucleon coupling. We also provide a parametrization of the axial form factor $G_A(Q^2)$ over the range 0 ≤ $Q^2$ ≤ 1 GeV 2 for use in phenomenology and a comparison with other lattice determinations. We find that the various lattice data agree within 10% but are significantly different from the extraction of $G_A(Q^2)$ from the $ν$-deuterium scattering data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Frontal Slice Approaches for Tensor Linear Systems

Inspired by the row and column action methods for solving large-scale linear systems, in this work, we explore the use of frontal slices for solving tensor linear systems. In particular, this paper presents a novel approach for using frontal slices of a tensor $\mathcal{A}$ to solve tensor linear systems $\mathcal{A} ∗\mathcal{X} = \mathcal{B}$ where ∗ denotes the $t$-product. In addition, we consider variations of this method, including cyclic, block, and randomized approaches, each designed to optimize performance in different operational contexts. Our primary contribution lies in the development and convergence analysis of these methods. Experimental results on synthetically generated and real-world data, including applications such as image and video deblurring, demonstrate the efficacy of our proposed approaches and validate our theoretical findings.

Luo, Hengrui↗

NASA NTRS1989

A mechanics framework for a progressive failure methodology for laminated composites

A laminate strength and life prediction methodology has been postulated for laminated composites which accounts for the progressive development of microstructural damage to structural failure. A damage dependent constitutive model predicts the stress redistribution in an average sense that accompanies damage development in laminates. Each mode of microstructural damage is represented by a second-order tensor valued internal state variable which is a strain like quantity. The mechanics framework together with the global-local strategy for predicting laminate strength and life is presented in the paper. The kinematic effects of damage are represented by effective engineering moduli in the global analysis and the results of the global analysis provide the boundary conditions for the local ply level stress analysis. Damage evolution laws are based on experimental results.

Harris, Charles E.↗

DOE OSTI2023

Search for flavor-changing neutral-current couplings between the top quark and the Z boson with proton-proton collisions at s = 13 TeV with the ATLAS detector

A search for flavor-changing neutral-current couplings between a top quark, an up or charm quark, and a $Z$ boson is presented, using proton-proton collision data at $\sqrt{s}$ = 13 TeV collected by the ATLAS detector at the Large Hadron Collider. The analyzed data set corresponds to an integrated luminosity of 139 fb -1 . The search targets both single-top-quark events produced as $gq$ → $tZ$ (with $q$ = $u, c$) and top-quark-pair events, with one top quark decaying through the $t$ → $Zq$ channel. The analysis considers events with three leptons (electrons or muons), a $b$-tagged jet, possible additional jets, and missing transverse momentum. The data are found to be consistent with the background-only hypothesis and 95% confidence-level limits on the $t$ → $Zq$ branching ratios, assuming only tensor operators of the Standard Model effective field theory framework contribute to the $tZq$ vertices. These are 6.2 x 10 -5 (13 x 10 -5 ) for $t$ → $Zu$ ($t$ → $Zc$) for a left-handed $tZq$ coupling, and 6.6 x 10 -5 (12 x 10 -5 ) in the case of a right-handed coupling. These results are interpreted as 95% CL upper limits on the strength of the corresponding couplings, yielding limits for |${C}^{(13)*}_{uW}$| and |${C}^{(13)*}_{uB}$| (|${C}^{(31)}_{uW}$| and |${C}^{(31)}_{uB}$|) of 0.15 (0.16), and limits for |${C}^{(23)*}_{uW}$| and |${C}^{(23)*}_{uB}$| (|${C}^{(32)}_{uW}$| and |${C}^{(32)}_{uB}$|) of 0.22 (0.21), assuming a new-physics energy scale Λ NP of 1 TeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Nucleon isovector axial form factors

We present results for the isovector axial vector form factors obtained using thirteen 2 + 1 + 1 -flavor highly improved staggered quark (HISQ) ensembles generated by the MILC collaboration. The calculation of nucleon two- and three-point correlation functions has been done using Wilson-clover fermions. In the analysis of these data, we quantify the sensitivity of the results to strategies used for removing excited state contamination and invoke the partially conserved axial current relation between the form factors to choose between them. Our data driven analysis includes removing contributions from multihadron N π states that make significant contributions. Our final results are g A = 1.292 ( 53 ) stat ( 24 ) sys for the axial charge; g S = 1.085 ( 50 ) stat ( 103 ) sys and g T = 0.991 ( 21 ) stat ( 10 ) sys for the scalar and tensor charges; ⟨ r A 2 ⟩ = 0.439 ( 56 ) stat ( 34 ) sys fm 2 for the mean squared axial charge radius, g P * = 9.03 ( 47 ) stat ( 42 ) sys for the induced pseudoscalar charge; and g π N N = 14.14 ( 81 ) stat ( 85 ) sys for the pion-nucleon coupling. We also provide a parametrization of the axial form factor G A ( Q 2 ) over the range 0 ≤ Q 2 ≤ 1 GeV 2 for use in phenomenology and a comparison with other lattice determinations. We find that the various lattice data agree within 10% but are significantly different from the extraction of G A ( Q 2 ) from the ν -deuterium scattering data. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Chapter 4: Physically informed deep learning networks for simulating microstructure evolution of 3D polycrystals

As discussed in the previous chapter, high energy diffraction microscopy (HEDM) is used to study the micromechanical evolution of a material during in situ loading. HEDM experiments have been used to verify crystal plasticity (CP) simulations [119, 91, 90, 120], for experimental planning, material design, and to further analyze experimental results. However, Fast Fourier transform-based CP (CP-FFT) or finite element-based CP (CP-FE) methods are often too slow to be used in real-time during an experiment. CP-FFT is faster than CP-FE simulations due to the absence of meshing, but can still take hours to simulate the response of a single volume depending on the size and number of strain steps [127]. Reducing computation time would create a larger exploration space in planning and design, and enable faster analysis of experimental results and real-time feedback during an experiment. This research expands upon previous works to develop a workflow for predicting the full-field evolution of a 3D polycrystal. The workflow is simplified from previous works to predict only orientation and elastic strain tensors (from which stress tensors are calculated). The network is physically informed through loss functions and network architecture for a more robust model. The orientation predictions are informed about the cubic crystal symmetry of the material by incorporating disorientation and misorientation information into the network architecture and loss. The Von Mises stress is used to enforce the correct stress-strain trends in the strain tensor predictions. Additional total strain steps from the elastic and elastoplastic region are included to better capture the stress-strain evolution at smaller total strain steps. Material and hardening parameters are additional inputs into the networks to further inform the network and to study the network’s ability to predict different materials other than those used for training.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Effective Resistivity in Relativistic Collisionless Reconnection

Magnetic reconnection can power spectacular high-energy astrophysical phenomena by producing nonthermal energy distributions in highly magnetized regions around compact objects. By means of two-dimensional fully kinetic particle-in-cell (PIC) simulations, we investigate relativistic collisionless plasmoid-mediated reconnection in magnetically dominated pair plasmas with and without a guide field. In X-points, where diverging flows result in a nondiagonal thermal pressure tensor, a finite residence time for particles gives rise to a localized collisionless effective resistivity. Here, for the first time for relativistic reconnection in a fully developed plasmoid chain, we identify the mechanisms driving the nonideal electric field using a full Ohm law by means of a statistical analysis based on our PIC simulations. We show that the nonideal electric field is predominantly driven by gradients of nongyrotropic thermal pressures. We propose a kinetic physics motivated nonuniform effective resistivity model that is negligible on global scales and becomes significant only locally in X-points. It captures the properties of collisionless reconnection with the aim of mimicking its essentials in nonideal magnetohydrodynamic descriptions. This effective resistivity model provides a viable opportunity to design physically grounded global models for reconnection-powered high-energy emission.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1984

An analysis of the high-latitude thermospheric wind pattern calculated by a thermospheric general circulation model. I - Momentum forcing

A diagnostic processor (DP) was developed for analysis of hydrodynamic and thermodynamic processes predicted by the NCAR thermospheric general circulation model (TGCM). The TGCM contains a history file on the projected wind, temperature and composition fields at each grid point for each hour of universal time. The DP assimilates the history file plus ion drag tensors and drift velocities, specific heats, coefficients of viscosity, and thermal conductivity and calculates the individual forcing terms for the momentum and energy equations for a given altitude. Sample momentum forcings were calculated for high latitudes in the presence of forcing by solar radiation and magnetospheric convection with a 60 kV cross-tail potential, i.e., conditions on Oct. 21, 1981. It was found that ion drag and pressure forces balance out at F region heights where ion drift velocities are small. The magnetic polar cap/auroral zone boundary featured the largest residual force or net acceleration. Diurnal oscillations were detected in the thermospheric convection, and geostrophic balance was dominant in the E layer.

Killeen, T. L.↗

NASA NTRS1989

Predictions of Poisson's ratio in cross-ply laminates containing matrix cracks and delaminations

A damage-dependent constitutive model for laminated composites has been developed for the combined damage modes of matrix cracks and delaminations. The model is based on the concept of continuum damage mechanics and uses second-order tensor valued internal state variables to represent each mode of damage. The internal state variables are defined as the local volume average of the relative crack face displacements. Since the local volume for delaminations is specified at the laminate level, the constitutive model takes the form of laminate analysis equations modified by the internal state variables. Model implementation is demonstrated for the laminate engineering modulus E(x) and Poisson's ratio nu(xy) of quasi-isotropic and cross-ply laminates. The model predictions are in close agreement to experimental results obtained for graphite/epoxy laminates.

Harris, Charles E.↗

NASA NTRS2024

Analysis of faults and pit chains in Noctis Labyrinthus: Implications for early extension and possible magmatic plumbing

Noctis Labyrinthus has been a region of disputed origin due to its complexity and poor understanding of how various processes and mechanisms may have combined to form it. The surface is an integrated record of intensive tectonic activity expressed by a multiple extended sets of dip-slip faults oriented in different directions, and thought to have acted on this region over its history. These faults are always coalescent to pits and pit chains displaying a complicated geological history in the region. To understand this geological history, we mapped the surface features in Noctis Labyrinthus using the High-Resolution Stereo Camera (HRSC) onboard Mars Express ND2 nadir channel basemaps, and we adapted the Digital Terrain Map (DTM) from the Mission Experiment Gridded Data Record (MEGDR) of Mars Orbiter Laser Altimeter (MOLA) onboard Mars Global Surveyor (MGS) for the topography. We have investigated the spatial distribution and trend of fault systems, the pit chains' morphology, and the correlation between these two types of features. Our results show three fault systems: i) NS and NNE-SSW, ii) EW and ENE-WSW, and iii) NNW-SSE and NW. The analysis of the faults trending, cross-cutting correlation and the superimposition led to identify multiple intersections between these faults that have been alongside with the reactivations of some inherited faults. We interpreted the first system of fault to be related to coeval lateral extension, generated by regional stress tensor, which is probably related to the slight bending of Valles Marineris within two phases of bidirectional extension. The second system of faults has been generated by the radial oblate stress tensor related to the formation of the small shield volcanoes in Syria Planum. However, the third system is likely related to the external driving process, probably in the Tharsis province. We classified pits in four evolutionary stages based on their morphometric attributes. We believe that the formation of the pit chains in Noctis Labyrinthus is related to a surface collapse after a pressure drop related to the magma chamber deflation associated with Syria Planum volcanic province. We propose a deformational model based on early extension and magmatic plumbing as driving processes for the formation of Noctis Labyrinthus.

Mayssa El Yazidi↗

NASA NTRS2015

Analysis of Fiber Clustering in Composite Materials Using High-Fidelity Multiscale Micromechanics

A new multiscale micromechanical approach is developed for the prediction of the behavior of fiber reinforced composites in presence of fiber clustering. The developed method is based on a coupled two-scale implementation of the High-Fidelity Generalized Method of Cells theory, wherein both the local and global scales are represented using this micromechanical method. Concentration tensors and effective constitutive equations are established on both scales and linked to establish the required coupling, thus providing the local fields throughout the composite as well as the global properties and effective nonlinear response. Two nondimensional parameters, in conjunction with actual composite micrographs, are used to characterize the clustering of fibers in the composite. Based on the predicted local fields, initial yield and damage envelopes are generated for various clustering parameters for a polymer matrix composite with both carbon and glass fibers. Nonlinear epoxy matrix behavior is also considered, with results in the form of effective nonlinear response curves, with varying fiber clustering and for two sets of nonlinear matrix parameters.

composite materials↗

DOE OSTI2024

A reduced-order model for nonlinear radiative transfer problems based on moment equations and POD-Petrov-Galerkin projection of the normalized Boltzmann transport equation

A data-driven projection-based reduced-order model (ROM) for nonlinear thermal radiative transfer (TRT) problems is presented. The TRT ROM is formulated by (i) a hierarchy of low-order quasidiffusion (aka variable Eddington factor) equations for moments of the radiation intensity and (ii) the normalized Boltzmann transport equation (BTE). The multilevel system of moment equations is derived by projection of the BTE onto a sequence of subspaces which represent elements of the phase space of the problem. Exact closure for the moment equations is provided by the Eddington tensor. A Petrov-Galerkin (PG) projection of the normalized BTE is formulated using a proper orthogonal decomposition (POD) basis representing the normalized radiation intensity over the whole phase space and time. The Eddington tensor linearly depends on the solution of the normalized BTE. By linear superposition of the POD basis functions, a low-rank expansion of the Eddington tensor is constructed with coefficients defined by the PG projected normalized BTE. The material energy balance (MEB) equation is coupled with the effective gray low-order equations which exist on the same dimensional scale as the MEB equation. The resulting TRT ROM is structure and asymptotic preserving. A detailed analysis of the ROM is performed on the classical Fleck-Cummings (F-C) TRT multigroup test problem in 2D geometry. Numerical results are presented to demonstrate the ROM's effectiveness in the simulation of radiation wave phenomena. Importantly, the ROM is shown to produce solutions with sufficiently high accuracy while using low-rank approximation of the normalized BTE solution. Essential physical characteristics of supersonic radiation wave are preserved in the ROM solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Tensor factorization for elucidating mechanisms of piezoresponse relaxation via dynamic Piezoresponse Force Spectroscopy

Spatially resolved time and voltage-dependent polarization dynamics in PbTiO 3 thin films is explored using dynamic piezoresponse force microscopy (D-PFM) in conjunction with interferometric displacement sensing. This approach gives rise to 4D data sets containing information on bias-dependent relaxation dynamics at each spatial location without long-range electrostatic artifacts. To interpret these data sets in the absence of defined physical models, we employ a non-negative tensor factorization method which clearly presents the data as a product of simple behaviors allowing for direct physics interpretation. Correspondingly, we perform phase-field modeling finding the existence of ‘hard’ and ‘soft’ domain wall edges. This approach can be extended to other multidimensional spectroscopies for which even exploratory data analysis leads to unsatisfactory results due to many components in the decomposition.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Scanning X-Ray Diffraction Microscopy for Diamond Quantum Sensing

An understanding of nano- and microscale crystal strain in chemical-vapor-deposition diamond is crucial to the advancement of diamond quantum technologies. In particular, the presence of such strain and its characterization presents a challenge to diamond-based quantum sensing and information applications-as well as for future dark-matter detectors, where the directional information about incoming particles is encoded in crystal strain. Here, we exploit nanofocused scanning x-ray diffraction microscopy to quantitatively measure crystal deformation from defects in diamond with high spatial and strain resolution. The combination of information from multiple Bragg angles allows stereoscopic three-dimensional modeling of strain-feature geometry; the diffraction results are validated via comparison to optical measurements of the strain tensor based on spin-state-dependent spectroscopy of ensembles of nitrogen-vacancy centers in the diamond. Finally, our results demonstrate both strain and spatial resolution sufficient for directional detection of dark matter via x-ray measurement of crystal strain and provide a promising tool for diamond growth analysis and improvement of defect-based sensing.

42 ENGINEERING↗

DOE OSTI2024

The Simons Observatory: Combining cross-spectral foreground cleaning with multitracer B -mode delensing for improved constraints on inflation

The Simons Observatory (SO), due to start full science operations in early 2025, aims to set tight constraints on inflationary physics by inferring the tensor-to-scalar ratio r from measurements of cosmic microwave background (CMB) polarization B-modes. Its nominal design including three small-aperture telescopes (SATs) targets a precision σ⁡(r = 0) ≤ 0.003 without delensing. Achieving this goal and further reducing uncertainties requires a thorough understanding and mitigation of other large-scale B-mode sources such as Galactic foregrounds and weak gravitational lensing. We present an analysis pipeline aiming to estimate r by including delensing within a cross-spectral likelihood, and demonstrate it for the first time on SO-like simulations accounting for various levels of foreground complexity, inhomogeneous noise and partial sky coverage. As introduced in an earlier SO delensing paper, lensing B-modes are synthesized using internal CMB lensing reconstructions as well as Planck-like cosmic infrared background maps and LSST-like galaxy density maps. We then extend SO’s power-spectrum-based foreground-cleaning algorithm to include all auto- and cross-spectra between the lensing template and the SAT B-modes in the likelihood function. This allows us to constrain r and the parameters of our foreground model simultaneously. Within this framework, we demonstrate the equivalence of map-based and cross-spectral delensing and use it to motivate an optimized pixel-weighting scheme for power spectrum estimation. We start by validating our pipeline in the simplistic case of uniform foreground spectral energy distributions. In the absence of primordial B-modes, we find that the 1⁢σ statistical uncertainty on r, σ⁡(r), decreases by 37% as a result of delensing. Tensor modes at the level of r = 0.01 are successfully detected by our pipeline. Even when using more realistic foreground models including spatial variations in the dust and synchrotron spectral properties, we obtain unbiased estimates of r both with and without delensing by employing the moment-expansion method. In this case, uncertainties are increased due to the higher number of model parameters, and delensing-related improvements range between 27% and 31%. These results constitute the first realistic assessment of the delensing performance at SO’s nominal sensitivity level.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

A Visual Understanding of Circular Dichroism Spectroscopy

Mapping chemical and structural properties to electronic and magnetic responses is critical to many applications such as quantum information science, where the precise storage and transmission of unique information is paramount. Specifically, constructing molecules and materials that provide strong polarized responses at tunable frequencies and with large anisotropies is key to optical processing of quantum information. Chiral molecules provide chiroptical response to circularly polarized light, making them attractive for quantum information science and other applications related to sensing, polarized photodetectors, and spintronics. Predicting a molecular design, a priori, with large anisotropies to circularly polarized light is challenging due to the complex interplay between electric and magnetic components of the optical response. In this work, we explore a visual representation of the electronic chiroptical response by decomposing the rotary strength into its constituent components. Here, we make use of the intuitive electronic oscillator framework to develop classical intuition regarding the rotary strength and its constituents. We explore three model chemical systems that exhibit local and global chirality. Our analysis reveals that local chirality necessarily exhibits competition between the local chiral center and chirality induced in other fragments of the molecule, resulting in both unexpected nonmonotonic trends and sign flips in chemically adjacent geometries. Furthermore, we can visually distinguish between local and global chirality via examination of the transition chiral tensor. Interestingly, we make strong connections to ferromagnetic and antiferromagnetic spin systems in that chiroptically inactive transitions exhibit antiferromagnetic-like alternating orbital patterns while active transitions show domain formation in an ferromagnetic-like alignment that produces a net chiroptical response.

36 MATERIALS SCIENCE↗