Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix approximation”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

A Piloted Simulator Evaluation of Transport Aircraft Rudder Pedal Force/Feel Characteristics

A piloted simulation study has been conducted in a fixed-base research simulator to assess the directional handling qualities for various rudder pedal feel characteristics for commercial transport airplanes. That is, the effects of static pedal force at maximum pedal travel, breakout force, and maximum pedal travel on handling qualities were studied. An artificial maneuver with a severe lateral wind shear and requiring runway tracking at an altitude of 50 feet in a crosswind was used to fully exercise the rudder pedals. Twelve active airline pilots voluntarily participated in the study and flew approximately 500 maneuvers. The pilots rated the maneuver performance with various rudder pedal feel characteristics using the Cooper- Harper rating scale. The test matrix had 15 unique combinations of the 3 static pedal feel characteristics. A 10-term, second-order equation for the Cooper-Harper pilot rating as a function of the 3 independent pedal feel parameters was fit to the data. The test matrix utilized a Central Composite Design that is very efficient for fitting an equation of this form. The equation was used to produce contour plots of constant pilot ratings as a function of two of the parameters with the third parameter held constant. These contour plots showed regions of good handling qualities as well as regions of degraded handling qualities. In addition, a numerical equation solver was used to predict the optimum parameter values (those with the lowest pilot rating). Quantitative pilot performance data were also analyzed. This analysis found that the peak values of the cross power spectra of the pedal force and heading angle could be used to quantify the tendency toward directional pilot induced oscillations (PIO). Larger peak values of the cross power spectra were correlated with larger (degraded) Cooper-Harper pilot ratings. Thus, the subjective data (Cooper-Harper pilot ratings) were consistent with the objective data (peak values of the cross power spectra).

Stewart, Eric C.↗
NASA NTRS2011

A Multiple Sphere T-Matrix Fortran Code for Use on Parallel Computer Clusters

A general-purpose Fortran-90 code for calculation of the electromagnetic scattering and absorption properties of multiple sphere clusters is described. The code can calculate the efficiency factors and scattering matrix elements of the cluster for either fixed or random orientation with respect to the incident beam and for plane wave or localized- approximation Gaussian incident fields. In addition, the code can calculate maps of the electric field both interior and exterior to the spheres.The code is written with message passing interface instructions to enable the use on distributed memory compute clusters, and for such platforms the code can make feasible the calculation of absorption, scattering, and general EM characteristics of systems containing several thousand spheres.

Mackowski, D. W.↗
NASA NTRS1990

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.↗
NASA NTRS2019

Materials and Processes for New Propulsion Systems with Reduced Environmental Impact

Aeronautics research at NASA Glenn Research Center includes development, characterization and modeling of high temperature, lightweight materials and fabrication processes for aircraft propulsion systems with increased efficiency and reduced emissions, fuel burn and noise. Current propulsion materials research includes Ceramic Matrix Composites and Environmental Barrier Coatings, Polymer Matrix Composites and Additive Manufacturing processes. This presentation will summarize recent progress and plans in these areas. Ceramic Matrix Composites for Turbine Components: As part of NASA's Aeronautics research, Glenn Research Center has developed Ceramic Matrix Composites for 2700 degrees Fahrenheit turbine engine applications in the next generation of ultra-efficient aircraft. In this presentation, the development of advanced fiber and matrix constituents that enabled this advancement will be reviewed, and characterization of the resulting improvements in mechanical properties and durability will be summarized. Progress toward the development and validation of models predicting the effects of the engine environment on durability of Ceramic Matrix Composites and Environmental Barrier Coatings will be summarized. Progress and plans for collaborative research with industry and other government agencies will be reviewed. Polymer Matrix Composites for Powertrain and Impact Protection: A lightweight, hybrid polymer matrix composite/steel gear concept for rotorcraft power transmission applications was evaluated at realistic speed, torque, and power conditions. Dynamic testing was performed at 5000 horsepower and 5400 revolutions per minute with no performance degradation. Test results demonstrated the potential of hybrid composite / metal gears to reduce gear weight by approximately 15 percent, enabling the implementation of multi-speed drive systems, which would otherwise have a weight penalty, for increased speed and efficiency. Results of post-test inspection and endurance tests will be presented. Separately, ballistic impact tests demonstrated improved damage tolerance in polymer matrix composites by incorporating lightweight, thermoplastic veil materials between selected composite plies during the composite fabrication process. Results of impact tests and post-test inspection will be presented. Additive Manufacturing for Electric Propulsion: Additive Manufacturing processes offer the potential to fabricate new, high power density electric motor designs with complex geometries, multi-material components and optimally designed components that would not be feasible with traditional manufacturing processes. High performance stator components using advanced 3-phase conductive coils were fabricated using direct printing with optimized silver pastes. Additional additive manufacturing methods are being used to fabricate rotors and motor housing components. Systems studies show that the higher performance and lower weight motors offer improved energy efficiency and reduced emissions. Improvements to motor performance resulting from optimized materials and component designs will be presented.

Additive Manufacturing↗
DOE OSTI2024

A high-order Shifted Interface Method for Lagrangian shock hydrodynamics

Here, we present a new method for two-material Lagrangian hydrodynamics, which combines the Shifted Interface Method (SIM) with a high-order Finite Element Method. Our approach relies on an exact (or sharp) material interface representation, that is, it uses the precise location of the material interface. The interface is represented by the zero level-set of a continuous high-order finite element function that moves with the material velocity. This strategy allows to evolve curved material interfaces inside curved elements. By reformulating the original interface problem over a surrogate (approximate) interface, located in proximity of the true interface, the SIM avoids cut cells and the associated problematic issues regarding implementation, numerical stability, and matrix conditioning. Accuracy is maintained by modifying the original interface conditions using Taylor expansions. We demonstrate the performance of the proposed algorithms on established numerical benchmarks in one, two and three dimensions.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1972

Power series evaluation of transition and covariance matrices.

Reexamination power series solutions to the matrix covariance differential equation and the transition differential equation. Truncation error bounds are derived which are computationally attractive and which extend previous results. Polynomial approximations are obtained by exploiting the functional equations satisfied by the transition and covariance matrices. The series-functional equation propagation technique represents a fast and accurate alternative to the numerical integration of the time-invariant transition and covariance equations.

Bierman, G. J.↗
NASA NTRS1977

Comment on the relation between the nonadiabatic coupling and the complex intersection of potential energy curves

Simple relations are discussed that provide a correspondence between the complex intersection of two potential surfaces and the nonadiabatic coupling matrix element between those surfaces. These are key quantities in semiclassical and quantum mechanical theories of collision induced electronic transitions. Within the two state approximation, the complex intersection is shown to be directly related to the location and magnitude of the peak in the nonadiabatic coupling. Two cases are considered: the avoided crossing between two potential surfaces; and the spin orbit interaction due to a P-2 halogen atom. Comparisons are made between the results of the two-state model and the results of ab initio quantum chemical calculations.

Jaffe, R. L.↗
NASA NTRS1978

Steady, Oscillatory and Unsteady, Subsonic and Supersonic Aerodynamics (SOUSSA) for complex aircraft configurations

The Green's function method and the computer program SOUSSA (Steady Oscillatory and Unsteady Subsonic and Supersonic Aerodynamics) are reviewed. The Green's function method is applied to the fully unsteady potential equation yielding an integro-differential-delay equation. This equation is approximated by a set of differential-delay equations in time using the finite element method. The Laplace transform is used to yield a matrix relating the velocity potential to the normal wash. The matrix of the generalized aerodynamic forces is obtained by premultiplying and postmultiplying the matrices relating generalized forces to the potential and the normal wash by the generalized coordinates. The program SOUSSA is compared with existing numerical results. Results indicate that the program is not only general, flexible, and easy to use, but also accurate and fast.

Morino, L.↗
NASA NTRS1977

Hamilton's principle - Finite-element methods and flexible body dynamics

A variational formulation is given for the equations of motion for an unconstrained elastic body, and Hamilton's principle is used to derive the equations of motion and deformation of the body. Finite-element approximations are developed for these dynamical equations with respect to a body axis system satisfying the mean axis conditions. The free-body influence matrix for the body then is developed in terms of the finite-element model parameters.

Cavin, R. K., III↗
NASA NTRS1991

Solution of three-dimensional afterbody flow using reduced Navier-Stokes equations

The flow over afterbody geometries was investigated using the reduced Navier-Stokes (RNS) approximation. Both pressure velocity flux-split and composites velocity primitive variable formulations were considered. Pressure or pseudopotential relaxation procedures are combined with sparse matrix or coupled strongly implicit algorithms to form a three-dimensional solver for general non-orthogonal coordinates. Three-dimensional subsonic and transonic viscous/inviscid interacting flows were evaluated. Solutions with and without regions of recirculation were obtained.

Almahroos, H. M. H.↗
NASA NTRS1992

Asteroid orbital error analysis: Theory and application

We present a rigorous Bayesian theory for asteroid orbital error estimation in which the probability density of the orbital elements is derived from the noise statistics of the observations. For Gaussian noise in a linearized approximation the probability density is also Gaussian, and the errors of the orbital elements at a given epoch are fully described by the covariance matrix. The law of error propagation can then be applied to calculate past and future positional uncertainty ellipsoids (Cappellari et al. 1976, Yeomans et al. 1987, Whipple et al. 1991). To our knowledge, this is the first time a Bayesian approach has been formulated for orbital element estimation. In contrast to the classical Fisherian school of statistics, the Bayesian school allows a priori information to be formally present in the final estimation. However, Bayesian estimation does give the same results as Fisherian estimation when no priori information is assumed (Lehtinen 1988, and reference therein).

Muinonen, K.↗
NASA NTRS1995

Explicit and implicit compact high-resolution shock-capturing methods for multidimensional Euler equations 1: Formulation

Two classes of explicit compact high-resolution shock-capturing methods for the multidimensional compressible Euler equations for fluid dynamics are constructed. Some of these schemes can be fourth-order accurate away from discontinuities. For the semi-discrete case their shock-capturing properties are of the total variation diminishing (TVD), total variation bounded (TVB), total variation diminishing in the mean (TVDM), essentially nonoscillatory (ENO), or positive type of scheme for 1-D scalar hyperbolic conservation laws and are positive schemes in more than one dimension. These fourth-order schemes require the same grid stencil as their second-order non-compact cousins. One class does not require the standard matrix inversion or a special numerical boundary condition treatment associated with typical compact schemes. Due to the construction, these schemes can be viewed as approximations to genuinely multidimensional schemes in the sense that they might produce less distortion in spherical type shocks and are more accurate in vortex type flows than schemes based purely on one-dimensional extensions. However, one class has a more desirable high-resolution shock-capturing property and a smaller operation count in 3-D than the other class. The extension of these schemes to coupled nonlinear systems can be accomplished using the Roe approximate Riemann solver, the generalized Steger and Warming flux-vector splitting or the van Leer type flux-vector splitting. Modification to existing high-resolution second- or third-order non-compact shock-capturing computer codes is minimal. High-resolution shock-capturing properties can also be achieved via a variant of the second-order Lax-Friedrichs numerical flux without the use of Riemann solvers for coupled nonlinear systems with comparable operations count to their classical shock-capturing counterparts. The simplest extension to viscous flows can be achieved by using the standard fourth-order compact or non-compact formula for the viscous terms.

Yee, H. C.↗
NASA NTRS2013

Highly Reduced Forsterite and Enstatite from Stardust Track 61: Implications for Radial Transport of E Asteroid Material

The Stardust Mission returned a large fraction of high-temperature, crystalline material that was radially transported from the inner solar system to the Kuiper Belt [1,2]. The mineralogical diversity found in this single cometary collection points to an even greater number of source materials than most primitive chondrites. In particular, the type II olivine found in Wild 2 includes the three distinct Fe/Mn ratios found in the matrix and chondrules of carbonaceous chondrites (CCs) and unequilibrated ordinary chondrites (UOCs) [3]. We also find that low-Ca pyroxene is quite variable (approximately Fs3-29) and is usually indistinguishable from CC, UOC, and EH3 pyroxene as well. However, occasional olivine and pyroxene compositions are found in Wild 2 that are inconsistent with chondrites. The Stardust track 61 terminal particle (TP) is one such example and is the focus of this study. It s highly reduced forsterite and enstatite is consistent only with that in Aubrites, in which FeO is essentially absent from these phases (less than approximately 0.1 wt.% FeO) [4].

Frank, David R.↗
NASA NTRS1985

Linear and nonlinear dynamic analysis of redundant load path bearingless rotor systems

The bearingless rotorcraft offers reduced weight, less complexity and superior flying qualities. Almost all the current industrial structural dynamic programs of conventional rotors which consist of single load path rotor blades employ the transfer matrix method to determine natural vibration characteristics because this method is ideally suited for one dimensional chain like structures. This method is extended to multiple load path rotor blades without resorting to an equivalent single load path approximation. Unlike the conventional blades, it isk necessary to introduce the axial-degree-of-freedom into the solution process to account for the differential axial displacements in the different load paths. With the present extension, the current rotor dynamic programs can be modified with relative ease to account for the multiple load paths without resorting to the equivalent single load path modeling. The results obtained by the transfer matrix method are validated by comparing with the finite element solutions. A differential stiffness matrix due to blade rotation is derived to facilitate the finite element solutions.

Murthy, V. R.↗
DOE OSTI2021

Analytic techniques for solving the transport equations in electroweak baryogenesis

We develop an efficient method for solving transport equations, particularly in the context of electroweak baryogenesis. It provides fully-analytical results under mild approximations and can also test semi-analytical results, which are applicable in more general cases. Key elements of our method include the reduction of the second-order differential equations to first order, representing the set of coupled equations as a block matrix of the particle densities and their derivatives, identification of zero modes, and block decomposition of the matrix. We apply our method to calculate the baryon asymmetry of the Universe (BAU) in a Standard Model effective field theory framework of complex Yukawa couplings to determine the sensitivity of the resulting BAU to modifications of various model parameters and rates, and to estimate the effect of the commonly-used thin-wall approximation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Entanglement perspective on the quantum approximate optimization algorithm

Many quantum algorithms seek to output a specific bitstring solving the problem of interest—or a few if the solution is degenerate. It is the case for the quantum approximate optimization algorithm (QAOA) in the limit of large circuit depth, which aims to solve quadratic unconstrained binary optimization problems. Hence, the expected final state for these algorithms is either a product state or a low-entangled superposition involving a few bitstrings. What happens in between the initial N -qubit product state | 0 〉 ⊗ N and the final one regarding entanglement? Here, we consider the QAOA algorithm for solving the paradigmatic MaxCut problem on different types of graphs. We study the entanglement growth and spread resulting from randomized and optimized QAOA circuits and find that there is a volume-law entanglement barrier between the initial and final states. We also investigate the entanglement spectrum in connection with random matrix theory. In addition, we compare the entanglement production with a quantum annealing protocol aiming to solve the same MaxCut problems. Finally, we discuss the implications of our results for the simulation of QAOA circuits with tensor network-based methods relying on low-entanglement for efficiency, such as matrix product states.

Dupont, Maxime↗
DOE OSTI2022

Vibrational Entropy of Crystalline Solids from Covariance of Atomic Displacements

The vibrational entropy of a solid at finite temperature is investigated from the perspective of information theory. Ab initio molecular dynamics (AIMD) simulations generate ensembles of atomic configurations at finite temperature from which we obtain the N-body distribution of atomic displacements, ρN. We calculate the information-theoretic entropy from the expectation value of lnρN. At a first level of approximation, treating individual atomic displacements independently, our method may be applied using Debye–Waller B-factors, allowing diffraction experiments to obtain an upper bound on the thermodynamic entropy. At the next level of approximation we correct the overestimation through inclusion of displacement covariances. We apply this approach to elemental body-centered cubic sodium and face-centered cubic aluminum, showing good agreement with experimental values above the Debye temperatures of the metals. Below the Debye temperatures, we extract an effective vibrational density of states from eigenvalues of the covariance matrix, and then evaluate the entropy quantum mechanically, again yielding good agreement with experiment down to low temperatures. Our method readily generalizes to complex solids, as we demonstrate for a high entropy alloy. Further, our method applies in cases where the quasiharmonic approximation fails, as we demonstrate by calculating the HCP/BCC transition in Ti.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Entanglement and many-body effects in collective neutrino oscillations

Collective neutrino oscillations play a crucial role in transporting lepton flavor in astrophysical settings, such as supernovae, where the neutrino density is large. In this regime, neutrino-neutrino interactions are important and simulations in the mean-field approximation show evidence for collective oscillations occurring at timescales much shorter than those associated with vacuum oscillations. In this work, we study the out-of-equilibrium dynamics of a corresponding spin model using matrix product states and show how collective bipolar oscillations can be triggered by many-body correlations if appropriate initial conditions are present. We find entanglement entropies scaling at most logarithmically in the system size suggesting that classical tensor network methods could be efficient in describing collective neutrino dynamics more generally. These observation provide a clear path forward, not only to increase the accuracy of current simulations, but also to elucidate the mechanism behind collective flavor oscillations without resorting to the mean-field approximation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗