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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The Transient Multi-Level method for Monte Carlo reactor statics calculations

The Transient Multi-Level (TML) method is applied to a time-dependent Monte Carlo transport solver to offload some of the computational burden of the expensive Monte Carlo solve to lower-order Coarse Mesh Finite Difference (CMFD) and Exact Point Kinetics Equations (EPKE) solvers via factorization of the neutron flux at the transport and CMFD levels using the Predictor Corrector Quasi-Static Method (PCQM). The Monte Carlo transient is solved by a modified fission source iteration scheme that introduces a single transient source bank. The method is implemented in the production-level Monte Carlo code, Shift, and verified with prescribed reactivity ramps from the two-dimensional version of the C5G7-TD reactor benchmark. The results show that, as compared to other quasi-static methods, the TML reduces the stochastic noise inherent to the transient Monte Carlo solver by factors of ~2 to 6 for various norm comparisons of the reactor power amplitude. Finally, the TML additionally reduces the number of Monte Carlo evaluations needed to simulate the transient, leading to roughly an order of magnitude improvement in CPU time relative to the standard PCQM for the problems tested.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Determination of latent dimensionality in international trade flow

Currently, high-dimensional data is ubiquitous in data science, which necessitates the development of techniques to decompose and interpret such multidimensional (aka tensor) datasets. Finding a low dimensional representation of the data, that is, its inherent structure, is one of the approaches that can serve to understand the dynamics of low dimensional latent features hidden in the data. Moreover, decomposition methods with non-negative constraints are shown to extract more insightful factors. Nonnegative RESCAL is one such technique, particularly well suited to analyze self-relational data, such as dynamic networks found in international trade flows. Particularly, non-negative RESCAL computes a low dimensional tensor representation by finding the latent space containing multiple modalities. Furthermore, estimating the dimensionality of this latent space is crucial for extracting meaningful latent features. Here, to determine the dimensionality of the latent space with non-negative RESCAL, we propose a latent dimension determination method which is based on clustering of the solutions of multiple realizations of non-negative RESCAL decompositions. We demonstrate the performance of our model selection method on synthetic data. We then apply our method to decompose a network of international trade flows data from International Monetary Fund and shows that with a correct latent dimension determination, the resulting features are able to capture relevant empirical facts from economic literature.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Suppressing the sample variance of DESI-like galaxy clustering with fast simulations

Ongoing and upcoming galaxy redshift surveys, such as the Dark Energy Spectroscopic Instrument (DESI) survey, will observe vast regions of sky and a wide range of redshifts. In order to model the observations and address various systematic uncertainties, N-body simulations are routinely adopted, however, the number of large simulations with sufficiently high mass resolution is usually limited by available computing time. Therefore, achieving a simulation volume with the effective statistical errors significantly smaller than those of the observations becomes prohibitively expensive. In this study, we apply the Convergence Acceleration by Regression and Pooling (CARPool) method to mitigate the sample variance of the DESI-like galaxy clustering in the AbacusSummit simulations, with the assistance of the quasi-N-body simulations FastPM. Based on the halo occupation distribution (HOD) models, we construct different FastPM galaxy catalogs, including the luminous red galaxies (LRGs), emission line galaxies (ELGs), and quasars, with their number densities and two-point clustering statistics well matched to those of AbacusSummit. We also employ the same initial conditions between AbacusSummit and FastPM to achieve high cross-correlation, as it is useful in effectively suppressing the variance. Our method of reducing noise in clustering is equivalent to performing a simulation with volume larger by a factor of 5 and 4 for LRGs and ELGs, respectively. We also mitigate the standard deviation of the LRG bispectrum with the triangular configurations k 2 = 2k 1 = 0.2 h Mpc -1 by a factor of 1.6. With smaller sample variance on galaxy clustering, we are able to constrain the baryon acoustic oscillations (BAO) scale parameters to higher precision. The CARPool method will be beneficial to better constrain the theoretical systematics of BAO, redshift space distortions (RSD) and primordial non-Gaussianity (NG).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A direct detection method of galaxy intrinsic ellipticity-gravitational shear correlation in non-linear regimes using self-calibration

Intrinsic alignment (IA) of galaxies is a challenging source of contamination in the Cosmic shear (GG) signals. The galaxy intrinsic ellipticity-gravitational shear (IG) correlation is generally the most dominant component of such contamination for cross-correlating redshift bins. One of the most effective techniques to mitigate such contamination is the self-calibration (SC) method which extracts the IG correlation and allows for its removal from the GG signal. In a photometric survey, the SC method first extracts the galaxy number density-galaxy intrinsic ellipticity (gI) correlation from the observed galaxy-galaxy lensing correlation using the redshift dependence of lens-source pairs. The IG correlation is computed through a scaling relation using the gI correlation and other lensing observables. The applicability of the SC method has so far been focused on the linear IA scales and the linear galaxy bias. We extend the SC method beyond the linear regime by modifying its scaling relation which can account for the non-linear galaxy bias model and various IA models. In this study, we provide a framework to detect the IG correlation for the redshift bins for source galaxies for the proposed year 1 survey of the Rubin Legacy Survey of Space and Time (LSST Y1). We tested the method for the tidal alignment and tidal torquing (TATT) model of IA and we found that the scaling relation is accurate within 10% and 20% for cross-correlating and auto-correlating redshift bins, respectively. Hence the suppression of IG contamination in observed GG correlation can be accomplished with a factor of 10 and 5, for cross-correlating and auto-correlating redshift bins, respectively. We tested the method's robustness and found that the suppression of IG contamination by a factor of 5 is still achievable for all combinations of cross-correlating bins even with the inclusion of a moderate amount of uncertainties on IA and bias parameters, respectively. We also make available, a branch of the code FAST-PT to provide gI correlations up to 1-loop order term used by the new SC method.

gravitational lensing↗
DOE OSTI2023

The Association of Radiation Exposure with Stable Chromosome Aberrations in Atomic Bomb Survivors Based on DS02R1 Dosimetry and FISH Methods

The frequency of stable chromosome aberrations (sCA) in lymphocytes is a recognized radiation biological dosimeter. Its analysis can provide insights into factors that affect individual susceptibility as well as into the adequacy of radiation dose estimates used in studies of atomic bomb survivors. We analyzed the relationship between atomic bomb radiation exposure using the most recent DS02R1 dose estimates and the frequency of sCA as determined by FISH in 1,868 atomic bomb survivors. We investigated factors that may affect the background sCA rate and the shape and magnitude of the dose response. As in previous analyses of sCA in atomic bomb survivors that were based on Giemsa staining methods and used older DS86 dose estimates, the relationship between radiation dose and sCA rate was significant (P < 0.0001) with a linear-quadratic relationship at lower doses that did not persist at higher doses. As before, age at the time of the bombing and type of radiation shielding were significant dose-effect modifiers (P < 0.0001), but in contrast the difference in dose response by city was not so pronounced (P = 0.026) with a city effect not evident at doses below 1.25Gy. Background sCA rate increased with age at the time of examination (P < 0.0001), but neither sex, city, nor smoking was significantly associated with background rate. Based on FISH methods and recent dosimetry, the relationship between radiation dose and sCA frequency is largely consistent with previous findings, although the lesser importance of city as an effect modifier may reflect better dosimetry as well as more reproducible scoring of sCA. The persisting difference in sCA dose response by shielding category points to remaining problems with the accuracy or precision of radiation dose estimates in some A-bomb survivors.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2020

Energy Efficient Thermoplastic Composite Manufacturing

The objective of the project was to establish an effective and affordable method to lay-up and consolidate/join large thermoplastic composite aerospace structure with cycle times measured in minutes rather than hours. Composite airplane designs have proven efficient and effective however future potential product production rates are challenging what the current systems can efficiently produce due to material lay-down constraints and extended thermal cycle times. The ability to lay-up then rapidly heat, consolidate, and cool large complex composite structures plus very accurately tool them (i.e. matching CTE of composite materials) along with very precise thermal control is a difficult challenge. Current systems such as autoclave processing of thermoset materials require long cycle times due to method of heating and the large associated thermal masses. These extended cycle times inhibit the ability to meet higher rate production scenarios due to the need for multiple sets of equipment and tools. Thermoplastic composite materials were used to facilitate more rapid cycle times via the elimination of a need for a cure dwell at temperature. Also, utilized induction heating along with smart susceptors to enable the quick cycle times needed while providing precise intrinsic thermal control. The intrinsic control of the induction heating process via the smart susceptors along with the use of laser assisted fiber placement were the key innovations areas developed. The ability to rapidly heat and then precisely control the temperature of the consolidation component was paramount to the success of this project. It is this unique processing attribute coupled with the thermoplastic material characteristics along with rapid lay-up methods that are novel and advantageous. While many requirements and influencing factors decide the materials and processes utilized for future airplane construction, the forecast of accelerated production rates and the recent performance successes of composites in airplane construction provide an opportunity for this processing technology to have significant influence. During the course of this project a number of key accomplishments were completed. They are as follows: The needed laser assisted fiber placement capabilities have been developed and validated on a component of meaningful scale and complexity. Subsequently, a thermoplastic preform of the part scale-up part design was laser assisted fiber placed and made availalbe for consolidation. An induction consolidation system was sized, designed, and fabricated resulting in the full capability to consolidate the scale-up component. The tool for induction consolidating of the scale-up component was designed and fabricated. An aluminum pressure bladder was welded and formed for use as the pressurization membrane to perform the needed application of pressure on the part at temperature to achieve consolidation. The induction consolidation system consisting of the restraint, induction power supply, induction tool with integrated induction coils and smart susceptor liners, along with the aluminum bladder were successfully used to consolidate a large thermoplastic skin. The capability of accomplishing rapid heat-up rates and cool down rates for consolidation of large thermoplastic skins with precise thermal control and even pressure application at the consolidation temperature was validated. The general conclusion is that the use of induction consolidation with smart susceptors is an effective method to rapidly and efficiently consolidate large thermoplastic composite components and enable affordable high rate manufacturing of composites. More specific conclusions generated by this project are as follows: A restraint type induction consolidation system provides an affordable option over the standard press when using the induction consolidation processing method. Aluminum bladders can be formed to the complex shapes needed and supply uniform pressure for consolidation and co-consolidation of thermoplastic composite materials. Straightforward modifications can be made to existing fiber placement robotic systems to enable their use for thermoplastic lay-up applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Overview of ASDEX upgrade results in view of ITER and DEMO

Experiments on ASDEX Upgrade (AUG) in 2021 and 2022 have addressed a number of critical issues for ITER and EU DEMO. A major objective of the AUG programme is to shed light on the underlying physics of confinement, stability, and plasma exhaust in order to allow reliable extrapolation of results obtained on present day machines to these reactor-grade devices. Concerning pedestal physics, the mitigation of edge localised modes (ELMs) using resonant magnetic perturbations (RMPs) was found to be consistent with a reduction of the linear peeling-ballooning stability threshold due to the helical deformation of the plasma. Conversely, ELM suppression by RMPs is ascribed to an increased pedestal transport that keeps the plasma away from this boundary. Candidates for this increased transport are locally enhanced turbulence and a locked magnetic island in the pedestal. The enhanced D-alpha (EDA) and quasi-continuous exhaust (QCE) regimes have been established as promising ELM-free scenarios. Here, the pressure gradient at the foot of the H-mode pedestal is reduced by a quasi-coherent mode, consistent with violation of the high-n ballooning mode stability limit there. This is suggestive that the EDA and QCE regimes have a common underlying physics origin. In the area of transport physics, full radius models for both L- and H-modes have been developed. These models predict energy confinement in AUG better than the commonly used global scaling laws, representing a large step towards the goal of predictive capability. A new momentum transport analysis framework has been developed that provides access to the intrinsic torque in the plasma core. In the field of exhaust, the X-Point Radiator (XPR), a cold and dense plasma region on closed flux surfaces close to the X-point, was described by an analytical model that provides an understanding of its formation as well as its stability, i.e., the conditions under which it transitions into a deleterious MARFE with the potential to result in a disruptive termination. With the XPR close to the divertor target, a new detached divertor concept, the compact radiative divertor, was developed. Here, the exhaust power is radiated before reaching the target, allowing close proximity of the X-point to the target. No limitations by the shallow field line angle due to the large flux expansion were observed, and sufficient compression of neutral density was demonstrated. With respect to the pumping of non-recycling impurities, the divertor enrichment was found to mainly depend on the ionisation energy of the impurity under consideration. In the area of MHD physics, analysis of the hot plasma core motion in sawtooth crashes showed good agreement with nonlinear 2-fluid simulations. This indicates that the fast reconnection observed in these events is adequately described including the pressure gradient and the electron inertia in the parallel Ohm’s law. Concerning disruption physics, a shattered pellet injection system was installed in collaboration with the ITER International Organisation. Thanks to the ability to vary the shard size distribution independently of the injection velocity, as well as its impurity admixture, it was possible to tailor the current quench rate, which is an important requirement for future large devices such as ITER. Progress was also made modelling the force reduction of VDEs induced by massive gas injection on AUG. The H-mode density limit was characterised in terms of safe operational space with a newly developed active feedback control method that allowed the stability boundary to be probed several times within a single discharge without inducing a disruptive termination. Regarding integrated operation scenarios, the role of density peaking in the confinement of the ITER baseline scenario (high plasma current) was clarified. The usual energy confinement scaling ITER98(p,y) does not capture this effect, but the more recent H20 scaling does, highlighting again the importance of developing adequate physics based models. Advanced tokamak scenarios, aiming at large non-inductive current fraction due to non-standard profiles of the safety factor in combination with high normalised plasma pressure were studied with a focus on their access conditions. A method to guide the approach of the targeted safety factor profiles was developed, and the conditions for achieving good confinement were clarified. Based on this, two types of advanced scenarios (‘hybrid’ and ‘elevated’

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Structural Uncertainty in the Sensitivity of Urban Temperatures to Anthropogenic Heat Flux

One key source of uncertainty for weather and climate models is structural uncertainty arising from the fact that these models must simplify or approximate complex physical, chemical, and biological processes that occur in the real world. However, structural uncertainty is rarely examined in the context of simulated effects of anthropogenic heat flux in cities. Using the Weather Research and Forecasting (WRF) model coupled with a single-layer urban canopy model, it is found that the sensitivity of urban canopy air temperature to anthropogenic heat flux can differ by an order of magnitude depending on how anthropogenic heat flux is released to the urban environment. Moreover, varying model structures through changing the treatment of roof-air interaction and the parameterization of convective heat transfer between the canopy air and the atmosphere can affect the sensitivity of urban canopy air temperature by a factor of 4. Urban surface temperature and 2-m air temperature are less sensitive to the methods of anthropogenic heat flux release and the examined model structural variants than urban canopy air temperature, but their sensitivities to anthropogenic heat flux can still vary by as much as a factor of 4 for surface temperature and 2 for 2-m air temperature. Our study recommends using temperature sensitivity instead of temperature response to understand how various physical processes (and their representations in numerical models) modulate the simulated effects of anthropogenic heat flux.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Identifying spatiotemporal patterns in opioid vulnerability: investigating the links between disability, prescription opioids and opioid-related mortality

Background: The opioid crisis remains one of the most daunting and complex public health problems in the United States. This study investigates the national epidemic by analyzing vulnerability profiles of three key factors: opioid-related mortality rates, opioid prescription dispensing rates, and disability rank ordered rates. Methods: This study utilizes county level data, spanning the years 2014 through 2020, on the rates of opioid-related mortality, opioid prescription dispensing, and disability. To successfully estimate and predict trends in these opioid-related factors, we augment the Kalman Filter with a novel spatial component. To define opioid vulnerability profiles, we create heat maps of our filter’s predicted rates across the nation’s counties and identify the hotspots. In this context, hotspots are defined on a year-by-year basis as counties with rates in the top 5% nationally. Results: Our spatial Kalman filter demonstrates strong predictive performance. From 2014 to 2018, these predictions highlight consistent spatiotemporal patterns across all three factors, with Appalachia distinguished as the nation’s most vulnerable region. Starting in 2019 however, the dispensing rate profiles undergo a dramatic and chaotic shift. Conclusions: The initial primary drivers of opioid abuse in the Appalachian region were likely prescription opioids; however, it now appears that abuse is sustained by illegal drugs. Additionally, we find that the disabled subpopulation may be more at risk of opioid-related mortality than the general population. Public health initiatives must extend beyond controlling prescription practices to address the transition to and impact of illicit drug use.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Defect-specific strength factors and superposition model for predicting strengthening of ion irradiated Fe18Cr alloy

Here, a high-purity binary alloy Fe18Cr was subjected to heavy-ion irradiation in order to provide improved understanding of the irradiation effect on radiation hardening associated with dislocation loop and network formation, α’ precipitation. The specimens were irradiated with 8 MeV Fe ions (~2 μm ion range) to midrange doses of 0.37 and 3.7 dpa at 300, 350 and 450 °C using dose rates of ~10 –5 –10 –3 dpa/s. Nanoindentation testing was performed to extract the bulk equivalent hardness at low depth from these ion irradiated specimens. High-quality transmission electron microscopy (TEM) images were acquired by utilizing a flash electropolishing method on the Focused Ion Beam (FIB) prepared liftouts. An accurate calculation of the strength factor was provided based on detailed TEM characterization of the irradiated microstructures. A newly refined hardening superposition method was applied to combine the strengthening components from each microstructure. The good agreement between microstructure predicted strength and measured strength further demonstrated the fidelity to use hardening model to quantify the mechanical properties of ion irradiated materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Toward Establishing Uniqueness of Experimentally Determined Transference Numbers

The passage of current through a battery results in the development of concentration gradients in the electrolytic phase. For a fully characterized binary electrolyte, where the conductivity, salt diffusion coefficient, cation transference number, and the thermodynamic factor are known, concentration and potential gradients in the electrolytic phase can be modeled using Newman’s concentrated solution theory. We report two methods for measuring the transference number: the standard method based on electrochemical measurements ( t + , echem 0 ) and electrophoretic NMR ( t + , eNMR 0 ). The electrochemical approach requires combining measurements from multiple experiments; the equations used to determine the cation transference number and the thermodynamic factor are coupled, nonlinear algebraic equations. In the electrophoretic-NMR-based approach, however, the equations used to determine the cation transference number and the thermodynamic factor are decoupled. We find for a liquid electrolyte comprised of a lithium salt dissolved in tetraglyme, the values of the transference numbers obtained by these two methods are distinct. For example, at 30 °C, t + , echem 0 = −1.02 ± 1.11 and t + , eNMR 0 = 0.25 ± 0.04. The corresponding thermodynamic factors are also different. While the magnitude of the predicted concentration gradients based on the two sets of parameters are different, the predicted current-voltage relationships are similar.

Hickson, Darby T. (ORCID:0000000251339755)↗
DOE OSTI2022

Complete characterization of a lithium battery electrolyte using a combination of electrophoretic NMR and electrochemical methods

Improving transport properties of electrolytes is important for developing lithium-ion batteries for future energy storage applications. In Newman's concentrated solution theory, electrolytes are characterized by three transport parameters, conductivity, diffusion coefficient, and transference number, in addition to the thermodynamic factor. In this work, these parameters are all determined for an exemplar liquid electrolyte, lithium bis(trifluoromethanesulfonyl)imide mixed in tetraethylene glycol dimethyl ether, using electrochemical methods. The intrinsic coupling between parameters obtained by electrochemical methods results in large error bars in the transference number that obscure the transport behavior of the electrolyte. Here, we use electrophoretic NMR (eNMR) to measure the electric-field-induced ion and solvent velocities to obtain the transference number directly, which enables determination of the thermodynamic factor with greater certainty. In conclusion, our work indicates that the combination of eNMR and electrochemical methods provides a robust approach for complete characterization of battery electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2020

Comparative Analysis of Thrombin Calibration Algorithms and Correction for Thrombin-α2macroglobulin Activity

Background: The thrombin generation (TG) test is useful for characterizing global hemostasis potential, but fluorescence substrate artifacts, such as thrombin-α2macroglobulin (T-α2MG) signal, inner filter effect (IFE), substrate consumption, and calibration algorithms have been suggested as sources of intra- and inter-laboratory variance, which may limit its clinical utility. Methods: Effects of internal vs. external normalization, IFE and T-α2MG on TG curves in normal plasma supplemented with coagulation factors, thrombomodulin, and tissue factor were studied using the Calibrated Automated Thrombinography (CAT; Diagnostica Stago, Parsippany, NJ, USA) and in-house software. Results: The various calibration methods demonstrated no significant difference in producing TG curves, nor increased the robustness of the TG assay. Several TG parameters, including thrombin peak height (TPH), produced from internal linear calibration did not differ significantly from uncalibrated TG parameters. Further, TPH values from internal linear and nonlinear calibration with or without T-α2MG correction correlated well with TPH from external calibration. Higher coefficients of variation (CVs) for TPH values were observed in both platelet-free and platelet-rich plasma with added thrombomodulin. Conclusions: Our work suggests minimal differences between distinct computational approaches toward calibrating and correcting fluorescence signals into TG levels, with most samples returning similar or equivalent TPH results.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

A new calibration method for charm jet identification validated with proton-proton collision events at $\sqrt{s}$ =13 TeV

Many measurements at the LHC require efficient identification of heavy-flavour jets, i.e. jets originating from bottom (b) or charm (c) quarks. An overview of the algorithms used to identify c jets is described and a novel method to calibrate them is presented. This new method adjusts the entire distributions of the outputs obtained when the algorithms are applied to jets of different flavours. It is based on an iterative approach exploiting three distinct control regions that are enriched with either b jets, c jets, or light-flavour and gluon jets. Results are presented in the form of correction factors evaluated using proton-proton collision data with an integrated luminosity of 41.5 fb$^{-1}$ at √s = 13 TeV, collected by the CMS experiment in 2017. The closure of the method is tested by applying the measured correction factors on simulated data sets and checking the agreement between the adjusted simulation and collision data. Furthermore, a validation is performed by testing the method on pseudodata, which emulate various mismodelling conditions. The calibrated results enable the use of the full distributions of heavy-flavour identification algorithm outputs, e.g. as inputs to machine-learning models. Thus, they are expected to increase the sensitivity of future physics analyses.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Pass-efficient methods for compression of high-dimensional turbulent flow data

The future of high-performance computing, specifically on future Exascale computers, will presumably see memory capacity and bandwidth fail to keep pace with data generated, for instance, from massively parallel partial differential equation (PDE) systems. Current strategies proposed to address this bottleneck entail the omission of large fractions of data, as well as the incorporation of in situ compression algorithms to avoid overuse of memory. To ensure that post-processing operations are successful, this must be done in a way that a sufficiently accurate representation of the solution is stored. Moreover, in situations where the input/output system becomes a bottleneck in analysis, visualization, etc., or the execution of the PDE solver is expensive, the number of passes made over the data must be minimized. In the interest of addressing this problem, this work focuses on the utility of pass-efficient, parallelizable, low-rank, matrix decomposition methods in compressing high-dimensional simulation data from turbulent flows. Additionally, a particular emphasis is placed on using coarse representation of the data – compatible with the PDE discretization grid – to accelerate the construction of the low-rank factorization. This includes the presentation of a novel single-pass matrix decomposition algorithm for computing the so-called interpolative decomposition. The methods are described extensively and numerical experiments on two turbulent channel flow data are performed. In the first (unladen) channel flow case, compression factors exceeding 400 are achieved while maintaining accuracy with respect to first- and second-order flow statistics. In the particle-laden case, compression factors of 100 are achieved and the compressed data is used to recover particle velocities. These results show that these compression methods can enable efficient computation of various quantities of interest in both the carrier and disperse phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Robust Automatic EXAFS First-Shell Fits

Extended X-ray absorption fine structure (EXAFS) is a widely used technique for atomic structure determination. Fourier transformation connects EXAFS in k space and R space. However, determining the appropriate k-range for the transformation can be challenging, but critical for the first-shell fit. In this study, we present an automatic method to determine the k-range using the Larch package and a Python program. The first step is to estimate spectral noise across a series of k-ranges with a fixed minimum value and identify the optimal maximum value in the k-range (k max ). The k max is determined by an empirical noise threshold that marks the point where the noise level in the Fourier transformed spectrum changes dramatically. Using the obtained k max value, the first shell is modeled to determine the minimum k value (k min ) by optimizing the background function through alignment of the spectrum with theory. The optimal k min corresponds to the point of the minimum R-factor, which quantifies the difference between the experimental and fitted spectrum. Our method was tested on various typical datasets and yielded suitable k-ranges for Fourier transformation and accurate first-shell fits. This approach helps avoid unreliable, irreproducible data analysis, especially for noisy data from diluted samples, and enables robust automatic first-shell EXAFS fitting.

EXAFS analysis↗
DOE OSTI2020

Enhancing Oxygen Reduction Activity of Pt-based Electrocatalysts: From Theoretical Mechanisms to Practical Methods

Pt-based electrocatalysts are considered as one of the most promising choices to facilitate the oxygen reduction reaction (ORR), and the key factor enabling their success is to reduce the required amount of platinum. Herein, we focus on illuminating both the theoretical mechanisms which enable enhanced and sustained ORR activity and the practical methods to achieve them in catalysts. The various multi-step pathways of ORR are firstly reviewed and the rate-determining steps based on the reaction intermediates and their binding energies are analyzed. We then explain the critical aspects of Pt-based electrocatalysts to tune oxygen reduction properties from the viewpoints of active sites exposure and altering the surface electronic structure, and further summarize representative research progress towards practically achieving these activity enhancements with a focus on platinum size reduction, shape control and core Pt elimination methods. Finally, we outline the remaining challenges and provide our perspectives with regard to further enhancing their activities.

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

Multiplexing and Demultiplexing Signals for Radiography Application Using the Discrete Fourier Transform

Our goal is to develop an X-ray phase-contrast imaging system that can provide excellent soft tissue contrast of phase, attenuation, and small-angle scatter. We propose to replace the common system of G0, G1, and G2 gradings with a biprism array to replace the G1 grading and introduce a novel X-ray tube designed to replace the motion of the phase stepping grading G2. The proposed X-ray tube uses temporal multiplexing to provide simultaneous virtual “electronic phase stepping.” In this work the discrete Fourier transform is used to separate from the composite measurement individual X-ray phase contrast measurements sampled at different frequencies. The method performs a discrete Fourier transform of a composite refence sequence to obtain using the frequency amplitudes calibration factors needed to extract the X-ray phase contrast measurement amplitudes from the composite image. The composite reference sequence is the sum of the individual sequences, at different frequencies, with amplitudes of one. The method takes the discrete Fourier transform of this composite reference sequence; whereby, the amplitude of each frequency component is compared with the total sum of its stand-alone sequence amplitude. A calibration factor is determined so that the amplitude of this composite reference frequency times the calibration factor must equal the total sum of the sequence amplitude—the zero-frequency amplitude of the discrete Fourier transform of its stand-alone sequence. To demultiplex the composite measured signal these calibration factors are multiplied by the amplitudes of the frequency components of the discrete Fourier transform of the composite X-phase-contrast measurement to obtain the amplitude of each frequency encoded measurement. Using these calibration factors, we demonstrate with the discrete Fourier transform in Mathematica the extraction of individual images from a composite image that one would expect obtaining from our proposed new X-ray phase contrast imaging system. We then demonstrate as an example how using images from X-ray phase contrast data one can calculate phase, attenuation and the dark field images using grading phase step data supplied to use from Microworks, GmbH in Karlsruhe, Germany.

42 ENGINEERING↗