Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix elements”

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,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

On the Singularity in the Estimation of the Quaternion-of-Rotation

It has been claimed in the archival literature that the covariance matrix of a Kalman filter, which is designed to estimate the quaternion-of-rotation, is necessarily rank deficient because the normality constraint of the quaternion produces dependence between the quaternion elements. In reality, though, this phenomenon does not occur. The covariance matrix is not singular, and the filter is well behaved. Several simple examples are presented that demonstrate the regularity of the covariance matrix. First, estimation cases are presented where a relationship exists between the estimated variables, and yet the covariance matrix is not singular. Then the particular problem of quaternion estimation is analyzed. It is shown that the discrepancy stems from the fact that a functional relationship exists between the elements of the true quaternion but not between its estimated elements.

Bar-Itzhack, Itzhack Y.↗
NASA NTRS1991

Analysis of delamination onset and growth in curved laminates

Delamination onset and growth from a matrix crack in a curved composite laminate are studied. Finite element analysis is used to determine values of G which can be compared to the interlaminar fracture toughness of the composite for delamination predictions. The delamination is assumed to grow from the matrix crack and to remain straight and perpendicular to the edge. G is larger at the free edge than in the center and increases with delamination growth. Thus, once delamination begins from the matrix crack it should continue around the curved section with no increase in applied load. In the local vicinity of the matrix crack the rate of change of G decreases initially and then begins to increase with a point of inflection in between. Measured failure loads and the interlaminar fracture toughness of the composite are used to show that the point of inflection corresponds to the critical value of G.

Martin, Roderick H.↗
DOE OSTI2022

Phenomenological R -Matrix parameterization of direct, doorway, and compound nuclear reactions [Abstract]

Although formal expressions for scattering matrix accounting for direct, doorway, and compound nuclear (CN) resonant reactions have been derived several decades ago in both the transition ( T -)matrix formalism and the reactance ( K -)matrix formalism, the absence of corresponding expressions in phenomenological R -matrix formalism has limited the application of the latter to CN resonant reactions only. We remove this limitation by parameterizing direct, doorway, and CN resonant reactions in a phenomenological R -matrix scattering matrix, and provide a parameterization for a corresponding Reich-Moore approximation of eliminated capture channels. Direct reactions induce (previously neglected) mixing among the incoming or outgoing R -matrix channel wave functions, parameterized by real and orthonormal channel-rotation matrix, M , whereby the original scattering matrix U is transformed into M T UM . Any real and orthonormal matrix, M , can be equivalently expressed as e η , where η is a real and skew-symmetric 2 rotation-generating matrix that subsequently yields a more intuitive parameterization of eliminated direct capture reactions in Reich-Moore approximation. A phenomenological R -matrix parameterization of doorway reactions is inferred by equating the expression for reactance ( K -)matrix, given in terms of Brune’s alternative R -matrix parameterization, to a corresponding expression derived using Feshbach’s projection operator formalism. Assuming that all doorway states, just like CN states, are confined within spheres defined by R -matrix channel radii, a new R -matrix-like term induced by doorway states is gleaned, wherein each doorway state is parameterized by its energy, width, and the strength of its coupling to each CN state. Since a Reich-Moore approximation for retained-channel scattering matrix ought to approximate the effect of eliminated capture channels taking place via direct, doorway, or CN reactions, each of the three kinds of reactions contributing to the capture entails a corresponding Reich-Moore parameterization in a first-order approximation: direct contribution is parameterized by introducing finite diagonal elements of a retained-channel rotation-generating matrix, doorway contribution is parameterized by doorway capture widths, while CN contribution is parameterized by conventional Reich-Moore capture widths. We will present evidence of direct and doorway reactions observed in recent measurements of resolved resonance cross sections at the Gaerttner LINAC Center at Rensselaer Polytechnic Institute, and will outline a path for implementing this new R -matrix parameterization into the SAMMY nuclear data evaluation code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Matrix-free preconditioning for high-order H (curl) discretizations

The greater arithmetic intensity of high-order finite element discretizations makes them attractive for implementation on next-generation hardware, but assembly of high-order finite element operators as matrices is prohibitively expensive. As a result, the development of general algebraic solvers for such operators has been an open research challenge. Fast matrix-free application of high-order operators has received significant attention in the literature in the context of Poisson-type problems, but preconditioners and solvers for inverting more general operators are not very well-developed. In this paper, we consider the problem of preconditioning a definite Maxwell operator at high polynomial order without assembling a matrix. We show that given efficient preconditioners for high-order H 1 finite element problems on the same mesh, efficient H(curl) preconditioners can be constructed in an auxiliary space framework. We demonstrate the resulting preconditioners in a practical setting with tensor-product basis functions on an unstructured mesh of quadrilaterals. Overall, our approach uses a sparsified H 1 solver constructed on a low-order mesh of the nodal points of the underlying high-order space, and we show that the resulting H(curl) preconditioner is effective at very high polynomial orders for two-dimensional model problems with complicated geometry, varying piecewise constant coefficients, and curved elements. The resulting preconditioner scales with nearly optimal O(p d+1 ) floating point operation count and optimal O(p d ) memory transfer requirements, outperforming existing Maxwell preconditioners in the high-order regime.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2003

Mineralogy and Textural Characteristics of Fine-grained Rims in the Yamato 791198 CM2 Carbonaceous Chondrite: Constraints on the Location of Aqueous Alteration

Carbonaceous chondrites provide important clues into the nature of physical and chemical processes in the early solar system. A question of key importance concerns the role of water in solar nebular and asteroidal processes. The effects of water on primary mineral assemblages have been widely recognized in chondritic meteorites, especially the CI and CM carbonaceous chondrites. These meteorites have undergone extensive aqueous alteration that occurred prior to their arrival on Earth. In the case of the CM chondrites, this alteration has resulted in the partial to complete replacement of the primary nebular phases with secondary alteration phases. Considerable controversy exists as to the exact location where the alteration of the CM chondrites occurred. Several textural lines of evidence have been cited in support of aqueous alteration prior to the accretion of the final parent asteroid. An important line of evidence to support this hypothesis is the dis-equilibrium nature of fine-grained rims and matrix materials. [2] also noted the juxtaposition of micron-sized Fe-Ni metal grains and apparently unaltered chondrule glass against hydrated rim silicates. Conversely, there is a large body of evidence in favor of parent body alteration such as the occurrence of undisturbed Fe-rich aureoles and the systematic redistribution of elemental components over millimeters, e.g., Mg(+2) into the matrix and Fe(+2) into chondrules etc.

Chizmadia, Lysa J.↗
NASA NTRS1984

Modeling material failure with a vectorized routine

The computational aspects of modelling material failure in structural wood members are presented with particular reference to vector processing aspects. Wood members are considered to be highly orthotropic, inhomogeneous, and discontinuous due to the complex microstructure of wood material and the presence of natural growth characteristics such as knots, cracks and cross grain in wood members. The simulation of strength behavior of wood members is accomplished through the use of a special purpose finite element/fracture mechanics routine, program STARW (Strength Analysis Routine for Wood). Program STARW employs quadratic finite elements combined with singular crack tip elements in a finite element mesh. Vector processing techniques are employed in mesh generation, stiffness matrix formation, simultaneous equation solution, and material failure calculations. The paper addresses these techniques along with the time and effort requirements needed to convert existing finite element code to a vectorized version. Comparisons in execution time between vectorized and nonvectorized routines are provided.

Cramer, S. M.↗
NASA NTRS2002

On the Singularity in the Estimation of the Quaternion-of-Rotation

It has been claimed in the archival literature that the covariance matrix of a Kalman filter, which is designed to estimate the quaternion-of-rotation, is necessarily rank, deficient because the normality constraint of the quaternion produces dependence between the quaternion elements. In reality, though, this phenomenon does not occur. The covariance matrix is not singular, and the filter is well behaved. Several simple examples are presented th at demonstrate the regularity of the covariance matrix. First, a Kalman filter is designed to estimate variables subject to a functional relationship. Then the particular problem of quaternion estimation is analyzed. It is shown that the discrepancy stems from the fact that the functional relationship exists between the elements of the quaternion but not between its estimate elements.

Bar-Itzhack, Itzhack Y.↗
NASA NTRS2016

Using Tabulated Experimental Data to Drive an Orthotropic Elasto-Plastic Three-Dimensional Model for Impact Analysis

An orthotropic elasto-plastic-damage three-dimensional model with tabulated input has been developed to analyze the impact response of composite materials. The theory has been implemented as MAT 213 into a tailored version of LS-DYNA being developed under a joint effort of the FAA and NASA and has the following features: (a) the theory addresses any composite architecture that can be experimentally characterized as an orthotropic material and includes rate and temperature sensitivities, (b) the formulation is applicable for solid as well as shell element implementations and utilizes input data in a tabulated form directly from processed experimental data, (c) deformation and damage mechanics are both accounted for within the material model, (d) failure criteria are established that are functions of strain and damage parameters, and mesh size dependence is included, and (e) the theory can be efficiently implemented into a commercial code for both sequential and parallel executions. The salient features of the theory as implemented in LS-DYNA are illustrated using a widely used composite - the T800S/3900-2B[P2352W-19] BMS8-276 Rev-H-Unitape fiber/resin unidirectional composite. First, the experimental tests to characterize the deformation, damage and failure parameters in the material behavior are discussed. Second, the MAT213 input model and implementation details are presented with particular attention given to procedures that have been incorporated to ensure that the yield surfaces in the rate and temperature dependent plasticity model are convex. Finally, the paper concludes with a validation test designed to test the stability, accuracy and efficiency of the implemented model.

Polymer Matrix Composites↗
NASA NTRS1989

Optical laboratory solution and error model simulation of a linear time-varying finite element equation

The use of simplified error models to accurately simulate and evaluate the performance of an optical linear-algebra processor is described. The optical architecture used to perform banded matrix-vector products is reviewed, along with a linear dynamic finite-element case study. The laboratory hardware and ac-modulation technique used are presented. The individual processor error-source models and their simulator implementation are detailed. Several significant simplifications are introduced to ease the computational requirements and complexity of the simulations. The error models are verified with a laboratory implementation of the processor, and are used to evaluate its potential performance.

Taylor, B. K.↗
NASA NTRS1996

Lunar Meteorite Queen Alexandra Range 93069 and the Iron Concentration of the Lunar Highlands Surface

Lunar meteorite Queen Alexandra Range 93069 is a clast-rich, glassy-matrix regolith breccia of ferroan, highly aluminous bulk composition. It is similar in composition to other feldspathic lunar meteorites but differs in having higher concentrations of siderophile elements and incompatible trace elements. Based on electron microprobe analyses of the fusion crust, glassy matrix, and clasts, and instrumental neutron activation analysis of breccia fragments, QUE 93069 is dominated by nonmare components of ferroan, noritic- anorthosite bulk composition. Thin section QUE 93069,31 also contains a large, impact-melted, partially devitrified clast of magnesian, anorthositic-norite composition. The enrichment in Fe, Sc, and Cr and lower Mg/Fe ratio of lunar meteorites Yamato 791197 and Yamato 82192/3 compared to other feldspathic lunar meteorites can be attributed to a small proportion (5-10%) of low-Ti mare basalt. It is likely that the non- mare components of Yamato 82192/3 are similar to and occur in similar abundance to those of Yamato 86032, with which it is paired. There is a significant difference between the average FeO concentration of the lunar highlands surface as inferred from the feldspathic lunar meteorites (mean: approx. 5.0%; range: 4.3-6.1 %) and a recent estimate based on data from the Clementine mission (3.6%).

Korotev, Randy L.↗
NASA NTRS1976

A comparison of matrix methods for calculating eigenvalues in acoustically lined ducts

Three approximate methods - finite differences, weighted residuals, and finite elements - were used to solve the eigenvalue problem which arises in finding the acoustic modes and propagation constants in an absorptively lined two-dimensional duct without airflow. The matrix equations derived for each of these methods were solved for the eigenvalues corresponding to various values of wall impedance. Two matrix orders, 20 x 20 and 40 x 40, were used. The cases considered included values of wall admittance for which exact eigenvalues were known and for which several nearly equal roots were present. Ten of the lower order eigenvalues obtained from the three approximate methods were compared with solutions calculated from the exact characteristic equation in order to make an assessment of the relative accuracy and reliability of the three methods. The best results were given by the finite element method using a cubic polynomial. Excellent accuracy was consistently obtained, even for nearly equal eigenvalues, by using a 20 x 20 order matrix.

Watson, W.↗
NASA NTRS1994

A new look at the simultaneous analysis and design of structures

The minimum weight optimization of structural systems, subject to strength and displacement constraints as well as size side constraints, was investigated by the Simultaneous ANalysis and Design (SAND) approach. As an optimizer, the code NPSOL was used which is based on a sequential quadratic programming (SQP) algorithm. The structures were modeled by the finite element method. The finite element related input to NPSOL was automatically generated from the input decks of such standard FEM/optimization codes as NASTRAN or ASTROS, with the stiffness matrices, at present, extracted from the FEM code ANALYZE. In order to avoid ill-conditioned matrices that can be encountered when the global stiffness equations are used as additional nonlinear equality constraints in the SAND approach (with the displacements as additional variables), the matrix displacement method was applied. In this approach, the element stiffness equations are used as constraints instead of the global stiffness equations, in conjunction with the nodal force equilibrium equations. This approach adds the element forces as variables to the system. Since, for complex structures and the associated large and very sparce matrices, the execution times of the optimization code became excessive due to the large number of required constraint gradient evaluations, the Kreisselmeier-Steinhauser function approach was used to decrease the computational effort by reducing the nonlinear equality constraint system to essentially a single combined constraint equation. As the linear equality and inequality constraints require much less computational effort to evaluate, they were kept in their previous form to limit the complexity of the KS function evaluation. To date, the standard three-bar, ten-bar, and 72-bar trusses have been tested. For the standard SAND approach, correct results were obtained for all three trusses although convergence became slower for the 72-bar truss. When the matrix displacement method was used, correct results were still obtained, but the execution times became excessive due to the large number of constraint gradient evaluations required. Using the KS function, the computational effort dropped, but the optimization seemed to become less robust. The investigation of this phenomenon is continuing. As an alternate approach, the code MINOS for the optimization of sparse matrices can be applied to the problem in lieu of the Kreisselmeier-Steinhauser function. This investigation is underway.

Striz, Alfred G.↗
NASA NTRS2011

A Bi-Modal Distribution of ALHA77307 Matrix Olivine: Evidence for Fine-Grained Mixing from Multiple Reservoirs in the CO Formation Zone

The CO 3.0 chondrite ALHA77307 is thought to be the least metamorphosed of all the CO chondrites [1]. As such, the fine-grained (<30 m) olivine found in its matrix is a valuable resource for investigating the CO formation environment since its compositions should be primary. In the CO matrix, we indeed find a wide range of major element compositions (Fa(0.5-71)). However, more importantly, we find that the olivines make up two compositionally distinct populations (Fa(0.5-5) and Fa(21-71)). Grains from both populations are found within an extremely close proximity and we see no obvious evidence of two distinct lithologies within our samples. Therefore, we conclude that the olivine grains found in the ALHA77307 matrix must have crystallized within two unique formation conditions and were later mixed at a very fine scale during the accretion epoch. Here, we propose a possible explanation based on Cr and Mn concentrations in the olivine.

Frank, D.↗
NASA NTRS2014

Tailorable Dielectric Material with Complex Permittivity Characteristics

A dielectric material includes a network of nanosubstrates, such as but not limited to nanotubes, nanosheets, or other nanomaterials or nanostructures, a polymer base material or matrix, and nanoparticles constructed at least partially of an elemental metal. The network has a predetermined nanosubstrate loading percentage by weight with respect to a total weight of the dielectric material, and a preferential or predetermined longitudinal alignment with respect to an orientation of an incident electrical field. A method of forming the dielectric material includes depositing the metal-based nanoparticles onto the nanosubstrates and subsequently mixing these with a polymer matrix. Once mixed, alignment can be achieved by melt extrusion or a similar mechanical shearing process. Alignment of the nanosubstrate may be in horizontal or vertical direction with respect to the orientation of an incident electrical field.

Dudley, Kenneth L.↗
NASA NTRS1976

On a finite dynamic element method for free vibration analysis of structures

This paper explores the concept of finite dynamic elements involving higher order dynamic correction terms in the associated stiffness and mass matrices. Such matrices are then developed for a rectangular prestressed membrane element. Next, efficient analysis techniques for the eigenproblem solution of the resulting quadratic matrix equations are described in detail. These are followed by suitable numerical examples which indicate that employment of such dynamic elements in conjunction with an efficient quadratic matric solution technique will result in a most significant economy in the free vibration analysis of structures.

Gupta, K. K.↗
DOE OSTI2024

Simulation of Multiphase Flow and Poromechanical Effects Around Injection Wells in CO 2 Storage Sites

In geological CO 2 storage operations, wellbore deformations and leakage pathways formations can occur around injection and abandoned wells subjected to high rates and long-term CO 2 injection. To guide engineering design and prevent CO 2 leakage risks, a full understanding of the underlying physics and robust numerical models is necessary to evaluate the response of underground formations in the near wellbore region and in the reservoir. In this study, a multi-scale and multi-physics open-source simulator (GEOS) is used to simulate multiphase flow and poromechanical deformations over time in three dimensions. The governing equations for mechanical deformations of the rock body and multiphase compositional fluid flow within the rock matrix are solved with a fully coupled finite element and finite volume approach. The Drucker–Prager model with friction hardening is applied to simulate elastoplastic deformation and a multiphase fluid model with power-law correlations for relative permeability is used to model the migration of CO 2 plume, which are coupled with numerical implicit scheme. Simulation results are verified against multiple analytical solutions for multiphase flow and wellbore problems, thus demonstrating the accuracy of this advanced simulator. In two engineering applications, here we highlight the impact of elastoplastic deformation and coupled modeling for assessing induced displacements and stress perturbations, which are more pronounced in the near wellbore regions. This work focuses on short-term processes in the vicinity of injection wells where stress evolutions, rock deformations and multiphase compositional flow and transport are simulated jointly to ensure wellbore stability and prevent damage. This fully coupled geomechanical model can simulate multiphase flow and any associated poromechanical effects within the CO 2 storage site and in the surrounding formations. Such a large-scale, long-term, multi-physics simulation model is useful in many ways: it can guide operational decisions for CO 2 injection, assess the containment potential and risks of a site, and analyze the wellbore stability and integrity during and after CO 2 injection.

58 GEOSCIENCES↗
DOE OSTI2021

Applications of Strain-Coupled Magnetoelectric Composites

This article deals with research, development and future directions of magnetoelectric composites for practical devices applications. In the past 20 years there has been a surge of research in the area of multiferroics (MF) and magnetoelectrics (ME) due to their potential to replace existing technologies based only on ferroelectric or ferromagnetic materials. Some of the magnetoelectric composites show exceptionally high potential in the area of magnetic field sensors, however, work remains before commercialization can be realized. The cross coupling among various ferroic parameters in magnetoelectric composites is several orders higher than single phase magetoelectrics, which make its favorable for low detection (nT or pT) magnetic field sensors. Advances in both layered structures or controlled three-dimensional matrix composites for applications as magnetoelectric nonvolatile memory elements are both required. Robust cross-coupling among various parameters with more than four logic states and its compatibility with complementary-symmetry metal–oxide–semiconductor (CMOS) technology are the main requirements for heterostructure magnetoelectric thin films. The major hurdles in the area of magnetoelectric nonvolatile memory elements are poor interfacial properties and weak magnetoelectric coupling for high density fast read and write processes. Another potential area is strained coupled magneto-electric composites where magnetostriction mediated dimensional change in the magnetic layer effectively modulate the change in the dimension of piezoelectric layers via piezostriction, which leads to a strong ME coupling where their coupling magnitude is sufficient for magnetic field sensors. Energy harvesters based on magnetoelectric composites are also intriguing concepts to capture various types of waste energy, in the form mechanical vibration, pressure, wind energy, hydrothermal and waste temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Characterization of γ′/γ″ compact and sandwich type precipitates during long-term high-temperature exposure in Ni-based superalloys

Here, the formation of γ′/γ″ co-precipitates is investigated in Ni-based superalloys with a varying Ti/Al ratio and Ta content and their stability is studied using long-term high-temperature exposure. Transmission electron microscopy and atom probe tomography analyses demonstrate that both higher Ti/Al ratios and increased Ta promote γ″ phase formation, leading to sandwich and compact structures. The compact co-precipitation significantly reduces γ′ precipitate coarsening during 10,000 h exposure at 700°C by restricting elemental diffusion, particularly of aluminum, from the γ matrix to the γ′ phase. For the alloy without a compact structure, and only γʹ precipitates, at the beginning of the exposure, the coarsening rate over 10,000 h was 3.5 times faster than for the alloy with compact γʹ/γʺ precipitates. Evidence of destabilization of the compact morphology was found to occur between 5,000 and 10,000 h exposure and originated from the extensive formation and growth of δ platelets that extended throughout the grains. Thus, the outer layer of the compact, which consists of γʺ, was subjected to the γ″ to δ phase transformation.

gamma double prime↗