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,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Scalable semidefinite programming approach to variational embedding for quantum many-body problems

In quantum embedding theories, a quantum many-body system is divided into localized clusters of sites which are treated with an accurate ‘high-level’ theory and glued together self-consistently by a less accurate ‘low-level’ theory at the global scale. The recently introduced variational embedding approach for quantum many-body problems combines the insights of semidefinite relaxation and quantum embedding theory to provide a lower bound on the ground-state energy that improves as the cluster size is increased. The variational embedding method is formulated as a semidefinite program (SDP), which can suffer from poor computational scaling when treated with black-box solvers. Here, we exploit the interpretation of this SDP as an embedding method to develop an algorithm which alternates parallelizable local updates of the high-level quantities with updates that enforce the low-level global constraints. Moreover, we show how translation invariance in lattice systems can be exploited to reduce the complexity of projecting a key matrix to the positive semidefinite cone.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1983

Energetics and optimum motion of oscillating lifting surfaces

Low-frequency, unsteady, lifting-line theory is used to characterize the energetics and optimum motion of an unswept rigid wing oscillating harmonically in an inviscid, incompressible flow. The energetics calculations account for the leading edge suction force, the power absorbed in the wing oscillations, and the energy loss rate produced by vortex shedding. Optimization is achieved by minimizing the average energy loss rate in relation to a given thrust, and a unique solution is found in the three dimensional case for low, reduced frequencies. The two-dimensional solution is nonunique, a condition which is examined in terms of the normal modes of the energy loss rate matrix. An invisible mode with a hydrodynamic efficiency of 100 pct is obtained in the two-dimensional case, causing the nonuniqueness of the solution by yielding no fixed positive thrust through perfect unsteady feathering.

Ahmadi, A. R.↗
NASA NTRS1992

FONO - A difficult case for theory and experiment

High levels of ab initio theory are used to investigate the equilibrium structures, vibrational spectra, and relative energetics of FNO2, cis-FONO, and trans-FONO isomers. FNO2 is determined to be 36.9 +/- 2.5 kcal/mol (0 K) more stable than cis-FONO, which is more stable than trans-FONO by 2.5 +/- 1.0 kcal/mol (0 K). The molecular structure of cis-FONO is shown to have typical F-O and central O-N single bond distances, indicating that cis-FONO has a structure analogous to HONO. The computed vibrational spectrum of cis-FONO is shown to be consistent with IR matrix isolation experiments, including isotopic shifts. The experimentally deduced structure of FONO is shown to be incorrect because one of the vibrational bands included in the normal coordinate analysis is either a combination band or an overtone.

Lee, Timothy J.↗
DOE OSTI2020

Improved band gaps and structural properties from Wannier–Fermi–Löwdin self-interaction corrections for periodic systems

The accurate prediction of band gaps and structural properties in periodic systems continues to be one of the central goals of electronic structure theory. However, band gaps obtained from popular exchange–correlation (XC) functionals (such as LDA and PBE) are severely underestimated partly due to the spurious self-interaction error (SIE) inherent to these functionals. In this work, we present a new formulation and implementation of Wannier function-derived Fermi–Löwdin (WFL) orbitals for correcting the SIE in periodic systems. Since our approach utilizes a variational minimization of the self-interaction energy with respect to the Wannier charge centers (WCC), it is computationally more efficient than the HSE hybrid functional and other self-interaction corrections that require a large number of transformation matrix elements. Here, calculations on several (17 in total) prototypical molecular solids, semiconductors, and wide-bandgap materials show that our WFL self-interaction correction approach gives better band gaps and bulk moduli compared to semilocal functionals, largely due to the partial removal of self-interaction errors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Developments in Performance and Portability for MadGraph5_aMC@NLO

Event generators simulate particle interactions using Monte Carlo techniques, providing the primary connection between experiment and theory in experimental high energy physics. These software packages, which are the first step in the simulation worflow of collider experiments, represent approximately 5 to 20% of the annual WLCG usage for the ATLAS and CMS experiments. With computing architectures becoming more heterogeneous, it is important to ensure that these key software frameworks can be run on future systems, large and small. In this contribution, recent progress on porting and speeding up the Madgraph5_aMC@NLO event generator on hybrid architectures, i.e. CPU with GPU accelerators, is discussed. The main focus of this work has been in the calculation of scattering amplitudes and "matrix elements", which is the computational bottleneck of an event generation application. For physics processes limited to QCD leading order, the code generation toolkit has been expanded to produce matrix element calculations using C++ vector instructions on CPUs and using CUDA for NVidia GPUs, as well as using Alpaka, Kokkos and SYCL for multiple CPU and GPU architectures. Performance is reported in terms of matrix element calculations per time on NVidia, Intel, and AMD devices. The status and outlook for the integration of this work into a production release usable by the LHC experiments, with the same functionalities and very similar user interfaces as the current Fortran version, is also described.

Valassi, Andrea↗
NASA NTRS2013

Optical Alignment of the Global Precipitation Measurement (GPM) Star Trackers

The optical alignment of the star trackers on the Global Precipitation Measurement (GPM) core spacecraft at NASA Goddard Space Flight Center (GSFC) was challenging due to the layout and structural design of the GPM Lower Bus Structure (LBS) in which the star trackers are mounted as well as the presence of the star tracker shades that blocked line-of-sight to the primary star tracker optical references. The initial solution was to negotiate minor changes in the original LBS design to allow for the installation of a removable item of ground support equipment (GSE) that could be installed whenever measurements of the star tracker optical references were needed. However, this GSE could only be used to measure secondary optical reference cube faces not used by the star tracker vendor to obtain the relationship information and matrix transformations necessary to determine star tracker alignment. Unfortunately, due to unexpectedly large orthogonality errors between the measured secondary adjacent cube faces and the lack of cube calibration data, we required a method that could be used to measure the same reference cube faces as originally measured by the vendor. We describe an alternative technique to theodolite auto-collimation for measurement of an optical reference mirror pointing direction when normal incidence measurements are not possible. This technique was used to successfully align the GPM star trackers and has been used on a number of other NASA flight projects. We also discuss alignment theory as well as a GSFC-developed theodolite data analysis package used to analyze angular metrology data.

Hetherington, Samuel↗
DOE OSTI2024

Systematic Uncertainties from Gribov Copies in Lattice Calculation of Parton Distributions in the Coulomb Gauge

Recently, a new method has been proposed to compute parton distributions using boosted correlators fixed in the Coulomb gauge (CG) within the framework of large-momentum effective theory. This approach, which does not involve Wilson lines, could greatly improve the efficiency and precision of lattice quantum chromodynamics calculations. However, concerns remain regarding whether systematic uncertainties from Gribov copies, which correspond to ambiguities in lattice gauge-fixing, are adequately controlled. This work assesses the effects of Gribov copies on Coulomb-gauge-fixed quark correlators. We utilize different strategies for Coulomb-gauge fixing, selecting two different groups of Gribov copies based on lattice gauge configurations. We examine the differences in the resulting spatial quark correlators in both vacuum and pion states. Our findings indicate that the statistical errors of the matrix elements from both Gribov copies, regardless of the correlation range, decrease proportionally to the square root of the number of gauge configurations. The difference between the strategies does not show statistical significance compared to the gauge noise, demonstrating that the effect of the Gribov copies can be neglected in practical lattice calculations of quark parton distributions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Tunable-fidelity wave functions for the ab initio description of scattering and reactions

Here, the no-core shell model (NCSM) is an ab initio method that solves the nuclear many-body problem by expanding the many-particle wave function into a (typically) harmonic oscillator basis and minimizing the energy to obtain the expansion coefficients. Extensions of the NCSM, such as its coupling with microscopic-cluster basis states, further allow for an ab initio treatment of light-ion nuclear reactions of interest for both astrophysics and nuclear technology applications. A downside of the method is the exponential scaling of the basis size with increasing number of nucleons and excitation quanta, which limits its applicability to mass A ≲ 16 nuclei, except for variants where the basis is further down-selected via some truncation scheme. We consider a basis selection method for the NCSM that was first introduced in the context of the large-scale shell model and captures the essential degrees of freedom of the nuclear wave function leading to a favorable complexity scaling for calculations and enabling ab initio reaction calculations in sd-shell nuclei. The particle configurations within the NCSM basis are ordered based on their contribution to the first moment of the Hamiltonian matrix that results from the projection onto the many-body basis. The truncation scheme then consists in retaining only the lowest-first-moment configurations, which typically contain only few many-body basis states (Slater determinants). As the energy threshold above which configurations are disregarded is increased, the size of the basis becomes an almost-continuous variable, allowing for tunable fidelity in the obtained wave functions. The resulting wave functions can then be used directly in ab initio reaction calculations. We present calculations for 7 Li and n + 12 C scattering using nucleon-nucleon interactions derived from chiral effective field theory and softened using the similarity renormalization group method. The obtained energy levels invariably demonstrate exponential convergence with the size of the basis, and we find improved convergence in scattering calculations. To demonstrate the possibilities enabled by the approach, we also present a first calculation for the scattering of neutrons from 24 Mg. The method presented in this work appears promising for future studies of nuclei with mass A > 16, opening multiple future research directions impacting both nuclear astrophysics and nuclear technology applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1989

Application of the pseudostate theory to the three-body Lambert problem

The pseudostate theory, which approximates three-body trajectories by overlapping the conic effects of both massive bodies on the third body, has been used to solve boundary-value problems. Frequently, the approach to the secondary is quite close, as in interplanetary gravity-assist or satellite-tour trajectories. In this case, the orbit with respect to the primary is radically changed so that perturbation techniques are time consuming, yet higher accuracy than point-to-point conics is necessary. This method reduces the solution of the three-body Lambert problem to solving two conic Lambert problems and inverting a 7 x 7 matrix, the components of which are all found analytically. Typically 90-95 percent of the point-to-point conic error, with respect to an integrated trajectory, is eliminated.

Byrnes, Dennis V.↗
NASA NTRS2008

Reliability and Confidence Interval Analysis of a CMC Turbine Stator Vane

High temperature ceramic matrix composites (CMC) are being explored as viable candidate materials for hot section gas turbine components. These advanced composites can potentially lead to reduced weight, enable higher operating temperatures requiring less cooling and thus leading to increased engine efficiencies. However, these materials are brittle and show degradation with time at high operating temperatures due to creep as well as cyclic mechanical and thermal loads. In addition, these materials are heterogeneous in their make-up and various factors affect their properties in a specific design environment. Most of these advanced composites involve two- and three-dimensional fiber architectures and require a complex multi-step high temperature processing. Since there are uncertainties associated with each of these in addition to the variability in the constituent material properties, the observed behavior of composite materials exhibits scatter. Traditional material failure analyses employing a deterministic approach, where failure is assumed to occur when some allowable stress level or equivalent stress is exceeded, are not adequate for brittle material component design. Such phenomenological failure theories are reasonably successful when applied to ductile materials such as metals. Analysis of failure in structural components is governed by the observed scatter in strength, stiffness and loading conditions. In such situations, statistical design approaches must be used. Accounting for these phenomena requires a change in philosophy on the design engineer s part that leads to a reduced focus on the use of safety factors in favor of reliability analyses. The reliability approach demands that the design engineer must tolerate a finite risk of unacceptable performance. This risk of unacceptable performance is identified as a component's probability of failure (or alternatively, component reliability). The primary concern of the engineer is minimizing this risk in an economical manner. The methods to accurately determine the service life of an engine component with associated variability have become increasingly difficult. This results, in part, from the complex missions which are now routinely considered during the design process. These missions include large variations of multi-axial stresses and temperatures experienced by critical engine parts. There is a need for a convenient design tool that can accommodate various loading conditions induced by engine operating environments, and material data with their associated uncertainties to estimate the minimum predicted life of a structural component. A probabilistic composite micromechanics technique in combination with woven composite micromechanics, structural analysis and Fast Probability Integration (FPI) techniques has been used to evaluate the maximum stress and its probabilistic distribution in a CMC turbine stator vane. Furthermore, input variables causing scatter are identified and ranked based upon their sensitivity magnitude. Since the measured data for the ceramic matrix composite properties is very limited, obtaining a probabilistic distribution with their corresponding parameters is difficult. In case of limited data, confidence bounds are essential to quantify the uncertainty associated with the distribution. Usually 90 and 95% confidence intervals are computed for material properties. Failure properties are then computed with the confidence bounds. Best estimates and the confidence bounds on the best estimate of the cumulative probability function for R-S (strength - stress) are plotted. The methodologies and the results from these analyses will be discussed in the presentation.

Murthy, Pappu L. N.↗
NASA NTRS1979

Are large concentration of atomic H storable in tritium-impregnated solid in H2 below 0.10 K

The storage and release of atomic hydrogen produced by the beta decay of tritium contained in a crystalline solid H2 matrix at concentrations greater than 2% and temperatures below 0.80 K are investigated. The temperature of a sample chamber containing tritium-impregnated H2 and placed in the mixing chamber of a dilution refrigerator was measured as the chamber was heated and cooled in order to determine the rates of energy storage and release. It is found that for samples containing 1.2 wt.% tritium, after storage at 0.054 K for 40 h, an increase in sample temperature to a trigger point of 0.17 K leads to an energy release due to the destabilization of atomic H in H2 as predicted by the phenomenological rate process theory. For a tritium weight fraction of 2.5%, energy releases were triggered at 0.54 and 0.82 K after storage at 0.080 K, indicating the trapping of H atoms at the sites of T2 and HT molecules in the sample. The application of a 15 kG magnetic field is shown to increase the storage capacity of T2 traps while reducing that of HT traps, and to lower the trigger temperatures of both. Results suggest that the direct conversion of nuclear energy to chemical energy may become technically feasible in the future.

Rosen, G.↗
DOE OSTI2024

Nonrelativistic nuclear reduction for tensor couplings in dark matter direct detection and μ → e conversion

The nonrelativistic effective field theory (NRET) is widely used in dark matter direct detection and charged-lepton flavor violation studies through μ → e conversion. However, existing literature has not fully considered tensor couplings. This study fills this gap by utilizing an innovative tensor decomposition method, extending NRET to incorporate previously overlooked tensor interactions. This development is expected to have a significant impact on ongoing experiments seeking physics beyond the Standard Model and on our understanding of the new-physics interactions. Notably, we identify additional operators in μ → e conversion that are absent in scalar and vector couplings. To support further research and experimental analyses, comprehensive tables featuring tensor matrix elements and their corresponding operators are provided. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
NASA NTRS1979

Application of the pseudostate theory to the three-body Lambert problem

The pseudostate theory, which approximates three-body trajectories by overlapping the conic effects of both massive bodies on the third body, has been used to solve boundary value problems. Frequently, the approach to the secondary is quite close, as in interplanetary gravity assist trajectories or satellite tour trajectories. In this case the orbit with respect to the primary is radically changed so that perturbation techniques are time consuming, yet higher accuracy than point-to-point conics (V-infinity matching) is necessary. This method reduces the solution of the three-body Lambert problem to solving two conic Lambert problems and inverting a 7 x 7 matrix, the components of which are all found analytically. Typically 90-95% of the point-to-point conic error with respect to an integrated trajectory is eliminated.

Byrnes, D. V.↗
NASA NTRSDate not supplied

Machine Learning Approaches for Rare-Earth Silicate Environmental Barrier Coating Thermochemical and Thermomechanical Property Predictions

Environmental barrier coatings (EBCs) are a necessary enabling technology for the transition from superalloys to silicon carbide (SiC) ceramic matrix composites (CMCs) in gas turbine engines for increased efficiency and decreased fuel costs. SiC-based CMCs are prone to oxidation-based degradation in the engine hot section, and rare-earth (RE) silicates are promising candidates for EBCs due to their close thermal expansion match to the composite substrate and oxidation resistance. However, the design of EBCs is hindered by the large chemical space of candidate materials and the difficulty in obtaining material properties for engineering optimization. This is especially difficult as research continues into mixed-cation or “high-entropy” RE silicates. First-principles computational methods such as density functional theory (DFT) are highly effective at calculating material properties to guide coating design but are limited by their computational cost. Atomistic simulations have the potential to both accelerate property calculations and expand the properties able to be calculated due to their lower computational compared to DFT. However, they require interatomic potentials (IAPs) specific to the material system of interest, and, to our knowledge, there are no suitable IAPs for RE silicates. Machine learning (ML) is a promising technique to accelerate material property predictions indirectly by generating IAPs for atomistic simulations or via direct prediction. In this work, we present two ML approaches to accelerate the calculation of RE silicate properties relevant to EBC design: 1) a ML-derived interatomic potential (IAP) for atomistic simulations of yttrium disilicate (Y2Si2O7) from DFT training data, and 2) a neural network (NN) model to directly predict thermochemical properties of RE silicates and oxides directly from easily obtainable unit cell parameters. Classical MD simulations using the IAP yield lattice properties and bond lengths in good agreement with both DFT and experimental results from x-ray diffraction. Thermodynamic properties calculated using the finite-displacement phonon method and quasi-harmonic approximation were orders of magnitude faster than DFT with good agreement to the DFT results. The IAP was also used to calculate properties such as coefficient of thermal expansion (CTE) that require large simulation supercells and are therefore difficult with DFT. The IAP correctly predicted the anisotropic nature of the CTE in three different phases of Y2Si2O7. The NN model predicts constant pressure heat capacity, Cp, orders of magnitude faster than DFT calculations, which can enable its use as a surrogate model for multiscale simulations. The two methods presented in this work demonstrate the utility of ML for accelerating the prediction of RE silicate properties, which can in turn accelerate EBC design and optimization.

machine learning↗
DOE OSTI2020

Elastic Antiproton-Nucleus Scattering from Chiral Forces

Elastic scattering of antiprotons off He 4 , C 12 , and O 16 , 18 is described for the first time with a consistent microscopic approach based on the calculation of an optical potential (OP) describing the antiproton-target interaction. The OP is derived using the recent antiproton-nucleon ( p ¯ N ) chiral interaction to calculate the p ¯ N t matrix, while the target densities are computed with the ab initio no-core shell model using chiral interactions as well. Here, our results are in good agreement with the existing experimental data and the results computed at different chiral orders of the p ¯ N interaction display a well-defined convergence pattern.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Remote detectability from entanglement bootstrap I: Kirby’s torus trick

Remote detectability is often taken as a physical assumption in the study of topologically ordered systems, and it is a central axiom of mathematical frameworks of topological quantum field theories. We show under the entanglement bootstrap approach that remote detectability is a necessary property; that is, we derive it as a theorem. Starting from a single wave function on a topologically-trivial region satisfying the entanglement bootstrap axioms, we can construct states on closed manifolds. The crucial technique is to immerse the punctured manifold into the topologically trivial region and then heal the puncture. This is analogous to Kirby’s torus trick. We then analyze a special class of such manifolds, which we call pairing manifolds. For each pairing manifold, which pairs two classes of excitations, we identify an analog of the topological S S -matrix. This pairing matrix is unitary, which implies remote detectability between two classes of excitations. These matrices are in general not associated with the mapping class group of the manifold. As a by-product, we can count excitation types (e.g., graph excitations in 3+1d). The pairing phenomenon occurs in many physical contexts, including systems in different dimensions, with or without gapped boundaries. We provide a variety of examples to illustrate its scope.

Shi, Bowen (ORCID:0000000206899964)↗
NASA NTRS2000

Ultrasonic Assessment of Impact-Induced Damage and Microcracking in Polymer Matrix Composites

The main objective of this NASA FAR project is to conduct ultrasonic assessment of impact-induced damage and microcracking in polymer matrix composites at various temperatures. It is believed that the proposed study of impact damage assessment on polymer matrix composites will benefit several NASA's missions and current interests, such as ballistic impact testing of composite fan containment and high strain rate deformation modeling of polymer matrix composites. Currently, impact-induced delamination and fracture in 6061-T6 aluminum/cast acrylic sandwich plates adhered by epoxy were generated in an instrumented drop-weight impact machine. Although only a small dent was produced on the aluminum side when a hemispherical penetrator tup was dropped onto it from a couple of inches, a large ring of delamination at the interface was observed. The delamination damage was often accompanied by severe shattering in the acrylic substratum. Damage patterns in the acrylic layer include radial and ring cracks and, together with delamination at the interface, may cause peeling-off of acrylic material from the sandwich plate. Theory of stress-wave propagation can be used to explain these damage patterns. The impact tests were conducted at various temperatures. The results also show clearly that temperature effect is very important in impact damage. For pure cast acrylic nil-ductile transition (NDT) occurs between 185-195 F Excessive impact energy was dissipated into fracture energy when tested at temperature below this range or through plastic deformation when tested at temperature above the NDT temperature. Results from this study will be used as baseline data for studying fiber-metal laminates, such as GLARE and ARALL for advanced aeronautical and astronautical applications.

Liaw, Benjamin↗
NASA NTRS2004

Lunar and Planetary Science XXXV: Big Dust, Little Dust, and Aerogel

The titles in this section include: 1) A First Look at Graphite Grains from Orgueil: Morphology, Carbon, Nitrogen and Neon Isotopic Compositions of Individual, Chemically Separated Grains; 2) Discovery of an In-Situ Presolar Silicate Grain with GEMS-Like Composition in the Bishumpur Matrix; 3) Revised Global Accretion Rates of Micrometeorites in the Last Glacial Period; 4) Isotopic Fractionation of Potassium in Stony Cosmic Spherules; 5) Crumbs from the Crust of Vesta: Achondritic Micrometeorites from the South Pole Water Well; 6) Preservation of Micrometeoroid Remains on Solar Cells Returned from Orbit; 7) High Input Rates of Micrometeoritic Sulfur, 'Smoke' Particles and Oligoelements on the Early Earth; 8) High Resolution Confocal Raman Imaging of an IDP; 9) A micro-Raman Survey of 10 IDPs and 6 Carbonaceous Chondrites; 10) Mid- and Far-Infrared Spectroscopy at the Advanced Light Source; 11) Hypervelocity Impact Energy Loss and Track Shape in Aerogels: Theory and Experiment; 12) Fiducial Marks for Location of Particles in Aerogel.

Source record↗