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 739 records · Page 41

System and method for focal-plane illuminator/detector (FASID) design for improved graded index lenses

The present disclosure relates to a detector system for imaging an optical signal received by a graded index (GRIN) optical element to account for known variations in a graded index distribution of the GRIN optical element. The detector system incorporates a plurality of optical detector elements configured to receive optical rays received by the GRIN optical element at specific locations on a plane of the GRIN optical element. Ray tracing software is configured to receive and map the optical rays to a plurality of additional specific locations on the plane based on the known variations in the graded index distribution of the GRIN optical element. A processor uses algorithms for diagonalization of a linear system matrix to determine information on both an intensity and an angle of the received optical rays at each one of the plurality of specific locations on the plane.

Bude, Jeffrey Devin↗

Nuclear waste reduction: Exploring new pathways one step at a time

In my home country of Venezuela, nuclear energy is not a topic that attracts much attention. The government briefly oversaw some nuclear energy programs during the 1950s, but currently there are no active nuclear power facilities in the country. In fact, the Venezuelan government signed and ratified the treaty of the prohibition of nuclear weapons in 2021, which states that Venezuela has never owned, possessed or controlled nuclear weapons or programs of any kind. When I moved to the United States, however, nuclear energy became an extremely relevant topic. In the 1940s, the U.S. government established and oversaw the Manhattan Project to build atomic bombs for use in World War II. After the war, the government encouraged scientists to use this information on nuclear reactions to develop nuclear energy for peaceful civilian purposes instead.1 During these early days of nuclear research, there were no formal regulatory standards for nuclear waste management. Policies usually were self-regulated and often created based on existing policies of disposal for non-nuclear waste.2 As a result, there were instances of nuclear waste leaching into the environment and affecting local communities. So, much research has been conducted since then to characterize and store nuclear waste safely and securely.3 I first became interested in nuclear energy during my undergraduate studies when I worked on a project involving ligand synthesis to help extract actinides from nuclear waste. I then studied electrochemistry in molten salt systems for nuclear energy applications during my Ph.D. As I approached graduation, I started looking into national laboratories that have programs involving nuclear energy and waste management. At Idaho National Laboratory (INL), the focus is more on applied processes and how nuclear energy can be innovated to realize next-generation reactor design and technologies. This focus led me to apply for a Seaborg distinguished postdoctoral position at INL, for which I was chosen based on my proposal of a way to improve nuclear waste recycling. To understand my proposal, we must familiarize ourselves with the makeup of nuclear waste. After uranium dioxide is used as nuclear fuel in a reactor, the fuel matrix is then characterized by various fission products, including rare earth elements, alkali and alkaline earths, and actinides. Some of these fission products can potentially be recovered through pyroprocessing, 4 which involves the electrochemical dissolution of the used nuclear fuel in a molten chloride salt mixture at high temperatures. Though some of the fission products can be easily recovered—for example, uranium is reduced onto an inert cathode by applied potentials—numerous other fission products such as rare earth elements are difficult to recover due to their multivalent oxidation states and side reactions.5 To improve the recovery efficiency of rare earth elements specifically, I proposed investigating the fundamental interactions between rare earth elements in the molten chloride salt and their metallic form (Figure 1). The kinetic pathways and the chemical reactions of these elements, which will be elucidated through spectro-electrochemistry at high temperatures, will give insights on how the recovery efficiency can be improved. Although my research focuses on fundamental science, it will benefit the applied process by generating new scientific knowledge and closing the gap for efficient recycling of the waste: one step at a time.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Probing for the Trace Estimation of a Permuted Matrix Inverse Corresponding to a Lattice Displacement

We report thatpProbing is a general technique that is used to reduce the variance of the Hutchinson stochastic estimator for the trace of the inverse of a large, sparse matrix A. The variance of the estimator is the sum of the squares of the off-diagonal elements of A -1 . Therefore, this technique computes probing vectors that when used in the estimator annihilate the largest off-diagonal elements. For matrices that display decay of the magnitude of |A$^{-1}_{ij}$| with the graph distance between nodes i and j, this is achieved through graph coloring of increasing powers A k . Equivalently, when a matrix stems from a lattice discretization, it is computationally beneficial to find a distance-k coloring of the lattice. Previously, a hierarchical coloring was proposed so that k can be increased at runtime as needed without discarding previous work. In this work, we study probing for the more general problem of computing the trace of a permutation of A -1 , say PA -1 . The motivation comes from lattice quantum chromodynamics (QCD), where we need to construct “disconnected diagrams” to extract flavor-separated generalized parton functions. In lattice QCD, where the matrix has a four-dimensional toroidal lattice structure, these nonlocal operators correspond to a PA -1 , where P is the permutation relating to some displacement $\vec{p}$ in one or more dimensions. We focus on a single dimension displacement (p), but our methods are general. We show that probing on A k or (PA) k does not annihilate the largest magnitude elements. To resolve this issue, our displacement-based probing works on PA k using a new coloring scheme that works directly on appropriately displaced neighborhoods on the lattice. We prove lower bounds on the number of colors needed and study the effect of this scheme on variance reduction, both theoretically and experimentally on a real-world lattice QCD calculation. We achieve orders of magnitude speedup over the unprobed or the naively probed methods.

97 MATHEMATICS AND COMPUTING↗

A predictive modeling tool for damage analysis and design of hydrogen storage composite pressure vessels

In this study, a predictive modeling tool is developed for damage analysis and design of hydrogen (H 2 ) storage composite pressure vessels. It integrates micromechanics of matrix cracking into a continuum damage mechanics (CDM) description for damage evolution, and three-dimensional (3D) finite element (FE) modeling of the vessel structural response. At the scale of the composite layer (mesoscale), the temperature-dependent stiffness reduction law in terms of the damage variable for transverse matrix cracking is computed using an Eshelby-Mori-Tanaka approach for the initial composite thermoelastic properties and a self-consistent model for the stiffness reduction as a function of the damage variable. While transverse matrix cracking obeying a damage evolution relation can progressively evolve from an initiation to a saturation state, fiber failure is predicted by a micromechanical fiber rupture criterion that accounts for the fiber strength and matrix stress. The implementation of this integrated multiscale modeling model into a 3D FE formulation enables damage analysis and design of H 2 storage composite pressure vessels. The developed tool is illustrated through 3D damage analyses of a cryogenically compressed H 2 storage vessel model subjected to thermomechanical loadings to investigate effects of the helical layer fiber orientation and loading scenario on damage development, vessel integrity and burst pressure.

08 HYDROGEN↗

A comparison between the one- and two-step spin-orbit coupling approaches based on the ab initio Density Matrix Renormalization Group

The efficient and reliable treatment of both spin-orbit coupling (SOC) and electron correlation is essential for understanding f-element chemistry. Here, we analyze two approaches to the problem, the one-step approach where both effects are treated simultaneously, and the two-step state interaction approach. We report an implementation of the ab initio density matrix renormalization group(DMRG) with a one-step treatment of the SOC effect which can be compared to prior two-step treatments on an equal footing. Using a dysprosium octahedral complex and bridged dimer as benchmark systems, we identify characteristics of problems where the one-step approach is beneficial for obtaining the low-energy spectrum.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ocean Energy Systems Wave Energy Modeling Task 10.4: Numerical Modeling of a Fixed Oscillating Water Column

This paper reports on an ongoing international effort to establish guidelines for numerical modeling of wave energy converters, initiated by the International Energy Agency Technology Collaboration Program for Ocean Energy Systems. Initial results for point absorbers were presented in previous work, and here we present results for a breakwater-mounted Oscillating Water Column (OWC) device. The experimental model is at scale 1:4 relative to a full-scale installation in a water depth of 12.8 m. The power-extracting air turbine is modeled by an orifice plate of 1–2% of the internal chamber surface area. Measurements of chamber surface elevation, air flow through the orifice, and pressure difference across the orifice are compared with numerical calculations using both weakly-nonlinear potential flow theory and computational fluid dynamics. Both compressible- and incompressible-flow models are considered, and the effects of air compressibility are found to have a significant influence on the motion of the internal chamber surface. Recommendations are made for reducing uncertainties in future experimental campaigns, which are critical to enable firm conclusions to be drawn about the relative accuracy of the numerical models. It is well-known that boundary element method solutions of the linear potential flow problem (e.g., WAMIT) are singular at infinite frequency when panels are placed directly on the free surface. This is problematic for time-domain solutions where the value of the added mass matrix at infinite frequency is critical, especially for OWC chambers, which are modeled by zero-mass elements on the free surface. A straightforward rational procedure is described to replace ad-hoc solutions to this problem that have been proposed in the literature.

16 TIDAL AND WAVE POWER↗

hr-Adaptivity for nonconforming high-order meshes with the target matrix optimization paradigm

In this study, we present an \(hr\) -adaptivity framework for optimization of high-order meshes. This work extends the r -adaptivity method by Dobrev et al. (Comput Fluids, 2020), where we utilized the Target-Matrix Optimization Paradigm (TMOP) to minimize a functional that depends on each element’s current and target geometric parameters: element aspect-ratio , size , skew , and rotation . Since fixed mesh topology limits the ability to achieve the target size and aspect-ratio at each position, in this paper, we augment the r -adaptivity framework with nonconforming adaptive mesh refinement to further reduce the error with respect to the target geometric parameters. The proposed formulation, referred to as \(hr\) -adaptivity, introduces TMOP-based quality estimators to satisfy the aspect-ratio target via anisotropic refinements and size target via isotropic refinements in each element of the mesh. The methodology presented is purely algebraic, extends to both simplices and hexahedra/quadrilaterals of any order, and supports nonconforming isotropic and anisotropic refinements in 2D and 3D. Using a problem with a known exact solution, we demonstrate the effectiveness of \(hr\) -adaptivity over both r - and \(h\) -adaptivity in obtaining similar accuracy in the solution with significantly fewer mesh nodes. We also present several examples that show that \(hr\) -adaptivity can help satisfy geometric targets even when \(r\) -adaptivity fails to do so, due to the topology of the initial mesh.

42 ENGINEERING↗

Xylan-cellulose core structure of oat water-extractable β-glucan macromolecule: Insight into interactions and organization of the cell wall complex

Water-extractable β-glucan with high molar mass (HM) determines health benefits of oat food. Oat β-glucan was extracted by a standardized in vitro digestion method and co-existing water-extractable polysaccharide (WEP) fraction and its HM-arabinoxylan (HM-AX) subfraction were isolated to identify their highly acid-resistant subunit and investigate molecular interactions between constituent polymers. The WEP and HM-AX samples consisted of arabinoxylans (AXs) (74 and 76 %, respectively), however, cellulose constituted the secondary component (6.6 and 12.8 %, respectively). Multi-detection HPSEC along with specific enzymatic hydrolysis of AXs revealed the presence of the HM-xylan domain (16 and 34 %, respectively) built of numerous single- and multi-component populations with random coil and rod-like conformations, which were embedded in a xylan matrix with spherical conformation and controlled the macromolecular shape. Unlike single-component populations, the multi-component ones were resistant to hydrolytic action of AX-hydrolyzing enzymes and represented the subunits that anchor matrix polysaccharides onto cellulose surface. These results indicate that water-extractable β-glucan macromolecule comprises as integral element a cellulose core with two linking populations, HM-xylan and low molar mass glucomannan, which are surrounded by a feruloylated AX-arabinan-arabinogalactan composite and next laminated by β-glucan matrix. Here, the stiff cellulose-xylan backbone is the basis of HM β-glucan organization, controlled by its cellulose-like segments.

59 BASIC BIOLOGICAL SCIENCES↗

Hierarchical morphologies in co-sputter deposited thin films

Co-depositions of immiscible alloy films at specific processing conditions have yielded hierarchical microstructures which consist of distinct features at multiple length scales, often agglomerates and concentration modulations on the sub-micrometer-scale and fine nanoprecipitates in a matrix on the nanoscale. The present work examined a series of immiscible alloy systems: Cu-Mo, Cu-Ag, Cu-Fe, Cu-Ta, Mo-Ag, Cu-Mo-Ag, to determine the kinetic conditions favorable for hierarchical organization and the formation mechanism of such structures. Thin films of six immiscible systems were sputter co-deposited over a range of deposition rates from 0.12 to 2 nm/s and various temperatures from 400 to 800°C. The resulting microstructures indicate that hierarchical structures form with sufficient disparity in kinetic energy between the constituent atoms, one species being highly mobile (A) and the other relatively immobile (B). This condition arises typically at elevated deposition temperatures and reduced deposition rates but is also alloy dependent. The hierarchical structures form during deposition via phase separation and self-organization processes across the multiple length scales. The adatoms diffuse on the film surface with the highly mobile species swiftly agglomerating into A-rich domains within which B-rich nanoprecipitates form, often self-organizing into periodic arrays. The smallest B-rich nanoprecipitates in the A-rich domains are found to be coherent and in a metastable crystal structure (B taking the structure of A), but coarser precipitates that exhibit the equilibrium structure of B element. The A-rich domains are surrounded by a B-rich matrix that phase-separates into a concentration modulated structure. In conclusion, the observations are interpreted via a model incorporating material properties and process parameters.

36 MATERIALS SCIENCE↗

Designing Nuclear Fuels with a Multi-Principal Element Alloying Approach

Previous research has shown that multi-principal element alloys (MPEAs) using chromium, molybdenum, niobium, tantalum, titanium, vanadium, and zirconium can form stable body-centered-cubic (BCC) structures across a large temperature region (25°C to 1000°C). This is the same crystal structure as γ-uranium (U), which has shown desirable thermal and irradiation behavior in previous alloy fuel research. It is hypothesized then that the MPEA alloying approach can be used to produce a stable BCC uranium-bearing alloy and to retain its stability throughout anticipated operating regimes of power-producing reactors. Candidate elements were assessed using Monte Carlo N-Particle (MCNP) analysis to determine uranium densities necessary to make the alloy an economically viable fuel compared to conventional fuel forms. Following neutronic considerations, materials property databases and empirical predictors were used to determine the compositions with a high potential to form a BCC solid solution alloy. The final four alloys were MoNbTaU 2 , MoNbTiU 2 , NbTaTiU 2 , and NbTaVU 2 , which were cast using arc melting of raw elemental foils and chunks. Characterization of the fabricated alloys included scanning electron microscopy, X-ray diffraction, energy dispersive X-ray spectroscopy, and transmission electron microscopy. The results showed a two-phase system with a U-rich matrix phase surrounding the refractory precipitates. The U phase was found to contain varying concentrations of the alloying elements and was a BCC γ-U phase. These results warrant further research to identify ideal compositions for use as an advanced alloy fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Evaluating Spatial Accelerator Architectures with Tiled Matrix-Matrix Multiplication.

There is a growing interest in custom spatial accelerators for machine learning applications. These accelerators employ a spatial array of processing elements (PEs) interacting via custom buffer hierarchies and networks-on-chip. The efficiency of these accelerators comes from employing optimized dataflow (i.e., spatial/temporal partitioning of data across the PEs and fine-grained scheduling) strategies to optimize data reuse. The focus of this work is to evaluate these accelerator architectures using a tiled general matrix-matrix multiplication (GEMM) kernel. To do so, we develop a framework that finds optimized mappings (dataflow and tile sizes) for a tiled GEMM for a given spatial accelerator and workload combination, leveraging an analytical cost model for runtime and energy. Finally, our evaluations over five spatial accelerators demonstrate that the tiled GEMM mappings systematically generated by our framework achieve high performance on various GEMM workloads and accelerators.

42 ENGINEERING↗

Evaluating Spatial Accelerator Architectures with Tiled Matrix-Matrix Multiplication

There is a growing interest in custom spatial accelerators for machine learning applications. These accelerators employ a spatial array of processing elements (PEs) interacting via custom buffer hierarchies and networks-on-chip. The efficiency of these accelerators comes from employing optimized dataflow (i.e., spatial/temporal partitioning of data across the PEs and fine-grained scheduling) strategies to optimize data reuse. The focus of this work is to evaluate these accelerator architectures using a tiled general matrix-matrix multiplication (GEMM) kernel. To do so, we develop a framework that finds optimized mappings (dataflow and tile sizes) for a tiled GEMM for a given spatial accelerator and workload combination, leveraging an analytical cost model for runtime and energy. Our evaluations over five spatial accelerators demonstrate that the tiled GEMM mappings systematically generated by our framework achieve high performance on various GEMM workloads and accelerators.

43 PARTICLE ACCELERATORS↗

Ordered Particle Packing in Dense TRISO/SiC Fuel Elements and Preliminary Assessment of Neutronic and Thermomechanical Characteristics

Detailed analysis of the particle distribution in Transformational Challenge Reactor fuel elements indicates that particle packing is not random; instead, it follows a relatively ordered structure near fuel element surfaces. Discrete particle neutronic simulations indicate that the core reactivity is not impacted when assuming homogenization of particles with the silicon carbide matrix. However, the neutronic power distribution resulting from the ordered packing structure indicates that the highest-power particles reside at the top and bottom of the fuel elements and nearest the YH 1.85 moderator rods. The power distribution results were applied to thermo-mechanical simulations using mesh-based power distributions. Previous results indicated high stress at the bottom of the fuel element, where packing is most ordered. Additionally, to reduce this stress concentration, additively manufactured protrusions were added to the bottom of a test fuel element to disrupt dense particle packing. These protrusions reduced the overall power peaking, but the thermomechanical simulations did not indicate a significant change in the fuel element’s maximum stress or failure probability.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Total metals & anion concentration data; Slate River floodplain, Crested Butte, CO; May 2020-September 2020

This data package includes processed and undiluted measurements for metal and anion concentrations from pore water (groundwater) samples from the Slate River floodplain of Crested Butte, CO, a focus field site for the SLAC Floodplain Hydro-Biogeochemistry SFA. The data was generated as part of the work targeting the overarching research question for the SLAC SFA: How do ubiquitous subsurface interfaces mediate molecular-scale biogeochemical processes and groundwater quality in floodplains and watersheds? Samples were collected between May and September of 2020. These measurements were all recorded at the Arizona Laboratory for Emerging Contaminants (ALEC) at the University of Arizona located in Tucson, AZ. Groundwater samples were extracted from a network of installed rhizon (Rhizosphere Research Products, part no. 19.60.21F, 0.6 micrometer mesh size) and piezometer wells within the river floodplain. All water samples were shaded from sun exposure during extraction from the subsurface and preserved at 4C until measured at ALEC.Analysis by ICP-MS:Measurements for total metals were made on the Agilent 7700x ICP-MS (for total metals) – Agilent Technologies, Santa Clara, CA.The analytical QA/QC protocol was adapted from US EPA Method 200.8 for analysis by ICP-MS. Calibration standards were prepared from multi-element stock solution (Sigma-Aldrich Multielement standard solution for ICP, St. Louis, MO) using matrix matched to sample solutions (either 2% HCl or HNO3 from AriStar Plus,grade acids from VWR Scientific). Calibration curves include at least 7 points with correlation coefficients > 0.995. The QC protocol includes a continuing calibration blank (CCB), a continuing calibration verification (CCV) solution and at least one quality control sample (QCS) to be analyzed just after calibration and again after every 12 samples and at the completion of the run. The QCS solutions are from an independent source, such as NIST SRM 1643e - Trace Elements in Water, or QCS solutions from High Purity Standards (Charleston, SC). Acceptable QC responses must be between 90 and 110% of the certified value. An internal standard (Rh) is added via on-line addition into the sample line using a mixing tee.Analysis by Ion Chromatography (Anions):The protocol follows Method 4110 in Standard Methods for Examination of Water and Wastewater.The instrument used is the Thermo Scientific Dionex ICS-6000 using AS+AG22 column set for anion analysis with isocratic method using sodium carbonate eluent. Detection is by chemical suppression of eluent conductivity. Quality control solutions and mixed analyte standards purchased from Inorganic Ventures, Christiansburg, VA.All files are in csv format.

54 ENVIRONMENTAL SCIENCES↗

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↗

Alloying of U-Al-SS as a Simulant for Pu-Al-SS Alloying

The Savannah River National Laboratory evaluated several options for disposition of stainless steel (SS)–clad plutonium metal, particularly Pu-10.6 at. % Al (Pu- 1.3 wt% Al) alloy fuel. One technology considered was alloying fuel with SS. The goal of the alloying would be to make a SS-Pu alloy that was a nonproliferable waste form with secondary Pu-rich microencapsulated regions distributed throughout the refractory SS. The microencapsulation of the Pu regions should therefore allow the waste form to meet the requirements for a low attractiveness waste as defined by the U.S. Department of Energy. Plutonium-bearing alloys at these levels could potentially be suitable for disposal at a waste isolation pilot plant. Four metal ingots were successfully fabricated using U and Al as a surrogate for Pu-Al. The U was distributed and microencapsulated by the alloy matrix, thereby setting the stage for subsequent tests using SS-clad fuel elements containing Pu-10.6Al.

36 MATERIALS SCIENCE↗