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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 451 records · Page 25

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Revisiting the Hail Radar Reflectivity–Kinetic Energy Flux Relation by Combining T-Matrix and Discrete Dipole Approximation Calculations to Size Distribution Observations

The retrieval of hail kinetic energy with weather radars or its simulation in numerical models is challenging because of the shape complexity and variable density of hailstones. We combine 3D scans of individual hailstones with measurements of the particle size distributions (PSD) and T-matrix calculations to understand how hail reflectivity Z changes when approximating hailstones as spheroids, as compared to the realistic shapes obtained by 3D scanning technology. Additionally, recent terminal velocity relations are used to compare Z to the hail kinetic energy flux E ˙ . We parameterize the hail backscattering cross sections at L, S, C, and X bands as a function of size between 0.5 and 5.0 cm, matching the range of the observed PSDs. The scattering calculations use the T-matrix method for size parameters below 1.0 and the discrete dipole approximation (DDA) method otherwise. The DDA calculations are done for 48 digital models of realistic hailstones of sizes between 1 and 5 cm. The DDA cross sections are calculated for multiple orientations and averaged assuming a fully random orientation distribution to provide a single value per hailstone. The T-matrix reflectivity assuming solid ice spheres presents negligible differences to DDA results for size parameters below 1.0. Therefore, T matrix was used to fill in the gaps left by the DDA calculations. The results are mapped to the same size bins of the observed PSDs, allowing the calculation of the radar reflectivity. This is then correlated to E ˙ , allowing a potential improvement of past retrieval methods of E ˙ from Z in multiple wavelengths.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Driving rapid atomic order in MnAl via low-magnitude magnetic field annealing

Application of a mild (60 mT), uniform magnetic field during short-term thermal treatment of kinetically retained, atomically disordered (paramagnetic) ε-MnAl was found to deliver a significant ~50 % increase in the formation of L1 0 atomically ordered (ferromagnetic) τ-MnAl product phase, compared to that produced by conventional (i.e., zero-field) annealing under identical thermal conditions. The magnetic field, applied in a passive closed-circuit configuration during annealing, induced significant changes in the structural, magnetic, and phase evolution of the material. Computational results based on electronic structure calculations demonstrate that the effective magnetic susceptibility of τ-MnAl is sensitive to the orientation, rather than the magnitude, of an applied magnetic field in the vicinity of the Curie temperature. The uniaxial magnetocrystalline anisotropy of the L1 0 structure is proposed to act as a filter for selective propagation of the population of τ-MnAl variants that are favorably aligned with the applied field. In this manner, crystallographic “gridlock” is alleviated that would otherwise arise from the coexistence of multiple, energetically equivalent τ-phase variants within the parent ε-phase matrix. These results confirm that static, low-magnitude magnetic field annealing is able to accelerate L1 0 atomic ordering in the MnAl system and likely can exert similar influences in relevant magnetic systems, facilitating efficient tailoring of structure-sensitive magnetic properties for the manufacture of magnetic materials.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Multi-objective Bayesian alloy design using multi-task Gaussian processes

In design applications, correlations among material properties (such as the tendency for stronger materials to be less ductile) are often neglected. This approach is echoed in multi-objective optimization techniques which treat each performance characteristic as an independent objective, aiming to optimize scalar functions and find optimal Pareto fronts. However, this overlooks the statistical relationships between performance characteristics inherent in a material system. To address this, we propose the use of Bayesian optimization, a highly efficient black-box optimization algorithm known for constructing Gaussian processes (GPs) – uncorrelated surrogates - to model objective functions. Rather than evaluating multiple GPs for each objective function separately, we argue for a shift towards jointly modeling these objective functions, considering their statistical correlations. This integrated approach utilizes naturally occurring relationships among material properties, providing additional information to enhance the performance of the design framework. This requires the replacement of multiple independent GPs with a single multi-task GP, employing a correlation matrix to construct a multi-task kernel function, wherein each task corresponds to a single objective function. Here, we anticipate this refined methodology will better leverage material correlations, improving design optimization results.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Pareto-optimal target definition for multi-axis random vibration testing

In random vibration testing with multiple control channels, existing control laws require specification of a complete spectral density matrix at each control frequency. Spectral density matrices include autospectral densities on the diagonal and cross-spectral densities on the off-diagonal. In practice, the off-diagonal terms are often unknown, and recent vibration testing research has focused on fixing the diagonal and specifying the off-diagonal to minimize the required control energy, subject to a constraint that the target matrix is positive semidefinite. This paper shows that, even with a fixed diagonal, off-diagonal terms strongly affect control residuals. This overlooked effect occurs in both square and rectangular systems. By jointly considering input energy and control residuals, open-loop inputs are derived directly from the diagonal without specifying the off-diagonal terms. Vibration targets that can be used in closed-loop control are then derived using the optimal inputs, with positive semidefinite constraints applied during the derivation. The result is a set of Pareto-optimal control solutions. For each solution in the set, any other possible solution produces greater control error, greater input energy, or both. A balanced solution is selected automatically, though others can be chosen based on test needs. Simulations and experiments show that the proposed method outperforms state-of-the-art energy-minimizing approaches, achieving significant reductions in both control error and input energy.

Autospectral density↗

DOE OSTI2020

Faster Johnson–Lindenstrauss transforms via Kronecker products

The Kronecker product is an important matrix operation with a wide range of applications in signal processing, graph theory, quantum computing and deep learning. In this work, we introduce a generalization of the fast Johnson–Lindenstrauss projection for embedding vectors with Kronecker product structure, the Kronecker fast Johnson–Lindenstrauss transform (KFJLT). The KFJLT reduces the embedding cost by an exponential factor of the standard fast Johnson–Lindenstrauss transform’s cost when applied to vectors with Kronecker structure, by avoiding explicitly forming the full Kronecker products. Here, we prove that this computational gain comes with only a small price in embedding power: consider a finite set of $p$ points in a tensor product of $d$ constituent Euclidean spaces $\bigotimes _{k=d}^{1}{\mathbb{R}}^{n_k}$, and let $N = \prod _{k=1}^{d}n_k$. With high probability, a random KFJLT matrix of dimension $m \times N$ embeds the set of points up to multiplicative distortion $(1\pm \varepsilon )$ provided $m \gtrsim \varepsilon ^{-2} \, \log ^{2d - 1} (p) \, \log N$. We conclude by describing a direct application of the KFJLT to the efficient solution of large-scale Kronecker-structured least squares problems for fitting the CP tensor decomposition.

Kronecker structure↗

DOE OSTI2022

𝛼-cluster structure in 19 F and 19 Ne in resonant scattering

This work investigates the nuclear structure of 19 F and 19 Ne, which is important for understanding α clustering in the A = 20 mass region and for astrophysical applications. The only high-resolution, broad angular- and energy-range study of the 19 F resonance structure in α + 15 N scattering was published over 60 years ago, when a detailed analysis of complex excitation functions with overlapping resonances and multiple decay channels was not feasible. We have performed a modern R-matrix analysis of these data to assign spins and determine resonance parameters for levels in 19 F up to an excitation energy of 8.2 MeV. Our R-matrix parameters were successfully tested by fitting recent α + 15 N data obtained with the Thick Target Inverse Kinematics (TTIK) method at 180°. The new 19 F parameters were then used to fit TTIK data for α + 15 O, the mirror resonant reaction. In conclusion, comparison of these isobaric mirror reactions provides valuable insight into the underlying nuclear structure.

6 ≤ A ≤ 19↗

DOE OSTI2026

A Novel Low-Profile High-Efficiency Three-Phase Matrix Transformer

High step-down isolated DC-DC conversion from an 800 V DC bus to low-voltage, high-current outputs is required in automotive auxiliary converters and data center power supplies. In such applications, conventional transformer-based converters require large turns ratios, which increase winding resistance, leakage inductance, and magnetic height. This paper proposes a novel low-profile three-phase matrix transformer that realizes a large effective voltage ratio through flux division among multiple secondary legs, without increasing the physical turns count of each winding. As a result, the proposed structure reduces copper usage and transformer height while preserving the voltage conversion capability of a conventional three-phase transformer. Finite element analysis shows that the proposed design reduces magnetic height by 27%, ferrite volume by 34%, and copper volume by 28%. Circuit-level simulations of an 800 V/12 V,3 kW CLLLC dual-active-bridge converter further show that the lower winding resistance reduces total system loss by 91% and increases DC-DC efficiency from 82.6% to 97.6% at 3 kW output.

Inoue, Shuntaro [ORNL] (ORCID:0000000262637627)↗

DOE OSTI2020

Mixed Material Scintillator Systems for Neutron Detection

Scintillators are a fundamental technology that enables radiation detection. Their use has historically enabled the detection of incoming gamma radiation, with recent advances in plastic scintillator development enabling detection and discrimination of neutron and gamma radiation through pulse shape discrimination. However, scintillators have so far consisted of a single material, typically one primary dye and one secondary dye embedded in a plastic matrix. Based on the simulation studies described previously, mixed material scintillator systems with multiple secondary dyes arranged in alternating, periodic micro-sized regions opens the door for scintillators with a wide variety of different capability, such as directional sensing, positional reconstruction, particle identification, and pulse shape discrimination. However, fabrication of a mixed material scintillator system by conventional manufacturing methods is difficult to manufacture and scale, creating a barrier to their wide spread implementation. On the other hand, additive manufacturing additive manufacturing presents itself as the optimal method for fabricating mixed material scintillator systems based on its flexibility, ability to fabricate fine features, and potential for multimaterial fabrication.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1976

Efficient reanalysis of structures by a direct modification method

A procedure for the local stiffness modifications of large structures is described. It enables structural modifications without an a priori definition of the changes in the original structure and without loss of efficiency due to multiple loading conditions. The solution procedure, implemented in NASTRAN, involved the decomposed stiffness matrix and the displacement vectors of the original structure. It solves the modified structure exactly, irrespective of the magnitude of the stiffness changes. In order to investigate the efficiency of the present procedure and to test its applicability within a design environment, several real and large structures were solved. The results of the efficiency studies indicate that the break-even point of the procedure varies between 8% and 60% stiffness modifications, depending upon the structure's characteristics and the options employed.

Raibstein, A. I.↗

NASA NTRS1984

Future remote sensing instruments and systems

Scientific requirements for terrestrial remote-sensing satellites in the 1990's and instrument packages to meet them are discussed, summarizing the results of a preliminary study conducted by the NASA System Z Science and Mission Requirements Working Group. Major questions regarding atmospheric chemistry and circulation, global climate, ocean dynamics, biogeochemical cycles, ice-budget and hydrological cycles, continental geology, biomass dynamics, and land cover and use are listed along with the observational frequency and resolution necessary to answer them, and an instrument use matrix is presented. Large essentially permanent man-tended arrays of instruments serving multiple functions and linked to the users via advanced information-processing networks and (possibly including expert systems) are recommended, and a basic payload and two supplemental packages are outlined.

Broome, D. R., Jr.↗

NASA NTRS1992

The link evaluation terminal for the advanced communications technology satellite experiments program

The experimental NASA satellite, Advanced Communications Technology Satellite (ACTS), introduces new technology for high throughput 30 to 20 GHz satellite services. Contained in a single communication payload is both a regenerative TDMA system and multiple 800 MHz 'bent pipe' channels routed to spot beams by a switch matrix. While only one mode of operation is typical during any experiment, both modes can operate simultaneously with reduced capability due to sharing of the transponder. NASA-Lewis instituted a ground terminal development program in anticipation of the satellite launch to verify the performance of the switch matrix mode of operations. Specific functions are built into the ground terminal to evaluate rain fade compensation with uplink power control and to monitor satellite transponder performance with bit error rate measurements. These functions were the genesis of the ground terminal's name, Link Evaluation Terminal, often referred to as LET. Connectors are included in LET that allow independent experimenters to run unique modulation or network experiments through ACTS using only the RF transmit and receive portions of LET. Test data indicate that LET will be able to verify important parts of ACTS technology and provide independent experimenters with a useful ground terminal. Lab measurements of major subsystems integrated into LET are presented. Bit error rate is measured with LET in an internal loopback mode.

May, Brian D.↗

NASA NTRS1992

The Link Evaluation Terminal for the Advanced Communications Technology Satellite Experiments Program

May, Brian D.↗

NASA NTRS2000

Adjoint Techniques for Topology Optimization of Structures Under Damage Conditions

The objective of this cooperative agreement was to seek computationally efficient ways to optimize aerospace structures subject to damage tolerance criteria. Optimization was to involve sizing as well as topology optimization. The work was done in collaboration with Steve Scotti, Chauncey Wu and Joanne Walsh at the NASA Langley Research Center. Computation of constraint sensitivity is normally the most time-consuming step of an optimization procedure. The cooperative work first focused on this issue and implemented the adjoint method of sensitivity computation (Haftka and Gurdal, 1992) in an optimization code (runstream) written in Engineering Analysis Language (EAL). The method was implemented both for bar and plate elements including buckling sensitivity for the latter. Lumping of constraints was investigated as a means to reduce the computational cost. Adjoint sensitivity computation was developed and implemented for lumped stress and buckling constraints. Cost of the direct method and the adjoint method was compared for various structures with and without lumping. The results were reported in two papers (Akgun et al., 1998a and 1999). It is desirable to optimize topology of an aerospace structure subject to a large number of damage scenarios so that a damage tolerant structure is obtained. Including damage scenarios in the design procedure is critical in order to avoid large mass penalties at later stages (Haftka et al., 1983). A common method for topology optimization is that of compliance minimization (Bendsoe, 1995) which has not been used for damage tolerant design. In the present work, topology optimization is treated as a conventional problem aiming to minimize the weight subject to stress constraints. Multiple damage configurations (scenarios) are considered. Each configuration has its own structural stiffness matrix and, normally, requires factoring of the matrix and solution of the system of equations. Damage that is expected to be tolerated is local and represents a small change in the stiffness matrix compared to the baseline (undamaged) structure. The exact solution to a slightly modified set of equations can be obtained from the baseline solution economically without actually solving the modified system.. Shennan-Morrison-Woodbury (SMW) formulas are matrix update formulas that allow this (Akgun et al., 1998b). SMW formulas were therefore used here to compute adjoint displacements for sensitivity computation and structural displacements in damaged configurations.

Akgun, Mehmet A.↗

NASA NTRS1999

Optimization of Aerospace Structure Subject to Damage Tolerance Criteria

The objective of this cooperative agreement was to seek computationally efficient ways to optimize aerospace structures subject to damage tolerance criteria. Optimization was to involve sizing as well as topology optimization. The work was done in collaboration with Steve Scotti, Chauncey Wu and Joanne Walsh at the NASA Langley Research Center. Computation of constraint sensitivity is normally the most time-consuming step of an optimization procedure. The cooperative work first focused on this issue and implemented the adjoint method of sensitivity computation in an optimization code (runstream) written in Engineering Analysis Language (EAL). The method was implemented both for bar and plate elements including buckling sensitivity for the latter. Lumping of constraints was investigated as a means to reduce the computational cost. Adjoint sensitivity computation was developed and implemented for lumped stress and buckling constraints. Cost of the direct method and the adjoint method was compared for various structures with and without lumping. The results were reported in two papers. It is desirable to optimize topology of an aerospace structure subject to a large number of damage scenarios so that a damage tolerant structure is obtained. Including damage scenarios in the design procedure is critical in order to avoid large mass penalties at later stages. A common method for topology optimization is that of compliance minimization which has not been used for damage tolerant design. In the present work, topology optimization is treated as a conventional problem aiming to minimize the weight subject to stress constraints. Multiple damage configurations (scenarios) are considered. Each configuration has its own structural stiffness matrix and, normally, requires factoring of the matrix and solution of the system of equations. Damage that is expected to be tolerated is local and represents a small change in the stiffness matrix compared to the baseline (undamaged) structure. The exact solution to a slightly modified set of equations can be obtained from the baseline solution economically without actually solving the modified system. Sherrnan-Morrison-Woodbury (SMW) formulas are matrix update formulas that allow this. SMW formulas were therefore used here to compute adjoint displacements for sensitivity computation and structural displacements in damaged configurations.

Akgun, Mehmet A.↗

NASA NTRS2004

Random Matrix Approach to Quantum Adiabatic Evolution Algorithms

We analyze the power of quantum adiabatic evolution algorithms (Q-QA) for solving random NP-hard optimization problems within a theoretical framework based on the random matrix theory (RMT). We present two types of the driven RMT models. In the first model, the driving Hamiltonian is represented by Brownian motion in the matrix space. We use the Brownian motion model to obtain a description of multiple avoided crossing phenomena. We show that the failure mechanism of the QAA is due to the interaction of the ground state with the "cloud" formed by all the excited states, confirming that in the driven RMT models. the Landau-Zener mechanism of dissipation is not important. We show that the QAEA has a finite probability of success in a certain range of parameters. implying the polynomial complexity of the algorithm. The second model corresponds to the standard QAEA with the problem Hamiltonian taken from the Gaussian Unitary RMT ensemble (GUE). We show that the level dynamics in this model can be mapped onto the dynamics in the Brownian motion model. However, the driven RMT model always leads to the exponential complexity of the algorithm due to the presence of the long-range intertemporal correlations of the eigenvalues. Our results indicate that the weakness of effective transitions is the leading effect that can make the Markovian type QAEA successful.

Boulatov, Alexei↗

NASA NTRS2016

Incorporation of Damage and Failure into an Orthotropic Elasto-Plastic Three-Dimensional Model with Tabulated Input Suitable for Use in Composite Impact Problems

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 LS-DYNA(Registered Trademark) is under development. In particular, the material model, which is being implemented as MAT 213 into a tailored version of LS-DYNA being jointly developed by the FAA and NASA, incorporates both plasticity and damage within 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), and is able to analyze the response of composites composed with a variety of fiber architectures. 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, in composites it has been found that loading in one direction can lead to damage in multiple coordinate directions. To account for this phenomena, 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. The onset of material failure, and thus element deletion, is being developed to be a function of the stresses and plastic strains in the various coordinate directions. Systematic procedures are being developed to generate the required input parameters based on the results of experimental tests.

Finite Element Method↗

NASA NTRS2022

Quantifying Cybersecurity Risk for NASA Missions

An end-to-end cyber risk assessment process is presented that is based on the combination of guidelines from the National Institute of Standards & Technology (NIST), the standard 5x5 risk matrix, and quantitative methods for generating loss exceedance curves. The NIST guidelines provide a framework for cyber risk assessment, and the standard 5x5 matrix is widely used across the industry for the representation of risk across multiple disciplines. Loss exceedance curves are a means of quantitatively assessing the loss that occurs due to a given risk profile. Combining these different techniques enables us to follow the guidelines, adhere to standard 5x5 risk management practices and develop quantitative metrics simultaneously. Our quantification process is based on the consideration of the NASA and JPL Cost Risk assessment modeling techniques as we define the cost associated with the cybersecurity risk profile of a mission as a function of the mission cost.

Miller, Robert L.↗

NASA NTRS2021

Polarimetric calibration of the Multi-Angle Imager for Aerosols (MAIA)

NASA’s Multi-Angle Imager for Aerosols (MAIA) mission, under development at the Jet Propulsion Laboratory, is designed to study the adverse health effects of different types of particulate air pollution. Planned for launch in 2023 for a 3-year mission, the MAIA satellite instrument will focus on a selected set of metropolitan target areas, where air quality monitors and health data are available. Aerosol concentration and speciation are inferred from multi-angle measurements of backscattered sunlight in 14 spectral bands from 350-2200 nm, with bands near 442, 645 and 1040 nm measuring the degree (DoLP) and angle of linear polarization (AoLP) in addition to radiance. The pushbroom camera has a ~240-km cross-track field of view with a nadir resolution of ~200 m, and is mounted onto a biaxial gimbal to provide along-track view angles within ±60°, to extend the field of regard to ±48°, and to view the instrument’s onboard calibrator (OBC) and dark target. The OBC consists of a sunlit transmissive diffuser, followed by 12 polarizers at different orientations. MAIA’s polarimetry is implemented using miniature wire grid polarizers on the focal plane array, and dual photoelastic modulators (PEMs) and achromatic quarter-wave plates to rapidly rotate the polarization. The resulting ~26-Hz intensity modulation encodes the linearly polarized and total radiance in each pixel, leaving the DoLP and AoLP insensitive to gain calibration. We report on the polarimetric calibration of the MAIA camera using a vacuum-compatible polarization state generator, consisting of a 1600W Xenon lamp, 12-inch integrating sphere, and rotating high-extinction polarizer. Mueller-matrix-based calibration coefficients for each detector pixel are derived from measurements at multiple polarizer angles, and are used to correct the measurements for instrumental polarization aberrations. Prior to flight, the calibrated MAIA camera is panned across the OBC to characterize its output, using uniform illumination with an irradiance similar to the Sun.

Werne, Thomas↗