Search NASA⌕ Search

SEARCH · Search NASA

Results for “Approximation”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Characterization and thermometry of dissipatively stabilized steady states

In this work we study the properties of dissipatively stabilized steady states of noisy quantum algorithms, exploring the extent to which they can be well approximated as thermal distributions, and proposing methods to extract the effective temperature T. We study an algorithm called the relaxational quantum eigensolver (RQE), which is one of a family of algorithms that attempt to find ground states and balance error in noisy quantum devices. In RQE, we weakly couple a second register of auxiliary ‘shadow’ qubits to the primary system in Trotterized evolution, thus engineering an approximate zero-temperature bath by periodically resetting the auxiliary qubits during the algorithm’s runtime. Balancing the infinite temperature bath of random gate error, RQE returns states with an average energy equal to a constant fraction of the ground state. We probe the steady states of this algorithm for a range of base error rates, using several methods for estimating both T and deviations from thermal behavior. In particular, we both confirm that the steady states of these systems are often well-approximated by thermal distributions, and show that the same resources used for cooling can be adopted for thermometry, yielding a fairly reliable measure of the temperature. These methods could be readily implemented in near-term quantum hardware, and for stabilizing and probing Hamiltonians where simulating approximate thermal states is hard for classical computers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Free-energy orbital-free density functional theory: recent developments, perspective, and outlook

By summarizing the constraint-based development of orbital-free free-energy density functional approximations, we provide a perspective on progress over the last 15 years, the limitations of existing functionals, and the challenges awaiting resolution. We outline the chronology of the development of noninteracting and exchange-correlation free-energy orbital-free functionals and summarize the theoretical basis of existing local density approximation, second-order approximation, generalized gradient approximation (GGA), and meta-GGAs. We discuss limitations and challenges such as problems with thermodynamic derivatives, free-energy nonadditivity and the closely related issue of all-electron versus valence-only local pseudo-potential performance.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Small-𝑥 asymptotics of the leading-twist flavor-singlet quark TMDs

In this paper, we investigate the small-𝑥 behavior of the flavor-singlet, leading-twist quark transverse-momentum-dependent parton distribution functions (TMDs) using the light-cone operator treatment. This formalism allows us to express TMD operators at small 𝑥 in terms of polarized dipole amplitudes, enabling a systematic approach to their small-𝑥 evolution. We derive the evolution equations for these TMDs and solve them within the large-𝑁 𝑐 approximation under the linearized, double-logarithmic approximation, where 𝑁 𝑐 represents the number of quark colors. Expanding on previous work on unpolarized and helicity TMDs, we present the small-𝑥 asymptotics for a comprehensive set of TMDs, including the Sivers function, helicity worm-gear, transversity, pretzelosity, Boer-Mulders, and transversity worm-gear distributions. Our results provide a complete picture of the small-𝑥 asymptotic behavior for all leading-twist flavor-singlet quark TMDs. We also discuss the implications of our findings for phenomenological applications and outline potential avenues for further research, particularly in understanding nonlinear effects and extending beyond the double-logarithmic approximation and large-𝑁 𝑐 approximations.

Adamiak, Daniel [Thomas Jefferson National Acceler↗
DOE OSTI2025

A High-Speed, High-Resolution Transition Edge Sensor Spectrometer for Soft X-Rays at the Advanced Photon Source

This project explores the design and development of a transition edge sensor (TES) spectrometer for resonant soft X-ray scattering (RSXS) measurements developed in collaboration between Argonne National Laboratory (ANL) and the National Institute of Standards and Technology (NIST). Soft X-ray scattering is a powerful technique for studying the electronic and magnetic properties of materials on a microscopic level. However, the lack of high-performance soft X-ray spectrometers has limited the potential of this technique. TES spectrometers have the potential to overcome these limitations due to their high energy resolution, high efficiency, and broad energy range. This project aims to optimize the design of a TES spectrometer for RSXS measurements and more generally soft X-ray spectroscopy at the Advanced Photon Source (APS) 29-ID, leading to improved understanding of advanced materials. We will present a detailed description of the instrument design and implementation. The spectrometer consists of a large array of approximately 250 high-speed and high-resolution pixels. The pixels have saturation energies of approximately 1 keV, sub-ms pulse duration and energy resolution of approximately 1 eV. The array is read out using microwave multiplexing chips with MHz bandwidth per channel, enabling efficient data throughput. To facilitate measurement of samples in situ under ultra-high vacuum conditions at the beamline, the spectrometer is integrated with an approximately 1 m long snout.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Efficient Neural Network Approaches for Conditional Optimal Transport with Applications in Bayesian Inference

In this work, we present two neural network approaches that approximate the solutions of static and dynamic conditional optimal transport (COT) problems. Both approaches enable conditional sampling and conditional density estimation, which are core tasks in Bayesian inference—particularly in the simulation-based (“likelihood-free”) setting. Our methods represent the target conditional distribution as a transformation of a tractable reference distribution. Obtaining such a transformation, chosen here to be an approximation of the COT map, is computationally challenging even in moderate dimensions. To improve scalability, our numerical algorithms use neural networks to parameterize candidate maps and further exploit the structure of the COT problem. Our static approach approximates the map as the gradient of a partially input convex neural network. It uses a novel numerical implementation to increase computational efficiency compared to state-of-the-art alternatives. Our dynamic approach approximates the conditional optimal transport via the flow map of a regularized neural ODE; compared to the static approach, it is slower to train but offers more modeling choices and can lead to faster sampling. We demonstrate both algorithms numerically, comparing them with competing state-of-the-art approaches, using benchmark datasets and simulation-based Bayesian inverse problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Large $$N_c$$ QCD phase diagram at $$\mu _B=0$$

Abstract Lattice studies suggest that at zero baryon chemical potential and increasing temperature there are three characteristic regimes in QCD that are connected by smooth analytical crossovers: a hadron gas regime at$$T < T_{ch}\sim 155$$ T < T ch ∼ 155 MeV, an intermediate regime, called stringy fluid, at$$T_{ch}< T < \sim 3 T_{ch}$$ T ch < T < ∼ 3 T ch , and a quark-gluon plasma regime at higher temperatures. These regimes have been interpreted to reflect different approximate symmetries and effective degrees of freedom. In the hadron gas the effective degrees of freedom are hadrons and the approximate chiral symmetry of QCD is spontaneously broken. The intermediate regime has been interpreted as lacking spontaneous chiral symmetry breaking along with the emergence of new approximate symmetry, chiral spin symmetry, that is not a symmetry of the Dirac Lagrangian, but is a symmetry of the confining part of the QCD Lagrangian. While the high temperature regime is the usual quark-gluon plasma which is often considered to reflect “deconfinement” in some way. This paper explores the behavior of these regimes of QCD as the number of colors in the theory,$$N_c$$ N c , gets large. In the large$$N_c$$ N c limit the theory is center-symmetric and notions of confinement and deconfinement are unambiguous. The energy density is$$\mathcal{O}(N_c^0)$$ O ( N c 0 ) in the meson gas,$${{\mathcal {O}}}(N_c^1)$$ O ( N c 1 ) in the intermediate regime and$${{\mathcal {O}}}(N_c^2)$$ O ( N c 2 ) in the quark-gluon plasma regime. In the large$$N_c$$ N c limit these regimes may become distinct phases separated by first order phase transitions. The intermediate phase has the peculiar feature that glueballs should exist and have properties that are unchanged from what is seen in the vacuum (up to$$1/N_c $$ 1 / N c corrections), while the ordinary dilute gas of mesons with broken chiral symmetry disappears and approximate chiral spin symmetry should emerge.

Physics↗
DOE OSTI2024

nuclear-score-maximization v1.0

This software library presents efficient and multithreaded implementations of matrix low rank approximation via column selection in C++17 code. The algorithms are described in Fornace, Mark, and Michael Lindsey. "Column and row subset selection using nuclear scores: algorithms and theory for Nystro m approximation, CUR decomposition, and graph Laplacian reduction." arXiv preprint arXiv:2407.01698 (2024). The presented methods are by-and-large ver novel, have provable approximation guarantees, multiple use-cases, and exhibit higher quality approximations on a variety of studied examples.

Fornace, Mark↗
DOE OSTI2026

Data for Multisite Field Evaluation of Oil Accumulation and Agronomic Performance in Grain and Sweet Sorghums Engineered for Lipid Hyperaccumulation

Oil sorghum (OS) has been developed by engineering grain (TX430) and sweet (Ramada) genetic backgrounds to accumulate triacylglycerols (TAG) in vegetative tissues as an energy-dense feedstock for sustainable aviation fuel (SAF) and other biofuels. This study evaluated two TX430 OS lines (TxHO-2, TxHO-3) and two Ramada OS lines (RmHO-1, RmHO-2) alongside wild-type (WT) lines in NE and IL over 2 years (2023–2024) to quantify genotype × environment effects on agronomic performance and TAG accumulation. Across four environments, TX430 OS lines showed average TAG concentrations of 15.0 g kg−1 in leaves and 12.8 g kg−1 in stems, approximately 19-fold higher than WT. Ramada OS lines accumulated 26.1 g kg−1 in leaves and 12.3 g kg−1 in stems, approximately 25-fold and 13-fold increases over WT, respectively. OS lines in TX430 exhibited an 18% reduction in biomass (8.4 vs. 9.9 Mg ha−1 for WT), while Ramada OS lines had similar WT biomass (18.3 vs. 19.9 Mg ha−1 for WT). Among TX430 OS lines, TxHO-2 achieved the highest TAG yield (190 kg ha−1), while RmHO-1 led the Ramada lines (335 kg ha−1) due to higher biomass and similar TAG concentration. Enhanced TAG accumulation increased N, P, and K removal in TX430 lines but not in Ramada lines. Structural carbohydrate and ash concentration were unaffected. Overall, results confirm vegetative lipid accumulation as a viable strategy for high-biomass sorghum, supporting its potential as a dual-purpose feedstock for SAF. Future work should focus on minimizing biomass yield penalties and improving nutrient use efficiency in oil sorghum systems.

Agronomy↗
DOE OSTI2025

Destructive PIE and Safety Testing of Six AGR-5/6/7 Capsule 2 Compacts

This study evaluates fission product retention and particle failure mechanisms in AGR-5/6/7 Capsule 2 uranium carbide and uranium oxide (UCO) tristructural isotropic (TRISO) fuel under high-temperature gas reactor accident-relevant conditions using high-temperature safety tests and destructive postirradiation examination. Three Capsule 2 compacts were held isothermally at 1600°C for approximately 300 hours and one compact at 1800°C for approximately 300 hours; two additional compacts were examined in the as-irradiated state. Post-test deconsolidation–leach–burn–leach (DLBL) quantified nuclide inventories in matrix and particles. Individual particles were surveyed for radioisotope inventories, and microanalytical approaches resolved microstructural evolution and fission product distributions within the coating layers. At 1600°C, no krypton was detected above the minimum detectable limit, and cesium releases were far below a single particle equivalent, indicating the absence of full TRISO failure or SiC failures. Silver releases were limited and primarily reflected depleted postirradiation inventories, consistent with prior compact-level exams indicating substantial in-pile 110mAg loss. At 1800°C, cumulative 134Cs release of approximately 2.5 particle equivalents and delayed 85Kr totaling approximately 0.53 particle equivalents were consistent with one full TRISO failure and two SiC failures. Europium and strontium releases were roughly one order of magnitude higher than at 1600°C and comparable to AGR-1/AGR-2 high-temperature tests, with sustained late-hold rates indicating diffusion through intact coatings coupled with matrix depletion. Overall, AGR-5/6/7 Capsule 2 UCO fuel demonstrated fission product retention during safety testing consistent with prior AGR campaigns, while distinctive in-pile 110mAg depletion and measurable 1600°C europium loss motivate targeted follow-on studies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Destructive PIE and Safety Testing of Six AGR-5/6/7 Capsule 2 Compacts

This study evaluates fission product retention and particle failure mechanisms in AGR-5/6/7 Capsule 2 uranium carbide and uranium oxide (UCO) tristructural isotropic (TRISO) fuel under high-temperature gas reactor accident-relevant conditions using high-temperature safety tests and destructive post-irradiation examination. Three Capsule 2 compacts were held isothermally at 1600°C for approximately 300 hours and one compact at 1800°C for approximately 300 hours; two additional compacts were examined in the as-irradiated state. Post-test deconsolidation–leach–burn–leach (DLBL) quantified nuclide inventories in matrix and particles. Individual particles were surveyed for radioisotope inventories, and microanalytical approaches resolved microstructural evolution and fission product distributions within the coating layers. At 1600°C, no krypton was detected above the minimum detectable limit, and cesium releases were far below a single particle equivalent, indicating the absence of full TRISO failure or SiC failures. Silver releases were limited and primarily reflected depleted post-irradiation inventories, consistent with prior compact-level exams indicating substantial in-pile 110m Ag loss. At 1800°C, cumulative 134 Cs release of approximately 2.5 particle equivalents and delayed 85 Kr totaling approximately 0.53 particle equivalents were consistent with one full TRISO failure and two SiC failures. Europium and strontium releases were roughly one order of magnitude higher than at 1600°C and comparable to AGR-1/AGR-2 high-temperature tests, with sustained late-hold rates indicating diffusion through intact coatings coupled with matrix depletion. Overall, AGR-5/6/7 Capsule 2 UCO fuel demonstrated fission product retention during safety testing consistent with prior AGR campaigns, while distinctive in-pile 110m Ag depletion and measurable 1600°C europium loss motivate targeted follow-on studies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Using Separation-Enhanced Isotope Ratio Mass Spectrometry to Enable Increased Renewable Carbon Content in Transportation Fuels (CRADA 525)

Stable isotope ratio measurements of carbon atoms using isotope ratio mass spectrometry (IRMS) can be an effective tool for quantifying biogenic carbon in co-processed fuels, with results approaching the precision and accuracy of accelerator mass spectrometry (AMS). The lower cost of an IRMS may enable deployment to refineries, improving access and analysis turnaround times (≤2 hours), and, by extension, provide data that can allow process optimization to maximize renewable carbon in desired refinery products. This project explored the integration of chemical separation with IRMS analyses to enable highly detailed tracking of biogenic carbon into fuel product streams separated by boiling point range, chemical class, or specific compound. Forty-nine fuels and fuel components of fossil and biogenic origin, spanning gasoline and diesel boiling point ranges, were received from three refiners and were analyzed for their δ 13 C values via IRMS. Results spanned a 13 C range from ca. 10‰ to 44‰ and reflect materials derived from sustainable sources (e.g., C4 or C3 plants, animal-based pathways, syngas) or from fossil-derived fuels. Common ranges are approximately 18‰ to 9‰ and approximately 30‰ to 20‰ for C4 and C3 plants, respectively, and approximately 34‰ to 24‰ and approximately 70‰ to 33‰ for petroleum-derived fuels and methane, respectively. Fuel-like standards were developed and tested using direct-injection elemental analyzer (EA) IRMS for liquid fuels. This method was compared with the published methods, yielding statistically similar results. Four blend curve sets were produced ranging from 0% to 100% of a fuel containing biogenic carbon, focusing on 0% to 10% biogenic carbon. Linear fits were the most applicable for two of the four blend curve sets; however, two sets were found to exhibit slightly quadratic behavior, which was more pronounced in low biogenic blend samples, necessitating second-order fits. The origin of the slight quadratic behavior remains unclear; however, the discussion points to possible interpretations. CanmetENERGY thoroughly characterized a majority of the samples using one- and two-dimensional gas chromatography (GC and GC×GC, respectively) and other analyses. Selected samples were subjected to solid phase extraction (SPE) for saturate, olefin, aromatic, and polar (SOAP) analysis, and the resulting solvent-diluted fractions containing saturates and aromatics were returned to Pacific Northwest National Laboratory (PNNL), where the solvent was removed via evaporation or physical separation using GC techniques. Characterization and separations provided an understanding of saturate and aromatic content, as well as boiling point ranges for each sample and sample fraction. Samples resulting from SPE were examined using EA-IRMS and gas chromatography combustion IRMS (GC-C-IRMS) analyses. Both approaches suggest that the range in values between end-members can be increased by selecting the paraffinic or aromatic fraction of the end-member or by selecting among individual compounds resulting from GC separation of the paraffinic fractions. Considerable work remains to put these approaches into practice and statistically validate the benefit for using a fraction or individual compound over bulk analysis of a sample. However, initial results suggest that separations provide advantages for samples having blend ratios of less than 10% biogenic blendstocks. 13 C results showed statistically similar biofuel blend results to those obtained at PNNL, although additional work is needed to obtain better reproducibility. Select samples were sent to Los Alamos National Laboratory (LANL) for IRMS measurements and Beta Analytics for AMS measurements. This work suggests that IRMS and AMS yield closely comparable results and in some circumstances, IRMS could serve as a surrogate for AMS. While additional work is needed to better resolve statistical advantages for separations and better show the comparable nature of IRMS and AMS in both the biogenic carbon analysis of bulk chemical classes, initial results from this study suggest that these should be pursued in order to proliferate this approach for quantifying biogenic carbon in transportation fuels to the refinery level, thereby potentially enabling process optimization in co-processing scenarios.

09 BIOMASS FUELS↗
DOE OSTI2026

Screen-Printed SHJ Solar Cells with Complex Silver Inks

Metallization using complex metal inks has gained significant research interest due to its cost-effectiveness and ability to achieve performance comparable to traditional nanoparticle pastes. This study introduces the use of complex silver (Ag) inks applied via industrial screen-printing for silicon heterojunction (SHJ) solar cell metallization. The printed Ag lines exhibit a contact resistivity on SHJ tin-doped indium oxide (ITO) surfaces as low as approximately 0.2-12 mO cm2. Photoluminescence imaging reveals minimal surface passivation degradation (iVoc < 3.5 mV), while scanning electron microscopy (SEM) shows a denser structure compared to Ag layer from nanoparticle pastes. The printed Ag grid features thin (approximately 1 micrometer), continuous fingers approximately 100-120 micrometer wide, significantly thinner than conventional approximately 20-30 micrometer fingers produced with nanoparticle-based pastes. Double printing achieves SHJ device efficiencies exceeding 20%, the highest reported for industrial solar cell precursors using this technology. These findings highlight the potential of complex Ag inks as a sustainable alternative to particle-based pastes, reducing Ag consumption and processing temperatures without compromising efficiency.

14 SOLAR ENERGY↗
DOE OSTI2025

Geophysical measurements of a demolition at NRF

On May 31 st , 2025, at approximately 14:40:00 UTC, the former Submarine First Westinghouse (SW1) structure housing a large crane was demolished at the Naval Reactors Facility (NRF) located within the perimeter of the Idaho National Laboratory (INL). Several targeted explosives were used to facilitate this demolition. The INL Seismic Monitoring Program (SMP) installed four temporary seismic instrument suites around the site of the demolition in order to monitor ground motion including measuring peak ground acceleration (PGA), peak ground velocity (PGV), and to calculate peak ground displacement (PGD). The demolition produced three types of waves: body waves, surface waves, and acoustic waves. The acoustic wave represents the largest signal received. The PGA, measured on the eastern component of the instrument located approximately 95m away from the demolition site, was 0.187g. The PGV measured was 0.0004 m/s, and the PGD was approximately 0.000006m. Rayleigh waves incident on an instrument located at the Advanced Test Reactor, about 7.6km away, were used to estimate a magnitude of 0.3 at the ATR, representing a yield of energy equivalent to approximately 3 grams of TNT.

58 - GEOSCIENCES↗
DOE OSTI2026

HDG-1 Experiment Irradiation Monitoring Data Qualification Final Report

SUMMARY The U.S. Department of Energy (DOE) Advanced Reactor Technologies (ART) Graphite Research and Development (GRD) Program is conducting a series of six experiments to quantify the effects of irradiation on nuclear-grade graphite. This report documents the qualification of irradiation monitoring data for the fifth experiment, High Dose Graphite-1 (HDG-1). Qualified monitoring data are required by the ART program to support the design and licensing of the first high-temperature reactor (HTR) nuclear plant. Data are classified as Qualified if they meet the usage requirements described in the experiment planning and quality assurance (QA) documents, Failed if they do not meet those requirements and provide no usable information, or Trend if they do not fully meet all requirements but still provide useful information subject to an assessment of how any deficiencies may affect a particular use of the data. HDG-1 irradiation began with Advanced Test Reactor (ATR) Cycle 168B on August 24, 2020, and concluded after Cycle 173C on January 27, 2025. The HDG-1 capsule was removed from the reactor core twice—during core internal change (CIC) Cycle 170A and powered axial locator mechanism (PALM) Cycle 172A—to prevent overheating of the graphite specimens during high-power PALM cycles. The capsule was therefore irradiated during a total of seven normal ATR cycles: 168B, 169A, 171A, 171B, 173A, 173B, and 173C. Irradiation monitoring data evaluated in this report include thermocouple (TC) temperature, gas flow rate, gas moisture, gas pressure, specimen load, and graphite stack displacement. Temperature. A total of 14,508,065 TC temperature records were captured. Of these, 13,901,785 (95.8%) are Qualified and 606,280 (4.2%) are Failed. The principal source of failed temperature data was the instrument failure of TC-9 (Zone 2) on June 24, 2024, and TC-10 (Zone 1) on July 5, 2024, near the end of Cycle 173A, which resulted in 595,554 Failed readings. An additional 379 missing values and 10,347 slightly negative values from TC-13 during ATR outages are also Failed. Neither TC-9 nor TC-10 was used as a temperature-control TC, and their failures did not compromise capsule condition monitoring. Correlation analysis of all 13 TCs found no evidence of virtual junction formation. Control chart analysis revealed clear downward drift of approximately 80°C for TC-6 (Zone 3) relative to other stable TCs, and possible downward drift of approximately 60°C for TC-13 relative to the Zone 5 control TC (TC-1), though TC-13 remained consistent with the Zone 2 control TC (TC-12). Gas flow. A total of 20,088,090 gas flow rate records were captured. Of these, 19,941,463 (99.3%) are Qualified and 146,627 (0.7%) are Failed due to missing values. All argon, helium, and total gas flow data were within expected ranges throughout the irradiation. Gas moisture. A total of 1,116,005 outlet gas moisture values were captured. Of these, 1,101,421 (98.7%) are Qualified and 14,584 (1.3%) are Failed, comprising 14,556 out-of-range values and 28 missing values. The out-of-range moisture values exceeded 22,000 ppmv for approximately 1 week at the beginning of Cycle 173A, when accumulated moisture evaporated after the capsule was retrieved from water storage during PALM Cycle 172A and reinserted into the east flux trap. Moisture levels returned to below 25 ppmv for the remaining three cycles, and the transient high-moisture event did not affect the integrity of specimen irradiation. Gas pressure. A total of 7,812,035 gas pressure values were captured. Of these, 6,642,048 (85.0%) are Qualified and 1,169,987 (15.0%) outlet pressure values are Failed, comprising 718,537 zero outlet pressure values due to sensor failure from Cycle 168B through Cycle 171B, 54,550 missing values, and 396,900 too-low outlet pressure values, ranging from 1.1 to 1.6 psia after sensor replacement during Cycle 173A. Load. A total of 6,696,030 load values were captured. Of these, 6,694,580 (99.98%) are Qualified and 1,450 (0.02%) are Failed due to missing values. Applied loads to the six specimen stacks were stable throughout the irradiation. Stack displacement. A total of 6,696,030 displacement values were captured. Of these, 5,713,297 (85.32%) are Qualified and 3,781 (0.06%) are Failed due to missing values. Stack displacement increased consistently throughout the irradiation, reaching approximately 3.08 in. for Channels 5 and 6 by the end of irradiation. 978,952 (14.62%) substantially elevated displacements observed for Channel 6 beginning in Cycle 171A and for Channel 5 beginning in Cycle 173A are assigned Trend status. Raising pressure. A total of 1,115,999 raising pressure values were captured. Of these, 1,115,430 (99.95%) are Qualified and 569 (0.05%) are Failed due to missing values. Ram pressure. A total of 6,696,030 ram pressure values were captured. Of these, 6,692,249 (99.95%) are Qualified and 3,484 (0.05%) are Failed due to missing values. Stack raising was perf

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Image-based modeling of coupled electro-chemo-mechanical behavior of Li-ion battery cathode using an interface-modified reproducing kernel particle method

Abstract An interface-modified reproducing kernel particle method (IM-RKPM) is introduced in this work to allow for a direct model construction from image pixels of heterogeneous polycrystalline Li-ion battery microstructures. The interface-modified reproducing kernel (IM-RK) approximation is constructed through scaling of a kernel function by a regularized distance function in conjunction with strategic placement of interface node locations. This leads to RK shape functions with either weak or strong discontinuities across material interfaces, suitable for modeling various interface mechanics. With the placement of a triple junction node and distance-based scaling of kernel functions, the resulting IM-RK shape function also possesses proper discontinuities at the triple junctions. This IM-RK approximation effectively remedies the well-known Gibb’s oscillation in the smooth approximation of discontinuities. Different from the conventional meshfree approaches for interface discontinuities, this IM-RK approach is done without additional degrees of freedom associated with the enrichment functions, and it is formulated with the standard procedures in the RK shape function construction. This work focuses on identifying the accuracy and convergence properties of IM-RKPM for modeling the coupled electro-chemo-mechanical system. A linear patch test is formulated and numerically tested for the electro-chemo-mechanical coupled problem with a Butler–Volmer boundary condition representing the physical conditions in Li-ion battery microstructures. This is followed by verification of the optimal rates of convergence of IM-RKPM for solving the coupled problem with higher order solutions. The image-based modeling of Li-ion battery microstructures in the numerical examples demonstrates the applicability of the proposed method to realistic Li-ion battery materials modeling.

25 ENERGY STORAGE↗
DOE OSTI2025

McCormick envelopes in mixed-integer PDE-constrained optimization

McCormick envelopes are a standard tool for deriving convex relaxations of optimization problems that involve polynomial terms. Such McCormick relaxations provide lower bounds, for example, in branch-and-bound procedures for mixed-integer nonlinear programs but have not gained much attention in PDE-constrained optimization so far. This lack of attention may be due to the distributed nature of such problems, which on the one hand leads to infinitely many linear constraints (generally state constraints that may be difficult to handle) in addition to the state equation for a pointwise formulation of the McCormick envelopes and renders bound-tightening procedures that successively improve the resulting convex relaxations computationally intractable. We analyze McCormick envelopes for a model problem class that is governed by a semilinear PDE involving a bilinearity and integrality constraints. We approximate the nonlinearity and in turn the McCormick envelopes by averaging the involved terms over the cells of a partition of the computational domain on which the PDE is defined. This yields convex relaxations that underestimate the original problem up to an a priori error estimate that depends on the mesh size of the discretization. These approximate McCormick relaxations can be improved by means of an optimization-based bound-tightening procedure. We show that their minimizers converge to minimizers to a limit problem with a pointwise formulation of the McCormick envelopes when driving the mesh size to zero. We provide a computational example, for which we certify all of our imposed assumptions. The results point to both the potential of the methodology and the gaps in the research that need to be closed. Our methodology provides a framework first for obtaining pointwise underestimators for nonconvexities and second for approximating them with finitely many linear inequalities in an infinite-dimensional setting.

Approximations and Expansions↗
DOE OSTI2025

Slender-body approach for computing second-order wave loads in the frequency domain

This work presents a slender-body approach to evaluate the second-order wave loads acting on a floating structure in the frequency domain. The approach is in the same spirit as the common use of Morison’s equation to approximate the wave loads without solving the radiation/diffraction problem. To do so, we employ Rainey’s equation, which can be seen as an extension of the inertial part of Morison’s equation to include nonlinear effects. We introduce modifications to Rainey’s formulation in order to evaluate wave kinematics at the mean body position instead of the original approach of considering instantaneous displacements. We also propose a simple approximation to partially account for wave scattering effects on the second-order loads based on the analytical solution of a surface-piercing bottom-mounted vertical circular cylinder. Though limited to structures composed of cylinders, this slender-body approach is orders of magnitude faster than computing second-order wave coefficients with a radiation/diffraction code. We implemented this approach for difference-frequency (slow drift) loads in an open-source frequency-domain floating wind turbine model. We present comparisons against results obtained with radiation/diffraction theory for three reference floating wind turbine designs: the OC3-Hywind spar, the OC4-DeepCwind semisubmersible, and the VolturnUS-S semisubmersible. In general, the results show that the proposed slender-body approach with the correction to approximate wave scattering effects provides useful estimations of the difference-frequency wave loads and the resulting motions of the floater.

17 WIND ENERGY↗
DOE OSTI2024

Monte Carlo Explicitly Correlated Second-Order Many-Body Green’s Function Calculations of Semiconductor Band Gaps

A systematically converging series of ab initio, post-density-functional, size-consistent, electron-correlated approximations is desired for predictive computing of felectronic band structures of insulating, semiconducting, and metallic solids. A series that meets all of these desiderata (except the applicability to metals) is ab initio many-body Green's function theory based on Gaussian-type-orbital (GTO) basis sets. Here, its leading-order approximation, the second-order Green's function (GF2) method in the diagonal and frequency-independent approximations with the aug-cc-pVDZ basis set, is applied to the fundamental band gaps of three semiconductors (diamond, silicon, and silicon carbide in the zincblende structure) using cluster models. Corrections are made to the basis-set-incompleteness errors by the explicit-correlation (F12) ansatz (GF2-F12) for the valence band edges. The crystals are modeled as surface-passivated clusters of increasing sizes, whose wave functions are expanded by up to 2709 GTO basis functions. Immense computational costs of these calculations are overcome by the highly scalable stochastic algorithm of the Monte Carlo GF2-F12 method, whose operation cost per state increases only as a cubic power of system size, which has a tiny memory footprint and easily achieves near-perfect parallel efficiency on thousands of CPUs or on hundreds of GPUs. The correlated, F12-corrected highest-occupied and lowest-unoccupied molecular-orbital energy (HOMO-LUMO) gap is 5.78 ± 0.07 eV for C 87 H 76 as compared with the experimental value of the fundamental (indirect) band gap of bulk diamond at 5.48 eV. The correlated, F12-corrected HOMO-LUMO gaps for Si 75 H 76 and Si 32 C 43 H 76 are 2.56 ± 0.15 eV and 3.50 ± 0.12 eV, respectively, which are expected to decrease further with increasing cluster sizes. As a result, the experimental fundamental (indirect) band gaps of bulk silicon and silicon carbide are 1.17 eV and 2.42 eV, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗