Search NASA⌕ Search

SEARCH · Search NASA

Results for “general tensor structure”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Spin-orbit correlations in the nucleon in the large- N c limit

We study the twist-3 spin-orbit correlations of quarks described by the nucleon matrix elements of the parity-odd rank-2 tensor QCD operator (the parity-odd partner of the QCD energy-momentum tensor). Our treatment is based on the effective dynamics emerging from the spontaneous breaking of chiral symmetry and the mean-field picture of the nucleon in the large- N c limit. The twist-3 QCD operators are converted to effective operators, in which the QCD interactions are replaced by spin-flavor-dependent chiral interactions of the quarks with the pion field. We compute the nucleon matrix elements of the twist-3 effective operators and discuss the role of the chiral interactions in the spin-orbit correlations. We derive the first-quantized representation in the mean-field picture and develop a quantum-mechanical interpretation. The chiral interactions give rise to new spin-orbit couplings and qualitatively change the correlations compared to the quark model picture. We also derive the twist-3 matrix elements in the topological soliton picture where the quarks are integrated out (skyrmion). The methods used here can be extended to other QCD operators describing higher-twist nucleon structure and generalized parton distributions. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Multiplexing Focusing Analyzer for Efficient Stress-Strain Measurements

Statement of the problem or situation that is being addressed. Although thermal and cold neutron scattering is widely used and is critical for success in many areas of materials science and engineering, relatively low neutron fluxes severely limit applications of not only laboratory neutrons generators, but also large national neutron facilities. State-of-the-art thermal and cold neutron sources are large expensive national facilities, which serve diverse community of scientific and industrial users. The constant need to improve the instruments performance, stems from the fact that neutron methods are gaining in popularity, and becoming more and more powerful, while new neutron sources are not being constructed to keep pace with the developments and needs of the scientific community. Small research reactors at universities and National Labs, and laboratory-based neutron generators, are necessary not only for education and training, but also when samples cannot be transported to other facilities. However, the standard neutron techniques, which were developed for high-flux facilities, require much higher efficiencies to be used effectively with the low fluxes of small sources. Thus, the efficient use of neutron sources, such as with our proposed analyzer, is important for the progress and broader use of these neutron techniques. General statement of how this problem is being addressed. We propose to design and demonstrate novel diffractive optical device, which will enable very efficient residual stress neutron diffractometers. The proposed device will be a multi-foil analyzer, where each foil is constructed of focusing bent single crystals of Si. Such device will enable polychromatic residual stress neutron diffraction. At large national facilities, such as at Oak Ridge National Laboratory, these analyzers would enable very fast measurements for determining residual stress tensors, raster large samples or screen multiple samples. Commercial Applications and Other Benefits The outcome of this project would be the demonstration of commercial devices, novel neutron optical components, which could be utilized to improve the performance of existing instruments or build novel neutron scattering instruments at DOE neutron facilities and commercial laboratory neutron sources. These new devices will widen the scope of research conducted using neutrons and enable measurements not feasible at present. Summary for Members of Congress Thermal and cold neutron beams are a powerful materials science probe, which provide unique information about the structure of matter. The proposed innovations expand the reach of neutron-based investigations to new materials and industries by enabling new instrumentation capabilities, thereby greatly enhancing and expanding the role of small, laboratory-based neutron instrumentation, and improving education and training of neutron users.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Angle-resolved polarized Raman study of layered b-As x P x-1 alloys: Identification of As-P vibrational modes

In this study, the polarization-resolved Raman spectra of b-As x P x-1 flakes with varying arsenic concentration (x = 0, 0.4, and 0.8) are systematically investigated. A clear polarization dependence is observed in the Raman intensity of all the P-P, As-P, and As-As modes and confirmed by the first-principles calculations. The observed angle-dependence can be related to the structure of the Raman tensors of each peak. In particular, the Raman intensities corresponding to As-P vibration modes are studied in detail. The out of plane As-P mode ($A$$^{1}_{g}$) is clearly identified, but the butterfly-shaped behavior observed in the polar plot around 350 cm -1 is interpreted as the combination of the two in-plane vibrations, one along the armchair and the other along the zigzag direction. This work allows unambiguous assignment of the $A$$^{1}_{g}$, B 2g , and $A$$^{2}_{g}$ modes of the As-P vibrational modes and the importance of the consideration of the superposition of the in-plane B 2g and $A$$^{2}_{g}$ modes. Such study can be utilized to further characterize the anisotropic structural properties and analyze the local distribution of the As and P in b-As x P x-1 alloys which will play key roles in electronic and optical properties. Furthermore, this methodology can be further extended to alloyed two-dimensional material systems with in-plane anisotropy in general.

2D materials↗
DOE OSTI2025

Large Hyperfine Coupling Arising from Pseudo- 2 S Ground States in a Series of Lutetium(II) Metallocene Complexes

The synthesis of molecules with strong coupling between electronic and nuclear spins represents an important challenge in molecular quantum information science. Here, we report the synthesis and characterization of the divalent lutetium metallocene complexes Lu(Cp Me 5 )(Cp iPr 5 ) (Cp Me 5 = pentamethylcyclopentadienyl; Cp iPr 5 = pentaisopropylcyclopentadienyl), Lu(Cp iPr 4 Et ) 2 (Cp iPr 4 Et = ethyltetraisopropylcyclopentadienyl), and Lu(Cp iPr 4 ) 2 (Cp iPr 4 = tetraisopropylcyclopentadienyl). The molecular structures of these complexes, as determined through singlecrystal X-ray diffraction, feature a common bent sandwich geometry, with average Cp–Lu–Cp angles ranging from 159.9° to 152.6°. Analysis of continuous-wave electron paramagnetic resonance (EPR) spectra for the complexes reveals nearly isotropic g tensors with only a slight deviation from that of a free electron. Moreover, an extremely large splitting of the eight-line spectra indicates the presence of strong hyperfine coupling, and simulations provide isotropic hyperfine coupling constants of A iso = 4.38, 4.30, and 4.17 GHz across the series, where the value of A iso is found to decrease as the Cp–Lu–Cp angle becomes more acute. Notably, these values are the largest yet observed for any lanthanide complex. Moreover, EPR and computational analysis show that the large values of A iso stem from large s-orbital character up to 41.2% in the corresponding singly occupied molecular orbitals. To our knowledge, this degree of s-character in a molecular orbital is the largest yet reported for an open-shell isolable complex. These results outline a general strategy toward the isolation of paramagnetic molecules with strong hyperfine coupling and highly isotropic doublet electronic ground states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Quantum mereology: Factorizing Hilbert space into subsystems with quasiclassical dynamics

We study the question of how to decompose Hilbert space into a preferred tensor-product factorization without any preexisting structure other than a Hamiltonian operator, in particular the case of a bipartite decomposition into “system” and “environment.” Such a decomposition can be defined by looking for subsystems that exhibit quasiclassical behavior. The correct decomposition is one in which pointer states of the system are relatively robust against environmental monitoring (their entanglement with the environment does not continually and dramatically increase) and remain localized around approximately classical trajectories. We present an in-principle algorithm for finding such a decomposition by minimizing a combination of entanglement growth and internal spreading of the system. Both of these properties are related to locality in different ways. Furthermore, this formalism is relevant to questions in the foundations of quantum mechanics and the emergence of spacetime from quantum entanglement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1993

On the large eddy simulation of turbulent flows in complex geometry

Application of the method of Large Eddy Simulation (LES) to a turbulent flow consists of three separate steps. First, a filtering operation is performed on the Navier-Stokes equations to remove the small spatial scales. The resulting equations that describe the space time evolution of the 'large eddies' contain the subgrid-scale (sgs) stress tensor that describes the effect of the unresolved small scales on the resolved scales. The second step is the replacement of the sgs stress tensor by some expression involving the large scales - this is the problem of 'subgrid-scale modeling'. The final step is the numerical simulation of the resulting 'closed' equations for the large scale fields on a grid small enough to resolve the smallest of the large eddies, but still much larger than the fine scale structures at the Kolmogorov length. In dividing a turbulent flow field into 'large' and 'small' eddies, one presumes that a cut-off length delta can be sensibly chosen such that all fluctuations on a scale larger than delta are 'large eddies' and the remainder constitute the 'small scale' fluctuations. Typically, delta would be a length scale characterizing the smallest structures of interest in the flow. In an inhomogeneous flow, the 'sensible choice' for delta may vary significantly over the flow domain. For example, in a wall bounded turbulent flow, most statistical averages of interest vary much more rapidly with position near the wall than far away from it. Further, there are dynamically important organized structures near the wall on a scale much smaller than the boundary layer thickness. Therefore, the minimum size of eddies that need to be resolved is smaller near the wall. In general, for the LES of inhomogeneous flows, the width of the filtering kernel delta must be considered to be a function of position. If a filtering operation with a nonuniform filter width is performed on the Navier-Stokes equations, one does not in general get the standard large eddy equations. The complication is caused by the fact that a filtering operation with a nonuniform filter width in general does not commute with the operation of differentiation. This is one of the issues that we have looked at in detail as it is basic to any attempt at applying LES to complex geometry flows. Our principal findings are summarized.

Ghosal, Sandip↗
DOE OSTI2022

Polar liquids at charged interfaces: A dipolar shell theory

The structure of polar liquids and electrolytic solutions, such as water and aqueous electrolytes, at interfaces underlies numerous phenomena in physics, chemistry, biology, and engineering. In this work, we develop a continuum theory that captures the essential features of dielectric screening by polar liquids at charged interfaces, including decaying spatial oscillations in charge and mass, starting from the molecular properties of the solvent. The theory predicts an anisotropic dielectric tensor of interfacial polar liquids previously studied in molecular dynamics simulations. We explore the effect of the interfacial polar liquid properties on the capacitance of the electrode/electrolyte interface and on hydration forces between two plane-parallel polarized surfaces. In the linear response approximation, we obtain simple formulas for the characteristic decay lengths of molecular and ionic profiles at the interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Nanoscale lattice strains in self-ion implanted tungsten

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

36 MATERIALS SCIENCE↗
DOE OSTI2022

Elastic Solutions to 2D Plane Strain Problems: Nonlinear Contact and Settlement Analysis for Shallow Foundations

The classical Neumann boundary value problem of an isotropic, homogeneous elastic half-plane under plane strain conditions is readdressed as the limiting case of the fully three-dimensional problem. Analytical solutions of the stress and strain tensors are obtained by taking the limit from known three-dimensional solutions. It is shown that the displacement fields for the plane strain problem are not well defined. A small number of simple expressions are developed, which provide a general solution for linearly-varying traction over arbitrary regions on the boundary. A simple, efficient, and rapidly convergent algorithm is developed which uses these solutions as analytic elements and provides a solution approach to the general boundary value problem. The method is verified against known solutions for Hertzian contact between parallel cylinders. Two numerical examples are presented for the analysis of shallow foundation systems. In the first, the boundary conditions are informed by analytical elastoplastic calculations and a strain influence analysis is performed and compared with the Schmertmann method. Subsequently, empirical laboratory contact traction distributions measured by Bauer et al., in both the normal and tangential directions are employed as boundary conditions for an analysis of the underlying stress field.

42 ENGINEERING↗
DOE OSTI2026

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)↗
DOE OSTI2022

Models of $J/\varPsi$ photo-production reactions on the nucleon

The J/Ψ photo-production reactions on the nucleon can provide information on the roles of gluons in determining the J/Ψ-nucleon (J/Ψ-N) interactions and the structure of the nucleon. The information on the J/Ψ-N interactions is needed to test lattice QCD (LQCD) calculations and to understand the nucleon resonances such as N*(P c ) recently reported by the LHCb Collaboration. In addition, it is also needed to investigate the production of nuclei with hidden charms and to extract the gluon distributions in nuclei. The main purpose of this article is to review six models of the γ+p→J/Ψ+p reaction which have been and can be applied to analyze the data from Thomas Jefferson National Accelerator Facility (JLab). Here, the formulae for each model are given and used to obtain the results to show the extent to which the available data can be described. The models presented include the Pomeron-exchange model of Donnachie and Landshoff (Pom-DL) and its extensions to include J/Ψ-N potentials extracted from LQCD (Pom-pot) and to also use the constituent quark model (CQM) to account for the quark substructure of J/Ψ (Pom-CQM). The other three models are developed from applying the perturbative QCD approach to calculate the two-gluon exchange using the generalized parton distribution (GPD) of the nucleon (GPD-based), two- and three-gluon exchanges using the parton distribution of the nucleon (2g+3g), and the exchanges of scalar (0 ++ ) and tensor (2 ++ ) glueballs within the holographic formulation (holog). The results of investigating the excitation of the nucleon resonances N*(P c ) in the γ+p→J/Ψ+p reactions are also given. We demonstrate that the differences between these six models can be unambiguously distinguished and the N* can be better studied by using the forthcoming JLab data at large |t| and at energies very near the J/Ψ production threshold. Possible improvements of the considered models are discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

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↗
DOE OSTI2022

Analytic error analysis of cross section interpolation methods in nodal diffusion codes - I : Theory

This paper discusses two cross section interpolation methods commonly found in popular nodal codes; the partial derivatives and multiple tables models. The motivation for choosing a model, and thus a case matrix structure, is a trade off between accuracy and computational cost. Due to decades of experience, there are default structures that are sufficient for current light water reactor analysis. However, this is not necessarily the case for advanced reactor designs. Therefore, it is advantageous to understand the sources of error in cross section interpolation models so that the quality of a case matrix may be improved. A mathematical framework for these models is presented in this work that provides a more rigorous connection between the nuclear engineering field's cross section interpolation methods and the broader mathematical field of function approximation. The two cross section models examined in this paper were found to utilize Lagrange interpolation and are a subset of Lagrange tensor products. Classical results of Lagrange polynomial error analysis were then applied to the partial derivative and multiple tables models to derive expressions for the total point-wise error. The analytical results classify the total error into two parts: the model form error and interpolation error. Finally, based on our observations, a better foundation for improving the quality of a case matrix is proposed. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Mapping the intrinsic photocurrent streamlines through micromagnetic heterostructure devices

Photocurrent in quantum materials is often collected at global contacts far away from the initial photoexcitation. This collection process is highly nonlocal. It involves an intricate spatial pattern of photocurrent flow (streamlines) away from its primary photoexcitation that depends sensitively on the configuration of current collecting contacts as well as the spatial nonuniformity and tensor structure of conductivity. Direct imaging to track photocurrent streamlines is challenging. Here, we demonstrate a microscopy method to image photocurrent streamlines through ultrathin heterostructure devices comprising platinum on yttrium iron garnet (YIG). We accomplish this by combining scanning photovoltage microscopy with a uniform rotating magnetic field. Here, local photocurrent is generated through a photo-Nernst type effect with its direction controlled by the external magnetic field. This enables the mapping of photocurrent streamlines in a variety of geometries that include conventional Hall bar-type devices, but also unconventional wing-shaped devices called electrofoils. In these, we find that photocurrent streamlines display contortion, compression, and expansion behavior depending on the shape and angle of attack of the electrofoil devices, much in the same way as tracers in a wind tunnel map the flow of air around an aerodynamic airfoil. This affords a powerful tool to visualize and characterize charge flow in optoelectronic devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Physics-informed machine learning of the Lagrangian dynamics of velocity gradient tensor

Reduced models describing the Lagrangian dynamics of the velocity gradient tensor (VGT) in homogeneous isotropic turbulence (HIT) are developed under the physics-informed machine learning (PIML) framework. We consider the VGT at both Kolmogorov scale and coarse-grained scale within the inertial range of HIT. Building reduced models requires resolving the pressure Hessian and subfilter contributions, which is accomplished by constructing them using the integrity bases and invariants of the VGT. The developed models can be expressed using the extended tensor basis neural network (TBNN) introduced by Ling et al. [J. Fluid Mech. 807, 155 (2016)]. Physical constraints, such as Galilean invariance, rotational invariance, and incompressibility condition, are thus embedded in the models explicitly. Our PIML models are trained on the Lagrangian data from a high-Reynolds number direct numerical simulation (DNS). To validate the results, we perform a comprehensive out-of-sample test. We observe that the PIML model provides an improved representation for the magnitude and orientation of the small-scale pressure Hessian contributions. Statistics of the flow, as indicated by the joint PDF of second and third invariants of the VGT, show good agreement with the “ground-truth” DNS data. A number of other important features describing the structure of HIT are reproduced by the model successfully. We have also identified challenges in modeling inertial range dynamics, which indicates that a richer modeling strategy is required. This helps us identify important directions for future research, in particular towards including inertial range geometry into the TBNN.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Universal Tripartite Entanglement in One-Dimensional Many-Body Systems

Motivated by conjectures in holography relating the entanglement of purification and reflected entropy to the entanglement wedge cross section, we introduce two related non-negative measures of tripartite entanglement g and h. We prove structure theorems which show that states with nonzero g or h have nontrivial tripartite entanglement. We then establish that in one dimension these tripartite entanglement measures are universal quantities that depend only on the emergent low-energy theory. For a gapped system, we argue that either g≠0 and h=0 or g=h=0, depending on whether the ground state has long-range order. For a critical system, we develop a numerical algorithm for computing g and h from a lattice model. We compute g and h for various CFTs and show that h depends only on the central charge whereas g depends on the whole operator content.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Entanglement structures in quantum field theories: Negativity cores and bound entanglement in the vacuum

Here, the many-body entanglement between two finite (size-d) disjoint vacuum regions of noninteracting lattice scalar field theory in one spatial dimension, i.e., a (d A × d B ) mixed Gaussian continuous variable system, is locally transformed into a tensor-product core of (1 A × 1 B ) mixed entangled pairs. Accessible entanglement within these core pairs exhibits an exponential hierarchy and as such identifies the structure of dominant region modes from which vacuum entanglement could be extracted into a spatially separated pair of quantum detectors. Beyond the core, the remaining modes of the halo are determined to be AB separable in isolation, as well as separable from the core. However, state preparation protocols that distribute entanglement in the form of (1 A × 1 B ) mixed core pairs are found to require additional entanglement in the halo that is obscured by classical correlations. This inaccessible (bound) halo entanglement is found to mirror the accessible entanglement, but with a step behavior as the continuum is approached. It remains possible that alternate initialization protocols that do not utilize the exponential hierarchy of core-pair entanglement may require less inaccessible entanglement. Entanglement consolidation is expected to persist in higher dimensions and may aid classical and quantum simulations of asymptotically free gauge field theories, such as quantum chromodynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Geothermal Fault Zone and Fluid Imaging through Joint Airborne ZTEM and Ground MT Data Inversion Analysis

This project has aimed to achieve detailed electrical resistivity resolution at geothermal reservoir scales by combining airborne natural electromagnetic (EM) field surveying (ZTEM) with ground magnetotelluric (MT) measurements to approximate an airborne MT geophysical method. MT alone is relatively expensive and may have permitting challenges in sensitive areas. Airborne ZTEM field data contains only the magnetic field, requires a background assumption, and has been limited to relatively high frequencies, thus suffering uniqueness problems. Based on proto-type 2D simulations, ZTEM ambiguities may be reduced through formal incorporation with possibly sparse ground MT soundings, which we pursued in full 3D for this project. The methodology was tested at the high-temperature Roosevelt Hot Springs geothermal system, Utah, which was considered advantageous given the near total exposure of crystalline reservoir rocks across the project area. ZTEM and ground MT survey data were acquired in 2017, subcontracted to outside parties with which we have worked in the past. These included 80 remote-referenced tensor MT soundings over the Mineral Mountains and adjacent Roosevelt Hot Spring producing geothermal system. These MT stations abut later coverage of a similar number of MT stations taken for the Utah FORGE project providing excellent total data aperture to re-solve structure beneath both project areas better than either set alone. The airborne ZTEM survey covered 704 line kilometers in E-W flight lines with a 250 m line spacing. Although this survey was timed during a maintenance-related shutdown of power production at the Roosevelt Hot Springs, other noise sources difficult to identify but including two high-voltage state-scale transmission lines compromised the ZTEM survey badly leading to unusable responses. Thus, with DOE management concurrence, the project proceeded to emphasize inversion and interpretation of the joint SubTER-FORGE MT data sets with regard to the Roosevelt Hot Springs reservoir recharge and to deep heat sources for both it and the Utah FORGE EGS project area. We also investigated the joint ZTEM-MT sampling concept with data sets from the Eleven Mile Canyon prospect area donated by the U.S. Navy (A. Sabin, PoC). Inversion of the SubTER-FORGE MT data using the HexMT 3D finite element algorithm reveals a large, low-resistivity anomaly extending sub-vertically through the depth range of the crust beneath the western Mineral Mountains. The steep conductive zone connects in the lower crust to a more tabular conductor characteristic of much of the Great Basin that generally is ascribed to current mafic magmatic underplating, hybridization and fluid release. The location of the resolved anomaly relative to the recent (0.5-0.8 Ma) eruptive centers of the Mineral Mountains implicates it as remnants of the magma body which fed these centers. This structure appears to be currently feeding heat and fluids upward into the Roosevelt Hot Springs hydrothermal system, as well as heat laterally to the FORGE project area. Separate and joint inversion models were carried out for the donated Eleven Mile Canyon MT-ZTEM data set to demonstrate concept. ZTEM only inversion showed two main alteration zones in the western portion of the project area known from geological mapping. Joint inversion including an E-W profile of MT soundings sharpened these features considerably. It also resolved in much greater detail the graben related normal faulting structure of the central project area which lies at depths exceeding the sensitivity of ZTEM alone. The sparse number of MT da-ta relative to the ZTEM required upweighting the former by a factor of several, but an exact procedure awaits future research. Our final impression is that sparse MT data can improve resolution of the subsurface over that of ZTEM alone. However, well sampled MT data are to be preferred and offer the simplicity of interpreting just one data type, and possess the superior resolution capability coming with the electric field everywhere, and from their high bandwidth.

15 GEOTHERMAL ENERGY↗