Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix completion”

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 145 records · Page 8

Study of design and control of remote manipulators. Part 2: Vibration considerations in manipulator design

An investigation is reported involving a mathematical procedure using 4 x 4 transformation matrices for analyzing the vibrations of flexible manipulators. Previous studies with the procedure are summarized and the method is extended to include flexible joints as well as links, and to account for the effects of various power transmission schemes. A systematic study of the allocation of structural material and the placement of components such as motors and gearboxes was undertaken using the analytical tools developed. As one step in this direction the variables which relate the vibration parameters of the arm to the task and environment of the arm were isolated and nondimensionalized. The 4 x 4 transformation matrices were also used to develop analytical expressions for the terms of the complete 6 x 6 compliance matrix for the case of two flexible links joined by a rotating joint, flexible about its axis of rotation.

Book, W. J.↗

Theoretical studies of molecular collisions

The following subject areas are covered: (1) total integral reactive cross sections and vibrationally resolved reaction probabilities for F + H2 = HF + H; (2) a theoretical study of inelastic O + N2 collisions; (3) body frame close coupling wave packet approach to gas phase atom-rigit rotor inelastic collisions; (4) wave packet study of gas phase atom-rigit motor scattering; (5) the application of optical potentials for reactive scattering; (6) time dependent, three dimensional body frame quantal wave packet treatment of the H + H2 exchange reaction; (7) a time dependent wave packet approach to atom-diatom reactive collision probabilities; (8) time dependent wave packet for the complete determination of s-matrix elements for reactive molecular collisions in three dimensions; (9) a comparison of three time dependent wave packet methods for calculating electron-atom elastic scattering cross sections; and (10) a numerically exact full wave packet approach to molecule-surface scattering.

Kouri, Donald J.↗

Cratering and penetration experiments in teflon targets at velocities from 1 to 7 km/s

Approximately 20 sq m of protective thermal blankets, largely composed of Teflon, were retrieved from the Long Duration Exposure Facility after the spacecraft spent approximately 5.7 years in space. Examination of these blankets revealed that they contained thousands of hypervelocity impact features ranging from micron-sized craters to penetration holes several millimeters in diameter. We conducted impact experiments to reproduce such features and to understand the relationships between projectile size and the resulting crater or penetration hole diameter over a wide range of impact velocities. Such relationships are needed to derive the size and mass frequency distribution and flux of natural and man-made particles in low-earth orbit. Powder propellant and light-gas guns were used to launch soda-lime glass spheres into pure Teflon targets at velocities ranging from 1 to 7 km/s. Target thickness varied over more than three orders of magnitude from finite halfspace targets to very thin films. Cratering and penetration of massive Teflon targets is dominated by brittle failure and the development of extensive spall zones at the target's front and, if penetrated, the target's rear side. Mass removal by spallation at the back side of Teflon targets may be so severe that the absolute penetration hole diameter can become larger than that of a standard crater. The crater diameter in infinite halfspace Teflon targets increases, at otherwise constant impact conditions, with encounter velocity by a factor of V (exp 0.44). In contrast, the penetration hole size in very thin foils is essentially unaffected by impact velocity. Penetrations at target thicknesses intermediate to these extremes will scale with variable exponents of V. Our experimental matrix is sufficiently systematic and complete, up to 7 km/s, to make reasonable recommendations for velocity-scaling of Teflon craters and penetrations. We specifically suggest that cratering behavior and associated equations apply to all impacts in which the shock-pulse duration of the projectile is shorter than that assigned a unique projectile size, provided an impact velocity is known or assumed. This calibration seems superior to the traditional ballistic-limit approach.

Horz, Friedrich↗

Cratering and penetration experiments in Teflon targets at velocities from 1 to 7 km/s

Approximately 20 sq m of protective thermal blankets, largely composed of Teflon, were retrieved from the Long Duration Exposure Facility (LDEF) after the spacecraft had spent approximately 5.7 years in space. Examination of these blankets revealed that they contained thousands of hypervelocity impact features ranging from micron-sized craters to penetration holes several millimeters in diameter. We conducted impact experiments in an effort to reproduce such features and to -- hopefully -- understand the relationships between projectile size and the resulting crater or penetration-hole diameter over a wide range of impact velocity. Such relationships are needed to derive the size- and mass-frequency distribution and flux of natural and man-made particles in low-Earth orbit. Powder propellant and light-gas guns were used to launch soda-lime glass spheres of 3.175 mm (1/8 inch) nominal diameter (Dp) into pure Teflon FEP targets at velocities ranging from 1 to 7 km/s. Target thickness (T) was varied over more than three orders of magnitude from infinite halfspace targets (Dp/T less than 0.1) to very thin films (Dp/T greater than 100). Cratering and penetration of massive Teflon targets is dominated by brittle failure and the development of extensive spall zones at the target's front and, if penetrated, the target's rear side. Mass removal by spallation at the back side of Teflon targets may be so severe that the absolute penetration-hole diameter (Dh) can become larger than that of a standard crater (Dc) at relative target thicknesses of Dp/T = 0.6-0.9. The crater diameter is infinite halfspace Teflon targets increases -- at otherwise constant impact conditions -- with encounter velocity by a factor of V0.44. In contrast, the penetration-hole size is very thin foils (Dp/T greater than 50) is essentially unaffected by impact velocity. Penetrations at target thicknesses intermediate to these extremes will scale with variable exponents of V. Our experimental matrix is sufficiently systematic and complete, up to 7 km/s, to make reasonable recommendations for the velocity-scaling of Teflon craters and penetrations. We specifically suggest that cratering behavior and associated equations dominate all impacts in which the shock-pulse duration of the projectile (tp) is shorter than that of the target (tt). We also demonstrate that each penetration hole from space-retrieved surfaces may be assigned a unique projectile size, provided an impact velocity is known or assumed. This calibration seems superior to the traditional ballistic-limit approach.

Hoerz, Friedrich↗

Radioprotection of Human Cell Nuclear DNA by Polyamines: Radiosensitivity of Chromatin is Influenced by Tightly Bound Spermine

The polyamines putrescine (PUT) and spermine (SPM) were examined for their ability to protect human cell Deoxyribonucleic Acid (DNA) against the formation of radiation-induced double-strand breaks (DSBs). As observed previously, under conditions where polyamines were shown to be almost completely absent, association with nuclear matrix protein into a nucleoid, and organization into chromatin structure, protected DNA from induction of DSBs by factors of 4.5 and 95, respectively. At concentrations below 1 mM, PUT or SPM provided equivalent levels of protection to deproteinized nuclear DNA, consistent with their capacity to scavenge radiation-induced radicals. At constant ionic strength, 5 mM SPM protected deproteinized DNA and nucleoid DNA and DNA in nuclear chromatin by factors of 100 and 26, respectively. At 5 mM, SPM provided 15 times greater protection of deproteinized DNA than did PUT. Under physiologically relevant conditions, 5 mM SPM protected DNA in the intact nucleus from the induction of DSBs by a factor of 2 relative to DNA in the absence of SPM. Studies of SPM binding during cellular fractionation revealed that a significant fraction of the cellular SPM is tightly bound in the nucleus but can be removed by extended washing. Thus the association of SPM with nuclear chromatin appears to be a significant contributor to the resistance of the cell's DNA to the induction of DSBs.

Warters, Raymond L.↗

Influence of Finite Element Size in Residual Strength Prediction of Composite Structures

The sensitivity of failure load to the element size used in a progressive failure analysis (PFA) of carbon composite center notched laminates is evaluated. The sensitivity study employs a PFA methodology previously developed by the authors consisting of Hashin-Rotem intra-laminar fiber and matrix failure criteria and a complete stress degradation scheme for damage simulation. The approach is implemented with a user defined subroutine in the ABAQUS/Explicit finite element package. The effect of element size near the notch tips on residual strength predictions was assessed for a brittle failure mode with a parametric study that included three laminates of varying material system, thickness and stacking sequence. The study resulted in the selection of an element size of 0.09 in. X 0.09 in., which was later used for predicting crack paths and failure loads in sandwich panels and monolithic laminated panels. Comparison of predicted crack paths and failure loads for these panels agreed well with experimental observations. Additionally, the element size vs. normalized failure load relationship, determined in the parametric study, was used to evaluate strength-scaling factors for three different element sizes. The failure loads predicted with all three element sizes provided converged failure loads with respect to that corresponding with the 0.09 in. X 0.09 in. element size. Though preliminary in nature, the strength-scaling concept has the potential to greatly reduce the computational time required for PFA and can enable the analysis of large scale structural components where failure is dominated by fiber failure in tension.

Satyanarayana, Arunkumar↗

Updates to Implicit Edge-Based Gradient Methods

In this paper, we report updates to the implicit edge-based gradient methods originally introduced in [H. Nishikawa, AIAA Paper 2020-3048, 2020]. First, we clarify the relationship between gradient accuracy and truncation error and show that the quadratic method involves a free parameter. Then, we provide a complete description and a simplified matrix form of the implicit gradient systems including a consistent boundary treatment, and derive a set of parameters for achieving fourth-order gradient accuracy on regular tetrahedral grids. A stability analysis is performed for a relaxation scheme used to solve the implicit gradient systems, and the result serves as a guide for choosing parameters. Numerical results are shown for accuracy verification and also for realistic inviscid flow problems in three dimensions, including flows with shock waves.

Weighted Least-Squares↗

Alternate Inversion Paradigm for Spectroheliogram Data

Over the past five years, new methods to reconstruct spectrally pure maps of the Sun from spectroheliogram data have emerged, essentially unlocking this long-abandoned method of obtaining both spatial and spectral information over a large field of view simultaneously. The original inversion method determined the plasma’s emission measure distribution as a function of temperature at every spatial location in the field of view. To complete this inversion, a response matrix had to be created mapping the emission measure in different (temperature, space) bins to detector, requiring assumptions on the thermal and ionization equilibrium and abundance state of the plasma. We have since derived a new method of the inversion that does not require these atomic assumptions to be made. Instead, we use only the different locations of spectral lines from the same ion species and the possible ratios of those single-species spectral lines, removing the need for a priori knowledge on the state of the emitting plasma. In this presentation, we demonstrate this method using observed data from the Marshall Grazing Incidence X-ray Spectrometer (MaGIXS) and simulated data from the EUV CME and Coronal Connectivity Observatory (ECCCO) investigation.

Amy Winebarger↗

Library reuse in a rapid development environment

The Aeroscience and Flight Mechanics Division (AFMD) established a Rapid Development Laboratory (RDL) to investigate and improve new 'rapid development' software production processes and refine the use of commercial, off-the-shelf (COTS) tools. These tools and processes take an avionics design project from initial inception through high fidelity, real-time, hardware-in-the-loop (HIL) testing. One central theme of a rapid development process is the use and integration of a variety of COTS tools: This paper discusses the RDL MATRIX(sub x)(R) libraries, as well as the techniques for managing and documenting these libraries. This paper also shows the methods used for building simulations with the Advanced Simulation Development System (ASDS) libraries, and provides metrics to illustrate the amount of reuse for five complete simulations. Combining ASDS libraries with MATRIX(sub x)(R) libraries is discussed.

Uhde, JO↗

User and technical documentation

The program LIBRATE calculates velocities for trajectories from low earth orbit (LEO) to four of the five libration points (L2, L3, L4, and L5), and from low lunar orbit (LLO) to libration points L1 and L2. The flight to be analyzed departs from a circular orbit of any altitude and inclination about the Earth or Moon and finishes in a circular orbit about the Earth at the desired libration point within a specified flight time. This program produces a matrix of the delta V's needed to complete the desired flight. The user specifies the departure orbit, and the maximum flight time. A matrix is then developed with 10 inclinations, ranging from 0 to 90 degrees, forming the columns, and 19 possible flight times, ranging from the flight time (input) to 36 hours less than the input value, in decrements of 2 hours, forming the rows. This matrix is presented in three different reports including the total delta V's, and both of the delta V components discussed. The input required from the user to define the flight is discussed. The contents of the three reports that are produced as outputs are also described. The instructions are also included which are needed to execute the program.

Source record↗

Quantum-classical hybrid algorithm for the simulation of all-electron correlation

While chemical systems containing hundreds to thousands of electrons remain beyond the reach of quantum devices, hybrid quantum-classical algorithms present a promising pathway toward a quantum advantage. Hybrid algorithms treat the exponentially scaling part of the calculation-the static correlation-on the quantum computer and the non-exponentially scaling part-the dynamic correlation-on the classical computer. While a variety of algorithms have been proposed, the dependence of many methods on the total wave function limits the development of easy-to-use classical post-processing implementations. Here, we present a novel combination of quantum and classical algorithms, which computes the all-electron energy of a strongly correlated molecular system on the classical computer from the 2-electron reduced density matrix (2-RDM) evaluated on the quantum device. Significantly, we circumvent the wave function in the all-electron calculations by using density matrix methods that only require input of the statically correlated 2-RDM. Although the algorithm is completely general, we test it with two classical density matrix methods, the anti-Hermitian contracted Schrödinger equation (ACSE) and multiconfiguration pair-density functional theories, using the recently developed quantum ACSE method for simulating the statically correlated 2-RDM. Furthermore, we obtain experimental accuracy for the relative energies of all three benzyne isomers and thereby demonstrate the ability of the developed algorithm to achieve chemically relevant and accurate results on noisy intermediate-scale quantum devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Filled Elastomers: Mechanistic and Physics-Driven Modeling and Applications as Smart Materials

Elastomers are made of chain-like molecules to form networks that can sustain large deformation. Rubbers are thermosetting elastomers that are obtained from irreversible curing reactions. Curing reactions create permanent bonds between the molecular chains. On the other hand, thermoplastic elastomers do not need curing reactions. Incorporation of appropriated filler particles, as has been practiced for decades, can significantly enhance mechanical properties of elastomers. However, there are fundamental questions about polymer matrix composites (PMCs) that still elude complete understanding. This is because the macroscopic properties of PMCs depend not only on the overall volume fraction (ϕ) of the filler particles, but also on their spatial distribution (i.e., primary, secondary, and tertiary structure). This work aims at reviewing how the mechanical properties of PMCs are related to the microstructure of filler particles and to the interaction between filler particles and polymer matrices. Overall, soft rubbery matrices dictate the elasticity/hyperelasticity of the PMCs while the reinforcement involves polymer–particle interactions that can significantly influence the mechanical properties of the polymer matrix interface. For ϕ values higher than a threshold, percolation of the filler particles can lead to significant reinforcement. While viscoelastic behavior may be attributed to the soft rubbery component, inelastic behaviors like the Mullins and Payne effects are highly correlated to the microstructures of the polymer matrix and the filler particles, as well as that of the polymer–particle interface. Additionally, the incorporation of specific filler particles within intelligently designed polymer systems has been shown to yield a variety of functional and responsive materials, commonly termed smart materials. We review three types of smart PMCs, i.e., magnetoelastic (M-), shape-memory (SM-), and self-healing (SH-) PMCs, and discuss the constitutive models for these smart materials.

36 MATERIALS SCIENCE↗

Local stress within a granular molecular solvent matrix, a mechanism for individual ion hydration

The activity of water is important in many natural systems and directed processes. To better understand water activity, the hydration of individual ions is approximated for a recently developed speciation-based solution model. Here a matrix of values for hydration of complete salts (anion and cation) was used as source data to calculate individual ion hydration values. Some electrolytes were excluded from calculations of individual ion hydration due to the prevalence of ion pairing and resulting liberation of water upon ion pairing. Evaluation of the periodic trends emerging from individual ion hydration values suggests that ionic volume, relative to the volume of molecular water, dictates both modeled hydration requirements (related to vapor-liquid-equilibrium-derived activity) as well as properties historically described by hydrodynamic radius (ion diffusion and solution viscosity). When water is conceptualized as a granular three-dimensional matrix, solutes may function as substitutional defects (dilation or constriction) inducing local compressive or tensile stress. These results have implications on the mechanism of the Hofmeister series and series reversals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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

The goal of this research is to develop the transfer matrix method to treat nonlinear autonomous boundary value problems with multiple branches. The application is the complete nonlinear aeroelastic analysis of multiple-branched rotor blades. Once the development is complete, it can be incorporated into the existing transfer matrix analyses. There are several difficulties to be overcome in reaching this objective. The conventional transfer matrix method is limited in that it is applicable only to linear branch chain-like structures, but consideration of multiple branch modeling is important for bearingless rotors. Also, hingeless and bearingless rotor blade dynamic characteristics (particularly their aeroelasticity problems) are inherently nonlinear. The nonlinear equations of motion and the multiple-branched boundary value problem are treated together using a direct transfer matrix method. First, the formulation is applied to a nonlinear single-branch blade to validate the nonlinear portion of the formulation. The nonlinear system of equations is iteratively solved using a form of Newton-Raphson iteration scheme developed for differential equations of continuous systems. The formulation is then applied to determine the nonlinear steady state trim and aeroelastic stability of a rotor blade in hover with two branches at the root. A comprehensive computer program is developed and is used to obtain numerical results for the (1) free vibration, (2) nonlinearly deformed steady state, (3) free vibration about the nonlinearly deformed steady state, and (4) aeroelastic stability tasks. The numerical results obtained by the present method agree with results from other methods.

Murthy, V. R.↗

Seamless metal-clad fiber-reinforced organic matrix composite structures and process for their manufacture

A metallic outer sleeve is provided which is capable of enveloping a hollow metallic inner member having continuous reinforcing fibers attached to the distal end thereof. The inner member is then introduced into outer sleeve until inner member is completely enveloped by outer sleeve. A liquid matrix member is then injected into space between inner member and outer sleeve. A pressurized heat transfer medium is flowed through the inside of inner member, thereby forming a fiber reinforced matrix composite material. The wall thicknesses of both inner member and outer sleeve are then reduced to the appropriate size by chemical etching, to adjust the thermal expansion coefficient of the metal-clad composite structure to the desired value. thereby forming a fiber reinforced matrix composite material. The wall thicknesses of both inner member and outer sleeve are then reduced to the appropriate size by chemical etching, to adjust the thermal expansion coefficient of the metal-clad composite structure to the desired value. The novelty of this invention resides in the development of a efficient method of producing seamless metal clad fiber reinforced organic matrix composite structures.

Bluck, Raymond M.↗

Mechanical Properties of Continuous Fiber Reinforced Zirconium Diboride Matrix Composites

Continuous fiber reinforced zirconium diboride matrix composites, SCS-9a-(RBSiCZrB2)matrix, are being developed for leading edge, rocket nozzle and turbine engine applications. Recently, the composite materials have been characterized for tensile properties to 1250 C, the highest temperature tested. The tensile properties are fiber dominated as the matrix is microcracked on fabrication, but favorable failure characteristic are observed. Compression and shear mechanical testing results will be reported if completed. The effects of fiber volume fraction and matrix density on mechanical properties will be discussed. The target applications of the materials will be discussed. Specific testing being performed towards qualification for these applications will be included.

Stuffle, Kevin↗

Convection equation modeling: A non-iterative direct matrix solution algorithm for use with SINDA

The determination of the boundary conditions for a component-level analysis, applying discrete finite element and finite difference modeling techniques often requires an analysis of complex coupled phenomenon that cannot be described algebraically. For example, an analysis of the temperature field of a coldplate surface with an integral fluid loop requires a solution to the parabolic heat equation and also requires the boundary conditions that describe the local fluid temperature. However, the local fluid temperature is described by a convection equation that can only be solved with the knowledge of the locally-coupled coldplate temperatures. Generally speaking, it is not computationally efficient, and sometimes, not even possible to perform a direct, coupled phenomenon analysis of the component-level and boundary condition models within a single analysis code. An alternative is to perform a disjoint analysis, but transmit the necessary information between models during the simulation to provide an indirect coupling. For this approach to be effective, the component-level model retains full detail while the boundary condition model is simplified to provide a fast, first-order prediction of the phenomenon in question. Specifically for the present study, the coldplate structure is analyzed with a discrete, numerical model (SINDA) while the fluid loop convection equation is analyzed with a discrete, analytical model (direct matrix solution). This indirect coupling allows a satisfactory prediction of the boundary condition, while not subjugating the overall computational efficiency of the component-level analysis. In the present study a discussion of the complete analysis of the derivation and direct matrix solution algorithm of the convection equation is presented. Discretization is analyzed and discussed to extend of solution accuracy, stability and computation speed. Case studies considering a pulsed and harmonic inlet disturbance to the fluid loop are analyzed to assist in the discussion of numerical dissipation and accuracy. In addition, the issues of code melding or integration with standard class solvers such as SINDA are discussed to advise the user of the potential problems to be encountered.

Schrage, Dean S.↗