Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “factorization methods”

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 OSTI2025

Precision Measurement of the Neutron Magnetic Form Factor via the Ratio Method at Jefferson Lab Hall A

Protons and neutrons, collectively known as nucleons, are composed of quarks and gluons. The Sachs electromagnetic form factors encode information about the spatial distributions of charge and magnetization in the nucleon, particularly at low momentum transfer. In particular, the neutron magnetic form factor (GMn) provides crucial information about the distribution of magnetization inside the neutron and helps constrain theoretical models of nucleon structure. Quasi-elastic electron scattering from deuterium was measured up to Q^2=13.5 GeV^2 using the Super BigBite Spectrometer in Hall A at Jefferson Lab. In this work, the neutron magnetic form factor GMn was extracted at Q^2 = 3.0 GeV^2 and Q^2=4.5 GeV^2 using the Ratio Method. These results represent a subset of the full dataset collected in this experiment, which extended to significantly higher Q^2. The extracted GMn values agree with the existing global fit within approximately two standard deviations at Q^2=3.0 and show excellent agreement at Q^2=4.5. The measurements achieved systematic uncertainties of about 2% and statistical uncertainties below 0.5%, among the most precise determinations of GMn at these kinematics. These results demonstrate the robustness of the experimental technique and provide an important validation point for future extractions at higher Q^2, where data remain scarce. In addition, the GRINCH heavy gas Cherenkov detector—a key component of the experimental apparatus—was commissioned and achieved an electron detection efficiency of approximately 97%, supporting reliable particle identification. Together, the analysis presented here advances both our understanding of nucleon structure and the validation of the experimental methods and instrumentation used to access it.

Satnik, Maria [College of William and Mary, Willia↗
DOE OSTI2023

Sparse Approximate Multifrontal Factorization with Composite Compression Methods

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

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Different flavors of exact-factorization-based mixed quantum-classical methods for multistate dynamics

Here, the exact factorization approach has led to the development of new mixed quantum-classical methods for simulating coupled electron–ion dynamics. We compare their performance for dynamics when more than two electronic states are occupied at a given time, and analyze: (1) the use of coupled versus auxiliary trajectories in evaluating the electron–nuclear correlation terms, (2) the approximation of using these terms within surface-hopping and Ehrenfest frameworks, and (3) the relevance of the exact conditions of zero population transfer away from nonadiabatic coupling regions and total energy conservation. Dynamics through the three-state conical intersection in the uracil radical cation as well as polaritonic models in one dimension are studied.

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

Improvement of the AC orbit bump method of local kick factor measurement using reference bunches

Here we have improved the AC orbit bump technique of local kick factor measurement. This technique, developed at NSLS-II a few years ago, is based on local sine-wave excitation of the electron beam by four fast correctors adjacent to the impedance location. By implementing reference bunches and the time gate function of beam position monitors, we are now able to measure the AC orbits of bunches with different intensities simultaneously. This improvement allows us to obtain the impedance-driven orbit distortion in a single measurement, significantly improving the measurement accuracy and reducing the measurement time.

36 MATERIALS SCIENCE↗
DOE OSTI2023

GeoThermalCloud: A Machine Learning Tool for Discovery, Exploration, and Development of Hidden Geothermal Resources

In this 25 minute presentation, we showcase our open source “GeoThermalCloud” tool for identifying hidden geothermal resources using a publicly available dataset for southwestern New Mexico. The presenters include Bulbul Ahmmed and Luke Frash. All of the visuals use source material from LA-UR approved publications and this work falls under the Earth Sciences DUSA. The code shown in this video is already released with LANL approval in open source format on GitHub and DockerHub. The audio in this video includes only material on the topics of geothermal energy and machine learning applied to geothermal energy. The primary machine learning method used is LANL’s Non-negative Matrix Factorization “NMFk” method. Modeling work also mentions LANL’s Geothermal Design Tool “GeoDT” which is another approved open source code that has been released by LANL. This work was performed for DOE Geothermal Technologies Office (DE-EE-3.1.8.1). The host for the released video is intended to be YouTube or a suitable perpetual data repository such as GDR.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

New constraints on sodium production in globular clusters from the Na 23 ( He 3 , d ) Mg 24 reaction

The star-to-star anticorrelation of sodium and oxygen is a defining feature of globular clusters, but, to date, the astrophysical site responsible for this unique chemical signature remains unknown. Sodium enrichment within these clusters depends sensitively on reaction rate of the sodium destroying reactions 23 Na(p, γ) and 23 Na(p,α). In this paper, we report the results of a 23 Na( 3 He,d) 24 Mg transfer reaction carried out at Triangle Universities Nuclear Laboratory using a 21 MeV 3 He beam. Astrophysically relevant states in 24 Mg between 11 < E x < 12 MeV were studied using high-resolution magnetic spectroscopy, thereby allowing the extraction of excitation energies and spectroscopic factors. Bayesian methods are combined with the distorted wave Born approximation to assign statistically meaningful uncertainties to the extracted spectroscopic factors. For the first time, these uncertainties are propagated through to the estimation of proton partial widths. Our experimental data are used to calculate the reaction rate. The impact of the new rates are investigated using asymptotic giant branch star models. Furthermore, it is found that while the astrophysical conditions still dominate the total uncertainty, intramodel variations on sodium production from the 23 Na(p, γ) and 23 Na(p,α) reaction channels are a lingering source of uncertainty.

20 ≤ A ≤ 38↗
DOE OSTI2022

Case studies in bias reduction and inference for electronic health record data with selection bias and phenotype misclassification

Electronic health records (EHR) are not designed for population‐based research, but they provide easy and quick access to longitudinal health information for a large number of individuals. Many statistical methods have been proposed to account for selection bias, missing data, phenotyping errors, or other problems that arise in EHR data analysis. However, addressing multiple sources of bias simultaneously is challenging. We developed a methodological framework (R package, SAMBA ) for jointly handling both selection bias and phenotype misclassification in the EHR setting that leverages external data sources. These methods assume factors related to selection and misclassification are fully observed, but these factors may be poorly understood and partially observed in practice. As a follow‐up to the methodological work, we demonstrate how to apply these methods for two real‐world case studies, and we evaluate their performance. In both examples, we use individual patient‐level data collected through the University of Michigan Health System and various external population‐based data sources. In case study (a), we explore the impact of these methods on estimated associations between gender and cancer diagnosis. In case study (b), we compare corrected associations between previously identified genetic loci and age‐related macular degeneration with gold standard external summary estimates. These case studies illustrate how to utilize diverse auxiliary information to achieve less biased inference in EHR‐based research.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Liquid solution centrifugation for safe, scalable, and efficient isotope separation

A general method of separating isotopes by centrifuging dissolved chemical compounds in a liquid is introduced. This technique can be applied to almost all elements and leads to large separation factors. The method has been demonstrated in several isotopic systems including Ca, Mo, O, and Li with single-stage selectivities of 1.046 to 1.067 per neutron mass difference (e.g., 1.43 in 40 Ca/ 48 Ca), which are beyond the capabilities of various conventional methods. Equations are derived to model the process, and the results agree with those of the experiments. The scalability of the technique has been demonstrated by a three-stage enrichment of 48 Ca with a total 40 Ca/ 48 Ca selectivity of 2.43, and the scalability is more broadly supported through analogies to gas centrifuge, whereby countercurrent centrifugation can further multiply the separation factor by 5 to 10 times per stage in a continuous process. Optimal centrifuge conditions and solutions can achieve both high-throughput and highly efficient isotope separation.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2024

Acceleration of Solvation Free Energy Calculation via Thermodynamic Integration Coupled with Gaussian Process Regression and Improved Gelman–Rubin Convergence Diagnostics

The determination of the solvation free energy of ions and molecules holds profound importance across a spectrum of applications spanning chemistry, biology, energy storage, and the environment. Molecular dynamics simulations are powerful tools for computing this critical parameter. Nevertheless, the accurate and efficient calculation of the solvation free energy becomes a formidable endeavor when dealing with complex systems characterized by potent Coulombic interactions and sluggish ion dynamics and, consequently, slow transition across various metastable states. Here, in the present study, we expose limitations stemming from the conventional calculation of the statistical inefficiency g in the thermodynamic integration method, a factor that can hinder the determination of convergence of the solvation free energy and its associated uncertainty. Instead, we propose a robust scheme based on Gelman–Rubin convergence diagnostics. We leverage this improved estimation of uncertainties to introduce an innovative accelerated thermodynamic integration method based on the Gaussian Process regression. This methodology is applied to the calculation of the solvation free energy of trivalent rare-earth elements immersed in ionic liquids, a scenario in which the aforementioned challenges render standard approaches ineffective. The proposed method proves to be effective in computing solvation free energy in situations where traditional thermodynamic integration methods fall short.

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

New Methods for Predicting Non-Born-Oppenheimer Chemistry

Current methods for modeling non-adiabatic molecular dynamics face fundamental limitations when treating geometric phase effects: quantum mechanical phenomena where nuclear wavepackets acquire phase shifts when encircling conical intersections. Existing approaches either neglect these effects entirely or rely on potential energy surfaces arising from the Born-Oppenheimer approximation, which introduce artificial singularities and can overestimate geometric phase contributions. This project developed a new theoretical framework based on exact factorization (XF) methods to overcome these limitations. We derived mathematical formulations for hybrid quantum-classical XF dynamics that selectively treat critical nuclear degrees of freedom quantum mechanically while propagating others classically. This approach addresses the computational intractability that has previously limited exact methods to toy systems. Key innovations include a new approach to systematically identifying nuclear coordinates requiring quantum treatment, as well as novel implementation strategies that interface with existing quantum chemistry codes. The project also developed a proof-of-concept code for treating Jahn-Teller systems and creation of educational materials on non-adiabatic dynamics geared at the graduate level. The theoretical framework developed will enable future systematically improvable calculations of nuclear quantum effects in realistic molecular systems, filling a critical gap in non-adiabatic dynamics methods. This foundation supports future development of predictive tools for designing energy-relevant photochemical processes where quantum coherence effects may be exploited to control reaction outcomes.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Controlled Synthesis of SnO 2 Nanocrystals with Tunable Band Gaps

Tin(IV) oxide nanocrystals (SnO 2 NCs) have significant potential in various applications, with their performance closely related to their band gap. The band gap is influenced by the size and shape of the NCs, which can be precisely controlled by adjusting reaction conditions. In this study, we present deliberately designed synthesis protocols to produce high-quality SnO 2 NCs with tunable band gaps using different methods. Key factors affecting the synthesis include control of the oxidizing agent, reaction temperature, solvent selection, and reaction time optimization. The resulting NCs were characterized by using TEM, XRD, XPS, and optical spectroscopy. Notably, SnO 2 NCs synthesized by controlling the oxidizing agent (air injection) in a hot organic solution were smaller in size and exhibited abundant oxygen vacancies. In contrast, extending the reaction time or using ethanol as a solvent in hydrothermal systems facilitated larger spherical or rod-like SnO 2 NCs with fewer oxygen vacancies. Further analysis of the band gap and valence band maximum energy revealed that the abundant vacancies in SnO 2 NCs synthesized with the air-controlled hot organic solution method resulted in a narrower band gap and an upshifted valence band. These synthetic strategies illustrate the potential for deliberately designing SnO 2 NCs with optimized electronic structures for various applications.

band gap↗
DOE OSTI2022

Demonstration and Evaluation of the Human-Technology Integration Guidance for Plant Modernization

The significance of nuclear power in its role producing carbon-free electricity to the U.S. cannot be overstated. However, with changes in the energy market coupled with changes in incentives given to certain resources like solar and wind, the operating and maintenance costs for these sources have seen a significant reduction, which has consequently negatively impacted the economic viability of the existing U.S. nuclear power plant fleet. Digital instrumentation and control (I&C) and control room modernization is a major critical work domain to reduce operating and maintenance costs. Existing nuclear power plants are commonly configured with mostly legacy analog I&C as well as isolated pockets of digital I&C (a plant process computer, digital recorders, etc.). One challenge with this analog I&C is that replacement parts are becoming prohibitively more expensive and difficult to obtain. Moreover, a significant challenge with the existing analog I&C is that the way in which plants are currently operated and maintained is no longer competitive with other electricity generating sources, like natural gas, where advanced digital I&C technologies are commonplace. This gap between the waynuclear power plants are operated compared to other electricity generating sources significantly challenges the economic viability of the nuclear industry. Indeed, digital I&C systems can fundamentally change the way the plant is operated (i.e., the concept of operation). The introduction of digital I&C technologies offers a wealth of benefits to the nuclear industry. However, it is important to emphasize that, to realize these benefits, a careful understanding of how to integrate technology in a collaborative way that leverages the capabilities of people and technologies is necessary. Human-technology integration applies human factors engineering methods and tools to ensure the safe and reliable use of these technologies while ensuring that the inherent features of the technologies that provide economic value are not missed. The scope of this work documents the demonstration of the recently developed human and technology integration methodology, as applied to developing a new vision and concept of operations for a major U.S. nuclear power plant fleet. This report focuses on the methodological aspects of developing a vision and concept of operations. This report shares the tools, activities, and lessons learned during a modernization currently underway for industry as a whole to consider when planning any significant digital modification and developing a new vision and concept of operations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Investigating the universality of five-point QCD scattering amplitudes at high energy

We investigate 2 → 3 QCD scattering amplitudes in multi-Regge kinematics, i.e. where the final partons are strongly ordered in rapidity. In this regime amplitudes exhibit intriguing factorisation properties which can be understood in terms of effective degrees of freedom called reggeons. Working within the Balitsky/JIMWLK framework, we predict these amplitudes for the first time to next-to-next-to-leading logarithmic order, and compare against the limit of QCD scattering amplitudes in full colour and kinematics. We find that the latter can be described in terms of universal objects, and that the apparent non-universality arising at NNLL comes from well-defined and under-control contributions that we can predict. Thanks to this observation, we extract for the first time the universal vertex that controls the emission of the central-rapidity gluon, both in QCD and $\mathcal{N}$ = 4 super Yang-Mills.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Speciation mapping of the oxidation layer on aged uranium dioxide using scanning transmission x-ray spectromicroscopy

In this study, UO 2 was aged in humid air and prepared as a thin section using a focused ion beam (FIB) instrument. The specimen was measured using synchrotron radiation spectromicroscopy techniques at the Beamline 11.0.2 STXM end station of the Advanced Light Source (ALS). Non-negative matrix factorization (NMF) methods were used to identify and map three component x-ray absorption near-edge structure (XANES) spectra in the oxygen K-edge data, revealing a surface layer of U 4 O 9 with a thickness of 206 ± 21 nm, and the bulk of the sample remaining as UO 2 . Uranium N 4,5 -edge XANES spectromicroscopy supports these results. The diffusion-controlled parabolic rate constant for UO 2 oxidation to U 4 O 9 was calculated from the observed layer thickness and compared to literature values. Complementary transmission electron microscopy (TEM) was used to image the sample and identify the phases present in various regions, confirming the STXM results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Interpreting accelerated tests on perovskite modules using photooxidation of MAPbI 3 as an example

Solar panels (modules) based on metal halide perovskites are following a fast track to commercialization. Unlike more established solar cell materials, there are not yet decades-long field observations to increase consumer confidence. The "physics and chemistry of failure" approach is used in other industries and estimates product degradation based on laboratory accelerated tests integrated with an understanding of degradation mechanisms. This work uses that approach to quantify the relationship between accelerated tests and projected product behavior for metal halide perovskite modules. Degradation involving photooxidation of methylammonium lead iodide is used to illustrate the method. Acceleration factors in common accelerated tests are found to be low. Conclusions emphasize that the accelerated tests on photovoltaics should not be interpreted as equivalent across module types or as a green light for commercialization unless supported by the appropriate field data or physics and chemistry of failure analysis.

14 SOLAR ENERGY↗
DOE OSTI2023

Robust design of semi-automated clustering models for 4D-STEM datasets

Materials discovery and design require characterizing material structures at the nanometer and sub-nanometer scale. Four-Dimensional Scanning Transmission Electron Microscopy (4D-STEM) resolves the crystal structure of materials, but many 4D-STEM data analysis pipelines are not suited for the identification of anomalous and unexpected structures. This work introduces improvements to the iterative Non-Negative Matrix Factorization (NMF) method by implementing consensus clustering for ensemble learning. We evaluate the performance of models during parameter tuning and find that consensus clustering improves performance in all cases and is able to recover specific grains missed by the best performing model in the ensemble. The methods introduced in this work can be applied broadly to materials characterization datasets to aid in the design of new materials.

Bruefach, Alexandra (ORCID:0000000209323477)↗
DOE OSTI2023

Finite-temperature many-body perturbation theory for anharmonic vibrations: Recursions, algebraic reduction, second-quantized reduction, diagrammatic rules, linked-diagram theorem, finite-temperature self-consistent field, and general-order algorithm

A unified theory is presented for finite-temperature many-body perturbation expansions of the anharmonic vibrational contributions to thermodynamic functions, i.e., the free energy, internal energy, and entropy. The theory is diagrammatically size-consistent at any order, as ensured by the linked-diagram theorem proved in this study, and, thus, applicable to molecular gases and solids on an equal footing. It is also a basis-set-free formalism, just like its underlying Bose–Einstein theory, capable of summing anharmonic effects over an infinite number of states analytically. It is formulated by the Rayleigh–Schrödinger-style recursions, generating sum-over-states formulas for the perturbation series, which unambiguously converges at the finite-temperature vibrational full-configuration-interaction limits. Two strategies are introduced to reduce these sum-over-states formulas into compact sum-over-modes analytical formulas. One is a purely algebraic method that factorizes each many-mode thermal average into a product of one-mode thermal averages, which are then evaluated by the thermal Born–Huang rules. Canonical forms of these rules are proposed, dramatically expediting the reduction process. The other is finite-temperature normal-ordered second quantization, which is fully developed in this study, including a proof of thermal Wick’s theorem and the derivation of a normal-ordered vibrational Hamiltonian at finite temperature. The latter naturally defines a finite-temperature extension of size-extensive vibrational self-consistent field theory. These reduced formulas can be represented graphically as Feynman diagrams with resolvent lines, which include anomalous and renormalization diagrams. Two order-by-order and one general-order algorithms of computing these perturbation corrections are implemented and applied up to the eighth order. The results show no signs of Kohn–Luttinger-type nonconvergence.

74 ATOMIC AND MOLECULAR PHYSICS↗