Engineering Papers⌕ Search

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Tough and circular glass fiber composites via a tailored dynamic boronic ester interface

Glass fiber reinforced polymer (GFRP) composites are valued for their strength and cost-effectiveness. However, traditional GFRPs often face challenges for end-of-life recycling due to their non-depolymerizable thermoset matrices, and long-term performance due to inadequate interfacial adhesion, which can lead to fiber–matrix delamination. Here, in this work, we have designed dynamic fiber–matrix interfaces to allow tough and closed-loop recyclable GFRPs by utilizing a vitrimer, derived from upcycled polystyrene-b-poly(ethylene-co-butylene)-b-polystyrene (SEBS) with boronic ester (S-Bpin) and amine-based diol crosslinker. The boronic ester groups in S-Bpin form dynamic covalent bonds with the naturally present hydroxyl groups on the unsized GF surface, which eliminates the need for fiber sizing and enables facile closed-loop recyclability of both the fibers and the vitrimer matrix. The resulting strong fiber–matrix interface, depicted by the Raman mapping, leads to a 552% increase in-plane shear toughness (6.2 ± 0.3 MJ m -3 ) and 27% ultimate tensile strength (361 ± 89.2 MPa) compared to those of the conventional epoxy-based matrix (0.95 ± 0.05 MJ m -3 and 264 ± 59.7 MPa, respectively). The network rearrangement through dynamic boronic ester exchange enables fast thermoformability and repairability of micro-cracks at elevated temperatures. Additionally, both the matrix and composite demonstrate strong adhesion to various surfaces including steel and glasses exhibiting ≥6 MPa lap shear strength, which expands their suitability for diverse industrial applications. The readily created dynamic interface between boronic ester functionalized vitrimer and neat GFs presents a promising strategy for developing closed-loop recyclable, multifunctional structural materials, offering a sustainable alternative to non-recyclable thermoset GFRPs and contributes to a circular economy in composite materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The impact of nanoparticle softness on its tracer diffusion coefficient in all polymer nanocomposites

The diffusion of nanoparticles in a polymer matrix is an area of current interest. However, a complete understanding is still limited as it is often difficult to quantify the much slower motion of nanoparticles in a polymer matrix. To combat this problem, we have developed a protocol to measure the diffusion coefficient of soft nanoparticles in a linear polymer matrix. Recently developed synthetic control over soft nanoparticle structures combined with this protocol provides a pathway to separately elucidate the effects of the molecular weight and nanoparticle softness on its diffusive behavior. These results indicate that the nanoparticle softness and deformability dictate its motion. Increasing the cross-linking density of the nanoparticle for all molecular weights increases its hardness and suppresses its motion in the linear matrix. Additionally, the nanoparticle molecular weight dependence deviates from the exponential dependence for star polymers suggesting that these nanoparticles benefit from the cooperative motion of the matrix to open pathways for the nanoparticle. Finally, comparison of these experimentally determined values to the Stokes–Einstein theory demonstrates that the nanoparticles diffuse much slower than a hard sphere. This is interpreted to indicate that there exist additional interactions between the nanoparticle and polymer matrix that are not captured by Stokes–Einstein, including threading or entanglement of the linear chain with the nanoparticle.

Rostom, Sahar↗
DOE OSTI2025

Particle hit clustering and identification using point set transformers in liquid argon time projection chambers

Liquid argon time projection chambers are often used in neutrino physics and dark-matter searches because of their high spatial resolution. The images generated by these detectors are extremely sparse, as the energy values detected by most of the detector are equal to 0, meaning that despite their high resolution, most of the detector is unused in a particular interaction. Instead of representing all of the empty detections, the interaction is usually stored as a sparse matrix, a list of detection locations paired with their energy values. Traditional machine learning methods that have been applied to particle reconstruction such as convolutional neural networks (CNNs), however, cannot operate over data stored in this way and therefore must have the matrix fully instantiated as a dense matrix. Operating on dense matrices requires a lot of memory and computation time, in contrast to directly operating on the sparse matrix. We propose a machine learning model using a point set neural network that operates over a sparse matrix, greatly improving both processing speed and accuracy over methods that instantiate the dense matrix, as well as over other methods that operate over sparse matrices. Compared to competing state-of-the-art methods, our method improves classification performance by 14%, segmentation performance by more than 22%, while taking 80% less time and using 66% less memory. Compared to state-of-the-art CNN methods, our method improves classification performance by more than 86%, segmentation performance by more than 71%, while reducing runtime by 91% and reducing memory usage by 61%.

calibration and fitting methods↗
DOE OSTI2023

Using angular two-point correlations to self-calibrate the photometric redshift distributions of DECaLS DR9

Calibrating the redshift distributions of photometric galaxy samples is essential in weak lensing studies. The self-calibration method combines angular auto- and cross-correlations between galaxies in multiple photometric redshift (photo-z) bins to reconstruct the scattering rates matrix between redshift bins. In this paper, we test a recently proposed self-calibration algorithm using the DECaLS Data Release 9 and investigate to what extent the scattering rates are determined. We first mitigate the spurious angular correlations due to imaging systematics by a machine learning based method. We then improve the algorithm for χ 2 minimization and error estimation. Finally, we solve for the scattering matrices, carry out a series of consistency tests, and find reasonable agreements: (1) finer photo-z bins return a high-resolution scattering matrix, and it is broadly consistent with the low-resolution matrix from wider bins; (2) the scattering matrix from the Northern Galactic Cap is almost identical to that from Southern Galactic Cap; (3) the scattering matrices are in reasonable agreement with those constructed from the power spectrum and the weighted spectroscopic subsample. We also evaluate the impact of cosmic magnification. Although it changes little the diagonal elements of the scattering matrix, it affects the off-diagonals significantly. The scattering matrix also shows some dependence on scale cut of input correlations, which may be related to a known numerical degeneracy between certain scattering pairs. In conclusion, this work demonstrates the feasibility of the self-calibration method in real data and provides a practical alternative to calibrate the redshift distributions of photometric samples.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Quantum Joule expansion of one-dimensional systems

Here, we investigate the Joule expansion of nonintegrable quantum systems that contain bosons or spinless fermions in one-dimensional lattices. A barrier initially confines the particles to be in half of the system in a thermal state described by the canonical ensemble and is removed at time t = 0. We investigate the properties of the time-evolved density matrix, the diagonal ensemble density matrix, and the corresponding canonical ensemble density matrix with an effective temperature determined by the total energy conservation using exact diagonalization. The weights for the diagonal ensemble and the canonical ensemble match well for high initial temperatures that correspond to negative effective final temperatures after the expansion. At long times after the barrier is removed, the time-evolved Rényi entropy of subsystems bigger than half can equilibrate to the thermal entropy with exponentially small fluctuations. The time-evolved reduced density matrix at long times can be approximated by a thermal density matrix for small subsystems. Few-body observables, like the momentum distribution function, can be approximated by a thermal expectation of the canonical ensemble with strongly suppressed fluctuations. The negative effective temperatures for finite systems go to non-negative temperatures in the thermodynamic limit for bosons, but there is a true thermodynamic effect for fermions, which is confirmed by finite-temperature density matrix renormalization group calculations. We propose the Joule expansion as a way to dynamically create negative temperature states for fermion systems with repulsive interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Reconstructing thermal quantum quench dynamics from pure states

Simulating the nonequilibrium dynamics of thermal states is a fundamental problem across scales from high-energy to condensed-matter physics. Quantum computers may provide a way to solve this problem efficiently. Preparing a thermal state on a quantum computer is challenging, but there exist methods to circumvent this by computing a weighted sum of time-dependent matrix elements in a convenient basis. Further, while the number of basis states can be large, in this paper we show that it can be reduced by simulating only the largest density matrix elements by weight, capturing the density matrix to a specified precision. Leveraging Hamiltonian symmetries enables further reductions. This approach paves the way to more accurate thermal-state dynamics simulations on near-term quantum hardware.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Assessing the Feasibility of Bordered Block Diagonal Reordering in Power System Matrices using Fully Convolutional Network

In electromagnetic transient (EMT) simulations for power systems and inverter-based resources (IBRs), the arrangement of states within the system's linear equations, represented by matrix A in Ax=b, is critical. The state ordering in matrix A can highlight distinct characteristics of the system's graph, and identifying an optimal state ordering is crucial for efficient computation. The choice of state ordering, however, is dependent on the solver used, as each solver may perform optimally with different matrix patterns. With a wide array of matrix reordering algorithms available, selecting the most suitable one becomes challenging without insights into the matrix's ideal configuration. To address this, the paper proposes a fully convolutional network (FCN) to evaluate the reordering potential of the A matrix into a bordered block diagonal (BBD) pattern, which is commonly observed in power system and IBR modeling. The FCN's assessment aims to streamline the solver's operation, which in turn could substantially reduce the computational time required to find a solution.

Xia, Qianxue↗
DOE OSTI2023

Parallel Memory-Independent Communication Bounds for SYRK

In this paper, we focus on the parallel communication cost of multiplying a matrix with its transpose, known as a symmetric rank-k update (SYRK). SYRK requires half the computation of general matrix multiplication because of the symmetry of the output matrix. Recent work (Beaumont et al., SPAA '22) has demonstrated that the sequential I/O complexity of SYRK is also a constant factor smaller than that of general matrix multiplication. Inspired by this progress, we establish memory-independent parallel communication lower bounds for SYRK with smaller constants than general matrix multiplication, and we show that these constants are tight by presenting communication-optimal algorithms. The crux of the lower bound proof relies on extending a key geometric inequality to symmetric computations and analytically solving a constrained nonlinear optimization problem. Here, the optimal algorithms use a triangular blocking scheme for parallel distribution of the symmetric output matrix and corresponding computation.

Communication costs↗
DOE OSTI2024

Flax fiber-reinforced fatty acid vitrimer biocomposite with enhanced chemical recyclability

Here, this study unveils a sustainable, easily recyclable biocomposite, leveraging the dynamic nature of covalently adaptive bonds in a vitrimer matrix. The fabrication involved a fatty acid-derived vitrimer as the polymer matrix and multi-layered, nonwoven flax mat as reinforcing scaffold. The incorporation of these fibers significantly improved the mechanical performance of the vitrimer matrix uniformly. The ester-based covalently adaptive network plays a crucial role in enabling exceptional fiber-matrix bonding, as well as recyclability. The vitrimer matrix dissolves in ethylene glycol through transesterification, facilitating complete material recovery and biocomposite recycling without compromising the original properties of the matrix and reinforcing fibers.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Wettability Alteration with Weak Acid-Assisted Surfactant Flood

Oil-wetness and heterogeneity are two key reasons for low oil recovery by waterflooding in carbonate reservoirs. Surfactants have been effective in altering the oil-wet matrix to a more water-wet condition and initiating spontaneous imbibition. Because it takes time for the surfactant to alter wettability, oil recovery from the tight matrix is slow and sometimes not economically feasible. Acids have the potential of dissolving minerals, which may alter wettability. In this study, the enhanced oil recovery (EOR) performance of an acid-assisted surfactant solution, a novel technique, was evaluated for low-temperature applications. A set of acids and their acetates were tested. Bulk rock-acid reaction, wettability alteration (WA) tests, and spontaneous imbibition measurements were conducted at reservoir temperature (35°C) to identify effective candidates. Coreflood tests were then performed to evaluate the selected acid-surfactant formulations. Before and after a coreflood test, the core was scanned using micro-computed tomography (CT) to investigate pore structure alteration. Bulk reaction measurements showed that weak acids, especially acetic acid (AA), have the desired low reaction rates at 35°C. WA tests showed that AA can remove the crude oil off the rock surface and alter wettability through mineral dissolution. The surfactant can reduce contact angles from 160° to 58°; adding acid into the surfactant can further reduce it to 52°. Spontaneous imbibition experiments showed the synergy between the acid and the surfactant; the AA-surfactant solution had the highest oil recovery (62.6%) among acid-surfactant formulations. The acid improves the WA efficiency by the surfactant through surface mineral dissolution and lower ζ-potential. The imbibition transports the acid-surfactant solutions into the matrix, which minimizes face dissolution. Coreflood tests show that the AA-surfactant flood can increase the oil recovery rate and recover about 8% more oil compared to the surfactant flood. Micro-CT showed that a few mineral particles were transported along the core and partially plugged pore throats, which reduced permeability and diverted flow leading to improved oil recovery. Here, the transport of the acid in reservoirs and the potential plugging issues have to be carefully evaluated in future studies.

58 GEOSCIENCES↗
DOE OSTI2025

Designing Physicochemically‐Ordered Interphases for High‐Performance Composites

To enhance the mechanical properties of carbon fiber‐reinforced polymer composites, a physicochemical scaffold is designed incorporating microscopically architected chemically reactive nanofibers that act as a multiscale bridge between the carbon fibers and the matrix. Thermally activated nanofibers leverage their morphologically driven mechanochemical properties to form covalent bonds with adjacent polymer molecules, creating a co‐continuous network that dramatically enhances fiber‐matrix load transfer. By meticulously controlling the nanofiber architecture through variable surface area, functional group availability, and polymer chain alignment effects, the extent of covalent bonding between nanofibers and the matrix is manipulated ultimately resulting in improved carbon fiber‐matrix adhesion. Further, the concept was validated using polyacrylonitrile nanofibers within an acrylonitrile butadiene styrene matrix in a discontinuous carbon fiber‐reinforced composite system. Nanomechanical studies using atomic force microscopy and low‐field nuclear magnetic resonance spectroscopy confirmed immobilized, chemically transferred, and ordered nanostructures at the interphase. The resulting composites demonstrated ≈56% and ≈175% improvements in tensile strength and toughness, respectively, compared to composites without nanofiber. Comprehensive thermal, rheological, and X‐ray scattering analyses, along side all‐atomic molecular dynamics simulations, revealed the fundamental mechanisms behind these improvements in mechanical behavior. The versatility and efficacy of the approach have the potential to address longstanding interphase challenges in the composite industry.

36 MATERIALS SCIENCE↗
DOE OSTI2020

ELSI — An open infrastructure for electronic structure solvers

Routine applications of electronic structure theory to molecules and periodic systems need to compute the electron density from given Hamiltonian and, in case of non-orthogonal basis sets, overlap matrices. System sizes can range from few to thousands or, in some examples, millions of atoms. Different discretization schemes (basis sets) and different system geometries (finite non-periodic vs. infinite periodic boundary conditions) yield matrices with different structures. The ELectronic Structure Infrastructure (ELSI) project provides an open-source software interface to facilitate the implementation and optimal use of high-performance solver libraries covering cubic scaling eigensolvers, linear scaling density-matrix-based algorithms, and other reduced scaling methods in between. In this paper, we present recent improvements and developments inside ELSI, mainly covering (1) new solvers connected to the interface, (2) matrix layout and communication adapted for parallel calculations of periodic and/or spin-polarized systems, (3) routines for density matrix extrapolation in geometry optimization and molecular dynamics calculations, and (4) general utilities such as parallel matrix I/O and JSON output. The ELSI interface has been integrated into four electronic structure code projects (DFTB+, DGDFT, FHI-aims, SIESTA), allowing us to rigorously benchmark the performance of the solvers on an equal footing. Based on results of a systematic set of large-scale benchmarks performed with Kohn–Sham density-functional theory and density-functional tight-binding theory, we identify factors that strongly affect the efficiency of the solvers, and propose a decision layer that assists with the solver selection process. As a result, we describe a reverse communication interface encoding matrix-free iterative solver strategies that are amenable, e.g., for use with planewave basis sets.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Predicting ductility of Mg/SiCp nanocomposite under multiaxial loading conditions based on unit cell modeling

This article presents an investigation on the ductility of Mg/SiCp under nine different multiaxial loading conditions based on a unit cell (or called representative volumetric element, RVE) finite element model. The constitutive relationship of magnesium matrix is characterized by a model combined Lode angle dependent plasticity model with the CPB06 plasticity model, which is able to describe the tension-compression asymmetry in magnesium matrix. The debonding of the matrix-particle interface and matrix fracture are considered by the modified Mohr-Coulomb (MMC) fracture model. Simulation results show that the interface debonding and subsequent matrix fracture are responsible for the ductility and fracture of nanocomposite. The computational results with parametric studies are compared and well correlated with experimental data, and it finally gives a stress-based full fracture locus of Mg/SiCp nanocomposite.

42 ENGINEERING↗
DOE OSTI2023

A fast, dense Chebyshev solver for electronic structure on GPUs

Matrix diagonalization is almost always involved in computing the density matrix needed in quantum chemistry calculations. In the case of modest matrix sizes (≲4000), performance of traditional dense diagonalization algorithms on modern GPUs is underwhelming compared to the peak performance of these devices. This motivates the exploration of alternative algorithms better suited to these types of architectures. We newly derive, and present in detail, an existing Chebyshev expansion algorithm whose number of required matrix multiplications scales with the square root of the number of terms in the expansion. Focusing on dense matrices of modest size, our implementation on GPUs results in large speed ups when compared to diagonalization. Additionally, we improve upon this existing method by capitalizing on the inherent task parallelism and concurrency in the algorithm. Furthermore, this improvement is implemented on GPUs by using CUDA and HIP streams via the MAGMA library and leads to a significant speed up over the serial-only approach for smaller (≲1000) matrix sizes. Finally, we apply our technique to a model system with a high density of states around the Fermi level, which typically presents significant challenges.

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

Online Voltage Event Detection Using Synchrophasor Data with Structured Sparsity-Inducing Norms

This paper develops an accurate and computationally efficient data-driven framework to detect voltage events from PMU data streams. It develops an innovative Proximal Bilateral Random Projection (PBRP) algorithm to quickly decompose the PMU data matrix into a low-rank matrix, a row-sparse event-pattern matrix and a noise matrix. Here, the row-sparse pattern matrix significantly distinguishes events from normal behavior. These matrices are then fed into a clustering algorithm to separate voltage events from normal operating conditions. Large-scale numerical study results on real-world PMU data show that the proposed algorithm is computationally more efficient and achieves higher F scores than state-of-the-art benchmarks.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Microstructure of Neutron-Irradiated Al 3 Hf-Al Thermal Neutron Absorber Materials

A thermal neutron-absorbing metal matrix composite (MMC) comprised of Al 3 Hf particles in an aluminum matrix was developed to filter out thermal neutrons and create a fast flux environment for material testing in a mixed-spectrum nuclear reactor. Intermetallic Al 3 Hf particles capture thermal neutrons and are embedded in a highly conductive aluminum matrix that provides conductive cooling of the heat generated due to thermal neutron capture by the hafnium. These Al 3 Hf-Al MMCs were fabricated using powder metallurgy via hot pressing. The specimens were neutron-irradiated to between 1.12 and 5.38 dpa and temperatures ranging from 286 °C to 400 °C. The post-irradiation examination included microstructure characterization using transmission electron microscopy (TEM) and energy-dispersive X-ray spectroscopy. This study reports the microstructural observations of four irradiated samples and one unirradiated control sample. All the samples showed the presence of oxide at the particle–matrix interface. The irradiated specimens revealed needle-like structures that extended from the surface of the Al 3 Hf particles into the Al matrix. An automated segmentation tool was implemented based on a YOLO11 computer vision-based approach to identify dislocation lines and loops in TEM images of the irradiated Al-Al 3 Hf MMCs. This work provides insight into the microstructural stability of Al 3 Hf-Al MMCs under irradiation, supporting their consideration as a novel neutron absorber that enables advanced spectral tailoring.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Normalization of ZZ instanton amplitudes in type 0B minimal superstring theory

Abstract We study ZZ instanton corrections in the (2,4k)$$ \mathcal{N} $$ N = 1 minimal superstring theory with the type 0B GSO projection, which becomes the type 0B$$ \mathcal{N} $$ N = 1 super-JT gravity in thek→ ∞ limit. Each member of the (2,4k) family of theories has two phases distinguished by the sign of the Liouville bulk cosmological constant. The worldsheet method for computing the one-loop normalization constant multiplying the instanton corrections gives an ill-defined answer in both phases. We fix these divergences using insights from string field theory and find finite, unambiguous results. Each member of the (2,4k) family of theories is dual to a double-scaled one-matrix integral, where the double-scaling limit can be obtained starting either from a unitary matrix integral with a leading one-cut saddle point, or from a hermitian matrix integral with a leading two-cut saddle point. The matrix integral exhibits a gap-closing transition, which is the same as the double-scaled Gross-Witten-Wadia transition whenk= 1. We also compute instanton corrections in the double-scaled matrix integral for allkand in both phases, and find perfect agreement with the string theory results.

Physics↗
DOE OSTI2024

Distributed strain sensing using Bi-metallic coated fiber optic sensors embedded in stainless steel

Silica fiber optic sensors are resistant to corrosive environments and high temperatures, making them attractive candidates for harsh conditions found in nuclear and aerospace industries. Moreover, fibers can be deployed remotely for continuous measuring of spatially distributed temperatures and strains. Here, this study investigated embedding a Ni/Cu bi-metallic coated fiber in a stainless-steel 316 (SS316) matrix using laser powder bed fusion towards functionalizing metal components for site-specific health monitoring. The embedded fiber was continuously interrogated during controlled heating to 1000°C. The measured fiber strains were similar to the expected differential thermal strains between the fiber and the SS316 matrix, until divergent behavior was observed at temperatures >500°C. No debonding at the matrix–coating–fiber interfaces was observed during microscopy, but significant interactions between the coatings and matrix resulted in diffusion-driven chemistry variations and Kirkendall void formation. Applying the strain-lag theory revealed plastic behavior in the Ni coating at temperatures >500°C, limiting the strain transfer to the fiber at higher temperatures. It was estimated that the elastic modulus in the Ni coating had decreased from ~200 GPa at room temperature to below 40 GPa, starting at 600°C. The low elastic modulus above 600°C is within the margin of what the tangent modulus would be in the case of bilinear isotropic hardening. Regardless of the divergent strain transfer at higher temperatures, the fiber was exposed to the equivalent of 1.9 % engineering strain at 1000°C, but measured only a 0.7 % engineering strain due to the poor strain transfer. Although compensating for the plastic behavior of Ni proved challenging, the bonding of a brittle silica fiber to a metal matrix surviving to 1000°C invites potential iterations on coating material for future application. For example, the embedded fiber is sufficient for acoustic energy transfer, realizing high temperature distributed acoustic sensing.

36 MATERIALS SCIENCE↗