Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “LU factorization”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
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

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 OSTI2023

Combining Sparse Approximate Factorizations with Mixed-precision Iterative Refinement

The standard LU factorization-based solution process for linear systems can be enhanced in speed or accuracy by employing mixed-precision iterative refinement. Most recent work has focused on dense systems. We investigate the potential of mixed-precision iterative refinement to enhance methods for sparse systems based on approximate sparse factorizations. In doing so, we first develop a new error analysis for LU- and GMRES-based iterative refinement under a general model of LU factorization that accounts for the approximation methods typically used by modern sparse solvers, such as low-rank approximations or relaxed pivoting strategies. We then provide a detailed performance analysis of both the execution time and memory consumption of different algorithms, based on a selected set of iterative refinement variants and approximate sparse factorizations. Our performance study uses the multifrontal solver MUMPS, which can exploit block low-rank factorization and static pivoting. We evaluate the performance of the algorithms on large, sparse problems coming from a variety of real-life and industrial applications showing that mixed-precision iterative refinement combined with approximate sparse factorization can lead to considerable reductions of both the time and memory consumption.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Multifrontal Non-negative Matrix Factorization

Non-negative matrix factorization (Nmf) is an important tool in high-performance large scale data analytics with applications ranging from community detection, recommender system, feature detection and linear and non-linear unmixing. While traditional Nmf works well when the data set is relatively dense, however, it may not extract sufficient structure when the data is extremely sparse. Specifically, traditional Nmf fails to exploit the structured sparsity of the large and sparse data sets resulting in dense factors. We propose a new algorithm for performing Nmf on sparse data that we call multifrontal Nmf (Mf-Nmf) since it borrows several ideas from the multifrontal method for unconstrained factorization (e.g. LU and QR). We also present an efficient shared memory parallel implementation of Mf-Nmf and discuss its performance and scalability. We conduct several experiments on synthetic and realworld datasets and demonstrate the usefulness of the algorithm by comparing it against standard baselines. We obtain a speedup of 1.2x to 19.5x on 24 cores with an average speed up of 10.3x across all the real world datasets.

Sao, Piyush↗
DOE OSTI2025

A two-level GPU-accelerated incomplete LU preconditioner for general sparse linear systems

This paper presents a parallel preconditioning approach based on incomplete LU (ILU) factorizations in the framework of Domain Decomposition (DD) for general sparse linear systems. We focus on distributed memory parallel architectures, specifically, those that are equipped with graphic processing units (GPUs). In addition to block-Jacobi, we present general purpose two-level ILU Schur complement-based approaches, where different strategies are presented to solve the coarse-level reduced system. These strategies are combined with modified ILU methods in the construction of the coarse-level operator, in order to effectively remove smooth errors by targeting an algebraically smooth vector. We leverage available GPU-based sparse matrix kernels to accelerate the setup and the solve phases of the proposed ILU preconditioner. We evaluate the efficiency of the proposed methods as a smoother for algebraic multigrid (AMG) and as a preconditioner for Krylov subspace methods on challenging anisotropic diffusion problems and a collection of general sparse matrices.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

High-throughput characterization of Lu-doped zirconia

Lu-doped zirconia (LuDZ) was reported in the 1960s to be a poor oxygen ion conductor. However, we identified this system as worthy of reexamination based on the results of a machine learning-based screen for promising oxygen ion conductors, further supported by a physical argument: because Lu is the heaviest and smallest of the lanthanides, it should offer the possibility of high conductivity due to its high polarizability and size compatibility with the zirconia structure. Here the relationship between dopant concentration and conductivity in the LuDZ system was determined using a high throughput experimental methodology in which spatially resolved conductivity was measured across a compositionally graded film, spanning from Zr 0.93 Lu 0.07 O 1.965 to Zr 0.82 Lu 0.18 O 1.91 . The conductivity broadly peaks with composition at 15 cation % Lu, coinciding with the approximate phase boundary of the fully stabilized cubic phase. The composition-dependent conductivity measurements reveal a surprising increase in activation energy for ion transport with increasing Lu content with no change in this trend at the composition of the conductivity maximum. The pre-exponential factor increases sharply with Lu content, but displays a distinct plateau once the conductivity maximum composition is attained. At 500 °C, the conductivity of Zr 0.85 Lu 0.15 O 1.925 is 2 × 10 –3 S cm –1 , approaching that of optimally doped Sc-stabilized zirconia, which has the highest reported conductivity among zirconia compounds. In comparison to other oxide ion conductors, Lu-stabilized zirconia offers the advantage over ceria of negligible electronic conductivity and is easier to process than (La,Sr)(Ga,Mg)O 3-δ .

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

Achieving high-ZT in Te/Se-based thermoelectric materials via volume-constrained shear-induced plastic deformation

This study investigated the correlation between severe plastic deformation and microstructural evolution, along with the corresponding response in the thermoelectric performance. Specifically, severe plastic deformation introduced into intrinsically brittle thermoelectric materials like Bi 2 Te 3 - x Se x leads to microstructural changes ranging from nano- to micro- scales, including phase transformation, phase segregation, and stacking faults. The resultant residual strain and lattice defects in the Bi 2 Te 3-x Se x materials led to engineering of electronic bands, which increases the effective mass of electrons while controlling their concentration, thereby enhancing the power factor.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Measurements and computational analysis of the natural decay of 176 Lu

Background: Mainly because of its long half-life and despite its scientific relevance, spectroscopic measurements of 176 Lu forbidden β decays are very limited and lack formulation of shape factors. A direct precise measurement of its Q value is also presently unreported. In addition, the description of forbidden decays provides interesting challenges for nuclear theory. The comparison of precise experimental results with theoretical calculations for these decays can help to test underlying models and can aid the interpretation of data from other experiments. Purpose: Perform the first precision measurements of 176 Lu β-decay spectra and attempt the observation of its electron capture decays, as well as perform the first precision direct measurement of the 176 Lu β-decay Q value. Compare the shape of the precisely determined experimental β spectra to theoretical calculations, and compare the end point energy to that obtained from an independent Q value measurement. Method: Here, the 176 Lu β-decay spectra measurements and the search for electron capture decays were performed with an experimental setup that employed lutetium-containing scintillator crystals and a NaI(Tl) spectrometer for coincidence counting. The β decay Q value was determined via high-precision Penning trap mass spectrometry (PTMS) with the LEBIT facility at the National Superconducting Cyclotron Laboratory. The β-spectrum calculations were performed within the Fermi theory formalism with nuclear structure effects calculated using a shell model approach. Results: Both β transitions of 176 Lu were experimentally observed and corresponding shape factors formulated in their entire energy ranges. The search for electron capture decay branches led to an experimental upper limit of 6.3×10 –6 relative to its β decays. The 176 Lu β-decay and electron capture Q values were measured using PTMS to be 1193.0(6) and 108.9(8) keV, respectively. This enabled precise β end point energies of 596.2(6) and 195.3(6) keV to be determined for the primary and secondary β decays, respectively. The conserved vector current hypothesis was applied to calculate the relativistic vector matrix elements. The β-spectrum shape was shown to significantly depend on the Coulomb displacement energy and on the value of the axial vector coupling constant gA, which was extracted according to different assumptions. Conclusion: The implemented self-scintillation method has provided unmatched observations of 176 Lu, independently validated by the first direct measurements of its β-decay Q value by Penning trap mass spectrometry. Theoretical study of the main β transition led to the extraction of very different effective gA and log10f values, showing that a high-precision description of this transition would require a realistic nuclear structure with nucleus deformation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Understanding the Adsorption of Rare-Earth Elements in Oligo-Grafted Mesoporous Carbon

Rare-earth elements (REEs) are 17 elements of the periodic table primarily consisting of lanthanides. In modern society, the usage of REEs is ubiquitous in almost all modern gadgets and therefore efficient recovery and separation of REEs are of high importance. Selective adsorption and chelation of REEs in solid sorbents is a unique and sustainable process for their recovery. In this work, single-stranded oligos with 100 units of thymine were grafted onto carboxylated mesoporous carbon to synthesize a sorbent with phosphorus and oxygen functionalities. Additionally, the sorbent was characterized by X-ray photoelectron spectroscopy, Fourier-transform infrared spectroscopy, and scanning electron microscopy-energy-dispersive X-ray spectroscopy. Three different REEs with varying atomic radii and densities, Lu, Dy, and La, were adsorbed onto the carbon from aqueous solutions. It was observed that the adsorbed amounts increased with the increase in the atomic radius or decrease in the atomic density. Calculation of the distribution coefficients for all the equilibrium adsorption amounts suggested that adsorption is more effective in the lower concentration region. The L 3 -edge X-ray absorption near-edge structure confirmed a 3+ oxidation state of REEs in the adsorbed phase. Extended X-ray absorption fine structure (EXAFS) confirmed the binding of REEs with oxygen functionalities in the adsorbed phase. The radial distribution functions calculated from the EXAFS data suggest a longer RE–O distance for La compared to those for Lu and Dy. The coordination numbers and Debye–Waller factors have typical values of about 8–9 atoms and 0.01–0.02 Å 2 , respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Structural Complexity and Tuned Thermoelectric Properties of a Polymorph of the Zintl Phase Ca 2 CdSb 2 with a Non-centrosymmetric Monoclinic Structure

The Zintl phase Ca 2 CdSb 2 was found to be dimorphic. Besides the orthorhombic Ca 2 CdSb 2 (-o), here we report on the synthesis, the structural characterization, and the thermoelectric transport properties of its monoclinic form, Ca 2 CdSb 2 (-m), and its Lu-doped variant Ca 2–x Lu x CdSb 2 (x ≈ 0.02). The monoclinic structure exhibits complex structural characteristics and constitutes a new structure type with the non-centrosymmetric space group Cm (Z = 30). The electrical resistivity ρ(T) measured on single crystals of both phases portrays a transition from a semiconductor to a degenerate p-type semiconductor upon doping with Lu and with an attendant change in the Hall carrier concentration nH from 7.15 × 10 18 to 2.30 × 10 19 cm –3 at 300 K. The Seebeck coefficient S(T) of both phases are comparable and indicate a hole-dominated carrier transport mechanism with magnitudes of 133 and 116 μV/K at 600 K for Ca 2 CdSb 2 (-m) and Ca 2–x Lu x CdSb 2 , respectively. The convoluted atomic bonding with an attendant large unit cell volume of ~4365 Å 3 drives a putative low thermal conductivity in these materials resulting in a power factor PF of 1.63 μW/cm K 2 and an estimated thermoelectric figure of merit zT of ~0.5 for Ca 2–x Lu x CdSb 2 at 600 K. Furthermore, differential scanning calorimetry results reveal the stability of these phases up to about 960 K, making them candidates for moderate temperature thermoelectric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Connecting Geomechanical Properties with Potential for Proppant Embedment and Production Decline for the Emerging Caney Shale, Oklahoma

The Caney Shale is emerging as a target for hydrocarbon production, creating opportunities to study rock mechanical origins of observed/anticipated challenges to effective and sustained stimulation. Here we examine five subunits within the Caney, two of which are dubbed “ductile” and three of which are dubbed “reservoir” rock types. The “ductile” versus “reservoir” identification is initially based on elastic properties ascertained from well logs. By comparing core-based mechanical tests, it is found nominally ductile and reservoir layers do not differ in terms of brittleness, but instead the nominally ductile layers tend to be weaker and more prone to creep deformation compared to the nominally more brittle reservoir layers. This difference in mechanical properties is shown to generate higher susceptibility to proppant embedment in the weaker and more creep-prone layers. Specifically, over a 5-year period, the creep is expected to have little impact on reservoir layers but is predicted to reduce propped fracture aperture by a factor of 2 in the ductile layers.

Benge, Margaret↗
DOE OSTI2024

Energy Analysis of Combi Heat Pump System Configurations for Space Conditioning and Domestic Hot Water Heating in Residential Buildings

Combi heat pump systems, also referred to multifunctional variable refrigerant flow heat recovery (MF-VRFHR) systems, are specifically designed for residential applications to manage both space conditioning and domestic hot water (DHW). They have attracted attention due to their potential for energy conservation through heat recovery. The incorporation of a hot water tank introduces various system configurations, each characterized by distinct pros and cons related to energy efficiency, system stability, and maintenance. Despite this, a critical gap exists as the specific energy performance remains unquantified under diverse operational modes (e.g., heating mode and heat recovery mode). This paper aims to bridge this gap by conducting a comprehensive comparative analysis of two prevalent system configurations while considering feasible proposed control logics. Configuration 1 integrates a separate hot water tank and a refrigerant-to-water heat exchanger (HEX), also known as a Hydro Kit while Configuration 2 incorporates a refrigerant-wrapped hot water tank. To facilitate this analysis, we developed high-fidelity system models for both configurations in Modelica, capturing system dynamics and detailed control sequences effectively. These system models were built upon the TIL library for HVAC equipment components and the Buildings library for residential building thermal load calculations. The validation of the simulation testbed utilized data from experiments conducted in the PNNL lab home for Configuration 1. To establish the simulation testbed for Configuration 2, we extended the modeling setup derived from Configuration 1. This extension specifically involved substituting the separate hot water tank and Hydro Kit with a refrigerant-wrapped hot water tank of similar sizing sourced from an actual product. The simulation analysis of heating-only and heat recovery modes reveals that Configuration 2 not only saves energy and maintains warmer tank temperatures but also demonstrates faster water heating capabilities. This is attributed to decreased energy loss and improved heat transfer. The study encompasses a wide range of scenarios, considering diverse thermal loads and water usage patterns across heating and heat recovery modes. Overall, the comprehensive results indicate that Configuration 2 achieves energy savings ranging from 3.5% to 12.2% compared to Configuration 1, depending on factors such as water usage patterns, thermal loads, and operational modes.

Configuration, Comparison, Multi-functional, Resid↗
DOE OSTI2022

Comparative Small-Signal Stability Analysis of Grid-Forming and Grid-Following Inverters in Low-Inertia Power Systems: Preprint

Inverter-based resources (IBRs) are gradually replacing conventional fossil energy in modern power systems, which could make the system susceptible due to insufficient inertia, especially when moving towards a 100% inverter-based power grid. Grid-forming (GFM) and grid-following (GFL) controls are two common inverter control algorithms. However, the corresponding stability issues along with a comparative study on low-inertia power systems should still be investigated. This paper focuses on small-signal stability and conducts a comprehensive analysis of GFM and GFL controls, with emphases on the interactions between synchronous generators and IBRs. The detailed modeling of GFM and GFL inverters, and synchronous generators (SGs) are established. The coupling mechanism of IBRs and SGs is then conducted through participation factor analysis. The impacts of critical parameters on system stability are finally identified. The theoretical analyses are validated through time-domain simulation in MATLAB/Simulink.

grid-following↗