Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix compression”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Slow plastic strain rate compressive flow in binary CoAl intermetallics

Constant-velocity elevated temperature compression tests have been conducted on a series of binary CoAl intermetallics produced by hot extrusion of blended prealloyed powders. The as-extruded materials were polycrystalline, and they retained their nominal 10-micron grain size after being tested between 1100 and 1400 K at strain rates ranging from 2 x 10 to the -4th to 2 x 10 to the -7th per sec. Significant plastic flow was obtained in all cases; while cracking was observed, much of this could be due to failure at matrix-oxide interfaces along extrusion stringers rather than to solely intergranular fracture. A maximum in flow strength occurs at an aluminum-to-cobalt ratio of 0.975, and the stress exponent appears to be constant for aluminum-to-cobalt ratios of 0.85 or more. It is likely that very aluminum-deficient materials deform by a different mechanism than do other compositions.

Whittenberger, J. D.↗
NASA NTRS1991

Unified aeroacoustics analysis for high speed turboprop aerodynamics and noise. Volume 4: Computer user's manual for UAAP turboprop aeroacoustic code

The Unified AeroAcoustic Program (UAAP) code calculates the airloads on a single rotation prop-fan, or propeller, and couples these airloads with an acoustic radiation theory, to provide estimates of near-field or far-field noise levels. The steady airloads can also be used to calculate the nonuniform velocity components in the propeller wake. The airloads are calculated using a three dimensional compressible panel method which considers the effects of thin, cambered, multiple blades which may be highly swept. These airloads may be either steady or unsteady. The acoustic model uses the blade thickness distribution and the steady or unsteady aerodynamic loads to calculate the acoustic radiation. The users manual for the UAAP code is divided into five sections: general code description; input description; output description; system description; and error codes. The user must have access to IMSL10 libraries (MATH and SFUN) for numerous calls made for Bessel functions and matrix inversion. For plotted output users must modify the dummy calls to plotting routines included in the code to system-specific calls appropriate to the user's installation.

Menthe, R. W.↗
NASA NTRS1977

A finite element algorithm for sound propagation in axisymmetric ducts containing compressible mean flow

An accurate mathematical model for sound propagation in axisymmetric aircraft engine ducts with compressible mean flow is reported. The model is based on the usual perturbation of the basic fluid mechanics equations for small motions. Mean flow parameters are derived in the absence of fluctuating quantities and are then substituted into the equations for the acoustic quantities which were linearized by eliminating higher order terms. Mean swirl is assumed to be zero from the restriction of axisymmetry. A linear rectangular serendipity element is formulated from these equations using a Galerkin procedure and assembled in a special purpose computer program in which the matrix map for a rectangular mesh was specifically coded. Representations of the fluctuating quantities, mean quantities and coordinate transformations are isoparametric. The global matrix is solved by foreward and back substitution following an L-U decomposition with pivoting restricted internally to the blocks. Results from the model were compared with results from several alternative analyses and yielded satisfactory agreement.

Abrahamson, A. L.↗
NASA NTRS2006

Fuel Cells Utilizing Oxygen From Air at Low Pressures

A fuel cell stack has been developed to supply power for a high-altitude aircraft with a minimum of air handling. The fuel cell is capable of utilizing oxygen from ambient air at low pressure with no need for compression. For such an application, it is advantageous to take oxygen from the air (in contradistinction to carrying a supply of oxygen onboard), but it is a challenging problem to design a fuel-cell stack of reasonable weight that can generate sufficient power while operating at reduced pressures. The present fuel-cell design is a response to this challenge. The design features a novel bipolar plate structure in combination with a gas-diffusion structure based on a conductive metal core and a carbon gas-diffusion matrix. This combination makes it possible for the flow fields in the stack to have a large open fraction (ratio between open volume and total volume) to permit large volumes of air to flow through with exceptionally low backpressure. Operations at reduced pressure require a corresponding increase in the volume of air that must be handled to deliver the same number of moles of oxygen to the anodes. Moreover, the increase in the open fraction, relative to that of a comparable prior fuel-cell design, reduces the mass of the stack. The fuel cell has been demonstrated to operate at a power density as high as 105 W/cm2 at an air pressure as low as 2 psia (absolute pressure 14 kPa), which is the atmospheric pressure at an altitude of about 50,000 ft ( 15.2 km). The improvements in the design of this fuel cell could be incorporated into designs of other fuel cells to make them lighter in weight and effective at altitudes higher than those of prior designs. Potential commercial applications for these improvements include most applications now under consideration for fuel cells.

Cisar, Alan↗
DOE OSTI2024

Understanding interfacial crystallization dynamics on carbon fiber reinforced polypropylene composite manufacturing

Reinforcing polymers with discontinuous fibers improves mechanical properties, such as strength and stiffness, and in some cases achieve isotropic properties, rendering them suitable for various engineering applications. Matrix materials are generally highly engineered thermosets (e.g. crosslinked epoxies), bonded to the fiber periphery by proprietary surface and sizing chemistries. Semicrystalline thermoplastic matrices are less utilized due to poor fiber-matrix bonding resulting in inefficient interfacial load-transfer in reinforced composites. However, flexibility with melt-processing or molding conditions can be leveraged to promote non-covalent interfacial bonding between matrix and fiber via crystallization of the matrix onto fiber surface. In the present study, we utilize a co-mingle chopped carbon and isotactic polypropylene fibers to form isotropic composites, tailoring interfacial immobilized matrix or interphase morphology to optimize performance through precise control of thermal processing/molding windows. Calorimetry and optical microscopy were employed to investigate the impact of carbon fiber at various volume fractions (10, 20, and 30 %) on isotactic polypropylene crystallization and mechanical performance. Variations in mechanical properties correspond to the structural evolution of the interfacial region and are correlated to underlying microstructural attributes using wide-angle X-ray scattering, thermal analysis, and low-field nuclear magnetic resonance spectroscopy. These results provide a practical framework for the manufacturing of thermoplastic matrix composites. In conclusion, the results presented provide a guide for the strategic optimization of interphase design, showcasing tailorable tensile strengths which outperform any isotactic polypropylene carbon fiber composites previously reported in literature.

36 MATERIALS SCIENCE↗
NASA NTRS2012

Some Examples of the Relations Between Processing and Damage Tolerance

Most structures made of laminated polymer matrix composites (PMCs) must be designed to some damage tolerance requirement that includes foreign object impact damage. Thus from the beginning of a part s life, impact damage is assumed to exist in the material and the part is designed to carry the required load with the prescribed impact damage present. By doing this, some processing defects may automatically be accounted for in the reduced design allowable due to these impacts. This paper will present examples of how a given level of impact damage and certain processing defects affect the compression strength of a laminate that contains both. Knowledge of the impact damage tolerance requirements, before processing begins, can broaden material options and processing techniques since the structure is not being designed to pristine properties.

Nettles, Alan T.↗
DOE OSTI2021

Ocean Energy Systems Wave Energy Modeling Task 10.4: Numerical Modeling of a Fixed Oscillating Water Column

This paper reports on an ongoing international effort to establish guidelines for numerical modeling of wave energy converters, initiated by the International Energy Agency Technology Collaboration Program for Ocean Energy Systems. Initial results for point absorbers were presented in previous work, and here we present results for a breakwater-mounted Oscillating Water Column (OWC) device. The experimental model is at scale 1:4 relative to a full-scale installation in a water depth of 12.8 m. The power-extracting air turbine is modeled by an orifice plate of 1–2% of the internal chamber surface area. Measurements of chamber surface elevation, air flow through the orifice, and pressure difference across the orifice are compared with numerical calculations using both weakly-nonlinear potential flow theory and computational fluid dynamics. Both compressible- and incompressible-flow models are considered, and the effects of air compressibility are found to have a significant influence on the motion of the internal chamber surface. Recommendations are made for reducing uncertainties in future experimental campaigns, which are critical to enable firm conclusions to be drawn about the relative accuracy of the numerical models. It is well-known that boundary element method solutions of the linear potential flow problem (e.g., WAMIT) are singular at infinite frequency when panels are placed directly on the free surface. This is problematic for time-domain solutions where the value of the added mass matrix at infinite frequency is critical, especially for OWC chambers, which are modeled by zero-mass elements on the free surface. A straightforward rational procedure is described to replace ad-hoc solutions to this problem that have been proposed in the literature.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Online randomized interpolative decomposition with a posteriori error estimator for temporal PDE data reduction

Traditional low-rank approximation is a powerful tool for compressing large data matrices that arise in simulations of partial differential equations (PDEs), but suffers from high computational cost and requires several passes over the PDE data. The compressed data may also lack interpretability thus making it difficult to identify feature patterns from the original data. Here, to address these issues, we present an online randomized algorithm to compute the interpolative decomposition (ID) of large-scale data matrices in situ. Compared to previous randomized IDs that used the QR decomposition to determine the column basis, we adopt a streaming ridge leverage score-based column subset selection algorithm that dynamically selects proper basis columns from the data and thus avoids an extra pass over the data to compute the coefficient matrix of the ID. In particular, we adopt a single-pass error estimator based on the non-adaptive Hutch++ algorithm to provide real-time error approximation for determining the best coefficients. As a result, our approach only needs a single pass over the original data and thus is suitable for large and high-dimensional matrices stored outside of core memory or generated in PDE simulations. A strategy to improve the accuracy of the reconstructed data gradient, when desired, within the ID framework is also presented. We provide numerical experiments on turbulent channel flow and ignition simulations, and on the NSTX Gas Puff Image dataset, comparing our algorithm with the offline ID algorithm to demonstrate its utility in real-world applications.

Column subset selection↗
DOE OSTI2020

Simulation-driven optimization of high-order meshes in ALE hydrodynamics

Here we propose tools for high-order mesh optimization and demonstrate their benefits in the context of multi-material Arbitrary Lagrangian-Eulerian (ALE) compressible shock hydrodynamic applications. The mesh optimization process is driven by information provided by the simulation which uses the optimized mesh, such as shock positions, material regions, known error estimates, etc. These simulation features are usually represented discretely, for instance, as finite element functions on the Lagrangian mesh. The discrete nature of the input is critical for the practical applicability of the algorithms we propose and distinguishes this work from approaches that strictly require analytical information. Our methods are based on node movement through a high-order extension of the Target-Matrix Optimization Paradigm (TMOP). The proposed formulation is fully algebraic and relies only on local Jacobian matrices, so it is applicable to all types of mesh elements, in 2D and 3D, and any order of the mesh. We discuss the notions of constructing adaptive target matrices and obtaining their derivatives, reconstructing discrete data in intermediate meshes, node limiting that enables improvement of global mesh quality while preserving space-dependent local mesh features, and appropriate normalization of the objective function. The adaptivity methods are combined with automatic ALE triggers that can provide robustness of the mesh evolution and avoid excessive remap procedures. The benefits of the new high-order TMOP technology are illustrated on several simulations performed in the high-order ALE application BLAST.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Tuning the temperature range of superelastic Ni-Ti alloys for elastocaloric cooling via thermal processing

Abstract Caloric cooling enlisting solid-state refrigerants is potentially a promising eco-friendly alternative to conventional cooling based on vapor compression. The most common refrigerant materials for elastocaloric cooling to date are Ni-Ti based superelastic shape memory alloys. Here, we have explored tuning the operation temperature range of Ni 50.8 Ti 49.2 for elastocaloric cooling. In particular, we have studied the effect of thermal treatments (a.k.a. aging) on the transformation temperature, superelasticity, and elastocaloric effects of Ni 50.8 Ti 49.2 shape memory alloy tubes. The isothermal compressive test revealed that the residual strain of thermally-treated Ni-Ti tubes at room temperature approaches zero as aging time is increased. Short-time aging treatment at 400 °C resulted in good superelasticity and elastocaloric cooling performance with a large tunable austenite finish ( A f ) temperature range of 24.7 °C, as determined from the A f temperature of the samples that were aged 5–120 min. The main reason of the property change is the formation of a different amount of Ni 4 Ti 3 precipitates in the NiTi matrix. Our findings show that it is possible to tailor the A f temperature range for development of cascade elastocaloric cooling systems by thermally treating a starting single composition Ni-Ti alloy.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Post-impact fatigue of cross-plied, through-the-thickness reinforced carbon/epoxy composites

An experimental investigation of the post-impact fatigue response of integrally woven carbon/epoxy composites was conducted. Five different through-the-thickness (TTT) reinforcing fibers were used in an experimental textile process that produced an integrally woven (0/90/0/90/0/90/0/90/0)(sub T) ply layup with 21K AS4 carbon tow fiber. The resin was Hercules 3501-6, and the five TTT reinforcing fibers were Kevlar, Toray carbon, AS4 carbon, glass, and IM6 carbon. The purpose of this investigation was to study the post-impact fatigue response of these material systems and to identify the optimum TTT fiber. Samples were impacted with one half inch diameter aluminum balls with an average velocity of 543 ft/sec. Post-impact static compression and constant amplitude tension-compression fatigue tests were conducted. Fatigue tests were conducted with a loading ratio of R=-5, and frequency of 4 Hz. Damage growth was monitored using x-radiographic and sectioning techniques and by examining the stress-strain response (across the impact site) throughout the fatigue tests. The static compressive stress versus far-field strain response was nearly linear for all material groups. All the samples had a transverse shear failure mode. The average compressive modulus (from far-field strain) was about 10 Msi. The average post-impact static compressive strength was about 35.5 Ksi. The IM6 carbon sample had a strength of over 40 Ksi, more than 16 percent stronger than average. There was considerable scatter in the S-N data. However, the IM6 carbon samples clearly had the best fatigue response. The response of the other materials, while worse than IM6 carbon, could not be ranked definitively. The initial damage zones caused by the impact loading and damage growth from fatigue loading were similar for all five TTT reinforcing materials. The initial damage zones were circular and consisted of delaminations, matrix cracks and ply cracks. Post-impact fatigue loading caused delamination growth, ply cracking and fiber bundle failures, typically 45 deg from impact load direction. During the initial 97 percent of fatigue life, delaminations, ply cracks and fiber bundle failures primarily grew at and near the impact site. During the final 3 percent of life, damage grew rapidly transverse to the loading direction as a through-the-thickness transverse shear failure. The stress-strain response was typically linear during the initial 50 percent of life, and stiffness dropped about 20 percent during this period. During the next 47 percent of life, stiffness dropped about 34 percent, and the stress-strain response was no longer linear. The stiffness decreased about 23 percent during the final 3 percent of life. These trends were typical of all the materials tested. Therefore, by monitoring stiffness loss, fatigue failure could be accurately anticipated.

Serdinak, Thomas E.↗
DOE OSTI2024

Thermal field reconstruction and compressive sensing using proper orthogonal decomposition

Model order reduction allows critical information about sensor placement and experiment design to be distilled from raw fluid mechanics simulation data. In many cases, sensed information in conjunction with reduced order models can also be used to regenerate full field variables. In this paper, a proper orthogonal decomposition (POD) inferencing method is extended to the modeling and compressive sensing of temperature, a scalar field variable. The method is applied to a simulated, critically stable, incompressible flow over a heated cylinder (Re = 1000) with Prandtl number varying between 0.001 and 50. The model is trained on pressure and temperature data from simulations. Field reconstructions are then generated using data from selected sensors and the POD model. Finally, the reconstruction error is evaluated across all Prandtl numbers for different numbers of retained modes and sensors. The predicted trend of increasing reconstruction accuracy with decreasing Prandtl number is confirmed and a Prandtl number/sensor count error matrix is presented.

Matulis, John↗
DOE OSTI2020

Boundary-consistent B-spline filtering schemes and application to high-fidelity simulations of turbulence

A filtering operation, based on B-spline discretizations, is introduced to target weakly growing mesh-scale oscillations that can arise in high-fidelity turbulence simulations. This is a spectral regularization that can be described using the singular values of a banded matrix operator, with the filtering strength set by a scalar- or vector-valued penalty parameter. The penalty parameter can be specified though it can also be advantageously selected to minimize the generalized cross validation (GCV) measure of distance between the pre- and post-filtered solutions. Efficient algorithms are developed to compute both the scalar and vector penalty parameters. The B-spline filter has a sharper localization to high-wavenumber than compact or explicit filters of the same stencil width and is demonstrated for solutions of the Burgers' equation, decaying Burgers' turbulence, and compressible Navier–Stokes turbulent channel flow. Furthermore, these simulations confirm the scheme's numerical stability and ability to narrowly target the high wavenumber components of numerical solutions. An advantage over finite-difference filters is that these B-spline filters are stable on bounded domains and even preserve formal order of accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Emergent viscous hydrodynamics from a single quantum particle

We investigate an explicit example of how spatial decoherence can lead to hydrodynamic behavior in the late-time, long-wavelength regime of open quantum systems. We focus on the case of a single nonrelativistic quantum particle linearly coupled to a thermal bath of noninteracting harmonic oscillators at temperature T , a la Caldeira and Leggett. Taking advantage of decoherence in the position representation, we expand the reduced density matrix in powers of the off-diagonal spatial components, so that high-order terms are suppressed at late times. Truncating the resulting power series at second order leads to a set of dissipative transient hydrodynamic equations similar to the nonrelativistic limit of equations widely used in simulations of the quark-gluon plasma formed in ultrarelativistic heavy-ion collisions. Transport coefficients are directly determined by the damping constant γ , which quantifies the influence of the environment. The asymptotic limit of our hydrodynamic equations reduces to the celebrated Navier-Stokes equations for a compressible fluid in the presence of a drag force. Furthermore, our results shed new light on the onset of hydrodynamic behavior in open quantum systems where a system with few degrees of freedom is coupled to a large thermal environment.

Hydrodynamics↗
DOE OSTI2024

Brittle failure analysis and modeling of high-burnup PWR fuel cladding alloys

The aim of this research is the development of methods for predicting mechanical behavior and identification of limiting conditions to prevent brittle failure of high-burnup (HBU) pressure water reactor (PWR) fuel cladding alloys. A finite element (FE) model of the ring compression test (RCT) was created to analyze the failure behavior of zirconium-based alloys with radial hydrides during the RCT. An elastic-plastic material model describes the zirconium alloy. The stress-strain curve needed for the elastic-plastic material model was derived by inverse finite element analyses. Cohesive zone modeling is used to reproduce sudden load drops during RCT loading. Based on the failure mechanism in non-irradiated ZIRLO (R) claddings, a micro-mechanical model was developed that distinguishes between brittle failure along hydrides and ductile failure of the zirconium matrix. Two different cohesive laws representing these types of failure are present in the same cohesive interface. The key differences between these constitutive laws are the cohesive strength, the stress at which damage initiates, and the cohesive energy, which is the damage energy dissipated by the cohesive zone. Statistically generated matrix-hydride distributions were mapped onto the cohesive elements and simulations with focus on the first load drop were performed. Computational results are in good agreement with the RCT results conducted on high-burnup M5 (R) samples. It could be shown that crack initiation and propagation strongly depend on the specific configuration of hydrides and matrix material in the fracture area.

Simbruner, Kai↗
NASA NTRS1999

Tensegrity and mechanoregulation: from skeleton to cytoskeleton

OBJECTIVE: To elucidate how mechanical stresses that are applied to the whole organism are transmitted to individual cells and transduced into a biochemical response. DESIGN: In this article, we describe fundamental design principles that are used to stabilize the musculoskeletal system at many different size scales and show that these design features are embodied in one particular form of architecture that is known as tensegrity. RESULTS: Tensegrity structures are characterized by use of continuous tension and local compression; architecture, prestress (internal stress prior to application of external force), and triangulation play the most critical roles in terms of determining their mechanical stability. In living organisms, use of a hierarchy of tensegrity networks both optimizes structural efficiency and provides a mechanism to mechanically couple the parts with the whole: mechanical stresses applied at the macroscale result in structural rearrangements at the cell and molecular level. CONCLUSION: Due to use of tensegrity architecture, mechanical stress is concentrated and focused on signal transducing molecules that physically associate with cell surface molecules that anchor cells to extracellular matrix, such as integrins, and with load-bearing elements within the internal cytoskeleton and nucleus. Mechanochemical transduction may then proceed through local stress-dependent changes in molecular mechanics, thermodynamics, and kinetics within the cell. In this manner, the entire cellular response to stress may be orchestrated and tuned by altering the prestress in the cell, just as changing muscular tone can alter mechanical stability and structural coordination throughout the whole musculoskeletal system.

NASA Discipline Cell Biology↗
DOE OSTI2022

Recent Advances toward Efficient Calculation of Higher Nuclear Derivatives in Quantum Chemistry

In this article, we provide an overview of state-of-the-art techniques that are being developed for efficient calculation of second and higher nuclear derivatives of quantum mechanical (QM) energy. Calculations of nuclear Hessians and anharmonic terms incur high costs and memory and scale poorly with system size. Three emerging classes of methods—machine learning (ML), automatic differentiation (AD), and matrix completion (MC)—have demonstrated promise in overcoming these challenges. We illustrate studies that employ unsupervised ML methods to reduce the need for multiple Hessian calculations in dynamics simulations and those that utilize supervised ML to construct approximate potential energy surfaces and estimate Hessians and anharmonic terms at reduced cost. By extension, if electronic structure operations could be written in a manner similar to functions underlying ML methods, rapid differentiation or AD routines can be employed to inexpensively calculate higher arbitrary-order derivatives. While ML approaches are typically black-box, we describe methods such as compressed sensing (CS) and MC, which explicitly leverage problem-specific mathematical properties of higher derivatives such as sparsity and low-rank, to complete higher derivative information using only a small, incomplete sample. The three classes of methods facilitate reliable predictions of observables ranging from infrared spectra to thermal conductivity and constitute a promising way forward in accurately capturing otherwise intractable higher-order responses of QM energy to nuclear perturbations.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Scalable quantum simulations of scattering in scalar field theory on 120 qubits

Simulations of collisions of fundamental particles on a quantum computer are expected to have an exponential advantage over classical methods and promise to enhance searches for new physics. Furthermore, scattering in scalar field theory has been shown to be bounded-error quantum polynomial time (BQP) complete, making it a representative problem for which quantum computation is efficient. As a step toward large-scale quantum simulations of collision processes, scattering of wave packets in one-dimensional scalar field theory is simulated using 120 qubits of IBM’s Heron superconducting quantum computer ibm_fez. Variational circuits compressing vacuum preparation, wave packet initialization, and time evolution are determined using classical resources. By leveraging physical properties of states in the theory, such as symmetries and locality, the variational quantum algorithm constructs scalable circuits that can be used to simulate arbitrarily large system sizes. A new strategy is introduced to mitigate errors in quantum simulations, which enables the extraction of meaningful results from circuits with up to 4924 two-qubit gates and two-qubit gate depths of 103. The effect of interactions is clearly seen, and is found to be in agreement with classical matrix product state simulations. Finally, the developments that will be necessary to simulate high-energy inelastic collisions on a quantum computer are discussed.

quantum circuits↗