Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix elements”

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

Why Is Attention Sparse In Particle Transformer?

Transformer-based models have achieved state-of-the-art performance in jet tagging at the CERN Large Hadron Collider (LHC), with the Particle Transformer (ParT) representing a leading example of such models. A striking feature of ParT is its sparse, nearly binary, attention structure, raising questions about the origin of this behavior and whether it encodes physically meaningful correlations. In this work, we investigate the source of ParT's sparse attention by comparing models trained on multiple benchmark datasets and examine the relative contributions of the attention term and the physics-inspired interaction matrix before softmax. We find that binary sparsity arises primarily from the attention mechanism itself, with the interaction matrix playing a secondary role. Moreove, we show that ParT is able to identify key jet substructure elements, such as leptons in semileptonic top decays, even without explicit particle identification inputs. These results provide new insight into the interpretability of transformer-based jet taggers and clarify the conditions under which sparse attention patterns emerge in ParT.

Legge, Timothy [UC, San Diego]↗

Impact of Molten Gallium on the Microstructure and Corrosion Behavior of Aluminum and Uranium-Aluminum Alloys for Used Nuclear Fuel Reprocessing

Test reactors around the world utilize highly enriched uranium fuel to achieve high neutron fluxes for materials testing. Once spent, the remaining uranium is a valuable resource for subsequent fuel fabrication. However, some of these test reactor cores consist of curved plate-type fuel elements, fabricated using aluminum alloy 6061 (AA6061) cladding to encapsulate a uranium-aluminum alloy (UAlx) fuel matrix. These assemblies require non-standard reprocessing approaches for uranium recovery, as aluminum dissolves readily in acidic solutions, generating large volumes of waste and complicating downstream chemical separations. In this work, we investigate a novel chemical decladding strategy based on the interaction between AA6061/UAlx and molten gallium (Ga). Ga is known to induce severe degradation of aluminum metal through liquid metal embrittlement (LME), even at relatively low Ga concentrations. By penetrating the aluminum crystal lattice, Ga disrupts grain cohesion and facilitates fracture or dissolution of the aluminum matrix. Thermodynamic analysis of the Al–Ga binary phase diagram suggests that Ga may offer a viable pathway to selectively weaken or dissolve the AA6061 cladding, and potentially the aluminum component of the UAlx fuel matrix within. To this end, parametric experiments were performed at 50 °C and 100 °C across a range of Al–Ga atomic fractions. At lower Al fractions, the AA6061 was completely molten after 2 hours of exposure to the Ga metal. In contrast, samples with higher Al fractions (0.9 Al, 0.1 Ga) contained residual solids after 2 hours, which were characterized by microstructural examination using electron backscatter diffraction (EBSD) and transmission electron microscopy (TEM). These Al-Ga compositions were also evaluated using FactSage thermodynamic modeling to further elucidate the relationship between phase diagram behavior and LME.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Simulation-driven optimization of high-order meshes in ALE hydrodynamics

Here we propose tools for high-order mesh optimization and demonstrate their benefits in the context of multi-material Arbitrary Lagrangian-Eulerian (ALE) compressible shock hydrodynamic applications. The mesh optimization process is driven by information provided by the simulation which uses the optimized mesh, such as shock positions, material regions, known error estimates, etc. These simulation features are usually represented discretely, for instance, as finite element functions on the Lagrangian mesh. The discrete nature of the input is critical for the practical applicability of the algorithms we propose and distinguishes this work from approaches that strictly require analytical information. Our methods are based on node movement through a high-order extension of the Target-Matrix Optimization Paradigm (TMOP). The proposed formulation is fully algebraic and relies only on local Jacobian matrices, so it is applicable to all types of mesh elements, in 2D and 3D, and any order of the mesh. We discuss the notions of constructing adaptive target matrices and obtaining their derivatives, reconstructing discrete data in intermediate meshes, node limiting that enables improvement of global mesh quality while preserving space-dependent local mesh features, and appropriate normalization of the objective function. The adaptivity methods are combined with automatic ALE triggers that can provide robustness of the mesh evolution and avoid excessive remap procedures. The benefits of the new high-order TMOP technology are illustrated on several simulations performed in the high-order ALE application BLAST.

97 MATHEMATICS AND COMPUTING↗

Extending a 10‐node composite tetrahedral finite element for solid mechanics

Abstract We propose to extend the composite tetrahedral finite element first introduced by Thoutireddy et al. (2002) and recently reformulated by Ostien et al. (2016). We generalize the gradient operator and mass matrix to curved domains through analytical expressions weighted by subtetrahedra Jacobians. Optimal integration weights are constructed to increase the accuracy of Gaussian quadrature. We preserve the variational structure of the formulation through a new five‐field functional with additional, independent fields for the Jacobian and the pressure. This approach not only obviates volumetric locking but also yields symmetry. A deleterious soft mode, common to both the quadratic and composite tetrahedral element with constant pressure formulations is effectively stabilized through a novel convex energy penalty function. Numerous numerical examples spanning a patch test to the impact of a Taylor bar demonstrate the accuracy, robustness, and convergence of the extended composite tetrahedral element for application to structural metals.

Foulk III, James W.↗

Collaborative Study for Certification of Trace Elements in Uranium Ore Concentrate CRMs UCLO-1, UCHI-1, and UPER-1

Trace impurity patterns are important nuclear forensic signatures in uranium ore concentrates (UOCs) and Certified Reference Materials (CRMs) are used to validate the analysis methods employed by end users. Herein, we discuss the certification campaign for three new UOC CRMs from the National Research Council Canada: UCLO-1 (https://doi.org/10.4224/crm.2020.uclo-1), UCHI-1 (https://doi.org/10.4224/crm.2020.uchi-1), and UPER-1 (https://doi.org/10.4224/crm.2020.uper-1). This study involved 15 laboratories from 10 countries, using sector-field and (triple) quadrupole inductively-coupled plasma mass spectrometry to analyze 64 trace element impurities. We discuss the importance of the acids used for sample digestion, difficulties analyzing in a high uranium matrix, and data combination and uncertainty evaluation for this large dataset.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Draft Experimental Plan for Commercial SNF Degradation in Repository Environments with a Focus on Fracture Matrix Degradation

This report documents work performed supporting US Department of Energy (DOE) Nuclear Energy Spent Fuel and Waste Disposition, Spent Fuel and Waste Science and Technology, under work breakdown structure element 1.08.01.03.04, “Inventory and Waste Form Characteristics and Performance.” In particular, this report fulfills the M3 milestone, M3SF-21OR010309072 "Draft Experimental Plan for Commercial SNF Degradation in Repository Environments” as described in work package SF-21OR01030907 "SNF Degradation Testing-ORNL".This experimental plan is a collaboration between ORNL, SNL, ANL, and PNNL. It is intended to support source term model development and spent fuel degradation testing in a range of generic repository conditions. Potential testing and experiments to support validation of the Fuel Matrix Dissolution Model (FMDM) and the overall development/validation of repository source term modeling are discussed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Multiscale characterization and representation of variability in ceramic matrix composites

Low density, high strength, and high creep and oxidation resistance properties of ceramic matrix composites (CMCs) make them an ideal choice for use in extreme environments in space and military applications. This paper presents a detailed characterization study of structural and manufacturing flaws in Carbon fiber Silicon-Carbide-Nitride matrix (C/SiNC) CMCs at different length-scales. Energy-dispersive spectroscopy (EDS) is used for the chemical characterization of the material’s elemental constituents. High-resolution multiscale graphs obtained from scanning electron microscope (SEM) and confocal laser scanning microscope (LSM) are used to characterize the distribution and morphology of defects at different length scales. This is followed by the classification and quantification of the common manufacturing defects. An image processing algorithm based on the image segmentation process is developed to quantify the variability of various scale-dependent architectural parameters. Finally, a three-dimensional stochastic representative volume element (SRVE) generation algorithm is developed to provide precise representations of material textures at multiple length scales. The developed algorithm accurately accounts for material features and flaws based on a range of multiscale structural and defects characterization results.

36 MATERIALS SCIENCE↗

Neutron Irradiation Induced Changes in Isotopic Abundance of 6Li and 3D Nanoscale Distribution of Tritium in LiAlO2 Pellets Analyzed by Atom Probe Tomography

Tritium, a radioactive isotope of hydrogen is produced by neutron irradiation of 6Li enriched LiAlO2 pellets in Tritium producing burnable absorber rods (TPBARs). Three-dimensional nanoscale mapping of Tritium in irradiated pellets is critical to understand its spatial variation and its association with various microstructural features in the pellets. Spatially resolved analysis of 6Li isotopic enrichment in such ceramic pellets before and after irradiation can provide insights to heterogeneities in Li isotopic distribution. However, Li being a light isotope, its spatially resolved compositional or isotopic analysis is a challenging task to most analytical methods. Here we used atom probe tomography, to evaluate the Li isotopic enrichment in both the LiAlO¬2 matrix and the Li deficient secondary phase LiAl5O8 before irradiation. These results were then compared with 6Li isotopic ratio measured by APT analysis of neutron irradiated pellets to clearly demonstrate the unique capability of APT to quantitatively analyze the isotopic enrichment of 6Li as well as for other light elements present in these materials. Evidence for heterogenous nanoscale distribution of 3H and O3H within irradiated pellet microstructure is also provided, which is attributed to 3H trapping in irradiation induced vacancy clusters or voids within the LiAlO2 matrix. The characterization results given here, now prove the feasibility for using APT for analyzing nanoscale spatially resolved isotopic enrichment of Li as well as other light elements such as Tritium in TPBARs. This work also paves way to adopt APT for analyzing light elements such as Li and tritium in materials used for nuclear fusion reactor applications, geochemistry, astrophysics as well as for energy storage applications.

Devaraj, Arun↗

Design and experimental testing of a 150 kWh thermal battery using thermosiphons embedded in a concrete matrix for power plant flexible operation

One of the options for achieving the global temperature limitation of 1.5 °C target, for the mitigation of global warming, is based on the better penetration of renewables into the electrical grid. This has imposed a burden to fossil fuel fired power plants since they are required to operate away from their baseload mode to compensate for the inherent intermittence of the renewable power. Integrating energy storage with fossil plants is an option to achieve their needed flexibility. A cost competitive energy storage option for the solution is based on storing sensible heat in concrete. Here, this paper reports research results and development of a thermal battery cell (TBC) capable of operating at temperatures up to 425 °C. A novel concept consisting of a concrete matrix for sensible heat storage, engineered to provide enhanced thermal and mechanical properties, and twenty-two thermosiphon elements, engineered for dual action were designed and fabricated into a single thermal energy storage (TES) module. Research for the development of the components for the TBC was performed in the laboratory. Efficient heat transfer, to/from the storage media, was demonstrated under several charging and discharging conditions with a thermal storage capacity of 150 kWh th and a rapid discharge, making the TBC suitable for fast ramping when integrated with a fossil fuel fired power plant. Efficient radial heat transfer to the concrete was observed due to the well designed spacing and location of thermosiphons in the radial direction. A minimal temperature difference of 2 °C, between the thermosiphons bottom and top was obtained, demonstrating the isothermicity of those elements. An overall end-to-end TBC energy-to-energy round trip efficiency of 70% was achieved.

25 ENERGY STORAGE↗

SparseLU, A Novel Algorithm and Math Library for Sparse LU Factorization

Decomposing sparse matrices into lower and upper triangular matrices (sparse LU factorization) is a key operation in many computational scientific applications. We developed SparseLU, a sparse linear algebra library that implements a new algorithm for LU factorization on general sparse matrices. The new algorithm divides the input matrix into tiles to which OpenMP tasks are created for factorization computation, where only tiles that contain nonzero elements are computed. For comparative performance analysis, we used the reference library SuperLU. Testing was performed on synthetically generated matrices which replicate the conditions of the real-world matrices. SparseLU is able to reach a mean speedup of ~29× compared to SuperLU.

Valero Lara, Pedro↗

Sparse Approximate Multifrontal Factorization with Butterfly Compression for High-Frequency Wave Equations

In this work, we present a fast and approximate multifrontal solver for large-scale sparse linear systems arising from finite-difference, finite-volume or finite-element discretization of high-frequency wave equations. The proposed solver leverages the butterfly algorithm and its hierarchical matrix extension for compressing and factorizing large frontal matrices via graph-distance guided entry evaluation or randomized matrix-vector multiplication-based schemes. Complexity analysis and numerical experiments demonstrate $\mathcal{O}(N\log^2 N)$ computation and $\mathcal{O}(N)$ memory complexity when applied to an $N\times N$ sparse system arising from 3D high-frequency Helmholtz and Maxwell problems.

97 MATHEMATICS AND COMPUTING↗

INL CRMO XRF Special Report No. 3: INL CRMO XRF Calibration Methods for the Analysis of Ceramics and Brick

This report details the analytical methods used by the Idaho National Laboratory Cultural Resource Management Office to characterize ceramic and brick artifacts using X-ray fluorescence spectrometry. More specifically, this report provides details on instrumentation, analysis times, and matrix-specific calibration protocols for ceramic materials. Calibration curves were developed to adjust raw instrumental output for 13 elements using a suite of reference standards made from historic brick borrowed from the Yale Peabody Museum. Calibration curves were validated by plotting predicted values against published reference values for the calibration set as well as a certified soil reference standard (SRM 2711a).

36 MATERIALS SCIENCE↗

A stable potential-based time-domain method for wideband elec- tromagnetic analysis

In previous research, the frequency-domain A-ϕ formulation has been validated using the finite element method for electromagnetic simulations of low-frequency and multi- scale problems, demonstrating excellent numerical accuracy, good matrix condition, and high computational efficiency. Time- domain simulations provide significant advantages for modeling wideband problems and are crucial for multiphysics applications. In this paper, the frequency-domain A-ϕ formulation is extended to the time domain. The central difference scheme is employed for temporal discretization to ensure both accuracy and stability. A numerical example is presented to demonstrate the capability of the proposed time-domain method in wideband electromagnetic analysis.

Mekonnen, Minyichil↗

An initial exploration of Bayesian model calibration for estimating the composition of rocks and soils on Mars

The Mars Curiosity rover carries an instrument, ChemCam, designed to measure the composition of surface rocks and soil using laser-induced breakdown spectroscopy (LIBS). The measured spectra from this instrument must be analyzed to identify the component elements in the target sample, as well as their relative proportions. This process, which we call disaggregation, is complicated by so-called matrix effects, which describe nonlinear changes in the relative heights of emission lines as an unknown function of composition due to atomic interactions within the LIBS plasma. In this work, we explore the use of the plasma physics code ATOMIC, developed at Los Alamos National Laboratory, for the disaggregation task. ATOMIC has recently been used to model LIBS spectra and can robustly reproduce matrix effects from first principles. The ability of ATOMIC to predict LIBS spectra presents an exciting opportunity to perform disaggregation in a manner not yet tried in the LIBS community, namely via Bayesian model calibration. However, using it directly to solve our inverse problem is computationally intractable due to the large parameter space and the computation time required to produce a single output. Therefore, we also explore the use of emulators as a fast solution for this analysis. We discuss a proof of concept Gaussian process emulator for disaggregating two-element compounds of sodium and copper. The training and test datasets were simulated with ATOMIC using a Latin hypercube design. After testing the performance of the emulator, we successfully recover the composition of 25 test spectra with Bayesian model calibration.

97 MATHEMATICS AND COMPUTING↗

Constitutive Modeling of the Cross-Linked Polymers During the Trio-Aging Mechanism: Effects of Water, Oxygen, UV

Abstract The synergized effects of diverse environmental conditions such as humidity, temperature, UV, and oxygen, are considered through a generic model by combining the concepts of network alteration, hygrothermal, and Photo-thermo oxidation Models [1–4]. Hygrothermal and Photo-oxidation are considered dual aging mechanisms which are due to the synergized effects of moisture-oxygen, and UV-oxygen, respectively at elevated temperatures. Here Trio-aging is considered as the synergized effects of Hygrothermal, and Photo-oxidation which the main environmental elements involved in polymer degradation are oxygen, water, and UV radiation where each of which has a different micro-structural impact on the polymer matrix. While oxygen imposes chain scission and cross-link formation, water, and UV radiation mainly cause chain scission. Similar to our recent models, the dual network hypothesis of Topolosky is utilized and as a result, the strain energy of the polymer matrix is written in terms of two states namely unaged and infinity. The infinity state network, itself, is decomposed into three networks, namely brittle, flexible, and dead networks. While the brittle network models the cross-link formation as sub-structural phenomena, the flexible and dead networks are utilized to represent the cross-link reduction and chain scission. The contribution of these networks in the infinity state is based on the amount of relevant environmental elements. The model is developed for a domain with plenty of oxygen and holds true for a gradual aging process that happens in super-thin samples.

Bahrololoumi, Amir↗

Discrete Green’s functions and spectral graph theory for computationally efficient thermal modeling

Here, this work concerns solutions of the heat equation with the spectral graph method, for which the temperature is defined at discrete points in the domain and the spatial relationship among the points is described by a graph. The heat equation on the graph is solved using matrix techniques involving the eigenvectors and eigenvalues of the Laplacian matrix. The spectral graph approach precludes the computationally intensive meshing and numerous time-integration steps of the finite element method. In the present work, the spectral graph method is extended to include heat loss at the boundaries with a generalized boundary condition, and physics-based edge weights are introduced which simplify the calibration process. From this approach a discrete Green’s function is defined which allows for solutions under a variety of heating conditions including: space-varying initial conditions; time-and-space varying internal heating; and, time-and-space-varying heating at boundaries of type 1 (Dirichlet), type 2 (Neumann) and type 3 (Robin). Results are provided for benchmark heat transfer problems in one spatial dimension and in three spatial dimensions, and verification is provided by comparison with exact analytical solutions and finite difference solutions. The spectral graph method converges within 0.4% error of the analytical solution. The practical utility of the approach is demonstrated by thermal simulation of a multilayer additive manufacturing process. The spectral graph results are compared to experimentally-obtained temperature data for two metal parts, with error less than 5% of the experimental measurements, with computation time less than one minute on a desktop computer.

36 MATERIALS SCIENCE↗