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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Accelerating high-order mesh optimization using finite element partial assembly on GPUs

In this paper we present a new GPU-oriented mesh optimization method based on high order finite elements. Our approach relies on node movement with fixed topology, through the Target-Matrix Optimization Paradigm (TMOP) and uses a global nonlinear solve over the whole computational mesh, i.e., all mesh nodes are moved together. A key property of the method is that the mesh optimization process is recast in terms of finite element operations, which allows us to utilize recent advances in the field of GPU-accelerated high order finite element algorithms. For example, we reduce data motion by using tensor factorization and matrix-free methods, which have superior performance characteristics compared to traditional full finite element matrix assembly and offer advantages for GPU based HPC hardware. Furthermore, we describe the major mathematical components of the method along with their efficient GPU-oriented implementation. In addition, we propose an easily reproducible mesh optimization test that can serve as a performance benchmark for the mesh optimization community.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Uncertainty quantification for competing failure mechanisms in unidirectionally reinforced carbon–carbon composites

Microstructure-informed finite element models play a key role in the carbon–carbon composite design process. Variability in manufacturing process parameters and experimental limitations introduce model parameter uncertainty. This study quantifies the effect of model parameter uncertainty on transverse tensile fracture behavior and proposes a methodology to predict the failure mode based on competing microscale damage mechanisms. Finite element simulations incorporate fiber–matrix interface debonding with cohesive zones and matrix damage with a smeared crack band approach in a unidirectional carbon–carbon composite. Results from a variance-based global sensitivity analysis identifies interfacial and matrix damage parameters as the primary source of variability in fracture behavior. Sobol’ indices indicate that matrix and cohesive zone strengths contribute 94% of the variance in the effective ultimate stress. A local analysis elucidates the relationship between these constituent strength parameters and failure mode by estimating the probability of cohesive, matrix, and mixed-mode dominated failure. Based on the results for 4000 simulations, 93% exhibit mixed-mode or interfacial dominated failure, which underscores the crucial role of fiber–matrix interface debonding in the transverse tensile failure of carbon–carbon composites. These uncertainty quantification results facilitate more efficient model calibration and provide a framework for microstructure-informed failure predictions in the face of manufacturing-induced uncertainty.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Postirradiation characterization of palladium as an additive for fuel cladding chemical interaction mitigation in metallic fuel

This work describes the microstructural and elemental characterization of irradiated metallic fuels containing palladium as an additive. The use of additives has been proposed to control Fuel-Cladding Chemical Interaction (FCCI) and thus to promote higher fuel utilization (i.e., higher burnup). In this work, Pd has been investigated as a potential additive to metallic fuel to bind lanthanides, impeding their migration and attack on the cladding. The influence of Pd on the microstructure, chemistry and performance of metallic fuel has been characterized via scanning electron microscopy for two metallic fuel designs—namely, annular and solid fuel. Pd was observed to play an important role in the chemistry of the fuel. Indeed, the addition of Pd leads to the formation of new phases. Pd was detected to combine not only with the lanthanides, as intended, but also with Zr, a main element of the fuel matrix. While Pd proved to be effective in preventing lanthanide migration and their attack on the cladding, the Pd-Zr compound may potentially lead to other unexpected fuel-performance issues, such as the formation of low-melting point phases and increased unalloyed U available for FCCI interaction with Fe in the cladding. Even the increase of Zr to 13wt%. did not completely mitigate this adverse phenomenon generated by the Pd-Zr interaction. Furthermore, the efficacy of using this additive needs further investigation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

The j and k dependencies of the line coupling and line mixing effects: Theoretical studies of the relaxation matrices of N 2 -broadened CH 3 D

Line coupling and line mixing effects in parallel and perpendicular bands of CH 3 D perturbed by N 2 have been studied. The work focuses on exhibiting the j and k dependencies of these two processes. The calculations were based on a previously reported anisotropic intermolecular potential including both the long-range multipole, induction, and dispersion forces and a short-range atom-atom model. Throughout this work, it is shown that components with L 1 = 3 of the atom-atom model are dominant contributions to the diffusion operator. As a consequence, in comparison with other molecular systems such as the CH 3 Cl-N 2 and CH 3 I-N 2 , theoretically predicted line coupling and line mixing effects exhibit completely new j and k dependencies. In general, the theoretically calculated halfwidths and intra-doublets’ off-diagonal elements of the relaxation matrix are in reasonable agreement with measurements.

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

Off-line commissioning of the St. Benedict radiofrequency quadrupole cooler-buncher

The St. Benedict ion trapping system, which aims to measure the β - v angular correlation parameter in superallowed-mixed mirror transitions, is under construction at the University of Notre Dame. These measurements will provide much-needed data to improve the accuracy of the V ud element of the CKM matrix. One of the major components of this system is the radio frequency quadrupole cooler-buncher, which is necessary to create low-emittance ion bunches for injection into the measurement Paul trap. The off-line commissioning of the cooler-buncher, using a potassium ion source, determined that the device could produce cooled ion bunches characterized by a 50-ns full-width-half-maximum time width at its exit. The commissioning results also determined the trapping efficiency to be 93(1)% and the trapping half-life to be 20.0(5) s.

Paul trap↗
DOE OSTI2021

Numerical simulation of parametric resonance in point absorbers using a simplified model

Abstract Parametric resonance is a non‐linear phenomenon in which a system can oscillate at a frequency different from its exciting frequency. Some wave energy converters are prone to this phenomenon, which is usually detrimental to their performance. Here, a computationally efficient way of simulating parametric resonance in point absorbers is presented. The model is based on linear potential theory, so the wave forces are evaluated at the mean position of the body. However, the first‐order variation of the body's centres of gravity and buoyancy is taken into account. This gives essentially the same result as a more rigorous approach of keeping terms in the equation of motion up to second order in the body motions. The only difference from a linear model is the presence of non‐zero off‐diagonal elements in the mass matrix. The model is benchmarked against state‐of‐the‐art non‐linear Froude–Krylov and computational fluid dynamics models for free decay, regular wave, and focused wave group cases. It is shown that the simplified model is able to simulate parametric resonance in pitch to a reasonable accuracy even though no non‐linear wave forces are included. The simulation speed on a standard computer is up to two orders of magnitude faster than real time.

16 TIDAL AND WAVE POWER↗
DOE OSTI2020

Trace Elemental Analysis of Bulk Thorium Using an Automated Separation–Inductively Coupled Plasma Optical Emission Spectroscopy Methodology

Presented here is a novel automated method for determining the trace element composition of bulk thorium by inductively coupled plasma–optical emission spectroscopy (ICP-OES). ICP-OES is a universal approach for measuring the trace elemental impurities present in actinide-rich materials; however, due to the emission rich spectrum of the actinide, a separation from the trace elements is warranted for spectrochemical analysis. Here, AG MP-1 ion exchange resin was utilized for retention of the Th matrix, while allowing the trace element impurities to be separated prior to subsequent analysis using ICP-OES. After demonstrating the separation on traditional gravity-driven columns, the methodology was transitioned to an automated platform for comparison. This automated platform utilizes syringe-driven sample and solvent flow and can collect the trace element and thorium fractions in separate locations. Finally, while reducing the sample size (500 µL, 1.5 mg of Th), maintaining the overall separation efficiency (recoveries >95%), and illustrating the sample throughput ability (n = 10+), this automated methodology could be readily adopted to nuclear facilities in which the determination of trace elemental impurities in Th samples is warranted.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Structure of the polysaccharide sheath from the B race of the green microalga Botryococcus braunii

The green microalga Botryococcus braunii forms a colony of cells held together by an intricate extracellular matrix. This extracellular matrix also acts as a storage location for liquid hydrocarbons initially produced inside the cells. These hydrocarbons can be used as feedstocks to produce transportation fuels. In the B race of B. braunii, the extracellular matrix is composed of several elements: a cross-linked hydrocarbon network filled with the liquid hydrocarbons, a retaining wall to hold the hydrocarbons in the extracellular matrix, an arabinogalactan polysaccharide fibrillar sheath that extends outward from the retaining wall, and a single protein known as polysaccharide associated protein is found at the base of the polysaccharide fibers where they connect to the retaining wall. During growth, colonies of the B race of B. braunii shed fragments of the retaining wall with polysaccharide fibrils attached into the media. Here we describe the isolation and structural analysis of these polysaccharides using a series of chemical degradations in combination with glycosyl residue composition and linkage analysis, as well as 2-dimensional nuclear magnetic resonance spectroscopy. The results demonstrate that the polysaccharide is highly branched, consisting of a galactan backbone with short branches comprised of Dgalactopyranose, L-arabinofuranose, 4-O-methyl-D-glucopyranuronic acid, and 6-deoxy-D-altropyranose. To our knowledge, this is the first description of 6-deoxyaltrose in algae or plants.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

A matrix-free approach for finite-strain hyperelastic problems using geometric multigrid

This paper investigates matrix-free algorithms for problems in quasi-static finite-strain hyperelasticity. Iterative solvers with matrix-free operator evaluation have emerged as an attractive alternative to sparse matrices in the fluid dynamics and wave propagation communities because they significantly reduce the memory traffic, the limiting factor in classical finite element solvers. Specifically, we study different matrix-free realizations of the finite element tangent operator and determine whether generalized methods of incorporating complex constitutive behavior might be feasible. Furthermore, in order to improve the convergence behavior of iterative solvers, we also propose a method by which to construct level tangent operators and employ them to define a geometric multigrid preconditioner. Additionally, the performance of the matrix-free operator and the geometric multigrid preconditioner is compared to the matrix-based implementation with an algebraic multigrid (AMG) preconditioner on a single node for a representative numerical example of a heterogeneous hyperelastic material in two and three dimensions. We find that matrix-free methods for finite-strain solid mechanics are very promising, outperforming linear matrix-based schemes by two to five times, and that it is possible to develop numerically efficient implementations that are independent of the hyperelastic constitutive law.

42 ENGINEERING↗
DOE OSTI2021

Prompt Gamma Analysis of ARIES Materials and Updates to the 2015 Calibration Equations

Prompt gamma (PG) analysis is a nondestructive, nuclear, elemental analysis technique that uses charged particle reactions to interrogate a sample, and elements present in the sample matrix are identified through the characteristic gamma-rays emitted from the product nuclei in alpha-p and alpha-n nuclear reactions. This technique has been applied to plutonium oxide packaged in over 4,000 individual 3013 containers, and the concentrations of certain light elements were determined from the integrated peak areas based on a calibration that was published previously. This report provides the results for a new population of 3013 containers packaged with oxide materials produced by the conversion of metal by Advanced Recovery and Integrated Extraction System (ARIES) project using the direct metal oxidation (DMO) process and muffle furnaces. New PG and analytical chemistry data collected since 2015 were added to the existing calibration data sets to refine the calibration parameters. Calibration equations were also developed for determining beryllium and fluorine at low concentrations in high-purity ARIES product oxides. The new fluorine calibration provides an order of magnitude greater sensitivity and results in an additional 1,408 containers in the original population identified as having fluorine as an impurity. Additionally, equations for calculating the lower limits of detection (LLDs) as a function of the actual counting time (live time) were obtained using WLS regression technique for samples in the calibration data set. This resulted in changes to the LLDs published previously.

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

Reduction of the molecular hamiltonian matrix using quantum community detection

Abstract Quantum chemistry is interested in calculating ground and excited states of molecular systems by solving the electronic Schrödinger equation. The exact numerical solution of this equation, frequently represented as an eigenvalue problem, remains unfeasible for most molecules and requires approximate methods. In this paper we introduce the use of Quantum Community Detection performed using the D-Wave quantum annealer to reduce the molecular Hamiltonian matrix in Slater determinant basis without chemical knowledge. Given a molecule represented by a matrix of Slater determinants, the connectivity between Slater determinants (as off-diagonal elements) is viewed as a graph adjacency matrix for determining multiple communities based on modularity maximization. A gauge metric based on perturbation theory is used to determine the lowest energy cluster. This cluster or sub-matrix of Slater determinants is used to calculate approximate ground state and excited state energies within chemical accuracy. The details of this method are described along with demonstrating its performance across multiple molecules of interest and bond dissociation cases. These examples provide proof-of-principle results for approximate solution of the electronic structure problem using quantum computing. This approach is general and shows potential to reduce the computational complexity of post-Hartree–Fock methods as future advances in quantum hardware become available.

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

Majorana versus Dirac Constraints on the Neutrino Dipole Moments

Massive neutrinos are guaranteed to have nonzero electromagnetic moments and, since there are at least three neutrino species, these dipole moments define a matrix. Here, we estimate the current upper bounds on all independent neutrino electromagnetic moments, concentrating on Earth-bound experiments and measurements with solar neutrinos, including the very recent results reported by XENONnT. We make no simplifying assumptions and compare the hypotheses that neutrinos are Majorana fermions or Dirac fermions. In particular, we fully explore constraints in the Dirac-neutrino parameter space. Majorana and Dirac neutrinos are different; for example, the upper bounds on the magnitudes of the elements of the dipole moment matrix are weaker for Dirac neutrinos, relative to Majorana neutrinos. The potential physics reach of next-generation experiments also depends on the nature of the neutrino. We find that a next-generation experiment two orders of magnitude more sensitive to the neutrino electromagnetic moments via $\nu_{\mu}$ elastic scattering may discover that the neutrino electromagnetic moments are nonzero if the neutrinos are Dirac fermions. Instead, if the neutrinos are Majorana fermions, such a discovery is ruled out by existing solar neutrino data, unless there are more than three light neutrinos.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Microwave Cavity Simulation Using Ansys HFSS

The design of microwave cavity detectors for axion dark matter research is often accomplished using advanced full-wave electromagnetic simulation software tools. These tools provide a cost-effective approach to evaluate a wide variety of cavity configurations, frequency tuning mechanisms, and conductive or dielectric materials and coatings. One simulation software package used for this application is Ansys High Frequency Structure Simulator (HFSS), which is based upon the well-established finite element method. HFSS includes numerous features useful for microwave cavity design such as parametric geometry modeling, adaptive meshing algorithm, curvilinear mesh elements, driven modal and eigenmode matrix solvers, and optimization algorithms. This paper describes the use of the HFSS software to simulate microwave cavities for axion haloscope detectors, with an example tutorial for a cylindrical cavity. Excellent agreement between the simulated and analytical results is shown for the resonant frequency, quality factor, and form factor.

Ansys HFSS, Cavity simulation, finite element mode↗
DOE OSTI2021

Method and apparatus for shape-based energy analysis of solids

A computer-readable medium stores instructions including storing a model of a physical structure and defining a mesh for the model. The mesh includes finite elements, each defined by a respective set of edges. The instructions include, for each finite element, identifying a governing differential equation and a set of complementary functions that exactly satisfy the governing differential equation. The instructions include determining an applied physical stimulus for the physical structure. The instructions include generating an energy optimization model that minimizes a difference between internal energy and external energy of the finite elements in response to the applied physical stimulus. The model includes a matrix of respective scalar multipliers for the complementary functions for each finite element. The instructions include transforming the matrix and calculating a physical parameter of interest. The instructions include, in response to the physical parameter not satisfying a design parameter, updating the model of the physical structure.

Spears, Robert E.↗
DOE OSTI2026

Hourglass control in staggered-grid hydrodynamics using virtual element stabilization techniques

Numerical simulations using the staggered-grid hydrodynamics (SGH) discretization suffer from hourglass instabilities. In this work, we develop a stabilization method to suppress the hourglass instabilities using techniques from the virtual element method (VEM). The stiffness matrix of the VEM consists of two terms: the consistency matrix which is rank deficient and the stability matrix. Here, we first show that in two dimensions and on general polygons, the stiffness matrix of the SGH is identical to the consistency matrix of the linear VEM for both the diffusion equation and the linear elasticity equation. These analyses explain the origin of the hourglass instabilities of the SGH discretization method, and establish a theoretical foundation for our proposed stabilization method by augmenting the stiffness matrix of the SGH discretization using the VEM stability matrix. Then, we present numerical examples using Lagrangian SGH simulations. The numerical experiments demonstrate that the proposed VEM stabilization method is effective at eliminating hourglass modes in the SGH discretization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Simulation of coupled multiphase flow and geomechanics in porous media with embedded discrete fractures

In fractured natural formations, the equations governing fluid flow and geomechanics are strongly coupled. Hydrodynamical properties depend on the mechanical configuration, and they are therefore difficult to accurately resolve using uncoupled methods. In recent years, significant research has focused on discretization strategies for these coupled systems, particularly in the presence of complicated fracture network geometries. In this work, we explore a finite-volume discretization for the multiphase flow equations coupled with a finite-element scheme for the mechanical equations. Fractures are treated as lower dimensional surfaces embedded in a background grid. Interactions are captured using the embedded discrete fracture model (EDFM) and the embedded finite element method (EFEM) for the flow and the mechanics, respectively. This nonconforming approach significantly alleviates meshing challenges. EDFM considers fractures as lower dimension finite volumes that exchange fluxes with the rock matrix cells. The EFEM method provides, instead, a local enrichment of the finite-element space inside each matrix cell cut by a fracture element. Both the use of piecewise constant and piecewise linear enrichments are investigated. They are also compared to an extended finite element approach. One key advantage of EFEM is the element-based nature of the enrichment, which reduces the geometric complexity of the implementation and leads to linear systems with advantageous properties. Synthetic numerical tests are presented to study the convergence and accuracy of the proposed method. It is also applied to a realistic scenario, involving a heterogeneous reservoir with a complex fracture distribution, to demonstrate its relevance for field applications.

58 GEOSCIENCES↗
DOE OSTI2023

FPGA Acceleration of GCN in Light of the Symmetry of Graph Adjacency Matrix

Graph Convolutional Neural Networks (GCNs) are widely used to process large-scale graph data. Different from deep neural networks (DNNs), GCNs are sparse, irregular, and unstructured, posing unique challenges to hardware acceleration with regular processing elements (PEs). In particular, the adjacency matrix of a GCN is extremely sparse, leading to frequent but irregular memory access, low spatial/temporal data locality and poor data reuse. Furthermore, a realistic graph usually consists of unstructured data (e.g., unbalanced distributions), creating significantly different processing times and imbalanced workload for each node in GCN acceleration. To overcome these challenges, we propose an end-to-end hardware-software co-design to accelerate GCNs on resource-constrained FPGAs with the features including: (1) A custom dataflow that leverages symmetry along the diagonal of the adjacency matrix to accelerate feature aggregation for undirected graphs. We utilize either the upper or the lower triangular matrix of the adjacency matrix to perform aggregation in GCN to improve data reuse. (2) Unified compute cores for both aggregation and transform phases, with full support to the symmetry-based dataflow. These cores can be dynamically reconfigured to the systolic mode for transformation or as individual accumulators for aggregation in GCN processing. (3) Preprocessing of the graph in software to rearrange the edges and features to match the custom dataflow. This step improves the regularity in memory access and data reuse in the aggregation phase. Moreover, we quantize the GCN precision from FP32 to INT8 to reduce the memory footprint without losing the inference accuracy. We implement our accelerator design in Intel Stratix10 MX FPGA board with HBM2, and demonstrate 1.3x-110.5x improvement in end-to-end GCN latency as compared to the state-of the-art FPGA implementations, on the graph datasets of Cora, Pubmed, Citeseer and Reddit.

Nair, Gopikrishnan R.↗
DOE OSTI2020

Perspectives on Quenching and Tempering 4340 Steel

Steels are ubiquitous due to their affordability and the landscape of useful properties that can be generated for engineering applications. But to further expand the performance envelope, one must be able to understand and control microstructure development by alloying and processing. In this paper we use multiscale, advanced characterization to better understand the structural and chemical evolution of AISI 4340 steel after quenching and tempering (Q&T), including the role of quench rate and short-time, isothermal tempering below 573 K (300 °C), with an emphasis on carbide formation. We compare the microstructure and/or property changes produced by conventional tempering to those produced by higher temperature, short-time “rapid” tempering. We underscore that no single characterization technique can fully capture the subtle microstructure changes like carbon redistribution, transition carbide and/or cementite formation, and retained austenite decomposition that occur during Q&T. Only the use of multiple techniques begins to unravel these complexities. After controlled fast or slow quenching, η transition carbides clearly exist in the microstructure, likely associated with autotempering of this high martensite start temperature (Ms) steel. Isothermal tempering below 598 K (325 °C) results in the relief of carbon supersaturation in the martensite, primarily by the formation of η transition carbides that exhibit a range of carbon levels, seemingly without substitutional element partitioning between the carbide and matrix phases. Hägg transition carbide is present between 300 °C and 325 °C. After conventional tempering at or above 598 K (325 °C) for 2 h, cementite is predominant, but small amounts of cementite are also present in other conditions, even after quenching. Previous work has indicated that silicon (Si) and substitutional elements partition between the cementite, which initially forms under paraequilibrium conditions, and the matrix. Phosphorous (P) may also be preferentially located at cementite/matrix interfaces after high temperature tempering. Slower quench rates result in greater amounts of retained austenite compared to those after fast quenching, which we attribute to increased austenite stability resulting from “autopartitioning”. Rapid, high temperature tempering is also found to diminish tempered martensite embrittlement (TME) believed to be associated with the extent of austenite decomposition, resulting in mechanical properties not attainable by conventional tempering, which may have important implications with respect to industrial heat treatment processes like induction tempering. Controlling the amount and stability of retained austenite is not only relevant to the properties of Q&T steels, but also next-generation advanced high strength steels (AHSS) with austenite/martensite mixtures.

36 MATERIALS SCIENCE↗