SEARCH · Engineering Papers

Results for “Matrix”

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 577 records · Page 32

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Modeling Aerial Gamma-Ray Backgrounds Using Non-negative Matrix Factorization

Airborne gamma-ray surveys are useful for many applications, ranging from geology and mining to public health and nuclear security. In all these contexts, the ability to decompose a measured spectrum into a linear combination of background source terms can provide useful insights into the data and lead to improvements in the techniques that use spectral energy windows. Multiple methods for the linear decomposition of spectra exist but are subject to various drawbacks, such as allowing negative photon fluxes or requiring detailed Monte Carlo modeling. In this work, we propose using non-negative matrix factorization (NMF) as a data-driven approach to spectral decomposition. Using aerial surveys that include flights over water, we demonstrate that the mathematical approach of NMF finds physically relevant structure in the aerial gamma-ray background, namely, that measured spectra can be expressed as the sum of nearby terrestrial emission, distant terrestrial emission, and radon and cosmic emission. Furthermore, these NMF background components are compared with the background components obtained by noise-adjusted singular value decomposition (NASVD), which contain negative photon fluxes and, thus, do not represent the emission spectra in as straightforward a way. Finally, we comment on the potential areas of research that are enabled by NMF decompositions, such as new approaches to spectral anomaly detection and data fusion.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Identifying Entangled Physics Relationships Through Sparse Matrix Decomposition to Inform Plasma Fusion Design

We report a sustainable burn platform through inertial confinement fusion (ICF) has been an ongoing challenge for over 50 years. Mitigating engineering limitations and improving the current design involves an understanding of the complex coupling of physical processes. While sophisticated simulation codes are used to model ICF implosions, these tools contain necessary numerical approximation but miss physical processes that limit predictive capability. Identification of relationships between controllable design inputs to ICF experiments and measurable outcomes (e.g., neutron yield, neutron velocity, areal density) from performed experiments can help guide the future design of experiments and development of simulation codes, to potentially improve the accuracy of the computational models used to simulate ICF experiments. We use sparse matrix decomposition methods to identify clusters of a few related design variables. Sparse principal component analysis (SPCA) identifies groupings that are related to the physical origin of the variables (laser, hohlraum, and capsule). A variable importance analysis finds that in addition to variables highly correlated with neutron yield, such as picket power and laser energy, variables that represent a dramatic change of the ICF design, such as number of pulse steps, are also very important. The obtained sparse components are then used to train a random forest (RF) regression surrogate for predicting total yield. The RF performance on the training and testing data compares with the performance of the RF trained using all the design variables considered. This work is intended to inform design changes in future ICF experiments by augmenting the expert intuition and simulation results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Algorithms for Non-Negative Matrix Factorization on Noisy Data With Negative Values

Non-negative matrix factorization (NMF) is a dimensionality reduction technique that has shown promise for analyzing noisy data, especially astronomical data. For these datasets, the observed data may contain negative values due to noise even when the true underlying physical signal is strictly positive. Prior NMF work has not treated negative data in a statistically consistent manner, which becomes problematic for low signal-to-noise data with many negative values. In this paper we present two algorithms, Shift-NMF and Nearly-NMF, that can handle both the noisiness of the input data and also any introduced negativity. Both of these algorithms use the negative data space without clipping or masking and recover non-negative signals without any introduced positive offset that occurs when clipping or masking negative data. We demonstrate this numerically on both simple and more realistic examples, and prove that both algorithms have monotonically decreasing update rules.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

In Situ Oxygen Isotope Determination in Serpentine Minerals by SIMS: Addressing Matrix Effects and Providing New Insights on Serpentinisation at Hole BA1B (Samail ophiolite, Oman)

The ability to constrain the petrogenesis of multiple serpentine generations recorded at the microscale is crucial for estimating the extent and conditions of modern versus fossil serpentinisation in ophiolites. To address matrix bias effects during oxygen isotope analysis by SIMS, we present the first investigation analysing antigorite in the compositional range Mg# = 77.5–99.5 mole %, using a CAMECA IMS–1280 secondary ion mass spectrometer. Spot–to–spot homogeneity is ≤ 0.5‰ (2s) for the new antigorite reference materials. The relative bias between antigorite reference materials with different Mg/Fe ratios is described by a second order polynomial, and a maximum difference in bias of ~ 1.8‰ was measured for Mg# ~ 78 to 100. We observed a bias up to ~ 1.0‰ between lizardite and antigorite attributed to their different crystal structure. Orientation effects up to ~ 1‰ were observed in chrysotile. The new analytical protocol allowed the identification of oxygen isotope zoning up to ~ 7‰ in serpentine minerals from two serpentinites recovered from an area of active serpentinisation in the Samail ophiolite. Furthermore in situ analysis is capable of resolving isotopic heterogeneity that may directly reflect changes in the physical and chemical conditions of multiple serpentinisation events in the Samail ophiolite.

47 OTHER INSTRUMENTATION↗

DOE OSTI2024

Posterior Covariance Matrix Approximations

Here, the Davis equation of state (EOS) is commonly used to model thermodynamic relationships for high explosive (HE) reactants. Typically, the parameters in the EOS are calibrated, with uncertainty, using a Bayesian framework and Markov Chain Monte Carlo (MCMC) methods. However, MCMC methods are computationally expensive, especially for complex models with many parameters. This paper provides a comparison between MCMC and less computationally expensive Variational methods (Variational Bayesian and Hessian Variational Bayesian) for computing the posterior distribution and approximating the posterior covariance matrix based on heterogeneous experimental data. All three methods recover similar posterior distributions and posterior covariance matrices. This study demonstrates that for this EOS parameter calibration application, the assumptions made in the two Variational methods significantly reduce the computational cost but do not substantially change the results compared to MCMC.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Force Reconstruction Using the Inverse of the Mode-Shape Matrix [Book Chapter]

A new technique for force reconstruction is developed. To estimate the externally applied force, this technique sums the weight-scaled acceleration signals, and is referred to as the Sum of Weighted Accelerations Technique (SWAT). To obtain the scalar weights the inverse of the mode shape matrix is used. Application of this technique is illustrated with both numerical calculations using a mass-spring model and experimental data from a structure impacting a rigid barrier.

055001 -- Nuclear Fuels-- Safeguards, Inspection, ↗

DOE OSTI2025

Environmental matrix and moisture influence soil microbial phenotypes in a simplified porous media incubation

Soil moisture and porosity regulate microbial metabolism by influencing factors, such as system chemistry, substrate availability, and soil connectivity. However, accurately representing the soil environment and establishing a tractable microbial community that limits confounding variables is difficult. Here, we use a reduced-complexity microbial consortium grown in a glass bead porous media amended with chitin to test the effects of moisture and a structural matrix on microbial phenotypes. Leveraging metagenomes, metatranscriptomes, metaproteomes, and metabolomes, we saw that our porous media system significantly altered microbial phenotypes compared with the liquid incubations, denoting the importance of incorporating pores and surfaces for understanding microbial phenotypes in soils. These phenotypic shifts were mainly driven by differences in expression of Streptomyces and Ensifer, which included a significant decrease in overall chitin degradation between porous media and liquid. Our findings suggest that the success of Ensifer in porous media is likely related to its ability to repurpose carbon via the glyoxylate shunt amidst a lack of chitin degradation byproducts while potentially using polyhydroxyalkanoate granules as a C source. We also identified traits expressed by Ensifer and others, including motility, stress resistance, and carbon conservation, that likely influence the metabolic profiles observed across treatments. Together, these results demonstrate that porous media incubations promote structure-induced microbial phenotypes and are likely a better proxy for soil conditions than liquid culture systems. Furthermore, they emphasize that microbial phenotypes encompass not only the multi-enzyme pathways involved in metabolism but also include the complex interactions with the environment and other community members.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

SiPM-matrix readout of two-phase argon detectors using electroluminescence in the visible and near infrared range

Proportional electroluminescence (EL) in noble gases is used in two-phase detectors for dark matter searches to record (in the gas phase) the ionization signal induced by particle scattering in the liquid phase. The “standard” EL mechanism is considered to be due to noble gas excimer emission in the vacuum ultraviolet (VUV). In addition, there are two alternative mechanisms, producing light in the visible and near infrared (NIR) ranges. The first is due to bremsstrahlung of electrons scattered on neutral atoms (“neutral bremsstrahlung”, NBrS). The second, responsible for electron avalanche scintillation in the NIR at higher electric fields, is due to transitions between excited atomic states. In this work, we have for the first time demonstrated two alternative techniques of the optical readout of two-phase argon detectors, in the visible and NIR range, using a silicon photomultiplier matrix and electroluminescence due to either neutral bremsstrahlung or avalanche scintillation. The amplitude yield and position resolution were measured for these readout techniques, which allowed to assess the detection threshold for electron and nuclear recoils in two-phase argon detectors for dark matter searches. To the best of our knowledge, this is the first practical application of the NBrS effect in detection science.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2023

Spin density matrix elements in exclusive $\rho ^0$ meson muoproduction

We report on a measurement of Spin Density Matrix Elements (SDMEs) in hard exclusive $\rho ^0$ meson muoproduction at COMPASS using 160 GeV/c polarised μ + and μ – beams impinging on a liquid hydrogen target. The measurement covers the kinematic range 5.0 GeV/c 2 < W < 17.0 GeV/c 2 , 1.0 (GeV/c) 2 < $\mathrm{Q}$ 2 < 10.0 (GeV/c) 2 and 0.01 (GeV/c) 2 < $p_{\textrm{T}}^2$ < 0.5 (GeV/c) 2 . Here, W denotes the mass of the final hadronic system, $\mathrm{Q}$ 2 the virtuality of the exchanged photon, and $p_{\textrm{T}}$ the transverse momentum of the $\rho ^0$ meson with respect to the virtual-photon direction. The measured non-zero SDMEs for the transitions of transversely polarised virtual photons to longitudinally polarised vector mesons ($\gamma ^*_T \rightarrow V^{ }_L$) indicate a violation of s-channel helicity conservation. Additionally, we observe a dominant contribution of natural-parity-exchange transitions and a very small contribution of unnatural-parity-exchange transitions, which is compatible with zero within experimental uncertainties. The results provide important input for modelling Generalised Parton Distributions (GPDs). In particular, they may allow one to evaluate in a model-dependent way the role of parton helicity-flip GPDs in exclusive $\rho ^0$ production.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Matrix element corrections in the PYTHIA8 parton shower in the context of matched simulations at next-to-leading order

We discuss the role of matrix element corrections (MEC) to parton showers in the context of MC@NLO-type matchings for processes that feature unstable resonances, where MEC are liable to result in double-counting issues, and are thus generally not employed. By working with PYTHIA8 , we show that disabling all MEC is actually unnecessary in computations based on the narrow-width approximation, and we propose alternative MEC settings which, while still avoiding double counting, allow one to include hard-recoil effects in the simulations of resonance decays. We illustrate our findings by considering $t\overline{t}$ production at the LHC, and by comparing MadGraph5-aMC@NLO predictions with those of POWHEG-BOX and standalone PYTHIA8 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Improving inference with matrix elements and machine learning

Particle physics processes bring together a high-energy amplitude described by quantum field theory and nonperturbative effects and detector interactions described by complex computer simulations. We review some recently developed multivariate inference techniques that leverage this structure and combine matrix-element information with machine learning. Automated by the MadMiner package, the new techniques have been applied to multiple problems in particle physics, allowing for stronger limits than traditional analysis methods and showing their potential to improve the sensitivity of the LHC legacy measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Sampling-based Sublinear Low-rank Matrix Arithmetic Framework for Dequantizing Quantum Machine Learning

We present an algorithmic framework for quantum-inspired classical algorithms on close-to-low-rank matrices, generalizing the series of results started by Tang’s breakthrough quantum-inspired algorithm for recommendation systems [STOC’19]. Motivated by quantum linear algebra algorithms and the quantum singular value transformation (SVT) framework of Gilyén et al. [STOC’19], we develop classical algorithms for SVT that run in time independent of input dimension, under suitable quantum-inspired sampling assumptions. Our results give compelling evidence that in the corresponding QRAM data structure input model, quantum SVT does not yield exponential quantum speedups. Since the quantum SVT framework generalizes essentially all known techniques for quantum linear algebra, our results, combined with sampling lemmas from previous work, suffice to generalize all prior results about dequantizing quantum machine learning algorithms. In particular, our classical SVT framework recovers and often improves the dequantization results on recommendation systems, principal component analysis, supervised clustering, support vector machines, low-rank regression, and semidefinite program solving. We also give additional dequantization results on low-rank Hamiltonian simulation and discriminant analysis. Our improvements come from identifying the key feature of the quantum-inspired input model that is at the core of all prior quantum-inspired results: ℓ 2 -norm sampling can approximate matrix products in time independent of their dimension. We reduce all our main results to this fact, making our exposition concise, self-contained, and intuitive.

Computer Science↗

DOE OSTI2023

Tiling Framework for Heterogeneous Computing of Matrix based Tiled Algorithms

Tiling matrix operations can improve the load balancing and performance of applications on heterogeneous computing resources. Writing a tile-based algorithm for each operation with a traditional, hand-tuned tiling approach that uses for loops in C/C++ is cumbersome and error prone. Moreover, it must enable and support the heterogeneous memory management of data objects and also explore architecture-supported, native, tiled-data transfer APIs instead of copying the tiled data to continuous memory before the data transfer. The tiling framework provides a tiled data structure for heterogeneous memory mapping and parameterization to a heterogeneous task specification API. We have integrated our tiled framework into MatRIS (Math kernels library using IRIS). IRIS is a heterogeneous run-time framework with a heterogeneous programming model, memory model, and task execution model. Experiments reveal that the tiled framework for BLAS operations has improved the programmability of tiled BLAS and improved performance by ~20% when compared against the traditional method that copies the data to continuous memory locations for heterogeneous computing.

Miniskar, Narasinga Rao↗

DOE OSTI2025

A Study of Performance Portability of Low-bit Fused Matrix-Vector Multiplication Kernels in SYCL

Understanding the causes of performance gaps between a portable programming model and a vendor-specific programming model is important for improving performance portability. This paper studies performance portability of low-bit fused general matrix-vector multiplication kernels in SYCL on vendors’ graphics processing units (GPUs). This work introduces the use case, explains the kernel implementations in detail, evaluates the performance of the CUDA, HIP, and SYCL kernels on datacenter, desktop, and laptop GPUs, and investigates the causes of performance gaps. The results show that loop unrolling, kernel dispatch overhead, and sum reduction contribute to the gaps.

Jin, Zheming [ORNL] (ORCID:000000027197780X)↗

DOE OSTI2021

Time-Resolved X-ray Operando Observations of Lithiation Gradients across the Cathode Matrix and Individual Oxide Particles during Fast Cycling of a Li-Ion Cell

Lithiated transition metal oxides serve as active materials in the positive electrode (cathode) of lithium-ion cells. During electrochemical cycling, lithium ions intercalate and deintercalate into these oxide particles. This behavior causes two types of lithiation gradients to emerge: (i) a bulk gradient across the depth of the cathode matrix (averaged over individual oxide particles) and (ii) a microscopic gradient across the particles themselves, which also depends on their location in the electrode. Here we show how both gradients can be studied using operando X-ray diffraction during 4C charge and 4C discharge. The oxide (de)lithiation is estimated from the unit cell parameters by indexing the X-ray diffraction spectra. By fitting the lithiation profiles with orthogonal polynomials, the bulk gradients across the electrode thickness are quantified. These gradients develop as the current flows through the cell and dissipate during open-circuit and potentiostatic-hold periods. Further details of lithiation dynamics can be obtained through shape analysis of the Bragg peaks. In particular, from electrochemical model simulations, we show that the width and skewness of the (003) peak track (de)lithiation fronts moving across the individual oxide particles.

(003) peak↗

DOE OSTI2024

Unraveling Morphology and Chemistry Dynamics in Fluoroethylene Carbonate Generated Silicon Anode Solid Electrolyte Interphase Across Delithiated and Lithiated States: Relative Cycling Stability Enabled by an Elastomeric Polymer Matrix

The silicon solid electrolyte interphase (SEI) faces cyclical cracking and reconstruction due to the ∼350% volume expansion. Understanding the SEI dynamic morphology and chemistry evolution from delithiated to lithiated states is thereby paramount to engineering a stable Si anode. Fluoroethylene carbonate (FEC) is a preferred additive with widely demonstrated enhancement of the Si cycling. Thus, insights into the dynamics of the FEC-SEI may provide hints toward engineering the Si interface. Herein, complementary ATR-FTIR, AFM, tip IR, and XPS probing reveal the presence of an elastomeric polycarbonate-like matrix in the FEC-generated SEI which is absent from the FEC-free SEI. Adding FEC to the baseline 1 M LiPF 6 in EC:EMC (1:1) electrolyte promotes formation of a thinner and more conformal SEI, and subdues morphology and chemistry changes between consecutive half-cycles. From AFM, morphological stabilization of the FEC-SEI occurs earlier. Furthermore, conventional SEI biproducts such as Li 2 CO 3 and LiEDC appear in reduced quantities in the FEC-SEI implying a reduced quantity of Li-consuming species. The thin polymeric FEC-SEI enables deeper (de)lithiation of silicon. In conclusion, the enhanced mechanical compliance, chemical invariance, and reduced Li inventory consumption of the FEC-SEI are logically the key features underlying the Si cycling enhancement by FEC.

25 ENERGY STORAGE↗

DOE OSTI2020

CSBCOO SPMM: Compressed Sparse Block Coordinate format sparse matrix dense vector multiply benchmark (CSBCOO SPMM) v1.0

The purpose of this microbenchmark is to provide a means to explore different programming methodologies using a simple, but not trivial, mathematical kernel. The kernel is based on the matrix vector product in the MFDn nuclear configuration interaction code's LOBPCG eigensolver. This particular mathematical operation is a key element in the solution of systems of equations and block eigensolvers.

Cook, Brandon↗

DOE OSTI2024

Matrix-based Omnidirectional Pressure Integratio

This software enables the calculation of the pressure field based on experimental velocity measurements utilizing the matrix-based omnidirectional pressure integration algorithm developed by Dr. John Charonko and Dr. Fernando Zigunov, which is several orders of magnitude faster than the current state of the art.

Zigunov, Fernando↗