Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix Theory”

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 1,783 records · Page 99
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

A time dependent difference theory for sound propagation in ducts with flow

A time dependent numerical solution of the linearized continuity and momentum equation is developed for sound propagation in a two-dimensional straight hard or soft wall duct with a sheared mean flow. The time dependent governing acoustic-difference equations and boundary conditions are developed along with a numerical determination of the maximum stable time increments. The analysis begins with a harmonic noise source radiating into a quiescent duct. This explicit iteration method then calculates stepwise in real time to obtain the transient as well as the 'steady' state solution of the acoustic field. Example calculations are presented for sound propagation in hard and soft wall ducts, with no flow and with plug flow. Although the problem with sheared flow has been formulated and programmed, sample calculations have not yet been examined. So far, the time dependent finite difference analysis has been found to be superior to the steady state finite difference and finite element techniques because of shorter solution times and the elimination of large matrix storage requirements.

Baumeister, K. J.↗
DOE OSTI2025

Kaon gluon parton distribution and momentum fraction from 2+1+1 lattice QCD with high statistics

We present a high-statistics lattice-QCD determination of the kaon gluon parton distribution function and gluon momentum fraction. We use clover valence fermion action to take 1,296,640 kaon-correlator measurements on a highly improved staggered quark ensemble with 𝑎 ≈ 0.12 fm and 310-MeV pion mass generated by the MILC Collaboration. A detailed investigation into the impact of gauge-link smearing on the gluonic matrix elements indicates that five steps of hypercubic smearing offer an effective balance between signal quality and preservation of long-distance physics. We report a nonperturbatively renormalized kaon gluon momentum fraction of Math output error at 𝜇 = 2 GeV in the Math output error scheme. Using reduced pseudo-Ioffe-time-distribution matrix elements and pseudo-parton-distribution-function (PDF) matching, we extract the kaon gluon PDF and compare with the prediction from the Dyson-Schwinger equation and with the pion PDF obtained from the same ensemble.

First-principles calculations↗
DOE OSTI2020

Piezoelectric Electron-Phonon Interaction from Ab Initio Dynamical Quadrupoles: Impact on Charge Transport in Wurtzite GaN

First-principles calculations of $\textit{e}$-ph interactions are becoming a pillar of electronic structure theory. However, the current approach is incomplete. The piezoelectric (PE) $\textit{e}$-ph interaction, a long-range scattering mechanism due to acoustic phonons in noncentrosymmetric polar materials, is not accurately described at present. Current calculations include short-range $\textit{e}$-ph interactions (obtained by interpolation) and the dipolelike Frölich long-range coupling in polar materials, but lack important quadrupole effects for acoustic modes and PE materials. Here we derive and compute the long-range $\textit{e}$-ph interaction due to dynamical quadrupoles, and apply this framework to investigate $\textit{e}$-ph interactions and the carrier mobility in the PE material wurtzite GaN. In this work, we show that the quadrupole contribution is essential to obtain accurate $\textit{e}$-ph matrix elements for acoustic modes and to compute PE scattering. Our work resolves the outstanding problem of correctly computing $\textit{e}$-ph interactions for acoustic modes from first principles, and enables studies of $\textit{e}$-ph coupling and charge transport in PE materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Deformed Fredkin model for the ν = 5 / 2 Moore-Read state on thin cylinders

We propose a frustration-free model for the Moore-Read quantum Hall state on sufficiently thin cylinders with circumferences ≲ 7 magnetic lengths. While the Moore-Read Hamiltonian involves complicated long-range interactions between triplets of electrons in a Landau level, our effective model is a simpler one-dimensional chain of qubits with deformed Fredkin gates. We show that the ground state of the Fredkin model has high overlap with the Moore-Read wave function and accurately reproduces the latter's entanglement properties. Moreover, we demonstrate that the model captures the dynamical response of the Moore-Read state to a geometric quench, induced by suddenly changing the anisotropy of the system. We elucidate the underlying mechanism of the quench dynamics and show that it coincides with the linearized bimetric field theory. The minimal model introduced here can be directly implemented as a first step towards quantum simulation of the Moore-Read state, as we demonstrate by deriving an efficient circuit approximation to the ground state and implementing it on an IBM quantum processor. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Low-dimensional de Sitter quantum gravity

We study aspects of Jackiw-Teitelboim (JT) quantum gravity in two-dimensional nearly de Sitter (dS) spacetime, as well as pure de Sitter quantum gravity in three dimensions. These are each theories of boundary modes, which include a reparameterization field on each connected component of the boundary as well as topological degrees of freedom. In two dimensions, the boundary theory is closely related to the Schwarzian path integral, and in three dimensions to the quantization of coadjoint orbits of the Virasoro group. Using these boundary theories we compute loop corrections to the wavefunction of the universe, and investigate gravitational contributions to scattering. Along the way, we show that JT gravity in dS2 is an analytic continuation of JT gravity in Euclidean AdS2, and that pure gravity in dS3 is a continuation of pure gravity in Euclidean AdS3. We define a genus expansion for de Sitter JT gravity by summing over higher genus generalizations of surfaces used in the Hartle-Hawking construction. Assuming a conjecture regarding the volumes of moduli spaces of such surfaces, we find that the de Sitter genus expansion is the continuation of the recently discovered AdS genus expansion. Then both may be understood as coming from the genus expansion of the same double-scaled matrix model, which would provide a non-perturbative completion of de Sitter JT gravity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2007

Lightning Charge Retrievals: Dimensional Reduction, LDAR Constraints, and a First Comparison with LIS Satellite Data

A "dimensional reduction" ("DR") method is introduced for analyzing lightning field changes ((Delta)Es) whereby the number of unknowns in a discrete two-charge model is reduced from the standard eight (x, y, z, Q, x', y', z', Q') to just four (x, y, z, Q). The four unknowns (x, y, z, Q) are found by performing a numerical minimization of a chi-square function. At each step of the minimization, an overdetermined fixed matrix (OFM) method is used to immediately retrieve the best "residual source" (x', y', z', Q'), given the values of (x, y, z, Q). In this way, all eight parameters Ix. y, z, Q, x', y', z', Q') are found, yet a numerical search of only four parameters (x, y, z, Q) is required. The DR method has been used to analyze lightning-caused (Delta)Es derived from multiple ground-based electric field measurements at the NASA Kennedy Space Center (KSC) and U.S. Air Force Eastern Range (ER). The accuracy of the DR method has been assessed by comparing retrievals with data provided by the lightning detection and ranging (LDAR) system at the KSC-ER. and from least squares error estimation theory, and the method is shown to be a useful "stand alone" charge retrieval tool. Since more than one charge distribution describes a finite set of (Delta)Es (i.e., solutions are nonunique), and since there can be appreciable differences in the physical characteristics of these solutions, not all DR solutions are physically acceptable. Hence. an alternative and more accurate method of analysis is introduced that uses LDAR data to constrain the geometry of the charge solutions. thereby removing physically unacceptable retrievals. The charge solutions derived from this method are shown to compare well with independent satellite- and ground-based observations of lightning m several Florida storms.

Koshak, W. J.↗
DOE OSTI2023

Two-dimensional patterning of mesoscale fibers using acoustophoresis

The performance of functional composites can rely critically on the arrangement of secondary phases; for example, patterned networks of conductive particles can impart anisotropic thermal, electric or ionic conductivity while preserving flexibility in the matrix. We demonstrate the use of standing acoustic waves to generate periodic patterns of short fibers. We extend the range of possible patterns with the first demonstration of both rectangular grids and arrays of octagons interspersed with rectangles. These newly demonstrated patterns are rationalized using theoretical models of acoustic forces and torques on fibers that account for two-dimensional spatial variations arising from applied acoustic fields. The models enable simulations of fiber motion, which are used to (i) map out final fiber positions as a function of initial position and orientation, and (ii) corroborate experiments visualizing fiber motion and final patterns. This approach provides a fast and accurate way to predict emergent fiber patterns as a function of excitation modality and fiber length. The theory and experiments clearly indicate strong coupling between the length of the fibers and the spacing of the acoustic nodes. Furthermore, this coupling is used to estimate reductions in percolation thresholds associated with the ratio of fiber length and acoustic wavelength.

25 ENERGY STORAGE↗
DOE OSTI2021

Robust d-Wave Superconductivity in the Square-Lattice t–J Model

Unravelling competing orders emergent in doped Mott insulators and their interplay with unconventional superconductivity is one of the major challenges in condensed matter physics. Here, to explore the possible superconducting state in a doped Mott insulator, we study the square-lattice t-J model with both the nearest-neighbor and next-nearest-neighbor electron hoppings and spin interactions. By using the state-of-the-art density matrix renormalization group calculation with imposing charge U(1) and spin SU(2) symmetries on the six-leg cylinders, we establish a quantum phase diagram including three phases: a stripe charge density wave phase, a superconducting phase without static charge order, and a superconducting phase coexistent with a weak charge stripe order. Crucially, we demonstrate that the superconducting phase has a power-law pairing correlation that decays much slower than the charge density and spin correlations, which is a quasi-1D descendant of the uniform d-wave superconductor in two dimensions. These findings reveal that enhanced charge and spin fluctuations with optimal doping is able to produce robust d-wave superconductivity in doped Mott insulators, providing a foundation for connecting theories of superconductivity to models of strongly correlated systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2003

Flutter and Divergence Analysis using the Generalized Aeroelastic Analysis Method

The Generalized Aeroelastic Analysis Method (GAAM) is applied to the analysis of three well-studied checkcases: restrained and unrestrained airfoil models, and a wing model. An eigenvalue iteration procedure is used for converging upon roots of the complex stability matrix. For the airfoil models, exact root loci are given which clearly illustrate the nature of the flutter and divergence instabilities. The singularities involved are enumerated, including an additional pole at the origin for the unrestrained airfoil case and the emergence of an additional pole on the positive real axis at the divergence speed for the restrained airfoil case. Inconsistencies and differences among published aeroelastic root loci and the new, exact results are discussed and resolved. The generalization of a Doublet Lattice Method computer code is described and the code is applied to the calculation of root loci for the wing model for incompressible and for subsonic flow conditions. The error introduced in the reduction of the singular integral equation underlying the unsteady lifting surface theory to a linear algebraic equation is discussed. Acknowledging this inherent error, the solutions of the algebraic equation by GAAM are termed 'exact.' The singularities of the problem are discussed and exponential series approximations used in the evaluation of the kernel function shown to introduce a dense collection of poles and zeroes on the negative real axis. Again, inconsistencies and differences among published aeroelastic root loci and the new 'exact' results are discussed and resolved. In all cases, aeroelastic flutter and divergence speeds and frequencies are in good agreement with published results. The GAAM solution procedure allows complete control over Mach number, velocity, density, and complex frequency. Thus all points on the computed root loci can be matched-point, consistent solutions without recourse to complex mode tracking logic or dataset interpolation, as in the k and p-k solution methods.

Edwards, John W.↗
DOE OSTI2025

Neutrino properties from muonium-antimuonium mixing

The nature of neutrino mass — whether neutrinos are Dirac or Majorana particles — remains one of the central open questions in particle physics. While observation of neutrinoless double beta decay would confirm the Majorana case, the absence of a signal offers no definitive insight. In light of this, we investigate muonium-antimuonium mixing — proposed for further study at the MACE experiment — as an alternative probe of neutrino properties. We compute the mixing amplitude in the Standard Model minimally extended to include massive Dirac or Majorana neutrinos, and correct previous calculations by properly treating the relevant infrared scales. As the GIM mechanism strongly suppresses the Dirac contribution, we explore whether relaxing unitarity of the PMNS matrix can enhance the mixing without obscuring neutrino properties. Surprisingly, the answer to this question is negative. We also examine the pseudo-Dirac case — predominantly Dirac neutrinos with small Majorana masses — and find that this scenario can significantly enhance the mixing compared to the pure Dirac case, especially for normal mass ordering.

Neutrino interactions↗
DOE OSTI2024

Systematic improvement of x -dependent unpolarized nucleon generalized parton distributions in lattice-QCD calculation

We present a first study of the effects of renormalization-group resummation (RGR) and leading-renormalon resummation (LRR) on the systematic errors of the unpolarized isovector nucleon generalized parton distribution in the framework of large-momentum effective theory. This work is done using lattice gauge ensembles generated by the MILC Collaboration, consisting of 2 + 1 + 1 flavors of highly improved staggered quarks with a physical pion mass at lattice spacing a ≈ 0.09 fm and a box width L ≈ 5.76 fm . We present results for the nucleon H and E generalized parton distributions (GPDs) with average boost momentum P z ≈ 2 GeV at momentum transfers Q 2 = [ 0 , 0.97 ] GeV 2 at skewness ξ = 0 as well as Q 2 ∈ 0.23 GeV 2 at ξ = 0.1 , renormalized in the modified minimal subtraction ( MS ¯ ) scheme at scale μ = 2.0 GeV , with two- and one-loop matching, respectively. We demonstrate that the simultaneous application of RGR and LRR significantly reduces the systematic errors in renormalized matrix elements and distributions for both the zero and nonzero skewness GPDs, and that it is necessary to include both RGR and LRR at higher orders in the matching and renormalization processes. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
NASA NTRS1983

Effects of cobalt in nickel-base superalloys

The role of cobalt in a representative wrought nickel-base superalloy was determined. The results show cobalt affecting the solubility of elements in the gamma matrix, resulting in enhanced gamma' volume fraction, in the stabilization of MC-type carbides, and in the stabilization of sigma phase. In the particular alloy studied, these microstructural and microchemistry changes are insufficient in extent to impact on tensile strength, yield strength, and in the ductilities. Depending on the heat treatment, creep and stress rupture resistance can be cobalt sensitive. In the coarse grain, fully solutioned and aged condition, all of the alloy's 17% cobalt can be replaced by nickel without deleteriously affecting this resistance. In the fine grain, partially solutioned and aged condition, this resistance is deleteriously affected only when one-half or more of the initial cobalt content is removed. The structure and property results are discussed with respect to existing theories and with respect to other recent and earlier findings on the impact of cobalt, if any, on the performance of nickel-base superalloys.

Tien, J. K.↗
DOE OSTI2020

Machine learning a molecular Hamiltonian for predicting electron dynamics

We develop a computational method to learn a molecular Hamiltonian matrix from matrix-valued time series of the electron density. As we demonstrate for three small molecules, the resulting Hamiltonians can be used for electron density evolution, producing highly accurate results even when propagating 1,000 time steps beyond the training data. As a more rigorous test, we use the learned Hamiltonians to simulate electron dynamics in the presence of an applied electric field, extrapolating to a problem that is beyond the field-free training data. We find that the resulting electron dynamics predicted by our learned Hamiltonian are in close quantitative agreement with the ground truth. Our method relies on combining a reduced-dimensional, linear statistical model of the Hamiltonian with a time-discretization of the quantum Liouville equation within time-dependent Hartree Fock theory. In conclusion, we train the model using a least-squares solver, avoiding numerous, CPU-intensive optimization steps. For both field-free and field-on problems, we quantify training and propagation errors, highlighting areas for future development.

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

Microcarbonation of Naphthalene: An Experimental and Computational Study of Photoionization in Naphthalene-Carbon Dioxide Clusters

The photoionization of naphthalene (N)-carbon dioxide (CO 2 ) clusters was studied using tunable vacuum ultraviolet (VUV) radiation from a synchrotron in the photon range of 8.0 to 13.7 eV, in combination with time-of-flight mass spectrometry. Clusters of monomer, dimer, and trimer naphthalene with CO 2 (N­(CO 2 ) 0–6 , N 2 ­(CO 2 ) 0–3 , N3) were observed. The lowest-energy conformers were obtained via a conformer search, followed by geometry optimizations at the ωB97X-V2/aug-cc-pVTZ (monomer) and ωB97X-V2/aug-cc-pVDZ (dimer) levels of theory. Carbon dioxide was found to preferentially cluster on top of the naphthalene molecule (in an out-of-plane configuration). From the mass spectra, photoionization intensity curves (PICs) were constructed, and appearance energies (AEs) were determined. No substantial trend in AE was observed with increasing size of the naphthalene-carbon dioxide clusters; rather, AE oscillations around the value for pure naphthalene were observed. These AE oscillations are also observed in a recently studied naphthalene-water cluster system, though in this system a slight downward trend in AE for pure naphthalene clusters (N1/2/3/4) was observed. The differences between the two systems are attributed to differing interaction strengths. In conclusion, understanding these differences may aid in determining how photoprocessing can proceed differently depending on the dominant matrix component in interstellar ices.

Wannenmacher, Anna [Lawrence Berkeley National Lab↗
DOE OSTI2021

Convergence of Eigenvector Continuation

Eigenvector continuation is a computational method that finds the extremal eigenvalues and eigenvectors of a Hamiltonian matrix with one or more control parameters. It does this by projection onto a subspace of eigenvectors corresponding to selected training values of the control parameters. The method has proven to be very efficient and accurate for interpolating and extrapolating eigenvectors. However, almost nothing is known about how the method converges, and its rapid convergence properties have remained mysterious. In this Letter, we present the first study of the convergence of eigenvector continuation. In order to perform the mathematical analysis, we introduce a new variant of eigenvector continuation that we call vector continuation. We first prove that eigenvector continuation and vector continuation have identical convergence properties and then analyze the convergence of vector continuation. Our analysis shows that, in general, eigenvector continuation converges more rapidly than perturbation theory. The faster convergence is achieved by eliminating a phenomenon that we call differential folding, the interference between nonorthogonal vectors appearing at different orders in perturbation theory. From our analysis we can predict how eigenvector continuation converges both inside and outside the radius of convergence of perturbation theory. Further, while eigenvector continuation is a nonperturbative method, we show that its rate of convergence can be deduced from power series expansions of the eigenvectors. Our results also yield new insights into the nature of divergences in perturbation theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1991

Unified aeroacoustics analysis for high speed turboprop aerodynamics and noise. Volume 4: Computer user's manual for UAAP turboprop aeroacoustic code

The Unified AeroAcoustic Program (UAAP) code calculates the airloads on a single rotation prop-fan, or propeller, and couples these airloads with an acoustic radiation theory, to provide estimates of near-field or far-field noise levels. The steady airloads can also be used to calculate the nonuniform velocity components in the propeller wake. The airloads are calculated using a three dimensional compressible panel method which considers the effects of thin, cambered, multiple blades which may be highly swept. These airloads may be either steady or unsteady. The acoustic model uses the blade thickness distribution and the steady or unsteady aerodynamic loads to calculate the acoustic radiation. The users manual for the UAAP code is divided into five sections: general code description; input description; output description; system description; and error codes. The user must have access to IMSL10 libraries (MATH and SFUN) for numerous calls made for Bessel functions and matrix inversion. For plotted output users must modify the dummy calls to plotting routines included in the code to system-specific calls appropriate to the user's installation.

Menthe, R. W.↗
DOE OSTI2022

Application of the bRapid fission matrix burnup methodology to the JSI TRIGA Mark II research reactor

The novel burnup methodology, bRAPID, of the RAPID code system has been applied to the TRIGA Mark II research reactor at the 'Jozef Stefan' Institute. This paper presents the first step in our goal for on-line determination of isotope inventory and burnup, that the RAPID code system with its Fission Matrix (FM) method enables. Burnup analysis of changes in FM coefficients due to TRIGA fuel burnup at different irradiation times and specific power distribution in the reactor core is presented. Our results demonstrate that primary contribution for changes of FM coefficients is the burnup of the destination s, while the spectrum differences, that arise due to fuel's core location, have negligible effect on FM coefficients. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

First-Principles Quantum Transport Modeling of Spin-Transfer and Spin-Orbit Torques in Magnetic Multilayers

A unified approach for computing (i) spin-transfer torque in magnetic trilayers like spin valve and magnetic tunnel junction, where injected charge current flows perpendicularly to interfaces, and (ii) spin-orbit torque in magnetic bilayers of the type ferromagnet/spin-orbit-coupled material, where injected charge current flows parallel to the interface, was reviewed. The experimentally explored and technologically relevant spin-orbit-coupled materials include 5d heavy metals, topological insulators, Weyl semimetals, and transition metal dichalcogenides. This approach requires to construct the torque operator for a given Hamiltonian of the device and the steady-state nonequilibrium density matrix, where the latter is expressed in terms of the nonequilibrium Green’s functions and split into three contributions. Tracing these contributions with the torque operator automatically yields field-like and damping-like components of spin-transfer torque or spin-orbit torque vector, which is particularly advantageous for spin-orbit torque where the direction of these components depends on the unknown-in-advance orientation of the current-driven nonequilibrium spin density in the presence of spin-orbit coupling. Illustrative examples are provided by computing spin-transfer torque in a one-dimensional toy model of a magnetic tunnel junction and realistic Co/Cu/Co spin valve, both of which are described by first-principles Hamiltonians obtained from noncollinear density functional theory calculations, as well as by computing spin-orbit torque in a ferromagnetic layer described by a tight-binding Hamiltonian which includes spin-orbit proximity effect within ferromagnetic monolayers assumed to be generated by the adjacent monolayer transition metal dichalcogenide. In addition, it is shown here how spin-orbit proximity effect, quantified by computing (via first-principles retarded Green’s function) spectral functions and spin textures on monolayers of realistic ferromagnetic material like Co in contact with heavy metal or monolayer transition metal dichalcogenide, can be tailored to enhance the magnitude of spin-orbit torque. Furthermore, errors made in the calculation of spin-transfer torque are quantified when using Hamiltonian from collinear density functional theory, with rigidly rotated magnetic moments to create noncollinear magnetization configurations, instead of proper (but computationally more expensive) self-consistent Hamiltonian obtained from noncollinear density functional theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗