Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix 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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Phase-field modeling of thermally-grown oxide and damage evolution in environmental barrier coatings

Silicon carbide-based ceramic matrix composites protected by environmental barrier coatings (EBCs) present a promising materials solution for next-generation gas turbines. Developing more robust and efficient EBCs is therefore of significant technological importance. During the service in high-temperature oxidative environments, there is a thermally grown oxide (TGO) layer, spontaneously formed in the EBC system. TGO is recognized as a critical factor for the degradation and failure of EBCs, yet the detailed mechanisms of TGO growth and its effect on EBC failure remain unclear. In this study we develop a comprehensive chemo-mechano-phase-field model to simulate growth of the TGO in EBCs, factoring in creep and deformation, and especially the cracking behaviors. The volume expansion due to TGO growth and the resulting large inelastic deformation are addressed by using our recently developed, so-called incremental realization of inelastic deformation (IRID) algorithm, in combination with an adapted Hu-Chen spectral solver for elasticity. Simulations of TGO growth are performed considering different growth modes of TGOs determined mainly by the ratio of oxidant permeability in the topcoat to that in the TGO itself. Large-scale three-dimensional (3D) simulations are performed to model the formation of interconnecting vertical/channel cracks (often called ‘mud cracks’). The simulated crack morphology are in excellent agreement with the experimental observations from the literature. The simulations also provide insights into the cracking of EBCs and its dependence on the structure and constituent properties of the coating system. Furthermore, these results demonstrate the developed damage model can be a useful tool for design of more durable EBCs.

Cracking↗
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 OSTI2025

Enhancing cold spray coatings: Microstructural dynamics and performance attributes of Inconel 625 with chromium carbide incorporation for hydropower applications

The incorporation of chromium carbide (CrC) particles into the cold spray (CS) process is known to mitigate nozzle clogging, although at the expense of deposition efficiency. This study explores the intricate microstructural changes induced by varying amounts of CrC powder (12.5 % and 6 %) in conjunction with Inconel-625 (Inc-625) powder. The deposition was carried out onto A27 cast steel under different CS parameters. Microstructural characterization, including detailed electron microscopy studies, reveals a complex yet structurally stable coating. Noteworthy features include grain fragmentation and a cellular structure enriched with Nb and Mo, with minimal plastic deformation of CrC in the matrix. The cold-sprayed coatings exhibit a significant (~4 times) increase in microhardness compared to the A27 substrate. Mechanical and cavitation erosion properties were systematically investigated. Coatings subjected to higher particle energy conditions with a gas pressure of 600 psi and gas temperature of 650 °C, demonstrated superior resistance to cavitation erosion. This resistance is attributed to a combination of factors, including microstructural characteristics and porosity. Altogether, the study provides valuable insights into the structural dynamics and performance of CS coatings enriched with CrC particles.

A27 cast steel↗
DOE OSTI2024

Graph network heterogeneity predicts interplant wake losses

Wind plants generate large-scale wakes, which can affect the performance of neighboring installations. Such wakes are challenging to estimate due to the inherent complexity in modeling wake interactions between large quantities of turbines at various distances. Weighted directed graph networks can inform complex models by linking turbine pairs into chains of upstream and downstream neighbors for a given wind direction. A novel interpretation of the graph network adjacency matrix is proposed where each element of the matrix represents the cumulative impact of upstream turbines on an individual. In this study, wake losses were estimated with an engineering wake model across a range of inflow conditions for nine parametric variations of a system containing two neighboring wind plants. The parametric nature of the study isolates turbine spacing within the plant, separation distance between plants, and wind direction as the main drivers of wake losses. Spatial heterogeneity is computed from the weighted average adjacency matrix of each plant arrangement. The proposed method is orders of magnitude faster than wake modeling and does not require detailed turbine information or atmospheric conditions. Furthermore, the weighted average adjacency matrix provides insight on the spatial organization of wake losses at various scales. Plant heterogeneity is correlated with wake losses within and among plants. Framing wind plant wake interaction in terms of graph network spatial heterogeneity provides an efficient approach for predicting wake losses within and among neighboring wind plants with applications to other complex systems where wake interactions are key factor.

17 WIND ENERGY↗
DOE OSTI2025

An Adaptive Newton-Based Free-Boundary Grad–Shafranov Solver

Equilibria in magnetic confinement devices result from force balancing between the Lorentz force and the plasma pressure gradient. In an axisymmetric configuration like a tokamak, such an equilibrium is described by an elliptic equation for the poloidal magnetic flux, commonly known as the Grad–Shafranov equation. It is challenging to develop a scalable and accurate free-boundary Grad–Shafranov solver, since it is a fully nonlinear optimization problem that simultaneously solves for the magnetic field coil current outside the plasma to control the plasma shape. In this work, we develop a Newton-based free-boundary Grad–Shafranov solver using adaptive finite elements and preconditioning strategies. The free-boundary interaction leads to the evaluation of a domain-dependent nonlinear form of which its contribution to the Jacobian matrix is achieved through shape calculus. The optimization problem aims to minimize the distance between the plasma boundary and specified control points while satisfying two nontrivial constraints, which correspond to the nonlinear finite element discretization of the Grad–Shafranov equation and a constraint on the total plasma current involving a nonlocal coupling term. The linear system is solved by a block factorization, and AMG is called for subblock elliptic operators. The unique contributions of this work include the treatment of a global constraint, preconditioning strategies, nonlocal reformulation, and the implementation of adaptive finite elements. Furthermore, it is found that the resulting Newton solver is robust, successfully reducing the nonlinear residual to 1e-6 and lower in a small handful of iterations while addressing the challenging case to find a Taylor state equilibrium where conventional Picard-based solvers fail to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Palladium-Percolated Networks Enabled by Low Loadings of Branched Nanorods for Enhanced H 2 Separations

Nanoparticles (NPs) at high loadings are often used in mixed matrix membranes (MMMs) to improve gas separation properties, but they can lead to defects and poor processability that impede membrane fabrication. Herein, it is demonstrated that branched nanorods (NRs) with controlled aspect ratios can significantly reduce the required loading to achieve superior gas separation properties while maintaining excellent processability, as demonstrated by the dispersion of palladium (Pd) NRs in polybenzimidazole for H 2 /CO 2 separation. Increasing the aspect ratio from 1 for NPs to 40 for NRs decreases the percolation threshold volume fraction by a factor of 30, from 0.35 to 0.011. An MMM with percolated networks formed by Pd NRs at a volume fraction of 0.039 exhibits H 2 permeability of 110 Barrer and H 2 /CO 2 selectivity of 31 when challenged with simulated syngas at 200 °C, surpassing Robeson's upper bound. In conclusion, this work highlights the advantage of NRs over NPs and nanowires and shows that right-sizing nanofillers in MMMs is critical to construct highly sieving pathways at minimal loadings. In conclusion, this work paves the way for this general feature to be applied across materials systems for a variety of chemical separations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Thermal conductivity of multilayer polymer-nanocomposite thin films

The development of electrical insulators that are thermally conducting is critical for thermal management applications in many advanced electronics and electrical devices. Here, we synthesized polymer nanocomposite (PNC) films composed of polymers [polyethylenimine, poly(vinylamine), poly(acrylic acid), and poly(ethylene oxide)] and dielectric fillers (montmorillonite clay and hexagonal boron nitride) by layer-by-layer technique. The cross-plane thermal conductivity (k⊥) of the film was measured by the 3ω method. The effect of various factors such as film growth, filler type, filler volume fraction, polymer chemical structures, and temperature on the thermal conductivity is reported. The k⊥ of PNCs with thickness from 37 nm to 1.34 μm was found to be in the range of 0.11 to 0.21 ± 0.02 W m−1 K−1. The k⊥ values were found to be lower than the constituent polymer matrix. The experimental result is compared with existing theoretical models of nanocomposite systems to get insight into heat transfer behavior in such layered films composed of dielectrics and polymers.

Physics↗
DOE OSTI2022

Discovery top-down proteomics in symbiotic soybean root nodules

Proteomic methods have been widely used to study proteins in complex biological samples to understand biological molecular mechanisms. Most well-established methods (known as bottom-up proteomics, BUP) employ an enzymatic digestion step to cleave intact proteins into smaller peptides for liquid chromatography (LC) mass spectrometry (MS) detection. In contrast, top-down proteomics (TDP) directly characterizes intact proteins including all possible post-translational modifications (PTMs), thus offering unique insights into proteoform biology where combinations of individual PTMs may play important roles. We performed TDP on soybean root nodules infected by the symbiotic Bradyrhizobium japonicum in both the wildtype bacterium and a nifH- mutant, which lacks the ability to fix nitrogen in the soybean root nodule. TDP captured 1648 proteoforms derived from 313 bacterial genes and 178 soybean genes. Leghemoglobin, the most abundant protein in the sample, existed in many truncated proteoforms. Interestingly, these truncated proteoforms were considerably more abundant in the wildtype relative to the nifH- mutant, implicating protease activity as an important factor in nitrogen fixation. Proteoforms with various PTMs and combinations thereof were identified using an unrestricted open modification search. This included less common PTMs such as myristoylation, palmitoylation, cyanylation, and sulfation. In parallel, we collected high resolution MS imaging (MSI) data of intact proteins and biopolymers (<20 kDa due to current technical limitations) from sections of the soybean root nodules using matrix-assisted laser desorption/ionization (MALDI) coupled to high resolution Orbitrap. Several detected proteoforms exhibited unique spatial distributions inside the infection zone and cortex, suggesting functional compartmentalization in these regions. A subset of peaks from the MALDI-MSI were assigned to proteoforms detected in TDP LCMS data based on matching accurate masses. Many of the proteins detected in both LCMS and MALDI-MSI are currently uncharacterized in UniProt: the PTM and spatial information presented here will be valuable in understanding their biological functions. Taken together, our study demonstrates how untargeted TDP approach can provide unique insights into plant proteoform biology. On-going technology developments are expected to further improve TDP coverage for more comprehensive high-throughput analysis of proteoforms.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Scaling of Imbibition Front Dynamics in Heterogeneous Porous Media

Here we report a quantitative study on the dynamics of countercurrent imbibition in fractured rock samples with the application of CO 2 geosequestration, as the rock matrix is saturated with CO 2 and brine exists in the adjacent microfracture. Direct pore-scale simulations of spontaneous and forced imbibition in a three-dimensional heterogeneous porous matrix enable us to capture the effect of microscopic properties on the imbibition dynamics as these pore-level descriptions are ignored in the existing coarse-grained models. We have established new scaling classes for the propagation of the center of mass of the imbibition front and fluid-fluid interface broadening. It has become evident that the wettability is a major factor which leads to anomalies in these scaling laws. While the governing capillarity under the water-wet condition induces fast roughening dynamics of the interface, the adverse effect of capillary pressure under the intermediate-wet condition leads to a damping effect on the imbibition front fluctuations.

58 GEOSCIENCES↗
DOE OSTI2021

Intrinsic quality factor extraction of multi-port cavity with arbitrary coupling

In this work, we derived S-parameter based expressions for the intrinsic quality factor of an arbitrary coupled multiport microwave cavity and a non-ideal test fixture. Practical accuracy limitations of the obtained expressions specifically for superconducting accelerator cavities were evaluated both analytically and with the simulation software. The resulting formulas can be used to extract the intrinsic quality factor of normal-conducting and superconducting cavities directly from the calibrated S-parameter measurements.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Bottomonium suppression from the three-loop QCD potential

We compute the suppression of bottomonium in the quark-gluon plasma using the three-loop QCD static potential. The potential describes the spin-averaged bottomonium spectrum below threshold with a less than 1% error. Within potential nonrelativistic quantum chromodynamics and an open quantum systems framework, we compute the evolution of the bottomonium density matrix. The values of the quarkonium transport coefficients are obtained from lattice QCD measurements of the bottomonium in-medium width and thermal mass shift; we additionally include for the first time a vacuum contribution to the dispersive coefficient γ . Using the three-loop potential and the values of the heavy quarkonium transport coefficients, we find that the resulting bottomonium nuclear modification factor is consistent with experimental observations, while at the same time reproducing the lattice measurements of the in-medium width. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
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 OSTI2022

The response of accident tolerant fuel cladding to LOCA burst testing: A comparative study of leading concepts

Accident tolerant fuel claddings seek to improve safety margins during loss-of-coolant accident (LOCA) scenarios by reducing cladding steam oxidation rates, hydrogen production, delaying or mitigating cladding burst, and reducing cladding deformation. After rapid and extensive progress in recent years, Cr coated Zr-based alloys, ferritic FeCrAl, and SiC/SiC ceramic matrix composites (CMCs) have emerged as the most promising cladding concepts, but there are few studies comparing the accident tolerance of all three concepts simultaneously. Here, in this study, LOCA burst testing has been conducted on Cr high impulse power magnetron sputtered Zry-4 (Cr/Zry-4), C26M (FeCrAl), and SiC/SiC SiGA (CMC). Experimental observations indicate Cr/Zry-4 burst at higher temperatures than Zry-4 however the effectiveness of the coating on mitigating ZrO 2 formation via Cr 2 O 3 formation depended on physical proximity to the burst opening. No dramatic improvement in diametrical strain nor opening geometry was observed. C26M burst at higher temperatures than Zry-4 at lower pressures, exhibiting small amounts of strain and smaller openings, and formed slow-growing Al 2 O 3 regardless of location. The CMCs did not burst or show any macroscopic signs of deformation, all the way up to an internal cladding pressure of 15.2 MPa. High temperature expanding Nb plug testing indicated that CMC room temperature strength may persist up to 1650°C, with an approximately 40% loss of strength at 1900°C. CMC steam oxidation was also compared to reference chemically vapor deposited SiC up to 1700°C, and CMCs were found to oxidize at slightly accelerated rates that were still orders of magnitude lower than bare Zry-4. Projected H 2 generation during LOCA burst testing showed that when accounting for cladding burst exposing inner diameters to steam, Cr coated Zry-4 reduced H 2 generation by a factor of two compared with bare Zry-4, while both C26M and SiC/SiC reduced H 2 generation by orders of magnitude relative to bare Zry-4.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Insights into pulmonary phosphate homeostasis and osteoclastogenesis emerge from the study of pulmonary alveolar microlithiasis

Pulmonary alveolar microlithiasis is an autosomal recessive lung disease caused by a deficiency in the pulmonary epithelial Npt2b sodium-phosphate co-transporter that results in accumulation of phosphate and formation of hydroxyapatite microliths in the alveolar space. The single cell transcriptomic analysis of a pulmonary alveolar microlithiasis lung explant showing a robust osteoclast gene signature in alveolar monocytes and the finding that calcium phosphate microliths contain a rich protein and lipid matrix that includes bone resorbing osteoclast enzymes and other proteins suggested a role for osteoclast-like cells in the host response to microliths. While investigating the mechanisms of microlith clearance, we found that Npt2b modulates pulmonary phosphate homeostasis through effects on alternative phosphate transporter activity and alveolar osteoprotegerin, and that microliths induce osteoclast formation and activation in a receptor activator of nuclear factor-κB ligand and dietary phosphate dependent manner. This work reveals that Npt2b and pulmonary osteoclast-like cells play key roles in pulmonary homeostasis and suggest potential new therapeutic targets for the treatment of lung disease.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Accelerating LHC event generation with simplified pilot runs and fast PDFs

High-precision calculations are an indispensable ingredient for the success of the LHC physics programme, yet their poor computing efficiency has been a growing cause for concern, threatening to become a paralysing bottleneck in the coming years. We present solutions to eliminate the apprehension by focussing on two major components of generalpurpose Monte Carlo event generators: the evaluation of parton distribution functions, and the generation of perturbative matrix elements. We show that for the cost-driving event samples employed by the ATLAS experiment to model omnipresent, irreducible Standard Model backgrounds, such as weak boson or top-quark pair production in association with jets, these computational components dominate the overall run time by up to 80 %. We demonstrate that a reduction of the computing footprint of LHAPDF and SHERPA by factors of around 40 can be achieved for multi-leg NLO event generation.

Bothmann, Enrico [Gottingen U.]↗
DOE OSTI2024

Polarized 𝐽/𝜓 production in semi-inclusive DIS at large 𝑄 2 : Comparing quark fragmentation and photon-gluon fusion

We compare the relative importance of different mechanisms for polarized 𝐽/𝜓 production in semi-inclusive deep inelastic scattering processes at large 𝑄 2 . The transverse momentum dependent (TMD) factorization framework and nonrelativistic quantum chromodynamics are used to study the leading contributions from light quark fragmentation to polarized 𝐽/𝜓, and compared to direct production via photon-gluon fusion, which can proceed through color-singlet as well as color-octet mechanisms. We identify kinematic regimes where light quark fragmentation dominates, allowing for the extraction of the 3 𝑆$^{[8]}_{1}$ matrix element, as well as regimes where photon gluon fusion dominates, suggesting that the gluon TMD parton distribution function can be probed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Gauge-invariant renormalization of four-quark operators

We study the renormalization of four-quark operators in one-loop perturbation theory. We employ a coordinate-space gauge-invariant renormalization scheme (GIRS), which can be advantageous compared to other schemes, especially in nonperturbative lattice investigations. From our perturbative calculations, we extract the conversion factors between GIRS and the modified minimal subtraction scheme ( MS ¯ ) at the next-to-leading order. As a by-product, we also obtain the relevant anomalous dimensions in the GIRS scheme. A formidable issue in the study of the four-quark operators is that operators with different Dirac matrices mix among themselves upon renormalization. We focus on both parity-conserving and parity-violating four-quark operators, which change flavor numbers by two units ( Δ F = 2 ). The extraction of the conversion factors entails the calculation of two-point Green’s functions involving products of two four-quark operators, as well as three-point Green’s functions with one four-quark and two bilinear operators. The significance of our results lies in their potential to refine our understanding of QCD phenomena, offering insights into the precision of Cabibbo-Kobayashi-Maskawa (CKM) matrix elements and shedding light on the nonperturbative treatment of complex mixing patterns associated with four-quark operators. Published by the American Physical Society 2024

Constantinou, M. (ORCID:0000000269881745)↗