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 595 records · Page 33

All sources NASA NTRS DOE OSTI

NASA NTRS1994

Optimization of Residual Stresses in MMC's Using Compensating/Compliant Interfacial Layers. Part 2: OPTCOMP User's Guide

A user's guide for the computer program OPTCOMP is presented in this report. This program provides a capability to optimize the fabrication or service-induced residual stresses in uni-directional metal matrix composites subjected to combined thermo-mechanical axisymmetric loading using compensating or compliant layers at the fiber/matrix interface. The user specifies the architecture and the initial material parameters of the interfacial region, which can be either elastic or elastoplastic, and defines the design variables, together with the objective function, the associated constraints and the loading history through a user-friendly data input interface. The optimization procedure is based on an efficient solution methodology for the elastoplastic response of an arbitrarily layered multiple concentric cylinder model that is coupled to the commercial optimization package DOT. The solution methodology for the arbitrarily layered cylinder is based on the local-global stiffness matrix formulation and Mendelson's iterative technique of successive elastic solutions developed for elastoplastic boundary-value problems. The optimization algorithm employed in DOT is based on the method of feasible directions.

Pindera, Marek-Jerzy↗

DOE OSTI2021

Fault isolation and fault-tolerant control for Takagi-Sugeno fuzzy time-varying delay stochastic distribution systems

A fault isolation, estimation, and fault-tolerant control (FTC) scheme for nonlinear time-varying delay stochastic distribution control systems was presented in this paper. The Takagi-Sugeno fuzzy model was adopted to approach the nonlinear dynamics of time-varying delay systems. According to the output equivalence principle and Laplace transformation, an augmented state vector was given to solve the time-varying delay problem. When multiple actuator faults and interference occur simultaneously, fault detection, isolation and fault estimation was designed to obtained the fault information. To decouple faults and obtain the value and location information of the fault, the system was separated into two parts through the designed multiple conversion matrices, in which one subsystem was only affected by one actuator fault. This has simplified the design of fault isolation and estimation. A adaptive observer for fault estimation was given. Then, fault information such as the time, location, and size was determined. The observer gain matrices were calculated using linear matrix inequality (LMI). When a fault was detected and diagnosed, a FTC algorithm was devised using the proportional-integral control scheme to compensate the fault as much as possible. It has been shown that even if multiple faults actuator occurred simultaneously, the FTC controller still ensured the output probability density function of the system traced the desired probability density function when a fault occurred. Finally, the expected results were obtained through the simulation example, which confirmed the effectiveness of the method.

42 ENGINEERING↗

NASA NTRS2006

An Analysis of Nondestructive Evaluation Techniques for Polymer Matrix Composite Sandwich Materials

Structural sandwich materials composed of triaxially braided polymer matrix composite material face sheets sandwiching a foam core are being utilized for applications including aerospace components and recreational equipment. Since full scale components are being made from these sandwich materials, it is necessary to develop proper inspection practices for their manufacture and in-field use. Specifically, nondestructive evaluation (NDE) techniques need to be investigated for analysis of components made from these materials. Hockey blades made from sandwich materials and a flat sandwich sample were examined with multiple NDE techniques including thermographic, radiographic, and shearographic methods to investigate damage induced in the blades and flat panel components. Hockey blades used during actual play and a flat polymer matrix composite sandwich sample with damage inserted into the foam core were investigated with each technique. NDE images from the samples were presented and discussed. Structural elements within each blade were observed with radiographic imaging. Damaged regions and some structural elements of the hockey blades were identified with thermographic imaging. Structural elements, damaged regions, and other material variations were detected in the hockey blades with shearography. Each technique s advantages and disadvantages were considered in making recommendations for inspection of components made from these types of materials.

Cosgriff, Laura M.↗

DOE OSTI2025

Machine Learning–Guided Boolean Matrix Inference for Real-Time O-RAN Conflict Detection

Open Radio Access Networks (O-RAN) are emerging, software-driven cellular architectures that promote flexibility by enabling components from different vendors to interoperate. Multiple control applications called xApps can independently adjust network parameters in near real time, often without awareness of each other's actions. This creates a system highly prone to unintended conflicts and performance degradation due to the inherent complexity of such openness. To model such systems and ultimately prevent or mitigate xApp conflicts, it is essential to understand the dynamic relationships between xApps (A), the control parameters they adjust (P), and the resulting KPI responses (K). While the mappings from A to P and from K to A can often be derived from xApp specifications, the relationship from P to K is typically hidden within the system’s dynamics and must be inferred from observed data. We propose a novel data-driven Boolean inference framework that uncovers the hidden P?K dependencies using machine learning and interpretable rule induction. Continuous parameters and KPIs are first binarized using decision tree classifiers, and a binary influence matrix L is then inferred by solving Boolean matrix equations over time. This compact representation improves interpretability and enables real-time tracking of dynamically evolving parameter-KPI dependencies. We demonstrate the effectiveness of our method in a realistic mobile handover scenario, where it accurately recovers the underlying logic and enables proactive conflict detection.

42 - ENGINEERING↗

DOE OSTI2020

Nuclear Reaction Data Evaluations for Nuclear Applications Physics and modeling [Slides]

This presentation provides an evaluation workflow scheme for the project. It also discusses nuclear data libraries, including an overview of worldwide nuclear data libraries and current funding agencies. The presentation also discusses the theoretical framework commonly used in the evaluated work, including the R-matrix theory (≈ up to a few keVs for heavy nuclei) as well as pre-equilibrium (PE) and Hauser-Feshbach (HF) theory (⪆1 MeV for medium/heavy nuclei). Inclusion of experimental effects are also discussed. Examples of evaluated work, including theory, experiments, and applications are discussed, specifically n+ 28,29,30 Si evaluation (direct capture component), n+ 53 Cr evaluation (strong multiple scattering effects), n+ 239 Pu(n,γf) reaction (fluctuations in neutron multiplicities), and α+ 17,18 O evaluations (emerging neutron spectra).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Quantum-classical embedding via ghost Gutzwiller approximation for enhanced simulations of correlated electron systems

Simulating correlated materials on present-day quantum hardware remains challenging due to limited quantum resources. Quantum embedding methods offer a promising route by reducing computational complexity through the mapping of bulk systems onto effective impurity models, allowing more feasible simulations on pre- and early-fault-tolerant quantum devices. Here, this work develops a quantum-classical embedding framework based on the ghost Gutzwiller approximation to enable quantum-enhanced simulations of ground-state properties and spectral functions of correlated electron systems. Circuit complexity is analyzed using an adaptive variational quantum algorithm on a statevector simulator, applied to the infinite-dimensional Hubbard model with increasing ghost mode numbers from 3 to 5, resulting in circuit depths growing from 16 to 104. Noise effects are examined using a realistic error model, revealing significant impact on the spectral weight of the Hubbard bands. To mitigate these effects, the Iceberg quantum error detection code is employed, achieving up to 40% error reduction in simulations. Finally, the accuracy of the density matrix estimation and the derived spectral function is benchmarked on IBM and Quantinuum quantum hardware, featuring distinct qubit-connectivity and employing multiple levels of error mitigation techniques.

Chen, I-Chi [Ames Laboratory (AMES), Ames, IA (Uni↗

DOE OSTI2024

U(1) fields from qubits: An approach via D-theory algebra

A new quantum link microstructure was proposed for the lattice quantum chromodynamics (QCD) Hamiltonian, replacing the Wilson gauge links with a bilinear of fermionic qubits, later generalized to D-theory. This formalism provides a general framework for building lattice field theory algorithms for quantum computing. We focus mostly on the simplest case of a quantum rotor for a single compact U(1) field. We also make some progress for non-Abelian setups, making it clear that the ideas developed in the U(1) case extend to other groups. These in turn are building blocks for 1 + 0 -dimensional ( 1 + 0 -D) matrix models, 1 + 1 -D sigma models and non-Abelian gauge theories in 2 + 1 and 3 + 1 dimensions. By introducing multiple flavors for the U(1) field, where the flavor symmetry is gauged, we can efficiently approach the infinite-dimensional Hilbert space of the quantum O(2) rotor with increasing flavors. The emphasis of the method is on preserving the symplectic algebra exchanging fermionic qubits by sigma matrices (or hard bosons) and developing a formal strategy capable of generalization to a SU ( 3 ) field for lattice QCD and other non-Abelian 1 + 1 -D sigma models or 3 + 1 -D gauge theories. For U(1), we discuss briefly the qubit algorithms for the study of the discrete 1 + 1 -D sine-Gordon equation. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

DOE OSTI2020

LDA + $U_{sc}$ calculations of phase relations in FeO

Using the LDA + $U_{sc}$ method, we present calculation phase relations of iron monoxides involving five polytypes in multiple spin-state configurations. In this work, the Hubbard parameter $\textit{U}$ is determined self-consistently simultaneously with the occupation matrix and structures at arbitrary pressures. The Hubbard parameter strongly depends on pressure, structure, and spin state. Comparison with experimental structural data indicates the LDA + $U_{sc}$ can predict structure, compression curves, phase relations, and transition pressures very well for the insulating $\textit{B}$1 and $\textit{iB}$8 states. However, it requires additional calculations using the Mermin functional that includes the electronic entropic contribution to the free energy to obtain an $\textit{nB}$8 metallic state and a consistent $\textit{iB}$8 to $\textit{nB}$8 insulator to metal transition pressure.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Quantum Simulators and Applications on Quantum Framework

Simulating quantum circuits is essential for validating quantum algorithms. However, no single simulator consistently performs best - efficiency depends on circuit structure, entanglement, and depth. In this work, we integrate Qiskit-Aer (state-vector and matrix product state) and QTensor, a tree-tensor-network based simulator, into the Quantum Framework (QFw), a modular platform that supports multiple quantum backends via a unified interface. We also enable distributed quantum approximate optimization algorithm (DQAOA) application compatibility with QFw, allowing sub-problems to be solved in parallel at scale. We then benchmark DQAOA and TFIM (transverse field Ising model) circuits across supported simulators, showing how performance varies significantly with problem type. All simulations are deployed on the Frontier supercomputer using QFw's MPI-based orchestration for distributed, multinode execution. These results underscore the need for simulatoragnostic infrastructure to enable systematic evaluation and highperformance scaling of quantum workloads. QFw provides a practical and extensible path toward reproducible quantum algorithm development across diverse application domains.

Chundury, Srikar [ORNL] (ORCID:0009000183359259)↗

DOE OSTI2024

Optimizing the sextupole configuration for simultaneous correction of third order resonances at the recycler ring

For the Recycler Ring at Fermilab, space charge tune shifts of almost 0.1 will have to be dealt with under the Proton Improvement Plan (PIP-II) framework. This will lead to the excitation of third order resonances. The minimization of Resonance Driving Terms (RDTs) allows to mitigate the harmful effect of these betatron resonances. Past work has shown that previously-installed sextupoles can compensate the RDTs of individual third order resonance lines, thus reducing particle losses in these operational regimes. Nevertheless, trying to compensate multiple resonances of the same order simultaneously with these existing sextupoles is limited due to current constraints in the magnets. The following study showcases the procedure to install additional sextupoles in order to aid the compensation of multiple resonances. This includes the optimization of the new sextupoles' locations in order to cancel out multiple RDTs while minimizing the currents needed. This is followed by a verification of their effectiveness by means of the RDT response matrix.

43 PARTICLE ACCELERATORS↗

DOE OSTI2022

Blast Effects on Buildings (Final Report)

Lawrence Livermore National Laboratory (LLNL) has conservatively reduced the explosive safety standards suggested by the Small Quantities in Research Laboratories (SQRL) program testing by 60%. For this reason, further research into the detailed effects of small amounts of explosives in typically constructed rooms is needed in order to improve their factor of safety. The team was tasked with designing an experiment to investigate the effect of different variables on drywall under explosive blasts. In order to meet this objective, the team conducted a comprehensive literature review to gain an understanding of industry-standard construction practices and review previous tests conducted by the Army corp of engineers. Since the team used PBXN-5 rather than C4, as in SQRL, an initial shot was conducted to compare damage levels. From those results and the physical constraints of the testing chamber, the team redesigned multiple single-panel drywall frames to capture the entirety of the incurred damage. Proposed designs were narrowed down using a decision matrix. From the study of previous tests and literature review, variables were chosen that the team hypothesized to have an impact on drywall strength. The variables that were tested were chosen from the results of that work and specific variables the sponsor was interested in, and they were paint, humidity/moisture content, and explosive positioning relative to the studs. Detailed plans were made for each variable according to what conditions the team wanted to investigate. For humidity, this involved testing low, ambient, and high conditions by treating the panels in a chamber. Preliminary shots were performed to test the structural integrity of the frame and streamline the test diagnostics which involved a high-speed camera placed behind the drywall, outside of the chamber, and a pressure probe placed behind the drywall, inside the chamber. Once the instrumentation, diagnostics, and frame design were finalized, a quantitative damage criteria matrix was created to categorize the results of the main shot series. In conjunction with evidence from the high-speed video, the achieved damage levels indicate that high moisture content drywall is better able to withstand explosive blasts. Larger stud damage and lower drywall damage occurred when the explosive was located directly in front of a stud. Paint had no noticeable impact on strength. Ultimately, the team conducted a total of 17 tests, leaving the door open for future in-depth research into the impact of humidity.

36 MATERIALS SCIENCE↗

NASA NTRS1972

The k-space formulation of the n-dimensional scattering problem

The n-dimensional scattering problem is solved by means of a k-space formulation of the field equations, thereby replacing the conventional integral equation formulation by a set of two algebraic equations in two unknowns in two spaces (the constitutive equation being an algebraic equation in x-space). These equations are solved by an iterative method with the aid of the fast Fourier transform (FFT) algorithm connecting the two spaces, requiring very simple initial approximations. Since algebraic and FFT equations are used, the number of arithmetic multiple-add operations and storage allocations required for a numerical solution are reduced from the order of N sq (for solving the matrix equations resulting from the conventional integral equations) to the order of N(log base 2 of N) and N, respectively (where N is the number of data points required for the specification of the problem). The advantage gained in speed and storage is thus of the order of N/log base 2 of N and N, respectively. This method is thus considerably more efficient than the conventional matrix method, and permits exact numerical solutions for much larger problems. Arguments are presented toward the view that the field equations are more fundamental in k-space. The details and some numerical results of the application of this method to the three-dimensional electromagnetic scattering problems are presented as an example.

Bojarski, N. N.↗

NASA NTRS1983

Evaluation of the effect of vibration nonlinearity on convergence behavior of adaptive higher harmonic controllers

Effect of nonlinearity on convergence of the local linear and global linear adaptive controllers is evaluated. A nonlinear helicopter vibration model is selected for the evaluation which has sufficient nonlinearity, including multiple minimum, to assess the vibration reduction capability of the adaptive controllers. The adaptive control algorithms are based upon a linear transfer matrix assumption and the presence of nonlinearity has a significant effect on algorithm behavior. Simulation results are presented which demonstrate the importance of the caution property in the global linear controller. Caution is represented by a time varying rate weighting term in the local linear controller and this improves the algorithm convergence. Nonlinearity in some cases causes Kalman filter divergence. Two forms of the Kalman filter covariance equation are investigated.

Molusis, J. A.↗

NASA NTRS1987

Effective propagation constants in dense random media under effective medium approximation

The effective medium approximation is applied to study the effective propagation constants in a dense random medium. The dyadic Green's function is introduced to establish the effective medium approximation formalism for electromagnetic waves. The multiple scattering equations and the Lippmann-Schwinger equations for the transition operator are obtained in the configuration average form. The dispersion equations of multiple scattering is derived by using a standard method in quantum mechanics. To obtain an expression for the effective propagation constants the matrix elements of the configuration average dyadic transiton operator are calculated in momentum representation. Numerical illustations are carried out to demonstrate the difference in the effective propagation constants between the use of this approximation and the well-known quasi-crystalline approximation. A comparison is made with measured loss tangent in dry snow.

Zhu, P. Y.↗

NASA NTRS1988

An empirical model for polarized and cross-polarized scattering from a vegetation layer

An empirical model for scattering from a vegetation layer above an irregular ground surface is developed in terms of the first-order solution for like-polarized scattering and the second-order solution for cross-polarized scattering. The effects of multiple scattering within the layer and at the surface-volume boundary are compensated by using a correction factor based on the matrix doubling method. The major feature of this model is that all parameters in the model are physical parameters of the vegetation medium. There are no regression parameters. Comparisons of this empirical model with theoretical matrix-doubling method and radar measurements indicate good agreements in polarization, angular trends, and k sub a up to 4, where k is the wave number and a is the disk radius. The computational time is shortened by a factor of 8, relative to the theoretical model calculation.

Liu, H. L.↗

NASA NTRS2004

Microsystem Cooler Development

A patented microsystem Stirling cooler is under development with potential application to electronics, sensors, optical and radio frequency (RF) systems, microarrays, and other microsystems. The microsystem Stirling cooler is most suited to volume-limited applications that require cooling below the ambient or sink temperature. Primary components of the planar device include: two diaphragm actuators that replace the pistons found in traditional-scale Stirling machines; and a micro-regenerator that stores and releases thermal energy to the working gas during the Stirling cycle. The use of diaphragms eliminates frictional losses and bypass leakage concerns associated with pistons, while permitting reversal of the hot and cold sides of the device during operation to allow precise temperature control. Three candidate microregenerators were custom fabricated for initial evaluation: two constructed of porous ceramic, and one made of multiple layers of nickel and photoresist in an offset grating pattern. An additional regenerator was prepared with a random stainless steel fiber matrix commonly used in existing Stirling machines for comparison to the custom fabricated regenerators. The candidate regenerators were tested in a piezoelectric-actuated test apparatus designed to simulate the Stirling refrigeration cycle. In parallel with the regenerator testing, electrostatically-driven comb-drive diaphragm actuators for the prototype device have been designed for deep reactive ion etching (DRIE) fabrication.

Moran, Matthew E.↗

NASA NTRS2018

Electrical Properties of Unidirectional Carbon Fiber Reinforced Epoxy-Graphene Composites

This paper aims at establishing connections between electrical resistivity and stresses in carbon fiber reinforced epoxy-graphene composites. To avoid the interconnection of multiple parameters, we examined composite specimens consisting of a single strand or tow (12,000 carbon fiber filaments) embedded in an epoxy matrix containing up to 1.0 wt% of graphene (TG679) particles. The specimens show resistance increase with the tensile stress. The characteristic of the dependence varies with the content of graphene particles, tending to be linear as the percentage of graphene content increases. We did not observe any significant changes in the ultimate strength of the composite with an increase of the graphene particles concentration.

Daniel Wentzel↗

NASA NTRS2022

Deep Interacting Multiple Model Filtering

In this paper, a deep learning-based multiple model estimation framework is presented for the state estimation of hybrid dynamical systems from high dimensional observations such as camera images. A low dimensional vector which represents the measurement of the latent dynamical system and its corresponding variance are learned using a deep encoder neural network. An Interacting Multiple Model (IMM) filter is used to generate the latent state estimates and covariances using multiple dynamical models, which can be learned using backpropagation through time. The state estimates of the dynamical system and the corresponding covariance matrix are generated from the latent state estimates and covariance using a deep decoder neural network. The whole network is trained in an end-to-end manner using a loss function which minimizes the negative log-likelihood of the neural network parameters. Simulation results are presented using a 2D bouncing ball example and estimation error statistics are computed which demonstrates the accuracy and consistency of the estimation.

Ghananeel Rotithor↗