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.

SourceResearch connectionFrom yearThrough year

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

At least 1,333 records · Page 74

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Verification of Coupled Hydraulic Fracturing Simulators Using Laboratory-Scale Experiments

In this work, we aim to verify the predictions of the numerical simulators, which are used for designing field-scale hydraulic stimulation experiments. Although a strong theoretical understanding of this process has been gained over the past few decades, numerical predictions of fracture propagation in low-permeability rocks still remains a challenge. Against this background, we performed controlled laboratory-scale hydraulic fracturing experiments in granite samples, which not only provides high-quality experimental data but also a well-characterized experimental set-up. Using the experimental pressure responses and the final fracture sizes as benchmark, we compared the numerical predictions of two coupled hydraulic fracturing simulators—CSMP and GEOS. Both the simulators reproduced the experimental pressure behavior by implementing the physics of Linear Elastic Fracture Mechanics (LEFM) and lubrication theory within a reasonable degree of accuracy. The simulation results indicate that even in the very low-porosity (1–2 %) and low-permeability (10 -18 m 2 - 10 -19 m 2 ) crystalline rocks, which are usually the target of EGS, fluid-loss into the matrix and unsaturated flow impacts the formation breakdown pressure and the post-breakdown pressure trends. Therefore, underestimation of such parameters in numerical modeling can lead to significant underestimation of breakdown pressure. The simulation results also indicate the importance of implementing wellbore solvers for considering the effect of system compressibility and pressure drop due to friction in the injection line. The varying injection rate as a result of decompression at the instant of fracture initiation affects the fracture size, while the entry friction at the connection between the well and the initial notch may cause an increase in the measured breakdown pressure.

15 GEOTHERMAL ENERGY↗

DOE OSTI2024

Signatures of the attractive interaction in spin spectra of one-dimensional cuprate chains

Identifying the minimal model for cuprates is crucial for explaining the high- T c pairing mechanism. Recent photoemission experiments have suggested a significant near-neighbor attractive interaction V in cuprate chains, favoring pairing instability. To determine its strength, we systematically investigate the dynamical spin structure factors S ( q , ω ) using the density matrix renormalization group. Our analysis quantitatively reveals a notable softening in the two-spinon continuum, particularly evident in the intense spectrum at large momentum. This softening is primarily driven by the renormalization of the superexchange interaction, as determined by a comparison with the slave-boson theory. We also demonstrate the feasibility of detecting this spectral shift in thin-film samples using resonant inelastic x-ray scattering. Therefore, this provides a distinctive fingerprint for the attractive interaction, motivating future experiments to unveil essential ingredients in cuprates. Published by the American Physical Society 2024

Shen, Zecheng↗

NASA NTRS1978

Finite element concepts in computational aerodynamics

Finite element theory was employed to establish an implicit numerical solution algorithm for the time averaged unsteady Navier-Stokes equations. Both the multidimensional and a time-split form of the algorithm were considered, the latter of particular interest for problem specification on a regular mesh. A Newton matrix iteration procedure is outlined for solving the resultant nonlinear algebraic equation systems. Multidimensional discretization procedures are discussed with emphasis on automated generation of specific nonuniform solution grids and accounting of curved surfaces. The time-split algorithm was evaluated with regards to accuracy and convergence properties for hyperbolic equations on rectangular coordinates. An overall assessment of the viability of the finite element concept for computational aerodynamics is made.

Baker, A. J.↗

NASA NTRS1998

Predictive Models for Semiconductor Device Design and Processing

The device feature size continues to be on a downward trend with a simultaneous upward trend in wafer size to 300 mm. Predictive models are needed more than ever before for this reason. At NASA Ames, a Device and Process Modeling effort has been initiated recently with a view to address these issues. Our activities cover sub-micron device physics, process and equipment modeling, computational chemistry and material science. This talk would outline these efforts and emphasize the interaction among various components. The device physics component is largely based on integrating quantum effects into device simulators. We have two parallel efforts, one based on a quantum mechanics approach and the second, a semiclassical hydrodynamics approach with quantum correction terms. Under the first approach, three different quantum simulators are being developed and compared: a nonequlibrium Green's function (NEGF) approach, Wigner function approach, and a density matrix approach. In this talk, results using various codes will be presented. Our process modeling work focuses primarily on epitaxy and etching using first-principles models coupling reactor level and wafer level features. For the latter, we are using a novel approach based on Level Set theory. Sample results from this effort will also be presented.

Meyyappan, Meyya↗

NASA NTRS2018

The Effect of Thin Film Adhesives on Mode I Interlaminar Fracture Toughness in Carbon Fiber Composites with Shape Memory Alloy Inserts

Shape Memory Alloy (SMA) was placed within Polymer Matrix Composite (PMC) panels alongside film adhesives to examine bonding. Double cantilever beam (DCB) testing was performed using ASTM D5528. C-scanning was performed before testing, modal acoustic emissions (MAE) were monitored during testing, and microscopy performed post-test. Data was analyzed using modified beam theory (MBT), compliance calibration (CC) and modified compliance calibration (MCC) methods. Fracture toughness for control specimens was higher than previously reported due to fiber-bridging. Specimens with SMAs and adhesives stabilized crack propagation. Results revealed SMA-bridging; a phenomenon mimicking fiber-bridging which increased the load and fracture toughness of SMA specimens.

Derek J Quade↗

DOE OSTI2020

Eigenmode analysis of the sheared-flow Z-pinch

Experiments have demonstrated that a Z-pinch can persist for thousands of times longer than the growth time of global magnetohydrodynamic (MHD) instabilities such as the m=0 sausage and m=1 kink modes. These modes have growth times on the order of ta=a/vi, where vi is the ion thermal speed and a is the pinch radius. Axial flows with duz/dr ≲ vi/a have been measured during the stable period, and the commonly accepted theory is that this amount of shear is sufficient to stabilize these modes as predicted by numerical studies using the ideal MHD equations. However, these studies only consider specific equilibrium profiles that typically have a modest magnitude for the logarithmic pressure gradient, qP≡d ln P/d ln r, and may not represent experimental conditions. Linear stability of the sheared-flow Z-pinch is studied here via a direct eigen-decomposition of the matrix operator obtained from the linear ideal MHD equations. Several equilibrium profiles with a large variation of qP are examined. Considering a practical range of k, 1/3 ≲ ka ≲ 10, it is shown that the shear required to stabilize m=0 modes can be expressed as duz/dr≥Cγ0/(ka)α. Here, γ0=γ0(ka) is the profile-specific growth rate in the absence of shear, which scales approximately with |qP|. Both C and α are profile-specific constants, but C is order unity and α≈1. It is further demonstrated that even a large value of shear, duz/dr=3vi/a, is not sufficient to provide linear stabilization of the m=1 kink mode for all profiles considered. This result is in contrast to the currently accepted theory predicting stabilization at much lower shear, duz/dr=0.1vi/a, and suggests that the experimentally observed stability cannot be explained within the linear ideal-MHD model.

Angus, J. R. (ORCID:0000000314740002)↗

NASA NTRS2016

The Relaxation Matrix for Symmetric Tops with Inversion Symmetry. I. Effects of Line Coupling on Self-Broadened v (sub 1) and Pure Rotational Bands of NH3

The Robert-Bonamy formalism has been commonly used to calculate half-widths and shifts of spectral lines for decades. This formalism is based on several approximations. Among them, two have not been fully addressed: the isolated line approximation and the neglect of coupling between the translational and internal motions. Recently, we have shown that the isolated line approximation is not necessary in developing semi-classical line shape theories. Based on this progress, we have been able to develop a new formalism that enables not only to reduce uncertainties on calculated half-widths and shifts, but also to model line mixing effects on spectra starting from the knowledge of the intermolecular potential. In our previous studies, the new formalism had been applied to linear and asymmetric-top molecules. In the present study, the method has been extended to symmetric-top molecules with inversion symmetry. As expected, the inversion splitting induces a complete failure of the isolated line approximation. We have calculated the complex relaxation matrices of selfbroadened NH3. The half-widths and shifts in the ν1 and the pure rotational bands are reported in the present paper. When compared with measurements, the calculated half-widths match the experimental data very well, since the inapplicable isolated line approximation has been removed. With respect to the shifts, only qualitative results are obtained and discussed. Calculated off-diagonal elements of the relaxation matrix and a comparison with the observed line mixing effects are reported in the companion paper (Paper II).

Trajectory models↗

DOE OSTI2022

Easy-plane anisotropic-exchange magnets on a honeycomb lattice: Quantum effects and dealing with them

We provide analytical and numerical insights into the phase diagram and other properties of the extended Kitaev-Heisenberg model on the honeycomb lattice in the easy-plane limit, in which interactions are only between spin components that belong to the plane of magnetic ions. This parameter subspace allows for a much-needed systematic quantitative investigation of spin excitations in the ordered phases and of their generic features. Specifically, we demonstrate that in this limit one can consistently take into account magnon interactions in both zero-field zigzag and field-polarized phases. For the nominally polarized phase, we propose a regularization of the unphysical divergences that occur at the critical field and are plaguing the 1/S approximation in this class of models. For the explored parameter subspace, all symmetry-allowed terms of the standard parametrization of the extended Kitaev-Heisenberg model, such as K, J, and Γ, are significant, making the offered consideration relevant to a much wider parameter space. Furthermore, the dynamical structure factor near the paramagnetic critical point illustrates this relevance by showing features that are reminiscent of the ones observed in α–RuCl 3 , underscoring that they are not unique and should be common to a wide range of parameters of the model and, by extension, to other materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Isospin 0 and 2 two-pion scattering at physical pion mass using all-to-all propagators with periodic boundary conditions in lattice QCD

A study of two-pion scattering for the isospin channels, 𝐼 = 0 and 𝐼 = 2, using lattice QCD is presented. Möbius domain-wall fermions, on top of the Iwasaki-DSDR gauge action for gluons with periodic boundary conditions, are used for the lattice computations, which are carried out on two ensembles of gauge field configurations generated by the RBC and UKQCD Collaborations with physical masses, inverse lattice spacings of 1.023 and 1.378 GeV, and spatial extents of 𝐿 = 4.63 and 4.58 fm, respectively. The all-to-all propagator method is employed to compute a matrix of correlation functions of two-pion operators. The generalized eigenvalue problem (GEVP) is solved for a matrix of correlation functions to extract phase shifts with multiple states—two pions with a nonzero relative momentum, as well as two pions at rest. Our results for phase shifts for both the 𝐼 = 0 and 𝐼 = 2 channels are consistent with the Roy equation and chiral perturbation theory, though at this preliminary stage our errors for 𝐼 = 0 are large. An important outcome of this work is that we are successful in extracting two-pion excited states, which are useful for studying 𝐾 → 𝜋⁢𝜋 decay, on physical-mass ensembles using the GEVP.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Rapid Design of Cost-Effective Refractory High Entropy Alloys Strengthened by Precipitation

Funded by the ARPA-E ULTIMATE program, this project aims to develop cost-effective refractory high entropy alloys (RHEAs) for gas turbine blades applications at 1300 degree Celsius and above. The alloy design focuses on carbides precipitation strengthening in a strong and ductile solid solution BCC matrix. The goal is to achieve balanced mechanical properties including room temperature ductility & fracture toughness and high temperature strength and creep resistance, while maintaining comparable oxidation resistance and densities to Ni-base superalloys. To accelerate alloy design, multiscale modeling including CALPHAD and first-principles density functional theory (DFT) calculations and machine learning are performed. Specifically, high throughput phase diagram calculations and screening are carried out using CALPHAD; DFT calculations are performed to predict intrinsic ductility, grain boundary strengthening potential, coefficient of thermal expansion (CTE) and temperature-dependent elastic constants. Temperature dependent yield stress, creep strain and CTE are also predicted using machine learning. Down-selected alloys are synthesized in small buttons of about 250g using arc melting for rapid evaluation on microstructure and mechanical properties before employing plasma arc melting and additive manufacturing for producing large ingots. Preliminary computational and experimental results of this project will be presented and future direction in RHEA development will be discussed.

high temperature materials↗

NASA NTRS1985

Viscoelastic behavior and life-time predictions

Fiber reinforced plastics were considered for many structural applications in automotive, aerospace and other industries. A major concern was and remains the failure modes associated with the polymer matrix which serves to bind the fibers together and transfer the load through connections, from fiber to fiber and ply to ply. An accelerated characterization procedure for prediction of delayed failures was developed. This method utilizes time-temperature-stress-moisture superposition principles in conjunction with laminated plate theory. Because failures are inherently nonlinear, the testing and analytic modeling for both moduli and strength is based upon nonlinear viscoelastic concepts.

Dillard, D. A.↗

NASA NTRS1994

Coupled 2-dimensional cascade theory for noise and unsteady aerodynamics of blade row interaction in turbofans. Volume 1: Theory development and parametric studies

Typical analytical models for interaction between rotor and stator in a turbofan analyze the effect of wakes from the rotor impinging on the stator, producing unsteady loading, and thereby generating noise. Reflection/transmission characteristics of the rotor are sometimes added in a separate calculation. In those models, there is a one-to-one relationship between wake harmonics and noise harmonics; that is, the BPF (blade passing frequency) wake harmonic causes only the BPF noise harmonic, etc. This report presents a more complete model in which flow tangency boundary conditions are satisfied on two cascades in relative motion for several harmonics simultaneously. By an extension of S.N. Smith's code for two dimensional flat plate cascades, the noise generation/frequency scattering/blade row reflection problem is solved in a single matrix inversion. It is found that the BPF harmonic excitation of the stator scatters considerable energy in the higher BPF harmonics due to relative motion between the blade rows. Furthermore, when swirl between the rotor and stator is modeled, a 'mode trapping' effect occurs which explains observations on fans operating at rotational speeds below BFP cuton: the BPF mode amplifies between blade rows by multiple reflections but cannot escape to the inlet and exit ducts. However, energy scattered into higher harmonics does propagate and dominates the spectrum at two and three times BPF. This report presents the complete derivation of the theory, comparison with a previous (more limited) coupled rotor/stator interaction theory due to Kaji and Okazaki, exploration of the mode trapping phenomenon, and parametric studies showing the effects of vane/blade ratio and rotor/stator interaction. For generality, the analysis applies to stages where the rotor is either upstream or downstream of the stator and to counter rotation stages. The theory has been coded in a FORTRAN program called CUP2D, documented in Volume 2 of this report. It is concluded that the new features of this analysis - unsteady coupling, frequency scattering, and flow turning between rotor and stator - have a profound effect on noise generation caused by rotor/stator interaction. Treating rotors and stators as isolated cascades is not adequate for noise analysis and prediction.

Hanson, Donald B.↗

DOE OSTI2026

Proton isovector helicity PDF at NNLO and the twist-3 moment $\tilde{d}$ 2 from lattice QCD at physical quark masses

We present a lattice quantum chromodynamics calculation of the 𝑥-dependent isovector quark helicity parton distribution function (PDF) of the proton in the large momentum effective theory (LaMET) framework. Through operator product expansion (OPE) we also extract the $\tilde{d}$ 2 moment of the twist-3 PDF 𝑔 𝑇 ⁡(𝑥) for the first time in the $\overline{MS}$ scheme, which is proportional to the average color Lorentz force experienced by the quark in the proton. This calculation is performed on a lattice of spacing 𝑎 =0.076 fm at physical quark masses. The quasi-PDF matrix elements are measured in proton states boosted to momenta 𝑃 𝑧 ={0,0.25,1.02,1.53} GeV. We first extract the lowest few helicity PDF moments from the renormalization-group (RG) invariant ratios of the matrix elements with OPE. Combined with the matrix elements relevant for 𝑔 𝑇 ⁡(𝑥), we obtain $\tilde{d}$$^{u-d}_2$⁡(2 GeV) =0.0024⁢(46) at next-to-leading order in $\overline{MS}$. Then, the helicity quasi-PDF matrix elements are renormalized in the hybrid scheme with linear renormalon resummation and Fourier transformed to the 𝑥-space after an asymptotic extrapolation. The quasi-PDF is perturbatively matched to the $\overline{MS}$ PDF with RG and threshold resummations at next-to-leading power and next-to-next-to-leading logarithmic accuracies. After resummations, we determine the PDF in the region 𝑥 ∈[0.25,0.75]. The end-point regions are then parametrized, combined with the LaMET prediction at moderate 𝑥, and fitted to the short-distance matrix elements in coordinate space.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Spin–orbit effects on the electronic and optical properties of lead iodide

Lead iodide (PbI2) has gained much interest due to its direct electronic gap in the visible range and layered crystal structure. It has thereby been considered as a promising material for applications in atomically thin optoelectronic devices. In this work, we present a detailed investigation of the effect of spin–orbit coupling (SOC) that arises from the presence of heavy atoms on the electronic and optical properties of PbI2 using first-principles calculations based on density-functional theory and many-body perturbation theory. We find that SOC not only alters the bandgap but also induces the mixing of orbital characters, resulting in a significant change in the overall band structure and charge carrier effective masses. Moreover, the band orbital mixing caused by SOC results in the dramatic change in optical transition matrix elements and, correspondingly, the absorption spectrum. Our experimentally measured absorption spectra validate the calculation results and demonstrate the importance of SOC in the optical processes of PbI2. Our findings provide insights that are important for the potential use of PbI2 as a material platform for visible optoelectronic devices.

Physics↗

DOE OSTI2023

Subspace recursive Fermi-operator expansion strategies for large-scale DFT eigenvalue problems on HPC architectures

Quantum mechanical calculations for material modeling using Kohn–Sham density functional theory (DFT) involve the solution of a nonlinear eigenvalue problem for N smallest eigenvector-eigenvalue pairs, with N proportional to the number of electrons in the material system. Here, these calculations are computationally demanding and have asymptotic cubic scaling complexity with the number of electrons. Large-scale matrix eigenvalue problems arising from the discretization of the Kohn–Sham DFT equations employing a systematically convergent basis traditionally rely on iterative orthogonal projection methods, which are shown to be computationally efficient and scalable on massively parallel computing architectures. However, as the size of the material system increases, these methods are known to incur dominant computational costs through the Rayleigh–Ritz projection step of the discretized Kohn–Sham Hamiltonian matrix and the subsequent subspace diagonalization of the projected matrix. This work explores the potential of polynomial expansion approaches based on recursive Fermi-operator expansion as an alternative to the subspace diagonalization of the projected Hamiltonian matrix to reduce the computational cost. Subsequently, we perform a detailed comparison of various recursive polynomial expansion approaches to the traditional approach of explicit diagonalization on both multi-node central processing unit and graphics processing unit architectures and assess their relative performance in terms of accuracy, computational efficiency, scaling behavior, and energy efficiency.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Low-energy physics of isotropic spin-1 chains in the critical and Haldane phases

Using a matrix product state algorithm with infinite boundary conditions, we compute high-resolution dynamic spin and quadrupolar structure factors in the thermodynamic limit to explore the low-energy excitations of isotropic bilinear-biquadratic spin-1 chains. Haldane mapped the spin-1 Heisenberg antiferromagnet to a continuum field theory, the nonlinear sigma model (NLσM). In this work, we find that the NLσM fails to capture the influence of the biquadratic term and provides only an unsatisfactory description of the Haldane phase physics. But several features in the Haldane phase can be explained by noninteracting multimagnon states. The physics at the Uimin-Lai-Sutherland point is characterized by multisoliton continua. Moving into the extended critical phase, we find that these excitation continua contract, which we explain using a field-theoretic description. New excitations emerge at higher energies and, in the vicinity of the purely biquadratic point, they show simple cosine dispersions. Using block fidelities, we identify them as elementary one-particle excitations and relate them to the integrable Temperley-Lieb chain.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Three-dimensional dispersion in the type-II Dirac semimetals PtTe 2 and PdTe 2 revealed through circular dichroism in angle-resolved photoemission spectroscopy

PtTe 2 and PdTe 2 are among the first transition metal dichalcogenides that were predicted to host type-II Dirac fermions, exotic particles prohibited in free space. These materials are layered and air stable, which makes them top candidates for technological applications that take advantage of their anisotropic magnetotransport properties. Here, in this work, we provide a detailed characterization of the electronic structure of PtTe 2 and PdTe 2 using angle-resolved photoemission spectroscopy (ARPES) and density functional theory calculations, offering an alternative interpretation for one of the Dirac-like dispersions in these materials. Through the use of circularly polarized light, we report a different behavior of such dispersion in PdTe 2 compared to PtTe 2 , that we relate to a symmetry analysis of the dipole matrix element. Such analysis reveals a link between the observed circular dichroism and the different momentum-dependent terms in the dispersion of these two compounds, despite their close similarity in crystal structure. Additionally, our data show a clear difference in the circular dichroic signal for the type-II Dirac cones characteristic of these materials, compared to their topologically protected surface states. Our paper provides a useful reference for the ARPES characterization of other transition metal dichalcogenides with topological properties and illustrates the use of circular dichroism as a guide to identify the topological character and attributes of two otherwise equivalent band dispersions.

angle-resolved photoemission spectroscopy↗

DOE OSTI2022

Robust ab initio prediction of nuclear electric quadrupole observables by scaling to the charge radius

Meaningful predictions for electric quadrupole (E2) observables from ab initio nuclear theory are necessary, if the ab initio description of collective correlations is to be confronted with experiment, as well as to provide predictive power for unknown E2 observables. However, converged results for E2 observables are notoriously challenging to obtain in ab initio no-core configuration interaction (NCCI) approaches. Matrix elements of the E2 operator are sensitive to the large-distance tails of the nuclear wave function, which converge slowly in an oscillator basis expansion. Similar convergence challenges beset ab initio prediction of the nuclear charge radius. We demonstrate that the convergence patterns of the E2 and radius observables are strongly correlated, and that meaningful predictions for the absolute scale of E2 observables may be made by calibrating to the experimentally-known ground-state charge radius. Here, we illustrate by providing robust ab initio predictions for several E2 transition strengths and quadrupole moments in p-shell nuclei, in cases where experimental results are available for comparison.

6 ≤ A ≤ 19↗