Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix models”

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 793 records · Page 44

Determining g A / g V with High-Resolution Spectral Measurements Using a LiInSe 2 Bolometer

Neutrinoless double beta decay (0νββ) processes sample a wide range of intermediate forbidden nuclear transitions, which may be impacted by quenching of the axial vector coupling constant (g A /g V ), the uncertainty of which plays a pivotal role in determining the sensitivity reach of 0νββ experiments. Here we present measurements performed on a high-resolution LiInSe 2 bolometer in a “source = detector” configuration to measure the spectral shape of the fourfold forbidden β decay of 115 In . The value of g A /g V is determined by comparing the spectral shape of theoretical predictions to the experimental β spectrum taking into account various simulated background components as well as a variety of detector effects. We find evidence of quenching of g A /g V at >5 σ with a model-dependent quenching factor of 0.655 ± 0.002 as compared to the free-nucleon value for the interacting shell model. We also measured the 115 In half-life to be [5.18 ± 0.06(stat) $^{+0.005}_{- 0.015}$(sys)] x 10 14 yr within the interacting shell model framework. This Letter demonstrates the power of the bolometeric technique to perform precision nuclear physics single-β decay measurements, which along with improved nuclear modeling can help reduce the uncertainties in the calculation of several decay nuclear matrix elements including those used in 0νββ sensitivity calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A Statistical Interpolation Code for Ocean Analysis and Forecasting

Abstract We present a data assimilation package for use with ocean circulation models in analysis, forecasting, and system evaluation applications. The basic functionality of the package is centered on a multivariate linear statistical estimation for a given predicted/background ocean state, observations, and error statistics. Novel features of the package include support for multiple covariance models, and the solution of the least squares normal equations either using the covariance matrix or its inverse—the information matrix. The main focus of this paper, however, is on the solution of the analysis equations using the information matrix, which offers several advantages for solving large problems efficiently. Details of the parameterization of the inverse covariance using Markov random fields are provided and its relationship to finite-difference discretizations of diffusion equations are pointed out. The package can assimilate a variety of observation types from both remote sensing and in situ platforms. The performance of the data assimilation methodology implemented in the package is demonstrated with a yearlong global ocean hindcast with a 1/4° ocean model. The code is implemented in modern Fortran, supports distributed memory, shared memory, multicore architectures, and uses climate and forecasts compliant Network Common Data Form for input/output. The package is freely available with an open source license from www.tendral.com/tsis/ .

Srinivasan, Ashwanth↗

Roadmap for thermal property measurements of Molten Salt Reactor systems

Thermal properties of molten salts, both with and without fissile and fertile material, are fundamental to modeling and simulation that supports design, construction, operation, and accident progression evaluation of molten salt reactors. This roadmap provides guidance to the research community on the types and amounts of required data for base, or carrier, salt systems, recommends simulated burnup fuel salt compositions to evaluate, and provides a technical foundation for a focused salt property evaluation matrix. The property values generated from both computation and experimental methods are used to develop models that constitute the Molten Salt Thermal Properties Database (MSTDB). While the building blocks of MSTDB are fundamental in nature, its data facilitates constructing engineering models that allow stakeholders to investigate the behavior of specific compositions of interest and supports broader modeling and simulation of MSRs through coupling to multi-physics, multi-scale mass accountancy tools.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Numerical characterization of support recovery in sparse regression with correlated design

Sparse regression is employed in diverse scientific settings as a feature selection method. A pervasive aspect of scientific data is the presence of correlations between predictive features. These correlations hamper both feature selection and estimation and jeopardize conclusions drawn from estimated models. On the other hand, theoretical results on sparsity-inducing regularized regression have largely addressed conditions for selection consistency via asymptotics, and disregard the problem of model selection, whereby regularization parameters are chosen. In this numerical study, we address these issues through exhaustive characterization of the performance of several regression estimators, coupled with a range of model selection strategies. These estimators and selection criteria were examined across correlated regression problems with varying degrees of signal to noise, distributions of non-zero model coefficients, and model sparsity. Our results reveal a fundamental tradeoff between false positive and false negative control in all regression estimators and model selection criteria examined. Additionally, we numerically explore a transition point modulated by the signal-to-noise ratio and spectral properties of the design covariance matrix at which the selection accuracy of all considered algorithms degrades. Overall, we find that SCAD coupled with BIC or empirical Bayes model selection performs the best feature selection across the regression problems considered.

97 MATHEMATICS AND COMPUTING↗

Robust design of semi-automated clustering models for 4D-STEM datasets

Materials discovery and design require characterizing material structures at the nanometer and sub-nanometer scale. Four-Dimensional Scanning Transmission Electron Microscopy (4D-STEM) resolves the crystal structure of materials, but many 4D-STEM data analysis pipelines are not suited for the identification of anomalous and unexpected structures. This work introduces improvements to the iterative Non-Negative Matrix Factorization (NMF) method by implementing consensus clustering for ensemble learning. We evaluate the performance of models during parameter tuning and find that consensus clustering improves performance in all cases and is able to recover specific grains missed by the best performing model in the ensemble. The methods introduced in this work can be applied broadly to materials characterization datasets to aid in the design of new materials.

Bruefach, Alexandra (ORCID:0000000209323477)↗

Distributed out-of-memory NMF on CPU/GPU architectures

We propose an efficient distributed out-of-memory implementation of the non-negative matrix factorization (NMF) algorithm for heterogeneous high-performance-computing systems. The proposed implementation is based on prior work on NMFk, which can perform automatic model selection and extract latent variables and patterns from data. In this work, we extend NMFk by adding support for dense and sparse matrix operation on multi-node, multi-GPU systems. The resulting algorithm is optimized for out-of-memory problems where the memory required to factorize a given matrix is greater than the available GPU memory. Memory complexity is reduced by batching/tiling strategies, and sparse and dense matrix operations are significantly accelerated with GPU cores (or tensor cores when available). Input/output latency associated with batch copies between host and device is hidden using CUDA streams to overlap data transfers and compute asynchronously, and latency associated with collective communications (both intra-node and inter-node) is reduced using optimized NVIDIA Collective Communication Library (NCCL) based communicators. Benchmark results show significant improvement, from 32X to 76x speedup, with the new implementation using GPUs over the CPU-based NMFk. Good weak scaling was demonstrated on up to 4096 multi-GPU cluster nodes with approximately 25,000 GPUs when decomposing a dense 340 Terabyte-size matrix and an 11 Exabyte-size sparse matrix of density 10 -6 .

97 MATHEMATICS AND COMPUTING↗

Dimer description of the SU(4) antiferromagnet on the triangular lattice

In systems with many local degrees of freedom, high-symmetry points in the phase diagram can provide an important starting point for the investigation of their properties throughout the phase diagram. In systems with both spin and orbital (or valley) degrees of freedom such a starting point gives rise to SU(4)-symmetric models. Here we consider SU(4)-symmetric "spin" models, corresponding to Mott phases at half-filling, i.e. the six-dimensional representation of SU(4). This may be relevant to twisted multilayer graphene. In particular, we study the SU(4) antiferromagnetic "Heisenberg" model on the triangular lattice, both in the classical limit and in the quantum regime. Carrying out a numerical study using the density matrix renormalization group (DMRG), we argue that the ground state is non-magnetic. We then derive a dimer expansion of the SU(4) spin model. An exact diagonalization (ED) study of the effective dimer model suggests that the ground state breaks translation invariance, forming a valence bond solid (VBS) with a 12-site unit cell. Finally, we consider the effect of SU(4)-symmetry breaking interactions due to Hund's coupling, and argue for a possible phase transition between a VBS and a magnetically ordered state.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Bootstrapping 2d ϕ 4 theory with Hamiltonian truncation data

We combine the methods of Hamiltonian Truncation and the recently proposed generalisation of the S-matrix bootstrap that includes local operators to determine the two-particle scattering amplitude and the two-particle form factor of the stress tensor at s > 0 in the 2d Φ 4 theory. We use the form factor of the stress tensor at s ≤ 0 and its spectral density computed using Lightcone Conformal Truncation (LCT), and inject them into the generalized S-matrix bootstrap set-up. The obtained results for the scattering amplitude and the form factor are fully reliable only in the elastic regime. We independently construct the “pure” S-matrix bootstrap bounds (bootstrap without including matrix elements of local operators), and find that the sinh-Gordon model and its analytic continuation the “staircase model” saturate these bounds. Surprisingly, the Φ 4 two-particle scattering amplitude also very nearly saturates these bounds, and moreover is extremely close to that of the sinh-Gordon/staircase model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Delamination-informed lifecycle decisions: A dielectric and machine learning framework for composite sorting and recycling

Composite materials are widely used in aerospace, marine, and automotive sectors due to their high strength-to-weight ratio and durability. However, their long-term reliability can be compromised by damage accumulation. Specifically, delamination initiation serves as a precursor to structural failure, which is often difficult to detect during damage inspection. Identifying and sorting delamination initiation in samples not only increases operational safety while providing critical information for end-of-life decisions, which influences both the service life extension value and the efficiency of fiber extraction during recycling. This research addresses two challenges: (1) developing a nondestructive, ex-situ framework to sort composite materials based on damage severity, particularly delamination, and (2) understanding how damage in composites influences resin removal during pyrolysis. Both experimental work and finite element analysis were performed to predict critical stress levels that are associated with delamination onset. Based on these results, three loading levels 50 %, 75 %, and 90 % of maximum stress, were selected for controlled experiments, generating composite samples with varying extents of damage for machine learning model training. Microscopic imaging of these samples confirmed the damage progression from matrix cracking to delamination, validating the computational predictions. We explored supervised machine learning using dielectric measurements to classify damage states. Preliminary results show an artificial neural network can identify early delamination which is a potential precursor to failure, with 94.44 % accuracy on our dataset. A parallel investigation into the effect of damage severity on pyrolysis recycling showed that heavily delaminated samples required significantly less energy for comparable matrix removal than undamaged samples.

dielectric variables↗

Effects of micrometer-scale surface roughness on thermal infrared emittance spectra of silica glass

Surface roughness is known to decrease thermal infrared (TIR) absorption band intensity, but studies of the effect on geologically relevant samples are relatively limited. To determine the effect of surface roughness (with features smaller than ~2/3 of the wavelength) on TIR spectra, we investigated two glass compositions with prepared roughened surfaces: (1) high purity fused silica and (2) soda-lime glass (73 wt% SiO$_2$). We roughened the surfaces of the glasses by sandblasting and polishing with grit paper. The surfaces were characterized with scanning electron microscopy and stylus profilometry. We then analyzed the roughened glasses with TIR emittance spectroscopy. Micrometer-scale roughness causes a decrease in TIR absorption band intensity, relative to a specular surface. No significant changes in band shape or shifts in wavelength were detected. As roughness increases, empirical results show a logarithmic decrease in TIR absorption band intensity. The logarithmic trends of the two glass compositions are different; empirical roughness calibrations do not translate across different compositions. A linear, least-squares spectral deconvolution using two endmembers, specular and blackbody, predicts model spectra of roughened glass surfaces with relatively low error. This is of consequence to orbital TIR measurements of poorly constrained targets, such as the martian surface, because micrometer-scale roughness is adequately modeled by the addition of a blackbody spectrum to the deconvolution endmember matrix.

36 MATERIALS SCIENCE↗

Machine-Learning-Guided Insights into Solid-Electrolyte Interphase Conductivity: Are Amorphous Lithium Fluorophosphates the Key?

Despite decades of study, the identity of the dominant Li + -conducting phase within the inorganic SEI of Li-ion batteries remains unresolved. While the mosaic model describes LiF/Li 2 O/Li 2 CO 3 nanocrystallites within a disordered matrix, these crystalline phases inherently offer limited ionic conductivity. Growing evidence suggests that interfaces, grain boundaries, and amorphous phases may instead host the primary fast-ion pathways. Using diffusion-based generative structure prediction and machine-learning interatomic potentials (MLIPs), we investigate lithium difluorophosphate (LiPO 2 F 2 ), a key mixed-anion decomposition product of phosphorus- and fluorine-containing electrolytes. We identify a stable crystalline polymorph and demonstrate that the amorphous counterpart is conductive, with projected room-temperature σ ≈ 0.18 mS cm –1 and E a ≈ 0.40 eV. Here, this enhancement stems from structural disorder flattening the Li site-energy landscape and a low formation energy for Li-interstitial defects, which supplies additional mobile carriers. We propose amorphous mixed-anion Li-P-O-F phases as a promising conducting medium in the SEI, offering a specific target for engineering improved battery interfaces.

Zhong, Peichen [University of California, Berkeley↗

Edge Reconstruction in a Quantum Spin Hall Insulator

We study interaction-driven edge reconstruction in a quantum spin Hall insulator described by the Bernevig-Hughes-Zhang model with Kanamori-Hubbard interactions using the real-space density matrix renormalization group method in both the grand-canonical and canonical ensembles. For a two-dimensional cylinder with a smooth edge, we identify discrete particle-number transitions that lead to a spin-polarized edge state stabilized by an emergent ferromagnetic exchange interaction. The reconstruction is orbital-selective, occurring predominantly in the 𝑠-orbital channel. Our results reveal a microscopic mechanism for emergent fluctuating moments at the edge that could compromise the topological protection of helical edge states by time-reversal symmetry.

Soni, Rahul [ORNL] (ORCID:0000000317714299)↗

Ultrafast laser-driven many-body dynamics and Kondo coherence collapse

Ultrafast laser pulse has provided a systematic way to inspect the dynamics of electrons in condensed matter systems. In this paper, by means of time-dependent density matrix renormalization group, we study an ultrafast laser-driven Kondo lattice model, in which conduction electrons are strongly coupled with magnetically local moments. The single-particle spectral function due to strong correlation effects and photon emission in the nonequilibrium states under laser driving are calculated. We find laser field excited collective doublon-hole pairs and an associated transient melting of Kondo coherence phase, signifying the collapse of Kondo energy gap. Moreover, we show that the photon emission, induced by a strong laser field, exhibits a different intensity characteristics than in the equilibrium Kondo insulator, which could be explained by the Kondo collapse and related suppression of both intraband and interband contribution in Kondo melting liquid. This theoretical insight is accessible with time- and angle-resolved photoemission spectroscopy and high-harmonic generation spectroscopy and will stimulate the investigation of nonequilibrium dynamics and nonlinear phenomenon in heavy fermion systems.

36 MATERIALS SCIENCE↗

Quantum order by disorder in frustrated spin nanotubes

Here we investigate quantum order by disorder in a frustrated spin nanotube formed by wrapping a J1–J2 Heisenberg model at 45° around a cylinder. Using Schwinger boson theory and density matrix renormalization group (DMRG), we have computed the ground-state phase diagram to reveal a Z2 phase in which collinear spin stripes form a right- or left-handed helix around the nanotube. We have derived an analytic estimate for the critical η c = J 1 /2J 2 of the Z 2 -helical phase transition, which is in agreement with the DMRG results. By evaluating the entanglement spectrum and nonlocal string order parameters we discuss the topology of the Z 2 -helical phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Python Codebase and Jupyter Notebooks - Applications of Machine Learning Techniques to Geothermal Play Fairway Analysis in the Great Basin Region, Nevada

Git archive containing Python modules and resources used to generate machine-learning models used in the "Applications of Machine Learning Techniques to Geothermal Play Fairway Analysis in the Great Basin Region, Nevada" project. This software is licensed as free to use, modify, and distribute with attribution. Full license details are included within the archive. See "documentation.zip" for setup instructions and file trees annotated with module descriptions.

Brown, Stephen↗

HERO WEC V1.0 2024 - WEC-Sim Detailed Simulation Runs and Summary Data

This dataset includes results from simulations of NREL's hydraulic and electric reverse osmosis wave energy converter (HEREO WEC). Simulation runs include 135 wave cases that were based on the updated WEC-Sim model, which is linked below. The data represented in this repository is based on an updated WEC-Sim model using laboratory data to tune and refine the original WEC-Sim model for the V1.0 HERO WEC. The 135 wave cases represent waves with the following wave height and wave period ranges: - Significant Wave Height: 0.25 - 3.75m in 0.25m increments - Wave Period: 5 - 13 sec in 1 sec increments Each run was simulated using a Pierson-Moskowitz irregular wave spectrum with a 100 second ramp time, a total simulation time of 3,100 seconds, and a simulation time-step of 0.005s. A reference table has been included to map each multi condition run (MCR) case with each wave condition. Summary data set includes a spreadsheet and image files with matrices that are associated with data from simulation runs. All matrices cover the same significant wave height and wave periods from the simulation runs, in the same increments. The following matrices are included: - Power Abs: The average absorbed power from the WEC (calculated from anchor reaction force and heave velocity) - Power Hyd: The average hydraulic power output at pump (calculated from pump output flow and pressure) - Power - Hyd ROi: The average hydraulic power measured at the RO system inlet (calculated from RO system pressure and flow (pre-accumulator)) - Flow - Pump out: The average flowrate measured at the pump outlet - Flow - Perm: The average permeate (clean water) production - Flow - RO (pre): The average flowrate measured at the inlet of the RO system before the accumulators - Flow - RO (post): The average flowrate measured after the accumulator bank in the RO system - Pressure - RO: The average pressure measured at the inlet of the RO system This data set has been developed by the National Renewable Energy Laboratory, operated by the Alliance for Sustainable Energy, LLC, for the U.S. Department of Energy (DOE) under Contract No. DE-AC36-08GO28308. Funding provided by the U.S. Department of Energy Office of Energy Efficiency and Renewable Energy Water Power Technologies Office.

16 TIDAL AND WAVE POWER↗

Encapsulant Void Breakdown Voltage Modeling & Experiments

For high voltage electrical devices, prevention of high voltage breakdown is critical for device function. Use of polymeric encapsulation such as epoxies is common, but these may include air bubbles or other voids of varying size. The present work aimed to model and experimentally determine the size dependence of breakdown voltage for voids in an epoxy matrix, as a step toward establishing size criteria for void screening. Effects were investigated experimentally for both one-dimensional metal/epoxy/air/epoxy/metal gap sizes from 50 μm to 10 mm, as well as spherical voids of 250 μm, 500 μm, 1 mm and 2 mm sizes. These experimental results were compared to modified Paschen curve and particle-in-cell discharge models; minimum breakdown voltages of 6 - 8.5 kV appeared to be predicted by 1D models and experiments, with minimum breakdown voltage for void sizes of 0.2 - 1 mm. In a limited set of 3D experiments on 250 μm, 500 μm, 1 mm and 2 mm voids within epoxy, the minimum breakdown voltages observed were 18.5 - 20 kV, for 500 μm void sizes. These experiments and models are aimed at providing initial size and voltage criteria for tolerable void sizes and expected discharge voltages to support design of encapsulated high voltage components.

42 ENGINEERING↗