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 469 records · Page 26

metabCombiner 2.0: Disparate Multi-Dataset Feature Alignment for LC-MS Metabolomics

Liquid chromatography–high-resolution mass spectrometry (LC-HRMS), as applied to untargeted metabolomics, enables the simultaneous detection of thousands of small molecules, generating complex datasets. Alignment is a crucial step in data processing pipelines, whereby LC-MS features derived from common ions are assembled into a unified matrix amenable to further analysis. Variability in the analytical factors that influence liquid chromatography separations complicates data alignment. This is prominent when aligning data acquired in different laboratories, generated using non-identical instruments, or between batches from large-scale studies. Previously, we developed metabCombiner for aligning disparately acquired LC-MS metabolomics datasets. Here, we report significant upgrades to metabCombiner that enable the stepwise alignment of multiple untargeted LC-MS metabolomics datasets, facilitating inter-laboratory reproducibility studies. To accomplish this, a “primary” feature list is used as a template for matching compounds in “target” feature lists. We demonstrate this workflow by aligning four lipidomics datasets from core laboratories generated using each institution’s in-house LC-MS instrumentation and methods. We also introduce batchCombine, an application of the metabCombiner framework for aligning experiments composed of multiple batches. metabCombiner is available as an R package on Github and Bioconductor, along with a new online version implemented as an R Shiny App.

97 MATHEMATICS AND COMPUTING↗

Mass spectrometry imaging of natural carbonyl products directly from agar-based microbial interactions using 4-APEBA derivatization

Aliphatic carboxylic acids, aldehydes, and ketones play diverse roles in microbial adaptation to their microenvironment, from excretion as toxins to adaptive metabolites for membrane fluidity. However, the spatial distribution of these molecules throughout biofilms and how microbes in these environments exchange these molecules remain elusive for many of these bioactive species due to inefficient molecular imaging strategies. Herein, we apply on-tissue chemical derivatization (OTCD) using 4-(2-((4-bromophenethyl)dimethylammonio)ethoxy)benzenaminium dibromide (4-APEBA) on a co-culture of a soil bacterium (Bacillus subtilis NCIB 3610) and fungus (Fusarium sp. DS 682) grown on agar as our model system. Using matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI), we spatially resolved more than 300 different metabolites containing carbonyl groups within this model system. Various spatial patterns are observable in these species, which indicate possible extracellular or intercellular processes of the metabolites and their up- or downregulation during microbial interaction. The unique chemistry of our approach allowed us to bring additional confidence in accurate carbonyl identification, especially when multiple isomeric candidates were possible, and this provided the ability to generate hypotheses about the potential role of some aliphatic carbonyls in this B. subtilis/Fusarium sp. interaction. The results shown here demonstrate the utility of 4-ABEBA-based OTCD MALDI-MSI in probing interkingdom interactions directly from microbial co-cultures, and these methods will enable future microbial interaction studies with expanded metabolic coverage.

4-APEBA↗

Iodine Sorption on AgZ-PAN

The Department of Energy (DOE) Nuclear Fuel Cycle and Supply Chain Program Material Recovery and Waste Form Development (MRWFD) Campaign Off-gas Sigma Team has performed research and development on iodine control and iodine waste forms for the past several years. This research and development has included iodine adsorption tests using a laboratory-scale test system containing multiple-segmented fixed beds of iodine adsorbents. Iodine capture performance has been tested using (a) non-radioactive synthetic gas mixtures blended from air, nitrogen, NO, NO 2 , water, to simulate dissolver off-gas (DOG) and vessel off-gas (VOG) streams from used nuclear fuel aqueous reprocessing, (b) inorganic and organic iodine species, and (c) different sorbents. An initial test of silver-functionalized mordenite in polyacrylonitrile matrix (AgZ-PAN) was conducted to assess if this sorbent, initially designed for selective sorption of Xe that evolves into used fuel reprocessing off-gas streams, could also be a candidate for iodine sorption. This test was conducted using the same deep-bed test system and methodologies used in many prior iodine sorption tests using other sorbents. Diatomic iodine (I2) was used as the target iodide in a gas stream containing air, NO, NO 2 , and moisture designed to simulate an aqueous reprocessing dissolver off-gas (DOG) stream.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Strain hardening magnesium-silicate-hydrate composites (SHMSHC) reinforced with short and randomly oriented polyvinyl alcohol microfibers

Magnesium-silicate-hydrate (M-S-H) cement has emerged as an alternative binder with potentially lower energy requirements and emissions. Due to the intrinsically low pH values of its matrix, however, normal steel reinforcement is not appropriate for M-S-H system. To toughen the matrix and to overcome the brittle nature of the material, a new strain hardening magnesium-silicate-hydrate composite (SHMSHC) is developed for the first time by incorporating 2 vol% short and randomly oriented polyvinyl alcohol (PVA) microfibers. The resulting SHMSHC exhibits significant strain hardening with a tensile strain capacity of more than 3%, a compressive strength beyond 50 MPa and a tensile strength of around 3 MPa. Remarkably, saturated multiple cracking with a tight crack width less than 10 μm invisible to the human naked eyes was observed. The newly developed SHMSHC addresses the fundamental challenge of M-S-H binder and is expected to widen possible application areas of M-S-H system.

36 MATERIALS SCIENCE↗

A differentiable perturbation-based weak lensing shear estimator

Upcoming imaging surveys will use weak gravitational lensing to study the large-scale structure of the Universe, demanding sub-per cent accuracy for precise cosmic shear measurements. We present a new differentiable implementation of our perturbation-based shear estimator (fpfs), using jax, which is publicly available as part of a new suite of analytic shear algorithms called anacal. This code can analytically calibrate the shear response of any non-linear observable constructed with the fpfs shapelets and detection modes utilizing autodifferentiation (ad), generalizing the formalism to include a family of shear estimators with corrections for detection and selection biases. Using the ad capability of jax, it calculates the full Hessian matrix of the non-linear observables, which improves the previously presented second-order noise bias correction in the shear estimation. As an illustration of the power of the new anacal framework, we optimize the effective galaxy number density in the space of the generalized shear estimators using an LSST-like galaxy image simulation for the 10 yr LSST. For the generic shear estimator, the magnitude of the multiplicative bias |m| is below 3 × 10 –3 (99.7 per cent confidence interval), and the effective galaxy number density is improved by 5 per cent. We also discuss some planned future additions to the anacal software suite to extend its applicability beyond the fpfs measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

VAN-DAMME: GPU-accelerated and symmetry-assisted quantum optimal control of multi-qubit systems

We present an open-source software package, VAN-DAMME (Versatile Approaches to Numerically Design, Accelerate, and Manipulate Magnetic Excitations), for massively-parallelized quantum optimal control (QOC) calculations of multi-qubit systems. To enable large QOC calculations, the VAN-DAMME software package utilizes symmetry-based techniques with custom GPU-enhanced algorithms. This combined approach allows for the simultaneous computation of hundreds of matrix exponential propagators that efficiently leverage the intra-GPU parallelism found in high-performance GPUs. In addition, to maximize the computational efficiency of the VAN-DAMME code, we carried out several extensive tests on data layout, computational complexity, memory requirements, and performance. These extensive analyses allowed us to develop computationally efficient approaches for evaluating complex-valued matrix exponential propagators based on Padé approximants. To assess the computational performance of our GPU-accelerated VAN-DAMME code, we carried out QOC calculations of systems containing 10 - 15 qubits, which showed that our GPU implementation is 18.4× faster than the corresponding CPU implementation. Our GPU-accelerated enhancements allow efficient calculations of multi-qubit systems, which can be used for the efficient implementation of QOC applications across multiple domains.

97 MATHEMATICS AND COMPUTING↗

A non-isothermal breakage-damage model for plastic-bonded granular materials incorporating temperature, pressure, and rate dependencies

Plastic-bonded granular materials (PBM) are widely used in industrial sectors, including building construction, abrasive applications, and defense applications such as plastic-bonded explosives. The mechanical behavior of PBM is highly nonlinear, irreversible, rate dependent, and temperature sensitive governed by various micromechanical attributions such as grain crushing and binder damage. This paper presents a thermodynamically consistent, microstructure-informed constitutive model to capture these characteristic behaviors of PBM. Key features of the model include a breakage internal variable to upscale the grain-scale information to the continuum level and to predict grain size evolution under mechanical loading. In addition, a damage internal state variable is introduced to account for the damage, deterioration, and debonding of the binder matrix upon loading. Temperature is taken as a fundamental external state variable to handle non-isothermal loading paths. The proposed model is able to capture with good accuracy several important aspects of the mechanical properties of PBM, such as pressure-dependent elasticity, pressure-dependent yield strength, brittle-to-ductile transition, temperature dependency, and rate dependency in the post-yielding regime. Furthermore, the model is validated against multiple published datasets obtained from confined and unconfined compression tests, covering various PBM compositions, confining pressures, temperatures, and strain rates.

Breakage↗

Discrete Fracture Network Modeling to Estimate Upscaled Parameters for the Topopah Spring, Lava Flow, and Tiva Canyon Aquifers at Pahute Mesa, Nevada National Security Site

This report describes the results of Discrete Fracture Network (DFN) simulations for the Topopah Spring Aquifer (TSA), Lava Flow Aquifer, and Tiva Canyon Aquifer (TCA), at Pahute Mesa on the Nevada National Security Site (NNSS), formerly the Nevada Test Site. The research focuses on calculating upscaled groundwater flow and contaminant transport parameters using DFNs generated according to fracture characteristics observed in the TSA, LFA and TCA at Pahute Mesa. The highly fractured and heterogeneous nature of these aquifers makes them candidates for stochastic DFN modeling of radionuclide transport on a small scale with subsequent upscaling. One hundred independent DFN realizations are generated for each aquifer, and the upscaled parameters for continuum simulations of subsurface flow and transport in fractured media at Pahute Mesa are calculated. Our goal is to implement a modeling approach that can translate parameters to larger-scale models that account for local-scale flow and transport processes, such as channelization of flow and transport along a few well connected, large fractures. Additionally, to simulate advective and advective-diffusive transport through the fracture networks, the Time Domain Random Walk (TDRW) approach is applied to account for matrix diffusion into a finite half-space. Moreover, a novel approach to calculate dynamic (active) fracture surface area to reflect flow channeling is implemented. This work will improve the representation of radionuclide transport processes in largescale, regulatory-focused models by providing estimates of hard-to-measure flow and contaminant transport parameters at large scales. In this report, we (1) show recent results of flow and transport simulations on multiple DFN realizations of the TSA, LFA, TCA; (2) discuss the resulting distributions of estimated upscaled parameters; (3) describe the estimation of upscaled parameters for an equivalent parallel-plate continuum model and (4) present a comparison between simulated transport from the equivalent continuum model and an actual DFN.

54 ENVIRONMENTAL SCIENCES↗

A variational mimetic finite difference method for elliptic interface problems on non-matching polytopal meshes with geometric interface inconsistencies

A new variational mimetic finite difference method for elliptic interface problems with perfect and imperfect thermal contacts on non-matching polytopal meshes with geometric interface inconsistencies is developed and analyzed theoretically and numerically. The method is defined on multiple non-matching submeshes with gaps and overlaps along their interfaces. The discrete equations are derived from a minimization problem for the augmented Dirichlet functional. For a perfect thermal contact, the functional uses a modified mimetic gradient with extended stencil which couples unknowns from both sides of an interface, as well as penalty terms to enforce weak continuity of temperature across the interface. The method leads to a symmetric positive definite matrix for any scaling of the penalty terms. For an imperfect thermal contact, the Dirichlet functional is supplemented with a quadratic jump term along the interface related to the interface thermal resistance. We prove that the method conserves the total heat flux across each interface. In conclusion, the obtained results are verified with numerical experiments showing convergence in the discrete L 2 and L ∞ norms.

97 MATHEMATICS AND COMPUTING↗

FEA-aided investigation of the effective thermal conductivity in a medium with embedded spheres

For multiple applications in nuclear energy, the ability to accurately represent material behavior with a simplified model is important to facilitate practical engineering-scale simulations. Here, we focus on the homogenized thermal response of a medium containing spherical inclusions, similar to a fuel form (compact or pebble) containing (TRISO) particles. A review on effective thermal conductivity (ETC) modeling is performed considering a random distribution of mono-sized spherical inclusions in a continuous matrix, with a primary focus on the analytical models. Finite element simulations are performed to evaluate each analytical model in varying material conditions. The model predictions are compared with the expected results obtained from the finite element predictions in addition to the Wiener and Hashin-Shtrikman bounds. Lastly, we included a practical discussion on the homogenization applied to a TRISO fuel pebble.

42 ENGINEERING↗

Carbon Additive Manufacturing with a Near-Replica “Green-to-Brown” Transformation

Nanocomposites containing nanoscale materials offer exciting opportunities to encode nanoscale features into macroscale dimensions, which produces unprecedented impact in material design and application. However, conventional methods cannot process nanocomposites with a high particle loading, as well as nanocomposites with the ability to be tailored at multiple scales. A composite architected mesoscale process strategy that brings particle loading nanoscale materials combined with multiscale features including nanoscale manipulation, mesoscale architecture, and macroscale formation to create spatially programmed nanocomposites with high particle loading and multiscale tailorability is reported. The process features a low-shrinking (<10%) “green-to-brown” transformation, making a near-geometric replica of the 3D design to produce a “brown” part with full nanomaterials to allow further matrix infill. Here, this demonstration includes additively manufactured carbon nanocomposites containing carbon nanotubes (CNTs) and thermoset epoxy, leading to multiscale CNTs tailorability, performance improvement, and 3D complex geometry feasibility. The process can produce nanomaterial-assembled architectures with 3D geometry and multiscale features and can incorporate a wide range of matrix materials, such as polymers, metals, and ceramics, to fabricate nanocomposites for new device structures and applications.

36 MATERIALS SCIENCE↗

Perspectives on Quenching and Tempering 4340 Steel

Steels are ubiquitous due to their affordability and the landscape of useful properties that can be generated for engineering applications. But to further expand the performance envelope, one must be able to understand and control microstructure development by alloying and processing. In this paper we use multiscale, advanced characterization to better understand the structural and chemical evolution of AISI 4340 steel after quenching and tempering (Q&T), including the role of quench rate and short-time, isothermal tempering below 573 K (300 °C), with an emphasis on carbide formation. We compare the microstructure and/or property changes produced by conventional tempering to those produced by higher temperature, short-time “rapid” tempering. We underscore that no single characterization technique can fully capture the subtle microstructure changes like carbon redistribution, transition carbide and/or cementite formation, and retained austenite decomposition that occur during Q&T. Only the use of multiple techniques begins to unravel these complexities. After controlled fast or slow quenching, η transition carbides clearly exist in the microstructure, likely associated with autotempering of this high martensite start temperature (Ms) steel. Isothermal tempering below 598 K (325 °C) results in the relief of carbon supersaturation in the martensite, primarily by the formation of η transition carbides that exhibit a range of carbon levels, seemingly without substitutional element partitioning between the carbide and matrix phases. Hägg transition carbide is present between 300 °C and 325 °C. After conventional tempering at or above 598 K (325 °C) for 2 h, cementite is predominant, but small amounts of cementite are also present in other conditions, even after quenching. Previous work has indicated that silicon (Si) and substitutional elements partition between the cementite, which initially forms under paraequilibrium conditions, and the matrix. Phosphorous (P) may also be preferentially located at cementite/matrix interfaces after high temperature tempering. Slower quench rates result in greater amounts of retained austenite compared to those after fast quenching, which we attribute to increased austenite stability resulting from “autopartitioning”. Rapid, high temperature tempering is also found to diminish tempered martensite embrittlement (TME) believed to be associated with the extent of austenite decomposition, resulting in mechanical properties not attainable by conventional tempering, which may have important implications with respect to industrial heat treatment processes like induction tempering. Controlling the amount and stability of retained austenite is not only relevant to the properties of Q&T steels, but also next-generation advanced high strength steels (AHSS) with austenite/martensite mixtures.

36 MATERIALS SCIENCE↗

Computing the QRPA level density with the finite amplitude method

Here, we describe a new algorithm to calculate the vibrational nuclear level density of an atomic nucleus. Fictitious perturbation operators that probe the response of the system are generated by drawing their matrix elements from some probability distribution function. We use the Finite Amplitude Method to explicitly compute the response for each such sample. With the help of the Kernel Polynomial Method, we build an estimator of the vibrational level density and provide the upper bound of the relative error in the limit of infinitely many random samples. The new algorithm can give accurate estimates of the vibrational level density. Since it is based on drawing multiple samples of perturbation operators, its computational implementation is naturally parallel and scales like the number of available processing units.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Advances in Cell Wall Matrix Research with a Focus on Mixed-Linkage Glucan

Mixed β(1,3;1,4)-linkage glucan (MLG) is commonly found in the monocot lineage, at particularly high levels in the Poaceae family, but also in the evolutionally distant genus, Equisetum. MLG has several properties that make it unique from other plant cell wall polysaccharides. It consists of β1,4-linked polymers of glucose interspersed with β1,3-linkages, but the presence of β1,3-linkages provides quite different physical properties compared to its closest form of the cell wall component, cellulose. The mechanisms of MLG biosynthesis have been investigated to understand whether single or multiple enzymes are required to build mixed linkages in the glucan chain. Currently, MLG synthesis by a single enzyme is supported by mutagenesis analyses of cellulose synthase-like F6, the major MLG synthase, but further investigation is needed to gather mechanistic insights. Because of transient accumulation of MLG in elongating cells and vegetative tissues, several hypotheses have been proposed to explain the role of MLG in the plant cell wall. Studies have been carried out to identify gene expression regulators during development and light cycles as well as enzymes involved in MLG organization in the cell wall. A role of MLG as a storage molecule in grains is evident, but the role of MLG in vegetative tissues is still not well understood. Characterization of a cell wall component is difficult due to the complex heterogeneity of the plant cell wall. Furthermore, as detailed in this review, recent exciting research has made significant impacts in the understanding of MLG biology in plants.

Transglucosylase/hydrolase (XTH)↗

Single-Molecule Kinetics of Styrene Hydrogenation on Silica-Supported Vanadium: The Role of Disorder for Single-Atom Catalysts

We report a theoretical approach for the study of supported atom catalysis is developed based on recent advances in the study of single-molecule kinetics. This view is particularly useful in exhibiting the role of disorder in single-atom and single-site catalysts on amorphous supports. The distribution of passage times (or waiting times) through a complex catalytic network originating from a set of coupled active sites is described by a probability distribution function (PDF), f(t), that reflects the local environment of the reaction center. An efficient algorithm is developed based on the linear algebra of the Markov transition matrix that produces f(t) or its moments. The kinetics of the hydrogenation reaction of styrene on an organovanadium(III) catalyst supported on amorphous silica is studied. A kinetic model consisting of three intertwined catalytic cycles emanating from three chemically distinct active sites is proposed to describe the chemistry. Density functional theory (DFT) calculations are employed to determine the free energy barriers of the reactions, which are used to construct the rate coefficient matrix. The disorder induced by the amorphous support material is divided into a low-dimensional short-range component reflecting the covalent structures near the reaction center and a weaker long-range component modeling the bulk randomness. The results are computed and analyzed for a wide range of concentration values and disorder scenarios. The unusual structure in the f(t) PDF is found to occur for certain cases that reveal the contribution of multiple catalytic pathways acting in concert.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-order harmonic generation in solid C 60

High-order harmonic generation (HHG) has unleashed the power of strong laser physics in solids. Here we investigate HHG from a large system, solid C 60 , with 240 valence electrons engaging harmonic generation at each crystal momentum. We employ the density functional theory and the time-dependent Liouville equation of the density matrix to compute HHG signals. We find that under a moderately strong laser pulse, HHG signals reach 15th order, consistent with the experimental results from C 60 plasma. The helicity dependence in solid C 60 is weak, due to the high symmetry. In contrast to the general belief, HHG is unsuitable for band structure mapping in C 60 . However, we find a window of opportunity using a long wavelength, where harmonics are generated through multiple-photon excitation. In particular, the fifth-order harmonic energies closely follow the transition energy dispersion between the valence and conduction bands. Finally, this finding is expected to motivate future experimental investigations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Laser-induced selective local patterning of vanadium oxide phases

The same elements can form different compounds with widely different physical properties. Synthesis of a single-phase material is commonly achieved by controlling experimental conditions. Synthesizing materials that incorporate multiple specific spatially distributed chemical phases is often challenging, especially if different phases must be organized into well-defined spatial patterns. Here, we present an efficient solid reaction laser annealing (SRLA) approach to directly write regions of different local chemical compositions. We demonstrate the practical utility of our approach by locally writing microscale patterns of distinct chemical phases in vanadium oxide thin films. Specifically, we achieved the controlled local recrystallization of a uniform V 2 O 3 matrix into VO 2 , V 3 O 5 , and V 4 O 7 regions exhibiting sharp 1st- and 2nd-order metal–insulator phase transitions over a wide range of critical temperatures, i.e., a characteristic feature of select vanadium oxides that is extremely sensitive to even minute structural or compositional imperfections. We utilized the local chemical phase writing to pattern spiking oscillators with distinct electrical behavior directly in the thin film sample without employing elaborate lithography fabrication. Our laser tuning local chemical composition opens a pathway to synthesize a wide range of artificially micropatterned composite materials, with precision and control unattainable in conventional material synthesis methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗