Search NASA⌕ Search

SEARCH · Search NASA

Results for “Factorization Techniques”

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 19 records

Separable physics-informed DeepONet: Breaking the curse of dimensionality in physics-informed machine learning

The deep operator network (DeepONet) has shown remarkable potential in solving partial differential equations (PDEs) by mapping between infinite-dimensional function spaces using labeled datasets. However, in scenarios lacking labeled data, the physics-informed DeepONet (PI-DeepONet) approach, which utilizes the residual loss of the governing PDE to optimize the network parameters, faces significant computational challenges, particularly due to the curse of dimensionality. This limitation has hindered its application to high-dimensional problems, making even standard 3D spatial with 1D temporal problems computationally prohibitive. Additionally, the computational requirement increases exponentially with the discretization density of the domain. Here, to address these challenges and enhance scalability for high-dimensional PDEs, we introduce the Separable physics-informed DeepONet (Sep-PI-DeepONet). This framework employs a factorization technique, utilizing sub-networks for individual one-dimensional coordinates, thereby reducing the number of forward passes and the size of the Jacobian matrix required for gradient computations. By incorporating forward-mode automatic differentiation (AD), we further optimize computational efficiency, achieving linear scaling of computational cost with discretization density and dimensionality, making our approach highly suitable for high-dimensional PDEs. We demonstrate the effectiveness of Sep-PI-DeepONet through three benchmark PDE models: the viscous Burgers’ equation, Biot’s consolidation theory, and a parameterized heat equation. Our framework maintains accuracy comparable to the conventional PI-DeepONet while reducing training time by two orders of magnitude. Notably, for the heat equation solved as a 4D problem, the conventional PI-DeepONet was computationally infeasible (estimated 289.35 h), while the Sep-PI-DeepONet completed training in just 2.5 h. These results underscore the potential of Sep-PI-DeepONet in efficiently solving complex, high-dimensional PDEs, marking a significant advancement in physics-informed machine learning.

Neural operator↗

Hybrid model predictive control techniques for safety factor profile and stored energy regulation while incorporating NBI constraints

Abstract A novel hybrid Model Predictive Control (MPC) algorithm has been designed for simultaneous safety factor ( q ) profile and stored energy ( w ) control while incorporating the pulse-width-modulation constraints associated with the neutral beam injection (NBI) system. Regulation of the q -profile has been extensively shown to be a key factor for improved confinement as well as non-inductive sustainment of the plasma current. Simultaneous control of w is necessary to prevent the triggering of pressure-driven magnetohydrodynamic instabilities as the controller shapes the q profile. Conventional MPC schemes proposed for q -profile control have considered the NBI powers as continuous-time signals, ignoring the discrete-time nature of these actuators and leading in some cases to performance loss. The hybrid MPC scheme in this work has the capability of incorporating the discrete-time actuator dynamics as additional constraints. In nonlinear simulations, the proposed hybrid MPC scheme demonstrates improved q -profile+ w control performance for NSTX-U operating scenarios.

Physics↗

A Measurement of the Neutron Electromagnetic Form Factor Ratio from a Rosenbluth Technique with Simultaneous Detection of Neutrons and Protons

The internal structure of protons and neutrons provides insight into both the dynamical behavior of the constitute quarks and gluons, and emergent properties of the nucleons (such as mass, spin, and electromagnetic distributions). Elastic electron-nucleon scattering can probe the elastic electromagnetic form factors of the nucleon. The electric and magnetic form factors, respectively, encode information about the internal charge and magnetization distributions within the nucleon. Precision data for these form factors, over a broad range of the four-momentum transfer squared, Q^2, can benchmark theoretical models describing the strong interaction of nuclear physics. The Super BigBite Spectrometer (SBS) program in Hall A at Jefferson Lab, is a series of high-precision experiments which seek to significantly extend the Q^2 reach of previous data for the nucleon electromagnetic form factors. The first two experiments of this program are known as G_M^n and the neutron Two Photon Exchange (nTPE) and the data were collected from October 2021 to February 2022. Both experiments were conducted with the simultaneous measurement of D(e,e'n) and D(e,e'p) reactions for quasi-elastic electron-deuteron scattering. The scattered electrons were detected in the BigBite Spectrometer, which features multiple large-acceptance Gas Electron Multiplier (GEM) detectors. The Super BigBite Spectrometer provided simultaneous detection of scattered nucleons, and utilized a large acceptance dipole magnet and Hadron Calorimeter (HCal). The G_M^n experiment provides precision measurements of the neutron magnetic form factor, via the ratio method, over a Q^2 range of 3.0 to 13.5 (GeV/c)2. From this data analysis, preliminary values for G_M^n/µ_n G_D are extracted. For Q^2=4.48 (GeV/c)2 we find G_M^n/µ_n G_D=0.9546±0.0132 and for Q^2=4.476 (GeV/c)2 we find G_M^n/µ_n G_D=0.9563±0.0110. These preliminary G_M^n/µ_n G_D values are more precise than existing world data in this Q^2 regime and are consistent with the most recent parameterization of the G_M^n/µ_n G_D world data. The nTPE experiment provides a first measurement of the neutron Rosenbluth Slope and seeks to quantify the two-photon exchange(TPE) contribution to elastic electron-neutron scattering at a fixed Q^2=4.5 (GeV/c)2 with two different beam energies and scattering angle values. For data of the proton form factor ratio, µ_p G_E^p/G_M^p, significant discrepancies exist between values obtained from Rosenbluth Separation and polarization transfer measurement, particularly at large Q^2, and TPE contributions are thought to resolve this discrepancy. The impacts of TPE contributions have not yet been experimentally established for the neutron. From the data analysis presented in this dissertation, a preliminary result for the neutron Rosenbluth Slope is found as S^n=(G_E^n )^2/t_n (G_M^n )^2=0.0916±0.0476 for Q^2=4.48 (GeV/c)2. This value of the neutron Rosenbluth Slope is consistent with the world data extrapolation and the absence of large TPE corrections.

Wertz, Ezekiel [Thomas Jefferson National Accelera↗

Answering old questions with new techniques: Understanding performance-limiting factors in transition metal dichalcogenide photoelectrochemical solar cells

In the late 1970s and early 1980s, several research groups reported high-efficiency and stable photoelectrochemical solar cells based on transition metal dichalcogenides (TMDS; e.g., MoS 2 , WSe 2 , and MoSe 2 ) immersed in iodide/tri-iodide electrolytes. A consensus emerged that smooth crystals were necessary for high efficiency based on significant evidence that rough crystals with exposed edge sites produced lower photocurrents and fill factors. However, anecdotal observations in the literature hinted at significant performance variation among apparently smooth crystals with the possibility of highly active “hot” edge sites. Furthermore, this mini-review article is a case study on how spatially resolved photoelectrochemical techniques developed in the 2020s are answering old questions regarding the origin of performance variation in high efficiency n-type TMD|I – ,I 3 – |Pt photoelectrochemical solar cells.

14 SOLAR ENERGY↗

GPU Accelerated Sparse Cholesky Factorization

The solution of sparse symmetric positive definite linear systems is an important computational kernel in large-scale scientific and engineering modeling and simulation. We will solve the linear systems using a direct method, in which a Cholesky factorization of the coefficient matrix is performed using a right-looking approach and the resulting triangular factors are used to compute the solution. Sparse Cholesky factorization is compute intensive. In this work we investigate techniques for reducing the factorization time in sparse Cholesky factorization by offloading some of the dense matrix operations on a GPU. We will describe the techniques we have considered. We achieved up to 4x speedup compared to the CPU-only version.

Karsavuran, M Ozan↗

Additive manufacturing of metal matrix composites

Although Metal matrix composites (MMCs) are superior to most sought-after metallic alloys, their challenging fabricability has limited their widespread use in bulk-form applications. Among the many advanced fabrication techniques, Additive Manufacturing (AM), owing to its unique capabilities to produce near-net shapes, has drawn significant traction in the past two decades, especially for materials that are difficult to process using traditional methods. However, unlike pure metal/alloy systems, MMCs are highly sensitive to the processing conditions prevailing in AM techniques due to factors such as the high melting point of reinforcement particles and the potential for in-situ reactions. Therefore, it may be a while before metal matrix composites are commercially produced via AM. This review will discuss the current state-of-the-art design, fabricability, and performance of various additively manufactured MMCs. A particular focus will be on microstructural evolution and microstructure-property relationships. The most employed AM techniques, such as directed energy deposition, powder bed fusion, binder jetting, sheet lamination, and solid-state friction stir processing, are fundamentally different in terms of thermo-kinetics, forming the perspective for this review. A detailed comparison of microstructural evolution and process parameter optimization, including feedstock preparation methods and the role of machine learning and modeling among the different AM processes, is also presented. Finally, a critical evaluation of emerging AM technologies for MMCs is also provided, highlighting their potential advantages and challenges.

36 - MATERIALS SCIENCE↗

The Frequency Shift and $Q$ of Disordered Superconducting RF Cavities

Niobium superconducting radio-frequency (SRF) cavities for high-energy accelerator applications have been greatly improved in terms of the quality factor $Q$ by techniques such as Nitrogen doping. However, the mechanisms leading to improvement in $Q$ are still not fully understood. Quite recently the SRF group at Fermilab measured anomalies in the frequency shift of N-doped SRF Niobium cavities near the transition temperature. Here we report our theoretical analysis of these results based on the microscopic theory of superconductivity that incorporates anisotropy of the superconducting gap and inhomogeneous disorder in the screening region of the SRF cavities. We are able to account for frequency shift anomalies very close to $T_c$ on the order of fractions of 1 kHz. Our results for the frequency shift and Q are in good agreement with the experimental data reported for all four N-doped Nb SRF cavities by Bafia et al. We also compare our theory with an earlier report on a Nb sample measured at 60 GHz. In addition, we show that the quality factor calculated theoretically has a peak of upper convexity with the largest $Q$ at intermediate levels of disorder. For strong disorder, i.e. the dirty limit, pair breaking in the presence of disorder and screening currents limits the $Q$.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Estimation of backgrounds from jets misidentified as τ-leptons using the Universal Fake Factor method with the ATLAS detector

Processes with τ$$\tau $$-leptons in the final state are important for Standard Model measurements and searches for physics beyond the Standard Model. The ATLAS experiment at the Large Hadron Collider observes τ$$\tau $$-leptons produced in proton–proton collisions only through their decay products. Data analyses involving hadronically decaying τ$$\tau $$-leptons face challenges due to backgrounds from jets misidentified as τ$$\tau $$-leptons that are not modelled reliably by Monte Carlo simulations. Data-driven methods such as the fake-factor method allow such misidentified backgrounds to be predicted by measuring transfer factors, known as fake factors, in data from dedicated regions. This paper describes a refined technique for determining the fake factors, the Universal Fake Factor method. It evaluates the fake factors for a signal region by using fake factors from samples enriched in different sources of jets misidentified as τ$$\tau $$-leptons (light-quark, gluon, b-quark, and pile-up jets). Each fake factor is calculated as a linear combination of fake factors measured in these different enriched samples. For the full Run 2 data set, the systematic uncertainty of the calculated fake factors, evaluated using W(μν)$$W(\mu u )$$ enriched event sample, ranges from 15 to 35% depending on the τ$$\tau $$-lepton’s transverse momentum and charged-particle decay multiplicity.

Aad, G↗

Analyzing the impact of design factors on solar module thermomechanical durability using interpretable machine learning techniques

Solar modules in utility-scale systems are expected to maintain decades of lifetime to rival conventional energy sources. However, cyclic thermomechanical loading often degrades their long-term performance, highlighting the importance of effective design to mitigate thermal expansion mismatches between module materials. Given the complex composition of solar modules, isolating the impact of individual components on overall durability remains a challenging task. In this work, we analyze a comprehensive data set that comprises bill-of-materials (BOM) and thermal cycling power loss from 251 distinct module designs to identify the predominant design factors and their impacts on the thermomechanical durability of modules. The methodology of our analysis combines machine learning modeling (random forest) and Shapley additive explanation (SHAP) to correlate design factors with power loss and interpret the model’s decision-making. The interpretation reveals that silicon type (monocrystalline or polycrystalline), encapsulant thickness, busbar numbers, and wafer thickness predominantly influence the degradation. With lower power loss of around 0.6% on average in the SHAP analysis, monocrystalline cells present better durability than polycrystalline cells. This finding is further substantiated by statistical testing on our raw data set. The SHAP analysis also demonstrates that while thicker encapsulants lead to reduced power loss, further increasing their thickness over around 0.6 to 0.7 mm does not yield additional benefits, particularly for the front side one. In addition, other important BOM features such as the number of busbars are analyzed. This study provides a blueprint for utilizing explainable machine learning techniques in a complex material system and can potentially guide future research on optimizing the design of solar modules.

14 SOLAR ENERGY↗

String-based parametrization of nucleon GPDs at any skewness: A comparison to lattice QCD

We introduce a string-based parametrization for nucleon quark and gluon generalized parton distributions (GPDs) valid at all skewness values. The conformal moments of the GPDs are expressed as sums of the spin-j nucleon A-form factor and the skewness-dependent spin-j nucleon D-form factor. This representation, which fulfills the polynomiality condition (due to Lorentz invariance) and does not rely on model-specific assumptions, is derived from t-channel string exchanges in anti-de Sitter spaces. The spin-j nucleon D-form factor is closely related to the spin-j nucleon A-form factor. We use the Mellin moments from empirical parton distributions to model the spin-j nucleon A-form factors. Using only five Regge slope parameters, fixed from the electromagnetic and gravitational form factors, our string-based parametrization generates accurate singlet, nonsinglet, isovector, and flavor-separated nucleon quark GPDs, along with symmetric nucleon gluon GPDs from their Mellin-Barnes integral representations. Our isovector nucleon quark GPD is in agreement with existing lattice data. Our string-based parametrization should facilitate the empirical extraction and global analysis of nucleon GPDs in exclusive processes, bypassing the deconvolution challenge.

Electron-ion collisions↗

Randomized Algorithms for Symmetric Nonnegative Matrix Factorization

Symmetric Nonnegative Matrix Factorization (SymNMF) is a technique in data analysis and machine learning that approximates a matrix with a product of a nonnegative, low-rank matrix and it transpose. To design faster and more scalable algorithms for SymNMF we develop two randomized algorithms for its computation. The first method uses randomized matrix sketching to compute an initial low-rank approximation to the input matrix and proceeds to uses this as a low-rank input to rapidly compute a SymNMF. The second methods uses randomized leverage score sampling to approximately solve constrained least squares problems. Many successful methods for SymNMF rely on (approximately) solving sequences of constrained least squares problems. Here, we prove theoretically that leverage score sampling can approximately solve constrained least squares problems to e-accuracy. Finally we demonstrate both methods work in practice by applying them to graph clustering tasks on large real world data sets. These experiments show that our methods approximately maintain solution quality and achieve significant speed ups for both large dense and large sparse problems.

97 MATHEMATICS AND COMPUTING↗

Experimental Examination of Additively Manufactured Patterns on Structural Nuclear Materials for Digital Image Correlation Strain Measurements

Abstract Background There are a limited number of commercially available sensors for monitoring the deformation of materials in-situ during harsh environment applications, such as those found in the nuclear and aerospace industries. Such sensing devices, including weldable strain gauges, extensometers, and linear variable differential transformers, can be destructive to material surfaces being investigated and typically require relatively large surface areas to attach (> 10 mm in length). Digital image correlation (DIC) is a viable, non-contact alternative to in-situ strain deformation. However, it often requires implementing artificial patterns using splattering techniques, which are difficult to reproduce. Objective Additive manufacturing capabilities offer consistent patterns using programmable fabrication methods. Methods In this work, a variety of small-scale periodic patterns with different geometries were printed directly on structural nuclear materials (i.e., stainless steel and aluminum tensile specimens) using an aerosol jet printer (AJP). Unlike other additive manufacturing techniques, AJP offers the advantage of materials selection. DIC was used to track and correlate strain to alternative measurement methods during cyclic loading, and tensile tests (up to 1100 µɛ) at room temperature. Results The results confirmed AJP has better control of pattern parameters for small fields of view and facilitate the ability of DIC algorithms to adequately process patterns with periodicity. More specifically, the printed 100 μm spaced dot and 150 μm spaced line patterns provided accurate measurements with a maximum error of less than 2% and 4% on aluminum samples when compared to an extensometer and commercially available strain gauges. Conclusion Our results highlight a new pattern fabrication technique that is form factor friendly for digital image correlation in nuclear applications.

Novich, K. A. (ORCID:0000000204466022)↗

Accelerating the discovery of low-energy structure configurations: A computational approach that integrates first-principles calculations, Monte Carlo sampling, and Machine Learning

Finding Minimum Energy Configurations (MECs) is essential in fields such as physics, chemistry, and materials science, as they represent the most stable states of the systems. In particular, identifying such MECs in multi-component alloys considered candidate PFMs is key because it determines the most stable arrangement of atoms within the alloy, directly influencing its phase stability, structural integrity, and thermo-mechanical properties. However, since the search space grows exponentially with the number of atoms considered, obtaining such MECs using computationally expensive first-principles DFT calculations often results in a cumbersome task. To escape the above compromise between physical fidelity and computational efficiency, we have developed a novel physics-based data-driven approach that combines Monte Carlo sampling, first-principles DFT calculations, and Machine Learning to accelerate the discovery of MECs in multi-component alloys. More specifically, we have leveraged well-established Cluster Expansion (CE) techniques with Local Outlier Factor models to establish strategies that enhance the reliability of the CE method. In this work, we demonstrated the capabilities of the proposed approach for the particular case of a tungsten-based quaternary high-entropy alloy. However, the method is applicable to other types of alloys and enables a wide range of applications.

36 MATERIALS SCIENCE↗

Indirect Measurement of the 23 Na(p, γ) 24 Mg Direct Capture Reaction Rate via ( 3 He,d) Spectroscopy

The cross section of the $^{23}\text{Na}(p,γ)^{24}\text{Mg}$ reaction is dominated by direct capture at low energies relevant for stellar burning. Such cross sections can be constrained using spectroscopic factors($C^2S$) or asymptotic normalization coefficients(ANCs) from transfer reactions. In this work, the $^{23}\text{Na}(^3\text{He},d)^{24}\text{Mg}$ reaction was measured at $E_{lab}=21$ MeV to extract spectroscopic factors for $^{24}\text{Mg}$ states with excitation energies in $E_x=7 - 12$ MeV using the Enge split-pole spectrograph at the Triangle Universities Nuclear Laboratory. A new non-resonant astrophysical S factor and the direct capture reaction rate for the $^{23}\text{Na}(p,γ)$ reaction are calculated and presented based on this measurement. The new rate at $T<0.04$ GK is 43$\%$ smaller than in previous studies. Rigorous treatments of uncertainties are presented using a Bayesian Markov Chain Monte Carlo (MCMC) method. Finally, sources of uncertainties for computing the direct capture cross section are also discussed in detail.

electromagnetic moments↗

Learning the factors controlling mineral dissolution in three-dimensional fracture networks: applications in geologic carbon sequestration

We perform a set of high-fidelity simulations of geochemical reactions within three-dimensional discrete fracture networks (DFN) and use various machine learning techniques to determine the primary factors controlling mineral dissolution. The DFN are partially filled with quartz that gradually dissolves until quasi-steady state conditions are reached. At this point, we measure the quartz remaining in each fracture within the domain as our primary quantity of interest. We observe that a primary sub-network of fractures exists, where the quartz has been fully dissolved out. This reduction in resistance to flow leads to increased flow channelization and reduced solute travel times. However, depending on the DFN topology and the rate of dissolution, we observe substantial variability in the volume of quartz remaining within fractures outside of the primary subnetwork. This variability indicates an interplay between the fracture network structure and geochemical reactions. We characterize the features controlling these processes by developing a machine learning framework to extract their relevant impact. Specifically, we use a combination of high-fidelity simulations with a graph-based approach to study geochemical reactive transport in a complex fracture network to determine the key features that control dissolution. We consider topological, geometric and hydrological features of the fracture network to predict the remaining quartz in quasi-steady state. We found that the dissolution reaction rate constant of quartz and the distance to the primary sub-network in the fracture network are the two most important features controlling the amount of quartz remaining. This study is a first step towards characterizing the parameters that control carbon mineralization using an approach with integrates computational physics and machine learning.

54 ENVIRONMENTAL SCIENCES↗

Simultaneous control of the electron temperature and safety factor profiles in DIII-D using model-based optimal control techniques

Future tokamak power plants will likely operate using a single, well-defined plasma scenario, either in steady state or for very long pulse lengths. In order to enhance the robustness of the scenario, feedback controllers for a variety of plasma properties will be necessary to counteract any disturbances and ensure safe operation. However, only a limited set of actuators will be available to control many different quantities. Because of this, it is necessary to develop controllers that are able to regulate multiple plasma properties using a limited set of actuators. To this end, a controller has been developed for the simultaneous regulation of both the electron temperature and safety factor profiles in DIII-D. This algorithm uses a linear quadratic integral control synthesis approach based on a linearized model of the dynamics of the two profiles. Two neural network surrogate models, NubeamNet and MMMnet, are included to improve the fidelity of the model. Furthermore, the controller has been tested in simulation using COTSIM, and has demonstrated the ability to simultaneously track changes in both the electron temperature and safety factor targets, including changes in both the magnitude and the shape of the profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Measurement of the Neutron Magnetic Form Factor at Large Momentum Transfer Using the Super-Bigbite Apparatus in Jefferson Lab Hall-A

The nucleon elastic electromagnetic form factors help us study the electromagnetic structure of the nucleon, benchmark theoretical models, and improve our understanding of non-perturbative quantum chromodynamics and confinement. The Nobel Prize-winning electron-nucleon scattering experiments by Robert Hofstadter and collaborators in the 1950s at Stanford High Energy Physics Lab were the first nucleon form factor measurements performed using leptonic probes. The Super Bigbite Spectrometer (SBS) program at Hall-A of Jefferson Lab represents the latest efforts to measure nucleon form factors. This ambitious program aims to significantly extend the current data set in terms of square momentum transfer (Q2) with high precision. The advent of novel detector technologies, like Gas Electron Multipliers (GEM), which provide excellent position resolution (< 100 ?m) while withstanding high background particle rates (several hundred MHz/cm2) over a large active area, has paved the way for open-geometry, moderate solid angle spectrometers, which are central to all form factor experiments in the SBS program. The first experimental run group in SBS ran successfully between September 2021 and February 2022, collecting data for the measurement of the magnetic form factor of the neutron Gn M at five squared momentum-transfer values: 3.0, 4.5, 7.5, 9.8, and 13.5 (GeV /c)2. This extends the existing high-precision data for Gn M by about a factor of four. The ratio technique was used, which involved the simultaneous measurement of exclusive quasielastic scattering of D(e,e?n)p and D(e,e?p)n from a deuterium target. Pre-preliminary results for D(e,e?n)p and D(e,e?p)n quasi-elastic ratio, and the neutron magnetic form factor Gn M , for Q2 points 3.0, 9.8, and 13.5 (GeV /c)2 are presented.

Mudiyanselage, Anuruddha Rathnayake↗

Developing a Prototype Methodology to Rank CO2-EOR Wells and Assess Their Reuse Potential for Geologic Carbon Storage

This paper presents a prototype methodology to assess the possible transition of Class II carbon dioxide-enhanced oil recovery (CO2-EOR) wells to Class VI wells. The focus is on wellbore construction materials—casing, cement, tubing, and the packer—and includes comprehensive workflows to evaluate these materials, with primary emphasis on compliance with Environmental Protection Agency (EPA) Class VI well construction and conversion guidelines. These workflows systematically assess material properties and performance criteria to ensure regulatory compliance and optimize long-term wellbore integrity and functionality. Utilizing Python scripts and JavaScript Object Notation (JSON) representations, the study automates checks on digitized Texas Railroad Commission (TRRC) data to rank wells based on workflow criteria. By emphasizing critical factors such as casing integrity, cementing techniques, tubing compatibility, and packer selection, the methodology helps well owners and operators prioritize wells for potential reuse as CO2 injection wells. Given limitations in digitized data, manual user verification is required in some sections. Future improvements include integrating non-digitized data through web scraping and machine learning techniques. This research serves as a practical guide for stakeholders, supporting environmental compliance and sustainable well operations.

geologic carbon sequestration↗