Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MATRICES”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Hypergraph Random Walks, Laplacians, and Clustering

We propose a flexible framework for clustering hypergraph-structured data based on recently proposed random walks utilizing edge-dependent vertex weights. When incorporating edge-dependent vertex weights (EDVW), a weight is associated with each vertex-hyperedge pair, yielding a weighted incidence matrix of the hypergraph. Such weightings have been utilized in term-document representations of text data sets. We explain how random walks with EDVW serve to construct different hypergraph Laplacian matrices, and then develop a suite of clustering methods that use these incidence matrices and Laplacians for hypergraph clustering. Using 20Newsgroup, U.S. patent, Reuters' Corpus Volume 1, and genetics data sets, we compare the performance of these clustering algorithms experimentally against a variety of existing hypergraph clustering methods. We show that the proposed methods produce higher-quality clusters.

hypergraphs, clustering, laplacian, random walk, M↗
DOE OSTI2021

Structure Perception in 3D Point Clouds

Understanding human perception is critical to the design of effective visualizations. The relative benefits of using 2D versus 3D techniques for data visualization is a complex decision space, with varying levels of uncertainty and disagreement in both the literature and in practice. This study aims to add easily reproducible, empirical evidence on the role of depth cues in perceiving structures or patterns in 3D point clouds. We describe a method to synthesize a 3D point cloud that contains a 3D structure, where 2D projections of the data strongly resemble a Gaussian distribution. We performed a within-subjects structure identification study with 128 participants that compared scatterplot matrices (canonical 2D projections) and 3D scatterplots under three types of motion: rotation, xy-translation, and z-translation. We found that users could consistently identify three separate hidden structures under rotation, while those structures remained hidden in the scatterplot matrices and under translation. This work contributes a set of 3D point clouds that provide definitive examples of 3D patterns perceptible in 3D scatterplots under rotation but imperceptible in 2D scatterplots.

data analysis↗
DOE OSTI2023

A Verified Optimizer for Quantum Circuits

We present VOQC, the first verified optimizer for quantum circuits, written using the Coq proof assistant. Quantum circuits are expressed as programs in a simple, low-level language called SQIR, a small quantum intermediate representation, which is deeply embedded in Coq. Optimizations and other transformations are expressed as Coq functions, which are proved correct with respect to a semantics of SQIR programs. SQIR programs denote complex-valued matrices, as is standard in quantum computation, but we treat matrices symbolically to reason about programs that use an arbitrary number of quantum bits. SQIR’s careful design and our provided automation make it possible to write and verify a broad range of optimizations in VOQC, including full-circuit transformations from cutting-edge optimizers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Fully Homomorphic Encryption

This code implements a Fully Homomorphic Encryption (FHE) system, enabling secure computation on encrypted data without requiring decryption. It supports encryption, decryption, and homomorphic operations like matrix multiplication and addition. This code is adaptable for integrating FHE into linear-time invariant (LTI) systems, including digital control and filtering. With proper configuration from subject matter expertise, encrypted system parameters and signals can be manipulated to perform tasks like state updates, output calculations, and convolution in the encrypted domain. By preserving the structure of LTI systems while ensuring privacy, the framework facilitates secure applications in areas such as autonomous systems, signal processing, and industrial automation. The code initializes the encryption system using parameters provided in the env dictionary. These parameters include the ciphertext modulus, key dimension, plaintext fixed-point scaling factor, and noise bound. During initialization, a secret key is generated, which is essential for encrypting and decrypting data securely. The modular design allows users to tailor these parameters to specific use cases or security requirements. The code implements multiple cryptographic schemes. The learning with errors (LWE) encryption method encodes cleartext message to their plaintext fixed-point representation then encrypted into ciphertext space with additive noise. This noise ensures the security of the scheme, relying on the computational hardness of the LWE problem. The code also includes the Gentry-Sahai-Waters (GSW) scheme based off the LWE problem. Homomorphic matrix multiplication is performed between the LWE and GSW to encrypted data. This is achieved using a decomposition function on the LWE ciphertext during the multiplication operation. For higher-dimensional data, the code includes a method to encrypt entire matrices (GSWMat) using GSW encryption. These encrypted matrices can then be used for homomorphic matrix multiplications (MatMult). The decryption function uses the secret key to recover the original plaintext, removing the added noise and scaling that was originally applied during encryption.

Lois, Roberts [Idaho National Laboratory (INL), Id↗
DOE OSTI2025

Matrix-based Parallel Redistribution

MatRed is a parallel redistribution tool for HPC applications. It provides a simple approach that only requires a few relation matrices between entities to build redistribution matrices in parallel simulation codes. In particular, MatRed is well-suited for simulation codes based on finite element/volume methods.

Kalchev, DelyanZ [Lawrence Livermore National Labo↗
DOE OSTI2022

Dataset of Simulated Intracardiac Transmembrane Voltage Recordings and ECG Signals

The dataset consists of pairs of computationally simulated intracardiac transmembrane voltage recordings and ECG signals. In total, 16140 organ-level simulations were conducted to create this dataset, using a range of cardiac geometries and physiological parameters. Simulations were performed at LLNL's Lassen supercomputer, concurrently utilizing 4 GPUs and 40 CPU cores. Each simulation produced pairs of 500ms-by-10 ECG signals and 500ms-by-75 transmembrane voltage signals. For convenience, the signals are concatenated and saved as matrices. Each of these matrices is then stored as a numpy array. See the documentation for further details.

Landajuela, Mikel↗
DOE OSTI2021

CASTELO: clustered atom subtypes aided lead optimization—a combined machine learning and molecular modeling method

Background: Drug discovery is a multi-stage process that comprises two costly major steps: pre-clinical research and clinical trials. Among its stages, lead optimization easily consumes more than half of the pre-clinical budget. We propose a combined machine learning and molecular modeling approach that partially automates lead optimization workflow in silico, providing suggestions for modification hot spots. Results: The initial data collection is achieved with physics-based molecular dynamics simulation. Contact matrices are calculated as the preliminary features extracted from the simulations. To take advantage of the temporal information from the simulations, we enhanced contact matrices data with temporal dynamism representation, which are then modeled with unsupervised convolutional variational autoencoder (CVAE). Finally, conventional and CVAE-based clustering methods are compared with metrics to rank the submolecular structures and propose potential candidates for lead optimization. Conclusion: With no need for extensive structure-activity data, our method provides new hints for drug modification hotspots which can be used to improve drug potency and reduce the lead optimization time. It can potentially become a valuable tool for medicinal chemists.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Rotational symmetry relation for efficient response function generation in the coarse mesh transport method COMET

The coarse mesh transport code COMET is a continuous energy hybrid stochastic-deterministic neutronics solver with high fidelity and formidable computation speed in solving reactor core problems. Its method is based on the incident flux expansion theory. In this work, we take advantage of the local geometric symmetry in many reactor cores lattices (e.g., fuel lattices and reflector blocks) to develop relations among the flux response expansion coefficients for symmetric surfaces to further improve the computational efficiency of the COMET response function generation tool (method). This is done by a rigorous derivation of the transformation matrices for the angular and spatial expansion moments resulting from a rotation of a coarse mesh by an arbitrary angle. The relations for the response coefficients for the symmetric surfaces can be then written as the Kronecker product of those transformation matrices. The method is implemented into COMET and tested on two advanced high temperature reactor (AHTR) full-length single assembly benchmark problems. The COMET results using the response function library based on the symmetry relations were compared to those using the library directly generated by continuous energy Monte Carlo for all surfaces. It was found that the eigenvalues and stripe-wise fission densities using the two libraries are in statistical agreement as expected. This indicates that the new method maintains the high fidelity of the original COMET method while improving the computational efficiency in the response function generation by 270% to 400%, depending on the local geometric symmetry. This method also reduces the size of the response function library by the same magnitude (270% to 400%). (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Efficient ptychographic phase retrieval via a matrix-free Levenberg-Marquardt algorithm

The phase retrieval problem, where one aims to recover a complex-valued image from far-field intensity measurements, is a classic problem encountered in a range of imaging applications. Modern phase retrieval approaches usually rely on gradient descent methods in a nonlinear minimization framework. Calculating closed-form gradients for use in these methods is tedious work, and formulating second order derivatives is even more laborious. Additionally, second order techniques often require the storage and inversion of large matrices of partial derivatives, with memory requirements that can be prohibitive for data-rich imaging modalities. We use a reverse-mode automatic differentiation (AD) framework to implement an efficient matrix-free version of the Levenberg-Marquardt (LM) algorithm, a longstanding method that finds popular use in nonlinear least-square minimization problems but which has seen little use in phase retrieval. Furthermore, we extend the basic LM algorithm so that it can be applied for more general constrained optimization problems (including phase retrieval problems) beyond just the least-square applications. Since we use AD, we only need to specify the physics-based forward model for a specific imaging application; the first and second-order derivative terms are calculated automatically through matrix-vector products, without explicitly forming the large Jacobian or Gauss-Newton matrices typically required for the LM method. We demonstrate that this algorithm can be used to solve both the unconstrained ptychographic object retrieval problem and the constrained “blind” ptychographic object and probe retrieval problems, under the popular Gaussian noise model as well as the Poisson noise model. We compare this algorithm to state-of-the-art first order ptychographic reconstruction methods to demonstrate empirically that this method outperforms best-in-class first-order methods: it provides excellent convergence guarantees with (in many cases) a superlinear rate of convergence, all with a computational cost comparable to, or lower than, the tested first-order algorithms.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Analysis, identification and confirmation of synthetic opioids using chloroformate chemistry: Retrospective detection of fentanyl and acetylfentanyl in urine and plasma samples by EI-GC-MS and HR-LC-MS

Electron Impact Gas Chromatography-Mass Spectrometry (EI-GC-MS) and High Resolution Liquid Chromatography-Mass Spectrometry (HR-LC-MS) have been used in the analysis of products arising from the trichloroethoxycarbonylation of fentanyl and acetylfentanyl in urine and plasma matrices. The method involves the initial extraction of both synthetic opioids separately from the matrices followed by detection of the unique products that arise from their reaction with 2,2,2-trichloroethoxycarbonyl chloride (Troc-Cl), namely Troc-norfentanyl and Troc-noracetylfentanyl. The optimized protocol was successfully evaluated for its efficacy at detecting these species formed from fentanyl and acetylfentanyl when present at low and high levels in urine (fentanyl: 5 and 10 ng/mL and acetylfentanyl: 20 and 100 ng/mL) and plasma (fentanyl: 10 and 20 ng/mL and acetylfentanyl: 50 and 200 ng/mL), values that reflect levels reported in overdose victims. The HR-LC-MS method’s LOQ (limit of quantitation) for the Troc-norfentanyl and Troc-noracetylfentanyl products was determined to be ~10 ng/mL for both species. Even though the superiority in the detection of these species by HR-LC-MS over EI-GC-MS, the latter method proved to be important in the detection of the second product from the reaction, namely 2-phenylethyl chloride that is crucial in the determination of the original opioid. This observation highlights the importance of using complimentary analytical techniques in the analysis of a sample, whether biological or environmental in nature. The method herein serves as a complementary, qualitative confirmation for the presence of a fentanyl in collected urine, plasma and by extension other biological samples amenable to the common extraction procedures described for opioid analysis. More importantly, the method’s main strength comes from its ability to react with unknown fentanyls to yield products that can be not only detected by EI-GC-MS and HR-LC-MS but can then be used to retrospectively identify an unknown fentanyl.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Daily water stable isotopes, transpiration, and matrix potential data for an aspen and engelmann stand in the East River Watershed (version 2)

We provide daily stable isotope (2H & 18O) ratios in soil water and xylem (plant stem) water, as well as the sap flow (transpiration) and the soil's matric potential at a forested site near Gothic, Colorado, in the East River catchment. We measured the stable isotopic composition of the transpiration and the daily transpiration flux sum of three aspen and three engelmann spruce. In both forest stands, we installed a soil profile and measured the soil matric potential at 15, 30, and 60 cm depth as well as the stable isotopes of soil pore water at 5, 10, 30, 60, and 90 cm depths. All isotope measurements were done in situ via vapor probes connected to a cavity ring down spectrometer (Picarro L1240i).We further report the daily meteorological data observed at billy barr near our study site. We also provide for each tree the relative share of root water uptake derived from the isotope measurements via a Bayesian mixing model (MixSIAR).The daily data is provided as a time series in "Iso_MP_Sap_DataDaily_ESSDiveUpload.csv" and the units are provided in "dd.csv"; the location of the instrumented trees and soil profiles are given as latitude and longitude coordinates saved as CSV and KMZ files; and a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata.The data was gathered to investigate the short-term changes of the water sources (i.e., variation of root water uptake from different soil depths) of the studied subalpine trees.Update 07/16/2025: The relative and absolute plant water uptake depths were grouped to ensure that the MixSIAR model was applied with endmembers that differed in their d2H value by at least 3 permill and at least 1 permill for d18O. Whenever the difference between observed d2H values for two or more probes at neighboring depths was less than 3 permill, we used the average value for the source water endmember. For days at which probes that were not next to each other measurements did not differ at least 3 permill, the average of all probes between these two depths was used as the water source endmember representing the depth range between these two probes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Contrast Matching Biopolymers: A SANS-Based Approach to Structural Characterization of Chitosan

Immobilizing enzymes in polysaccharide-based matrices has been shown to improve both their stability and their catalytic efficiency. Among available materials, chitosan was selected for its abundance, biocompatibility, and versatility in a range of applications. While the broader aim of this work is to study chitosan as a matrix for immobilizing carbonic anhydrase (CA) for CO₂ capture, this study focused on determining the contrast match point (CMP) of hydrogenated and deuterated chitosan in D₂O:H₂O mixtures to enable future structural investigations using small-angle neutron scattering (SANS). Previous work successfully synthesized deuterated chitosan from Rhizopus oryzae mycelia and assessed the degree of deuteration using Fourier-transform infrared spectroscopy (FTIR) and nuclear magnetic resonance (NMR). Building on this foundation, the current study determined the CMPs of both H-chitosan and D-chitosan by performing SANS measurements across a range of D₂O:H₂O solvent contrasts. The extracted CMPs for both forms correlated well with their respective calculated scattering length densities (SLDs), validating the experimental approach. These results establish a critical reference for selective contrast matching in future SANS studies of carbonic anhydrase immobilized in chitosan matrices, enabling detailed structural analysis of enzyme–matrix interactions.

Kosgallana, Chathurika J [ORNL] (ORCID:00000002851↗
DOE OSTI2022

Electron beam treatment for the removal of 1,4-dioxane in water and wastewater

Electron beam (e-beam) treatment uses accelerated electrons to form oxidizing and reducing radicals when applied to water without the use of external chemicals. In this study, electron beam treatment was used to degrade 1,4-dioxane in several water matrices. Removal improved in the progressively cleaner water matrices and removals as high as 94% to 99% were observed at a dose of 2.3 kGy in secondary effluent. 1,4-dioxane removal was confirmed to be primarily through hydroxyl radical oxidation. The calculated electrical energy per order was found to be 0.53, 0.26, and 0.08 kWh/m 3 /order for secondary effluent (Avg. total organic carbon (TOC) 9.25 mg/L), granular activated carbon effluent (TOC 3.46 mg/L), and ultrapure water, respectively, with a 70% generation and transfer efficiency applied.

1,4-dioxane↗
DOE OSTI2020

Tailoring Charge Transport and Magnetism in Complex Half-Heusler/Full-Heusler Nanocomposites (Final Report)

This project exploits (1) the facile co-crystallization of the half-Heusler (HH) and full-Heusler (FH) structures, (2) the high coherency and stability of the HH/FH interfaces, and (3) the solid-state inter-conversion between HH and FH structures through local atomic diffusion, to investigate, for selected compositions of the HH and FH phases, the effect of size, dimensionality, dispersion and volume fraction of FH nanostructures on the electronic, phonon and magnetic behavior of the resulting HH/FH nanocomposites. Our goal is to identify and control key material and growth parameters governing phase formation, microstructural evolution, and to understand the mechanism by which these changes of the internal structure modify the material’s performance. To elucidate the mechanism by which the FH nanostructures regulate electronic charge transport and induce ferromagnetism within non-magnetic and magnetic semiconducting HH matrices, suitable series of HH/0D-FH and HH/2D-FH nanocomposites will be fabricated using synthesis conditions suggested by theoretical calculations of the formation energy and stability of various FH phases (varying M) in HH matrices with known chemical composition. The structure of the HH/FH interfaces (size, chemical composition, volume fraction and dispersion) will be investigated along with the evaluation of the electronic (carrier density, mobility, thermopower, electrical conductivity, effective mass), thermal and magnetic (susceptibility, magnetization etc.) properties of the resulting HH/FH nanocomposites.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Functional Nanomaterials with Tunable Plasmonic and Emittive Properties (Final Project Report)

The goal of the project is on understanding the principles of non-covalent assembling of tunable functional materials with various nanocrystals in order to control light-matter interactions by changing localized interfacial refractive properties, coupling parameters, or the ionic environment of surrounding media. We explore novel soft functional matrices with controlled physical properties via the application of light, mechanical stresses, ionic environment, and electrical potential. These functional matrices is used in conjunction with novel anisotropic inorganic nanostructures such as noble metal/metal alloy anisotropic nanostructures with complex optical properties (multiple LSPR modes) and quantum dots with light-emissive or absorbing properties in a broad spectral range from the UV to the NIR for efficient solar energy absorption and light conversion.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Crystalline Silicotitanate Batch Contact Testing with Ba, Ca, Pb, and Sr

Washington River Protection Solutions is working to support initial production of immobilized low-activity waste (LAW) by feeding Hanford tank supernate from tank farms to the Hanford Waste Treatment and Immobilization Plant (WTP) LAW Facility. This goal incorporates the design of a Tank Side Cesium Removal system, which filters tank waste supernate to remove suspended solids and then removes Cs by processing it through crystalline silicotitanate (CST) ion exchange media manufactured by Honeywell UOP, LLC. The 137 Cs-depleted product is intended to be sent to the WTP for vitrification. Processing of actual tank waste supernate showed effectively complete uptake of Sr and Ba by CST and significant uptake of Ca. Further, Campbell et al. (2019) analyzed CST post-column testing and found significant (>1E-2 mmoles/g) uptake of Ca and Pb along with some Ba, Cd, Fe, Sr, and U. This led to concern that selected metals, particularly the +2 cations, Ca, Sr, Ba, and Pb may be consuming Cs exchange sites and possibly reducing CST capacity for Cs. Exchange of +2 cations was assumed to be associated with the M(OH) + species for the metal (M) ion in the caustic solution. A series of batch contact testing was conducted to evaluate the exchange behavior of Ba, Ca, Pb, and Sr onto CST. The Cs exchange behavior was also tested as a benchmark for direct comparisons. The CST was provided in the sodium form by Honeywell UOP, as IONSIV TM R9140-B, Lot 2002009604, 18 x 50 mesh. A <30-mesh aliquot was collected to match the sieve fraction expected for use in upcoming small column test configurations. Kinetic exchange rate and isotherms were measured at metal concentrations benchmarked from the AP-107 tank waste feed composition and as limited by the metal solubility in the alkaline solution. Two simplified matrices were tested: 1) 1.0 M NaOH/4.6 M NaNO 3 and 2) 0.1 M NaOH/5.5 M NaNO 3 ; these matrices represented the expected 5.6 M Na concentration of process feed and served to address the hydroxide concentration effect on exchange behavior.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

LAW Simulant Recipes for Evaluation of Real-Time, In-Line Monitoring Instruments

Statistically designed recipes were generated for evaluation of in-line instrumentation for the Real-Time In-Line Monitoring (RTIM) program in support of the Hanford Direct-Feed Low Activity Waste (DFLAW) program. The composition matrices were prepared based on Low Activity Waste (LAW) feed vector, melter feed, and Effluent Management Facility (EMF) feed information provided by the flowsheet integration modelling for the DFLAW portion of the Waste Treatment Plant (WTP) mission. These composition matrices were converted to recipes for the LAW feed using the approaches for LAW simulant recipes developed by Savannah River National Laboratory (SRNL). Recipes for the glass forming chemicals (GFCs) were developed based on expected mineral forms to be used to for each element to be added to the recipe.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Randomized Cholesky Preconditioning for Graph Partitioning Applications

Graph partitioning has emerged as an area of interest due to its use in various applications in computational research. One way to partition a graph is to solve for the eigenvectors of the corresponding graph Laplacian matrix. This project focuses on the eigensolver LOBPCG and the evaluation of a new preconditioner: Randomized Cholesky Factorization (rchol). This proconditioner was tested for its speed and accuracy against other well-known preconditioners for the method. After experiments were run on several known test matrices, rchol appears to be a better preconditioner for structured matrices. This research was sponsored by National Nuclear Security Administration Minority Serving Institutions Internship Program (NNSA-MSIIP) and completed at host facility Sandia National Laboratories. As such, after discussion of the research project itself, this report contains a brief reflection on experience gained as a result of participating in the NNSA-MSIIP.

97 MATHEMATICS AND COMPUTING↗