SEARCH · Engineering Papers

Results for “matrix approximation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

DOE OSTI2022

The performance of additively manufactured Haynes 282 in supercritical CO 2

The use of supercritical carbon dioxide (sCO 2 ) as a working fluid is garnering interest in next generation power production systems due to the possibility of increased operational efficiencies and lower associated costs. Its implementation requires alloys with excellent high temperature strength and corrosion resistance, which makes Haynes 282® (H282) a suitable candidate. With increasing adoption of additive manufacturing (AM) in the energy industry, there is a need to investigate long-term sCO2 exposure on AM produced materials. Additively manufactured H282 (AM-H282) samples built using laser powder bed fusion in two different orientations were analyzed through imaging, gravimetric analysis, and room-temperature tensile testing before and after 1000-h exposure to CO 2 at 750 °C and 20 MPa. Performed imaging included observation of the sample surfaces and of the material bulk (cross sections) through means of Scanning Electron Microscopy (SEM) and optical microscopy. In comparison to wrought H282, it was found that the exposed material exhibited a similar Cr 2 O 3 oxide protective layer about 2 μm in thickness with additional top-surface TiO 2 oxide and carbon-rich precipitates. It was also similarly observed that internal oxidation was present but limited to a depth of 20 μm, and a 1–3 μm γ’ denuded region appeared to surround all internal and surface oxidation. Thermal aging effects were noted with the precipitation of needle-like structures in the γ/γ’ matrix and a coarsening of the γ’ precipitates from 28 nm to 73 nm. Mass measurements analyzing oxide precipitation revealed a larger increase compared to wrought, representing an approximate 20% difference. Tensile testing results showed similar behavior to wrought with a slight increase in yield strength and a large reduction in elongation in the exposed samples.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Neutron elastic scattering on calcium isotopes from chiral nuclear optical potentials

We formulate microscopic neutron-nucleus optical potentials from many-body perturbation theory based on chiral two- and three-body forces. The neutron self-energy is first calculated in homogeneous matter to second order in perturbation theory, which gives the central real and imaginary terms of the optical potential. Here, the real spin-orbit term is calculated separately from the density matrix expansion using the same chiral interaction as in the self-energy. Finally, the full neutron-nucleus optical potential is derived within the improved local density approximation utilizing mean-field models consistent with the chiral nuclear force employed. We compare the results of the microscopic calculations to phenomenological models and experimental data up to projectile energies of E = 200 MeV. Experimental elastic differential scattering cross sections and vector analyzing powers are generally well reproduced by the chiral optical potential, but we find that total cross sections are overestimated at high energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Magnus method for electronic structure calculations at extreme conditions

We present the application of Magnus based methods to the solution of first order coupled ordinary differential equations in High Energy Density (HED) physics applications. Our focus is on the application to quantum mechanical methods, specifically on the solution of the radial Dirac equation for real and complex energies. HED applications require accurate solutions across a wide range of spatial and energy domains, including regimes where the solutions exhibit pronounced oscillatory behavior. Such cases pose significant computational challenges. We demonstrate that Magnus-based integrators can efficiently and accurately address these challenges. We discuss the implementation of the Magnus method for the solution of the radial Dirac equation, including practical considerations such as the evaluation of matrix exponentials, numerical integration, error estimation, and adaptive step size control. We also discuss the application of these methods to complex energy Green’s function techniques and the efficient approximation of integrals of the solutions relevant to HED electronic structure calculations. Here, we demonstrate the accuracy and robustness of the resulting method in applications to the free-particle case, for which analytic solutions are available for comparison, as well as the challenging case of gold at HED conditions.

general physics↗

DOE OSTI2021

First-principles study of electron transport in ScN

We investigate the conduction-band structure and electron mobility in rocksalt ScN based on density functional theory. The first-principles band structure allows us to obtain band velocities and effective masses as a function of energy. Electron-phonon scattering is assessed by explicitly computing the q-dependent electron-phonon matrix elements, with the inclusion of the long-range electrostatic interaction. The influence of free-carrier screening on the electron transport is assessed using the random phase approximation. We find a notable enhancement of electron mobility when the carrier concentration exceeds 10 20 cm -3 . We calculate the room-temperature electron mobility in ScN to be 587 cm2 /Vs at low carrier concentrations. When the carrier concentration is increased, the electron mobility starts to decrease significantly around n = 10 19 cm -3 , and drops to 240 cm 2 /Vs at n = 10 21 cm -3 . We also explore the influence of strain in (111)- and (100)-oriented ScN films. For (111) films, we find that a 1.0% compressive epitaxial strain increases the in-plane mobility by 72 cm 2 /Vs and the out-of-plane mobility by 50 cm 2 /Vs. For (100) films, a 1.0% compressive epitaxial strain increases the out-of-plane mobility by as much as 172 cm2 /Vs, but has a weak impact on the in-plane mobility. Our study sheds light on electron transport in ScN at different electron concentrations and shows how strain engineering could increase the electron mobility.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Reduced-Order Models of Static Power Grids based on Spectral Clustering

For large-scale interconnected power systems that cover large geographical areas, certain electrical studies are required so that appropriate decisions ensure system reliability and low cost. For such studies, it is often neither practical nor necessary to model in detail the entire power system, which is increasingly complex due to a more diverse range of grid assets to choose from in both short and long-term planning. The goal of this paper is to present a methodology to reduce the order of large-scale power networks based on spectral graph theory given that current methods for static network reduction are not scalable. A brief analysis of some spectral clustering properties to determine which graph Laplacian matrix should be used and why is included. The analysis shows that the utilization of the normalized graph Laplacian is more advantageous for clustering purposes. Techniques are proposed to approximate cost functions for the aggregated generators. This is done via linear regression. The reduced-order model obtained with the proposed methodology has an accuracy above 94% and solves the scalability issue commonly present in other reduction methods. If the utilization of the reduced-order model is either constrained to load levels above mid-peak demand, or cost functions of aggregated units are approximated via a piecewise quadratic approach, then the error distribution is in the order of 10^-3. .

Baquedano-Aguilar, Mario D.↗

DOE OSTI2024

Development of Segregated Thermal-Hydraulics Solvers in MOOSE

The simulation of fluid flows is an essential part of the design and analysis of nuclear systems. Algorithms able to simulate flows at different fidelity levels are available in the Multiphysics Object-Oriented Simulation Environment (MOOSE) and MOOSE-based applications such as Pronghorn \cite{novak2018pronghorn}, Pronghorn-Subchannel, RELAP-7, and SAM. Currently, significant effort is being invested in the development of coarse-mesh Computational Fluid Dynamics (CFD) capabilities within MOOSE and Pronghorn for the simulation of Generation IV nuclear reactors. Traditionally, the solution algorithms in MOOSE have relied on Newton or quasi-Newton methods (such as the preconditioned Jacobian-free Newton-Krylov method) where residuals and Jacobians (or approximations thereof) are constructed. Both Newton and quasi-Newton methods require the solution of a linear system at each nonlinear Newton iteration with the Jacobian as the system matrix. The Jacobian contains blocks originating from all variables in the problem (i.e., for thermal-hydraulics at least pressure, velocities, and temperature). Due to the formulation of the problem in a general multiphysics setting on unstructured mesh, creating a good preconditioner for the linear system can be challenging, thus many fluid applications have utilized direct solver-based methods such as LU factorization. However, with increasing system size and complexity in multi-dimensional problems, the direct solution of linear systems becomes computationally expensive both in execution time and and memory. For this reason, recent effort has focused on adapting segregated solution algorithms for CFD problems in MOOSE. These algorithms use fixed-point iteration between segregated systems whose assembly and preconditioning are easier those of the monolithic system. Initial results show that the segregated solution algorithm outperforms the monolithic approach in terms of memory usage and for large 3D problems in terms of CPU time as well.

42 ENGINEERING↗

DOE OSTI2024

Development of Segregated Thermal-Hydraulics Solvers in MOOSE

42 ENGINEERING↗

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 OSTI2024

Investigating sol–gel matrix loading capacity toward producing surrogate nuclear explosive debris with realistic composition

Post detonation nuclear forensic materials which resemble the size, color, elemental composition, and radionuclide content of real nuclear debris would be valuable for developing and validating new nuclear forensic techniques. As nuclear fallout types vary significantly, the ability to tailor each of these parameters accurately is desired to produce materials capable of testing analytical methods under a wide array of forensic scenarios. Sol–gel synthesis techniques can provide tunability of size, shape and composition for producing a wide variety of solid nuclear forensics benchmarking materials. Further, the sol–gel process consists of forming a metal oxide material, often silica, through polymerization of a metal-alkoxy precursor. In this work, we characterize the ability to load sol–gel particles with secondary elemental components such as iron, aluminum, and calcium toward producing benchmarking materials approximating the elemental composition of historic nuclear debris from the Nevada National Security Site. We also demonstrate quantitative radionuclide encapsulation toward producing benchmarking materials with controllable radionuclide content. Finally, we employ these techniques to produce nuclear debris benchmarking materials with controllable elemental matrix composition and radionuclide content and compare these samples with the composition of a historic fallout sample previously reported from the Nevada National Security Site.

36 MATERIALS SCIENCE↗

DOE OSTI2024

ZeroCAL: Eliminating Carbon Dioxide Emissions from Limestone’s Decomposition to Decarbonize Cement Production

Limestone (calcite, CaCO 3 ) is an abundant and cost-effective source of calcium oxide (CaO) for cement and lime production. However, the thermochemical decomposition of limestone (~800 °C, 1 bar) to produce lime (CaO) results in substantial carbon dioxide (CO 2(g) ) emissions and energy use, i.e., ~1 tonne [t] of CO 2 and ~1.4 MWh per t of CaO produced. Here, we describe a new pathway to use CaCO 3 as a Ca source to make hydrated lime (portlandite, Ca(OH) 2 ) at ambient conditions (p, T) while nearly eliminating process CO 2(g) emissions (as low as 1.5 mol. % of the CO 2 in the precursor CaCO 3 , equivalent to 9 kg of CO 2(g) per t of Ca(OH) 2 ) within an aqueous flowelectrolysis/ pH-swing process that coproduces hydrogen (H 2(g) ) and oxygen (O 2(g) ). Because Ca(OH) 2 is a zero-carbon precursor for cement and lime production, this approach represents a significant advancement in the production of zero-carbon cement. The Zero CArbon Lime (ZeroCAL) process includes dissolution, separation/recovery, and electrolysis stages according to the following steps: (Step 1) chelator (e.g., ethylenediaminetetraacetic acid, EDTA)-promoted dissolution of CaCO 3 and complexation of Ca 2+ under basic (>pH 9) conditions, (Step 2a) Ca enrichment and separation using nanofiltration (NF), which allows separation of the Ca-EDTA complex from the accompanying bicarbonate (HCO 3 – ) species, (Step 2b) acidity-promoted decomplexation of Ca from EDTA, which allows near-complete chelator recovery and the formation of a Ca-enriched stream, and (Step 3) rapid precipitation of Ca(OH) 2 from the Ca-enriched stream using electrolytically produced alkalinity. These reactions can be conducted in a seawater matrix yielding coproducts including hydrochloric acid (HCl) and sodium bicarbonate (NaHCO 3 ), resulting from electrolysis and limestone dissolution, respectively. Careful analysis of the reaction stoichiometries and energy balances indicates that approximately 1.35 t of CaCO 3 , 1.09 t of water, 0.79 t of sodium chloride (NaCl), and ~2 MWh of electrical energy are required to produce 1 t of Ca(OH) 2 , with significant opportunity for process intensification. This approach has major implications for decarbonizing cement production within a paradigm that emphasizes the use of existing cement plants and electrification of industrial operations, while also creating approaches for alkalinity production that enable cost-effective and scalable CO 2 mineralization via Ca(OH) 2 carbonation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Sub-system self-consistency in coupled cluster theory

Here, in this article, we provide numerical evidence indicating that the single-reference coupled-cluster (CC) energies can be calculated alternatively to their copybook definition. We demonstrate that the CC energy can be reconstructed by diagonalizing the effective Hamiltonians describing correlated sub-systems of the many-body system. In the extreme case, we provide numerical evidence that the CC energy can be reproduced through the diagonalization of the effective Hamiltonian describing sub-system composed of a single electron. These properties of the CC formalism can be exploited to design protocols to define effective interactions in sub-systems used as probes to calculate the energy of the entire system and introduce a new type of self-consistency for approximate CC approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Quantum dynamics in strongly driven random dipolar magnets

The random dipolar magnet LiHo x Y 1– x F 4 enters a strongly frustrated regime for small Ho 3+ concentrations with x < 0.05. In this regime, the magnetic moments of the Ho 3+ ions experience small quantum corrections to the common Ising approximation of LiHo x Y 1– x F 4 , which lead to a $\mathbb{Z}_2$-symmetry breaking and small, degeneracy breaking energy shifts between different eigenstates. Here we show that destructive interference between two almost degenerate excitation pathways burns spectral holes in the magnetic susceptibility of strongly driven magnetic moments in LiHo x Y 1– x F 4 . Such spectral holes in the susceptibility, microscopically described in terms of Fano resonances, can already occur in setups of only two or three frustrated moments, for which the driven level scheme has the paradigmatic Λ shape. For larger clusters of magnetic moments, the corresponding level schemes separate into almost isolated many-body Λ schemes, in the sense that either the transition matrix elements between them are negligibly small or the energy difference of the transitions is strongly off-resonant to the drive. This enables the observation of Fano resonances, caused by many-body quantum corrections to the common Ising approximation also in the thermodynamic limit. Finally, we discuss its dependence on the driving strength and frequency as well as the crucial role that is played by lattice dissipation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Introduction of the Adding and Doubling Method for Solving Bateman Equations for Nuclear Fuel Depletion

This paper introduces and evaluates the Adding and Doubling Method (ADM) for solving the Bateman equations for depletion systems with varying numbers of nuclides and compares it to the Chebyshev Rational Approximation Method (CRAM), both implemented in the reactor physics analysis application Griffin. ADM, when applied to the Crank-Nicolson Finite Difference method, can produce results comparable in accuracy and precision to CRAM with comparable run times for systems with 35 or 297 nuclides. For systems with more than 300 nuclides, the matrix-matrix operations required by ADM are significantly more costly than the matrix-vector operations required by CRAM, making CRAM the more efficient method for systems with large numbers of nuclides. ADM is an accurate method that maintains other advantages over CRAM in that it does not depend on pre-generated coefficients or require complex number operations. ADM also manages to outperform CRAM by a factor of more than 250 in terms of run time for depletion systems that require multiple Bateman solves while the depletion matrix and time step size remain constant over all depletion intervals.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

UN-SiC TRISO Post Irradiation Examination Developmental Work

As part of efforts to strengthen INL?s post irradiation analysis capabilities with non-Advanced Gas Reactor (AGR) tristructural isotropic (TRISO) fuels, two developmental activities were conducted. The first activity was to determine how best to analyze uranium nitride TRISO fuel kernels using electron probe microanalysis, while the second activity focused on developing a method to deconsolidate TRISO fuel particles that have been encased in a silicon carbide matrix. Because these two activities were unrelated, they have been presented separately in this report. Initial EPMA analyses showed nitrogen contents that far exceeded the concentration expected for UN--a line compound. Further examination showed that current literature values for the mass absorption coefficient (MAC) for the N ka X-ray absorbed by U ranged from approximately 1600 to 9500, with most values tending toward 9500. Measuring UN with five different progressively increasing accelerating voltages followed by using the modeling program xMAC suggests the actual MAC is approximately 2115. Additional MAC modifications were required to produce reasonable analytical results. Because of the inaccuracies of necessary MAC coefficients, UN analysis via scanning electron microscopy (SEM) is likely to produce inaccurate results. This is because SEM software does not typically allow the user to alter MACs. Tests have been performed to examine the feasibility of an electrochemical technique to liberate irradiated TRISO fuel from a SiC matrix without damaging the outer pyrolytic carbon layer of the fuel particle. The method is performed by electrochemically exposing the SiC to magnesium metal forming Mg2Si and C. Following exposure, the small SiC samples showed slight mass increases with no evidence of conversion to Mg2Si and C nor obvious degradation of the SiC samples.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Substrate effects on spin relaxation in two-dimensional Dirac materials with strong spin-orbit coupling

Understanding substrate effects on spin dynamics and relaxation is of key importance for spin-based information technologies. However, the key factors that determine such effects, in particular for materials with strong spin-orbit coupling (SOC), have not been well understood. Here we performed first-principles real-time density-matrix dynamics simulations with SOC and the electron phonon and electron-impurity scattering for spin lifetimes (τ s ) of supported/free-standing germanene, a prototypical strong SOC 2D Dirac material. We show that the effects of different substrates on τ s can surprisingly differ by two orders of magnitude. We find that substrate effects on τ s are closely related to substrate-induced modifications of the SOC-field anisotropy, which changes the spin flip scattering matrix elements. We propose a new electronic quantity, named spin-flip angle θ ↑↓ , to characterize spin relaxation through intervalley spin-flip scattering. We find that $τ^{-1}_s$ is approximately proportional to the averaged value of sin 2 (θ ↑↓ /2), which serves as a guiding parameter of controlling spin relaxation.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Sparse Approximate Multifrontal Factorization with Composite Compression Methods

This article presents a fast and approximate multifrontal solver for large sparse linear systems. In a recent work by Liu et al., we showed the efficiency of a multifrontal solver leveraging the butterfly algorithm and its hierarchical matrix extension, HODBF (hierarchical off-diagonal butterfly) compression to compress large frontal matrices. The resulting multifrontal solver can attain quasi-linear computation and memory complexity when applied to sparse linear systems arising from spatial discretization of high-frequency wave equations. To further reduce the overall number of operations and especially the factorization memory usage to scale to larger problem sizes, in this article we develop a composite multifrontal solver that employs the HODBF format for large-sized fronts, a reduced-memory version of the nonhierarchical block low-rank format for medium-sized fronts, and a lossy compression format for small-sized fronts. This allows us to solve sparse linear systems of dimension up to 2.7 × larger than before and leads to a memory consumption that is reduced by 70% while ensuring the same execution time. The code is made publicly available in GitHub.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

MAGMA: Enabling exascale performance with accelerated BLAS and LAPACK for diverse GPU architectures

MAGMA (Matrix Algebra for GPU and Multicore Architectures) is a pivotal open-source library in the landscape of GPU-enabled dense and sparse linear algebra computations. With a repertoire of approximately 750 numerical routines across four precisions, MAGMA is deeply ingrained in the DOE software stack, playing a crucial role in high-performance computing. Notable projects such as ExaConstit, HiOP, MARBL, and STRUMPACK, among others, directly harness the capabilities of MAGMA. In addition, the MAGMA development team has been acknowledged multiple times for contributing to the vendors’ numerical software stacks. Looking back over the time of the Exascale Computing Project (ECP), we highlight how MAGMA has adapted to recent changes in modern HPC systems, especially the growing gap between CPU and GPU compute capabilities, as well as the introduction of low precision arithmetic in modern GPUs. We also describe MAGMA’s direct impact on several ECP projects. Maintaining portable performance across NVIDIA and AMD GPUs, and with current efforts toward supporting Intel GPUs, MAGMA ensures its adaptability and relevance in the ever-evolving landscape of GPU architectures.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Using Homogenized Models to Explore the Effect of Fracture Densities on Weathering

Despite its importance, only a few researchers have incorporated the effects of fracturing into models of reactive transport for rock weathering. Here we explore 2D simulations that describe weathering under conditions of diffusive and advective transport within heterogeneous media consisting of rocky blocks and fractures. In our simulations, the Darcy velocities vary in space and time and depend on weathering processes within the rock matrix. We explore simulations with saturated and unsaturated flow for weathering bedrock that consists of blocks separated by inert or weathered material. The simulations show that a simplified homogenized model can approximate exact solutions for some of the simulated columns and hills and can allow exploration of coupling between flow and reaction in fractured rock. These hillslope simulations document that, even in the presence of 2D water flow, i) an increase in fracture density results in faster weathering advance rates; and ii) the water table locates deeper for a rock system that is weathered and fractured rather than weathered and unfractured. Some of these patterns have also been observed for natural systems. But these simulations also highlight how simplified models that do not use appropriate averaging of heterogeneities can be inaccurate in predicting weathering rate for natural systems. For example, if water flows both vertically and laterally through the vadose zone of a hill, then a prediction of the depth of regolith that is based on modeling strictly unidirectional downward infiltration will be unrealistically large. Likewise, if the fracture density observed near the land surface is used in a model to predict depth of weathering for a system where the fracture density decreases downward, the model will overestimate regolith depth. Learning how to develop accurately homogenized models could thus enable better conceptual models and predictions of weathering advance in natural systems.

58 GEOSCIENCES↗