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 649 records · Page 36

In-Situ Scanning Electron Microscope Experiments for Microscale Mechanical Testing and Validated Modeling of Fiber Reinforced Thermoplastics

A novel, in-situ, scanning electron microscope (SEM) mechanical testing capability for materials at the microscale which provides experimental validation to a machine learning (ML) toolset for full-field validation of physics-based micromechanics models is being developed by researchers at NASA Glenn Research Center. These are enabling technologies for the integration of multiscale digital twins for materials into system level models which will result in the improved performance, material discovery, reduced production cost and time, rapid characterization, and prognostic structural health monitoring (SHM) for materials and structures for extreme environments in support of NASA space exploration missions. In order to bridge the material structure-to-system gap for digital twins, physics-based models must be experimentally validated at multiple length scales. Seminal microscale experiments, conducted at the Air Force Research Laboratory (AFRL), were limited to transverse compression of single-layer, unidirectional thermoset polymer matrix composite (PMC) micropillar specimens [1]. The early phases of the current project followed those initial results and setup to reproduce the compression testing of PMC material on the custom-built piezoelectric actuated micromechanical testing rig built by MicroTesting Solutions LLC. In this work, samples of thermoplastic PMC material were first machined into 3 mm cubes, and then further machining and final milling was done using a Focused Ion Beam (FIB). The initial experiment was done on a pillar roughly 20 µm x 20 µm x 40 µm tall. Additional pillars were milled with final sizes ranging from 20 µm x 20 µm x 40 µm tall to 40 µm x 40 µm x 65 µm tall. A speckle pattern for in-situ full-field measurements using Digital Image Correlation (DIC) was applied with platinum, which was coated on the surface, and then the FIB was used to mill away some of the coating to produce an irregular pattern of Pt on the pillar surface. The samples were loaded into the custom testing rig and placed into the SEM and loaded under compression until failure. Images were collected in the SEM during testing. Post-processing of the images was conducted using DIC to obtain full-field displacement and strain measurements elucidating the role of the matrix as well as fiber-fiber interaction at the microscale within the composite subjected to compression loading well into the non-linear regime of the material. Moreover, the evolution of fiber-matrix debonding and matrix cracking is observed in-situ at the microscale. This data, along with images segmented with a newly developed ML toolset [2], was used to create and validate physics-based micromechanics models. An image of the failed micropillar is shown in Figure 1. The techniques developed in the initial compression experiment was tailored to the validation needs of the models and expanded to include different sized samples as well as possibly tension and fatigue.

Laura Wilson↗

NXREFINE

NXRefine implements a complete workflow for both data acquisition and reduction of single crystal diffuse x-ray scattering collected as a series of area detector frames during continuous sample rotation. The data are stored in HDF5 files, which conform to the NeXus data format standard. These files contain a comprehensive set of metadata, including the incident wavelength, powder calibration scans used to define detector distances, beam centers, and yaw, pitch, and roll corrections. Automated peak searches of the three-dimensional data arrays define a set of Bragg peaks that are used to define an orientation matrix, defined according to Busing and Levy, which is then used to transform the data from angular coordinates to reciprocal space coordinates using the CCTW program. Multiple rotations at different detector translations and offset rotation axes are combined to ensure that there are no gaps in the reconstructed data, and allow the construction of masks to eliminate detector artifacts due to Compton scattering within the sensor layers. Finally, the data can be transformed into 3D-ΔPDF maps using the punch-and-fill method after symmetrization.

OSBORN, RAYMOND↗

Temporal Sampling of Enzymes from Live Cells by Localized Electroporation and Quantification of Activity by SAMDI Mass Spectrometry

Abstract Measuring changes in enzymatic activity over time from small numbers of cells remains a significant technical challenge. In this work, a method for sampling the cytoplasm of cells is introduced to extract enzymes and measure their activity at multiple time points. A microfluidic device, termed the live cell analysis device (LCAD), is designed, where cells are cultured in microwell arrays fabricated on polymer membranes containing nanochannels. Localized electroporation of the cells opens transient pores in the cell membrane at the interface with the nanochannels, enabling extraction of enzymes into nanoliter‐volume chambers. In the extraction chambers, the enzymes modify immobilized substrates, and their activity is quantified by self‐assembled monolayers for matrix‐assisted laser desorption/ionization (SAMDI) mass spectrometry. By employing the LCAD‐SAMDI platform, protein delivery into cells is demonstrated. Next, it is shown that enzymes can be extracted, and their activity measured without a loss in viability. Lastly, cells are sampled at multiple time points to study changes in phosphatase activity in response to oxidation by hydrogen peroxide. With this unique sampling device and label‐free assay format, the LCAD with SAMDI enables a powerful new method for monitoring the dynamics of cellular activity from small populations of cells.

Mukherjee, Prithvijit↗

Depletion capabilities in the OpenMC Monte Carlo particle transport code

A depletion solver has been implemented in OpenMC and is described herein. The depletion solver is implemented in Python and interfaces with OpenMC’s transport solver through a C++ application programming interface, which enables an in-memory transport-depletion coupling. Multiple integration methods for advancing in time have been implemented and exhibit tradeoffs in cost, accuracy, and memory use. For all time integration methods, evaluation of the matrix exponential is performed by using the incomplete partial fraction form of the Chebyshev rational approximation method. Simulations of a pressurized water reactor (PWR) pincell and a sodium-cooled fast reactor (SFR) assembly were carried out with OpenMC and Serpent. For both problems, the use of a high-fidelity depletion chain results in predictions of k eff that agree within 20–30 pcm between OpenMC and Serpent. Predicted actinide concentrations were found to agree to a fraction of a percent, and most fission product concentrations were found to agree within 1%. Here, the few cases where larger differences were observed can be attributed either to differences in how the energy dependence of fission product yields is handled or deficiencies in the nuclear data used. OpenMC simulations of the PWR and SFR problems using a simplified 228-nuclide depletion chain demonstrate that it achieves accuracy close to that of the full, high-fidelity depletion chain with respect to the studied responses.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The convergence of quadrupole rotational invariants from the nuclear shell model

Nuclei exhibit both single-particle and collective degrees of freedom, with the latter often subdivided into vibrational and rotational motions. Experimentally identifying the relative roles of these collective modes is extremely challenging, particularly in the face of possible shape coexistence. Model-independent, invariant quantities describing the deformation of a nucleus in the intrinsic frame have long been known but their determination potentially requires a large quantity of experimental data to achieve convergence. Through comparison with the nuclear shell model, the question of convergence is addressed. Shell-model calculations performed in the sd- and pf-shell model spaces are used to determine electric-quadrupole matrix elements for a multitude of low-lying states using the first 40 states of the relevant spins. Relative contributions to the rotationally invariant quantities from multiple states can therefore be determined. It is found that, on average, the inclusion of four intermediate states results in the leading-order invariant, $\langle\hat{Q^2}\rangle$, converging to within 10% of its true value and the triaxiality term, cos (3δ), converging to its true value, though some variance remains. Higher-order quantities relating to the softness of the nuclear shape are found to converge more slowly. The convergence of quadrupole rotationally invariant sum rules was quantified in the sd- and pf-shell model spaces and indicates the challenge inherent in a full determination of nuclear shape. The present study is limited to relatively small valence spaces. Finally, larger spaces, such as the rare-earth region, potentially offer faster convergence.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Measurement of the Higgs boson production rate in association with top quarks in final states with electrons, muons, and hadronically decaying tau leptons at $\sqrt{s} =$ 13 TeV

The rate for Higgs (${\mathrm{H}} $) bosons production in association with either one (${\mathrm{t}} {\mathrm{H}} $) or two (${\mathrm{t}} {{\overline{{{\mathrm{t}}}}}} {\mathrm{H}} $) top quarks is measured in final states containing multiple electrons, muons, or tau leptons decaying to hadrons and a neutrino, using proton–proton collisions recorded at a center-of-mass energy of $13\,\text {Te}\text {V} $ by the CMS experiment. The analyzed data correspond to an integrated luminosity of 137$\,\text {fb}^{-1}$. The analysis is aimed at events that contain ${\mathrm{H}} \rightarrow {\mathrm{W}} {\mathrm{W}} $, ${\mathrm{H}} \rightarrow {\tau } {\tau } $, or ${\mathrm{H}} \rightarrow {\mathrm{Z}} {\mathrm{Z}} $ decays and each of the top quark(s) decays either to lepton+jets or all-jet channels. Sensitivity to signal is maximized by including ten signatures in the analysis, depending on the lepton multiplicity. The separation among ${\mathrm{t}} {\mathrm{H}} $, ${\mathrm{t}} {{\overline{{{\mathrm{t}}}}}} {\mathrm{H}} $, and the backgrounds is enhanced through machine-learning techniques and matrix-element methods. The measured production rates for the ${\mathrm{t}} {{\overline{{{\mathrm{t}}}}}} {\mathrm{H}} $ and ${\mathrm{t}} {\mathrm{H}} $ signals correspond to $0.92 \pm 0.19\,\text {(stat)} ^{+0.17}_{-0.13}\,\text {(syst)} $ and $5.7 \pm 2.7\,\text {(stat)} \pm 3.0\,\text {(syst)} $ of their respective standard model (SM) expectations. The corresponding observed (expected) significance amounts to 4.7 (5.2) standard deviations for ${\mathrm{t}} {{\overline{{{\mathrm{t}}}}}} {\mathrm{H}} $, and to 1.4 (0.3) for ${\mathrm{t}} {\mathrm{H}} $ production. Assuming that the Higgs boson coupling to the tau lepton is equal in strength to its expectation in the SM, the coupling $y_{{\mathrm{t}}}$ of the Higgs boson to the top quark divided by its SM expectation, $\kappa _{{\mathrm{t}}}=y_{{\mathrm{t}}}/y_{{\mathrm{t}}}^{\mathrm {SM}}$, is constrained to be within $-0.9< \kappa _{{\mathrm{t}}}< -0.7$ or $0.7< \kappa _{{\mathrm{t}}}< 1.1$, at 95% confidence level. This result is the most sensitive measurement of the ${\mathrm{t}} {{\overline{{{\mathrm{t}}}}}} {\mathrm{H}} $ production rate to date.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Compensator improvement for multivariable control systems

A theory and the associated numerical technique are developed for an iterative design improvement of the compensation for linear, time-invariant control systems with multiple inputs and multiple outputs. A strict constraint algorithm is used in obtaining a solution of the specified constraints of the control design. The result of the research effort is the multiple input, multiple output Compensator Improvement Program (CIP). The objective of the Compensator Improvement Program is to modify in an iterative manner the free parameters of the dynamic compensation matrix so that the system satisfies frequency domain specifications. In this exposition, the underlying principles of the multivariable CIP algorithm are presented and the practical utility of the program is illustrated with space vehicle related examples.

Mitchell, J. R.↗

H2, fixed architecture, control design for large scale systems

The H2, fixed architecture, control problem is a classic linear quadratic Gaussian (LQG) problem whose solution is constrained to be a linear time invariant compensator with a decentralized processing structure. The compensator can be made of p independent subcontrollers, each of which has a fixed order and connects selected sensors to selected actuators. The H2, fixed architecture, control problem allows the design of simplified feedback systems needed to control large scale systems. Its solution becomes more complicated, however, as more constraints are introduced. This work derives the necessary conditions for optimality for the problem and studies their properties. It is found that the filter and control problems couple when the architecture constraints are introduced, and that the different subcontrollers must be coordinated in order to achieve global system performance. The problem requires the simultaneous solution of highly coupled matrix equations. The use of homotopy is investigated as a numerical tool, and its convergence properties studied. It is found that the general constrained problem may have multiple stabilizing solutions, and that these solutions may be local minima or saddle points for the quadratic cost. The nature of the solution is not invariant when the parameters of the system are changed. Bifurcations occur, and a solution may continuously transform into a nonstabilizing compensator. Using a modified homotopy procedure, fixed architecture compensators are derived for models of large flexible structures to help understand the properties of the constrained solutions and compare them to the corresponding unconstrained ones.

Mercadal, Mathieu↗

Error Analyses of the North Alabama Lightning Mapping Array (LMA)

Two approaches are used to characterize how accurately the North Alabama Lightning Mapping Array (LMA) is able to locate lightning VHF sources in space and in time. The first method uses a Monte Carlo computer simulation to estimate source retrieval errors. The simulation applies a VHF source retrieval algorithm that was recently developed at the NASA-MSFC and that is similar, but not identical to, the standard New Mexico Tech retrieval algorithm. The second method uses a purely theoretical technique (i.e., chi-squared Curvature Matrix theory) to estimate retrieval errors. Both methods assume that the LMA system has an overall rms timing error of 50ns, but all other possible errors (e.g., multiple sources per retrieval attempt) are neglected. The detailed spatial distributions of retrieval errors are provided. Given that the two methods are completely independent of one another, it is shown that they provide remarkably similar results, except that the chi-squared theory produces larger altitude error estimates than the (more realistic) Monte Carlo simulation.

Koshak, W. J.↗

Error Analyses of the North Alabama Lightning Mapping Array (LMA)

Two approaches are used to characterize how accurately the North Alabama Lightning Mapping Array (LMA) is able to locate lightning VHF sources in space and in time. The first method uses a Monte Carlo computer simulation to estimate source retrieval errors. The simulation applies a VHF source retrieval algorithm that was recently developed at the NASA-MSFC and that is similar, but not identical to, the standard New Mexico Tech retrieval algorithm. The second method uses a purely theoretical technique (i.e., chi-squared Curvature Matrix theory) to estimate retrieval errors. Both methods assume that the LMA system has an overall rms timing error of 50ns, but all other possible errors (e.g., multiple sources per retrieval attempt) are neglected. The detailed spatial distributions of retrieval errors are provided. Given that the two methods are completely independent of one another, it is shown that they provide remarkably similar results, except that the chi-squared theory produces larger altitude error estimates than the (more realistic) Monte Carlo simulation.

Koshak, W. J.↗

Systems Engineering and Integration for Advanced Life Support System and HST

Systems engineering (SE) discipline has revolutionized the way engineers and managers think about solving issues related to design of complex systems: With continued development of state-of-the-art technologies, systems are becoming more complex and therefore, a systematic approach is essential to control and manage their integrated design and development. This complexity is driven from integration issues. In this case, subsystems must interact with one another in order to achieve integration objectives, and also achieve the overall system's required performance. Systems engineering process addresses these issues at multiple levels. It is a technology and management process dedicated to controlling all aspects of system life cycle to assure integration at all levels. The Advanced Integration Matrix (AIM) project serves as the systems engineering and integration function for the Human Support Technology (HST) program. AIM provides means for integrated test facilities and personnel for performance trade studies, analyses, integrated models, test results, and validated requirements of the integration of HST. The goal of AIM is to address systems-level integration issues for exploration missions. It will use an incremental systems integration approach to yield technologies, baselines for further development, and possible breakthrough concepts in the areas of technological and organizational interfaces, total information flow, system wide controls, technical synergism, mission operations protocols and procedures, and human-machine interfaces.

Kamarani, Ali K.↗

Materials for Heated Head Automated Thermoplastic Tape Placement

NASA Langley Research Center (LaRC) is currently pursuing multiple paths to develop out of autoclave (OOA) polymeric composite materials and processes. Polymeric composite materials development includes the synthesis of new and/or modified thermosetting and thermoplastic matrix resins designed for specific OOA processes. OOA processes currently under investigation include vacuum bag only (VBO) prepreg/composite fabrication, resin transfer molding (RTM), vacuum assisted resin transfer molding (VARTM) and heated head automated thermoplastic tape placement (HHATP). This paper will discuss the NASA Langley HHATP facility and capabilities and recent work on characterizing thermoplastic tape quality and requirements for quality part production. Samples of three distinct versions of APC-2 (AS4/PEEK) thermoplastic dry tape were obtained from two materials vendors, TENCATE, Inc. and CYTEC Engineered Materials** (standard grade and an experimental batch). Random specimens were taken from each of these samples and subjected to photo-microscopy and surface profilometry. The CYTEC standard grade of APC-2 tape had the most voids and splits and the highest surface roughness and/or waviness. Since the APC-2 tape is composed of a thermoplastic matrix, it offers the flexibility of reprocessing to improve quality, and thereby improve final quality of HHATP laminates. Discussions will also include potential research areas and future work that is required to advance the state of the art in the HHATP process for composite fabrication.

Jensen, Brian J.↗

HyKKT: a hybrid direct-iterative method for solving KKT linear systems

Here, we propose a solution strategy for the large indefinite linear systems arising in interior methods for nonlinear optimization. The method is suitable for implementation on hardware accelerators such as graphical processing units (GPUs). The current gold standard for sparse indefinite systems is the LBLT factorization where L is a lower triangular matrix and B is 1×1 or 2×2 block diagonal. However, this requires pivoting, which substantially increases communication cost and degrades performance on GPUs. Our approach solves a large indefinite system by solving multiple smaller positive definite systems, using an iterative solver on the Schur complement and an inner direct solve (via Cholesky factorization) within each iteration. Cholesky is stable without pivoting, thereby reducing communication and allowing reuse of the symbolic factorization. We demonstrate the practicality of our approach on large optimal power flow problems and show that it can efficiently utilize GPUs and outperform LBL T factorization of the full system.

97 MATHEMATICS AND COMPUTING↗

Efficient Parallelization of Irregular Applications on GPU Architectures

With the enlarging computation capacity of general Graphics Processing Units (GPUs), leveraging GPUs to accelerate parallel applications has become a critical topic in academia and industry. However, a wide range of irregular applications with the computation-/memory-intensive nature cannot easily achieve high GPU utilization. The challenges mainly involve the following aspects: first, data dependence leads to coarse-grained kernel and inefficient parallelism; second, heavy GPU memory usage may cause frequent memory evictions and extra overhead of I/O; third, specific computation patterns produce memory redundancies; last, workload balance and data reusability conjunctly benefit the overall performance, but there may exist a dynamic trade-off between them. Targeting these challenges, this dissertation proposes multiple optimizations to accelerate two real-world applications: many-body correlation functions to simulate nuclear physics in a large-scale scientific system; the other is the eALS-based matrix factorization recommendation system. To accelerate the calculations of many-body correlation functions, this dissertation presents three frameworks in GPU memory management and multi-GPU scheduling. Firstly, an optimized systematic GPU memory management framework, MemHC, utilizes a series of new memory reduction designs in GPU memory allocation, CPU/GPU communications, and GPU memory oversubscription. Secondly, an enhanced multi-GPU scheduling framework, MICCO, particularly by taking both data dimension (e.g., data reuse and data eviction) and computation dimension into account. MICCO designs a heuristic scheduling algorithm and a machine learning-based regression model to generate the optimal settings of a proposed new concept to manage the trade-off. Thirdly, a locality-aware multi-GPU scheduling framework. This scheduler leverages pipeline batch generation with a looking-ahead strategy by building local dependency graphs for memory transfer reduction and better data reuse, achieving up to 79.92% memory cost reduction and 1.67x speedup. To parallelize the eALS-based recommendation system, this dissertation proposes an efficient CPU/GPU heterogeneous recommendation system, HEALS. HEALS employs newly designed architecture-adaptive data formats to achieve load balance and good data locality on CPU and GPU. To mitigate the data dependence, HEALS presents a CPU/GPU collaboration model for both task parallelism and data parallelism with multiple kernel computation optimizations. In summary, this dissertation efficiently accelerates two typical irregular applications on GPUs by building four frameworks, including CPU/GPU collaboration, GPU memory management, and multi-GPU scheduling.

Wang, Qihan↗

In Situ Scattering Studies of Crystallization Kinetics in a Phase-Separated Zr–Cu–Fe–Al Bulk Metallic Glass

In this work, thermal stability and the crystallization kinetics of a phase-separated Zr-Cu-Fe-Al bulk metallic glass were investigated using in situ high-energy synchrotron X-ray and neutron diffraction, as well as small-angle synchrotron X-ray scattering. It was revealed that this glass with excellent glass-forming ability possesses a two-step crystallization behavior. The crystalline products and their evolution sequence are more complicated than a homogeneous Zr-Cu-Al glass with average glass-forming ability. The experimental results indicate that a finely distributed nanometer-sized cubic Zr 2 Cu phase forms first and then transforms to a tetragonal Zr 2 Cu phase, while the matrix transforms to an orthorhombic Zr3Fe phase. The strength of the Zr-Cu-Fe-Al composite containing cubic Zr 2 Cu phase and glass matrix increases, and the plasticity also improves compared to the as-cast Zr-Cu-Fe-Al bulk metallic glass. Our results suggest that the formation of multiple and complex crystalline products would be the characteristics of the Zr-Cu-Fe-Al glass with better glass-forming ability. Our study may shed light on the synthesis of bulk-sized glass-nanocrystals composites of high strength and good plasticity.

36 MATERIALS SCIENCE↗

Testing and signal identification for two-sample high-dimensional covariances via multi-level thresholding

The paper considers testing and signal identification for covariance matrices from two populations of marginally sub-Gaussian distributed. A multi-level thresholding procedure is proposed for testing the equality of two high-dimensional covariance matrices, which is designed to detect sparse and faint differences between the covariances. A novel -statistic composition is developed to establish the asymptotic distribution of the thresholding statistics in conjunction with the matrix blocking and the coupling techniques. It is shown that the proposed test is more powerful than the existing tests in detecting sparse and weak signals in covariances. Multiple testing procedures are constructed to discover different covariances and the sub-groups of variables with different covariance structures between the two populations. The proposed procedures are based on the multi-level thresholding test, which are able to control the false discovery proportion () with high power. Simulation experiments and a case study on the returns of the S&P 500 stocks before and after the COVID-19 pandemic are conducted to demonstrate and compare the utilities of the proposed methods.

97 MATHEMATICS AND COMPUTING↗

Monte Carlo Event Generation with Continuous Normalizing Flows

We apply continuous normalizing flows trained with the flow matching method to the problem of phase-space sampling in Monte Carlo event generation for high-energy collider physics. Focusing on lepton-pair and top-quark pair production with multiple jets, the two computationally most expensive processes at the Large Hadron Collider, we train helicity-conditioned continuous normalizing flows to remap the random numbers used in matrix element evaluation. Compared to standard methods, we achieve unweighting efficiency improvements by factors of up to 184 and 25 for the two processes at their respective highest jet number, at the cost of an increased evaluation time. When combining the advantages of continuous normalizing flows with the fast evaluation times of coupling-layer-based flows, using the RegFlow approach, we find parton-level unweighted event generation walltime gains of about a factor of 10 at the highest jet numbers. These substantial gains highlight the promise of samplers based on machine learning for next-generation collider experiments.

Bothmann, Enrico [CERN; Gottingen U.] (ORCID:00000↗

Towards a Quantum Algorithm for the Incompressible Nonlinear Navier-Stokes Equations

In this work, we present novel concepts for quantum algorithms to solve transient, nonlinear partial differential equations (PDEs). The challenge lies in how to effectively represent, encode, process, and evolve the nonlinear system of PDEs on quantum computers. We will discuss the new techniques using the incompressible Navier-Stokes equations as an example, because it represents the fundamental nonlinear feature and yet removes certain complexity in physics, allowing us to focus on the design of quantum algorithms. Previous attempts solving nonlinear PDEs in quantum computation have often involved storing multiple copies of solutions or employing linearizations. Neither is practical due to exponential scaling with evolution time or insufficient solution accuracy. We propose a new framework based on matrix product states (MPSs) and matrix product operators (MPOs), in addition to the Krylov subspace methods. For example, the solution variables of the Navier-Stokes equations are represented by MPSs, and the linear and nonlinear terms are processed by MPOs. The time evolution of the operators is attained by a fast-forwarding algorithm using Krylov subspace methods. Furthermore, we discuss various techniques for efficient encoding of MPSs, measurement reduction for MPOs, and use of tensor operations to treat multi-variate, multi-physics characteristics of Navier-Stokes.

Gopalakrishnan Meena, Murali [ORNL] (ORCID:0000000↗