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,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Behaviour of plate elements based on the first-order shear deformation theory

A new analytical technique to assess the performance of shear deformable elements is presented, which makes it possible to determine a priori whether a given element will lock when used to model thin structures. The role that shear constraints play in determining the behavior of thin elements was established by comparing the results of key numerical tests with the predictions of element behavior made by studying the form of the shear constraints. Conclusions regarding locking behavior and the effects of reduced integration in thin shear deformable elements are presented, including the findings (1) that singularity of the shear stiffness matrix is not sufficient to avoid locking; (2) that the effect of mesh refinement on an element that contains spurious constraints is two-fold; and (3) that reduced integration does not remove spurious constraints but rather relaxes them. The results of the study are in agreement with previous studies of Mindlin plate elements in regarding Lagrangian elements as superior to serendipity elements when either full or reduced integration is employed.

Averill, R. C.↗
NASA NTRS1987

Nonlinear Beam Kinematics by Decomposition of the Rotation Tensor

A simple matrix expression is obtained for the strain components of a beam in which the displacements and rotations are large. The only restrictions are on the magnitudes of the strain and of the local rotation, a newly-identified kinematical quantity. The local rotation is defined as the change of orientation of material elements relative to the change of orientation of the beam reference triad. The vectors and tensors in the theory are resolved along orthogonal triads of base vectors centered along the undeformed and deformed beam reference axes, so Cartesian tensor notation is used. Although a curvilinear coordinate system is natural to the beam problem, the complications usually associated with its use are circumvented. Local rotations appear explicitly in the resulting strain expressions, facilitating the treatment of beams with both open and closed cross sections in applications of the theory. The theory is used to obtain the kinematical relations for coupled bending, torsion extension, shear deformation, and warping of an initially curved and twisted beam.

Danielson, D. A.↗
DOE OSTI2022

Weyl semimetal in the rare-earth hexaboride family supporting a pseudonodal surface and a giant anomalous Hall effect

Rare-earth hexaborides offer a rich tapestry for exploring the interplay of topological orders, magnetism, and electron-correlation effects associated with protected quantum phenomena. Using first-principles modeling, we identify in cerium hexaboride (CeB 6 ) a topological state that harbors a pseudonodal surface (PNS) along with double-Weyl fermions. An analysis of the electronic states in the vicinity of the Fermi level reveals the presence of M z -mirror-symmetry-protected band crossings that involve itinerant (d) and localized (f) states of Ce. These band crossings are found to remain nearly gapless as one goes away from the mirror plane, driving the formation of PNS states. We also demonstrate the presence of coexisting double-Weyl fermions with chiral charges of ±2 on the C 4z rotational axis. By analyzing the Berry curvature field associated with the topological states, we predict that CeB 6 supports a giant anomalous Hall conductivity of -1840Ω -1 cm -1 in the slightly electron-doped regime, which would be larger than the values reported to date in magnetic Weyl semimetals. Our study thus indicates that CeB 6 , which combines ferromagnetic order with topological physics in a heavy-fermion matrix, would provide a promising playground for investigating topological magnetic states with giant anomalous Hall responses.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A low-energy perspective on the minimal left-right symmetric model

We perform a global analysis of the low-energy phenomenology of the minimal left-right symmetric model (mLRSM) with parity symmetry. We match the mLRSM to the Standard Model Effective Field Theory Lagrangian at the left-right-symmetry breaking scale and perform a comprehensive fit to low-energy data including mesonic, neutron, and nuclear β -decay processes, Δ F = 1 and Δ F = 2 CP-even and -odd processes in the bottom and strange sectors, and electric dipole moments (EDMs) of nucleons, nuclei, and atoms. We fit the Cabibbo-Kobayashi-Maskawa and mLRSM parameters simultaneously and determine a lower bound on the mass of the right-handed W R boson. In models where a Peccei-Quinn mechanism provides a solution to the strong CP problem, we obtain $ {M}_{W_R} $ ≳ 5 . 5 TeV at 95% C.L. which can be significantly improved with next-generation EDM experiments. In the P -symmetric mLRSM without a Peccei-Quinn mechanism we obtain a more stringent constraint $ {M}_{W_R} $ ≳ 17 TeV at 95% C.L., which is difficult to improve with low-energy measurements alone. In all cases, the additional scalar fields of the mLRSM are required to be a few times heavier than the right-handed gauge bosons. We consider a recent discrepancy in tests of first-row unitarity of the CKM matrix. We find that, while TeV-scale W R bosons can alleviate some of the tension found in the V ud,us determinations, a solution to the discrepancy is disfavored when taking into account other low-energy observables within the mLRSM.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Closo ‐Borate Gel Polymer Electrolyte with Remarkable Electrochemical Stability and a Wide Operating Temperature Window

Abstract A major challenge in the pursuit of higher‐energy‐density lithium batteries for carbon‐neutral‐mobility is electrolyte compatibility with a lithium metal electrode. This study demonstrates the robust and stable nature of a closo ‐borate based gel polymer electrolyte (GPE), which enables outstanding electrochemical stability and capacity retention upon extensive cycling. The GPE developed herein has an ionic conductivity of 7.3 × 10 −4 S cm −2 at room temperature and stability over a wide temperature range from −35 to 80 °C with a high lithium transference number ( = 0.51). Multinuclear nuclear magnetic resonance and Fourier transform infrared are used to understand the solvation environment and interaction between the GPE components. Density functional theory calculations are leveraged to gain additional insight into the coordination environment and support spectroscopic interpretations. The GPE is also established to be a suitable electrolyte for extended cycling with four different active electrode materials when paired with a lithium metal electrode. The GPE can also be incorporated into a flexible battery that is capable of being cut and still functional. The incorporation of a closo ‐borate into a gel polymer matrix represents a new direction for enhancing the electrochemical and physical properties of this class of materials.

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

Electronegativity-Guided Site Differentiation in High-Entropy Alloy for pH-universal Hydrogen Evolution Reactions

Enhancing the intrinsic activity of transition metal catalysts for the hydrogen evolution reaction (HER) remains a critical challenge in sustainable energy conversion. Herein, we report an electronegativity-guided site differentiation strategy in a single-phase CoNiCuMoW high-entropy alloy (HEA) via electrodeposition by incorporating high-electronegativity 4d/5d orbital transition metals (Mo, W) into the face-centered cubic (fcc) matrix (CoNiCu). The as-synthesized HEA demonstrates exceptional HER performance in all pH conditions, delivering an outstanding overpotential of 65 mV (alkaline), 28 mV (acidic), and 155 mV (neutral) at a current density of 100 mA cm−2, showing performance comparable to commercial Pt/C and has excellent long-term stability at high current density (1 A cm−2, 1000 h). X-Ray absorption spectroscopy (XAS) and density functional theory (DFT) calculations reveal that the incorporation of Mo/W simultaneously alters the local coordination environment and induces element-dependent charge redistribution, accompanied by a system-level d-band center downshift, thereby optimizing the hydrogen binding strength across multimetallic sites. Meanwhile, oxophilic Mo/W sites lower the water dissociation energy barrier. These synergistic effects collectively enable efficient and durable pH-universal HER performance.

Wu, Yutong↗
DOE OSTI2023

Electronic heat generation in semiconductors: Non-equilibrium excitation and evolution of zone-edge phonons via electron–phonon scattering in photo-excited germanium

We investigate experimentally and using first-principles theory the generation of phonons and the relaxation of carriers on picosecond timescales across the Brillouin zone of photo-excited Ge by inter-valley electron–phonon scattering. The phonons generated are typical of those generated in semiconductor devices, contributing to the accumulation of heat within the material. We simulate the time-evolution of phonon populations, based on first-principles band structure and electron–phonon and phonon–phonon matrix elements, and compare them to data from time-resolved x-ray diffuse scattering experiments, performed at the Linac Coherent Light Source x-ray free-electron laser facility, following photo-excitation by a 50 fs near-infrared optical pulse. We show that the intensity of the non-thermal x-ray diffuse scattering signal, which is observed to grow substantially near the L-point of the Brillouin zone over 3–5 ps, is due to phonons generated by scattering of carriers between the Δ and L valleys. These phonons have low group velocities, resulting in a heat bottleneck. With the inclusion of phonon decay through 3-phonon processes, the simulations also account for other non-thermal features observed in the x-ray diffuse scattering intensity, which are due to anharmonic phonon–phonon scattering of the phonons initially generated by electron–phonon scattering.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

A flexible gyro-fluid system of equations

Gyro-fluid equations are velocity space moments of the gyrokinetic equations. Special gyro-Landau-fluid closures have been developed that include the damping due to kinetic resonances by fitting to the collisionless local plasma response functions. This damping allows for accurate linear eigenmodes to be computed with a relatively low number of velocity space moments compared to the number of velocity quadrature points in gyrokinetic codes. However, none of the published gyro-Landau-fluid closure schemes considers the Onsager symmetries of the resulting quasi-linear fluxes as a constraint. Onsager symmetry guarantees that the matrix of diffusivities is positive definite, an important property for the numerical stability of a transport solver. A two-parameter real closure for improving the accuracy of low-resolution gyro-fluid equations, which preserves the Onsager symmetry and allows higher velocity space moments, is presented in this paper. The new linear gyro-fluid system (GFS) is used to extend the TGLF quasi-linear transport model so that it can compute the energy and momentum fluxes due to parallel magnetic fluctuations, completing the transport matrix. The GFS equations do not use a bounce average approximation. The GFS equations are fully electromagnetic with general flux surface magnetic geometry, pitch angle scattering for electron collisions, and subsonic equilibrium toroidal rotation. Using GFS eigenmodes in the quasi-linear TGLF model will be shown to yield a more accurate match to fluxes computed by CGYRO turbulence simulations. In conclusion, prospects for future applications of a quasi-linear theory to new plasma transport regimes and magnetic confinement devices in addition to tokamaks are opened by the flexibility of the GFS eigensolver.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Modeling of Magnetic Properties of Magnetorheological Elastomers Using JA Hysteresis Model

Magnetorheological elastomers (MREs) are composite materials that consist of magnetically permeable particles in a nonmagnetic polymeric matrix. Under the influence of an external magnetic field, a reversible deformation change occurs in the mechanical properties of these materials. Due to their coupled magnetomechanical response, these materials have been found suitable for a range of applications including tunable vibration absorbers, sensors, and actuators. Notably, improvement of such devices are prerequisites to efficient energy conversion systems, hence the need to understand further the MRE technology. The Jiles-Atherton (JA) theory takes into consideration the magneto-coupling experienced by effective domains in a magnetic material. Algorithm based on the theory yields five model parameters; saturation magnetization (M s ), domain density (a), domain coupling (α), loss coefficient (k), and reversibility (c). Using JA theory, model parameters were calculated and linked to the physical attributes of Fe powder and isotropic MRE. The results show that the calculated Ms for the MRE is reasonably related to that of the Fe powder by a factor of the particle's volume fraction used in the MRE. Further, the calculated k, a, and α provided support for the reduced pinning factor, domain density, and increased domain coupling in the MRE due to the changes in the domain structure between the two materials. From the calculated JA parameters, finite-element modeling (FEM) of the MRE hysteresis loop was performed. The analysis showed that the modeled magnetic properties including coercivity, remanence, and coordinates of the hysteresis loop tip vary with geometric position.

42 ENGINEERING↗
NASA NTRS1974

Boron/aluminum-graphite/resin advanced fiber composite hybrids

An investigation was conducted to determine the fabrication feasibility and to assess the potential of adhesively-bonded metal and resin matrix fiber composite hybrids as an advanced material, for aerospace and other structural applications. The results of fabrication studies and of evaluation of physical and mechanical properties show that using this hybrid concept it is possible to design a composite which, when compared to nonhybrid composites, has improved transverse strength, transverse stiffness, and impact resistance with only a small penalty on density and longitudinal properties. The results also show that laminate theory is suitable for perdicting the structural response of such hybrids. The sequence of fracture modes indicates that these types of hybrids can be readily designed to meet fail-safe requirements.

Chamis, C. C.↗
NASA NTRS1974

Boron/aluminum-graphite/resin advanced fiber composite hybrids

An investigation was conducted to determine the fabrication feasibility and to assess the potential of adhesively-bonded metal and resin matrix fiber composite hybrids as an advanced material, for aerospace and other structural applications. The results of fabrication studies and of evaluation of physical and mechanical properties show that using this hybrid concept it is possible to design a composite which, when compared to nonhybrid composites, has improved transverse strength, transverse stiffness, and impact resistance with only a small penalty on density and longitudinal properties. The results also show that laminate theory is suitable for predicting the structural response of such hybrids. The sequence of fracture modes indicates that these types of hybrids can be readily designed to meet fail-safe requirements.

Chamis, C. C.↗
NASA NTRS1979

A high order accurate numerical solution algorithm for turbulent boundary layer flow

A fourth-order accurate numerical solution algorithm is derived using finite element interpolation theory for the non-linear parabolic equations governing turbulent boundary layer flow including a two-equation turbulence closure model. The results of carefully controlled numerical experiments firmly quantize for the first time performance differences between finite element and finite difference solution methodology for this type of equation. The developed algorithm takes advantage of the apparent semi-analytical formulational procedure, in establishment of a single, retarded-evaluation Jacobian matrix iterative solution algorithm. Numerical results document performance of solution economy features in terms of computer requirements and solution accuracy. The developed algorithm should find wide application in aerodynamics analysis.

Soliman, M. O.↗
DOE OSTI2021

Islands in de Sitter space

We consider black holes in 2d de Sitter JT gravity coupled to a CFT, and entangled with matter in a disjoint non-gravitating universe. Tracing out the entangling matter leaves the CFT in a density matrix whose stress tensor backreacts on the de Sitter geometry, lengthening the wormhole behind the black hole horizon. Naively, the entropy of the entangling matter increases without bound as the strength of the entanglement increases, but the monogamy property predicts that this growth must level off. We compute the entropy via the replica trick, including wormholes between the replica copies of the de Sitter geometry, and find a competition between conventional field theory entanglement entropy and the surface area of extremal “islands” in the de Sitter geometry. The black hole and cosmological horizons both play a role in generating such islands in the backreacted geometry, and have the effect of stabilizing the entropy growth as required by monogamy. We first show this in a scenario in which the de Sitter spatial section has been decompactified to an interval. Then we consider the compact geometry, and argue for a novel interpretation of the island formula in the context of closed universes that recovers the Page curve. Finally, we comment on the application of our construction to the cosmological horizon in empty de Sitter space.

2D gravity↗
DOE OSTI2025

Efficient use of quantum computers for collider physics

Most observables at particle colliders involve physics at a wide variety of distance scales. Due to asymptotic freedom of the strong interaction, the physics at short distances can be calculated reliably using perturbative techniques, while long distance physics is non-perturbative in nature. Factorization theorems separate the contributions from different scales, allowing to identify the pieces that can be determined perturbatively from those that require non-perturbative information, and if the non-perturbative pieces can be reliably determined, one can use experimental measurements to extract the short distance effects, sensitive to possible new physics. Without the ability to compute the non-perturbative ingredients from first principles one typically identifies observables for which the non-perturbative information is universal in the sense that it can be extracted from some experimental observables and then used to predict other observables. In this paper we argue that the future ability to use quantum computers to calculate non-perturbative matrix elements from first principles will allow to make predictions for observables with non-universal non-perturbative long-distance physics.

Algorithms and Theoretical Developments↗
DOE OSTI2026

An Algorithm for Atom-Centered Lossy Compression of the Atomic Orbital Basis in Density Functional Theory Calculations

Large atomic-orbital (AO) basis sets of at least triple and preferably quadruple-ζ (QZ) size are required to adequately converge Kohn–Sham density functional theory (DFT) calculations toward the complete basis set limit. However, incrementing the cardinal number by one nearly doubles the AO basis dimension, and the computational cost scales as the cube of the AO dimension, so this is very computationally demanding. Here, in this work, we develop and test a threshold-based natural atomic orbital (NAO) scheme in which ϵ-NAOs are obtained as eigenfunctions of atomic blocks of the density matrix in a one-center orthogonalized representation. This enables compression of the AO basis that is optimal for a given threshold, 10 –ϵ , by discarding NAOs with occupation numbers below that threshold. Extensive pilot test calculations using the Hartree–Fock functional and taking the converged density matrix as input suggest that a threshold of 10 –5 can yield a compression factor (ratio of AO to compressed ϵ-NAO dimension) between 2.5 and 4.5 for the QZ pc-3 basis. The errors in relative energies are typically less than 0.1 kcal/mol when the compressed basis is used instead of the uncompressed basis. Between 10 and 100 times smaller errors (i.e., usually less than 0.01 kcal/mol) can be obtained with a threshold 10 –7 , while the compression factor is typically between 2 and 2.5.

basis sets↗
NASA NTRS1988

Hot piston ring/cylinder liner materials: Selection and evaluation

In current designs of the automotive (kinematic) Stirling engine, the piston rings are made of a reinforced polymer and are located below the pistons because they cannot withstand the high temperatures in the upper cylinder area. Theoretically, efficiency could be improved if hot piston rings were located near the top of the pistons. Described is a program to select piston ring and cylinder coating materials to test this theory. Candidate materials were screened, then subjected to a pin or disk friction and wear test machine. Tests were performed in hydrogen at specimen temperatures up to 760 C to simulate environmental conditions in the region of the hot piston ring reversal. Based on the results of these tests, a cobalt based alloy, Stellite 6B, was chosen for the piston rings and PS200, which consists of a metal-bonded chromium carbide matrix with dispersed solid lubricants, was chosen as the cylinder coating. Tests of a modified engine and a baseline engine showed that the hot ring reduced specific fuel consumption by up to 7 percent for some operating conditions and averaged about 3 percent for all conditions evaluated. Related applications of high-temperature coatings for shaft seals and as back-up lubricants are also described.

Sliney, Harold E.↗
NASA NTRS1988

Failure of compression-loaded multi-directional composite laminates

An analytical and experimental investigation of the failure of selected compression-loaded composite laminates is described. A general nonlinear theory is presented for predicting a laminate's compressive strength and failure mode. The theory includes the effects of out-of-plane ply waviness, inplane fiber waviness, and fiber scissoring. A simple compressive test technique is used to obtain the experimental data. The analytical and experimental results show good agreement for theta less than 45 deg and show excellent agreement for theta not less than 45 deg. The dominant compression failure modes for the laminates in this study were found to be interlaminar shearing, inplane matrix shearing, and matrix compression.

Shuart, Mark J.↗
NASA NTRS1987

Spatially random models, estimation theory, and robot arm dynamics

Spatially random models provide an alternative to the more traditional deterministic models used to describe robot arm dynamics. These alternative models can be used to establish a relationship between the methodologies of estimation theory and robot dynamics. A new class of algorithms for many of the fundamental robotics problems of inverse and forward dynamics, inverse kinematics, etc. can be developed that use computations typical in estimation theory. The algorithms make extensive use of the difference equations of Kalman filtering and Bryson-Frazier smoothing to conduct spatial recursions. The spatially random models are very easy to describe and are based on the assumption that all of the inertial (D'Alembert) forces in the system are represented by a spatially distributed white-noise model. The models can also be used to generate numerically the composite multibody system inertia matrix. This is done without resorting to the more common methods of deterministic modeling involving Lagrangian dynamics, Newton-Euler equations, etc. These methods make substantial use of human knowledge in derivation and minipulation of equations of motion for complex mechanical systems.

Rodriguez, G.↗