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 163 records · Page 9

Short-range correlation physics at low renormalization group resolution

Recent experiments have succeeded in isolating processes for which short-range correlation (SRC) physics is dominant and well accounted for by SRC phenomenology. But an alternative and compelling picture emerges from renormalization group (RG) evolution to low RG resolution. At high RG resolution, SRCs are identified as components in the nuclear wave function with relative pair momenta greater than the Fermi momentum. Scale separation results in wave-function factorization that can be exploited with phenomenologies such as the generalized contact formalism or the low-order correlation operator approximation. Evolution to lower resolution shifts SRC physics from nuclear structure to the reaction operators without changing the measured observables. We show how the features of SRC phenomenology manifested at high RG resolution are cleanly identified at low RG resolution using simple two-body operators and local-density approximations with uncorrelated wave functions, all of which can be systematically generalized. We verify that the experimental consequences to date follow directly at low resolution from well-established properties of nucleon-nucleon interactions such as the tensor force. Thus the RG reconciles the contrasting pictures of the same experiment and shows how to get correct results using wave functions without SRC components. Finally, our demonstration has implications for the analysis of knockout reactions for which SRC physics is not cleanly isolated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Enhanced Microseismicity During Production Pumping Cessation at the San Emidio Geothermal Field (Nevada, USA) in December 2016

Abstract Tectonic activity, geothermal fluids, and microseismic events (MSEs) tend to occur in similar locations as a result of spatiotemporal changes in the subsurface stress state. To quantify this association, we analyze data from a dense seismic array deployed at the San Emidio geothermal field, Nevada for 1 week in December 2016 to coincide with a 19.45‐hr shutdown of all injection and production pumping operations. 123 MSEs were detected, of which 101 occurred during the shutdown. The spatial association of the MSEs with the production wells suggests a causal relationship between the production cessation and the MSEs. Here we performed a detailed analysis to investigate reservoir material properties, distribution of seismically activated faults, and local stress state. We determined the hypocenters, magnitudes, and focal mechanisms for the MSEs, P‐wave tomographic velocity model, and local stress tensor. The results show that most MSEs occurred near the production wells. Magnitudes fall between −2.2 and 0.0 with larger events located closer to the production wells. Most MSEs occurred within a westward‐dipping normal fault zone in the reservoir associated with anomalously low P‐wave velocity values. The focal mechanism and stress inversion results show predominantly normal faulting with the maximum horizontal stress oriented north‐south. We suggest that the MSEs during shutdown were triggered on pre‐existing, small‐scale, critically stressed fault patches in the reservoir as the pore pressure increased around the production wells when the production pumping ceased. We interpret the larger MSE magnitudes closer to the production wells as a result of higher pore pressure increase.

15 GEOTHERMAL ENERGY↗

Colossal transverse magnetoresistance due to nematic superconducting phase fluctuations in a copper oxide

Electronic anisotropy (“nematicity”) has been detected in cuprate superconductors by various experimental techniques. Using angle-resolved transverse resistance (ARTR) measurements, a very sensitive and background-free technique that can detect 0.5% anisotropy in transport, we have observed it also in La 2-x Sr x CuO 4 (LSCO) for 0.02 ≤ x ≤ 0.25. A central enigma in LSCO is the rotation of the nematic director (orientation of the largest longitudinal resistance) with temperature; this has not been seen before in any material. Here, we address this puzzle by measuring the angle-resolved transverse magnetoresistance (ARTMR) in LSCO. We report the discovery of colossal transverse magnetoresistance (CTMR)—an order-of-magnitude drop in the transverse resistivity in the magnetic field of 6 T. We show that the apparent rotation of the nematic director is caused by anisotropic superconducting fluctuations, which are not aligned with the normal electron fluid, consistent with coexisting bond-aligned and diagonal nematic orders. We quantify this by modeling the (magneto-)conductivity as a sum of normal (Drude) and paraconducting (Aslamazov–Larkin) channels but extended to contain anisotropic Drude and Cooper-pair effective mass tensors. Strikingly, the anisotropy of Cooper-pair stiffness is much larger than that of the normal electrons. It grows dramatically on the underdoped side, where the fluctuations become quasi-one-dimensional. Our analysis is general rather than model dependent. Still, we discuss some candidate microscopic models, including coupled strongly-correlated ladders where the transverse (interladder) phase stiffness is low compared with the longitudinal intraladder stiffness, as well as the anisotropic superconducting fluctuations expected close to the transition to a pair-density wave state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Geometry-aware training of factorized layers in tensor Tucker format

Reducing parameter redundancies in neural network architectures is crucial for achieving feasible computational and memory requirements during train and inference of large networks. Given its easy implementation and flexibility, one promising approach is layer factorization, which reshapes weight tensors into a matrix format and parameterizes it as the product of two rank-r matrices. However, this family of approaches often requires an initial full-model warm-up phase, prior knowledge of a feasible rank, and it is sensitive to parameter initialization.In this work, we introduce a novel approach to train the factors of a Tucker decomposition of the weight tensors. Our training proposal proves to be optimal in locally approximating the original unfactorized dynamics and stable for the initialization. Furthermore, the rank of each mode is dynamically updated during training.We provide a theoretical analysis of the algorithm, showing convergence, approximation and local descent guarantees. The method's performance is further illustrated through a variety of experiments, showing remarkable training compression rates and comparable or even better performance than the full baseline and alternative layer factorization strategies.

Zangrando, Emanuele [Gran Sasso Science Institute ↗

Moment expansion of polarized dust SED: A new path towards capturing the CMB B -modes with LiteBIRD

Accurate characterization of the polarized dust emission from our Galaxy will be decisive in the quest for the cosmic microwave background (CMB) primordial B-modes. An incomplete modeling of its potentially complex spectral properties could lead to biases in the CMB polarization analyses and to a spurious measurement of the tensor-to-scalar ratio r. It is particularly crucial for future surveys like the LiteBIRD satellite, the goal of which is to constrain the faint primordial signal leftover by inflation with an accuracy on the tensor-to-scalar ratio r of the order of 10 –3 . Variations of the dust properties along and between lines of sight lead to unavoidable distortions of the spectral energy distribution (SED) that cannot be easily anticipated by standard component-separation methods. This issue can be tackled using a moment expansion of the dust SED, an innovative parametrization method imposing minimal assumptions on the sky complexity. In the present paper, we apply this formalism to the B-mode cross-angular power spectra computed from simulated LiteBIRD polarization data at frequencies between 100 and 402 GHz that contain CMB, dust, and instrumental noise. The spatial variation of the dust spectral parameters (spectral index β and temperature T) in our simulations lead to significant biases on r (~21 σ r ) if not properly taken into account. Performing the moment expansion in β, as in previous studies, reduces the bias but does not lead to sufficiently reliable estimates of r. We introduce, for the first time, the expansion of the cross-angular power spectra SED in both β and T, showing that, at the sensitivity of LiteBIRD, the SED complexity due to temperature variations needs to be taken into account in order to prevent analysis biases on r. Thanks to this expansion, and despite the existing correlations between some of the dust moments and the CMB signal responsible for a rise in the error on r, we can measure an unbiased value of the tensor-to-scalar ratio with a dispersion as low as σ r = 8.8 × 10 –4 .

79 ASTRONOMY AND ASTROPHYSICS↗

Resummation of QED radiative corrections in a strong constant crossed field

By considering radiative corrections of up to 3rd-loop order, Ritus and Narozhny conjectured that the proper expansion parameter for QED in a strong constant crossed field is $g=\alpha\chi^{2/3}$, where the dynamical quantum parameter $\chi=e\sqrt{-(Fp)^2}/m^3$ combines the particle momentum $p$ with the external field strength tensor $F$. Here we present and discuss the first non-perturbative result in this context, the resummed bubble-type polarization corrections to the electron self-energy in a constant crossed field. Our analysis confirms the relevance of the scaling parameter $g$ to the enhancement of bubble-type radiative corrections. This parameter actually represents the characteristic value of the ratio of the 1-loop polarization bubble to the photon virtuality. After an all-order resummation we identify and discuss two contributions to the self-energy with different formation regions and asymptotic behavior for $g\gg1$. Whereas the breakdown of perturbation theory occurs already for $g\gtrsim1$, the leading-order result remains dominant until the asymptotic regime $g\gg 1$ is reached. However, the latter is specific to processes like elastic scattering or photon emission and does not have to remain true for general higher-order QED processes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Reduced order models for thermal radiative transfer problems based on moment equations and data-driven approximations of the Eddington tensor

Here a new group of structure and asymptotic preserving reduced-order models (ROMs) for multidimensional nonlinear thermal radiative transfer (TRT) problems is presented. They are formulated by means of the nonlinear projective approach and data compression techniques. The nonlinear projection is applied to the Boltzmann transport equation (BTE) to derive a hierarchy of low-order moment equations. Approximation of the Eddington tensor that provides exact closure for the system of moment equations is found with projection-based data-driven methodologies. These include the (i) proper orthogonal decomposition (POD), (ii) dynamic mode decomposition (DMD) and (iii) a variant of the DMD. A parameterization is derived for this ROM for the temperature of radiation incoming to the problem domain (the radiation drive temperature). This parameterization is informed from results of a dimensionless study of the TRT problem. Analysis of the ROMs is performed on the classical Fleck-Cummings TRT multigroup test problem in 2D geometry with a radiation-driven Marshak wave. Numerical results are presented to demonstrate the performance of these ROMs for the simulation of evolving radiation and heat waves. Results show these models to be sufficiently accurate for practical computations with rather low-rank representations of the Eddington tensor. As the rank of the approximation is increased, the errors of solutions generated by the ROMs gradually decreases.

42 ENGINEERING↗

Effective point-charge analysis of crystal fields: Application to rare-earth pyrochlores and tripod kagome magnets 𝑅 ⁢3⁢ Mg⁢ 2 ⁢Sb⁢ 3⁢ O⁢ 14

An indispensable step to understand collective magnetic phenomena in rare-earth compounds is the determination of spatially anisotropic single-ion properties resulting from spin-orbit coupling and crystal field (CF). The CF Hamiltonian has a discrete energy spectrum—accessible to spectroscopic probes such as neutron scattering—controlled by a number of independent parameters reflecting the point symmetry of the magnetic sites. Determining these parameters in low-symmetry systems is often challenging. Here, we describe a general method to analyze CF excitation spectra using adjustable effective point-charges. We benchmark our method to existing neutron-scattering measurements on pyrochlore rare-earth oxides and obtain a universal point-charge model that describes a large family of related materials. We adapt this model to the newly discovered tripod kagome magnets (𝑅 3 ⁢Mg 2 ⁢Sb 3 ⁢O 14 , 𝑅 = Tb, Ho, Er, Yb) for which we report broadband inelastic neutron-scattering spectra. Analysis of these data using adjustable point-charges yields the CF wave functions for each compound. From this, we calculate thermomagnetic properties that accurately reflect our measurements on powder samples and predict the effective gyromagnetic tensor for pseudospin degrees of freedom—a crucial step to understand the exotic collective properties of these kagome magnets at low temperature. We present further applications of our method to other tripod kagome materials and triangular rare-earth compounds 𝑅⁢MgGaO 4 (𝑅 =Yb, Tm). Overall, this study establishes a widely applicable methodology to predict CF and single-ion properties of rare-earth compounds based on interpretable and adjustable models of effective point charges.

Crystal-field theory↗

Watertight Tensor-Product Spline Reconstruction (KCNSC PDRD Technical Report)

The standard for representing geometric models in Computer Aided Design (CAD) systems is boundary representation with trimmed surfaces. However, the approximate nature of surface-surface intersections in CAD systems leaves many small gaps between mating surfaces. These gaps lead to problems for downstream analysis and manufacturing applications, requiring time consuming repair operations and further approximation of the geometry as polygonal meshes, which do not support design iteration. In this study, nVariate, Inc. evaluated its watertight reconstruction technology on a test part. The resulting watertight model has the same absolute model accuracy of the original model zero gaps, providing a representation that requires no repair and supports further design iterations in the CAD system.

42 ENGINEERING↗

Analytical and numerical modeling of optical second harmonic generation in anisotropic crystals using #SHAARP package

Electric-dipole optical second harmonic generation (SHG) is a second-order nonlinear process that is widely used as a sensitive probe to detect broken inversion symmetry and local polar order. Analytical modeling of the SHG polarimetry of a nonlinear optical material is essential to extract its point group symmetry and the absolute nonlinear susceptibilities. Current literature on SHG analysis involves numerous approximations and a wide range of (in)accuracies. We have developed an open-source package called the Second Harmonic Analysis of Anisotropic Rotational Polarimetry (#SHAARP.si) which derives analytical and numerical solutions of reflection SHG polarimetry from a single interface (.si) for bulk homogeneous crystals with arbitrary symmetry group, arbitrary crystal orientation, complex and anisotropic linear dielectric tensor with frequency dispersion, a general SHG tensor and arbitrary light polarization. #SHAARP.si enables accurate modeling of polarimetry measurements in reflection geometry from highly absorbing crystals or wedge-shaped transparent crystals. The package is extendable to multiple interfaces.

36 MATERIALS SCIENCE↗

IRMA

IRMA (In)elastic Representation of Materials As S(α,β) evaluations IRMA turns one phonon model into three outputs that usually require three separate tool chains: an evaluated nuclear-data file, predicted neutron-scattering spectra, and scattering kernels for Monte Carlo transport. The three outputs draw on a single, consistent description of the material, so the evaluation, the spectroscopy that can validate it, and the transport that uses it always agree about the physics. Nuclear data. IRMA writes ENDF-6 File 7 thermal scattering evaluations on automatically constructed (α, β) grids. This part reimplements and generalizes NJOY's LEAPR: the classic kernels reproduce freshly generated NJOY2016 tapes digit for digit and published reference tapes to about 1e-4, and the generalized paths add the exact coherent one-phonon term, anisotropic Debye-Waller tensors, coherent elastic for arbitrary crystals, and a per-species partition for polyatomic materials. The tapes feed NJOY, AMPX, FUDGE, and every transport code downstream of them. Neutron spectroscopy. The irma.spectra forward model projects the same physics onto an instrument's kinematics and resolution: INS spectra for VISION and generic indirect geometries, and 2-D S(Q,E) powder maps for direct-geometry spectrometers, from a phonopy model or straight from a phonon DOS. It can be used to predict a proposed measurement before beam time; in analysis, it supplies the calculated single-scattering counterpart of a measured spectrum, from the same material description the evaluation was built from. Monte Carlo transport. The irma.ncrystal exporter writes per-temperature scattering kernels for the companion NCrystal plugin, so McStas, OpenMC, and other NCrystal-aware codes sample the same physics. The exported kernels carry the per-site anisotropic Debye-Waller tensors, keeping directional coherent-elastic physics that NCrystal's standard scalar treatment does not represent. With the same physics inside a transport code, an entire beamline becomes a virtual experiment: IRMA's end-to-end validation ran a custom McStas implementation of the ARCS spectrometer, assembled from the existing McVine and McStas models, against measured data. From a bare crystal structure. The irma mlip front end builds the phonon model itself: a structure file and a choice of potential are enough. Nine pretrained machine-learned interatomic potentials are supported, on a laptop CPU, with no first-principles calculation; an approximate phonon model for a new material costs minutes, not a DFT campaign, and the build emits prefilled inputs for all three outputs. The result is a good starting point rather than a finished evaluation: survey-quality physics with every parameter exposed for review. A converged atomistic calculation enters the same way, as a phonopy model, when higher fidelity is needed.

Ramic, Kemal [Oak Ridge National Laboratory (ORNL)↗

3D simulations of spinlike flames in Co/Al multilayers with enhanced conduction losses

Reactive Co/Al multilayers are uniformly structured materials that may be ignited to produce rapid and localized heating. Prior studies varying the bilayer thickness (i.e., sum of two individual layers of Co and Al) have revealed different types of flame morphologies, including: (a) steady/planar, (b) wavy/periodic, and (c) transverse bands, originating in the flame front. These instabilities resemble the “spin waves” first observed in the early studies of solid combustion (i.e., Ti cylinder in a N 2 atmosphere), and are likewise thought to be due to the balance of heat released by reaction and heat conduction forward into the unreacted multilayer. However, the multilayer geometry and three-dimensional (3D) edge effects are relatively unexplored. In this work, a new diffusion-limited reaction model for Co/Al multilayers was implemented in large, novel 3D finite element analysis (FEA) simulations, in order to study the origins of these spinlike flames. This reaction model builds upon previous work by introducing three new phase-dependent property models for: (1) the diffusion coefficient, (2) anisotropic thermal conductivity tensor, and (3) bulk heat capacity, as well as one additional model for the bilayer-dependent heat of reaction. These novel 3D simulations are the first to predict both steady and unsteady flames in Co/Al multilayers. Moreover, two unsteady modes of flame propagation are identified, which depend on the enhanced conduction losses with slower flames, as well as flame propagation around notched edges. Finally, future work will consider the generality of the current modeling approach and also seek to define a more generalized set of stability criteria for additional multilayer systems.

42 ENGINEERING↗

Quantum field theory of topological spin dynamics

Here, we develop a field theory of quantum magnets and magnetic (semi)metals, which is suitable for the analysis of their universal and topological properties. The systems of interest include collinear, coplanar, and general noncoplanar magnets. At the basic level, we describe the dynamics of magnetic moments using smooth vector fields in the continuum limit. Dzyaloshinskii-Moriya interaction is captured by a non-Abelian vector gauge field, and chiral spin couplings related to topological defects appear as higher-rank antisymmetric tensor gauge fields. We distinguish type-I and type-II magnets by their equilibrium response to the non-Abelian gauge flux, and characterize the resulting lattices of skyrmions and hedgehogs, the spectra of spin waves, and the chiral response to external perturbations. The general spin-orbit coupling of electrons is similarly described by non-Abelian gauge fields, including higher-rank tensors related to the electronic Berry flux. Itinerant electrons and local moments exchange their gauge fluxes through Kondo and Hund interactions. Hence, by utilizing gauge fields, this theory provides a unifying physical picture of “intrinsic” and “topological” anomalous Hall effects, spin-Hall effects, and other correlations between the topological properties of electrons and moments. We predict “topological” magnetoelectric effect in materials prone to hosting hedgehogs. Links to experiments and model calculations are provided by deriving the couplings and gauge fields from generic microscopic models, including the Hubbard model with spin-orbit interactions. Much of the formal analysis is generalized to d spatial dimensions in order to access the $π_{d–1}(S^{d–1})$ homotopy classification of the magnetic hedgehog topological defects, and establish the possibility of novel quantum spin liquids that exhibit a fractional magnetoelectric effect. However, we emphasize the form of all results in the physically relevant $\textit{d}$ = 3 dimensions, and discuss a few applications to topological magnetic conductors like Mn 3 Sn and Pr 2 Ir 2 O 7 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Reconstructing f ( T ) gravity and exploring the torsion driven warm inflationary cosmology

The current paper reports an investigation of a warm inflationary scenario in the context of f(T) gravity for a spatially flat FLRW universe. In our model, inflation is driven purely by the torsional sector of f(T) gravity, without introducing any additional scalar fields. We focus on the high dissipative regime (R >> 1), reconstruct the Hubble parameter as a function of the e-folding number N, and derive the slow-roll parameters ε 1 (N) and ε 2 (N). The study has encapsulated the dynamics of inflation and its duration under strong dissipation. The dissipative coefficient Γ is modeled with a temperature-dependent power-law form, linking the inflationary dynamics to thermal corrections and the particle content of the early universe. The analysis has affirmed that the torsion-induced energy density ρ T successfully transitions to radiation energy density ρ rad , facilitating a graceful exit from inflation. Finally, we have validated our model by comparing the scalar spectral index and tensor-to-scalar ratio with Planck 2018 results, demonstrating consistency within observational bounds. Additionally, it is verified that the thermal domination condition T * /H > 1 and the torsion dominance condition ρ T /ρ rad > 1 are satisfied.

Ghosh, Moli [Amity University, Kolkata (India); Mr↗

Measuring residual stress in Ti-6Al-4V with HR-EBSD, using reference patterns from annealed material

In this study, elastic strain and residual stress values in an engineering material were determined by high-resolution EBSD (HR-EBSD) using reference patterns from annealed, recrystallized material. Specifically, rolled Ti-6Al-4V sheet was analyzed using HR EBSD in both an as-annealed and a strained state, the latter established via a static two-point bend fixture. The commercial software CrossCourt4 (CC4) was used to evaluate the relative intragranular strain tensor components of the alpha grains mapped by HR-EBSD. The intragranular reference patterns identified in each grain by the software were then replaced with similarly-oriented patterns from scans of the as-annealed, recrystallized material. The CC4 analysis was then repeated to provide calibrated strain and residual stress values. Elastic strains on the order of 10–3 were detected, with corresponding residual stress values on the order of 1 GPa. Subsequent analyses using high energy synchrotron x-ray diffraction showed good correlation with the elastic strains measured by HR-EBSD using the annealed reference patterns. Additionally, the maximum macroscopic stress analytically predicted for two-point bending was in agreement with the residual stress observed. Good repeatability of the methodology and insight into the state of individual grains were also observed. The findings provide encouraging evidence that accurate intragranular and long-range residual stress measurements of engineering significance can be made with HR-EBSD using a physical reference pattern library.

36 MATERIALS SCIENCE↗

EGS Stimulation Design with Uncertainty Quantification at the EGS Collab Site

Engineering a robust hydraulic connection between wells is one of the most difficult aspects of enhanced geothermal systems (EGS). Designing and constructing such hydraulic connections requires and understanding of the in situ state of stress and the heterogeneities and discontinuities that naturally exist and may control the stimulation. Even with comprehensive stress and formation characterization programs substantial uncertainty remains in these key parameters. This is especially the case in high-temperature EGS environments where drilling conditions are often difficult and a far fewer logging and testing options are available. This paper presents a new approach for explicitly quantifying the uncertainties in the state of stress using a Bayesian Markov Chain Monte Carlo method. This approach produces a probability distribution for the stress tensor, including a general 3D orientation, that reflects the uncertainties in all the observations or indicators used to constrain the stress state. This method is demonstrated on the characterization data for the EGS Collab Experiment 2 site. The output of the analysis is used to guide the design of the planned stimulations. In the case of research projects like EGS Collab, explicitly quantifying the uncertainties in the stress state allow for more rigorous hypothesis testing by allowing conclusions drawn from the experiments to be interpreted in the context of the uncertain knowledge about conditions in the test bed.

Burghardt, Jeffrey A.↗

Rheological Properties of Small-Molecular Liquids at High Shear Strain Rates

Molecular-scale understanding of rheological properties of small-molecular liquids and polymers is critical to optimizing their performance in practical applications such as lubrication and hydraulic fracking. We combine nonequilibrium molecular dynamics simulations with two unsupervised machine learning methods: principal component analysis (PCA) and t-distributed stochastic neighbor embedding (t-SNE), to extract the correlation between the rheological properties and molecular structure of squalane sheared at high strain rates (10 6 –10 10 s -1 ) for which substantial shear thinning is observed under pressures P ϵ 0.1–955 MPa at 293 K. Intramolecular atom pair orientation tensors of 435 × 6 dimensions and the intermolecular atom pair orientation tensors of 61 × 6 dimensions are reduced and visualized using PCA and t-SNE to assess the changes in the orientation order during the shear thinning of squalane. Dimension reduction of intramolecular orientation tensors at low pressures P = 0.1,100 MPa reveals a strong correlation between changes in strain rate and the orientation of the side-backbone atom pairs, end-backbone atom pairs, short backbone-backbone atom pairs, and long backbone-backbone atom pairs associated with a squalane molecule. At high pressures P ≥ 400 MPa, the orientation tensors are better classified by these different pair types rather than strain rate, signaling an overall limited evolution of intramolecular orientation with changes in strain rate. Dimension reduction also finds no clear evidence of the link between shear thinning at high pressures and changes in the intermolecular orientation. The alignment of squalane molecules is found to be saturated over the entire range of rates during which squalane exhibits substantial shear thinning at high pressures.

36 MATERIALS SCIENCE↗