Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix multiplication”

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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Comparison of machine learning and electrical resistivity arrays to inverse modeling for locating and characterizing subsurface targets

Here, this study evaluates the performance of multiple machine learning (ML) algorithms and electrical resistivity (ER) arrays for inversion with comparison to a conventional Gauss-Newton numerical inversion method. Four different ML models and four arrays were used for the estimation of only six variables for locating and characterizing hypothetical subsurface targets. The combination of dipole-dipole with Multilayer Perceptron Neural Network (MLP-NN) had the highest accuracy. Evaluation showed that both MLP-NN and Gauss-Newton methods performed well for estimating the matrix resistivity while target resistivity accuracy was lower, and MLP-NN produced sharper contrast at target boundaries for the field and hypothetical data. Both methods exhibited comparable target characterization performance, whereas MLP-NN had increased accuracy compared to Gauss-Newton in prediction of target width and height, which was attributed to numerical smoothing present in the Gauss-Newton approach. MLP-NN was also applied to a field dataset acquired at U.S. DOE Hanford site.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Stability analysis of Monte Carlo burnup calculation based on PWR fuel rods

The existing Monte Carlo burnup algorithms may give some unstable results when using multiple burnable cells during burnup calculations. A very small time step is required to maintain the stability of the system, which pose some limits for the realistic applications. In this paper, we firstly repeat the instability issue by using a Monte Carlo code based on a typical PWR fuel rod model. Then, a theoretical stability analysis was carried out by using the perturbation theory in order to get a better understanding about the property of system stability. We determine a perturbation matrix which relates some main influence factors, such as the time step and the height of fuel rod, it can help people estimate the stability of burnup calculation system and the results are quantitatively consistent with the Monte Carlo code. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1995

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗
NASA NTRS1990

Unsteady transonic Navier-Stokes computations for an oscillating wing using single and multiple zones

Modern jet transports and maneuvering tactical fighters operating in the transonic regime often give rise to time-dependent fluid physics that interact with flexible structural components, e.g., vortical flow, shocks, and separation. Efficient computational fluid dynamic (CFD) methods are required to study such computationally intensive problems. In this work a numerical method is presented to address this problem. Time-dependent, compressible, Navier-Stokes equations are used to simulate unsteady transonic flow about a three-dimensional rigid wing undergoing a forced periodic motion in angle of attack. An efficient, implicit, diagonal algorithm is utilized because of its low operation count per time step compared to other methods that solve systems of block matrix equations. The formal time accuracy is addressed theoretically and demonstrated numerically by comparison of computational results with experimental data. A zonal grid approach, capable of treating complex geometries, is presented and its time accuracy is demonstrated by comparing a zonal computation with a single grid computation and experimental data.

Chaderjian, Neal M.↗
DOE OSTI2023

Effect of recycled fibers and shredded intermediates variation on the mechanical properties and energy absorption of fiber‐reinforced composite panels

Abstract The global composite industry generates large quantities of waste which mostly ends up in landfills due to a lack of established end‐use applications for multiple waste streams. The scrap from end‐of‐life (EoL) includes manufacturing waste such as dry chopped fiber tows, loose fibers, shredded fibers from fabric textile operations, cured/semi‐cured prepregs, and fully cured composite structure waste from aircraft, automobiles, wind blades, boats, and pressure vessels. In this work, different composite waste streams were reduced to shredded intermediates, followed by simple blending, and subjected to wet compression molding to produce composite panels. The panels/plaques were tested for mechanical properties (flexure and impact), fiber‐matrix wet‐out, and property bounds. It was found that wet‐compression molding was a viable and scale‐able approach to produce recycled panels from EoL composites shredded scrap. Furthermore, full‐scale size panels for use in truck bodies and intermodal shipping container flooring were manufactured and their impact resistance was tested using a drop weight impact test. They were tested both for high‐ and low‐velocity load. In the case of high‐velocity load, the average impact load was 14,673 N; the average absorbed energy was 101.6 J; the average elastic energy was 11.7 J and the impact resistance was 1065 J/m. In the case of the low‐velocity drop weight impact test, it was found that the average impact load was 7877.358 N; the average absorbed energy was 9.718 J; the average elastic energy was 9.14 J, and the impact resistance was 184.1 J/m. The shredded composite was shown to be a candidate material for the manufacture of truck bodies and intermodal containers’ flooring panels. Highlights By using shredded intermediates from different composite waste streams, it is possible to manufacture composite panels. Wet–compression process is an appropriate technique for manufacturing recycled fiber composite panels. Regardless of the source of scrap, the mechanical properties of the produced composite panels were improved. The recycling process technology can be transformed to commercial scale to produce full‐size transportation flooring panels. A product pathway is established in consideration of lower cost and improved recyclability.

Vaidya, Uday↗
DOE OSTI2024

MAGMA: Enabling exascale performance with accelerated BLAS and LAPACK for diverse GPU architectures

MAGMA (Matrix Algebra for GPU and Multicore Architectures) is a pivotal open-source library in the landscape of GPU-enabled dense and sparse linear algebra computations. With a repertoire of approximately 750 numerical routines across four precisions, MAGMA is deeply ingrained in the DOE software stack, playing a crucial role in high-performance computing. Notable projects such as ExaConstit, HiOP, MARBL, and STRUMPACK, among others, directly harness the capabilities of MAGMA. In addition, the MAGMA development team has been acknowledged multiple times for contributing to the vendors’ numerical software stacks. Looking back over the time of the Exascale Computing Project (ECP), we highlight how MAGMA has adapted to recent changes in modern HPC systems, especially the growing gap between CPU and GPU compute capabilities, as well as the introduction of low precision arithmetic in modern GPUs. We also describe MAGMA’s direct impact on several ECP projects. Maintaining portable performance across NVIDIA and AMD GPUs, and with current efforts toward supporting Intel GPUs, MAGMA ensures its adaptability and relevance in the ever-evolving landscape of GPU architectures.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1985

Characterization of erosion of metallic materials under cavitation attack in a mineral oil

Cavitation erosion and erosion rates of eight metallic materials representing three crystal structures were studied. The erosion experiments were conducted with a 20-kHz ultrasonic magnetostrictive oscillator in a viscous mineral oil. The erosion rates of the metals with an fcc matrix were 10 to 100 times higher than that of an hop-matrix titanium alloy. The erosion rates of iron and molybdenum, with bcc matrices, were higher than that of the titanium alloy but lower than those of those of the fcc materials. Studies with scanning electron microscopy indicated that the cavitation pits were initially formed at the grain boundaries and precipitates and that the pits formed at the junction of grain boundaries grew faster than the others. Transcrystalline craters formed by cavitation attack over the surface of grains and roughened the surfaces by multiple slip and twinning. Surface roughness measurements showed that the pits that formed over the grain boundaries deepened faster than pits. Computer analysis revealed that a geometric expression describes the nondimensional erosion curves during the time period 0.5 t (sub 0) t 2.5 t (sub 0), where t (sub 0) is the incubation period. The fcc metals had very short incubation periods; the titanium alloy had the longest incubation period.

Rao, B. C. S.↗
NASA NTRS2006

Exploring the Connection Between Sampling Problems in Bayesian Inference and Statistical Mechanics

The Bayesian and statistical mechanical communities often share the same objective in their work - estimating and integrating probability distribution functions (pdfs) describing stochastic systems, models or processes. Frequently, these pdfs are complex functions of random variables exhibiting multiple, well separated local minima. Conventional strategies for sampling such pdfs are inefficient, sometimes leading to an apparent non-ergodic behavior. Several recently developed techniques for handling this problem have been successfully applied in statistical mechanics. In the multicanonical and Wang-Landau Monte Carlo (MC) methods, the correct pdfs are recovered from uniform sampling of the parameter space by iteratively establishing proper weighting factors connecting these distributions. Trivial generalizations allow for sampling from any chosen pdf. The closely related transition matrix method relies on estimating transition probabilities between different states. All these methods proved to generate estimates of pdfs with high statistical accuracy. In another MC technique, parallel tempering, several random walks, each corresponding to a different value of a parameter (e.g. "temperature"), are generated and occasionally exchanged using the Metropolis criterion. This method can be considered as a statistically correct version of simulated annealing. An alternative approach is to represent the set of independent variables as a Hamiltonian system. Considerab!e progress has been made in understanding how to ensure that the system obeys the equipartition theorem or, equivalently, that coupling between the variables is correctly described. Then a host of techniques developed for dynamical systems can be used. Among them, probably the most powerful is the Adaptive Biasing Force method, in which thermodynamic integration and biased sampling are combined to yield very efficient estimates of pdfs. The third class of methods deals with transitions between states described by rate constants. These problems are isomorphic with chemical kinetics problems. Recently, several efficient techniques for this purpose have been developed based on the approach originally proposed by Gillespie. Although the utility of the techniques mentioned above for Bayesian problems has not been determined, further research along these lines is warranted

Pohorille, Andrew↗
NASA NTRS2016

Tribological Studies of Dynamic Thermal Seal Materials

Thermal seals are required on high-speed vehicles in many dynamic applications such as variable inlets in propulsion systems and control surfaces. These seals, often referred to as dynamic thermal seals, must not only mitigate inboard heat transfer, but must also exhibit sufficient durability when scrubbed against mating surfaces. For high-temperature high-speed vehicle applications, the mating surfaces are often made from thermal protection system (TPS) materials, which are typically rougher and more abrasive than TPS materials used at lower temperatures. The high-temperature TPS materials used can include non-ablative (e.g., lightweight porous oxides, ceramic matrix composites) andor ablative systems (e.g., phenolic systems). Due to the increased need for durable high-temperature dynamic seals, researchers working with the NASA Glenn Research Center embarked on an effort to (a) characterize the tribological performance of state-of-the-art thermal seal materials against a variety of TPS materials and (b) develop approaches for improved wear resistance. Tests were conducted using a recently upgraded high-temperature tribometer to assess wear resistance for a variety of tribopairs under multiple conditions. This data will begin to frame the challenges of using these materials and eventually permit an improved ability to design and implement these critical TPS components.

thermal protection↗
NASA NTRS1980

A new method for optimizing multiple flyby trajectories

A new procedure has been developed which minimizes total impulsive Delta V for multiple flyby trajectories with constraints on flyby parameters and maneuver times. The method involves solving a bounds-constrained parameter optimization problem with a Newton algorithm utilizing analytic first and second derivatives. Each trajectory segment connecting consecutive maneuver points is found by first targeting from the preceding maneuver point to the parameters of the upcoming flyby and then propagating the resulting trajectory to the next maneuver point. Multi-conic techniques are used for trajectory propagation and for computation of the state transition matrix. This procedure has successfully optimized Galileo satellite tours containing up to 11 flybys.

Damario, L. A.↗
NASA NTRS2000

Integrin activation and focal complex formation in cardiac hypertrophy

Cardiac hypertrophy is characterized by both remodeling of the extracellular matrix (ECM) and hypertrophic growth of the cardiocytes. Here we show increased expression and cytoskeletal association of the ECM proteins fibronectin and vitronectin in pressure-overloaded feline myocardium. These changes are accompanied by cytoskeletal binding and phosphorylation of focal adhesion kinase (FAK) at Tyr-397 and Tyr-925, c-Src at Tyr-416, recruitment of the adapter proteins p130(Cas), Shc, and Nck, and activation of the extracellular-regulated kinases ERK1/2. A synthetic peptide containing the Arg-Gly-Asp (RGD) motif of fibronectin and vitronectin was used to stimulate adult feline cardiomyocytes cultured on laminin or within a type-I collagen matrix. Whereas cardiocytes under both conditions showed RGD-stimulated ERK1/2 activation, only collagen-embedded cells exhibited cytoskeletal assembly of FAK, c-Src, Nck, and Shc. In RGD-stimulated collagen-embedded cells, FAK was phosphorylated only at Tyr-397 and c-Src association occurred without Tyr-416 phosphorylation and p130(Cas) association. Therefore, c-Src activation is not required for its cytoskeletal binding but may be important for additional phosphorylation of FAK. Overall, our study suggests that multiple signaling pathways originate in pressure-overloaded heart following integrin engagement with ECM proteins, including focal complex formation and ERK1/2 activation, and many of these pathways can be activated in cardiomyocytes via RGD-stimulated integrin activation.

NASA Discipline Cardiopulmonary↗
NASA NTRS2021

Nanostructure and compositional segregation of soft magnetic FeNi-based nanocomposites with multiple nanocrystalline phases

Soft magnetic metal amorphous nanocomposite alloys are produced through rapid solidification and thermal annealing yielding nanocrystals embedded within an amorphous precursor. Similar free energies in Co-rich and FeNi-based alloy systems result in multiple nanocrystalline phases being formed during devitrification. Studies of multi-phase crystallization processes have been reported for Co-rich alloys but relatively few have investigated FeNi-based systems. A detailed characterization of compositional partitioning and microstructure of an optimally annealed FeNi-based MANC (Fe70Ni30)80Nb4Si2B14 alloy is presented through complementary high-resolution transmission electron microscopy (HRTEM) and atom probe tomography (APT). HRTEM demonstrates orientation relationships between FCC and BCC nanocrystals, suggesting heterogeneous nucleation of nanocrystals in the amorphous matrix or a cooperative mechanism of nucleation between BCC and FCC nanocrystallites. APT results show evidence for (i) the segregation of Fe and Ni between nanocrystals of different phases, (ii) B partitioning to the amorphous phase, and (iii) an Nb-enriched shell surrounding nanocrystals.

soft magnetic materials↗
DOE OSTI2024

O-GlcNAc transferase regulates collagen deposition and fibrosis resolution in idiopathic pulmonary fibrosis

Idiopathic pulmonary fibrosis (IPF) is a chronic pulmonary disease that is characterized by an excessive accumulation of extracellular matrix (ECM) proteins (e.g. collagens) in the parenchyma, which ultimately leads to respiratory failure and death. While current therapies exist to slow the progression, no therapies are available to resolve fibrosis. We characterized the O-linked N-Acetylglucosamine (O-GlcNAc) transferase (OGT)/O-GlcNAc axis in IPF using single-cell RNA-sequencing (scRNA-seq) data and human lung sections and isolated fibroblasts from IPF and non-IPF donors. The underlying mechanism(s) of IPF were further investigated using multiple experimental models to modulate collagen expression and accumulation by genetically and pharmacologically targeting OGT. Furthermore, we hone in on the transforming growth factor-beta (TGF-β) effector molecule, Smad3, by co-expressing it with OGT to determine if it is modified and its subsequent effect on Smad3 activation. We found that OGT and O-GlcNAc levels are upregulated in patients with IPF compared to non-IPF. We report that the OGT regulates collagen deposition and fibrosis resolution, which is an evolutionarily conserved process demonstrated across multiple species. Co-expression of OGT and Smad3 showed that Smad3 is O-GlcNAc modified. Blocking OGT activity resulted in decreased phosphorylation at Ser-423/425 of Smad3 attenuating the effects of TGF-β1 induced collagen expression/deposition. OGT inhibition or knockdown successfully blocked and reversed collagen expression and accumulation, respectively. Smad3 is discovered to be a substrate of OGT and its O-GlcNAc modification(s) directly affects its phosphorylation state. These data identify OGT as a potential target in pulmonary fibrosis resolution, as well as other diseases that might have aberrant ECM/collagen accumulation.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1991

Generalized covariance analysis for partially autonomous deep space missions

A new covariance analysis method is presented that is suitable for the evaluation of multiple impulsive controllers acting on some stochastic process x. The method accommodates batch and sequential estimators with equal ease and accounts for time-delay effects in a natural manner. The formalism is developed in terms of a generalized state vector that is formed from the system state vector x, augmented by various fixed epoch estimates, and a data vector formed from discrete time observations of the system. Recursions are developed for time transition, measurement incorporation, and impulsive control updating of the generalized covariance matrix. Means of limiting the dimensional growth of the generalized state vector via the processes of estimator epoch adjustment and measurement vector deflation are described and the application of numerically stable matrix factorization methods to the generalized covariance recursions is outlined. The method is applied to the Magellan spacecraft to demonstrate the capability of ground-based optimal estimation and control of gyro/star scanner misalignment.

Boone, Jack N.↗
NASA NTRS2005

Application of a Constant Gain Extended Kalman Filter for In-Flight Estimation of Aircraft Engine Performance Parameters

An approach based on the Constant Gain Extended Kalman Filter (CGEKF) technique is investigated for the in-flight estimation of non-measurable performance parameters of aircraft engines. Performance parameters, such as thrust and stall margins, provide crucial information for operating an aircraft engine in a safe and efficient manner, but they cannot be directly measured during flight. A technique to accurately estimate these parameters is, therefore, essential for further enhancement of engine operation. In this paper, a CGEKF is developed by combining an on-board engine model and a single Kalman gain matrix. In order to make the on-board engine model adaptive to the real engine s performance variations due to degradation or anomalies, the CGEKF is designed with the ability to adjust its performance through the adjustment of artificial parameters called tuning parameters. With this design approach, the CGEKF can maintain accurate estimation performance when it is applied to aircraft engines at offnominal conditions. The performance of the CGEKF is evaluated in a simulation environment using numerous component degradation and fault scenarios at multiple operating conditions.

Kobayashi, Takahisa↗
NASA NTRS1990

Computing Scattering Matrices For Circular Waveguides

Scattering Matrix Program for Circular Waveguide Junctions computes scattering matrix for series of circular waveguide sections. Sections must have same axis, but radius and length of each section completely arbitrary. Devices analyzed include simple waveguide step discontinuity like that used in a dual-mode horn, stepped matching section, or corrugated waveguide section with constant or varying slot depth. Certain types of corrugated horns also analyzed with program. Mathematical model used in program accurately predicts reflection and transmission characteristics of such devices, taking into account excitation of modes of higher order as well as multiple reflections and energy stored at each discontinuity. Written in FORTRAN 77.

Hoppe, Daniel J.↗
NASA NTRS2012

The Spatial Distribution of Organic Matter and Mineralogical Relationships in Carbonaceous Chondrites

Organic matter present within primitive carbonaceous meteorites represents the complex conglomeration of species formed in a variety of physically and temporally distinct environments including circumstellar space, the interstellar medium, the Solar Nebula & Jovian sub-nebulae and asteroids. In each case, multiple chemical pathways would have been available for the synthesis of organic molecules. Consequently these meteorites constitute a unique record of organic chemical evolution in the Universe and one of the biggest challenges in organic cosmochemistry has been in deciphering this record. While bulk chemical analysis has provided a detailed description of the range and diversity of organic species present in carbonaceous chondrites, there is virtually no hard experimental data as to how these species are spatially distributed and their relationship to the host mineral matrix, (with one exception). The distribution of organic phases is nevertheless critical to understanding parent body processes. The CM and CI chondrites all display evidence of low temperature (< 350K) interaction with aqueous fluids, which based on O isotope data, flowed along thermal gradients within the respective parent bodies. This pervasive aqueous alteration may have led to aqueous geochromatographic separation of organics and synthesis of new organics coupled to aqueous mineral alteration. To address such issues we have applied the technique of microprobe two-step laser desorption / photoionization mass spectrometry (L2MS) to map in situ the spatial distribution of a broad range of organic species at the micron scale in the freshly exposed matrices of the Bells, Tagish Lake and Murchison (CM2) carbonaceous chondrites.

Clemett, S. J.↗