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 1,333 records · Page 74

FY24 Progress Report: Disposition of Cracks & Features Observed in the Inner Can Closure Weld Region (ICCWR) of the 3013 Package

The long-term integrity of the 3013 containers is of interest for the safe storage of Pu materials. Although the 3013 standard embodies multiple barrier concept, the integrity of the inner container is considered to be crucial since it protects the Safety Class outer container from its contents. The container is designed to withstand high pressures that could result from the complete radiolysis of the maximum water content permissible by the 3013 standard. Shelf-life studies and destructive examination have not found high gas pressures but have found that corrosive gases are generated. Pitting and stress corrosion cracking (SCC) are considered to be critical corrosion modes for the performance of the inner container and have been observed during destructive examinations. Considerable research has been conducted on the possibility of aqueous electrolytes condensing on the inner container that can promote SCC. These studies indicate the uncertainties surrounding the formation of corrosive environments and the corrosion behavior of container materials. In this initial report, a Bayesian network (BN) model is described that can consider the uncertainties and the causal connections between various factors influencing the corrosion modes of the inner container. The BN model is preliminary and provides an initial framework to identify the necessary information. The report also provides initial experimental results on the electrochemical behavior of stainless steels in anticipated condensed environments from gas phase migration of acidic gases. The experimental results are consistent with the corrosion model. Recommendation for further work on the BN model include assembling an expert group to provide input to the BN structure and quantification of the conditional probability matrix, experimental studies to characterize the microstructure of the container, electrochemical studies to identify critical potentials for localized corrosion and SCC, and crack growth rate studies

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Projection methods for the numerical solution of Markov chain models

Projection methods for computing stationary probability distributions for Markov chain models are presented. A general projection method is a method which seeks an approximation from a subspace of small dimension to the original problem. Thus, the original matrix problem of size N is approximated by one of dimension m, typically much smaller than N. A particularly successful class of methods based on this principle is that of Krylov subspace methods which utilize subspaces of the form span(v,av,...,A(exp m-1)v). These methods are effective in solving linear systems and eigenvalue problems (Lanczos, Arnoldi,...) as well as nonlinear equations. They can be combined with more traditional iterative methods such as successive overrelaxation, symmetric successive overrelaxation, or with incomplete factorization methods to enhance convergence.

Saad, Youcef↗

Scaling whole-chip QAOA for higher-order ising spin glass models on heavy-hex graphs

Abstract We show that the quantum approximate optimization algorithm (QAOA) for higher-order, random coefficient, heavy-hex compatible spin glass Ising models has strong parameter concentration across problem sizes from 16 up to 127 qubits for p = 1 up to p = 5, which allows for computationally efficient parameter transfer of QAOA angles. Matrix product state (MPS) simulation is used to compute noise-free QAOA performance. Hardware-compatible short-depth QAOA circuits are executed on ensembles of 100 higher-order Ising models on noisy IBM quantum superconducting processors with 16, 27, and 127 qubits using QAOA angles learned from a single 16-qubit instance using the JuliQAOA tool. We show that the best quantum processors find lower energy solutions up to p = 2 or p = 3, and find mean energies that are about a factor of two off from the noise-free distribution. We show that p = 1 QAOA energy landscapes remain very similar as the problem size increases using NISQ hardware gridsearches with up to a 414 qubit processor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Carbon Transport and Storage Planning and Viability Support Tools

The EDX disCO2ver Carbon Transport and Storage Planning and Viability Support Tools are made up of the Carbon Storage Planning Inquiry Tool (CS PlanIT, Justman et al. 2024) and the Carbon Storage Technical Viability Approach Support Tool (CS TVA). Together, these tools support data access to support understanding data availability to support planning efforts for carbon transport and storage. The Carbon Storage Planning Inquiry Tool (CS PlanIT) is an online web mapping application designed to help users explore, query, and evaluate multiple data layers to support and accelerate carbon storage resource and feasibility assessments and planning efforts. CS PlanIT currently contains a range of datasets associated with geologic, technical, and infrastructure factors. The data sets can be filtered geographically for an area of interest to update statistics and charts within the dashboard. The dashboard is divided into different sections called widgets, relating to different steps in the carbon storage planning process. The resources in this submission include a link to PlanIT, as well as a data catalog and link to user documentation. The original citation for the CS PlanIT tool, which has now been integrated into the toolset here, was: - Devin Justman, Scott Pantaleone, Maneesh Sharma, Lucy Romeo, Paige Morkner, CS PlanIT (Carbon Storage Planning Inquiry Tool) , 6/28/2024, https://edx.netl.doe.gov/dataset/cs-planit-carbon-storage-planning-inquiry-tool, DOI: 10.18141/2377953 The Carbon Storage Technical Viability Approach Support (CS TVA) Tool displays spatial data availability for the many components of Geologic Carbon Storage (GCS) technical viability assessment (Creason al 2025). Identifying sites suitable for GCS requires evaluating the intersection of myriad factors, including reservoir conditions, subsurface and surface hazards, infrastructure requirements, and energy community metrics. The technical viability of a site can only be confirmed for instances where all these factors have data available, and where those data support viability. Additional Resources related to the Technical Viability Assessment Tool: - Julia Mulhern, Casey White, Araceli Lara, Neyda Cordero Rodriguez, Zachary Jackson, Jacob Shay, Gabriel Creason, MacKenzie Mark-Moser, Paige Morkner, Kelly Rose, Carbon Storage Technical Viability Approach (CS TVA) Database, 3/26/2025, https://edx.netl.doe.gov/dataset/edx4ccs-carbon-storage-technical-viability-approach-database , DOI:10.18141/1984655 - Julia Mulhern, MacKenzie Mark-Moser, Gabriel Creason, Casey White, Araceli Lara, Neyda Cordero Rodriguez, Zach Jackson, Paige Morkner, Kelly Rose, Carbon Storage Technical Viability Approach (CS TVA) Matrix, 3/27/2025, https://edx.netl.doe.gov/dataset/carbon-storage-technical-viability-approach-cs-tva-matrix , DOI: 10.18141/2539979 - Gabriel Creason, Zach Jackson, Neyda Cordero Rodriguez, Julia Mulhern, Casey White, Araceli Lara, MacKenzie Mark-Moser, Paige Morkner, Kelly Rose, Carbon Storage Technical Viability Approach (CS TVA) Data Availability Results Database, 3/27/2025, https://edx.netl.doe.gov/dataset/carbon-storage-technical-viability-approach-cs-tva-data-availability-results-database, DOI:10.18141/2538557

Carbon storage↗

Inclusion of elastically connected members in exact buckling and frequency calculations

A standard stiffness matrix procedure which permits any combination of rigid, elastic, pinned or sliding connections of the degrees of freedom at the ends of a member to the nodes of its parent structure is described, in order to show how easily it can be extended to allow an existing algorithm to be used to ensure that no eigenvalues of the parent structure can be missed even when 'exact' member theory is used. The eigenvalues are the natural frequencies of undamped free vibration analyses or the critical load factors of buckling problems. The method preserves the exactness of the member theory and an efficient method for computer application is indicated. The theory also permits any combination of rigid, elastic, pinned or sliding connections between the freedoms of a substructure and those of its parent structure.

Williams, F. W.↗

Modeling and simulation of vascular tumors embedded in evolving capillary networks

Here in this work, we present a coupled 3D–1D model of solid tumor growth within a dynamically changing vascular network to facilitate realistic simulations of angiogenesis. Additionally, the model includes erosion of the extracellular matrix, interstitial flow, and coupled flow in blood vessels and tissue. We employ continuum mixture theory with stochastic Cahn–Hilliard type phase-field models of tumor growth. The interstitial flow is governed by a mesoscale version of Darcy’s law. The flow in the blood vessels is controlled by Poiseuille flow, and Starling’s law is applied to model the mass transfer in and out of blood vessels. The evolution of the network of blood vessels is orchestrated by the concentration of the tumor angiogenesis factors (TAFs); blood vessels grow towards the increasing TAFs concentrations. This process is not deterministic, allowing random growth of blood vessels and, therefore, due to the coupling of nutrients in tissue and vessels, makes the growth of tumors stochastic. We demonstrate the performance of the model by applying it to a variety of scenarios. Numerical experiments illustrate the flexibility of the model and its ability to generate satellite tumors. Simulations of the effects of angiogenesis on tumor growth are presented as well as sample-independent features of cancer.

3D–1D coupled blood flow models↗

Jacobian-based model diagnostics and application to equation oriented modeling of a carbon capture system

It can be difficult to identify the specific variables or equations responsible for convergence issues in large mathematical programming models. The Institute for the Design of Advanced Energy Systems Integrated Platform (IDAES-IP) contains a tool to identify poorly scaled constraints and variables by searching for rows and columns of the Jacobian matrix with small L2-norms. A singular value decomposition is then performed to identify degenerate sets of equations and remaining scaling issues. Here, this work presents a flowsheet developed for post-combustion carbon capture using a monoethanolamine (MEA) solvent system as a case study. This work takes the reader through the entire process of model diagnostics and reformulation, from a basic introduction to the mathematics behind these model diagnostics to the reformulations necessary to make the model numerically robust, including a significantly modified enhancement factor model.

IDAES↗

Standalone Block Copolymer Nanoballoons: Decoupling Self-Assembly from Implementation in Nanomanufacturing

Here we report a facile method to produce isolatable hollow-core elastomeric vesicles, “nanoballoons”, prepared via block copolymer self-assembly in a polymer blend. Poly(isoprene-block-dimethylsiloxane) (PI-PDMS) diblock copolymers are blended with PDMS homopolymers (h-PDMS) as a “solvent” phase to template the self-assembly of PDMS-tethered vesicles with PI walls. The walls are subsequently crosslinked to yield mechanically stabilized elastomeric vesicles. The h-PDMS inner-core and matrix are separated from the vesicles by dialysis to yield the matrix-free nanoballoons. These objects, 0.3 – 1 μm in diameter, can be further reincorporated into a crosslinkable PDMS. Throughout the self-assembly, recovery, and reincorporation processes we apply several techniques including solvent dispersion/dynamic light scattering (DLS) measurements, transmission electron microscopy (TEM)/energy-dispersive X-ray spectroscopy (EDS), and small-angle X-ray scattering (SAXS) to provide a consilient body of evidence that the nanoballoon morphology is retained. Furthermore, this work presents advanced nanomanufacturing schema that illustrate the decoupling of the thermodynamic and dynamic factors that govern macromolecular self-assembly from the environment in which the self-assembled objects are deployed.

36 MATERIALS SCIENCE↗

Complex Fluid‐Driven Fractures Caused by Crack‐Parallel Stress

Abstract Managing fluid‐driven fracture networks is crucial for subsurface resource utilization, yet the current understanding of the key controlling factors remains insufficient. While geologic discontinuities have been shown to significantly influence fracture network complexity, this study identifies another major contributor. We conducted a new set of experiments using a transparent true triaxial cell, which enabled video recording of the temporal evolution of fluid‐driven fracture paths. Using pseudo‐2D samples without macroscale structural discontinuities, we observed multiple occurrences of hydraulic fracture curving and branching under anisotropic boundary stresses. We proposed a theoretical model demonstrating that the stress parallel to the crack line in the solid matrix near the crack tip (i.e., the T ‐stress) accounts for the observed fracture curving behavior. This finding suggests that T ‐stress is an additional mechanism contributing to the complexity of fluid‐driven fracture networks in the subsurface, besides the geologic discontinuities.

58 GEOSCIENCES↗

Probing anisotropies of the Stochastic Gravitational Wave Background with LISA

Here, we investigate the sensitivity of the Laser Interferometer Space Antenna (LISA) to the anisotropies of the Stochastic Gravitational Wave Background (SGWB). We first discuss the main astrophysical and cosmological sources of SGWB which are characterized by anisotropies in the GW energy density, and we build a Signal-to-Noise estimator to quantify the sensitivity of LISA to different multipoles. We then perform a Fisher matrix analysis of the prospects of detectability of anisotropic features with LISA for individual multipoles, focusing on a SGWB with a power-law frequency profile. We compute the noise angular spectrum taking into account the specific scan strategy of the LISA detector. We analyze the case of the kinematic dipole and quadrupole generated by Doppler boosting an isotropic SGWB. We find that β Ω GW ~ 2 × 10 -11 is required to observe a dipolar signal with LISA. The detector response to the quadrupole has a factor ~ 10 3 β relative to that of the dipole. The characterization of the anisotropies, both from a theoretical perspective and from a map-making point of view, allows us to extract information that can be used to understand the origin of the SGWB, and to discriminate among distinct superimposed SGWB sources.

79 ASTRONOMY AND ASTROPHYSICS↗

Dynamical correlations and nonequilibrium sum rules in photodoped Hubbard ladders

Using matrix product state techniques we study the nonequilibrium dynamical response of the half-filled Hubbard ladder when subject to an optical pump. Optical pumping offers a way of producing and manipulating new strongly correlated phenomena by suppressing existing magnetic correlations. The ladder allows the effects of pump directionality to be investigated, and compared to a single chain it has strong spin-charge coupling and a fully gapped excitation spectrum, promising different nonequilibrium physics. We compute time-dependent correlations, including the nonequilibrium dynamical structure factors for spin and charge. By deriving a combined spin-charge sum rule that applies both in and out-of-equilibrium, we show that spectral weight is pumped directly from the antiferromagnetic spin response into a low energy 𝜔 ∼ 0 charge response below the Mott gap. The transfer of weight is pump direction dependent: pumping directed along the legs disrupts magnetic correlations more than pumping in the rung-direction, even if the post-pump energy density is similar. The charge correlation length is dramatically enhanced by the pump, whilst the spin correlations are most strongly suppressed at nearest- and next-nearest-neighbor spacings. After the pump the system is in a nonthermal correlated metallic state, with gapless charge excitations and approximately equal spin and charge correlation lengths, emphasizing the importance of treating these degrees of freedom on an equal footing in nonequilibrium systems.

36 MATERIALS SCIENCE↗

Development of Implicit Methods in CFD NASA Ames Research Center 1970's - 1980's

The focus here is on the early development (mid 1970's-1980's) at NASA Ames Research Center of implicit methods in Computational Fluid Dynamics (CFD). A class of implicit finite difference schemes of the Beam and Warming approximate factorization type will be addressed. The emphasis will be on the Euler equations. A review of material pertinent to the solution of the Euler equations within the framework of implicit methods will be presented. The eigensystem of the equations will be used extensively in developing a framework for various methods applied to the Euler equations. The development and analysis of various aspects of this class of schemes will be given along with the motivations behind many of the choices. Various acceleration and efficiency modifications such as matrix reduction, diagonalization and flux split schemes will be presented.

Pulliam, Thomas H.↗

Determining the Most Influencing Medical Conditions in MEDPRAT’s SIN Directed Graph

INTRODUCTION: The Susceptibility Inference Network (SIN) is a network of medical conditions, part of the Medical Extensible Probabilistic Risk Assessment Tool (MEDPRAT) developed by NASA to assess human health and medical risk to space exploration missions. The SIN is subject matter expert informed and acts as a prototype that provides relationships and dependencies between events modeled by MEDPRAT. Each vertex in the SIN has a weight which evaluates the severity of having the condition regardless of the progression from or to that condition. In this presentation, we consider two statistics to measure that stand alone risk: Quality Time Lost (QTL) and Loss of Crew Life (LOCL). Our goal is to identify the medical conditions that contribute the most to crew members QTL and LOCL risks due to progression of conditions in the network. We investigate how different computation parameters result in different condition rankings and address the choice of parameters that allows appropriate interventions to ensure space mission success. METHODS: The Katz score, one of many centrality measures created for ranking purposes in network analysis, takes into account all possible walks through the network, penalizing each additional step in a walk by a factor α called the Katz parameter. The literature does not provide specific values for the choice of α. We derive an analytical relationship between α and the maximum path length which has influence on the Katz score and ranking. Based on the probability of progression of each condition in the SIN, we identify that maximum path length of interest and calculate α that is then used in the Katz formula to rank the conditions in the SIN. RESULTS AND CONCLUSION: The effective probabilities of the SIN matrix generally fall below 10−6, which is below the level of the least influencing condition in the set. This corresponds to the probability of at most six consecutive progressions of a condition. Consequently, we calculate the Katz Parameter α and get 0.32. We rank the medical conditions and find that Acute Radiation Symptom is the condition the most prone to contribute to quality time loss due to progression.

risk analysis↗

Improved Evaluation of Large Network Matrices for Linear Power Flow Within Optimization Problems: Preprint

This work discusses methods for evaluating the Power Transfer Distribution Factor (PTDF) and Line Outage Distribution Factor (LODF) matrices by employing sparse linear algebra for large-scale computing applications. These matrices are critical in many power systems applications, such as the Unit Commitment Problem (UC), pre- and post-contingency power flow analysis, and transmission expansion. These matrices are typically dense, which means they require a significant amount of time and memory to be computed for large networks. However, by analyzing the structure of the matrices and their computation method, it is possible to use reduced memory methods based on sparse matrix operations. This paper shows that sparse linear algebra algorithms are faster and require less memory and time than traditional dense approaches. Additionally, we explore the effect of matrix sparsification by eliminating trailing digits on power flow calculations.

ENERGY PLANNING, POLICY, AND ECONOMY↗

Annual report for DOE VTO

Carbon fiber (CF)/polymer composites are a transformative class of high-performance, lightweight material, where high aspect-ratio CFs reinforce a polymer matrix and exceed the strength of steel alloys at a fraction of the density. Despite the advantages of such a class of material, the broader implementation of CF composites in a range of automotive, aerospace, and energy applications is hindered by limitations of current manufacturing methods. These current techniques (e.g., hand lay-up, wet filament winding) are costly and impose severe limitations on fiber placement, orientation, and angle, and thus a composite’s ultimate properties. Today’s CF composites are expensive to manufacture, limited in form factor, and utilize costly and sub-optimal continuous filament CF. Advanced additive manufacturing (AM) processes, combined with computational design optimization and new approaches to resin development, offer alternative design and manufacturing paradigms that have the realistic potential to lift these constraints. Such integrated AM approaches could thus help to realize the full potential of CF composite materials. One relevant application of CF composite materials where manufacturing constraints limit the cost-benefit ratio is in the manufacture of high-performance composite pressure vessels for onboard compressed natural gas (CNG) storage. Current CNG storage vessels (Types 3–5) are made from load-bearing filament-wound carbon-fiber composite and are ~3.5 times as expensive as an all-metallic Type-1 vessel. This cost is invariably tied to the complex and labor-intensive nature of conventional filament winding processes and the large volumes of expensive high tensile-strength CF tow feedstock required in manufacture. Our proposed approach to CNG storage vessel manufacture is based on a combination of AM technologies for CF composite printing and design optimization tools that were pioneered at Lawrence Livermore National Laboratory (LLNL), with advances in resin/composite formulation enabled by chemical and nano-material modification. Through the successful development of this technology, LLNL seeks to demonstrate the capability for advanced CNG storage vessel manufacture at reduced cost with no reduction in performance versus the most advanced, extant Type-5 designs.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Organic fiber/epoxy pressure vessels

We evaluated the performance of an organic fiber in an epoxy matrix by winding 20-cm diam spherical and cylindrical pressure vessels of various designs. For the spherical vessels, we used soft aluminum liners 0.76 mm thick for a double boss design and 2 mm thick for a single boss design. For the cylindrical vessels, we used both 0.5-mm rubber liners and 0.76-mm soft aluminum liners. Vessels of both types were tested for burst pressure and cyclic fatigue at room temperature and liquid hydrogen temperature. The effects of temperature and vessel shape on the vessel performance factor were negligible. Our vessel fatigue data were marred by premature failure of the liners.

Chiao, T. T.↗