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,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Progressive Failure Studies of Composite Panels with and without Cutouts

Progressive failure analyses results are presented for composite panels with and without a cutout and subjected to in-plane shear loading and compression loading well into their postbuckling regime. Ply damage modes such as matrix cracking, fiber-matrix shear, and fiber failure are modeled by degrading the material properties. Results from finite element analyses are compared with experimental data. Good agreement between experimental data and numerical results are observed for most structural configurations when initial geometric imperfections are appropriately modeled.

Jaunky, Navin↗
NASA NTRS2012

Modeling of Damage Initiation and Progression in a SiC/SiC Woven Ceramic Matrix Composite

The goal of an ongoing project at NASA Glenn is to investigate the effects of the complex microstructure of a woven ceramic matrix composite and its variability on the effective properties and the durability of the material. Detailed analysis of these complex microstructures may provide clues for the material scientists who `design the material? or to structural analysts and designers who `design with the material? regarding damage initiation and damage propagation. A model material system, specifically a five-harness satin weave architecture CVI SiC/SiC composite composed of Sylramic-iBN fibers and a SiC matrix, has been analyzed. Specimens of the material were serially sectioned and polished to capture the detailed images of fiber tows, matrix and porosity. Open source analysis tools were used to isolate various constituents and finite elements models were then generated from simplified models of those images. Detailed finite element analyses were performed that examine how the variability in the local microstructure affected the macroscopic behavior as well as the local damage initiation and progression. Results indicate that the locations where damage initiated and propagated is linked to specific microstructural features.

Mital, Subodh K.↗
NASA NTRS2013

Improved Subcell Model for the Prediction of Braided Composite Response

In this work, the modeling of triaxially braided composites was explored through a semi-analytical discretization. Four unique subcells, each approximated by a "mosaic" stacking of unidirectional composite plies, were modeled through the use of layered-shell elements within the explicit finite element code LS-DYNA. Two subcell discretizations were investigated: a model explicitly capturing pure matrix regions, and a novel model which absorbed pure matrix pockets into neighboring tow plies. The in-plane stiffness properties of both models, computed using bottom-up micromechanics, correlated well to experimental data. The absorbed matrix model, however, was found to best capture out-of- plane flexural properties by comparing numerical simulations of the out-of-plane displacements from single-ply tension tests to experimental full field data. This strong correlation of out-of-plane characteristics supports the current modeling approach as a viable candidate for future work involving impact simulations.

Cater, Christopher R.↗
NASA NTRS2012

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↗
NASA NTRS1981

Mineralogy, petrology, and trace element geochemistry of the Johnstown meteorite - A brecciated orthopyroxenite with siderophile and REE-rich components

The compositional and petrologic characteristics of the Johnstown meteorite show it to contain uncontaminated and unbrecciated orthopyroxenite clasts of cumulative origin that (1) must have undergone subsolidus recrystalization, (2) are parental to the brecciated matrix, and (3) show no evidence of a xenolithic, meteoritic contribution to the matrix except for contamination by the projectile which crushed it on impact. The trapped liquid was not introduced in the impact process. The variability of such trace elements as the light rare earth elements, and the presence of plagioclase and olivine in only one of the thin sections studied, demonstrates the heterogeneity of coarse-grained diogenites on a millimeter scale and the difficulty of obtaining representative samples of such meteorites. The data presented indicate that this meteorite is a monominct breccia.

Floran, R. J.↗
NASA NTRS2005

A Model for Siderophile Element Distribution in Planetary Differentiation

Planetary differentiation begins with partial melting of small planetesimals. At low degrees of partial melting, a sulfur-rich liquid segregates by physical mechanisms including deformation-assisted porous flow. Experimental studies of the physical mechanisms by which Fe-S melts segregate from the silicate matrix of a molten H chondrite are part of a companion paper. Geochemical studies of these experimental products revealed that metallic liquids were in equilibrium with residual metal in the H chondrite matrix. This contribution explores the geochemical signatures produced by early stages of core formation. Particularly, low-degree partial melt segregation of Fe-S liquids leaves residual metal in the silicate matrix. Some achondrites appear to be residues of partial melting, e.g., ureilites, which are known to contain metal. The metal in these achondrites may show a distinct elemental signature. To quantify the effect of sulfur on siderophile element contents of residual metal we have developed a model based on recent parametrizations of equilibrium solid metal-liquid metal partitioning experiments.

Humayun, M.↗
NASA NTRS2014

Origins and Distribution of Chondritic Olivine Inferred from Wild 2 and Chondrite Matrix

To date, only 180 particle impact tracks from Wild 2 have been extracted from the Stardust aerogel collector and even fewer have been thoroughly characterized. In order to provide a cohesive compositional dataset that can be compared to the meteorite record, we have made both major and minor element analyses (TEM/EDXS) of olivine and low-Ca pyroxene for 39 particles harvested from 26 tracks. However, the dearth of equivalent analyses for these phases in chondrite matrix hinders their comparison to the Wild 2 samples. To properly permit comparison of chondritic olivine and pyroxene to the Wild 2 samples, we have also provided a large, comprehensive EPMA dataset (greater than10(exp 3) analyses) of analogous grains (5-30 micrometers) isolated in L/LL3.0-4, CI, CM, CR, CH, CO, CV, Acfer 094, EH3, EL6, and Kakangari matrix

Frank, D. R.↗
NASA NTRSDate not supplied

Improved Plate and Beam Models for Thermoviscoelastic Constitutive Modeling of Composites

The effective properties of composites are influenced by the time-dependent behavior of polymer matrices very sensitive to changes in temperature. Improved plate and beam models are required to efficiently design, and simulate composite structures when the long-term performance of large anisotropic composite structures is the matter of interest. In this work, mechanics of structure genome (MSG) is used to con-struct linear thermoviscoelastic plate and beam models that can homogenize three-dimensional heterogeneous materials made of constituents with time- and temperature-dependent behavior. The formulation derives the transient strain energy based on integral formulation for thermorheologically simple materials subject to finite temperature changes with the restriction that the strain is small. The reduced time parameter is introduced to relate the time-temperature dependency of the anisotropic material by means of master curves at reference conditions. The new formulation has been implemented in SwiftCompTM, a general-purpose multiscale constitutive modeling code based on MSG. Experimental data and three-dimensional direct numerical simulations of thin-ply high-strain composites (TP-HSC) using a commercial finite element analysis (FEA) package are conducted to verify the accuracy of SwiftCompTM results. The paper also analyzes the relationship between the shift factor of the polymer matrix and the temperature dependencies of the effective beam properties.

Finite element analysis↗
DOE OSTI2023

A Study of a Surrogate Nuclear Thermal Propulsion Fuel Element

We report a surrogate fuel, Zirconium Carbide/Titanium Nitride has been prepared using a simple method to obtain a more uniform distribution of titanium nitride in the zirconium carbide matrix. Powders were sintered using Spark Plasma Sintering at 1750oC and 1850oC using a modified sintering profile. This resulted in a fairly uniform distribution of titanium nitride in the zirconium carbide matrix with a density of 95% of theoretical at 1850oC. Scanning electron microscopy, indentation hardness and x-ray diffraction results are given.

36 MATERIALS SCIENCE↗
DOE OSTI2026

PDE-constrained high-order mesh optimization

Here, we present a novel framework for PDE-constrained r-adaptivity of high-order meshes. The proposed method formulates mesh movement as an optimization problem, with an objective function defined as a convex combination of a mesh quality metric and a measure of the accuracy of the PDE solution obtained via finite element discretization. The proposed formulation achieves optimized, well-defined high-order meshes by integrating mesh quality control, PDE solution accuracy, and robust gradient regularization. We adopt the Target-Matrix Optimization Paradigm to control geometric properties across the mesh, independent of the PDE of interest. To incorporate the accuracy of the PDE solution, we introduce error measures that control the finite element discretization error. The implicit dependence of these error measures on the mesh nodal positions is accurately captured by adjoint sensitivity analysis. Additionally, a convolution-based gradient regularization strategy is used to ensure stable and effective adaptation of high-order meshes. We demonstrate that the proposed framework can improve mesh quality and reduce the error by up to 10 times for the solution of Poisson and linear elasto-static problems. The approach is general with respect to the dimensionality, the order of the mesh, the types of mesh elements, and can be applied to any PDE that admits well-defined adjoint operators.

Computer science↗
DOE OSTI2022

Krylov subspace recycling for evolving structures

Krylov subspace recycling is a powerful tool when solving a long series of large, sparse linear systems that change only slowly over time. In PDE constrained shape optimization, these series appear naturally, as typically hundreds or thousands of optimization steps are needed with only small changes in the geometry. In this setting, however, applying Krylov subspace recycling can be a difficult task. As the geometry evolves, in general, so does the finite element mesh defined on or representing this geometry, including the numbers of nodes and elements and element connectivity. This is especially the case if re-meshing techniques are used. As a result, the number of algebraic degrees of freedom in the system changes, and in general the linear system matrices resulting from the finite element discretization change size from one optimization step to the next. Changes in the mesh connectivity also lead to structural changes in the matrices. In the case of re-meshing, even if the geometry changes only a little, the corresponding mesh might differ substantially from the previous one. Obviously, this prevents any straightforward mapping of the approximate invariant subspace of the linear system matrix (the focus of recycling in this work) from one optimization step to the next; similar problems arise for other selected subspaces. In this paper, we present an algorithm to map an approximate invariant subspace of the linear system matrix for the previous optimization step to an approximate invariant subspace of the linear system matrix for the current optimization step, for general meshes. This is achieved by exploiting the map from coefficient vectors to finite element functions on the mesh, combined with interpolation or approximation of functions on the finite element mesh. We demonstrate the effectiveness of our approach numerically with several proof of concept studies for a specific meshing technique.

42 ENGINEERING↗
NASA NTRS2003

The Cytoskeleton and Force Response Mechanisms

The long term aim of this project was to define the mechanisms by which cells sense and respond to the physical forces experienced at 1g and missing in microgravity. Identification and characterization of the elements of the cells force response mechanism could provide pathways and molecules to serve as targets for pharmacological intervention to mitigate the pathologic effects of microgravity. Mechanical forces experienced by the organism can be transmitted to cells through molecules that allow cells to bind to the extracellular matrix and through other types of molecules which bind cells to each other. These molecules are coupled in large complexes of proteins to structural elements such as the actin cytoskeleton that give the cell the ability to sense, resist and respond to force. Application of small forces to tissue culture cells causes local elevation of intracellular calcium through stretch activated ion channels, increased tyrosine phosphorylation and a restructuring of the actin cytoskeleton. Using collagen coated iron oxide beads and strong magnets, we can apply different levels of force to cells in culture. We have found that force application causes the cells to polymerize actin at the site of mechanical deformation and unexpectedly, to depolymerize actin across the rest of the cell. Observations of GFP- actin expressing cells demonstrate that actin accumulates at the site of deformation within the first five minutes of force application and is maintained for many tens of minutes after force is removed. Consistent with the reinforcement of the cytoskeletal structures underlying the integrin-bead interaction, force also alters the motion of bound magnetic beads. This effect is seen following the removal of the magnetic field, and is only partially ablated by actin disruption with cytochalsin B. While actin is polymerizing locally at the site of force application, force also stimulates a global reduction in actin filament content within the cells. We have examined the roles of several actin filament disassembly factors in the global reduction of cellular actin filaments. The calcium regulated actin filament severing protein gelsolin is not necessary for the increased actin turnover, as cells derived from gelsolin null and wildtype mice still show a reduction in total actin filament content. Instead, our work suggests that the actin binding protein cofilin may be important for these changes in actin dynamics. Cofilin binds to and enhances the disassembly of actin filaments. Using immunological methods, we observe transient changes in the phosphorylation state of cofilin upon force application that suggests that cofilin may mediate actin filament turnover. Early after force application, cofilin is transiently dephosphorylated, activating its actin disassembly activity. Subsequently, we find a hyper-phosphorylation of cofilin, rendering it inactive. This reduction in cofilin activity may explain the stability of the force induced actin structuttes. In testing this hypothesis, we aimed to generate cells that express the constituitively active kinase (LIM-kinase) that phosphorylates cofilin. lnidial attempts in the cell lines used for the our previous studies proved unsuccessful. While we prepare this work for pubication, we are continuing to study other cell lines and tissue sources to determine whether they show a reduction in F-actin content after force application.

Allen, Philip Goodwin↗
NASA NTRS2019

Verification and Validation of a Three-Dimensional Composite Impact Model with Tabulated Input

A material model which incorporates several key capabilities which have been identified by the aerospace community as lacking in the composite impact models currently available in the commercial transient dynamic finite element code LS-DYNA® has been developed. The material model utilizes experimentally based tabulated input to define the evolution of plasticity and damage as opposed to specifying discrete input parameters (such as modulus and strength. The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. The capability to account for the rate and temperature dependent deformation response of composites has also been incorporated into the material model. For the damage model, a strain equivalent formulation is utilized to allow for the uncoupling of the deformation and damage analyses. In the damage model, a diagonal damage tensor is defined to account for the directionally dependent variation of damage. However, the terms in the damage matrix are semi-coupled such that the damage in a particular coordinate direction is a function of the stresses and plastic strains in all of the coordinate directions. For the failure model, a tabulated approach is utilized in which a stress or strain based invariant is defined as a function of the location of the current stress state in stress space to define the initiation of failure, which allows an arbitrarily shaped failure surface to be defined. A systematic series of validation and verification studies, at a variety of length scales ranging from single element simulations to simulations of a flat panel impact test, have been performed to fully exercise and evaluate the capabilities of the developed model.

Impact↗
DOE OSTI2025

ThinCurr: An open-source 3D thin-wall eddy current modeling code for the analysis of large-scale systems of conducting structures

In this paper we present a new thin-wall eddy current modeling code, ThinCurr, for studying inductively-coupled currents in 3D conducting structures -- with primary application focused on the interaction between currents flowing in coils, plasma, and conducting structures of magnetically-confined plasma devices. The code utilizes a boundary finite element method on an unstructured, triangular grid to accurately capture device structures. The new code, part of the broader Open FUSION Toolkit, is open-source and designed for ease of use without sacrificing capability and speed through a combination of Python, Fortran, and C/C++ components. Scalability to large models is enabled through use of hierarchical off-diagonal low-rank compression of the inductance matrix, which is otherwise dense. Ease of handling large models of complicated geometry is further supported by automatic determination of supplemental elements through a greedy homology approach. Here, a detailed description of the numerical methods of the code and verification of the implementation of those methods using cross-code comparisons against the VALEN code and Ansys commercial analysis software is shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1996

Some Experiences with Nonoverlapping Schur Complement Parallel Preconditioning for CFD Calculations

In this work we consider solving matrices which arise from the discretization of advection-diffusion field equations on arbitrary triangulated domains using stabilized numerical methods. The talk will discuss several candidate matrix preconditioning algorithms based on the 2 x 2 block factorization induced by an apriori partitioning of the triangulated domain. Application of the 2 x 2 block preconditioner requires the formation and inversion of the Schur complement submatrix. We consider several strategies for simplifying this task: incomplete Schur complement factorizations, drop tolerance element filling, Schur complement probing, and localized Schur complement inversion. Numerical results will be shown comparing performance and efficiency of these approximations. The matrix preconditioner has also been embedded into a Newton algorithm for solving the nonlinear Euler and Navier-Stokes equations governing compressible flow. The remainder of the talk will show numerous examples in CFD to demonstrate the efficiency and robustness of the techniques.

Barth, Timothy J.↗
NASA NTRS1977

Nebular condensation of moderately volatile elements and their abundances in ordinary chondrites

Two models accounting for the concentrations of moderately volatile elements in ordinary chondrites are discussed. The first relies on an apparent correlation between ordinary chondrite/type-I carbonaceous chondrite abundance ratios and ideal solution-condensation temperatures of such elements as Sn, Cu, Li, Na, and Mn to construct a model in which volatiles are lost as gases before nebular condensation or as finely-divided solids that are incompletely agglomerated, and the condensation/agglomeration efficiency declines as a function of time. The second theory largely discounts this gas-dust diffraction effect and instead proposed condensation on matrix and volatization from chondules and coarse metal grains as mechanisms explaining the apparent existence of three distinct classes of volatile-element abundance in ordinary chondrites.

Wai, C. M.↗
NASA NTRS2009

Economical Implementation of a Filter Engine in an FPGA

A logic design has been conceived for a field-programmable gate array (FPGA) that would implement a complex system of multiple digital state-space filters. The main innovative aspect of this design lies in providing for reuse of parts of the FPGA hardware to perform different parts of the filter computations at different times, in such a manner as to enable the timely performance of all required computations in the face of limitations on available FPGA hardware resources. The implementation of the digital state-space filter involves matrix vector multiplications, which, in the absence of the present innovation, would ordinarily necessitate some multiplexing of vector elements and/or routing of data flows along multiple paths. The design concept calls for implementing vector registers as shift registers to simplify operand access to multipliers and accumulators, obviating both multiplexing and routing of data along multiple paths. Each vector register would be reused for different parts of a calculation. Outputs would always be drawn from the same register, and inputs would always be loaded into the same register. A simple state machine would control each filter. The output of a given filter would be passed to the next filter, accompanied by a "valid" signal, which would start the state machine of the next filter. Multiple filter modules would share a multiplication/accumulation arithmetic unit. The filter computations would be timed by use of a clock having a frequency high enough, relative to the input and output data rate, to provide enough cycles for matrix and vector arithmetic operations. This design concept could prove beneficial in numerous applications in which digital filters are used and/or vectors are multiplied by coefficient matrices. Examples of such applications include general signal processing, filtering of signals in control systems, processing of geophysical measurements, and medical imaging. For these and other applications, it could be advantageous to combine compact FPGA digital filter implementations with other application-specific logic implementations on single integrated-circuit chips. An FPGA could readily be tailored to implement a variety of filters because the filter coefficients would be loaded into memory at startup.

Kowalski, James E.↗
DOE OSTI2022

Towards underwater additive manufacturing via additive friction stir deposition

Given the challenges in feed material supply and quality control, metal additive manufacturing has rarely been implemented in austere environments, especially underwater. This paper explores the underwater operation potential of an emerging solid-state additive technology: additive friction stir deposition, wherein material feeding and bonding are enabled by mechanical forces with minimal influences from water. It is demonstrated that additive friction stir deposition of 304 stainless steel can be successfully performed with the print head and substrate immersed in water. High temperature is reached in the deposition zone (>60% melting temperature); the material deposition behavior is similar to that in typical open-air operation. The as-deposited material is fully-dense, having fewer annealing twins and a substantially smaller grain size than the feed material (4.98 μm vs. 31.44 μm ). Such microstructural changes stem from dynamic recrystallization caused by the large strain and high temperature introduced during deposition. In addition to grain refinement, small equiaxed dispersoids (~2–3 μm or less) are formed and evenly distributed in the austenite steel matrix. Rich in Cr, Mn, and O, these particles likely result from the reaction between the elements in stainless steel and water at elevated temperatures.

16 TIDAL AND WAVE POWER↗