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,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

An Automated Meshing Framework for Progressive Damage Analysis of Fabrics Using CompDam

A meshing and model construction framework for progressive damage analysis of fabric composites is proposed. A mesoscale idealization is used such that the fiber tows are considered homogeneous and transversely isotropic continua, and the matrix is a homogeneous isotropic continuum. By using this idealization, relatively mature damage modeling techniques developed for laminates of unidirectional plies are transferrable to fabrics. Cohesive interface elements and a continuum damage mechanics (CDM) model implemented in the software CompDam are used to account for the key damage modes including cracking at the tow/matrix interface, transverse matrix cracking in the tows, and fiber breakage in the tows. Matrix degradation is accounted for by plasticity. The modeling approach is demonstrated for a 3D orthogonal woven fabric. The main novelty of the proposed model construction framework is in its flexibility resulting from a formulation that is agnostic to the fabric architecture and implementation using mesh based operations. The algorithm does not make any assumptions about the tow sizes, shape, path, or proximity to other tows, thus it is applicable to a broad range of fabric architectures. The matrix region mesh surrounding the tows is constructed from domain boundary and tow surface meshes (i.e. through a series of mesh-based operations) such that the model domain can be arbitrarily shaped.

Andrew Cole Bergan↗
DOE OSTI2025

Multistate resistance in TaN/(Hf,Zr)O 2 /Ta ferroelectric tunnel junctions

Ferroelectric tunnel junctions (FTJs) utilizing hafnium zirconium oxide (HZO) have emerged as promising non-volatile memory elements for microelectronics, compatible with back end of line (BEOL) complementary–metal–oxide semiconductor fabrication. This study investigates asymmetric electrode TaN/HZO/Ta devices with a 6 nm thick HZO layer as FTJs for multistate resistive memory applications. The individual FTJs exhibit a resistance ratio exceeding 10× when utilized as a binary state device, with pulsing between −1.7 and +1.4 V to set the high resistance state (HRS) and low resistance state (LRS), respectively. Following with reduced write voltage pulses allows the ferroelectric device to operate with a selection of over 32 distinct resistance states (2 5 bits) between the LRS and HRS. This work then explores the stability of the resistance states during write/read pulse cycling, along with the stability of the state after multiple read pulses. Accessing the multibit state shows stability within 50 reads with the binary state remaining stable for more than 4000 reads pulses. With their multistate tunability and versatility, FTJs hold promise as BEOL memory elements for compute-in-memory (CiM) arrays, binary digital memory, or weighted vector matrix multiplication applications with low power consumption during computations.

CMOS↗
DOE OSTI2026

A limit on the total lepton number in the Universe from BBN and the CMB

At temperatures below the QCD phase transition, any substantial lepton number in the Universe can only be present within the neutrino sector. In this work, we systematically explore the impact of a non-vanishing lepton number on Big Bang Nucleosynthesis (BBN) and the Cosmic Microwave Background (CMB). Relying on our recently developed framework based on momentum averaged quantum kinetic equations for the neutrino density matrix, we solve the full BBN reaction network to obtain the abundances of primordial elements. We find that the maximal primordial total lepton number L allowed by BBN and the CMB is -0.12 (-0.10) ≤ L ≤ 0.13 (0.12) for NH (IH), while specific flavor directions can be even more constrained. This bound is complementary to the limits obtained from avoiding baryon overproduction through sphaleron processes at the electroweak phase transition since, although numerically weaker, it applies at lower temperatures and is obtained completely independently. We publicly release the C++ code COFLASY-C on GitHub (https://github.com/mariofnavarro/COFLASY/tree/COFLASY-C) which solves for the evolution of the neutrino quantum kinetic equations numerically.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Improved precision in As speciation analysis with HERFD-XANES at the As K -edge: the case of As speciation in mine waste

High-energy-resolution fluorescence-detected (HERFD) X-ray absorption near-edge spectroscopy (XANES) is a spectroscopic method that allows for increased spectral feature resolution, and greater selectivity to decrease complex matrix effects compared with conventional XANES. XANES is an ideal tool for speciation of elements in solid-phase environmental samples. Accurate speciation of As in mine waste materials is important for understanding the mobility and toxicity of As in near-surface environments. In this study, linear combination fitting (LCF) was performed on synthetic spectra generated from mixtures of eight measured reference compounds for both HERFD-XANES and transmission-detected XANES to evaluate the improvement in quantitative speciation with HERFD-XANES spectra. The reference compounds arsenolite (As 2 O 3 ), orpiment (As 2 S 3 ), getchellite (AsSbS 3 ), arsenopyrite (FeAsS), kaňkite (FeAsO 4 ·3.5H 2 O), scorodite (FeAsO 4 ·2H 2 O), sodium arsenate (Na 3 AsO 4 ), and realgar (As 4 S 4 ) were selected for their importance in mine waste systems. Statistical methods of principal component analysis and target transformation were employed to determine whether HERFD improves identification of the components in a dataset of mixtures of reference compounds. LCF was performed on HERFD- and total fluorescence yield (TFY)-XANES spectra collected from mine waste samples. Arsenopyrite, arsenolite, orpiment, and sodium arsenate were more accurately identified in the synthetic HERFD-XANES spectra compared with the transmission-XANES spectra. In mine waste samples containing arsenopyrite and either scorodite or kaňkite, LCF with HERFD-XANES measurements resulted in fits with smaller R -factors than concurrently collected TFY measurements. The improved accuracy of HERFD-XANES analysis may provide enhanced delineation of As phases controlling biogeochemical reactions in mine wastes, contaminated soils, and remediation systems.

58 GEOSCIENCES↗
DOE OSTI2020

Multivariate Curve Resolution (MCR) using Principal Components Inputs and Rigorous Equality and Inequality Constraints

This MATLAB pseudocode perform multivariate curve resolution (MCR) using PCA scores & loadings of data as inputs. It employs rigorous least squares equality and inequality constraints for all elements in the solution factor matrices. It also can be used to perform nonnegative matrix factorization. SAND2020-12650 M Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525.

Van Benthem, Mark↗
NASA NTRS1975

Feasibility study for the implementation of NASTRAN on the ILLIAC 4 parallel processor

The ILLIAC IV, a fourth generation multiprocessor using parallel processing hardware concepts, is operational at Moffett Field, California. Its capability to excel at matrix manipulation, makes the ILLIAC well suited for performing structural analyses using the finite element displacement method. The feasibility of modifying the NASTRAN (NASA structural analysis) computer program to make effective use of the ILLIAC IV was investigated. The characteristics are summarized of the ILLIAC and the ARPANET, a telecommunications network which spans the continent making the ILLIAC accessible to nearly all major industrial centers in the United States. Two distinct approaches are studied: retaining NASTRAN as it now operates on many of the host computers of the ARPANET to process the input and output while using the ILLIAC only for the major computational tasks, and installing NASTRAN to operate entirely in the ILLIAC environment. Though both alternatives offer similar and significant increases in computational speed over modern third generation processors, the full installation of NASTRAN on the ILLIAC is recommended. Specifications are presented for performing that task with manpower estimates and schedules to correspond.

Field, E. I.↗
NASA NTRS1978

Finite element analysis of periodic transonic flow problems

Flow about an oscillating thin airfoil in a transonic stream was considered. It was assumed that the flow field can be decomposed into a mean flow plus a periodic perturbation. On the surface of the airfoil the usual Neumman conditions are imposed. Two computer programs were written, both using linear basis functions over triangles for the finite element space. The first program uses a banded Gaussian elimination solver to solve the matrix problem, while the second uses an iterative technique, namely SOR. The only results obtained are for an oscillating flat plate.

Fix, G. J.↗
NASA NTRS1982

Aeroelastic analysis for helicopter rotors with blade appended pendulum vibration absorbers. Mathematical derivations and program user's manual

Mathematical development is presented for the expanded capabilities of the United Technologies Research Center (UTRC) G400 Rotor Aeroelastic Analysis. This expanded analysis, G400PA, simulates the dynamics of teetered rotors, blade pendulum vibration absorbers and the higher harmonic excitations resulting from prescribed vibratory hub motions and higher harmonic blade pitch control. Formulations are also presented for calculating the rotor impedance matrix appropriate to these higher harmonic blade excitations. This impedance matrix and the associated vibratory hub loads are intended as the rotor blade characteristics elements for use in the Simplified Coupled Rotor/Fuselage Vibration Analysis (SIMVIB). Sections are included presenting updates to the development of the original G400 theory, and material appropriate to the user of the G400PA computer program. This material includes: (1) a general descriptionof the tructuring of the G400PA FORTRAN coding, (2) a detaild description of the required input data and other useful information for successfully running the program, and (3) a detailed description of the output results.

Bielawa, R. L.↗
NASA NTRS1981

Accretionary dark rims in unequilibrated chondrites

Textural and qualitative EDX investigations of dark-rimmed particles in six low petrologic type chondrites indicate that the rims accreted on host particles over a wide range of temperatures prior to initial accumulation and lithification of the meteorites in which the rimmed particles are now contained. Many dark rims are enriched in moderately volatile trace elements such as Na, Cl, P, and K, relative to the host particles and matrix. The range of physical/chemical environments associated with hypervelocity impacts may have offered the setting for the formation of dark-rimmed particles early in solar system history.

King, T. V. V.↗
NASA NTRS1985

Connection between variations of the atmosphere temperature profile and variations of the meson component intensity

The influence of temperature effects on intensity variations of the cosmic ray meson component were studied. The connection between the temperature variation delta T and the intensity variation delta I was established by using the temperature coefficient density technique. To realize how many devices are needed on the Earth for predicting the temperature variation of the atmosphere profile with a reasonable accuracy, IO isobaric levels and IO were calculated. The set of initial elements of the cosmic ray mesons are varied and it is shown that the matrix of the coefficients B sub ij is rather sensitive to their choice. It is found that if for the calculations of the atmospheric temperature variations the model is used, the number of meson components, essentially exceeding 3, should be considered.

Blokh, Y. L.↗
NASA NTRS1986

Preconditioned conjugate residual methods for the solution of spectral equations

Conjugate residual methods for the solution of spectral equations are described. An inexact finite-difference operator is introduced as a preconditioner in the iterative procedures. Application of these techniques is limited to problems for which the symmetric part of the coefficient matrix is positive definite. Although the spectral equation is a very ill-conditioned and full matrix problem, the computational effort of the present iterative methods for solving such a system is comparable to that for the sparse matrix equations obtained from the application of either finite-difference or finite-element methods to the same problems. Numerical experiments are shown for a self-adjoint elliptic partial differential equation with Dirichlet boundary conditions, and comparison with other solution procedures for spectral equations is presented.

Wong, Y. S.↗
NASA NTRS1986

Reduced conservatism in stability robustness bounds by state transformation

This note addresses the issue of 'conservatism' in the time domain stability robustness bounds obtained by the Liapunov approach. A state transformation is employed to improve the upper bounds on the linear time-varying perturbation of an asymptotically stable linear time-invariant system for robust stability. This improvement is due to the variance of the conservatism of the Liapunov approach with respect to the basis of the vector space in which the Liapunov function is constructed. Improved bounds are obtained, using a transformation, on elemental and vector norms of perturbations (i.e., structured perturbations) as well as on a matrix norm of perturbations (i.e., unstructured perturbations). For the case of a diagonal transformation, an algorithm is proposed to find the 'optimal' transformation. Several examples are presented to illustrate the proposed analysis.

Yedavalli, R. K.↗
NASA NTRS1986

The RP-2 sailplane

The all-composite design and fabrication methodology employed by Rensselaer Polytechnic Intitute's 'CAPGLIDE' program in producing such sailplanes as the RP-2 is discussed. The RP-2's production methods, which furnish hands-on experience for the students enrolled in the program, encompass wet lay-up and vacuum bag procedures, as well as autoclave curing of prepreg parts employing graphite, kevlar, and glass fiber reinforcements. Epoxy resins are the primary matrix materials, and are used in conjunction with polymer foams for the sandwich panel elements of the sailplane structure.

Winckler, Steven↗
NASA NTRS1993

Global dynamic modeling of a transmission system

The work performed on global dynamic simulation and noise correlation of gear transmission systems at the University of Akron is outlined. The objective is to develop a comprehensive procedure to simulate the dynamics of the gear transmission system coupled with the effects of gear box vibrations. The developed numerical model is benchmarked with results from experimental tests at NASA Lewis Research Center. The modal synthesis approach is used to develop the global transient vibration analysis procedure used in the model. Modal dynamic characteristics of the rotor-gear-bearing system are calculated by the matrix transfer method while those of the gear box are evaluated by the finite element method (NASTRAN). A three-dimensional, axial-lateral coupled bearing model is used to couple the rotor vibrations with the gear box motion. The vibrations between the individual rotor systems are coupled through the nonlinear gear mesh interactions. The global equations of motion are solved in modal coordinates and the transient vibration of the system is evaluated by a variable time-stepping integration scheme. The relationship between housing vibration and resulting noise of the gear transmission system is generated by linear transfer functions using experimental data. A nonlinear relationship of the noise components to the fundamental mesh frequency is developed using the hypercoherence function. The numerically simulated vibrations and predicted noise of the gear transmission system are compared with the experimental results from the gear noise test rig at NASA Lewis Research Center. Results of the comparison indicate that the global dynamic model developed can accurately simulate the dynamics of a gear transmission system.

Choy, F. K.↗
NASA NTRS1994

Accelerated panel methods using the fast multipole method

Panel methods are commonly used in computational fluid dynamics for the solution of potential flow problems. The methods are a numerical technique based on the surface distribution of singularity elements. The solution is the process of finding the strength of the singularity elements distributed over the body's surface. This process involves the solution of the matrix problem Pq = p' for a set of unknowns q. The Fast Multipole Method is used to directly compute q without using matrix solvers. The algorithm works in O(N) time for N points, a great improvement over standard matrix solvers. In panel methods, the surface of a body is divided into a series of quadrilateral panels. The methods involve the computation of the influence of all other panels on each individual panel. The influence is based on the surface distribution, though this can be approximated by the area for distant panels. An alternative approximation, though with arbitrary accuracy, is to develop a multipole expansion about the center of the panel to describe the effect of a given panel on distant points in space. The expansion is based on the moments of the panel, thus allow the use of various surface distributions without changing the basic algorithm, just the computation of the various moments. The expansions are then manipulated in a tree walk to develop Taylor series expansions about a point in space which describe the effect of all distant panels on any point within a volume of convergence. The effect of near panels then needs to be computed directly, but the effect of all distant panels can be computed by simply evaluating the resulting expansion. The Fast Multipole Method has been applied to panel methods for the solution of source and doublet distributions. A major feature of the algorithm is that the algorithm does not change to derive the potential and velocity for sources and doublets. The same expansions can be used for both sources and doublets. Since the velocity is related to the potential, and the doublet potential is related to the z-component of the source velocity, all values can be derived from the same expansion by taking a series of partial derivatives. This requires more expansion terms to be kept since terms are lost in the process of taking partial derivatives. Thus to maintain accuracy for the doublet computation, more terms are required than if just evaluating for sources. The resulting Fast Multipole code should then parallelize better than classical panel methods due to the locality of data dependencies found in the Fast Multipole Method. Theoretically the parallelized code should execute in O(log N) time with O(N) processors, though this is not practical. Ongoing work includes implementing the parallel accelerated panel method, including methods to improve the load balancing of the problem by taking advantage of the known geometry of panels, and to encorporate sensitivity analysis into the algorithm.

Leathrum, James F., Jr.↗
NASA NTRS1999

Development of GENOA Progressive Failure Parallel Processing Software Systems

A capability consisting of software development and experimental techniques has been developed and is described. The capability is integrated into GENOA-PFA to model polymer matrix composite (PMC) structures. The capability considers the physics and mechanics of composite materials and structure by integration of a hierarchical multilevel macro-scale (lamina, laminate, and structure) and micro scale (fiber, matrix, and interface) simulation analyses. The modeling involves (1) ply layering methodology utilizing FEM elements with through-the-thickness representation, (2) simulation of effects of material defects and conditions (e.g., voids, fiber waviness, and residual stress) on global static and cyclic fatigue strengths, (3) including material nonlinearities (by updating properties periodically) and geometrical nonlinearities (by Lagrangian updating), (4) simulating crack initiation. and growth to failure under static, cyclic, creep, and impact loads. (5) progressive fracture analysis to determine durability and damage tolerance. (6) identifying the percent contribution of various possible composite failure modes involved in critical damage events. and (7) determining sensitivities of failure modes to design parameters (e.g., fiber volume fraction, ply thickness, fiber orientation. and adhesive-bond thickness). GENOA-PFA progressive failure analysis is now ready for use to investigate the effects on structural responses to PMC material degradation from damage induced by static, cyclic (fatigue). creep, and impact loading in 2D/3D PMC structures subjected to hygrothermal environments. Its use will significantly facilitate targeting design parameter changes that will be most effective in reducing the probability of a given failure mode occurring.

Abdi, Frank↗
NASA NTRS2000

Validated Feasibility Study of Integrally Stiffened Metallic Fuselage Panels for Reducing Manufacturing Costs: Cost Assessment of Manufacturing/Design Concepts

The objective of the Integral Airframe Structures (IAS) program was to demonstrate, for an integrally stiffened structural concept, performance and weight equal to "built-up" structure with lower manufacturing cost. This report presents results of the cost assessment for several design configuration/manufacturing method combinations. The attributes of various cost analysis models were evaluated and COSTRAN selected for this study. A process/design cost evaluation matrix was developed based on material, forming, machining, and assembly of structural sub-elements and assembled structure. A hybrid design, made from high-speed machined extruded frames that are mechanically fastened to high-speed machined plate skin/stringer panels, was identified as the most cost-effective manufacturing solution. Recurring labor and material costs of the hybrid design are up to 61 percent less than the current built-up technology baseline. This would correspond to a total cost reduction of $1.7 million per ship set for a 777-sized airplane. However, there are important outstanding issues with regard to the cost of capacity of high technology machinery, and the ability to cost-effectively provide surface finish acceptable to the commercial aircraft industry. The projected high raw material cost of large extrusions also played an important role in the trade-off between plate and extruded concepts.

Metschan, S.↗
NASA NTRS2009

Pixelized Device Control Actuators for Large Adaptive Optics

A fully integrated, compact, adaptive space optic mirror assembly has been developed, incorporating new advances in ultralight, high-performance composite mirrors. The composite mirrors use Q-switch matrix architecture-based pixelized control (PMN-PT) actuators, which achieve high-performance, large adaptive optic capability, while reducing the weight of present adaptive optic systems. The self-contained, fully assembled, 11x11x4-in. (approx.= 28x28x10-cm) unit integrates a very-high-performance 8-in. (approx.=20-cm) optic, and has 8-kHz true bandwidth. The assembled unit weighs less than 15 pounds (=6.8 kg), including all mechanical assemblies, power electronics, control electronics, drive electronics, face sheet, wiring, and cabling. It requires just three wires to be attached (power, ground, and signal) for full-function systems integration, and uses a steel-frame and epoxied electronics. The three main innovations are: 1. Ultralightweight composite optics: A new replication method for fabrication of very thin composite 20-cm-diameter laminate face sheets with good as-fabricated optical figure was developed. The approach is a new mandrel resin surface deposition onto previously fabricated thin composite laminates. 2. Matrix (regenerative) power topology: Waveform correction can be achieved across an entire face sheet at 6 kHz, even for large actuator counts. In practice, it was found to be better to develop a quadrant drive, that is, four quadrants of 169 actuators behind the face sheet. Each quadrant has a single, small, regenerative power supply driving all 169 actuators at 8 kHz in effective parallel. 3. Q-switch drive architecture: The Q-switch innovation is at the heart of the matrix architecture, and allows for a very fast current draw into a desired actuator element in 120 counts of a MHz clock without any actuator coupling.

Knowles, Gareth J.↗