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 811 records · Page 45

Quantum Melting of Generalized Wigner Crystals in Transition Metal Dichalcogenide Moiré Systems

The generalized Wigner crystal (GWC) is a novel quantum phase of matter driven by further-range interaction at fractional fillings of a lattice. The role of further-range interaction as the driver for the incompressible state is akin to the Wigner crystal. On the other hand, the significant role of commensurate filling is akin to the Mott insulator. Recent progress in simulator platforms presents unprecedented opportunities to investigate quantum melting in the strongly interacting regime through synergy between theory and experiments. However, the earlier theory literature presents diverging predictions. We study the quantum freezing of GWC through large-scale density matrix renormalization group simulations of a triangular lattice extended Hubbard model. We find a single first-order phase transition between the Fermi liquid and the √3×√3 GWC state. Furthermore, the GWC state shows long-range antiferromagnetic 120° Néel order. Our results present the simplest answers to the question of the quantum phase transition into the GWC phase and the properties of the GWC phase.

Density matrix renormalization group↗

A Modified Network Alteration Model to Predict Quasi-Static Behavior of the Cross-Linked Polymers During Hydrolytic Aging

Abstract The network alteration theory is utilized to capture the effects of hydrolytic aging on the mechanical integrity of the polymer matrix. Hydrolytic degradation generally considered as a result of two simultaneous chemical phenomena; (i) reduction of the cross-links, and (ii) active chain loss. Assuming the alteration of the chain density is identical to the change of the number average molecular weight along the aging trajectory, the strain energy of polymer matrix can be rewritten as a function of deformation, deformation history, storage time, and aging temperature. The proposed model includes only five physically inspired material parameters. Thus, while it is computationally efficient, it shows good agreement with own experimental data, which performed on various range of accelerated aging temperatures and times. With respect to its computational efficiency, simplicity, accuracy, and interpret-ability, the model is a good choice for advanced implementations in FE programs.

Bahrololoumi, Amir↗

Minimum impulse three-body trajectories.

A rapid and accurate method of calculating optimal impulsive transfers in the restricted problem of three bodies has been developed. The technique combines a multi-conic method of trajectory integration with primer vector theory and an accelerated gradient method of trajectory optimization. A unique feature is that the state transition matrix and the primer vector are found analytical without additional integrations or differentiations. The method has been applied to the determination of optimal two and three impulse transfers between the L2 libration point and circular orbits about both the earth and the moon.

D'Amario, L.↗

Predicting the time-temperature dependent axial failure of B/A1 composites

Experimental and theoretical studies were conducted in order to understand and predict the effects of time, temperature, and stress on the axial failure modes of boron fibers and B/A1 composites. Due to the anelastic nature of boron fiber deformation, it was possible to determine simple creep functions which can be employed to accurately describe creep and fracture stress of as-produced fibers. Analysis of damping and strength data for B/6061 A1 composites indicates that fiber creep effects of creep on fiber fracture are measurably reduced by the composite fabrication process. The creep function appropriate for fibers with B/Al composites was also determined. A fracture theory is presented for predicting the time-temperature dependence of the axial tensile strength for metal matrix composites in general and B/A1 composites in particular.

Dicarlo, J. A.↗

Recursive Robot-Arm Dynamics via Filtering and Smoothing

Forward and inverse dynamics solved using Kalman filtering and Bryson-Frazier smoothing. Dynamics of serial-link robot arm solved by using recursive techniques from linear filtering and smoothing theory. Solutions of dynamical equations give forces, moments, and accelerations at joints between links, and multilink inertia matrix and its inverse. Theoretical developments lay foundation for use of filtering and smoothing techniques in design of robot controls.

Rodriguez, Guillermo↗

Stabilization of Taylor-Couette flow due to time-periodic outer cylinder oscillation

The linear stability of circular Couette flow between concentric infinite cylinders is considered for the case when the inner cylinder is rotated at a constant angular velocity and the outer cylinder is driven sinusoidally in time with zero mean rotation. This configuration was studied experimentally by Walsh and Donnelly. The critical Reynolds numbers calculated from linear stability theory agree well with the experimental values, except at large modulation amplitudes and small frequencies. The theoretical values are obtained using Floquet theory implemented in two distinct approaches: a truncated Fourier series representation in time, and a fundamental solution matrix based on a Chebyshev pseudospectral representation in space. For large amplitude, low frequency modulation, the linear eigenfunctions are temporally complex, consisting of a quiescent region followed by rapid change in the perturbed flow velocities.

Murray, B. T.↗

Statistical scaling relationships and size effects in the strength and creep rupture of fibrous composites

In this presentation we discuss a new theoretical model and supporting experimental results for the strength and lifetime in creep rupture of unidirectional, carbon fiber/epoxy matrix composites at ambient conditions. First we review the 'standard' Weibull/power-law methodology that has been standard practice. Then we discuss features of a recent model which build on the statistical aspects of fiber strength, micromechanical aspects of stress transfer around fiber breaks, and time-dependent creep of the matrix. The model is applied to 'microcomposites' consisting of seven fibers in a matrix for which strength and creep-rupture data are available. The model yields Weibull distributions in an envelope format for both strength and lifetime. The respective shape, scale and power-law parameters depend on such parameters as the Weibull shape parameter for fiber strength, the exponent for matrix creep, the effective load transfer length (which grows in time due to matrix creep) and the critical cluster size for failed fibers. The experimental results are consistent with the theory, though time-dependent debonding appears to be part of the failure process.

Phoenix, S. Leigh↗

CFD Computation of Broadband Fan Interaction Noise

In this study, a 3-D, unsteady, Reynolds Averaged Navier Stokes CFD code coupled to an acoustic calculation is used to predict the contribution of the exit guide vanes to broadband fan noise. The configuration investigated is that corresponding to the NASA Source Diagnostic Test (SDT) 22-in fan rig. Then an acoustic model introduced by Nallasamy which is based on 2-D strip theory is used to compute the broadband rotor-stator interaction noise. One configuration from the SDT matrix is considered here: the fan speed correlating to approach, and outlet guide vane count designed for cut-off of the blade passage frequency. Thus, in the chosen configuration, there are 22 rotor blades and 54 stator blades. The stators are located 2.5 tip chords downstream of the rotor trailing edge. The RANS computations are used to obtain the spectra of the unsteady surface pressure on the exit guide vanes. This surface pressure is then integrated together with the Green's function for and infinite cylindrical duct to obtain the acoustic field. The results from this investigation validate the use of the CFD code along with the acoustic model for broadband fan noise predictions. The validation enables future investigations such as the determination of rotor tip clearance and stator solidity effects on fan rotor-stator interaction noise.

Grace, Sheryl M.↗

Theoretical models for polarimetric radar clutter

The Mueller matrix and polarization covariance matrix are described for polarimetric radar systems. The clutter is modeled by a layer of random permittivity, described by a three-dimensional correlation function, with variance, and horizontal and vertical correlation lengths. This model is applied, using the wave theory with Born approximations carried to the second order, to find the backscattering elements of the polarimetric matrices. It is found that 8 out of 16 elements of the Mueller matrix are identically zero, corresponding to a covariance matrix with four zero elements. Theoretical predictions are matched with experimental data for vegetation fields.

Borgeaud, M.↗

Soft particles and infinite-dimensional geometry

Abstract In the sigma model, soft insertions of moduli scalars enact parallel transport ofS-matrix elements about the finite-dimensional moduli space of vacua, and the antisymmetric double-soft theorem calculates the curvature of the vacuum manifold. We explore the analogs of these statements in gauge theory and gravity in asymptotically flat spacetimes, where the relevant moduli spaces are infinite-dimensional. These models have spaces of vacua parameterized by (trivial) flat connections on the celestial sphere, and soft insertions of photons, gluons, and gravitons parallel transportS-matrix elements about these infinite-dimensional manifolds. We argue that the antisymmetric double-soft gluon theorem ind + 2 bulk dimensions computes the curvature of a connection on the infinite-dimensional space Map ( S d , G ) / G , whereGis the global part of the gauge group. The analogous metrics in abelian gauge theory and gravity are flat, as indicated by the vanishing of the antisymmetric double-soft theorems in those models. In other words, Feynman diagram calculations not only know about the vacuum manifold of Yang–Mills theory, they can also be used to compute its curvature. The results have interesting implications for flat space holography.

Astronomy & Astrophysics↗

Role of three-body dynamics in nucleon-deuteron correlation functions

Correlation functions of hadrons can be accessed in high-energy collisions of atomic nuclei, revealing information about the underlying interaction. This work complements experimental efforts to study nucleon-deuteron Nd—with N=p (proton) or N=n (neutron)—correlations with theory evaluations using different techniques. The correlation functions C nd and C pd are calculated based on a scattering wave function, extending previous benchmarks for the Nd scattering matrix to this new observable. We use hyperspherical harmonics and Faddeev techniques with one of the widely used nucleon-nucleon (NN) interactions, the Argonne v18 potential. Moreover, in the low-energy region we perform additional calculations in the framework of pionless effective field theory. The pd correlation function is computed in the large-energy region to make contact with a recent measurement by the ALICE Collaboration. We show that the scattering wave function has the proper dynamical input to describe an initial rise and subsequent oscillations of C pd as a function of the energy. Effects on the observables using different NN and three-nucleon potentials are evaluated with the conclusion that variations of around 2% are observed. Although these effects are small, future measurements can go beyond this accuracy, allowing for new detailed studies of strong interaction in light nuclear systems. Furthermore, the present study supports the current efforts devoted to the measurement of correlation functions in systems dominated by the strong interactions, such as pd, ppp, Λd, and ppΛ.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A static quantum embedding scheme based on coupled cluster theory

Here, we develop a static quantum embedding scheme that utilizes different levels of approximations to coupled cluster (CC) theory for an active fragment region and its environment. To reduce the computational cost, we solve the local fragment problem using a high-level CC method and address the environment problem with a lower-level Møller–Plesset (MP) perturbative method. This embedding approach inherits many conceptual developments from the hybrid second-order Møller–Plesset (MP2) and CC works by Nooijen [J. Chem. Phys. 111, 10815 (1999)] and Bochevarov and Sherrill [J. Chem. Phys. 122, 234110 (2005)]. We go beyond those works here by primarily targeting a specific localized fragment of a molecule and also introducing an alternative mechanism to relax the environment within this framework. We will call this approach MP-CC. We demonstrate the effectiveness of MP-CC on several potential energy curves and a set of thermochemical reaction energies, using CC with singles and doubles as the fragment solver, and MP2-like treatments of the environment. The results are substantially improved by the inclusion of orbital relaxation in the environment. Using localized bonds as the active fragment, we also report results for N=N bond breaking in azomethane and for the central C–C bond torsion in butadiene. We find that when the fragment Hilbert space size remains fixed (e.g., when determined by an intrinsic atomic orbital approach), the method achieves comparable accuracy with both a small and a large basis set. Additionally, our results indicate that increasing the fragment Hilbert space size systematically enhances the accuracy of observables, approaching the precision of the full CC solver.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Competing ionization and dissociation: Extension of the energy-dependent frame transformation to the gerade symmetry of H 2

This article solves two major tasks that frequently arise in the theory of electron collisions with a target molecular cation. First, it extends the energy-dependent frame transformation (EDFT) treatment, which is needed to map fixed-nuclei electron-molecule scattering matrices into an energy-dependent laboratory-frame scattering matrix with vibrational channel indices. The EDFT mapping can now be carried out even when the target molecule possesses multiple low-energy potential curves, significantly transcending previous applications. Second, it implements a method to extract the rest of the full laboratory-frame scattering matrix, i.e., the columns and rows describing input and/or output dissociation channels. The treatment is benchmarked in this article against the essentially exact solution of a refined two-dimensional model of the singlet gerade Σ symmetry of H 2 . Our tests demonstrate that the theory accurately maps fixed-nuclei scattering information, of the type provided by existing electron-molecule computer codes, into a laboratory-frame scattering matrix that includes both ionization and dissociation. Furthermore, this treatment can provide a general framework applicable to a broad class of electron collision processes involving diatomic target ions, suitable for an accurate description of challenging processes such as dissociative recombination.

74 ATOMIC AND MOLECULAR PHYSICS↗

Spin alignment of vector mesons by glasma fields

We explain how spin alignment of vector mesons can be induced by background fields, such as electromagnetic fields or soft gluon fields. Our study is based on the quantum kinetic theory of spinning quarks and antiquarks and incorporates the relaxation of the dynamically generated spin polarization. The spin density matrix of vector mesons is obtained by quark coalescence via the Wigner function and kinetic equation. Our approach predicts a local spin correlation that is distinct from the nonlocal expressions previously obtained in phenomenological derivations. We estimate the magnitude of such local correlations in the glasma model of the preequilibrium phase of relativistic heavy ion collisions. It is found that the resulting spin alignment could be greatly enhanced and may be comparable to the experimental measurement in order of magnitude. We further propose new phenomenological scenarios to qualitatively explain the transverse-momentum and centrality dependence of spin alignment in a self-consistent framework.

79 ASTRONOMY AND ASTROPHYSICS↗

Constituent-Based Life Models Being Developed for SiC/SiC Composites

For the successful utilization of ceramic matrix composites (CMC) as hot-section components in advanced aeropropulsion engines, the CMC constituents will need to be tailored and optimized to meet all the critical property demands of each component. Under the High-Speed Research (HSR) and Advanced High-Temperature Engine Materials Technology (HITEMP) Programs, the NASA Glenn Research Center at Lewis Field initiated research to develop mechanistic models for key CMC thermostructural properties. These models would describe the effects of different constituent factors (composition, geometry, and volume fraction) and of potential application conditions (stress, time, temperature, and environment) on these properties. Particular focus was placed on both analytical and numerical modeling of state-of-the-art SiC/SiC composites where the primary load-bearing constituents are stoichiometric SiC fibers in a complex multiphase SiC matrix produced by chemical vapor infiltration and melt infiltration. Recent studies have resulted in computer-based numerical models for the elastic modulus, thermal expansion, and thermal conductivity properties of the SiC/SiC system. Additional studies have generated analytical and empirical models for the time dependence of composite rupture strength at temperatures above 2200 F (1200 C), where CMC's have an important thermostructural advantage over current nickel-based superalloys. These life models utilize thermal activation theory and fiber stress-rupture results measured at Glenn to generate Larson-Miller (L M) plots of fiber rupture strength versus q, a single time- and temperature-dependent parameter. Assuming a worse case in which the SiC matrix is cracked, rupture is then controlled by the time-dependent fracture characteristics of the fiber bundles bridging the matrix cracks. With this as the controlling mechanism, one can then use simple composite theory and the fiber L M plots to predict CMC rupture strength versus the q parameter. The dashed line shows the predicted rupture strength of a SiC/SiC composite that is reinforced by a state-of-the-art stoichiometric SiC fiber. For the q parameter, time is in hours and temperature in degrees kelvin. To generate these predictions, a two-dimensional 0 /90 composite with approximately 16 percent fiber in the applied stress direction and an air test environment were assumed. As such, it is possible to compare the model predictions against limited stress-rupture data for this CMC as shown by the data points. The good agreement confirms the rupture model at least for the selected CMC and test conditions. Thus for this particular SiC/SiC composite, one can estimate a 1000-hr rupture strength of approximately 12 ksi at 2400 F (T = 1588 K and q = 39 700 K). At lower CMC application stresses, the SiC matrix is typically uncracked, so both the fiber and matrix constituents share the composite load. In this case, CMC rupture is controlled by the constituent with the longest rupture time based on the creep rate of the composite. Measured Monkman-Grant plots of rupture time versus creep rate for the two SiC constituents have been used to develop CMC life models for this important application condition. NASA and DOD are currently using this information to establish application and material goals for more advanced CMC's that can be used at even higher temperatures.

DiCarlo, James A.↗

Progressive Failure Analysis of 3D Woven Composites via Multiscale Recursive Micromechanics

Multiscale progressive failure simulations have been performed for a 3D woven composite considering six length scales, spanning from the woven composite repeating unit cell to the matrix material containing voids. The Multiscale Recursive Micromechanics approach, which enables micromechanics models to call other micromechanics models (or themselves recursively) to consider finer and finer length scales, has been employed. The multiscale model uses both the generalized method of cells and Mori-Tanaka micromechanics theories and considers failure and damage at the lowest length scale. A simple subvolume elimination damage model, as well as the crack band model, are employed for the matrix and their predicted composite responses compared with and without a preexisting binder tow disbond. Results are also compared with experimental data for an AS4 carbon fiber/RTM6 epoxy matrix 3D orthogonal woven PMC, with good agreement in terms of global stiffness and global failure stress.

3D woven↗

Theory of biaxial graded-index optical fiber

A biaxial graded-index fiber with a homogeneous cladding is studied. Two methods, wave equation and matrix differential equation, of formulating the problem and their respective solutions are discussed. For the wave equation formulation of the problem it is shown that for the case of a diagonal permittivity tensor the longitudinal electric and magnetic fields satisfy a pair of coupled second-order differential equations. Also, a generalized dispersion relation is derived in terms of the solutions for the longitudinal electric and magnetic fields. For the case of a step-index fiber, either isotropic or uniaxial, these differential equations can be solved exactly in terms of Bessel functions. For the cases of an istropic graded-index and a uniaxial graded-index fiber, a solution using the Wentzel, Krammers and Brillouin (WKB) approximation technique is shown. Results for some particular permittivity profiles are presented. Also the WKB solutions is compared with the vector solution found by Kurtz and Streifer. For the matrix formulation it is shown that the tangential components of the electric and magnetic fields satisfy a system of four first-order differential equations which can be conveniently written in matrix form. For the special case of meridional modes, the system of equations splits into two systems of two equations. A general iterative technique, asymptotic partitioning of systems of equations, for solving systems of differential equations is presented. As a simple example, Bessel's differential equation is written in matrix form and is solved using this asymptotic technique. Low order solutions for particular examples of a biaxial and uniaxial graded-index fiber are presented. Finally numerical results obtained using the asymptotic technique are presented for particular examples of isotropic and uniaxial step-index fibers and isotropic, uniaxial and biaxial graded-index fibers.

Kawalko, Stephen F.↗