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 1,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Pattern-matching indexing of Laue and monochromatic serial crystallography data for applications in materials science

Serial crystallography data can be challenging to index, as each frame is processed individually, rather than being processed as a whole like in conventional X-ray single-crystal crystallography. An algorithm has been developed to index still diffraction patterns arising from small-unit-cell samples. The algorithm is based on the matching of reciprocal-lattice vector pairs, as developed for Laue microdiffraction data indexing, combined with three-dimensional pattern matching using a nearest-neighbors approach. As a result, large-bandpass data (e.g. 5–24 keV energy range) and monochromatic data can be processed, the main requirement being prior knowledge of the unit cell. Furthermore, angles calculated in the vicinity of a few theoretical and experimental reciprocal-lattice vectors are compared, and only vectors with the highest number of common angles are selected as candidates to obtain the orientation matrix. Global matching on the entire pattern is then checked. Four indexing options are available, two for the ranking of the theoretical reciprocal-lattice vectors and two for reducing the number of possible candidates. Furthermore, the algorithm has been used to index several data sets collected under different experimental conditions on a series of model samples. Additionally, knowing the crystallographic structure of the sample and using this information to rank the theoretical reflections based on the structure factors helps the indexing of large-bandpass data for the largest-unit-cell samples. For small-bandpass data, shortening the candidate list to determine the orientation matrix should be based on matching pairs of reciprocal-lattice vectors instead of triplet matching.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Benchmark neutrinoless double- β decay matrix elements in a light nucleus

We compute nuclear matrix elements of neutrinoless double-β decay mediated by light Majorana-neutrino exchange in the $\textit{A}$ = 6 system. The goal for this work is to benchmark two many-body approaches, the no-core shell model and the multireference in-medium similarity renormalization group (SRG). We use a SRG-evolved next-to-next-to-next-to-leading chiral order two-body potential for the nuclear interaction, and make the approximation that the isospin is conserved. We compare the results of the two approaches as a function of the cutoff on the many-body basis space. Although differences are seen in the predicted nuclear radii, the ground-state energies and neutrinoless double-β decay matrix elements produced by the two approaches show significant agreement. We discuss the implications for calculations in heavier nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

A Linear-Complexity Tensor Butterfly Algorithm for Compressing High-Dimensional Oscillatory Integral Operators

This paper presents a multilevel tensor compression algorithm called tensor butterfly algorithm for efficiently representing large-scale and high-dimensional oscillatory integral operators, including Green's functions for wave equations and integral transforms such as Radon transforms and Fourier transforms. The proposed algorithm leverages a tensor extension of the so-called complementary low-rank property of existing matrix butterfly algorithms. The algorithm partitions the discretized integral operator tensor into subtensors of multiple levels and factorizes each subtensor at the middle level as a Tucker-type interpolative decomposition, whose factor matrices are formed in a multilevel fashion. For a d-dimensional (d > 1) integral operator discretized into a 2d-mode tensor with n2d entries, the overall CPU time and memory requirement scale as O(nd), in stark contrast to the O(nd log n) complexity of existing matrix algorithms such as matrix butterfly algorithms and fast Fourier transforms (FFTs), where n is the number of points per direction. When comparing with other tensor algorithms such as quantized tensor train (QTT), the proposed algorithm also shows superior CPU and memory performance for tensor contraction. Remarkably, the tensor butterfly algorithm can efficiently model high-frequency Green's function interactions between two unit cubes, each spanning 512 wavelengths per direction, which represents problems of scale over 512× larger than that existing butterfly algorithms can handle, with the same amount of computation resources. On the other hand, for a problem representing 64 wavelengths per direction, which is the largest size existing algebraic matrix algorithms can handle, our tensor butterfly algorithm exhibits 200x speedups and 30× memory reduction compared with existing ones. Moreover, the tensor butterfly algorithm also permits O(nd)-complexity FFTs and Radon transforms up to d = 6 dimensions.

Kielstra, P Michael↗
DOE OSTI2023

Simulation toolkit for digital material characterization of large image-based microstructures

In this paper, an efficient image-based simulation toolkit for material characterization is presented, which is scalable to work from personal computers to workstations. The effective thermal conductivity, elasticity, and permeability are evaluated employing a computational homogenization framework based on the Finite Element Method (FEM). Two complementary open-source packages are presented: one developed in Python, which can convert digital images into voxel meshes (pyTomoviewer); the other developed in Julia, that can run numerical simulations to compute effective material properties (chpack). Also, a CUDA C version of chpack is provided (chfem_gpu). They were designed to deal with large multi-phase models, so strategies were devised to minimize their memory footprint, while avoiding a high toll on execution time. The voxel-based approach significantly simplifies the FEM meshes and allows efficient matrix-free implementations. In that sense, to handle large linear systems of equations, the element-by-element (EBE) technique is adopted, in conjunction with a low-memory implementation of the Preconditioned Conjugate Gradient (PCG) method. Finally, the code was thoroughly tested on an artificial geometry made of a square array of cylinders, for which analytical solutions exist, as well as on a real micro-tomographic reconstruction of FiberForm TM , a carbon preform commonly used in thermal protection systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Matrix-Free High-Performance Saddle-Point Solvers for High-Order Problems in \(\boldsymbol{H}(\operatorname{\textbf{div}})\)

Here, this work describes the development of matrix-free GPU-accelerated solvers for high-order finite element problems in H(div). The solvers are applicable to grad-div and Darcy problems in saddle-point formulation, and have applications in radiation diffusion and porous media flow problems, among others. Using the interpolation–histopolation basis, efficient matrix-free preconditioners can be constructed for the (1, 1)-block and Schur complement of the block system. With these approximations, block-preconditioned MINRES converges in a number of iterations that is independent of the mesh size and polynomial degree. The approximate Schur complement takes the form of an M-matrix graph Laplacian and therefore can be well-preconditioned by highly scalable algebraic multigrid methods. High-performance GPU-accelerated algorithms for all components of the solution algorithm are developed, discussed, and benchmarked. Numerical results are presented on a number of challenging test cases, including the “crooked pipe” grad-div problem, the SPE10 reservoir modeling benchmark problem, and a nonlinear radiation diffusion test case.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Chirality transfer from chiral perovskite to molecular dopants via charge transfer states

Chiral perovskites are semiconductors with broken mirror symmetries. Their photo responses are often constrained in the UV range. In this work, we demonstrate that doping 2,3,5,6-Tetrafluoro-7,7,8,8-tetracyanoquinodimethane in the chiral perovskite matrix introduces a visible light absorption feature through the emerging charge-transfer electronic states. These charge-transfer states exhibits circular dichroism inherited from the chiral host, indicating effective chirality transfer from host to guest component via electronic coupling. Quantum-chemical modeling identifies a strong wave function overlap between an electron and a hole of the guest-host in a closely packed crystal configuration promoting the charge transfer state’s optical activity. We further integrate the doped chiral perovskite film into photodetectors and demonstrate a selective detection of circularly polarized light in both UV and visible regions. Our results suggest a universal approach of introducing visible photo absorption states to the chiral matrix to broaden the optical active range while enhancing the electrical conductivity.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Particle‐Size‐Dependent Lithium‐Ion Transport in PEO/LLZO Composite Electrolytes

Lithium-metal batteries with solid electrolytes can deliver higher energy density and improved safety than conventional Li-ion batteries. Among solid electrolyte candidates, polymer/ceramic composite electrolytes are attractive because they combine polymer flexibility with the high ionic conductivity of ceramics. However, whether ceramic fillers synergistically reduce polarization losses in the polymer matrix remains unclear. A central unknown is the critical polymer/ceramic interfacial resistance (Rint,crit), below which adding ceramics lowers electrolyte overpotential. Here, we present the first macroscale model framework to quantify R int,crit for composite electrolytes based on polyethylene oxide (PEO) and Ta-doped Li 7 La 3 Zr 2 O 12 (LLZO). A 1D model for DC-polarization of tri-layer cells (PEO-LiTFSI/LLZO/PEO-LiTFSI) shows that LLZO surface functionalization reduces the PEO/LLZO interfacial resistance, consistent with electrochemical impedance measurements. Extending to a 2D composite model, we show notably that Rint,crit scales linearly with LLZO particle diameter and shifts toward experimentally accessible values (e.g., 28.8 Ωcm 2 ) as particle size increases. At fixed ceramic volume fraction, larger LLZO particles reduce the number of interfacial crossings, driving more current through the ceramic phase and lowering concentration polarization. In contrast, R int,crit is largely independent of ceramic volume fraction. These results demonstrate that ceramic filler-size engineering can enable synergistic, energy-efficient transport in polymer/ceramic composite electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2026

Mitigation of polysulfide shuttle effect in Li-S batteries through catalytic disproportionation reaction

Polysulfides are poorly retained within porous cathodes and readily diffuse into the electrolyte over time, leading to the well-known shuttle effect that undermines the reversibility of Li-S batteries. Here, in this study, we demonstrate that catalytic disproportionation of polysulfides provides an effective pathway to suppress this process by rapidly converting dissolved species into solid sulfur and sulfides, thereby preventing their migration into the electrolyte. Fundamentally, the sluggish kinetics of sulfur redox reactions are responsible for the accumulation and redistribution of soluble polysulfides in the bulk electrolyte. By accelerating these kinetics, catalyzed disproportionation not only confines sulfur within the conductive cathode matrix but also promotes the homogeneous precipitation of Li₂S₂/Li₂S, which enhances electrochemical reversibility and cycling stability. Using nitrogen-doped carbon (NC800) as a model catalyst, we reveal its ability to drive a pseudo-16-electron reduction pathway, leading to a single dominant Li₂S product and uniform deposition within the porous framework. In contrast, a non-catalytic carbon (KB) yields multiple polysulfide intermediates and heterogeneous deposition. The mechanistic insights provided here highlight the pivotal role of catalytic disproportionation in reshaping sulfur redox pathways and offer a rational strategy for mitigating polysulfide shuttling in practical Li-S pouch cells.

25 ENERGY STORAGE↗
DOE OSTI2020

Modeling The Effects of Loading Scenario and Thermal Expansion Coefficient on Potential Failure of Cryo-compressed Hydrogen Vessels

A multiscale thermomechanical model for a simplified Type-3 cryogenic compressed hydrogen (H2) storage vessel is developed in this paper. The model accounts for the temperature-dependent elastic-plastic behavior of the vessel carbon/epoxy composite overwrap and aluminum alloy liner. The homogenized thermo-elastic-plastic behavior for the individual laminas of the vessel layup is obtained by an incremental Eshelby-Mori-Tanka approach associated with a micromechanical failure criterion to predict lamina failure while a standard elastic-plastic constitutive model is used to describe the behavior of a typical aluminum alloy assumed for the liner. The vessel response to external loadings is achieved by a finite element method. Four loading scenarios representing four thermomechanical cycles applied to the vessel are analyzed to evaluate constituent and lamina stresses as well as the associate failure criterion during the cycle according to these scenarios. The model can provide helpful guidance to mitigate thermal stresses by an adequate selection of loading scenario, optimizing the layup and by tailoring thermomechanical properties of the resin matrix.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Random insights into the complexity of two-dimensional tensor network calculations

Projected entangled pair states (PEPS) offer memory-efficient representations of some quantum many-body states that obey an entanglement area law and are the basis for classical simulations of ground states in two-dimensional (2d) condensed matter systems. However, rigorous results show that exactly computing observables from a 2d PEPS state is generically a computationally hard problem. Yet approximation schemes for computing properties of 2d PEPS are regularly used, and empirically seen to succeed, for a large subclass of (“not too entangled”) condensed matter ground states. Adopting the philosophy of random matrix theory, in this work, we analyze the complexity of approximately contracting a 2d random PEPS by exploiting an analytic mapping to an effective replicated statistical mechanics model that permits a controlled analysis at a large bond dimension. Through this statistical-mechanics lens, we argue that (i) although approximately sampling wave-function amplitudes of random PEPS faces a computational-complexity phase transition above a critical bond dimension, and (ii) one can generically efficiently estimate the norm and correlation functions for any finite bond dimension. Furthermore, these results are supported numerically for various bond-dimension regimes. It is an important open question whether the above results for random PEPS apply more generally also to PEPS representing physically relevant ground states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

New method for fitting coefficients in standard model effective theory

We present an alternative method for carrying out a principal-component analysis of Wilson coefficients in standard model effective field theory (SMEFT). The method is based on singular-value decomposition (SVD). The SVD method provides information about the sensitivity of experimental observables to physics beyond the standard model that is not accessible in the Fisher-information method. In principle, the SVD method can also have computational advantages over diagonalization of the Fisher information matrix. We demonstrate the SVD method by applying it to the dimension-6 coefficients for the process of top-quark decay to a b quark and a W boson and use this example to illustrate some pitfalls in widely used fitting procedures. We also outline an iterative procedure for applying the SVD method to dimension-8 SMEFT coefficients.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

𝐵 → 𝜌⁢ℓ⁢$\bar{v}$ Resonance Form Factors from 𝐵→ 𝜋⁢𝜋⁢ℓ⁢$\bar{v}$ in Lattice QCD

The decay 𝐵 → 𝜌⁢ℓ⁢$\bar{v}$ is an attractive process for determining the magnitude of the smallest Cabibbo-Kobayashi-Maskawa matrix element, |𝑉 𝑢⁢𝑏 |, and can provide new insights into the origin of the long-standing exclusive-inclusive discrepancy in determinations of this standard-model parameter. This requires a nonperturbative QCD calculation of the 𝐵 → 𝜌 form factors 𝑉, 𝐴 0 , 𝐴 1 , and 𝐴 12 . The unstable nature of the 𝜌 resonance has prevented precise lattice QCD calculations of these form factors to date. Here, we present the first lattice QCD calculation of the 𝐵 → 𝜌 form factors in which the 𝜌 is treated properly as a resonance in 𝑃-wave 𝜋⁢𝜋 scattering. To this end, we use the Lellouch-Lüscher finite-volume formalism to compute the 𝐵 → 𝜋⁢𝜋 form factors as a function of both momentum transfer and 𝜋⁢𝜋 invariant mass, and then analytically continue to the 𝜌 resonance pole. This calculation is performed with 2 + 1 dynamical quark flavors at a pion mass of approximately 320 MeV, and demonstrates a clear path toward results at the physical point.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

A hybrid machine-learning approach for analysis of methane hydrate formation dynamics in porous media with synchrotron CT imaging

Fast multi-phase processes in methane hydrate bearing samples pose a challenge for quantitative micro-computed tomography study and experiment steering due to complex tomographic data analysis involving time-consuming segmentation procedures. This is because of the sample's multi-scale structure, which changes over time, low contrast between solid and fluid materials, and the large amount of data acquired during dynamic processes. Here, a hybrid approach is proposed for the automatic segmentation of tomographic data from time-resolved imaging of methane gas-hydrate formation in sandy granular media, which includes a deep-learning 3D U-Net model. To prepare a training dataset for the 3D U-Net, a technique to automate data labeling based on sample-specific information about the mineral matrix immobility and occasional fluid movement in pores is proposed. Automatic segmentation allowed for studying properties of the hydrate growth in pores, as well as dynamic processes such as incremental flow and redistribution of pore brine. Results of the quantitative analysis showed that for typical gas-hydrate stability parameters (100 bar methane pressure, 7°C temperature) the rate of formation is slow (less than 1% per hour), after which the surface area of contact between brine and gas increases, resulting in faster formation (2.5% per hour). Hydrate growth reaches the saturation point after 11 h of the experiment. Finally, the efficacy of the proposed segmentation scheme in on-the-fly automatic data analysis and experiment steering with zooming to regions of interest is demonstrated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Controlling cantilevered adaptive X-ray mirrors

Modeling the behavior of a prototype cantilevered X-ray adaptive mirror (held from one end) demonstrates its potential for use on high-performance X-ray beamlines. Similar adaptive mirrors are used on X-ray beamlines to compensate optical aberrations, control wavefronts and tune mirror focal distances at will. Controlled by 1D arrays of piezoceramic actuators, these glancing-incidence mirrors can provide nanometre-scale surface shape adjustment capabilities. However, significant engineering challenges remain for mounting them with low distortion and low environmental sensitivity. Finite-element analysis is used to predict the micron-scale full actuation surface shape from each channel and then linear modeling is applied to investigate the mirrors' ability to reach target profiles. Using either uniform or arbitrary spatial weighting, actuator voltages are optimized using a Moore–Penrose matrix inverse, or pseudoinverse, revealing a spatial dependence on the shape fitting with increasing fidelity farther from the mount.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Exploration with Scalable Gaussian Process Reinforcement Learning

Exploration is a challenging problem in reinforcement learning (RL), especially in environments with sparse rewards. Quantifying and utilizing the parametric uncertainty has been shown to be paramount for successful exploration [Osband et al., 2018]. Bayesian, or approximately Bayesian, methods present a principled means of estimating the parametric uncertainty in RL problems. Gaussian processes, nonparametric Bayesian models, are often impractical due to poor scalability and computational bottlenecks. We introduce a scalable Gaussian process RL (GPRL) method which directly induces sparsity in the covariance matrix to facilitate faster computation. This is a departure from previous GPRL methods which instead rely on data reduction and subsampling. We compare various covariance-based exploration techniques (Thompson sampling, upper confidence bound, and probabilistic maximum variance) which leverage our scalable GP framework in sparse reward environments. Finally, we show favorable comparison against the bootstrapped deep Q-Network.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Development and demonstration of a BISON–Griffin modeling framework for the design of targeted TRISO transient experiments in the Transient Reactor Test Facility

Uranium oxycarbide (UCO)-bearing tri-structural isotropic (TRISO) particle fuels are expected to be used in numerous U.S. commercial reactor applications within the next decade. Here, in this work, we reviewed historical particle fuel transient experiments to identify gaps in TRISO fuel performance transient testing. A BISON–Griffin modeling framework was then developed to conduct preliminary TRISO transient analyses and begin to address these gaps. The framework was demonstrated using limiting-case transient conditions from a prototypic high-temperature gas-cooled reactor (HTGR). It was then applied to develop a matrix of experiments that could be performed in the Transient Reactor Test Facility (TREAT) to (1) evaluate UCO-fueled particle performance at moderate and high heat rates, (2) assess whether historical testing involving UO 2 -fueled particles is applicable to modern UCO-fueled particles, (3) deconvolute the impacts of temperature and heat rate on particle transient response, and (4) collect the data needed for fuel performance model validation and/or further development.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Phase-field simulations of fission gas bubble growth and interconnection in U-(Pu)-Zr nuclear fuel

Abstract The growth and interconnection of fission gas bubbles in the hotter central regions of U-(Pu)-Zr nuclear fuel has been simulated with a phase-field model. The Cahn-Hilliard equation was used to represent the two-phase microstructure, with a single defect species. The volume fraction of the bubble phase and surface area of the bubble-matrix interface were determined during growth and interconnection. Surface area increased rapidly during the initial stages of growth, then slowed and finally decreased as bubble interconnection began and coarsening acted to reduce surface area. The fraction of the bubbles vented to a simulation domain boundary, f V , was quantified as a measure of the microstructure’s interconnectivity and plotted as a function of porosity p . The defect species diffusivity was varied; although changes in diffusivity significantly affected the microstructure, the plots of f V vs. p did not change significantly. The percolation threshold p c was calculated to be approximately 0.26, depending on the assumed diffusivity and using an initial bubble number density based on experimental observations. This is slightly smaller than the percolation threshold for continuum percolation of overlapping 3D spheres. The simulation results were used to parameterize two different engineering-scale swelling models for U-(Pu)-Zr in the nuclear fuel performance code BISON.

Aagesen, Larry K. (ORCID:000000034936676X)↗